Combining finite element and finite difference methods for isotropic elastic wave simulations in an energy-conserving manner
KAUST DepartmentExtreme Computing Research Center
Computer, Electrical and Mathematical Sciences and Engineering (CEMSE) Division
Applied Mathematics and Computational Science Program
KAUST Grant NumberCCF-CAF/URF/1-2596
Permanent link to this recordhttp://hdl.handle.net/10754/627227
MetadataShow full item record
AbstractWe 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.
CitationGao L, Keyes D (2019) Combining finite element and finite difference methods for isotropic elastic wave simulations in an energy-conserving manner. Journal of Computational Physics 378: 665–685. Available: http://dx.doi.org/10.1016/j.jcp.2018.11.031.
SponsorsThe 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.
JournalJournal of Computational Physics