Gamra: Simple meshing for complex earthquakes
- Creators
- Landry, Walter
- Barbot, Sylvain
Abstract
The static offsets caused by earthquakes are well described by elastostatic models with a discontinuity in the displacement along the fault. A traditional approach to model this discontinuity is to align the numerical mesh with the fault and solve the equations using finite elements. However, this distorted mesh can be difficult to generate and update. We present a new numerical method, inspired by the Immersed Interface Method (Leveque and Li, 1994), for solving the elastostatic equations with embedded discontinuities. This method has been carefully designed so that it can be used on parallel machines on an adapted finite difference grid. We have implemented this method in Gamra, a new code for earth modeling. We demonstrate the correctness of the method with analytic tests, and we demonstrate its practical performance by solving a realistic earthquake model to extremely high precision.
Additional Information
© 2016 The Authors. Published by Elsevier Ltd. This is an open access article under the CC BY license(http://creativecommons.org/licenses/by/4.0/). Received 5 October 2015, Revised 18 February 2016, Accepted 19 February 2016, Available online 23 February 2016. We thank two anonymous reviewers for their comments that improved this paper. This research used resources of the Argonne Leadership Computing Facility, which is a DOE Office of Science User Facility supported under Contract DE-AC02-06CH11357. This research was supported by the Gordon and Betty Moore Foundation, the National Research Foundation of Singapore under the NRF Fellowship scheme (National Research Fellow Award No. NRF-NRFF2013-04) and by the Earth Observatory of Singapore and the National Research Foundation and the Singapore Ministry of Education under the Research Centres of Excellence initiative. This is EOS publication 113.Errata
The authors would like to acknowledge two more important sources of support for this research. This research was supported by NSF Award EAR-0949446 and the Tectonics Observatory. The authors would like to apologise for any inconvenience caused.Attached Files
Published - 1-s2.0-S0098300416300498-main.pdf
Submitted - 1603.01263v1.pdf
Supplemental Material - mmc1.zip
Erratum - 1-s2.0-S0098300417312086-main.pdf
Files
Additional details
- Alternative title
- Gamra: Simple Meshes for Complex Earthquakes
- Eprint ID
- 67489
- Resolver ID
- CaltechAUTHORS:20160531-100116966
- DE-AC02-06CH11357
- Department of Energy (DOE)
- Gordon and Betty Moore Foundation
- NRF-NRFF2013-04
- National Research Foundation (Singapore)
- Earth Observatory of Singapore
- Ministry of Education (Singapore)
- EAR-0949446
- NSF
- Caltech Tectonics Observatory
- Created
-
2016-05-31Created from EPrint's datestamp field
- Updated
-
2023-06-01Created from EPrint's last_modified field
- Caltech groups
- Infrared Processing and Analysis Center (IPAC)
- Other Numbering System Name
- Earth Observatory of Singapore
- Other Numbering System Identifier
- 113