Absorbing boundary conditions and geometric integration: A case study for the wave equation

This paper is concerned about the confluence of two subjects of the numerical solution of time evolution PDEs: numerical methods that preserve geometric properties of the flow and the use of absorbing boundary conditions to reduce the computation to a finite domain. This confluence is studied with special attention to the time stability of the resulting full discretization. For this, the stability regions of the time integrators are revisited. Since geometric methods are not always A-stable, it is necessary a suitable behavior of the real part of the eigenvalues of the spatially discretized problem to avoid in practice any time instability. A deep study is carried out for the case of the one dimensional wave equation when it is discretized with finite differences, showing that this suitable behavior happens. Numerical experiments confirming the previous results are included.

[1]  E. Hairer,et al.  Solving Ordinary Differential Equations II: Stiff and Differential-Algebraic Problems , 2010 .

[2]  Laurence Halpern,et al.  Absorbing boundary conditions for the discretization schemes of the one-dimensional wave equation , 1982 .

[3]  Isaías Alonso-Mallo,et al.  High order full discretizations of coupled wave equations with absorbing boundary conditions and geometric integration , 2014, J. Comput. Phys..

[4]  R. McLachlan.  Symplectic integration of Hamiltonian wave equations , 1993 .

[5]  J. M. Sanz-Serna,et al.  Order conditions for canonical Runge-Kutta schemes , 1991 .

[6]  Thomas Hagstrom,et al.  New Results on Absorbing Layers and Radiation Boundary Conditions , 2003 .

[7]  H. Fattorini.  Second Order Linear Differential Equations in Banach Spaces , 1985 .

[8]  T. Hagstrom.  Radiation boundary conditions for the numerical simulation of waves , 1999, Acta Numerica.

[9]  Lloyd N. Trefethen,et al.  Pseudospectra of Linear Operators , 1997, SIAM Rev..

[10]  Brian E. Moore.  A modified equations approach for multi-symplectic integration methods. , 2003 .

[11]  D. Givoli.  High-order local non-reflecting boundary conditions: a review☆ , 2004 .

[12]  E. Hairer,et al.  Geometric Numerical Integration: Structure Preserving Algorithms for Ordinary Differential Equations , 2004 .

[13]  Isaías Alonso-Mallo,et al.  A proof of the well posedness of discretized wave equation with an absorbing boundary condition , 2014, J. Num. Math..

[14]  S. Tsynkov.  Numerical solution of problems on unbounded domains. a review , 1998 .

[15]  S. Reich,et al.  Numerical methods for Hamiltonian PDEs , 2006 .

[16]  Siam Staff.  Weak Ill-Posedness of Spatial Discretizations of Absorbing Boundary Conditions for Schrödinger-Type Equations , 2002 .

[17]  Dan Givoli,et al.  High-order local absorbing conditions for the wave equation: Extensions and improvements , 2008, J. Comput. Phys..

[18]  Fermín S. Viloche Bazán,et al.  Chebyshev pseudospectral method for wave equation with absorbing boundary conditions that does not use a first order hyperbolic system , 2010, Math. Comput. Simul..

[19]  J. M. Sanz-Serna,et al.  An easily implementable fourth-order method for the time integration of wave problems , 1992 .

[20]  Isaías Alonso-Mallo,et al.  Discrete Absorbing Boundary Conditions for Schrödinger-Type Equations. Construction and Error Analysis , 2003, SIAM J. Numer. Anal..

[21]  J. M. Sanz-Serna,et al.  Numerical Hamiltonian Problems , 1994 .