Elastic wave propagation in variable media using a discontinuous Galerkin method.
Motivated by the needs of seismic inversion and building on our prior experience for fluid-dynamics systems, we present a high-order discontinuous Galerkin (DG) Runge-Kutta method applied to isotropic, linearized elasto-dynamics. Unlike other DG methods recently presented in the literature, our method allows for inhomogeneous material variations within each element that enables representation of realistic earth models - a feature critical for future use in seismic inversion. Likewise, our method supports curved elements and hybrid meshes that include both simplicial and nonsimplicial elements. We demonstrate the capabilities of this method through a series of numerical experiments including hybrid mesh discretizations of the Marmousi2 model as well as a modified Marmousi2 model with a oscillatory ocean bottom that is exactly captured by our discretization. A discontinuous Galerkin method for solving the equations of linear isotropic elasticity has been presented. The formulation is designed to accommodate variation of media parameters within elements, curved elements and unstructured heterogeneous meshes. We have demonstrated that each of these important features of the formulation can produce results that are significantly different from formulations that do not possess these capabilities suggesting that each of these capabilities may be important for effective full waveform inversion of elastic medium.
- Research Organization:
- Sandia National Laboratories
- Sponsoring Organization:
- USDOE
- DOE Contract Number:
- AC04-94AL85000
- OSTI ID:
- 992317
- Report Number(s):
- SAND2010-2166C
- Country of Publication:
- United States
- Language:
- English
Similar Records
CosmosDG: An hp -adaptive Discontinuous Galerkin Code for Hyper-resolved Relativistic MHD
A multiscale discontinuous Galerkin method.