Verification of a Python-based TRANsport Simulation Environment for density-driven fluid flow and coupled transport of heat and chemical species
GFZ German Research Centre for Geosciences, Telegrafenberg, 14473 Potsdam, Germany
Institute of Geosciences, University of Potsdam, Karl-Liebknecht-Str. 24–25,
14476 Potsdam, Germany
Related authors
Michael Kühn, Vaughan Stagpoole, Graham Paul D. Viskovic, and Thomas Kempka
Adv. Geosci., 65, 1–7, https://doi.org/10.5194/adgeo-65-1-2024, https://doi.org/10.5194/adgeo-65-1-2024, 2024
Short summary
Short summary
Waiwera is a small coastal village located on New Zealand's North Island above a geothermal reservoir. The origin of the warm water is not well understood. An inferred fault zone at the base of the reservoir is thought to channelise the uprising thermal water. The observed characteristic cold and hot water distribution in the system was simulated and the temperature profiles show an improved agreement in the near field around the centre of the reservoir.
Morgan Tranter, Svenja Steding, Christopher Otto, Konstantina Pyrgaki, Mansour Hedayatzadeh, Vasilis Sarhosis, Nikolaos Koukouzas, Georgios Louloudis, Christos Roumpos, and Thomas Kempka
Adv. Geosci., 58, 67–76, https://doi.org/10.5194/adgeo-58-67-2022, https://doi.org/10.5194/adgeo-58-67-2022, 2022
Short summary
Short summary
The quantification of the impacts on the environment and human health is a crucial prerequisite for geological sub-surface utilisation projects. With the presented approach, the shortcomings of using conceptually simplified models are substantially reduced, since subsurface complexities are accounted for. The transparency of the assessment basis should generally increase the acceptance of geoengineering projects, considered one of the crucial aspects for geological subsurface utilisation.
Christopher Otto, Svenja Steding, Morgan Tranter, Torsten Gorka, Mária Hámor-Vidó, Wioleta Basa, Krzysztof Kapusta, István Kalmár, and Thomas Kempka
Adv. Geosci., 58, 55–66, https://doi.org/10.5194/adgeo-58-55-2022, https://doi.org/10.5194/adgeo-58-55-2022, 2022
Short summary
Short summary
For a potential utilisation of coal resources located in Hungary, an assessment of groundwater pollution resulting from a water-borne contaminant pool has been undertaken. A sensitivity analysis was carried out by means of numerical simulations. Simulation results demonstrate that fluid flow via the regional faults is the main driver for a potential contamination of shallow groundwater aquifers. A parameter correlation analysis is presented.
Elena Chabab, Michael Kühn, and Thomas Kempka
Adv. Geosci., 58, 47–54, https://doi.org/10.5194/adgeo-58-47-2022, https://doi.org/10.5194/adgeo-58-47-2022, 2022
Short summary
Short summary
The present study, uses density-driven flow and transport models to evaluate mechanisms of saline water intrusion from deep aquifers into the freshwater column used for drinking water supply under different boundary conditions and for a specific site in the German Federal State of Brandenburg. Results show that mainly decreasing groundwater recharge leads to increased and earlier salinisation which highlights the need for waterworks to initiate effective countermeasures quickly and in time.
Michael Kühn, Melissa Präg, Ivy Becker, Christoph Hilgers, Andreas Grafe, and Thomas Kempka
Adv. Geosci., 58, 31–39, https://doi.org/10.5194/adgeo-58-31-2022, https://doi.org/10.5194/adgeo-58-31-2022, 2022
Short summary
Short summary
The geothermal hot water reservoir below the small town of Waiwera in New Zealand has been known to the indigenous Maori for many centuries. Overproduction by European immigrants led to a water level decrease and consequently artesian flow from the wells and the seeps on the beach ceased. The Te Kaunihera o Tāmaki Makaurau Auckland Council established the Waiwera Thermal Groundwater Allocation and Management Plan to allow the geothermal system to recover.
Thomas Kempka, Svenja Steding, and Michael Kühn
Adv. Geosci., 58, 19–29, https://doi.org/10.5194/adgeo-58-19-2022, https://doi.org/10.5194/adgeo-58-19-2022, 2022
Short summary
Short summary
The TRANSPORT Simulation Environment (TRANSPORTSE) was coupled with the geochemical reaction module PHREEQC, providing multiple new features that make it applicable to complex reactive transport problems in various geoscientific fields. Two computationally demanding and complex geochemical benchmarks were used in the present study to successfully verify the code implementation.
Maria Wetzel, Thomas Kempka, and Michael Kühn
Adv. Geosci., 58, 1–10, https://doi.org/10.5194/adgeo-58-1-2022, https://doi.org/10.5194/adgeo-58-1-2022, 2022
Short summary
Short summary
Porosity-permeability relations are simulated for a precipitation-dissolution cycle in a virtual sandstone. A hysteresis in permeability is observed depending on the geochemical process and dominating reaction regime, whereby permeability varies by more than two orders of magnitude. Controlling parameters for this hysteresis phenomenon are the closure and re-opening of micro-scale flow channels, derived from changes in pore throat diameter and connectivity of the pore network.
Michael Kühn, Natalie C. Nakaten, and Thomas Kempka
Adv. Geosci., 54, 173–178, https://doi.org/10.5194/adgeo-54-173-2020, https://doi.org/10.5194/adgeo-54-173-2020, 2020
Short summary
Short summary
Energy supply in Germany is subject to a profound change. The present paper addresses the German potential of storing excess energy from renewable power sources in the geological subsurface. Wind and solar electricity can be transformed into hydrogen, and with carbon dioxide subsequently into methane. The current potential for combined subsurface storage of methane and carbon dioxide allows to store far more than required to date and is estimated to provide the entire coverage in 2050.
Maria Wetzel, Thomas Kempka, and Michael Kühn
Adv. Geosci., 54, 33–39, https://doi.org/10.5194/adgeo-54-33-2020, https://doi.org/10.5194/adgeo-54-33-2020, 2020
Elena Tillner, Maria Langer, Thomas Kempka, and Michael Kühn
Hydrol. Earth Syst. Sci., 20, 1049–1067, https://doi.org/10.5194/hess-20-1049-2016, https://doi.org/10.5194/hess-20-1049-2016, 2016
Short summary
Short summary
The degree of shallow aquifer salinisation triggered by fluid injection into deeper brine-bearing aquifers and brine upward migration through hydraulically conductive faults strongly depends on the regional depth of the freshwater-saltwater boundary, since displaced brines originate only from the upper fault damage zones in the study area. The highest local salinity increase in shallow aquifers occurs in case of closed model boundaries and low fault damage zone volumes.
M. De Lucia, T. Kempka, and M. Kühn
Geosci. Model Dev., 8, 279–294, https://doi.org/10.5194/gmd-8-279-2015, https://doi.org/10.5194/gmd-8-279-2015, 2015
Michael Kühn, Vaughan Stagpoole, Graham Paul D. Viskovic, and Thomas Kempka
Adv. Geosci., 65, 1–7, https://doi.org/10.5194/adgeo-65-1-2024, https://doi.org/10.5194/adgeo-65-1-2024, 2024
Short summary
Short summary
Waiwera is a small coastal village located on New Zealand's North Island above a geothermal reservoir. The origin of the warm water is not well understood. An inferred fault zone at the base of the reservoir is thought to channelise the uprising thermal water. The observed characteristic cold and hot water distribution in the system was simulated and the temperature profiles show an improved agreement in the near field around the centre of the reservoir.
Morgan Tranter, Svenja Steding, Christopher Otto, Konstantina Pyrgaki, Mansour Hedayatzadeh, Vasilis Sarhosis, Nikolaos Koukouzas, Georgios Louloudis, Christos Roumpos, and Thomas Kempka
Adv. Geosci., 58, 67–76, https://doi.org/10.5194/adgeo-58-67-2022, https://doi.org/10.5194/adgeo-58-67-2022, 2022
Short summary
Short summary
The quantification of the impacts on the environment and human health is a crucial prerequisite for geological sub-surface utilisation projects. With the presented approach, the shortcomings of using conceptually simplified models are substantially reduced, since subsurface complexities are accounted for. The transparency of the assessment basis should generally increase the acceptance of geoengineering projects, considered one of the crucial aspects for geological subsurface utilisation.
Christopher Otto, Svenja Steding, Morgan Tranter, Torsten Gorka, Mária Hámor-Vidó, Wioleta Basa, Krzysztof Kapusta, István Kalmár, and Thomas Kempka
Adv. Geosci., 58, 55–66, https://doi.org/10.5194/adgeo-58-55-2022, https://doi.org/10.5194/adgeo-58-55-2022, 2022
Short summary
Short summary
For a potential utilisation of coal resources located in Hungary, an assessment of groundwater pollution resulting from a water-borne contaminant pool has been undertaken. A sensitivity analysis was carried out by means of numerical simulations. Simulation results demonstrate that fluid flow via the regional faults is the main driver for a potential contamination of shallow groundwater aquifers. A parameter correlation analysis is presented.
Elena Chabab, Michael Kühn, and Thomas Kempka
Adv. Geosci., 58, 47–54, https://doi.org/10.5194/adgeo-58-47-2022, https://doi.org/10.5194/adgeo-58-47-2022, 2022
Short summary
Short summary
The present study, uses density-driven flow and transport models to evaluate mechanisms of saline water intrusion from deep aquifers into the freshwater column used for drinking water supply under different boundary conditions and for a specific site in the German Federal State of Brandenburg. Results show that mainly decreasing groundwater recharge leads to increased and earlier salinisation which highlights the need for waterworks to initiate effective countermeasures quickly and in time.
Michael Kühn, Melissa Präg, Ivy Becker, Christoph Hilgers, Andreas Grafe, and Thomas Kempka
Adv. Geosci., 58, 31–39, https://doi.org/10.5194/adgeo-58-31-2022, https://doi.org/10.5194/adgeo-58-31-2022, 2022
Short summary
Short summary
The geothermal hot water reservoir below the small town of Waiwera in New Zealand has been known to the indigenous Maori for many centuries. Overproduction by European immigrants led to a water level decrease and consequently artesian flow from the wells and the seeps on the beach ceased. The Te Kaunihera o Tāmaki Makaurau Auckland Council established the Waiwera Thermal Groundwater Allocation and Management Plan to allow the geothermal system to recover.
Thomas Kempka, Svenja Steding, and Michael Kühn
Adv. Geosci., 58, 19–29, https://doi.org/10.5194/adgeo-58-19-2022, https://doi.org/10.5194/adgeo-58-19-2022, 2022
Short summary
Short summary
The TRANSPORT Simulation Environment (TRANSPORTSE) was coupled with the geochemical reaction module PHREEQC, providing multiple new features that make it applicable to complex reactive transport problems in various geoscientific fields. Two computationally demanding and complex geochemical benchmarks were used in the present study to successfully verify the code implementation.
Maria Wetzel, Thomas Kempka, and Michael Kühn
Adv. Geosci., 58, 1–10, https://doi.org/10.5194/adgeo-58-1-2022, https://doi.org/10.5194/adgeo-58-1-2022, 2022
Short summary
Short summary
Porosity-permeability relations are simulated for a precipitation-dissolution cycle in a virtual sandstone. A hysteresis in permeability is observed depending on the geochemical process and dominating reaction regime, whereby permeability varies by more than two orders of magnitude. Controlling parameters for this hysteresis phenomenon are the closure and re-opening of micro-scale flow channels, derived from changes in pore throat diameter and connectivity of the pore network.
Michael Kühn, Natalie C. Nakaten, and Thomas Kempka
Adv. Geosci., 54, 173–178, https://doi.org/10.5194/adgeo-54-173-2020, https://doi.org/10.5194/adgeo-54-173-2020, 2020
Short summary
Short summary
Energy supply in Germany is subject to a profound change. The present paper addresses the German potential of storing excess energy from renewable power sources in the geological subsurface. Wind and solar electricity can be transformed into hydrogen, and with carbon dioxide subsequently into methane. The current potential for combined subsurface storage of methane and carbon dioxide allows to store far more than required to date and is estimated to provide the entire coverage in 2050.
Maria Wetzel, Thomas Kempka, and Michael Kühn
Adv. Geosci., 54, 33–39, https://doi.org/10.5194/adgeo-54-33-2020, https://doi.org/10.5194/adgeo-54-33-2020, 2020
Elena Tillner, Maria Langer, Thomas Kempka, and Michael Kühn
Hydrol. Earth Syst. Sci., 20, 1049–1067, https://doi.org/10.5194/hess-20-1049-2016, https://doi.org/10.5194/hess-20-1049-2016, 2016
Short summary
Short summary
The degree of shallow aquifer salinisation triggered by fluid injection into deeper brine-bearing aquifers and brine upward migration through hydraulically conductive faults strongly depends on the regional depth of the freshwater-saltwater boundary, since displaced brines originate only from the upper fault damage zones in the study area. The highest local salinity increase in shallow aquifers occurs in case of closed model boundaries and low fault damage zone volumes.
M. De Lucia, T. Kempka, and M. Kühn
Geosci. Model Dev., 8, 279–294, https://doi.org/10.5194/gmd-8-279-2015, https://doi.org/10.5194/gmd-8-279-2015, 2015
Cited articles
Afanasyev, A.: Numerical modelling of solute flow dispersion in porous media
using simulator MUFITS, J. Phys. Conf. Ser., 1129, 012002,
https://doi.org/10.1088/1742-6596/1129/1/012002, 2018. a
Anderson, T. A., Liu, H., Kuper, L., Totoni, E., Vitek, J., and Shpeisman, T.:
Parallelizing Julia with a Non-Invasive DSL, in: 31st European Conference
on Object-Oriented Programming (ECOOP 2017), edited by: Müller, P.,
Vol. 74 of Leibniz International Proceedings in Informatics (LIPIcs),
4:1–4:29, Schloss Dagstuhl–Leibniz-Zentrum fuer Informatik, Dagstuhl,
Germany, https://doi.org/10.4230/LIPIcs.ECOOP.2017.4, 2017. a
Charlton, S. R. and Parkhurst, D. L.: Modules based on the geochemical model
PHREEQC for use in scripting and programming languages, Comput.
Geosci., 37, 1653–1663, https://doi.org/10.1016/j.cageo.2011.02.005, 2011. a
Childs, H., Brugger, E., Whitlock, B., Meredith, J., Ahern, S., Pugmire, D.,
Biagas, K., Miller, M., Harrison, C., Weber, G. H., Krishnan, H., Fogal, T.,
Sanderson, A., Garth, C., Bethel, E. W., Camp, D., Rübel, O., Durant, M.,
Favre, J. M., and Navrátil, P.: VisIt: An End-User Tool For Visualizing
and Analyzing Very Large Data, in: High Performance Visualization–Enabling
Extreme-Scale Scientific Insight, 357–372, Lawrence Berkeley National Lab. (LBNL), Berkeley, CA, United States, 2012. a
Crank, J. and Nicolson, P.: A practical method for numerical evaluation of
solutions of partial differential equations of the heat-conduction type,
Math. Proc. Cambridge, 43, 50–67,
https://doi.org/10.1017/S0305004100023197, 1947. a
De Lucia, M., Kempka, T., and Kühn, M.: A coupling alternative to reactive transport simulations for long-term prediction of chemical reactions in heterogeneous CO2 storage systems, Geosci. Model Dev., 8, 279–294, https://doi.org/10.5194/gmd-8-279-2015, 2015. a
De Lucia, M., Kempka, T., Jatnieks, J., and Kühn, M.: Integrating surrogate
models into subsurface simulation framework allows computation of complex
reactive transport scenarios, Enrgy. Proced., 125, 580–587,
https://doi.org/10.1016/j.egypro.2017.08.200, 2017. a
Elder, J.: Numerical experiments with free convection in a vertical slot,
J. Fluid Mech., 24, 823–843, https://doi.org/10.1017/S0022112066001022,
1966. a
Elder, J.: Transient convection in a porous medium, J. Fluid Mech.,
27, 609–623, https://doi.org/10.1017/S0022112067000576, 1967. a
Flemisch, B., Darcis, M., Erbertseder, K., Faigle, B., Lauser, A., Mosthaf, K.,
Müthing, S., Nuske, P., Tatomir, A., Wolff, M., and Helmig, R.: DuMux: DUNE
for multi-phase, component, scale, physics, … flow and transport in porous
media, Adv. Water Resour., 34, 1102–1112,
https://doi.org/10.1016/j.advwatres.2011.03.007, 2011. a
Goodwin, D. G., Moffat, H. K., and Speth, R. L.: Cantera: An Object-oriented
Software Toolkit for Chemical Kinetics, Thermodynamics, and Transport
Processes, version
2.3.0, available at: https://www.cantera.org (last access: 13 October 2020), Zenodo, https://doi.org/10.5281/zenodo.170284, 2017. a
Herrera, P.: PyEVTK: A self-contained Python module to write binary VTK files, version 0.2.0, available at:
https://github.com/paulo-herrera/PyEVTK (last access: 13 October 2020), 2017. a
Hunter, J. D.: Matplotlib: A 2D graphics environment, Compu. Sci.
Eng., 9, 90–95, 2007. a
Koch, T., Gläser, D., Weishaupt, K., Ackermann, S., Beck, M., Becker, B.,
Burbulla, S., Class, H., Coltman, E., Emmert, S., Fetzer, T., Grüninger, C.,
Heck, K., Hommel, J., Kurz, T., Lipp, M., Mohammadi, F., Scherrer, S.,
Schneider, M., Seitz, G., Stadler, L., Utz, M., Weinhardt, F., and Flemisch,
B.: DuMux 3 – an open-source simulator for solving flow
and transport problems in porous media with a focus on model coupling,
Comput. Math. Appl.,
https://doi.org/10.1016/j.camwa.2020.02.012, online first, 2020. a
Kolditz, O., Nagel, T., Shao, H., Wang, W., and Bauer, S.:
Thermo-Hydro-Mechanical-Chemical Processes in Fractured Porous Media:
Modelling and Benchmarking. From Benchmarking to Tutoring, Springer, Cham,
Springer International Publishing AG 2018, https://doi.org/10.1007/978-3-319-68225-9,
2018. a
Lam, S. K., Pitrou, A., and Seibert, S.: Numba: A LLVM-Based Python JIT
Compiler, in: Proceedings of the Second Workshop on the LLVM Compiler
Infrastructure in HPC, LLVM ’15, Association for Computing Machinery, New
York, NY, USA, https://doi.org/10.1145/2833157.2833162, 2015. a, b
Oliphant, T. E.: A guide to NumPy, Vol. 1, Trelgol Publishing USA, 2006. a
Otto, C. and Kempka, T.: Synthesis Gas Composition Prediction for Underground
Coal Gasification Using a Thermochemical Equilibrium Modeling Approach,
Energies, 13, 1171, https://doi.org/10.3390/en13051171, 2020. a
Parkhurst, D. L. and Appelo, C. A. J.: Description of input and examples for
PHREEQC version 3 – A computer program for speciation, batch-reaction,
one-dimensional transport, and inverse geochemical calculations,
available at: https://pubs.usgs.gov/tm/06/a43 (last access: 13 October 2020), 2013. a
Ramachandran, P. and Varoquaux, G.: Mayavi: 3D visualization of scientific
data, Comput. Sci. Eng., 13, 40–51,
https://doi.org/10.1109/MCSE.2011.35, 2011. a
Ségol, G.: Classic groundwater simulations: proving and improving numerical
models, Prentice Hall, New Jersey, United States, 1994. a
Steding, S., Zirkler, A., and Kühn, M.: Geochemical reaction models quantify
the composition of transition zones between brine occurrence and unaffected
salt rock, Chem. Geol., 532, 119349,
https://doi.org/10.1016/j.chemgeo.2019.119349, 2020. a
Steefel, C., Appelo, C., Arora, B., Jacques, D., Kalbacher, T., Kolditz, O.,
Lagneau, V., Lichtner, P., Mayer, K., Meeussen, J., Molins, S., Moulton, D.,
Shao, H., Simunek, J., Spycher, N., Yabusaki, S., and Yeh, G.: Reactive
transport codes for subsurface environmental simulation, Comput.
Geosci., 19, 445–478, https://doi.org/10.1007/s10596-014-9443-x, 2015. a
Sullivan, C. B. and Kaszynski, A.: PyVista: 3D plotting and mesh analysis
through a streamlined interface for the Visualization Toolkit (VTK),
Journal of Open Source Software, 4, 1450, https://doi.org/10.21105/joss.01450, 2019. a
Theis, C. V.: The relation between the lowering of the piezometric surface and
the rate and duration of discharge of a well using ground-water storage, EOS
T. Am. Geophys. Un., 16, 519–524, 1935. a
Trefry, M. G. and Muffels, C.: FEFLOW: A Finite-Element Ground Water Flow and
Transport Modeling Tool, Groundwater, 45, 525–528,
https://doi.org/10.1111/j.1745-6584.2007.00358.x, 2007. a
Utkarsh, A.: The ParaView Guide: A Parallel Visualization Application, Kitware Inc., United States, 2015. a
Van Rossum, G. and Drake, F. L.: Python 3 Reference Manual, CreateSpace, Scotts
Valley, CA, 2009. a
Voss, C. I. and Souza, W. R.: Variable density flow and solute transport
simulation of regional aquifers containing a narrow freshwater-saltwater
transition zone, Water Resour. Res., 23, 1851–1866,
https://doi.org/10.1029/WR023i010p01851, 1987. a, b
Young, D. M.: A bound for the optimum relaxation factor for the successive
overrelaxation method, Numer. Math., 16, 408–413,
https://doi.org/10.1007/BF02169150, 1971. a
Short summary
The TRANsport Simulation Environment (TRANSE) has been developed to improve the flexibility for coupling chemical libraries with fluid flow and the transport of heat and chemical species. The Python-based implementation of TRANSE enables users not experienced in low-level programming languages (e.g., C, C++ or FORTRAN) to undertake required code modifications and integrate chemical modules as required. TRANSE has been successfully verified against benchmarks on density-driven fluid flow.
The TRANsport Simulation Environment (TRANSE) has been developed to improve the flexibility for...