Numerical comparison of high-order absorbing boundary conditions and perfectly matched layers for a dispersive one-dimensional medium

Abstract

High-order absorbing boundary conditions (ABC) and perfectly matched layers (PML) are two powerful methods to numerically solve wave problems in unbounded domains. The aim of the proposed study is to analyze and compare the performance of these methods in the one-dimensional problem governed by the dispersive wave equation. The PMLs proposed in literature for time-harmonic dynamics are applied to time-dependent wave problems, and linear, quadratic and cubic polynomial stretching functions are considered. The resulting PMLs exhibit a double absorbing action: (i) they reduce the amplitude of the incident waves and (ii) they delay the wave propagation. Then the ABC proposed by Givoli and Neta and those proposed by Hagstrom, Mar-Or and Givoli are considered. The former are a reformulation of the Higdon high-order non-reflecting boundary conditions, the latter improve the Higdon conditions and extend them to take into account evanescent waves. The accuracy of the PMLs and ABCs, implemented in a finite element code, is first investigated with respect to the frequency of the incident wave, being it progressive or evanescent. Then the response to a wave train characterized by a broad frequency spectrum, resulting from an impulsive force, is studied. A detailed analysis is performed to detect the influence of the parameters of both the ABCs and the PMLs on the absorption of waves. The performances of PMLs and ABCs are compared, and merits and drawbacks of the two methods are pointed out.