A domain decomposition approach to implementing fault slip in finite-element models of quasi-static and dynamic crustal deformation
Journal of Geophysical Research B: Solid Earth
Links
- More information: Publisher Index Page (via DOI)
- Open Access Version: External Repository
- Download citation as: RIS | Dublin Core
Abstract
We employ a domain decomposition approach with Lagrange multipliers to implement fault slip in a finite-element code, PyLith, for use in both quasi-static and dynamic crustal deformation applications. This integrated approach to solving both quasi-static and dynamic simulations leverages common finite-element data structures and implementations of various boundary conditions, discretization schemes, and bulk and fault rheologies. We have developed a custom preconditioner for the Lagrange multiplier portion of the system of equations that provides excellent scalability with problem size compared to conventional additive Schwarz methods. We demonstrate application of this approach using benchmarks for both quasi-static viscoelastic deformation and dynamic spontaneous rupture propagation that verify the numerical implementation in PyLith.
| Publication type | Article |
|---|---|
| Publication Subtype | Journal Article |
| Title | A domain decomposition approach to implementing fault slip in finite-element models of quasi-static and dynamic crustal deformation |
| Series title | Journal of Geophysical Research B: Solid Earth |
| DOI | 10.1002/jgrb.50217 |
| Volume | 118 |
| Issue | 6 |
| Year Published | 2013 |
| Language | English |
| Publisher | American Geophysical Union |
| Contributing office(s) | Earthquake Science Center |
| Description | 21 p. |
| Larger Work Type | Article |
| Larger Work Subtype | Journal Article |
| Larger Work Title | Journal of Geophysical Research B: Solid Earth |
| First page | 3059 |
| Last page | 3079 |
| Google Analytic Metrics | Metrics page |