The double absorbing boundary method for elastodynamics in homogeneous and layered media

Abstract

Background: Recently the Double Absorbing Boundary (DAB) method was introduced as a new approach for solving wave problems in unbounded domains. It has common features to each of two types of existing techniques: local high-order Absorbing Boundary Conditions (ABC) and Perfectly Matched Layers (PML). However, it is different from both and enjoys relative advantages with respect to both. Methods: The DAB method is based on truncating the unbounded domain to produce a finite computational domain, and on applying a local high-order ABC on two parallel artificial boundaries, which are a small distance apart, and thus form a thin non-reflecting layer. Auxiliary variables are defined on the two boundaries and within the layer, and participate in the numerical scheme. In previous studies DAB was developed for acoustic waves which are solutions to the scalar wave equation. Here the approach is extended to time-dependent elastic waves in homogeneous and layered media. The equations are written in second-order form in space and time. Standard Finite Elements (FE) are used for space discretization and the damped Newmark scheme is used for time discretization. Results: The performance of the scheme is demonstrated via numerical examples. The DAB was applied to elastodynamics problems in conjunction with the FE method to demonstrate the performance of the method. Conclusions: DAB is a viable method for solving wave problems in unbounded domains.