Complete workflow for simulating a geological model with time and space dependent boundary conditions (advancing glacier)

In this three-part tutorial we look at the mechanical deformation of a sedimentary basin under an advancing glacier. In detail we show you a possibility how to

• create a mesh with Gmsh,
• implement space- and time-dependent boundary conditions,
• write an OGS input file,
• visualize the results with ParaView.

As a common programming language we use Python.

Part 1: Preprocessing

Part 2: Solving

Part 3: Postprocessing

Supplementary material

• mesh_basin.msh
• mesh_basin.py
• OGSinput_basin.prj
• timeBCs_glacier.py
• glacierclass.py
• history.sh
• msh2vtu -now integrated in ogstools
• Transient hydrodynamics within intercratonic sedimentary basins during glacial cycles, V. F. Bense M. A. Person, DOI: 10.1029/2007JF000969

NOTE: Due to ongoing development, msh2vtu.py was integrated in the ogs-repository for ogstools. In order to work with the tutorial and apply msh2vtu on mesh_basin.msh, ogstools need to be installed at first.

ogstools can be installed from PyPI using pip. It is recommended to set up a virtual environment beforehand:

python -m venv .venv
source .venv/bin/activate
pip install ogstools[ogs]

More detailed information about the usages and development of ogstools can be found in ogstools Documentation. After installation the tool can be applied as described in the first video-tutorial. To apply msh2vtu on the mesh_basin.msh the command should look like:

msh2vtu --ogs --rdcd mesh_basin.msh