On properties of compact 4th order finite-difference schemes for the variable coefficient wave equation

2021. No. ArXiv: 2101.10575v2[math.NA].

Zlotnik A., Čiegis R.

We consider an initial-boundary value problem for the $n$-dimensional wave equation with the variable sound speed, $n\geq 1$. We construct three-level implicit in time compact in space (three-point in each space direction) 4th order finite-difference schemes on the uniform rectangular meshes including their one-parameter (for $n=2$) and three-parameter (for $n=3$) families. They are closely connected to some methods and schemes constructed recently by several authors. In a unified manner, we prove the conditional stability of schemes in the strong and weak energy norms together with the 4th order error estimate under natural conditions on the time step. We also give an example of extending a compact scheme for non-uniform in space and time rectangular meshes. We suggest simple effective iterative methods based on FFT to implement the schemes whose convergence rate, under the stability condition, is fast and independent on both the meshes and variable sound speed. A new effective initial guess to start iterations is given too. We also present promising results of numerical experiments.

Keywords: Wave equation, stability, iterative methods, variable speed of sound, compact higher-order scheme

Publication based on the results of: Теоретический анализ численных методов повышенного порядка для решения задач волновой физики (2019)