We present a new method for near-source ground-motion calculations due to earthquake rupture on potentially geometrically complex faults. Following the recently introduced Discontinuous Galerkin approach with local time stepping on tetrahedral meshes, we use piecewise polynomial approximations of the unknown variables inside each element and achieve the same approximation order in time and space due to the new ADER time integration scheme that uses Arbitrary high-order DERivatives. We show how an external source term and its heterogeneous properties in space and time, given by a fine discretization of an extended rupture surface, can be included in much coarser tetrahedral meshes due to the subcell resolution of the high-order polynomial representation. Hereby, the rupture surface is represented very generally as a point cloud of the center locations of individual subfaults at which each polynomial test function is evaluated exactly inside an element and the space–time integration of the source term is accurately computed at each timestep. Besides the incorporation of complex source kinematics we also present the effects of model boundaries that can degrade the accuracy of seismograms due to weak artificial reflections. We propose an extended computational domain of a coarsely meshed buffer region and show that our scheme using the local timestepping completely avoids such boundary problems with only slightly increasing the computational cost. We validate the new approach against different test cases, comparing our results with analytic, quasianalytic, and a series of reference solutions. Our work shows that adding the functionality of accurately treating finite source-rupture models into the general framework of the ADER-Discontinuous Galerkin approach is an important contribution to modeling realistic earthquake scenarios, allowing the efficient inclusion of heterogeneous source kinematics and complex rupture-surface geometries in near-source ground-motion simulations.

Accurate Calculation of Fault-Rupture Models Using the High-Order Discontinuous Galerkin Method on Tetrahedral Meshes / M., Käser; P., Mai; Dumbser, Michael. - In: BULLETIN OF THE SEISMOLOGICAL SOCIETY OF AMERICA. - ISSN 0037-1106. - STAMPA. - v. 97:i. 5(2007), pp. 1570-1586. [10.1785/0120060253]

Accurate Calculation of Fault-Rupture Models Using the High-Order Discontinuous Galerkin Method on Tetrahedral Meshes

Dumbser, Michael
2007-01-01

Abstract

We present a new method for near-source ground-motion calculations due to earthquake rupture on potentially geometrically complex faults. Following the recently introduced Discontinuous Galerkin approach with local time stepping on tetrahedral meshes, we use piecewise polynomial approximations of the unknown variables inside each element and achieve the same approximation order in time and space due to the new ADER time integration scheme that uses Arbitrary high-order DERivatives. We show how an external source term and its heterogeneous properties in space and time, given by a fine discretization of an extended rupture surface, can be included in much coarser tetrahedral meshes due to the subcell resolution of the high-order polynomial representation. Hereby, the rupture surface is represented very generally as a point cloud of the center locations of individual subfaults at which each polynomial test function is evaluated exactly inside an element and the space–time integration of the source term is accurately computed at each timestep. Besides the incorporation of complex source kinematics we also present the effects of model boundaries that can degrade the accuracy of seismograms due to weak artificial reflections. We propose an extended computational domain of a coarsely meshed buffer region and show that our scheme using the local timestepping completely avoids such boundary problems with only slightly increasing the computational cost. We validate the new approach against different test cases, comparing our results with analytic, quasianalytic, and a series of reference solutions. Our work shows that adding the functionality of accurately treating finite source-rupture models into the general framework of the ADER-Discontinuous Galerkin approach is an important contribution to modeling realistic earthquake scenarios, allowing the efficient inclusion of heterogeneous source kinematics and complex rupture-surface geometries in near-source ground-motion simulations.
2007
i. 5
M., Käser; P., Mai; Dumbser, Michael
Accurate Calculation of Fault-Rupture Models Using the High-Order Discontinuous Galerkin Method on Tetrahedral Meshes / M., Käser; P., Mai; Dumbser, Michael. - In: BULLETIN OF THE SEISMOLOGICAL SOCIETY OF AMERICA. - ISSN 0037-1106. - STAMPA. - v. 97:i. 5(2007), pp. 1570-1586. [10.1785/0120060253]
File in questo prodotto:
File Dimensione Formato  
KaeserMaiDumbser.pdf

Solo gestori archivio

Tipologia: Post-print referato (Refereed author’s manuscript)
Licenza: Tutti i diritti riservati (All rights reserved)
Dimensione 504.75 kB
Formato Adobe PDF
504.75 kB Adobe PDF   Visualizza/Apri

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/11572/69697
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 14
  • ???jsp.display-item.citation.isi??? 14
  • OpenAlex ND
social impact