A Caltech Library Service

Three-dimensional dynamic rupture simulation with a high-order discontinuous Galerkin method on unstructured tetrahedral meshes

Pelties, Christian and de la Puente, Joseph and Ampuero, Jean-Paul and Brietzke, Gilbert B. and Käser, Martin (2012) Three-dimensional dynamic rupture simulation with a high-order discontinuous Galerkin method on unstructured tetrahedral meshes. Journal of Geophysical Research B, 117 (B2). Art. No. B02309. ISSN 0148-0227. doi:10.1029/2011JB008857.

PDF - Published Version
See Usage Policy.


Use this Persistent URL to link to this item:


Accurate and efficient numerical methods to simulate dynamic earthquake rupture and wave propagation in complex media and complex fault geometries are needed to address fundamental questions in earthquake dynamics, to integrate seismic and geodetic data into emerging approaches for dynamic source inversion, and to generate realistic physics-based earthquake scenarios for hazard assessment. Modeling of spontaneous earthquake rupture and seismic wave propagation by a high-order discontinuous Galerkin (DG) method combined with an arbitrarily high-order derivatives (ADER) time integration method was introduced in two dimensions by de la Puente et al. (2009). The ADER-DG method enables high accuracy in space and time and discretization by unstructured meshes. Here we extend this method to three-dimensional dynamic rupture problems. The high geometrical flexibility provided by the usage of tetrahedral elements and the lack of spurious mesh reflections in the ADER-DG method allows the refinement of the mesh close to the fault to model the rupture dynamics adequately while concentrating computational resources only where needed. Moreover, ADER-DG does not generate spurious high-frequency perturbations on the fault and hence does not require artificial Kelvin-Voigt damping. We verify our three-dimensional implementation by comparing results of the SCEC TPV3 test problem with two well-established numerical methods, finite differences, and spectral boundary integral. Furthermore, a convergence study is presented to demonstrate the systematic consistency of the method. To illustrate the capabilities of the high-order accurate ADER-DG scheme on unstructured meshes, we simulate an earthquake scenario, inspired by the 1992 Landers earthquake, that includes curved faults, fault branches, and surface topography.

Item Type:Article
Related URLs:
URLURL TypeDescription
Ampuero, Jean-Paul0000-0002-4827-7987
Additional Information:© 2012 American Geophysical Union. Received 8 September 2011; revised 30 December 2011; accepted 6 January 2012; published 18 February 2012. The authors thank the DFG (Deutsche Forschungsgemeinschaft), as the work was supported through the Emmy Noether-Programm (KA 2281/2-1). J.-P. A. was partially funded by NSF (grant EAR-0944288) and by the Southern California Earthquake Center (funded by NSF Cooperative Agreement EAR-0106924 and USGS Cooperative Agreement 02HQAG0008). The DFM data used for comparison were provided by Luis A. Dalguer and the SBIEM solutions where produced with the code of Eric M. Dunham (MDSBI: Multidimensional spectral boundary integral, version 3.9.10, 2008, available at http://pangea.stanford. edu/~edunham/codes/codes.html). Furthermore, we thank Luis A. Dalguer and Alan Schiemenz for very helpful and fruitful discussions. Cristóbal E. Castro gave valuable comments and advice on the solution of the Riemann problem and the parallelization. We also thank M. Mai for providing computational resources as many parallel tests, the convergence test, and the SCEC benchmark have been computed on the BlueGene/P Shaheen of the King Abdullah University of Science and Technology, Saudi Arabia. This paper is SCEC contribution 1526 and Caltech Seismological Lab contribution 10067. The reviews and comments by J.-P. Vilotte, S. M. Day, and the Associate Editor are appreciated and helped us to improve the manuscript.
Group:Seismological Laboratory
Funding AgencyGrant Number
Deutsche Forschungsgemeinschaft (DGF)KA 2281/2-1
Subject Keywords:computational seismology; dynamic rupture; earthquake physics; strong ground motion
Other Numbering System:
Other Numbering System NameOther Numbering System ID
Southern California Earthquake Center (SCEC)1526
Caltech Division of Geological and Planetary Sciences10067
Issue or Number:B2
Record Number:CaltechAUTHORS:20120328-074050760
Persistent URL:
Official Citation:Pelties, C., J. de la Puente, J.-P. Ampuero, G. B. Brietzke, and M. Käser (2012), Three-dimensional dynamic rupture simulation with a high-order discontinuous Galerkin method on unstructured tetrahedral meshes, J. Geophys. Res., 117, B02309, doi:10.1029/2011JB008857.
Usage Policy:No commercial reproduction, distribution, display or performance rights in this work are provided.
ID Code:29870
Deposited By: Ruth Sustaita
Deposited On:28 Mar 2012 15:07
Last Modified:09 Nov 2021 19:31

Repository Staff Only: item control page