Software

Primary development of

The FD-PDE framework uses finite difference staggered grids for solving partial differential equations (PDEs) for single-/two-phase flow magma dynamics and other applications in geodynamics.
The framework is based on PETSc. Governing equations are discretized as object PDEs (i.e., FDPDE Type), and the user only specifies coefficients, constitutive equations, and boundary conditions.
Models: mid-ocean ridge dynamics, melt extraction, small-scale convection.

Contributed to the development of

LaMEM (Lithosphere and Mantle Evolution Model) is a fast, parallel 3D numerical code for thermo-mechanical geodynamics problems with visco-elasto-plastic rheologies.
LaMEM is also based on PETSc.