Numerical simulation of magma intrusion

Based on the user_HydraulicBlock example, a two-dimensional magma intrusion model is established. The model contains 31,000 solid units and 47,000 fluid units. A weak interlayer is established in the middle and lower part, which also has a higher permeability. A circular area is excavated from the lower part as a "magma room".

Initially, the units are cemented with each other, and the permeability after the cement fracture is set to increase by 100 times. A constant pressure of 5MPa is applied to the upper part of the rock formation, and a constant pressure of 100MPa is applied in the "magma chamber". Simulate fracture formation and fluid migration. Using a laptop computer for calculation, the simulation takes 2 hours and the simulation real world time is 26.4 seconds.

Dynamic diagram of pore pressure field (red represents higher pressure)

Velocity field animation (red represents higher element velocity, higher velocity when the tip ruptures)

Motion picture of displacement field (red represents higher displacement)

Simulate multiple ejection processes