Abstract
We derive a method to locally change the order of accuracy of finite difference schemes that approximate the second derivative. The derivation is based on summation-by-parts operators, which are connected at interfaces using penalty terms. At such interfaces, the numerical solution has a double representation, with one representation in each domain. We merge this double representation into a single one, yielding a new scheme with unique solution values in all grid points. The resulting scheme is proven to be stable, accurate and dual consistent.
Original language | English |
---|---|
Pages (from-to) | 935-949 |
Number of pages | 15 |
Journal | Journal of Computational Physics |
Volume | 375 |
DOIs | |
Publication status | Published - 15 Dec 2018 |
Externally published | Yes |
Keywords
- Dual consistency
- Finite difference methods
- High order accuracy
- Interfaces
- Summation-by-parts
- Superconvergence
ASJC Scopus subject areas
- Numerical Analysis
- Modeling and Simulation
- Physics and Astronomy (miscellaneous)
- General Physics and Astronomy
- Computer Science Applications
- Computational Mathematics
- Applied Mathematics