Articles | Volume 13, issue 6
Method article
29 Jun 2022
Method article |  | 29 Jun 2022

An efficient partial-differential-equation-based method to compute pressure boundary conditions in regional geodynamic models

Anthony Jourdon and Dave A. May

Model code and software

Firedrake Firedrake team

FEniCS FEniCS team

Short summary
In this study we present a method to compute a reference pressure based on density structure in which we cast the problem in terms of a partial differential equation (PDE). We show in the context of 3D models of continental rifting that using the pressure as a boundary condition within the flow problem results in non-cylindrical velocity fields, producing strain localization in the lithosphere along large-scale strike-slip shear zones and allowing the formation and evolution of triple junctions.