Licentiate thesis 2015-005

Analysis of Boundary and Interface Closures for Finite Difference Methods for the Wave Equation

Siyang Wang

27 October 2015

Abstract:

We consider high order finite difference methods for the wave equations in the second order form, where the finite difference operators satisfy the summation-by-parts principle. Boundary conditions and interface conditions are imposed weakly by the simultaneous-approximation-term method, and non-conforming grid interfaces are handled by an interface operator that is based on either interpolating directly between the grids or on projecting to piecewise continuous polynomials on an intermediate grid.

Stability and accuracy are two important aspects of a numerical method. For accuracy, we prove the convergence rate of the summation-by-parts finite difference schemes for the wave equation. Our approach is based on Laplace transforming the error equation in time, and analyzing the solution to the boundary system in the Laplace space. In contrast to first order equations, we have found that the determinant condition for the second order equation is less often satisfied for a stable numerical scheme. If the determinant condition is satisfied uniformly in the right half plane, two orders are recovered from the boundary truncation error; otherwise we perform a detailed analysis of the solution to the boundary system in the Laplace space to obtain an error estimate. Numerical experiments demonstrate that our analysis gives a sharp error estimate.

For stability, we study the numerical treatment of non-conforming grid interfaces. In particular, we have explored two interface operators: the interpolation operators and projection operators applied to the wave equation. A norm-compatible condition involving the interface operator and the norm related to the SBP operator is essential to prove stability by the energy method for first order equations. In the analysis, we have found that in contrast to first order equations, besides the norm-compatibility condition an extra condition must be imposed on the interface operators to prove stability by the energy method. Furthermore, accuracy and efficiency studies are carried out for the numerical schemes.

Available as PDF (11.79 MB)

Download BibTeX entry.