Development of New Features for the PyLith Finite Element Code to Enable Innovative Research In Multiple Earthquake Cycle Simulations

来自这笔赠款的资金将支持一项为期两年的努力，以开发对现有的PyLith版本1.0地壳动力学建模包的代码增强。PyLith目前通过NSF/EAR支持的地球动力学计算基础设施(CIG)进行维护和支持，但这种支持解决了CIG目前没有针对未来开发的代码增强功能。改进措施将包括：1)耦合准静态和动态解，以实现对地震周期的完整模拟，包括动态破裂过程；2)开发非线性求解器，以模拟真实的地壳物质行为(例如，粘弹性、弹塑性、弹粘塑性)；3)重力预应力的自动计算；以及4)开发使用球面网格有限元框架模拟大空间区域动态的能力。这项工作将包括与CIG资助的人员(Matt Knepley，Argonne National Lab)、美国地质调查局(Brad Aagaard)和加州理工大学的一名博士生进行免费合作。代码开发将是开放源码的，并在cig存储库中维护。

Funding from this grant will support a two-year effort to develop code enhancements to the existing PyLith version 1.0 crustal dynamics modeling package. PyLith is currently maintained and supported through the NSF/EAR supported Computational Infrastructure for Geodynamics (CIG) facility but this support addresses code enhancements not currently targeted by CIG for future development. Enhancements will include: 1) coupling of quasi-static and dynamic solutions to allow for complete simulation of the earthquake cycle, including the dynamic rupture process; 2) development of nonlinear solvers to enable simulation of realistic crustal material behavior (e.g., viscoelastic, elastoplastic, elasto-viscoplastic); 3) automatic computation of gravitational prestresses; and 4) development of the capability to use a spherical mesh finite element framework for modeling the dynamics of large spatial areas. The work will involve no-cost collaborations with CIG funded personnel (Matt Knepley, Argonne National Lab), the USGS (Brad Aagaard) and a Caltech Ph.D. student. The code developments will be open-source and maintained in the CIG repository.