An Arbitrary High Order Discontinuous Galerkin Method for Elastic Waves on Unstructured Meshes II: The Three-Dimensional Isotropic Case

We present a new numerical method to solve the heterogeneous elastic wave equations formulated as a linear hyperbolic system using first order derivatives with arbitrary high order accuracy in space and time on three dimensional unstructured tetrahedral meshes. The method combines the Discontinuous Galerkin (DG) Finite Element method with the ADER approach using Arbitrary high order DERivatives for flux calculation. In the DG framework, in contrast to classical Finite Element methods, the numerical solution is approximated by piecewise polynomials which allow for discontinuities at element interfaces. Therefore, the well-established theory of numerical fluxes across element interfaces obtained by the solution of Riemann problems can be applied as in the finite volume framework. To define a suitable flux over the element surfaces, we solve so-called Generalized Riemann Problems (GRP) at the element interfaces. The GRP solution provides simultaneously a numerical flux function as well as a time-integration method. The main idea is a Taylor expansion in time in which all time derivatives are replaced by space derivatives using the so-called Cauchy-Kovalewski or Lax-Wendroff procedure which makes extensive use of the governing PDE. The numerical solution can thus be advanced for one time step without intermediate stages as typical e.g. for classical Runge-Kutta time stepping schemes. Due to the ADER time integration technique the same approximation order in space and time is achieved automatically. Furthermore, the projection of the tetrahedral elements in physical space onto a canonical reference tetrahedron allows for an efficient implementation, as many computations of three-dimensional integrals can be carried out analytically beforehand. Based on a numerical convergence analysis we demonstrate that the new schemes provide very high order accuracy even on unstructured tetrahedral meshes and computational cost and storage space for a desired accuracy can be reduced by higher order schemes. Moreover, due to the choice of the basis functions for the piecewise polynomial approximation, the new ADER-DG method shows spectral convergence on tetrahedral meshes. An application of the new method to a well-acknowledged test case and comparisons with analytical and reference solutions, obtained by different well-established methods, confirm the performance of the proposed method. Therefore, the development of the highly accurate ADER-DG approach for tetrahedral meshes provides a numerical technique to approach three-dimensional wave propagation problems in complex geometry with unforeseen accuracy.

@article{id593,
  author = {Dumbser, M. and K\"aser, M.},
  doi = {10.1111/j.1365-246X.2006.03120.x},
  journal = {Geophysical Journal International},
  language = {en},
  number = {1},
  pages = {319-336},
  title = {An Arbitrary High Order Discontinuous Galerkin Method for Elastic Waves on Unstructured Meshes II: The Three-Dimensional Isotropic Case},
  volume = {167},
  year = {2006},
}

%O Journal Article
%A Dumbser, M.
%A Käser, M.
%R 10.1111/j.1365-246X.2006.03120.x
%J Geophysical Journal International
%G en
%N 1
%P 319-336
%T An Arbitrary High Order Discontinuous Galerkin Method for Elastic Waves on Unstructured Meshes II: The Three-Dimensional Isotropic Case
%V 167
%D 2006