thumbnail

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

By:
, ,
DOI: 10.1002/jgrb.50217

Links

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.

Additional Publication Details

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