Abstract
We consider numerical simulation of the isotropic elastic wave equations arising from seismic applications with non-trivial land topography. The more flexible finite element method is applied to the shallow region of the simulation domain to account for the topography, and combined with the more efficient finite difference method that is applied to the deep region of the simulation domain. We demonstrate that these two discretization methods, albeit starting from different formulations of the elastic wave equation, can be joined together smoothly via weakly imposed interface conditions. Discrete energy analysis is employed to derive the proper interface treatment, leading to an overall discretization that is energy-conserving. Numerical examples are presented to demonstrate the efficacy of the proposed interface treatment.
| Original language | English (US) |
|---|---|
| Pages (from-to) | 665-685 |
| Number of pages | 21 |
| Journal | Journal of Computational Physics |
| Volume | 378 |
| DOIs | |
| State | Published - Nov 23 2018 |
Bibliographical note
KAUST Repository Item: Exported on 2020-10-01Acknowledged KAUST grant number(s): CCF-CAF/URF/1-2596
Acknowledgements: The authors gratefully acknowledge the support of KAUST's Office of Sponsored Research under CCF-CAF/URF/1-2596. The authors would also like to thank the anonymous reviewers for their thoughtful suggestions and comments that have led to significant improvements in this article.