Change search
ReferencesLink to record
Permanent link

Direct link
Simulation of Dynamic Earthquake Ruptures in Complex Geometries Using High-Order Finite Difference Methods
Department of Geophysics, Stanford University, USA.
Department of Geophysics, Stanford University, USA.
Linköping University, Department of Mathematics, Computational Mathematics. Linköping University, The Institute of Technology.ORCID iD: 0000-0002-7972-6183
2013 (English)In: Journal of Scientific Computing, ISSN 0885-7474, E-ISSN 1573-7691, Vol. 55, no 1, 92-124 p.Article in journal (Refereed) Published
Abstract [en]

We develop a stable and high-order accurate finite difference method for problems in earthquake rupture dynamics in complex geometries with multiple faults. The bulk material is an isotropic elastic solid cut by pre-existing fault interfaces that accommodate relative motion of the material on the two sides. The fields across the interfaces are related through friction laws which depend on the sliding velocity, tractions acting on the interface, and state variables which evolve according to ordinary differential equations involving local fields. The method is based on summation-by-parts finite difference operators with irregular geometries handled through coordinate transforms and multi-block meshes. Boundary conditions as well as block interface conditions (whether frictional or otherwise) are enforced weakly through the simultaneous approximation term method, resulting in a provably stable discretization. The theoretical accuracy and stability results are confirmed with the method of manufactured solutions. The practical benefits of the new methodology are illustrated in a simulation of a subduction zone megathrust earthquake, a challenging application problem involving complex free-surface topography, nonplanar faults, and varying material properties.

Place, publisher, year, edition, pages
Springer, 2013. Vol. 55, no 1, 92-124 p.
Keyword [en]
High-order finite difference – Nonlinear boundary conditions – Simultaneous approximation term method – Elastodynamics – Summation-by-parts – Friction – Wave propagation – Multi-block – Coordinate transforms – Weak boundary conditions
National Category
Computational Mathematics
URN: urn:nbn:se:liu:diva-80992DOI: 10.1007/s10915-012-9624-5ISI: 000316148400005OAI: diva2:549633
Available from: 2012-09-05 Created: 2012-09-05 Last updated: 2013-08-30

Open Access in DiVA

fulltext(2207 kB)289 downloads
File information
File name FULLTEXT01.pdfFile size 2207 kBChecksum SHA-512
Type fulltextMimetype application/pdf

Other links

Publisher's full text

Search in DiVA

By author/editor
Nordström, Jan
By organisation
Computational MathematicsThe Institute of Technology
In the same journal
Journal of Scientific Computing
Computational Mathematics

Search outside of DiVA

GoogleGoogle Scholar
Total: 289 downloads
The number of downloads is the sum of all downloads of full texts. It may include eg previous versions that are now no longer available

Altmetric score

Total: 210 hits
ReferencesLink to record
Permanent link

Direct link