An arbitrary high-order discontinuous Galerkin method for elastic waves on unstructured meshes -: I.: The two-dimensional isotropic case with external source terms

被引:308
作者
Kaser, Martin [1 ]
Dumbser, Michael [1 ]
机构
[1] Univ Trent, Dept Civil & Environm Engn, I-38100 Trento, Italy
关键词
arbitrary high order; discontinuous Galerkin method; elastic waves; unstructured mesh;
D O I
10.1111/j.1365-246X.2006.03051.x
中图分类号
P3 [地球物理学]; P59 [地球化学];
学科分类号
0708 ; 070902 ;
摘要
We present a new numerical approach to solve the elastic wave equation in heterogeneous media in the presence of externally given source terms with arbitrary high-order accuracy in space and time on unstructured triangular meshes. We combine a discontinuous Galerkin (DG) method with the ideas of the ADER time integration approach using Arbitrary high-order DERivatives. The time integration is performed via the so-called Cauchy-Kovalewski procedure using repeatedly the governing partial differential equation itself. In contrast to classical finite element methods we allow for discontinuities of the piecewise polynomial approximation of the solution at element interfaces. This way, we can use the well-established theory of fluxes across element interfaces based on the solution of Riemann problems as developed in the finite volume framework. In particular, we replace time derivatives in the Taylor expansion of the time integration procedure by space derivatives to obtain a numerical scheme of the same high order in space and time using only one single explicit step to evolve the solution from one time level to another. The method is specially suited for linear hyperbolic systems such as the heterogeneous elastic wave equations and allows an efficient implementation. We consider continuous sources in space and time and point sources characterized by a Delta distribution in space and some continuous source time function. Hereby, the method is able to deal with point sources at any position in the computational domain that does not necessarily need to coincide with a mesh point. Interpolation is automatically performed by evaluation of test functions at the source locations. The convergence analysis demonstrates that very high accuracy is retained even on strongly irregular meshes and by increasing the order of the ADER-DG schemes computational time and storage space can be reduced remarkably. Applications of the proposed method to Lamb's Problem, a problem of strong material heterogeneities and to an example of global seismic wave propagation finally confirm its accuracy, robustness and high flexibility.
引用
收藏
页码:855 / 877
页数:23
相关论文
共 58 条
[1]  
Aki K., 2002, QUANTITATIVE SEISMOL
[2]  
[Anonymous], 59 ANN INT M SOC EXP
[3]  
[Anonymous], FINITE ELEMENT METHO
[4]   Quadrature-free implementation of discontinuous Galerkin method for hyperbolic equations [J].
Atkins, HL ;
Shu, CW .
AIAA JOURNAL, 1998, 36 (05) :775-782
[5]  
Bedford A., 1994, INTRO ELASTIC WAVE P
[6]   A PERFECTLY MATCHED LAYER FOR THE ABSORPTION OF ELECTROMAGNETIC-WAVES [J].
BERENGER, JP .
JOURNAL OF COMPUTATIONAL PHYSICS, 1994, 114 (02) :185-200
[7]  
Berg P., 1994, MODELING EARTH OIL E, P421
[8]  
Butcher J. C., 1987, The Numerical Analysis of Ordinary Differential Equations: Runge-Kutta and General Linear Methods
[9]   THE WAVE-EQUATION IN GENERALIZED COORDINATES [J].
CARCIONE, JM .
GEOPHYSICS, 1994, 59 (12) :1911-1919
[10]   THE RUNGE-KUTTA LOCAL PROJECTION RHO-1-DISCONTINUOUS-GALERKIN FINITE-ELEMENT METHOD FOR SCALAR CONSERVATION-LAWS [J].
COCKBURN, B ;
SHU, CW .
ESAIM-MATHEMATICAL MODELLING AND NUMERICAL ANALYSIS-MODELISATION MATHEMATIQUE ET ANALYSE NUMERIQUE, 1991, 25 (03) :337-361