Multiphysics modelling in PyLith: Poroelasticity
Links
- More information: Publisher Index Page (via DOI)
- Download citation as: RIS | Dublin Core
Abstract
PyLith, a community, open-source code for modelling quasi-static and dynamic crustal deformation with an emphasis on earthquake faulting, has recently been updated with a flexible multiphysics implementation. We demonstrate the versatility of the multiphysics implementation by extending the code to model fully coupled continuum poromechanics. We verify the newly incorporated physics using standard benchmarks for a porous medium saturated with a slightly compressible fluid. The benchmarks include the 1-D consolidation problem as outlined by Terzaghi, Mandel’s problem for the 2-D case, and Cryer’s problem for the 3-D case. All three benchmarks have been added to the PyLith continuous integration test suite. We compare the closed form analytical solution for each benchmark against solutions generated by our updated code, and lastly, demonstrate that the poroelastic material formulation may be used alongside the existing fault implementation in PyLith.
| Publication type | Article |
|---|---|
| Publication Subtype | Journal Article |
| Title | Multiphysics modelling in PyLith: Poroelasticity |
| Series title | Geophysical Journal International |
| DOI | 10.1093/gji/ggad370 |
| Volume | 235 |
| Issue | 3 |
| Year Published | 2023 |
| Language | English |
| Publisher | Oxford Academic |
| Contributing office(s) | Geologic Hazards Science Center - Seismology / Geomagnetism |
| Description | 34 p. |
| First page | 2442 |
| Last page | 2475 |
| Google Analytic Metrics | Metrics page |