Dynamic Rupture Modeling on Unstructured Meshes Using a
Discontinuous Galerkin Method
J. de la Puente
Department of Earth and Environmental Sciences, Ludwig-Maximilians-University, Munich, Germany,
now at Institut de Ci`encies del Mar, CSIC, Barcelona, Spain.
Seismological Laboratory, California Institute of Technology, Pasadena, USA.
Department of Earth and Environmental Sciences, Ludwig-Maximilians-University, Munich, Germany.
Abstract. We introduce the application of an Arbitrary high-order DERivative (ADER) Dis-
continuous 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 using external automated tools. Discontinuous Galerkin
methods are well-suited for solving dynamic rupture problems in the velocity-stress formu-
lation as the variables are naturally discontinuous at the interface between two elements. There-
fore 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