Stable and high order accurate difference methods for the elastic wave equation in discontinuous media

Kenneth Duru, Kristoffer Virta

Research output: Contribution to journalArticlepeer-review

18 Scopus citations

Abstract

© 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.
Original languageEnglish (US)
Pages (from-to)37-62
Number of pages26
JournalJournal of Computational Physics
Volume279
DOIs
StatePublished - Dec 2014
Externally publishedYes

Fingerprint

Dive into the research topics of 'Stable and high order accurate difference methods for the elastic wave equation in discontinuous media'. Together they form a unique fingerprint.

Cite this