Convergence of finite difference methods for the wave equation in two space dimensions
Artikel i vetenskaplig tidskrift, 2018

When using a finite difference method to solve an initial-boundaryvalue problem, the truncation error is often of lower order at a few grid points near boundaries than in the interior. Normal mode analysis is a powerful tool to analyze the effect of the large truncation error near boundaries on the overall convergence rate, and has been used in many research works for different equations. However, existing work only concerns problems in one space dimension. In this paper, we extend the analysis to problems in two space dimensions. The two dimensional analysis is based on a diagonalization procedure that decomposes a two dimensional problem to many one dimensional problems of the same type. We present a general framework of analyzing convergence for such one dimensional problems, and explain how to obtain the result for the corresponding two dimensional problem. In particular, we consider two kinds of truncation errors in two space dimensions: the truncation error along an entire boundary, and the truncation error localized at a few grid points close to a corner of the computational domain. The accuracy analysis is in a general framework, here applied to the second order wave equation. Numerical experiments corroborate our accuracy analysis.

Accuracy

Normal mode analysis

Two space dimensions

Finite difference method

Second order wave equation

Convergence rate

Författare

Siyang Wang

Chalmers, Matematiska vetenskaper, Tillämpad matematik och statistik

Anna Nissen

Universitetet i Bergen

Kungliga Tekniska Högskolan (KTH)

Gunilla Kreiss

Uppsala universitet

Mathematics of Computation

0025-5718 (ISSN) 1088-6842 (eISSN)

Vol. 87 314 2737-2763

Ämneskategorier

Beräkningsmatematik

Annan matematik

Matematisk analys

DOI

10.1090/mcom/3319

Mer information

Senast uppdaterat

2018-09-18