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

Title: Optimal numerical parameterization of discontinuous Galerkin method applied to wave propagation problems

Abstract

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, we further show how to adequately choose elementary polynomial orders in order to equi-distribute a priori the discretization error. We also show a straightforward manner to allow variable polynomial orders in a DG scheme. We finally propose some numerical examples in the field of aero-acoustics.

Authors:
 [1];  [2];  [3];  [3];  [4]
  1. Universite catholique de Louvain, Department of Civil Engineering, Place du Levant 1, 1348 Louvain-la-Neuve (Belgium). E-mail: chevaugeon@gce.ucl.ac.be
  2. CENAERO CFD and Multiphysics Group, Batiment Mermoz 1, Av. J. Mermoz 30, b: 6041 Gosselies (Belgium)
  3. Free Field Technologies SA, Place de l'Universite, 1348 Louvain-la-Neuve (Belgium)
  4. Universite catholique de Louvain, Department of Civil Engineering, Place du Levant 1, 1348 Louvain-la-Neuve (Belgium) and Center for Systems Engineering and Applied Mechanics (CESAME), Universite catholique de Louvain, 1348 Louvain-la-Neuve (Belgium). E-mail: remacle@gce.ucl.ac.be
Publication Date:
OSTI Identifier:
20991571
Resource Type:
Journal Article
Resource Relation:
Journal Name: Journal of Computational Physics; Journal Volume: 223; Journal Issue: 1; Other Information: DOI: 10.1016/j.jcp.2006.09.005; PII: S0021-9991(06)00432-3; Copyright (c) 2006 Elsevier Science B.V., Amsterdam, The Netherlands, All rights reserved; Country of input: International Atomic Energy Agency (IAEA)
Country of Publication:
United States
Language:
English
Subject:
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS; ACCURACY; ACOUSTICS; CONVERGENCE; ERRORS; ONE-DIMENSIONAL CALCULATIONS; POLYNOMIALS; WAVE PROPAGATION

Citation Formats

Chevaugeon, Nicolas, Hillewaert, Koen, Gallez, Xavier, Ploumhans, Paul, and Remacle, Jean-Francois. Optimal numerical parameterization of discontinuous Galerkin method applied to wave propagation problems. United States: N. p., 2007. Web. doi:10.1016/j.jcp.2006.09.005.
Chevaugeon, Nicolas, Hillewaert, Koen, Gallez, Xavier, Ploumhans, Paul, & Remacle, Jean-Francois. Optimal numerical parameterization of discontinuous Galerkin method applied to wave propagation problems. United States. doi:10.1016/j.jcp.2006.09.005.
Chevaugeon, Nicolas, Hillewaert, Koen, Gallez, Xavier, Ploumhans, Paul, and Remacle, Jean-Francois. Tue . "Optimal numerical parameterization of discontinuous Galerkin method applied to wave propagation problems". United States. doi:10.1016/j.jcp.2006.09.005.
@article{osti_20991571,
title = {Optimal numerical parameterization of discontinuous Galerkin method applied to wave propagation problems},
author = {Chevaugeon, Nicolas and Hillewaert, Koen and Gallez, Xavier and Ploumhans, Paul and Remacle, Jean-Francois},
abstractNote = {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, we further show how to adequately choose elementary polynomial orders in order to equi-distribute a priori the discretization error. We also show a straightforward manner to allow variable polynomial orders in a DG scheme. We finally propose some numerical examples in the field of aero-acoustics.},
doi = {10.1016/j.jcp.2006.09.005},
journal = {Journal of Computational Physics},
number = 1,
volume = 223,
place = {United States},
year = {Tue Apr 10 00:00:00 EDT 2007},
month = {Tue Apr 10 00:00:00 EDT 2007}
}
  • 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 » 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.« less
  • We present an algorithm for the numerical simulation of seismic wave propagation in models with a complex near surface part and free surface topography. The approach is based on the combination of finite differences with the discontinuous Galerkin method. The discontinuous Galerkin method can be used on polyhedral meshes; thus, it is easy to handle the complex surfaces in the models. However, this approach is computationally intense in comparison with finite differences. Finite differences are computationally efficient, but in general, they require rectangular grids, leading to the stair-step approximation of the interfaces, which causes strong diffraction of the wavefield. Inmore » this research we present a hybrid algorithm where the discontinuous Galerkin method is used in a relatively small upper part of the model and finite differences are applied to the main part of the model.« less
  • Locally refined meshes impose severe stability constraints on explicit time-stepping methods for the numerical simulation of time dependent wave phenomena. To overcome that stability restriction, local time-stepping methods are developed, which allow arbitrarily small time steps precisely where small elements in the mesh are located. When combined with a discontinuous Galerkin finite element discretization in space, which inherently leads to a diagonal mass matrix, the resulting numerical schemes are fully explicit. Starting from the classical Adams-Bashforth multi-step methods, local time stepping schemes of arbitrarily high accuracy are derived. Numerical experiments validate the theory and illustrate the usefulness of the proposedmore » time integration schemes.« less
  • No abstract prepared.