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.

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

Integrated Two-Dimensional Surface and Three-Dimensional Subsurface Contaminant Transport Model Considering Soil Erosion and Sorption

Coupled Finite-Volume Model for 2D Surface and 3D Subsurface Flows

A Physically-Based Integrated Numerical Model for Flow, Upland Erosion, and Contaminant Transport in Surface-Subsurface Systems

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

Analytical Model for Steady-State Solute Diffusion in Non-Isothermal Fractured Porous Media

Numerical Simulation of Reactive Fluid Flow on Unstructured Meshes

Pore-scale simulation of multicomponent multiphase reactive transport with dissolution and precipitation

A benchmark study on reactive two-phase flow in porous media: Part I - model description

Modeling of Multiphysical–Chemical Coupling for Coordinated Mining of Coal and Uranium in a Complex Hydrogeological Environment

Combinatorial Model of Solute Transport in Porous Media.

Coupled Reactive Two-Phase Model Involving Dissolution and Dynamic Porosity for Deformable Porous Media Based on Mixture Coupling Theory

A Mixture Coupling Theory based model for solute transport in deformable dual-porosity media

Reactive transport modeling for assessing coupled hydrogeochemical processes at interfaces in deep geological repositories: from the laboratory to the real world

Modeling of Chemical Tracers for Two-Phase Flow in Advective-Dominated Porous Media at Core Scale

Phase-field modeling and effective simulation of non-isothermal reactive transport

Analytical Model for Volatile Organic Compound Transport in the Coupled Vadose Zone-Groundwater System

Macrocyclization by nickel-catalyzed, ester-promoted, epoxide-alkyne reductive coupling: total synthesis of (-)-gloeosporone.

Microcontinuum approach to multiscale modeling of multiphase reactive flow during mineral dissolution

In silico studies of the mechanism of methanol oxidation by quinoprotein methanol dehydrogenase.

A discrete fracture matrix framework for simulating single-phase flow and non-isothermal reactive transport