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