A Caltech Library Service

Dynamic rupture modeling on unstructured meshes using a discontinuous Galerkin method

de la Puente, J. and Ampuero, J.-P. and Käser, M. (2009) Dynamic rupture modeling on unstructured meshes using a discontinuous Galerkin method. Journal of Geophysical Research B, 114 (B10). Art. No. B10302. ISSN 0148-0227. doi:10.1029/2008JB006271.

PDF - Published Version
See Usage Policy.


Use this Persistent URL to link to this item:


We introduce the application of an arbitrary high-order derivative (ADER) discontinuous Galerkin (DG) method to simulate earthquake rupture dynamics. The ADER-DG method uses triangles as computational cells which simplifies the process of discretization of very complex surfaces and volumes by using external automated tools. Discontinuous Galerkin methods are well suited for solving dynamic rupture problems in the velocity-stress formulation as the variables are naturally discontinuous at the interface between two elements. Therefore, the fault has to be honored by the computational mesh. The so-called Riemann problem can be solved to obtain well defined values of the variables at the discontinuity itself. Fault geometries of high complexity can be modeled thanks to the flexibility of unstructured meshes, which solves a major bottleneck of other high-order numerical methods. Additionally, element refinement and coarsening are easily controlled in the meshing process to better resolve the near-fault area of the model. The fundamental properties of the method are shown, as well as a series of validating exercises with reference solutions and a comparison with the well-established finite difference, boundary integral, and spectral element methods, in order to test the accuracy of our formulation. An example of dynamic rupture on a nonplanar fault based upon the Landers 1992 earthquake fault system is presented to illustrate the main potentials of the new method.

Item Type:Article
Related URLs:
URLURL TypeDescription
Ampuero, J.-P.0000-0002-4827-7987
Additional Information:©2009 American Geophysical Union. Received 23 December 2008; accepted 6 July 2009; published 1 October 2009. The authors thank the Deutsche Forschungsgemeinschaft (DFG), as the work was supported through the Emmy Noether-program (KA 2281/2-1). We thank Gilbert Brietzke for providing us with the FD and SBIE solutions shown for the 2-D analog of the SCEC test case in section 6 and Cristóbal Castro for some fruitful discussions on the treatment of nonlinear Riemann problems. Steve M. Day and Peter Moczo are also acknowledged for their useful reviews which have helped to clarify and improve the original manuscript.
Group:Seismological Laboratory
Funding AgencyGrant Number
Deutsche Forschungsgemeinschaft (DFG)KA 2281/2-1
Subject Keywords:fault, rupture, modeling
Issue or Number:B10
Classification Code:7209 Seismology: Earthquake dynamics (1242); 7290 Seismology: Computational seismology; 7260 Seismology: Theory
Record Number:CaltechAUTHORS:20091022-113110359
Persistent URL:
Official Citation:de la Puente, J., J.-P. Ampuero, and M. Käser (2009), Dynamic rupture modeling on unstructured meshes using a discontinuous Galerkin method, J. Geophys. Res., 114, B10302, doi:10.1029/2008JB006271.
Usage Policy:No commercial reproduction, distribution, display or performance rights in this work are provided.
ID Code:16447
Deposited By: George Porter
Deposited On:22 Oct 2009 21:30
Last Modified:08 Nov 2021 23:27

Repository Staff Only: item control page