Hybrid perfectly-matched-layers for transient simulation of scalar elastic waves

Abstract

This paper presents a new formulation for forward scalar wave simulations in semi-infinite media. Perfectly-Matched-Layers (PMLs) are used as a wave absorbing boundary layer to surround a finite computational domain truncated from the semi-infinite domain. In this work, a hybrid formulation was developed for the simulation of scalar wave motion in two-dimensional PML-truncated domains. In this formulation, displacements and stresses are considered as unknowns in the PML domain, while only displacements are considered to be unknowns in the interior domain. This formulation reduces computational cost compared to fully-mixed formulations. To obtain governing wave equations in the PML region, complex coordinate stretching transformation was introduced to equilibrium, constitutive, and compatibility equations in the frequency domain. Then, equations were converted back to the time-domain using the inverse Fourier transform. The resulting equations are mixed (contain both displacements and stresses), and are coupled with the displacement-only equation in the regular domain. The Newmark method was used for the time integration of the semi-discrete equations.