Premium
Time‐domain hybrid formulations for wave simulations in three‐dimensional PML‐truncated heterogeneous media
Author(s) -
Fathi Arash,
Poursartip Babak,
Kallivokas Loukas F.
Publication year - 2014
Publication title -
international journal for numerical methods in engineering
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.421
H-Index - 168
eISSN - 1097-0207
pISSN - 0029-5981
DOI - 10.1002/nme.4780
Subject(s) - stability (learning theory) , perfectly matched layer , displacement (psychology) , field (mathematics) , displacement field , domain (mathematical analysis) , mathematics , scheme (mathematics) , computer science , mathematical analysis , finite element method , pure mathematics , engineering , boundary value problem , structural engineering , psychology , machine learning , psychotherapist
SUMMARY 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.