Abstract
We introduce solution dependent finite difference stencils whose coefficients adapt to the current numerical solution by minimizing the truncation error in the least squares sense. The resulting scheme has the resolution capacity of dispersion relation preserving difference stencils in under-resolved domains, together with the high order convergence rate of conventional central difference methods in well resolved regions. Numerical experiments reveal that the new stencils outperform their conventional counterparts on all grid resolutions from very coarse to very fine.
Original language | English |
---|---|
Article number | 109393 |
Journal | Journal of Computational Physics |
Volume | 410 |
DOIs | |
Publication status | Published - 1 Jun 2020 |
Externally published | Yes |
Keywords
- Accuracy
- Adaptivity
- Convergence
- Dispersion relation preserving
- Finite differences
- Least squares
ASJC Scopus subject areas
- Numerical Analysis
- Modeling and Simulation
- Physics and Astronomy (miscellaneous)
- General Physics and Astronomy
- Computer Science Applications
- Computational Mathematics
- Applied Mathematics