Provably Stable Full-Spectrum Dispersion Relation Preserving Schemes

Abstract: The dispersion error is often the dominant error for computed solutions of wave propagation problems with high-frequency components. In this paper, we define and give explicit examples of $\alpha$-dispersion-relation-preserving schemes. These are dual-pair finite-difference schemes for systems of hyperbolic partial differential equations which preserve the dispersion-relation of the continuous problem uniformly to an $\alpha \%$-error tolerance. We give a general framework to design provably stable finite difference operators that preserve the dispersion relation for hyperbolic systems such as the elastic wave equation. The operators we derive here can resolve the highest frequency ($\pi$-mode) present on any equidistant grid at a tolerance of $5\%$ error. This significantly improves on the current standard that have a tolerance of $100 \%$ error.