Difference Approximations of the Neumann Problem for the Second Order Wave Equation

Stability theory and numerical experiments are presented for a finite difference method that directly discretizes the Neumann problem for the second order wave equation. Complex geometries are discretized using a Cartesian embedded boundary technique. Both second and third order accurate approximati...

Full description

Saved in:
Bibliographic Details
Published inSIAM journal on numerical analysis Vol. 42; no. 3; pp. 1292 - 1323
Main Authors Kreiss, Heinz-Otto, Petersson, N. Anders, Yström, Jacob
Format Journal Article
LanguageEnglish
Published Philadelphia Society for Industrial and Applied Mathematics 01.01.2004
Subjects
Accuracy
Approximation
Boundary conditions
embedded boundary
Laboratories
Methods
Numerical analysis
stability
wave equation
Online AccessGet full text

Cover

More Information
Summary:Stability theory and numerical experiments are presented for a finite difference method that directly discretizes the Neumann problem for the second order wave equation. Complex geometries are discretized using a Cartesian embedded boundary technique. Both second and third order accurate approximations of the boundary conditions are presented. Away from the boundary, the basic second order method can be corrected to achieve fourth order spatial accuracy. To integrate in time, we present both a second order and a fourth order accurate explicit method. The stability of the method is ensured by adding a small fourth order dissipation operator, locally modified near the boundary to allow its application at all grid points inside the computational domain. Numerical experiments demonstrate the accuracy and long-time stability of the proposed method.
ISSN:0036-1429
1095-7170
1095-7170
DOI:10.1137/S003614290342827X