*** TOMORROW - MAY 12 @2 P PDT ***

Poroelastic Implementation in PyLith: Gateway to Multiphysics

Robert Walker, SUNY Buffalo

PyLith, a community, open-source code ( 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.