An Arbitrary High Order Discontinuous Galerkin Method for Elastic Waves on Unstructured Meshes II: The ThreeDimensional Isotropic Case
Dumbser, M., and M. Käser (2006),
An Arbitrary High Order Discontinuous Galerkin Method for Elastic Waves on Unstructured Meshes II: The ThreeDimensional Isotropic Case,
Geophysical Journal International, 167(1), 319336, doi:10.1111/j.1365246X.2006.03120.x.
 Abstract
 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 wellestablished 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 socalled Generalized Riemann Problems (GRP) at the element interfaces. The GRP solution provides simultaneously a numerical flux function as well as a timeintegration method. The main idea is a Taylor expansion in time in which all time derivatives are replaced by space derivatives using the socalled CauchyKovalewski or LaxWendroff 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 RungeKutta 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 threedimensional 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 ADERDG method shows spectral convergence on tetrahedral meshes. An application of the new method to a wellacknowledged test case and comparisons with analytical and reference solutions, obtained by different wellestablished methods, confirm the performance of the proposed method. Therefore, the development of the highly accurate ADERDG approach for tetrahedral meshes provides a numerical technique to approach threedimensional wave propagation problems in complex geometry with unforeseen accuracy.
 BibTeX

@article{id593,
author = {M. Dumbser and M. K{\"a}ser},
journal = {Geophysical Journal International},
number = {1},
pages = {319336},
title = {{An Arbitrary High Order Discontinuous Galerkin Method for Elastic Waves on Unstructured Meshes II: The ThreeDimensional Isotropic Case}},
volume = {167},
year = {2006},
doi = {10.1111/j.1365246X.2006.03120.x},
}
 EndNote

%0 Journal Article
%A Dumbser, M.
%A Käser, M.
%D 2006
%N 1
%V 167
%J Geophysical Journal International
%P 319336
%T An Arbitrary High Order Discontinuous Galerkin Method for Elastic Waves on Unstructured Meshes II: The ThreeDimensional Isotropic Case