High-order local absorbing conditions for the wave equation: Extensions and improvements

  • Authors:

  • Thomas Hagstrom;Assaf Mar-Or;Dan Givoli

  • Affiliations:

  • Department of Mathematics and Statistics, The University of New Mexico, Albuquerque, NM 87131, USA;Inter-Departmental Program for Applied Mathematics, Technion - Israel Institute of Technology, Haifa 32000, Israel;Department of Aerospace Engineering, Technion - Israel Institute of Technology, Haifa 32000, Israel

  • Venue:
  • Journal of Computational Physics
  • Year:
  • 2008

Quantified Score

Hi-index 31.49

Visualization

Abstract

The solution of the time-dependent wave equation in an unbounded domain is considered. An artificial boundary B is introduced which encloses a finite computational domain. On B an absorbing boundary condition (ABC) is imposed. A formulation of local high-order ABCs recently proposed by Hagstrom and Warburton and based on a modification of the Higdon ABCs, is further developed and extended in a number of ways. First, the ABC is analyzed in new ways and important information is extracted from this analysis. Second, The ABCs are extended to the case of a dispersive medium, for which the Klein-Gordon wave equation governs. Third, the case of a stratified medium is considered and the way to apply the ABCs to this case is explained. Fourth, the ABCs are extended to take into account evanescent modes in the exact solution. The analysis is applied throughout this paper to two-dimensional wave guides. Two numerical algorithms incorporating these ABCs are considered: a standard semi-discrete finite element formulation in space followed by time-stepping, and a high-order finite difference discretization in space and time. Numerical examples are provided to demonstrate the performance of the extended ABCs using these two methods.