This project aims to use realizations of heterogeneous hydraulic conductivity with a binary inclusion structures for numerical subsurface flow and transport simulations. Code is available to (i) generate and (ii) visualize random conductivity structures of binary inclusions. Furthermore, routines are provided (iii) to perform flow and transport simulations making use of the FEM-solver OpenGeoSys and (iv) to post-process results.

The project is organized as follows:

• README.md - description of the project
• results/ - folder with simulation results and plots
• src/ - folder containing the Python scripts of the project:
  • 01_plot_binary_inclusions.py - generate and visualize inclusion structures
  • 02_run_sim.py - run OGS 5 simulation on inclusion structures according to MADE setting
  • 03_run_ensemble.py - run OGS 5 ensemble of transport simulations
  • binary_inclusions.py - package containing inclusing structure classes
  • sim_processing.py - package containing post-processing routines

To make the example reproducible, we provide the following files:

• requirements.txt - requirements for pip to install all needed packages

You can contact us via a.zech@uu.nl.

MIT © 2020