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

doi:https://doi.org/10.1016/j.jcp.2010.09.008

Title: A high-order discontinuous Galerkin method for wave propagation through coupled elastic–acoustic media

Full Record
Other Related Research

Abstract

We introduce a high-order discontinuous Galerkin (dG) scheme for the numerical solution of three-dimensional (3D) wave propagation problems in coupled elastic-acoustic media. A velocity-strain formulation is used, which allows for the solution of the acoustic and elastic wave equations within the same unified framework. Careful attention is directed at the derivation of a numerical flux that preserves high-order accuracy in the presence of material discontinuities, including elastic-acoustic interfaces. Explicit expressions for the 3D upwind numerical flux, derived as an exact solution for the relevant Riemann problem, are provided. The method supports h-non-conforming meshes, which are particularly effective at allowing local adaptation of the mesh size to resolve strong contrasts in the local wavelength, as well as dynamic adaptivity to track solution features. The use of high-order elements controls numerical dispersion, enabling propagation over many wave periods. We prove consistency and stability of the proposed dG scheme. To study the numerical accuracy and convergence of the proposed method, we compare against analytical solutions for wave propagation problems with interfaces, including Rayleigh, Lamb, Scholte, and Stoneley waves as well as plane waves impinging on an elastic-acoustic interface. Spectral rates of convergence are demonstrated for these problems, which include a non-conforming mesh case.more »« less

Authors:: Wilcox, Lucas C.; Stadler, Georg; Burstedde, Carsten; Ghattas, Omar

Publication Date:: 2010-12-01

Research Org.:: Oak Ridge National Laboratory (ORNL), Oak Ridge, TN (United States). Oak Ridge Leadership Computing Facility (OLCF), Oak Ridge, TN (United States); Univ. of Texas, Austin, TX (United States)

Sponsoring Org.:: USDOE Office of Science (SC)

OSTI Identifier:: 1564731

DOE Contract Number:: FC02-06ER25782

Resource Type:: Journal Article

Journal Name:: Journal of Computational Physics

Additional Journal Information:: Journal Volume: 229; Journal Issue: 24; Journal ID: ISSN 0021-9991

Publisher:: Elsevier

Country of Publication:: United States

Language:: English

Subject:: Computer Science; Physics

Citation Formats


                    Wilcox, Lucas C., Stadler, Georg, Burstedde, Carsten, and Ghattas, Omar. A high-order discontinuous Galerkin method for wave propagation through coupled elastic–acoustic media.  United States: N. p., 2010. 
        Web.  doi:10.1016/j.jcp.2010.09.008.

Copy to clipboard


                    Wilcox, Lucas C., Stadler, Georg, Burstedde, Carsten, & Ghattas, Omar. A high-order discontinuous Galerkin method for wave propagation through coupled elastic–acoustic media.  United States.  https://doi.org/10.1016/j.jcp.2010.09.008

Copy to clipboard


                    Wilcox, Lucas C., Stadler, Georg, Burstedde, Carsten, and Ghattas, Omar. Wed .  
        "A high-order discontinuous Galerkin method for wave propagation through coupled elastic–acoustic media".  United States.  https://doi.org/10.1016/j.jcp.2010.09.008.

Copy to clipboard


                    
@article{osti_1564731,

  title        = {A high-order discontinuous Galerkin method for wave propagation through coupled elastic–acoustic media},

  author       = {Wilcox, Lucas C. and Stadler, Georg and Burstedde, Carsten and Ghattas, Omar},

  abstractNote = {We introduce a high-order discontinuous Galerkin (dG) scheme for the numerical solution of three-dimensional (3D) wave propagation problems in coupled elastic-acoustic media. A velocity-strain formulation is used, which allows for the solution of the acoustic and elastic wave equations within the same unified framework. Careful attention is directed at the derivation of a numerical flux that preserves high-order accuracy in the presence of material discontinuities, including elastic-acoustic interfaces. Explicit expressions for the 3D upwind numerical flux, derived as an exact solution for the relevant Riemann problem, are provided. The method supports h-non-conforming meshes, which are particularly effective at allowing local adaptation of the mesh size to resolve strong contrasts in the local wavelength, as well as dynamic adaptivity to track solution features. The use of high-order elements controls numerical dispersion, enabling propagation over many wave periods. We prove consistency and stability of the proposed dG scheme. To study the numerical accuracy and convergence of the proposed method, we compare against analytical solutions for wave propagation problems with interfaces, including Rayleigh, Lamb, Scholte, and Stoneley waves as well as plane waves impinging on an elastic-acoustic interface. Spectral rates of convergence are demonstrated for these problems, which include a non-conforming mesh case. Finally, we present scalability results for a parallel implementation of the proposed high-order dG scheme for large-scale seismic wave propagation in a simplified earth model, demonstrating high parallel efficiency for strong scaling to the full size of the Jaguar Cray XT5 supercomputer.},

  doi          = {10.1016/j.jcp.2010.09.008},

  url          = {https://www.osti.gov/biblio/1564731},
  journal      = {Journal of Computational Physics},

  issn         = {0021-9991},

  number       = 24,

  volume       = 229,

  place        = {United States},

  year         = {2010},

  month        = {12}

}

Copy to clipboard

Journal Article:

https://doi.org/10.1016/j.jcp.2010.09.008

Other availability

Find in Google Scholar

Search WorldCat to find libraries that may hold this journal

Save / Share:

Export Metadata

Save to My Library

Similar records in OSTI.GOV collections:

A high-order discontinuous Galerkin method for wave propagation through coupled elastic-acoustic media

Journal Article Wilcox, Lucas C., E-mail: lucas@hypercomp.ne ; Stadler, Georg ; Burstedde, Carsten ; ... - Journal of Computational Physics

We introduce a high-order discontinuous Galerkin (dG) scheme for the numerical solution of three-dimensional (3D) wave propagation problems in coupled elastic-acoustic media. A velocity-strain formulation is used, which allows for the solution of the acoustic and elastic wave equations within the same unified framework. Careful attention is directed at the derivation of a numerical flux that preserves high-order accuracy in the presence of material discontinuities, including elastic-acoustic interfaces. Explicit expressions for the 3D upwind numerical flux, derived as an exact solution for the relevant Riemann problem, are provided. The method supports h-non-conforming meshes, which are particularly effective at allowing localmore »« less
https://doi.org/10.1016/j.jcp.2010.09.008
Optimal numerical parameterization of discontinuous Galerkin method applied to wave propagation problems

Journal Article Chevaugeon, Nicolas ; Hillewaert, Koen ; Gallez, Xavier ; ... - Journal of Computational Physics

This paper deals with the high-order discontinuous Galerkin (DG) method for solving wave propagation problems. First, we develop a one-dimensional DG scheme and numerically compute dissipation and dispersion errors for various polynomial orders. An optimal combination of time stepping scheme together with the high-order DG spatial scheme is presented. It is shown that using a time stepping scheme with the same formal accuracy as the DG scheme is too expensive for the range of wave numbers that is relevant for practical applications. An efficient implementation of a high-order DG method in three dimensions is presented. Using 1D convergence results, wemore »« less
https://doi.org/10.1016/j.jcp.2006.09.005
Dynamic earthquake rupture simulations on nonplanar faults embedded in 3D geometrically complex, heterogeneous elastic solids

Journal Article Duru, Kenneth ; Dunham, Eric M. ; Institute for Computational and Mathematical Engineering, Stanford University, Stanford, CA - Journal of Computational Physics

Dynamic propagation of shear ruptures on a frictional interface in an elastic solid is a useful idealization of natural earthquakes. The conditions relating discontinuities in particle velocities across fault zones and tractions acting on the fault are often expressed as nonlinear friction laws. The corresponding initial boundary value problems are both numerically and computationally challenging. In addition, seismic waves generated by earthquake ruptures must be propagated for many wavelengths away from the fault. Therefore, reliable and efficient numerical simulations require both provably stable and high order accurate numerical methods. We present a high order accurate finite difference method for: a)more »« less
https://doi.org/10.1016/J.JCP.2015.10.021
Development of a Novel Accelerator for Neutron Transport Solution Using the Galerkin Spectral Element Methods (Final Report)

Technical Report Wang, Dean ; Downar, Thomas ; Xu, Yunlin ; ...

This report summarizes the results of a three-year research project sponsored by the U.S. Department of Energy (DOE) Nuclear Energy University Program (NEUP) to develop and implement advanced acceleration schemes for the DOE NEAMS neutronics code PROTEUSSN. The project team included the University of Massachusetts Lowell, The Ohio State University, University of Michigan, and Argonne National Laboratory. The PROTEUS code is a high-fidelity capable deterministic neutron transport code based on unstructured finite element meshes, which solves the steady-state and transient neutron transport problem using the 2nd-order discrete ordinate method (SN2ND), the method of characteristics (MOC), and the advanced nodal transportmore »« less
https://doi.org/10.2172/1511575

Full Text Available
High order ADER schemes for a unified first order hyperbolic formulation of continuum mechanics: Viscous heat-conducting fluids and elastic solids

Journal Article Dumbser, Michael ; Peshkov, Ilya ; Romenski, Evgeniy ; ... - Journal of Computational Physics

Highlights: • High order schemes for a unified first order hyperbolic formulation of continuum mechanics. • The mathematical model applies simultaneously to fluid mechanics and solid mechanics. • Viscous fluids are treated in the frame of hyper-elasticity as generalized visco-plastic solids. • Formal asymptotic analysis reveals the connection with the Navier–Stokes equations. • The distortion tensor A in the model appears to be well-suited for flow visualization. - Abstract: This paper is concerned with the numerical solution of the unified first order hyperbolic formulation of continuum mechanics recently proposed by Peshkov and Romenski [110], further denoted as HPR model. Inmore »« less
https://doi.org/10.1016/J.JCP.2016.02.015

Similar Records