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.

*.jou

Files used to construct the finite-element mesh using Cubit.

*.exo

Exodus II finite-element mesh files generated by Cubit.

*.spatialdb

Spatial database filesFiles associated with the spatial databases.

viz

Directory containing ParaView Python scripts and other files for visualizing results.

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.

Fig. 108 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

• Cubit Mesh
  • Geometry
  • Meshing using Journal Scripts
• Common Information
  • Metadata, Mesh, and Output
  • Physics
• Step 1: Magma inflation
  • Simulation parameter
  • Running the simulation
  • Visualizing the results
• Suggested Exercises