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

Title: An efficient high-order multiscale finite element method for frequency-domain elastic wave modeling

Abstract

Solving the frequency-domain elastic wave equation in highly heterogeneous and complex media is computationally challenging. Conventional methods for solving the elastic wave Helmholtz equation usually lead to a large-dimensional linear system that is difficult to solve without specialized and sophisticated techniques. Based on the multiscale finite element theory, we develop a novel method to solve the frequency-domain elastic wave equation in complex media. The key feature of our method is employing high-order multiscale basis functions defined by solving local linear problems to achieve model reduction, which eventually leads to a linear system with significantly reduced dimensions. Solving this reduced linear system therefore results in obvious computational time reduction. We use three 2D examples to verify the accuracy and efficiency of our high-order multiscale finite element method for solving the Helmholtz equation in complex isotropic and anisotropic elastic media. The results show that our new method can approximate the fine-scale reference solution on the coarse mesh with high accuracy and significantly reduced computational time at the linear system solving stage.

Authors:
 [1]; ORCiD logo [2];  [3];  [1]
  1. Chinese Univ. of Hong Kong (China)
  2. Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
  3. Texas A & M Univ., College Station, TX (United States)
Publication Date:
Research Org.:
Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
Sponsoring Org.:
USDOE Office of Energy Efficiency and Renewable Energy (EERE). Geothermal Technologies Office (EE-4G)
OSTI Identifier:
1558964
Alternate Identifier(s):
OSTI ID: 1558962
Report Number(s):
LA-UR-18-27125; LA-UR-18-24382
Journal ID: ISSN 1420-0597
Grant/Contract Number:  
89233218CNA000001
Resource Type:
Accepted Manuscript
Journal Name:
Computational Geosciences
Additional Journal Information:
Journal Volume: 23; Journal Issue: 5; Journal ID: ISSN 1420-0597
Publisher:
Springer
Country of Publication:
United States
Language:
English
Subject:
58 GEOSCIENCES; 97 MATHEMATICS AND COMPUTING; Earth Sciences

Citation Formats

Fu, Shubin, Gao, Kai, Gibson Jr., Richard, and Chung, Eric T. An efficient high-order multiscale finite element method for frequency-domain elastic wave modeling. United States: N. p., 2019. Web. doi:10.1007/s10596-019-09865-0.
Fu, Shubin, Gao, Kai, Gibson Jr., Richard, & Chung, Eric T. An efficient high-order multiscale finite element method for frequency-domain elastic wave modeling. United States. doi:10.1007/s10596-019-09865-0.
Fu, Shubin, Gao, Kai, Gibson Jr., Richard, and Chung, Eric T. Thu . "An efficient high-order multiscale finite element method for frequency-domain elastic wave modeling". United States. doi:10.1007/s10596-019-09865-0. https://www.osti.gov/servlets/purl/1558964.
@article{osti_1558964,
title = {An efficient high-order multiscale finite element method for frequency-domain elastic wave modeling},
author = {Fu, Shubin and Gao, Kai and Gibson Jr., Richard and Chung, Eric T.},
abstractNote = {Solving the frequency-domain elastic wave equation in highly heterogeneous and complex media is computationally challenging. Conventional methods for solving the elastic wave Helmholtz equation usually lead to a large-dimensional linear system that is difficult to solve without specialized and sophisticated techniques. Based on the multiscale finite element theory, we develop a novel method to solve the frequency-domain elastic wave equation in complex media. The key feature of our method is employing high-order multiscale basis functions defined by solving local linear problems to achieve model reduction, which eventually leads to a linear system with significantly reduced dimensions. Solving this reduced linear system therefore results in obvious computational time reduction. We use three 2D examples to verify the accuracy and efficiency of our high-order multiscale finite element method for solving the Helmholtz equation in complex isotropic and anisotropic elastic media. The results show that our new method can approximate the fine-scale reference solution on the coarse mesh with high accuracy and significantly reduced computational time at the linear system solving stage.},
doi = {10.1007/s10596-019-09865-0},
journal = {Computational Geosciences},
number = 5,
volume = 23,
place = {United States},
year = {2019},
month = {8}
}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record

Save / Share:

Works referenced in this record:

Multiscale modeling of acoustic wave propagation in 2D media
journal, March 2014


Massively parallel structured direct solver for equations describing time-harmonic qP-polarized waves in TTI media
journal, May 2012


An exact Riemann solver for wave propagation in arbitrary anisotropic elastic media with fluid coupling
journal, February 2018

  • Zhan, Qiwei; Ren, Qiang; Zhuang, Mingwei
  • Computer Methods in Applied Mechanics and Engineering, Vol. 329
  • DOI: 10.1016/j.cma.2017.09.007

A Galerkin least-squares finite element method for the two-dimensional Helmholtz equation
journal, February 1995

  • Thompson, Lonny L.; Pinsky, Peter M.
  • International Journal for Numerical Methods in Engineering, Vol. 38, Issue 3
  • DOI: 10.1002/nme.1620380303

High-order finite difference methods for the Helmholtz equation
journal, September 1998


A perfectly matched layer for the absorption of electromagnetic waves
journal, October 1994


A Nearly Analytic Discrete Method for Acoustic and Elastic Wave Equations in Anisotropic Media
journal, April 2003

  • Yang, D.
  • Bulletin of the Seismological Society of America, Vol. 93, Issue 2
  • DOI: 10.1785/0120020125

A staggered discontinuous Galerkin method for the simulation of seismic waves with surface topography
journal, July 2015


A Novel Multigrid Based Preconditioner For Heterogeneous Helmholtz Problems
journal, January 2006

  • Erlangga, Y. A.; Oosterlee, C. W.; Vuik, C.
  • SIAM Journal on Scientific Computing, Vol. 27, Issue 4
  • DOI: 10.1137/040615195

Spectral element method for elastic and acoustic waves in frequency domain
journal, December 2016


Fourth‐order finite‐difference P-SV seismograms
journal, November 1988


P-SV wave propagation in heterogeneous media: Velocity‐stress finite‐difference method
journal, April 1986


Wavenumber Explicit Convergence Analysis for Galerkin Discretizations of the Helmholtz Equation
journal, January 2011

  • Melenk, J. M.; Sauter, S.
  • SIAM Journal on Numerical Analysis, Vol. 49, Issue 3
  • DOI: 10.1137/090776202

Compact finite difference schemes of sixth order for the Helmholtz equation
journal, June 2007


An overview of SuperLU: Algorithms, implementation, and user interface
journal, September 2005


Sweeping preconditioners for elastic wave propagation with spectral element methods
journal, February 2014

  • Tsuji, Paul; Poulson, Jack; Engquist, Björn
  • ESAIM: Mathematical Modelling and Numerical Analysis, Vol. 48, Issue 2
  • DOI: 10.1051/m2an/2013114

Sweeping Preconditioner for the Helmholtz Equation: Moving Perfectly Matched Layers
journal, April 2011

  • Engquist, Björn; Ying, Lexing
  • Multiscale Modeling & Simulation, Vol. 9, Issue 2
  • DOI: 10.1137/100804644

A numerical homogenization method for heterogeneous, anisotropic elastic media based on multiscale theory
journal, July 2015


Plane Wave Discontinuous Galerkin Methods for the 2D Helmholtz Equation: Analysis of the p -Version
journal, January 2011

  • Hiptmair, R.; Moiola, A.; Perugia, I.
  • SIAM Journal on Numerical Analysis, Vol. 49, Issue 1
  • DOI: 10.1137/090761057

2-D non-periodic homogenization to upscale elastic media for P-SV waves: 2-D non-periodic homogenization, PSV case
journal, June 2010


On 3D modeling of seismic wave propagation via a structured parallel multifrontal direct Helmholtz solver
journal, August 2011


Multiscale finite element methods for high-contrast problems using local spectral basis functions
journal, February 2011

  • Efendiev, Yalchin; Galvis, Juan; Wu, Xiao-Hui
  • Journal of Computational Physics, Vol. 230, Issue 4
  • DOI: 10.1016/j.jcp.2010.09.026

Learning dominant wave directions for plane wave methods for high-frequency Helmholtz equations
journal, May 2017

  • Fang, Jun; Qian, Jianliang; Zepeda-Núñez, Leonardo
  • Research in the Mathematical Sciences, Vol. 4, Issue 1
  • DOI: 10.1186/s40687-017-0098-9

Using Plane Waves as Base Functions for Solving Time Harmonic Equations with the Ultra Weak Variational Formulation
journal, June 2003


Improved accuracy for the Helmholtz equation in unbounded domains
journal, April 2004

  • Turkel, Eli; Farhat, Charbel; Hetmaniuk, Ulrich
  • International Journal for Numerical Methods in Engineering, Vol. 59, Issue 15
  • DOI: 10.1002/nme.882

An Efficient Multiscale Finite‐Element Method for Frequency‐Domain Seismic Wave Propagation
journal, February 2018

  • Gao, Kai; Fu, Shubin; Chung, Eric T.
  • Bulletin of the Seismological Society of America, Vol. 108, Issue 2
  • DOI: 10.1785/0120170268

Earthquake ground motion modeling on parallel computers
conference, January 1996

  • Bao, Hesheng; Bielak, Jacobo; Ghattas, Omar
  • Proceedings of the 1996 ACM/IEEE conference on Supercomputing (CDROM) - Supercomputing '96
  • DOI: 10.1145/369028.369053

Preasymptotic Error Analysis of CIP-FEM and FEM for Helmholtz Equation with High Wave Number. Part II: $hp$ Version
journal, January 2013

  • Zhu, Lingxue; Wu, Haijun
  • SIAM Journal on Numerical Analysis, Vol. 51, Issue 3
  • DOI: 10.1137/120874643

Wavenumber explicit convergence analysis for finite element discretizations of general wave propagation problems
journal, May 2019


A rapidly converging domain decomposition method for the Helmholtz equation
journal, May 2013


Sweeping preconditioner for the Helmholtz equation: Hierarchical matrix representation
journal, February 2011

  • Engquist, Björn; Ying, Lexing
  • Communications on Pure and Applied Mathematics, Vol. 64, Issue 5
  • DOI: 10.1002/cpa.20358

Seismic imaging of complex onshore structures by 2D elastic frequency-domain full-waveform inversion
journal, November 2009

  • Brossier, Romain; Operto, Stéphane; Virieux, Jean
  • GEOPHYSICS, Vol. 74, Issue 6
  • DOI: 10.1190/1.3215771

3D finite-difference frequency-domain modeling of visco-acoustic wave propagation using a massively parallel direct solver: A feasibility study
journal, September 2007

  • Operto, Stéphane; Virieux, Jean; Amestoy, Patrick
  • GEOPHYSICS, Vol. 72, Issue 5
  • DOI: 10.1190/1.2759835

Second order homogenization of the elastic wave equation for non-periodic layered media
journal, August 2007


A dispersion minimizing finite difference scheme and preconditioned solver for the 3D Helmholtz equation
journal, October 2012

  • Chen, Zhongying; Cheng, Dongsheng; Wu, Tingting
  • Journal of Computational Physics, Vol. 231, Issue 24
  • DOI: 10.1016/j.jcp.2012.07.048

Generalized Multiscale Finite Element Methods for Wave Propagation in Heterogeneous Media
journal, January 2014

  • Chung, Eric T.; Efendiev, Yalchin; Leung, Wing Tat
  • Multiscale Modeling & Simulation, Vol. 12, Issue 4
  • DOI: 10.1137/130926675

Compact 2D and 3D sixth order schemes for the Helmholtz equation with variable wave number
journal, January 2013

  • Turkel, Eli; Gordon, Dan; Gordon, Rachel
  • Journal of Computational Physics, Vol. 232, Issue 1
  • DOI: 10.1016/j.jcp.2012.08.016

On high order methods for the heterogeneous Helmholtz equation
journal, November 2016


Introduction to the spectral element method for three-dimensional seismic wave propagation
journal, December 1999


A Parallel Sweeping Preconditioner for Heterogeneous 3D Helmholtz Equations
journal, January 2013

  • Poulson, Jack; Engquist, Björn; Li, Siwei
  • SIAM Journal on Scientific Computing, Vol. 35, Issue 3
  • DOI: 10.1137/120871985

A generalized plane-wave numerical method for smooth nonconstant coefficients
journal, October 2013

  • Imbert-Gerard, L. -M.; Despres, B.
  • IMA Journal of Numerical Analysis, Vol. 34, Issue 3
  • DOI: 10.1093/imanum/drt030

A high-order discontinuous Galerkin method for wave propagation through coupled elastic–acoustic media
journal, December 2010

  • Wilcox, Lucas C.; Stadler, Georg; Burstedde, Carsten
  • Journal of Computational Physics, Vol. 229, Issue 24
  • DOI: 10.1016/j.jcp.2010.09.008

Parallel Computation of Entries of ${A}^{-1}$
journal, January 2015

  • Amestoy, Patrick R.; Duff, Iain S.; L'Excellent, Jean-Yves
  • SIAM Journal on Scientific Computing, Vol. 37, Issue 2
  • DOI: 10.1137/120902616

A discontinuous Galerkin method with Lagrange multipliers for the solution of Helmholtz problems in the mid-frequency regime
journal, March 2003

  • Farhat, Charbel; Harari, Isaac; Hetmaniuk, Ulrich
  • Computer Methods in Applied Mechanics and Engineering, Vol. 192, Issue 11-12
  • DOI: 10.1016/S0045-7825(02)00646-1

An overview of full-waveform inversion in exploration geophysics
journal, November 2009


2D and 3D frequency-domain elastic wave modeling in complex media with a parallel iterative solver
journal, May 2015


Generalized multiscale finite element methods (GMsFEM)
journal, October 2013


A high-order multiscale finite-element method for time-domain acoustic-wave modeling
journal, May 2018


An unsplit convolutional perfectly matched layer technique improved at grazing incidence for the viscoelastic wave equation
journal, October 2009


Modeling the propagation of elastic waves using a modified finite-difference grid
journal, January 2000


Elastic full-waveform inversion application using multicomponent measurements of seismic data collection
journal, March 2014


The mean transmission properties of a fault with imperfect facial contact
journal, June 1997


Multiscale Finite Element Methods for Nonlinear Problems and Their Applications
journal, January 2004

  • Efendiev, Y.; Ginting, V.; Hou, T. Y.
  • Communications in Mathematical Sciences, Vol. 2, Issue 4
  • DOI: 10.4310/CMS.2004.v2.n4.a2