A Fourth-Order Accurate Embedded Boundary Method for the Wave Equation

A fourth-order accurate embedded boundary method for the scalar wave equation with Dirichlet or Neumann boundary conditions is described. The method is based on a compact Pade-type discretization of spatial derivatives together with a Taylor series method (modified equation) in time. A novel approach for enforcing boundary conditions is introduced which uses interior boundary points instead of exterior ghost points. This technique removes the small-cell stiffness problem for both Dirichlet and Neumann boundary conditions, is more accurate and robust than previous methods based on exterior ghost points, and guarantees that the solution is single-valued when slender bodies are treated. Numerical experiments are presented to illustrate the stability and accuracy of the method as well as its application to problems with complex geometries.

[1]  John R. Cary,et al.  Extracting degenerate modes and frequencies from time-domain simulations with filter-diagonalization , 2008, J. Comput. Phys..

[2]  Mark Lyon,et al.  High-order unconditionally stable FC-AD solvers for general smooth domains I. Basic elements , 2010, J. Comput. Phys..

[3]  Leslie Greengard,et al.  An Integral Evolution Formula for the Wave Equation , 2000 .

[4]  Mark Lyon,et al.  High-Order Unconditionally-Stable FC-AD PDE Solvers for General Domains , 2009 .

[5]  M. A. Dablain,et al.  The application of high-order differencing to the scalar wave equation , 1986 .

[6]  Leslie Greengard,et al.  High order marching schemes for the wave equation in complex geometry , 2004 .

[7]  Heinz-Otto Kreiss,et al.  Difference Approximations of the Neumann Problem for the Second Order Wave Equation , 2004, SIAM J. Numer. Anal..

[8]  Magnus Svärd,et al.  Stable and Accurate Artificial Dissipation , 2004, J. Sci. Comput..

[9]  Jean Virieux,et al.  Free and smooth boundaries in 2-D finite-difference schemes for transient elastic waves , 2007, 0706.3825.

[10]  HEINZ-OTTO KREISS,et al.  A Second Order Accurate Embedded Boundary Method for the Wave Equation with Dirichlet Data , 2005, SIAM J. Sci. Comput..

[11]  Pelle Olsson,et al.  The numerical behavior of high-order finite difference methods , 1994 .

[12]  Patrick Joly,et al.  Construction and analysis of higher order finite difference schemes for the 1D wave equation , 2000 .

[13]  Joël Piraux,et al.  Numerical treatment of two-dimensional interfaces for acoustic and elastic waves , 2004 .

[14]  Björn Sjögreen,et al.  A Cartesian Embedded Boundary Method for the Compressible Navier-Stokes Equations , 2008, J. Sci. Comput..

[15]  H. Kreiss,et al.  Finite Element and Finite Difference Methods for Hyperbolic Partial Differential Equations , 1974 .

[16]  S. Lele.  Compact finite difference schemes with spectral-like resolution , 1992 .

[17]  Randall J. LeVeque,et al.  Finite difference methods for ordinary and partial differential equations - steady-state and time-dependent problems , 2007 .

[18]  Mark Lyon,et al.  High-order unconditionally stable FC-AD solvers for general smooth domains II. Elliptic, parabolic and hyperbolic PDEs; theoretical considerations , 2010, J. Comput. Phys..

[19]  John B. Bell,et al.  A modified equation approach to constructing fourth order methods for acoustic wave propagation , 1987 .

[20]  Mats Viberg,et al.  Estimation of resonant frequencies and quality factors from time domain computations , 2003 .

[21]  Heinz-Otto Kreiss,et al.  Difference Approximations for the Second Order Wave Equation , 2002, SIAM J. Numer. Anal..

[22]  H. Kreiss,et al.  Comparison of accurate methods for the integration of hyperbolic equations , 1972 .

[23]  Stephen M. Wandzura,et al.  Stable, high-order discretization for evolution of the wave equation in 2 + 1 dimensions , 2004 .

[24]  Gunilla Kreiss,et al.  Application of a perfectly matched layer to the nonlinear wave equation , 2007 .

[25]  C. Peskin.  Flow patterns around heart valves: A numerical method , 1972 .

[26]  L. Collatz.  The numerical treatment of differential equations , 1961 .

[27]  B. Sjögreen,et al.  A Cartesian embedded boundary method for hyperbolic conservation laws , 2006 .

[28]  John B. Bell,et al.  Cartesian grid method for unsteady compressible flow in irregular regions , 1995 .

[29]  Amanda M. White,et al.  Stable, high-order discretization for evolution of the wave equation in 1 + 1 dimensions , 2004 .