toughio 1.2.1

Pre- and post-processing Python library for TOUGH

toughio

TOUGH (Transport Of Unsaturated Groundwater and Heat) is a general purpose numerical simulation software designed for fluid and heat flows of multiphase, multicomponent fluid mixtures in porous and fractured media developed at Lawrence Berkeley National Laboratory. It solves mass and energy balance equations that describe fluid and heat flow in multiphase and multicomponent systems. TOUGH handles all types of multiphase and multicomponent flow systems since the governing equations for fluid and heat flow have the same mathematical form. The nature and properties of fluid mixtures are described by thermophysical variables (e.g. density, viscosity, enthalpy) which are provided by an equation-of-state (EOS) module.

toughio is an open-source library that provides tools to facilitate pre- and post-processing for TOUGH using the latest Python standards. It aims to make setting up of a TOUGH simulation user-friendly by relying on existing well-established Python packages:

• numpy: vectorized calculation of N-dimensional arrays,
• meshio: input/output for many mesh formats,
• pyvista: 3D plotting and mesh analysis through a streamlined interface for the Visualization Toolkit (VTK).

Note that the results of a TOUGH simulation are sensitive to the quality of the mesh (ideally, it should satisfy the orthogonality condition). A mesh that contains too many ill-shaped cells can potentially lead to unexpected results although the simulation converged successfully. toughio does not verify the quality of the mesh which is left to the discretion of the user.

Simulation of CO₂ leakage along a fault.
Mesh generated with Gmsh and animation exported by PyVista.

Features

Meshing:

• Create simple 2D cylindric or 2D/3D structured meshes similarly to TOUGH's built-in MESHMAKER,
• Import mesh generated by external softwares (e.g. Abaqus, FLAC3D, Gmsh, LaGriT),
• Export imported or generated mesh to a MESH file for TOUGH assuming conformity (and optionally write initial condition file INCON).

Pre-processing:

• Easily add initial conditions, boundary conditions or other physical properties (e.g. porosity, permeability) using the convenient class toughio.Mesh,
• Import previous TOUGH simulation input file into Python or define simulation parameters in a human-readable and JSONable dictionary,
• Export simulation parameters dictionary to a TOUGH input file.

Post-processing:

• Import outputs of a TOUGH simulation into Python,
• Visualize results directly in Python using pyvista or export the results to another format supported by meshio (e.g. VTK, Tecplot...),
• Create animations (GIF or MP4) or export results for all time steps to a XDMF file to be visualized in ParaView.

Installation

The recommended way to install toughio and all its dependencies is through the Python Package Index:

pip install toughio[full] --user

Otherwise, clone and extract the package, then run from the package location:

pip install .[full]

Usage

In Python, to read a mesh and write the corresponding TOUGH MESH file (without any pre-processing), simply do

import toughio

mesh = toughio.read_mesh(
    filename,
    file_format="flac3d",   # Optional, inferred from file extension otherwise
)
mesh.to_tough()             # Write MESH file

Parameters of a TOUGH simulation can be defined as a dictionary with specific keywords following the JSON standard, for instance

parameters = {
    "title": "Sample title",
    "eos": "eco2n",
    "isothermal":, False,
    "default": {            # Default rock properties
        "density": 2600.0,
        "porosity": 0.1,
        # "permeability", "conductivity", "specific_heat"...
    },
    "rocks": {
        "shale": {          # To overwrite default rock properties
            "capillarity": {
                "id": 1,
                "parameters": [0.0, 0.0, 1.0],
            },
            # same keywords as in "default"
        },
        # other materials
    },
    "options": {
        "n_cycle": 100,
        "t_max": 3.0 * 365.25 * 24.0 * 3600.0,
        # "t_ini", "t_steps", "t_step_max", "gravity", "eps1", "eps2"...
    },
    # "extra_options", "selections", "solver", "generators"...
}
toughio.write_input("INFILE", parameters)

TOUGH simulation output can also be imported into Python as a list of namedtuple (time, labels, data)

output = toughio.read_output(
    filename,
    file_format="tough",    # Optional, "tough" or "tecplot"
)

toughio is mainly intended to be used as a Python scripting library for TOUGH. Nevertheless, several utility command line scripts are available for users who are not familiar with Python. From a console or terminal, the user can execute the following scripts:

• toughio-export: export TOUGH simulation results to a file for visualization (VTK, VTU, Tecplot or XDMF),
• toughio-extract: extract results from TOUGH main output file and reformat as a TOUGH3 element output file (mostly useful for TOUGH2 output before calling toughio-export),
• toughio-merge: merge input file, MESH and/or INCON into a single file (for storage or sharing).