Multiphysics Modeling in PyLith: Poroelasticity

Robert L Walker,Matthew G Knepley,Brad T Aagaard,Charles A Williams
DOI: https://doi.org/10.1093/gji/ggad370
IF: 3.352
2023-09-26
Geophysical Journal International
Abstract:Summary PyLith, a community, open-source code (https://geodynamics.org/resources/pylith/) for modeling 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 one-dimensional consolidation problem as outlined by Terzaghi, Mandel’s problem for the two-dimensional case, and Cryer’s problem for the three-dimensional 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.
geochemistry & geophysics
What problem does this paper attempt to address?