cianwilson / fenicsx_subduction

Overview

The thermal structure of subduction zones is fundamental to our understanding of the physical and chemical processes that occur at active convergent plate margins. These include magma generation and related arc volcanism, shallow and deep seismicity and metamorphic reactions that can release fluids. Computational models can predict the thermal structure to great numerical precision when models are fully described but this does not guarantee accuracy or applicability. In a trio of companion papers (I, II, and III), the construction of thermal subduction zone models, their use in subduction zone studies, and their link to geophysical and geochemical observations were explored. This Jupyter Book reproduces the examples presented in part II where the finite element techniques that can be used to predict the thermal structure were discussed in an introductory fashion along with their verification and validation. Unlike the examples presented in part II (github, zenodo), which mostly used the finite element model builder TerraFERMA and a legacy version of the finite element library FEniCS, here we present these models using the latest version of FEniCS, FEniCSx.

Interactive tutorials

As this book has been published as a Jupyter Book, we provide interactive notebooks that can be run in the browser. To start such a notebook click the -symbol in the top right corner of the relevant tutorial.

Obtaining the software

If you would like to work with FEniCSx outside of the binder-notebooks, you need to install the FEniCS software. The recommended way of installing DOLFINx for new users is by using Docker. Docker is a software that uses containers to supply software across different kinds of operating systems (Linux, Mac, Windows). The first step is to install docker, following the instructions at their webpage.

Docker images

Compatible images are available in the Github Packages.

To use the notebooks in this tutorial with DOLFINx on your own computer, you should (using a terminal) clone the repository and change into that directory:

git clone https://github.com/cianwilson/fenicsx_subduction.git
cd fenicsx_subduction

Then run the following docker command:

  docker run --init --rm -p 8888:8888 --workdir /root/shared -v "$(pwd)":/root/shared ghcr.io/cianwilson/fenicsx_subduction:v0.8.0

Finally, copy and paste the URL printed in the terminal into a web-browser.

Alternatively, this image can used as an interactive terminal by running:

  docker run -it --rm -p 8888:8888 --workdir /root/shared -v "$(pwd)":/root/shared  --entrypoint="/bin/bash" ghcr.io/cianwilson/fenicsx_subduction:v0.8.0

Jupyter lab can then be started from within the docker container:

  jupyter lab --ip 0.0.0.0 --port 8888 --no-browser --allow-root

again copying and pasting the resulting URL into a web browser.

Acknowledgments

This Jupyter Book is based on the FEniCSx Tutorial by Jørgen S. Dokken.