Abstract:Reactive transport simulation on unstructured meshes can provide fundamental insight into the effect that geometric complexity of geologic structures has on fluid flow and development of reaction fronts. When applied to conditions ranging from ambient to hydrothermal and combined with compressible flow, accounting for geometric complexity provides an advantage for applications such as enhanced geothermal systems, carbon dioxide sequestration, hydrothermal ore formation, and radioactive waste disposal. We introduce CSMP–GEMS, a thermo–hydro and chemical multicomponent reactive transport code based on coupling of the Complex System Modeling Platform (CSMP) transport modeling framework with the GEMS3K chemical speciation solver. GEMS3K features a comprehensive suite of non-ideal activity and equation-of-state models of solution phases (aqueous electrolyte, gas and fluid mixtures, solid solutions). Current features include transient, compressible, single-phase advective and/or dispersive fluid flow, mass transport, heat transport in saturated porous media, and geochemical reactions in subsurface hydrothermal systems. We present two one-dimensional numerical experiments to compare CSMP–GEMS with the reactive transport codes OpenGeoSys–GEM and TOUGHREACT. Each experiment simulates calcite dissolution and dolomite precipitation during advection and hydrodynamic dispersion. One experiment corresponds to an existing isothermal (25∘C)$$(25\,^{\circ }\mathrm{C})$$ benchmark; the second explores the applicability of the codes to non-isothermal problems. We also present a two-dimensional example that illustrates the application of CSMP–GEMS on unstructured meshes that can represent complex geologic relations. The results suggest that all three codes are well suited to predicting fluid circulation, heat transport, and mineral stability within hydrothermal systems relevant to enhanced geothermal systems and carbon dioxide sequestration in deep aquifers. Self-consistent accounting for kinetic processes is a major advantage of TOUGHREACT, but published applications are restricted to orthogonal meshes, potentially limiting the applicability of TOUGHREACT to geometrically less complex natural systems. OpenGeoSys–GEM can operate on unstructured meshes that may include multiple element types, facilitating the examination of non-orthogonal domains. However, due to its reliance on the groundwater equations, OpenGeoSys–GEM may be best suited for application to systems in which flow includes dispersion/diffusion and is not compressible. CSMP–GEMS does not currently calculate reaction kinetics, but may be useful for application to geometrically complex systems.

What problem does this paper attempt to address?

The main problem that this paper attempts to solve is to provide fundamental insights through reactive - transport simulations on unstructured grids regarding the influence of the geometric complexity of geological structures on fluid flow and the development of reaction fronts. Specifically, the paper introduces CSMP - GEMS, a platform based on thermo - hydrological and chemical multi - component reactive - transport codes, aiming to combine the transport modeling framework of the Complex System Modeling Platform (CSMP) with the GEMS3K chemical species analytical solver to solve reactive - transport problems from environmental to hydrothermal conditions. By considering geometric complexity, this method has advantages in applications such as enhanced geothermal systems, carbon dioxide sequestration, hydrothermal ore deposit formation, and radioactive waste disposal. The paper compares the performance of CSMP - GEMS with other reactive - transport codes (such as OpenGeoSys - GEM and TOUGHREACT) through two one - dimensional numerical experiments that simulate the processes of calcite dissolution and dolomite precipitation during convection and hydrodynamic dispersion. One of the experiments corresponds to an existing isothermal (25 °C) benchmark test; the other experiment explores the applicability of these codes in non - isothermal problems. In addition, the paper also provides a two - dimensional example to demonstrate the application of CSMP - GEMS on unstructured grids, which can represent complex geological relationships. Through these experiments, the paper evaluates the performance of different codes in predicting fluid circulation, heat transfer, and mineral stability in hydrothermal systems, especially for carbon dioxide sequestration in enhanced geothermal systems and deep aquifers. The research results show that all three codes are suitable for these applications, but each has different advantages and limitations when dealing with specific types of problems. For example, TOUGHREACT has an advantage in handling kinetic processes, but its application is mainly limited to orthogonal grids, which may limit its application in geometrically complex natural systems. OpenGeoSys - GEM can operate unstructured grids on non - orthogonal domains, but may have limitations when dealing with compressible flow. Although CSMP - GEMS does not currently calculate reaction kinetics, it may be suitable for the simulation of geometrically complex systems.

Numerical Simulation of Reactive Fluid Flow on Unstructured Meshes

Using Toughreact to Model Reactive Fluid Flow and Geochemical Transport in Hydrothermal Systems

TOUGHREACT - A Simulation Program for Non-Isothermal Multiphase Reactive Geochemical Transport in Variably Saturated Geologic Media: Applications to Geothermal Injectivity and CO2 Geological Sequestration

TOUGHREACT: a New Code of the TOUGH Family for Non-Isothermal Multiphase Reactive Geochemical Transport in Variably Saturated Geologic Media

A PARALLEL SECOND-ORDER ADAPTIVE MESH ALGORITHM FOR REACTIVE FLOW IN GEOCHEMICAL SYSTEMS

Coupled Modeling of Non-Isothermal Multiphase Flow, Solutetransport and Reactive Chemistry in Porous and Fractured Media: 1. ModelDevelopment and Validation

Benchmarking of Reactive Transport Codes for 2D Simulations with Mineral Dissolution–precipitation Reactions and Feedback on Transport Parameters

TOUGHREACT User's Guide: A Simulation Program for Non-isothermal Multiphase Reactive Geochemical Transport in Variable Saturated Geologic Media

Implementation and Evaluation of Permeability-Porosity and Tortuosity-Porosity Relationships Linked to Mineral Dissolution-Precipitation

Core2d V4: A Code for Water Flow, Heat and Solute Transport, Geochemical Reactions, and Microbial Processes

Numerical Modelling of Fluids Mixing, Heat Transfer and Non-Equilibrium Redox Chemical Reactions in Fluid-Saturated Porous Rocks

TOUGHREACT Version 2.0: A simulator for subsurface reactive transport under non-isothermal multiphase flow conditions

Lawrence Berkeley National Laboratory Lawrence Berkeley National Laboratory Title TOURGHREACT : A Simulation Program for Non-isothermal Multiphase Reactive Geochemical Transport in Variably Saturated Geologic

Modeling Multiphase Non-Isothermal Fluid Flow and Reactive Geochemical Transport in Variably Saturated Fractured Rocks: 1. Methodology

Multiphase Reactive Transport Modeling – the Impacts of Two-Phase Flow on Mineral Reaction Rate

A General Method for Simulating Reactive Dissolution in Carbonate Rocks with Arbitrary Geometry

A Hybrid Immersed Boundary-Lattice Boltzmann/finite Difference Method for Coupled Dynamics of Fluid Flow, Advection, Diffusion and Adsorption in Fractured and Porous Media

An Approach on Parallel Simulation of Coupled Solute Transport and Geochemical Reaction

IORSim: A Mathematical Workflow for Field-Scale Geochemistry Simulations in Porous Media

Modeling Reactive Multiphase Flow and Transport of Concentrated Solutions

An overview of computational methods for chemical equilibrium and kinetic calculations for geochemical and reactive transport modeling