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

We 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.