Time‐domain hybrid formulations for wave simulations in three‐dimensional PML‐truncated heterogeneous media

Abstract

SUMMARYWe are concerned with the numerical simulation of wave motion in arbitrarily heterogeneous, elastic, perfectly‐matched‐layer‐(PML)‐truncated media. We extend in three dimensions a recently developed two‐dimensional formulation, by treating the PML via an unsplit‐field, but mixed‐field, displacement‐stress formulation, which is then coupled to a standard displacement‐only formulation for the interior domain, thus leading to a computationally cost‐efficient hybrid scheme. The hybrid treatment leads to, at most, third‐order in time semi‐discrete forms. The formulation is flexible enough to accommodate the standard PML, as well as the multi‐axial PML.We discuss several time‐marching schemes, which can be used à la carte, depending on the application: (a) an extended Newmark scheme for third‐order in time, either unsymmetric or fully symmetric semi‐discrete forms; (b) a standard implicit Newmark for the second‐order, unsymmetric semi‐discrete forms; and (c) an explicit Runge–Kutta scheme for a first‐order in time unsymmetric system. The latter is well‐suited for large‐scale problems on parallel architectures, while the second‐order treatment is particularly attractive for ready incorporation in existing codes written originally for finite domains.We compare the schemes and report numerical results demonstrating stability and efficacy of the proposed formulations. Copyright © 2014 John Wiley & Sons, Ltd.

Non-Contact Ultrasonic Guided Wave Inspection of Rails: Next Generation Approach

“…The PML technique was introduced in 1994 by Berenger for electromagnetism [92], and later has been extended to other fields such as acoustics [93] and other elastic waves [94,95,96]. PML elements perfectly match the impedance of the material of study, so that waves enter in PML without reflections (at least in theory; in practice, studies have shown that some numerical reflections occur, mostly due to insufficient numbers of element per wavelength employed in the model [97]).…”

Source parameter inversion for wave energy focusing to a target inclusion embedded in a three‐dimensional heterogeneous halfspace

“…We are concerned with stress wave propagation in a three‐dimensional, heterogeneous, elastic halfspace containing a target inclusion. In order to obtain a finite computational model, we truncate the domain of interest ( ) using hybrid PMLs (Ω PML ) . Note that in Figure , Ω a represents the target inclusion and Ω b represents the heterogeneous elastic solid surrounding the target inclusion.…”

Hybrid Asynchronous Perfectly Matched Layer for seismic wave propagation in unbounded domains

“…From the frequency-domain equations of Basu and Chopra, Kucukcoban and Kallivokas derived an unsplit mixed approach of the PML, by retaining the displacement and stress fields as unknowns in the time domain [17]. Next, in order to couple their PML to a displacement-only field formulation in the physical interior domain, the authors extended their PML to a mixed hybrid approach [18,19].…”

“…The main benefit is a higher versatility and numerical efficiency of the PML which can be implemented using a more appropriate time integrator associated with a larger time step than the one employed in the soil medium imposed by the CFL condition for ensuring the algorithm stability [20]. The proposed PMLs can be viewed as a hybrid and asynchronous PML version because different time integrators can be adopted as in the work of Kucukcoban and Kallivokas [18,19], with the adding desirable properties of dealing with different time scales in the same dynamic simulation. Moreover, the proposed approach enables the number of unknowns to be reduced in comparison to the previous mixed displacement-stress formulation.…”