2022 Poroelastic Implementation in PyLith: Gateway to Multiphysics

By Robert L Walker

University at Buffalo



Published on


Thursday May 12, 2022 @ 2P PDT Poroelastic Implementation in PyLith: Gateway to Multiphysics Robert Walker, SUNY Buffalo 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 Terzaghi’s consolidation problem for the one-dimensional case, 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.