EGU22-5367, updated on 27 Mar 2022
EGU General Assembly 2022
© Author(s) 2022. This work is distributed under
the Creative Commons Attribution 4.0 License.

Verification of TRANsport Simulation Environment coupling with PHREEQC for reactive transport modelling

Svenja Steding1,2, Michael Kühn1,2, and Thomas Kempka1,2
Svenja Steding et al.
  • 1GFZ German Research Centre for Geosciences, Fluid Systems Modelling, Potsdam, Germany
  • 2University of Potsdam, Institute of Geosciences, Potsdam, Germany

Many types of geological subsurface utilisation are associated with fluid and heat flow as well as simultaneously occurring chemical reactions. For that reason, reactive transport models are required to understand and reproduce the governing processes. In this regard, reactive transport codes must be highly flexible to cover a wide range of applications, while being applicable by users without extensive programming skills at the same time. In this context, the TRANsport Simulation Environment (Kempka, 2020) was coupled with the geochemical reaction module PHREEQC (Parkhurst & Appelo, 2013), providing multiple features that make it applicable to complex reactive transport problems in various fields. Code readability is ensured by the applied high-level programming language Python which is relatively easy to learn compared to low-level programming languages. In the present study, common geochemical benchmarks are used to verify the numerical code implementation.

Currently, the coupled simulator can be used to investigate 3D single-phase fluid and heat flow as well as multicomponent solute transport in porous media. In addition to that, a wide range of equilibrium and nonequilibrium reactions can be considered. Chemical feedback on fluid flow is provided by adapting porosity and permeability of the porous media as well as fluid properties. Thereby, users are in full control of the underlying functions and equations of state. Both, the solution of the system of the partial differential equations and PHREEQC module, can be easily parallelised to increase computational efficiency.

The benchmarks used in the present study include density-driven flow as well as advective and diffusive reactive transport of solutes. Furthermore, porosity, permeability and diffusivity changes caused by kinetically controlled dissolution-precipitation reactions are considered to verify the main features of our reactive transport code. In future, the code implementation may be used to quantify processes encountered in different types of subsurface utilisation, such as geothermal energy production, geological energy, CO2 and nuclear waste storage.


Kempka, T. (2020). Verification of a Python-based TRANsport Simulation Environment for density-driven fluid flow and coupled transport of heat and chemical species. Adv. Geosci. 54, 67–77. (

Parkhurst, D.L.; Appelo, C.A.J. (2013). Description of Input and Examples for PHREEQC Version 3 - a Computer Program for Speciation, Batch-reaction, One-dimensional Transport, and Inverse Geochemical Calculations. In Techniques and Methods; Publisher: U.S. Geological Survey; Book 6, 497 pp. (

How to cite: Steding, S., Kühn, M., and Kempka, T.: Verification of TRANsport Simulation Environment coupling with PHREEQC for reactive transport modelling, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-5367,, 2022.