2D Magma Reservoir Using Poroelasticity

The files are in the directory examples/magma-2d. The files and directories for this set of examples includes:

README.md:

README file containing a brief description of the various examples.

*.cfg:

PyLith parameter files.

generate_gmsh.py:

Python script to generate mesh using Gmsh.

generate_cubit.py:

Python script to generate mesh using Cubit.

*.msh:

Gmsh mesh files generated by Gmsh.

*.exo:

Exodus II finite-element mesh files generated by Cubit.

*.spatialdb:

Spatial database filesFiles associated with the spatial databases.

output:

Directory containing simulation output. It is created automatically when running the simulations.

Overview

This example demonstrates the use of poroelasticity to model magma flow in a conduit and reservoir in the crust.

Step 1:

Magma influx with displacement and pressure boundary conditions.

Step 2:

Same as Step 1 with evolution of porosity state variable.

Fig. 157 Diagram of geometry for magma reservoir. The domain extends from 0 km to +20 km in the x direction and from -20 km to 0 km in the y direction. We refer to the domain boundaries using the names shown in the diagram.

Example Workflow

• Gmsh Mesh
  • Geometry
  • Meshing using Python Script
• Cubit Mesh
  • Geometry
  • Meshing using Python script
• Common Information
  • Metadata, Mesh, and Output
  • Physics
• Step 1: Magma inflation
  • Simulation parameters
  • Running the simulation
  • Visualizing the results
• Step 2: Magma inflation with evolution of porosity
  • Simulation parameters
  • Running the simulation
  • Visualizing the results
• Suggested Exercises