skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Super-Grid Modeling of the Elastic Wave Equation in Semi-Bounded Domains

Abstract

Abstract We develop a super-grid modeling technique for solving the elastic wave equation in semi-bounded two- and three-dimensional spatial domains. In this method, waves are slowed down and dissipated in sponge layers near the far-field boundaries. Mathematically, this is equivalent to a coordinate mapping that transforms a very large physical domain to a significantly smaller computational domain, where the elastic wave equation is solved numerically on a regular grid. To damp out waves that become poorly resolved because of the coordinate mapping, a high order artificial dissipation operator is added in layers near the boundaries of the computational domain. We prove by energy estimates that the super-grid modeling leads to a stable numerical method with decreasing energy, which is valid for heterogeneous material properties and a free surface boundary condition on one side of the domain. Our spatial discretization is based on a fourth order accurate finite difference method, which satisfies the principle of summation by parts. We show that the discrete energy estimate holds also when a centered finite difference stencil is combined with homogeneous Dirichlet conditions at several ghost points outside of the far-field boundaries. Therefore, the coefficients in the finite difference stencils need only be boundary modifiedmore » near the free surface. This allows for improved computational efficiency and significant simplifications of the implementation of the proposed method in multi-dimensional domains. Numerical experiments in three space dimensions show that the modeling error from truncating the domain can be made very small by choosing a sufficiently wide super-grid damping layer. The numerical accuracy is first evaluated against analytical solutions of Lamb’s problem, where fourth order accuracy is observed with a sixth order artificial dissipation. We then use successive grid refinements to study the numerical accuracy in the more complicated motion due to a point moment tensor source in a regularized layered material.« less

Authors:
;
Publication Date:
Research Org.:
Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
Sponsoring Org.:
USDOE
OSTI Identifier:
1224410
Report Number(s):
LLNL-JRNL-610212
Journal ID: ISSN 1815-2406
DOE Contract Number:  
AC52-07NA27344
Resource Type:
Journal Article
Resource Relation:
Journal Name: Communications in Computational Physics; Journal Volume: 16; Journal Issue: 04
Country of Publication:
United States
Language:
English
Subject:
58 GEOSCIENCES; 97 MATHEMATICS, COMPUTING, AND INFORMATION SCIENCE

Citation Formats

Petersson, N. Anders, and Sjögreen, Björn. Super-Grid Modeling of the Elastic Wave Equation in Semi-Bounded Domains. United States: N. p., 2014. Web. doi:10.4208/cicp.290113.220514a.
Petersson, N. Anders, & Sjögreen, Björn. Super-Grid Modeling of the Elastic Wave Equation in Semi-Bounded Domains. United States. doi:10.4208/cicp.290113.220514a.
Petersson, N. Anders, and Sjögreen, Björn. Wed . "Super-Grid Modeling of the Elastic Wave Equation in Semi-Bounded Domains". United States. doi:10.4208/cicp.290113.220514a. https://www.osti.gov/servlets/purl/1224410.
@article{osti_1224410,
title = {Super-Grid Modeling of the Elastic Wave Equation in Semi-Bounded Domains},
author = {Petersson, N. Anders and Sjögreen, Björn},
abstractNote = {Abstract We develop a super-grid modeling technique for solving the elastic wave equation in semi-bounded two- and three-dimensional spatial domains. In this method, waves are slowed down and dissipated in sponge layers near the far-field boundaries. Mathematically, this is equivalent to a coordinate mapping that transforms a very large physical domain to a significantly smaller computational domain, where the elastic wave equation is solved numerically on a regular grid. To damp out waves that become poorly resolved because of the coordinate mapping, a high order artificial dissipation operator is added in layers near the boundaries of the computational domain. We prove by energy estimates that the super-grid modeling leads to a stable numerical method with decreasing energy, which is valid for heterogeneous material properties and a free surface boundary condition on one side of the domain. Our spatial discretization is based on a fourth order accurate finite difference method, which satisfies the principle of summation by parts. We show that the discrete energy estimate holds also when a centered finite difference stencil is combined with homogeneous Dirichlet conditions at several ghost points outside of the far-field boundaries. Therefore, the coefficients in the finite difference stencils need only be boundary modified near the free surface. This allows for improved computational efficiency and significant simplifications of the implementation of the proposed method in multi-dimensional domains. Numerical experiments in three space dimensions show that the modeling error from truncating the domain can be made very small by choosing a sufficiently wide super-grid damping layer. The numerical accuracy is first evaluated against analytical solutions of Lamb’s problem, where fourth order accuracy is observed with a sixth order artificial dissipation. We then use successive grid refinements to study the numerical accuracy in the more complicated motion due to a point moment tensor source in a regularized layered material.},
doi = {10.4208/cicp.290113.220514a},
journal = {Communications in Computational Physics},
number = 04,
volume = 16,
place = {United States},
year = {Wed Oct 01 00:00:00 EDT 2014},
month = {Wed Oct 01 00:00:00 EDT 2014}
}