Combining finite element and finite difference methods for isotropic elastic wave simulations in an energy-conserving manner

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.

Analysis of wastewater facultative pond using advection-diffusion model based on explicit finite difference method

“…Brewitt-Taylor et al [14] discussed about the use of finite difference method for micro-perforated panel absorber predicting, whereas Patil et al [15] applied the finite difference method to analyze bio-heat transfer in human breast cyst. Furthermore, Gao et al [16] applied finite difference method for isotropic elastic wave simulations. Ben-romdhane [17] solved the Bratu's problem by using finite difference method.…”

A Novel Meshfree Strategy for a Viscous Wave Equation With Variable Coefficients

“…The earliest numerical methods for wave equations include the finite difference method [8], the finite element method [9], the boundary element method [10] and the finite volume method [11]. It should be noted that most of the existed numerical methods are still based on the finite difference method (FDM) [12,13], which lead to two-step finite difference approximations. More specifically, the finite difference method is used to deal with the time variable, and the rest procedures are finished by the other numerical methods.…”

Energy conservative SBP discretizations of the acoustic wave equation in covariant form on staggered curvilinear grids

“…SBP methods have also been extend beyond finite differences, e.g., to finite volume, discontinuous Galerkin, and flux reconstruction schemes [18,5,29]. Methods have also been devised for coupling different SBP discretizations [19,14,12,4].…”