Convergence of Summation-by-Parts Finite Difference Methods for the Wave Equation

Abstract

When using a finite difference method to solve a time dependent partial differential equation, the truncation error is often larger at a few grid points near a boundary or grid interface than in the interior. In computations, the observed convergence rate is often higher than the order of the large truncation error. In this paper, we develop techniques for analyzing this phenomenon, and particularly consider the second order wave equation. The equation is discretized by a finite difference operator satisfying a summation by parts property, and the boundary and grid interface conditions are imposed weakly by the simultaneous approximation term method. It is well-known that if the semi-discretized wave equation satisfies the determinant condition, that is the boundary system in Laplace space is nonsingular for all Re(s) ≥ 0 , two orders are gained from the large truncation error localized at a few grid points. By performing a normal mode analysis, we show that many common discretizations do not satisfy the determinant condition at s= 0. We then carefully analyze the error equation to determine the gain in the convergence rate. The result shows that stability does not automatically imply a gain of two orders in the convergence rate. The precise gain can be lower than, equal to or higher than two orders, depending on the boundary condition and numerical boundary treatment. The accuracy analysis is verified by numerical experiments, and very good agreement is obtained.