WP1-Numerical Models


Trajectories calculation in 2D for any stress field

Trajectories calculation in 2D based on Maccaferri et al, 2011 had been only performed for analytically-described stress field. We generalized the 2-D dike trajectory calculation for stress fields known only at the nodes of a mesh. This was applied for calculating magma intrusions trajectories below impact craters on the Moon (Michaut et al, 2020). This work correspond to deliverable D1.1.

Read more

2D numerical model coupling magma flow and fracturing

We provide a 2D numerical model that captures the coupling of fracturing at the crack tip and the transport of a viscous fluid. Built with the boundary element technique, our model allows for computation of the shape and velocity of a growing fluid-filled crack accounting for the viscosity of the fluid: the fluid flow induces a viscous pressure drop acting at the crack walls, and modifies the shape of the crack. The energy conservation equation provides the constraints to solve for the (...)

Read more

SAM: Simplified Analytical Model of Dyke Propagation in Three Dimensions

We provide a fast three-dimensional analytical model of dyke pathways based on principles of fracture mechanics. Dykes are driven by the interplay of an external stress field and the buoyancy force of the magma filling them.
The model is described in Mantiloni, L., Rivalta, E., & Davis, T. (2023). Mechanical modeling of pre-eruptive magma propagation scenarios at calderas. Journal of Geophysical Research: Solid Earth, 128, e2022JB025956. https://doi.org/10.1029/2022JB025956
You can (...)

Read more