Towards a staggered-grid finite difference code for modelling magmatic systems

The in-depth behaviour of magmatic systems is still poorly constrained due to their lack of accessibility and the difficulty of finding good analogue representations. However, progresses in the field of geological numerical modelling can allow to better understand and interpret those constraints. The MAGMA project aims on developing tools and software for studying a range of magmatic processes in the lithosphere. On the way to building a general framework able to model the behaviour of a fluid-solid chemically coupled magmatic system, we present here the current development of a mechanical staggered-grid finite difference 2D code using robust analytical linear and non-linear solvers via the PETSc infrastructure, able to run in parallel on highly performant computers. The mesh is assembled using the recently developed DMStag framework, which is part of PETSc. This code solves the Stokes equations for elasto-visco-plastic rheologies, including tensile plasticity essential to the development of dyke structures. Here, we will present the equations and implementations used, and show initial results and benchmarks.