TY - JOUR
T1 - Stable and high order accurate difference methods for the elastic wave equation in discontinuous media
AU - Duru, Kenneth
AU - Virta, Kristoffer
N1 - KAUST Repository Item: Exported on 2020-10-01
Acknowledgements: The work of the first author was supported by King Abdullah University of Science and Technology (KAUST) through a joint KAUST Academic Excellence Alliance (AEA) grant with Stanford. The first author also acknowledges the support of Eric M. Dunham during this work.
This publication acknowledges KAUST support, but has no KAUST affiliated authors.
PY - 2014/12
Y1 - 2014/12
N2 - © 2014 Elsevier Inc. In this paper, we develop a stable and systematic procedure for numerical treatment of elastic waves in discontinuous and layered media. We consider both planar and curved interfaces where media parameters are allowed to be discontinuous. The key feature is the highly accurate and provably stable treatment of interfaces where media discontinuities arise. We discretize in space using high order accurate finite difference schemes that satisfy the summation by parts rule. Conditions at layer interfaces are imposed weakly using penalties. By deriving lower bounds of the penalty strength and constructing discrete energy estimates we prove time stability. We present numerical experiments in two space dimensions to illustrate the usefulness of the proposed method for simulations involving typical interface phenomena in elastic materials. The numerical experiments verify high order accuracy and time stability.
AB - © 2014 Elsevier Inc. In this paper, we develop a stable and systematic procedure for numerical treatment of elastic waves in discontinuous and layered media. We consider both planar and curved interfaces where media parameters are allowed to be discontinuous. The key feature is the highly accurate and provably stable treatment of interfaces where media discontinuities arise. We discretize in space using high order accurate finite difference schemes that satisfy the summation by parts rule. Conditions at layer interfaces are imposed weakly using penalties. By deriving lower bounds of the penalty strength and constructing discrete energy estimates we prove time stability. We present numerical experiments in two space dimensions to illustrate the usefulness of the proposed method for simulations involving typical interface phenomena in elastic materials. The numerical experiments verify high order accuracy and time stability.
UR - http://hdl.handle.net/10754/599717
UR - https://linkinghub.elsevier.com/retrieve/pii/S0021999114006226
UR - http://www.scopus.com/inward/record.url?scp=84907546178&partnerID=8YFLogxK
U2 - 10.1016/j.jcp.2014.08.046
DO - 10.1016/j.jcp.2014.08.046
M3 - Article
SN - 0021-9991
VL - 279
SP - 37
EP - 62
JO - Journal of Computational Physics
JF - Journal of Computational Physics
ER -