Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reactive transport model”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

A Multicomponent Reactive Transport Model for Integrated Surface‐Subsurface Hydrology Problems

Abstract Despite the widespread use of integrated hydrology models in a variety of applications, consideration of multicomponent reactive transport is still not common. The implementation of these processes requires coupling transport at the surface‐subsurface interface and efficient solution of the non‐linear geochemical model that is consistent with the integrated hydrology solution. The Advanced Terrestrial Simulator provides a flexible multiphysics framework that facilitated this process. In this work, the integrated reactive transport process kernel (PK) was weakly coupled to the integrated hydrology PK. In turn, integrated transport and reactions were coupled using an operator splitting approach. This splitting enabled an explicit solution of the integrated transport problem, including a novel algorithm to calculate exchange fluxes across the surface‐subsurface interface and a point‐by‐point solution of the geochemical problem. Geochemical capabilities were added using well‐established external codes, but rather than using a custom interface to each, a generic interface was used that clearly specifies the variables and operations used by the chemistry PK. The implementation is demonstrated with two example simulations: transport of a tracer in a soil column as it saturates over time and water ponds on the surface and reactive transport in a hillslope driven by successive wet‐dry cycles that result in infiltration, runoff and exfiltration processes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Effects of Dissolution Regimes on Flow Channelization and Solute Transport in 3D Fracture Networks: Insights From Graph‐Based Reactive Transport Modeling

We investigate how mineral dissolution reshapes flow pathways and solute transport in three‐dimensional discrete fracture networks using a computationally efficient graph‐based reactive transport model. The DFNs are inspired by field‐site observations of fractured carbonate and represent realistic connectivity and structural heterogeneity. Flow is simulated with the Reynolds equation, and dissolution follows first‐order kinetics with diffusive limitations captured through an effective mass‐transfer coefficient. By systematically varying two key dimensionless parameters, the effective Damköhler number (Da), governing reaction versus advection rates, and a transport parameter (Da), analogous to the Thiele modulus, distinct flow channelization regimes emerge: mildly channelized at low G, highly channelized at intermediate Da, and extreme wormhole formation at high Da and low G. Eulerian and Lagrangian analyses, including breakthrough curves, particle tortuosity, dispersivity, and flow channeling indicators quantitatively characterize the progression of dissolution‐driven network restructuring. Across all regimes, initial fracture heterogeneity persists. The results underscore how the interplay between this initial structure, advection, reaction, and diffusion critically shapes subsurface flow pathways, with implications for applications ranging from groundwater remediation to enhanced geothermal systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Reactive Transport Modeling of Aquifer Thermal Energy Storage System at Stockton, NJ

This is the modeling data (input/output files of TOUGHREACT 4.10) used to simulate the reactive transport processes of the Aquifer Thermal Energy Storage (ATES) operations at Stockton University, NJ. Readme.txt lists all the files. TOUGHREACT 4.10 requires to reproduce the modeling output. The modeling data in this submission is related to the Aquifer Injection for Energy Storage purposes outlined in "Reactive Transport Modeling of Aquifer Thermal Energy Storage System at Stockton, NJ During Seasonal Operations".

15 GEOTHERMAL ENERGY↗
DOE OSTI2025

Reactive Transport Modeling of Hydrogen Production from Serpentinization of Olivine

Hydrogen production from serpentinization of ultramafic rocks represents a promising natural pathway for generating carbon-free energy, yet its kinetics and controlling factors remain incompletely understood. A key challenge in advancing serpentinization research lies in the heterogeneity of porosity and permeability in rocks, which leads to nonuniform fluid velocity fields, as well as uncertainties in estimating reactive surface area and identifying appropriate mineral reaction equilibria. Additional complexities arise from the role of dissolved SiO 2 , Fe 2+ /Fe 3+ partitioning, and the limited effect of pH variations within the strongly alkaline regime on hydrogen yields. These challenges hinder straightforward extrapolation from laboratory tests to practical applications of hydrogen production from natural rocks. Here, in this work, we address these questions using a simulation-based reactive transport modeling framework calibrated against controlled laboratory experiments reported elsewhere. The model couples geochemical kinetics, multiphase flow, and mineralogical feedbacks, enabling systematic evaluation of how surface area, dissolved silica concentration, Fe redox state, temperature, and pressure govern serpentinization and H2 generation. We find that surface area exerts the strongest control on reaction rates and hydrogen yields, while Fe 2+ /Fe 3+ ratios act as secondary modulators. Elevated dissolved silica concentrations suppress hydrogen production but accelerate serpentine precipitation, whereas increasing pH beyond 12 within the strongly alkaline regime produces only marginal gains. Finally, we demonstrate that integrating targeted experiments with calibrated simulations offers a powerful and efficient approach for predicting hydrogen yields and assessing parameter trade-offs in industrial-scale applications. This integration can substantially reduce the experimental burden while improving predictive capability, thereby enhancing both the mechanistic understanding and the practical feasibility of hydrogen production from serpentinization.

08 HYDROGEN↗
NASA NTRS2017

Using Reactive Transport Modeling to Understand Formation of the Stimson Sedimentary Unit and Altered Fracture Zones at Gale Crater, Mars

Water flowing through sediments at Gale Crater, Mars created environments that were likely habitable, and sampled basin-wide hydrological systems. However, many questions remain about these environments and the fluids that generated them. Measurements taken by the Mars Science Laboratory Curiosity of multiple fracture zones can help constrain the environments that formed them because they can be compared to nearby associated parent material (Figure 1). For example, measurements of altered fracture zones from the target Greenhorn in the Stimson sandstone can be compared to parent material measured in the nearby Big Sky target, allowing constraints to be placed on the alteration conditions that formed the Greenhorn target from the Big Sky target. Similarly, CheMin measurements of the powdered < 150 micron fraction from the drillhole at Big Sky and sample from the Rocknest eolian deposit indicate that the mineralogies are strikingly similar. The main differences are the presence of olivine in the Rocknest eolian deposit, which is absent in the Big Sky target, and the presence of far more abundant Fe oxides in the Big Sky target. Quantifying the changes between the Big Sky target and the Rocknest eolian deposit can therefore help us understand the diagenetic changes that occurred forming the Stimson sedimentary unit. In order to interpret these aqueous changes, we performed reactive transport modeling of 1) the formation of the Big Sky target from a Rocknest eolian deposit-like parent material, and 2) the formation of the Greenhorn target from the Big Sky target. This work allows us to test the relationships between the targets and the characteristics of the aqueous conditions that formed the Greenhorn target from the Big Sky target, and the Big Sky target from a Rocknest eolian deposit-like parent material.

Hausrath, E. M.↗
DOE OSTI2023

Reactive Transport Modeling of Anthropogenic Carbon Mineralization in Stacked Columbia River Basalt Reservoirs

Numerical simulation of CO2 storage in basalts and related reactive lithologies requires modeling complex, coupled hydrologic and chemical processes, including multi-phase flow and transport, partitioning of CO2 into the aqueous phase, and chemical interactions with aqueous fluids and rock minerals. We conducted reactive transport simulations of the Wallula pilot-scale CO2 injection into the flow tops of the Grande Ronde Basalt using our PNNL STOMP-CO2 simulator with the ECKEChem reactive module. Our mineralization simulation of the ~1,000 tons of injected CO2 into the interflow zones was based on the hydrologic transport model we previously developed. For this work, the simulations considered geochemical reactions involving the basalt components, precipitates, formation brine, and injected CO2. In our benchmark case, carbonate minerals precipitated, resulting in ~20% of the CO2 being mineralized in 10 years. Increasing the reaction rate of a single primary mineral phase (clinopyroxene) by an order of magnitude resulted in a carbon mineralization reaction extent of ~90% over the same time interval. Based on these initial sensitivity analysis results, it is clear that a thorough understanding of primary mineral dissolution rates is required for accurately predicting long-term fate and transport of injected CO2 into basalt formations. Our reactive transport numerical simulations will be key components of commercial-scale CO2 storage operation permitting, de-risking, and optimization in mafic and ultramafic reservoirs.

Cao, Ruoshi↗
DOE OSTI2022

A radioisotope - enabled reactive transport model for deep vadose zone carbon

In mountainous regions, which constitute the principle source of recharge to major rivers and regional aquifers, infiltration occurs through fractured, partially saturated, weathered bedrock that acts as a boundary layer between saturated aquifers and surface soil. Commonly this deep vadose zone (DVZ) is many meters thick, and yet its role in regulating the generation, retention and mobility of reactive solutes, including nutrients, contaminants and weathering products, is largely unknown. In particular, many of the key reactions that drive the formation of the weathered DVZ and the quality of water moving through it are redox processes, regulated by the availability of organic carbon and oxygen below the soil layer. The role of the DVZ is thus also poorly constrained in the context of carbon stocks and mobility, particularly in lithologies that are naturally high in organic carbon, such as shales. The overarching hypothesis of this study is that upland regions developed in geologic settings with abundant petrogenic carbon store and actively cycle carbon in the weathered DVZ below the soil and above the water table at rates that are significant and currently unconstrained. In order to quantify this cycling, the current study combines novel instrumentation techniques allowing new direct sampling of DVZ systems with advanced numerical reactive transport simulations of carbon transport and transformation. Critically, these simulations will explicitly treat the three isotopes of carbon (the abundant 12C, the stable rare 13C and the radioactive 14C) in a unified framework, thus clearly parsing between the contributions of modern surface derived carbon and lithologic carbon sources in integrated measurements of fluid and gas phase fluxes. This novel model capability will be applied to test the role of DVZ carbon cycling as a regulator of water quality and geological weathering in two complementary field sites both located in organic carbon rich shale lithologies. The first is the Eel River Critical Zone Observatory (ERCZO) in Mendocino County, California, and the second is the Lawrence Berkeley National Laboratory Watershed Function Scientific Focus Area (SFA) in the East River watershed, near Crested Butte, Colorado. At the ERCZO, a novel vadose zone monitoring system has been installed in a 20 m thick, partially saturated, weathered shale hillslope, and preliminary data already indicate substantial CO2 flux generated many meters below the soil surface. At the SFA field site, an instrumented hillslope transect indicates a more complex multi-dimensional fluid and solute transport regime, which will serve as a key test of the calibrated models. Collectively, this project will advance understanding of the cycling of carbon belowground and in relation to transport pathways across the poorly constrained DVZ characteristic of primary water recharge areas. The key product of this work will be enhanced isotope simulation capabilities that are robust and publicly available for application across a broad diversity of systems.

58 GEOSCIENCES↗
DOE OSTI2025

Geologic hydrogen as an emerging fuel: experimental insights, thermodynamics, kinetics, and reactive transport modeling

Geologic hydrogen (GeoH 2 ) is emerging as a viable clean energy source. It is largely produced through serpentinization, a geological process in which ultramafic rocks react with water under suitable temperature and pressure. Here, this review synthesizes the current understanding of H 2 generation by serpentinization, with an emphasis on reaction mechanisms, kinetics, and thermodynamics, as well as on modeling flow and transport of reacting fluids in geological formations. We describe the role of mineral assemblages, such as olivine and pyroxene, fluid-rock interactions, and catalytic surfaces, in influencing GeoH 2 yield and reaction rates. By integrating models of reaction kinetics, subsurface reactive flow and transport, and the serpentinization process, and by accounting for the thermodynamic state of the system, this review aims to guide future GeoH 2 research and to evaluate the potential of natural hydrogen as a sustainable clean energy source.

Moradi, Rasoul [Univ. of Southern California, Los ↗
DOE OSTI2025

Permeate fluxes from desalination of brines and produced waters: A reactive transport modeling study

The increasing interest in the use of membrane systems to desalinate inland brackish water, agricultural drainage, and industrially produced wastewater demands improved means of predicting desalination system performance under variable feedwater compositions. The interaction among water flow, solute transport, and chemical composition in these systems impacts permeate flux evolution. Here, an established multicomponent reactive transport simulator that accounts for these coupled processes is applied to compute osmotic pressure and permeate fluxes in reverse osmosis (RO) systems. The model is first validated by predicting permeate fluxes for a set of benchtop crossflow experiments subject to a range of feed flow rates and compositions, under fouling and non-fouling conditions. Results compare favorably with measured data that show that solutions with similar total dissolved solids concentrations but different compositions result in different permeate fluxes. The model is then applied to predict permeate fluxes from the desalination of produced waters using a commercial spiral wound RO module. For NaCl-dominant brines, at total dissolved salt concentrations (TDS) below about 70 g/L, permeate fluxes are inversely proportional to water mole fraction as the latter is a reasonable approximation of water activity (i.e. ideal mixing). In the case of Ca–Cl-, Na–CO3- and Na–SO4-dominant brines below about 70 g/L TDS, this relationship does not hold as well and tends to overpredict osmotic pressure and thus underpredict permeate fluxes. However, the opposite becomes true at higher TDS values for typical produced waters. The scaling potential of these waters is also computed by allowing the precipitation of minerals above their saturation limit on the RO membrane. This work demonstrates how reactive transport models developed for the analysis of waters from geological systems can be extended to improve process design, optimization, and control in desalination systems from produced waters and beyond.

Molins, Sergi↗
DOE OSTI2022

Reactive Transport Modeling of Mineral Precipitation and Carbon Trapping in Discrete Fracture Networks

In this study we use numerical experiments to analyze reactive flow and transport behavior in discrete intersecting fracture networks, focusing on (a) how reaction-induced changes in physical and chemical properties affect flow connectivity and (b) how fracture networks developed in the Earth's critical zone contribute to carbon sequestration via mineral weathering reactions. In the first part of the study, we used two-dimensional reactive flow and transport simulations to analyze the impacts of mixing in a natural discrete fracture network. We concluded that reaction-induced changes can substantially alter the flow connectivity, especially at fracture intersections. The second set of simulations considered the problem of natural weathering of fractured mafic and ultramafic rocks in the partially saturated Earth's critical zone as a function of infiltration rates, fracture permeability, and partially saturated flow parameters. As a model system, we considered an incongruent reaction network with dissolution of forsterite and precipitation of magnesite. The behavior is complex in terms of the rate-controlling processes because of the multicomponent nature of the system as shown by the grid Peclet number: the CO 2 behavior is gas diffusion-controlled in the partially saturated zone, while the rate of water flow via the Damkӧhler number controls Mg 2+ transport through the fracture network. The amounts of carbon that can be trapped are modest, but the naturally fractured domain considered here provides a useful “base case” against which various engineered solutions can be compared.

58 GEOSCIENCES↗
DOE OSTI2022

Understanding the Hydrogeochemical Response of a Mountainous Watershed Using Integrated Surface‐Subsurface Flow and Reactive Transport Modeling

Abstract Climate change and other disturbances significantly impact hydrogeochemical exports from mountainous headwater catchments such as the Upper Colorado River Basin. Developing a mechanistic understanding of how the physical and chemical processes interact in time and space in an integrated manner is key to quantifying the future impacts of such disturbances. The hydrogeochemical response of a mountainous catchment in the 2010–2019 period is evaluated quantitatively using a high‐resolution model that simulates integrated hydrology, and transport and reactions for selected solutes and minerals. The model assumes that pyrite is present only at depth while calcite is distributed uniformly, and captures the observed C‐Q reasonably well. Distinct C‐Q dynamics are observed in an average (WY16), a wet (WY17), and a dry (WY18) water year. The model also quantifies the water fraction from surface, shallow and deep groundwater compartments using tracers, and suggests greater groundwater contributions to peak stream discharge in the dry WY18. Results demonstrate that calcium concentrations do not change significantly from year to year, while sulfate shows significant temporal variability. Pyrite dissolution is affected by the changing hydrological drivers where it is enhanced in the dry WY18; calcite dissolution supplements calcium dilution under high flow conditions. The model simulates the reaction hotspots controlled by hydrological conditions, and the spatially‐resolved results show that higher soil saturation and less snowpack occur earlier on the south‐facing side than on the north‐facing side. This is a first‐of‐its‐kind demonstration of a model that integrates hydrologic processes, including evapotranspiration, and reactive transport to enable a predictive understanding of hydrogeochemical exports.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Development of a molecularly informed biogeochemical framework for reactive transport modeling of subsurface carbon inventories, transformations and fluxes (Final Report)

The overall objective of the project was to combine new molecular-level characterization strategies with soil carbon flux measurements to develop and evaluate model representations of subsurface carbon cycling. We expanded on existing studies in the East River watershed, Colorado, in collaboration with the Berkeley Lab Watershed Function Scientific Focus Area (SFA), SLAC Groundwater Quality SFA, and Rocky Mountain Biological Laboratory (RMBL) to develop an elevation and vegetation gradient that is now the subject of long-term monitoring by the USGS. To achieve the overall objective, we combined field studies of soil respiration with laboratory analyses, ranging from spectroscopy to incubation studies. This combination of techniques enabled us to develop a new understanding of the drivers of high-elevation soil respiration. In the process, we developed three new modeling approaches that improve our ability to conceptualize and ultimately to represent soil respiration in numerical models. The first approach is a plot-scale transient inverse model that can be used to determine in situ CO 2 production rates from measured concentration profiles and surface fluxes. Application of this method revealed the importance of plant phenology and deep CO 2 production in moderating CO 2 fluxes to the atmosphere. The second modeling approach is a molecular-scale tool that enables spectroscopic and elemental data for carbon speciation to be transformed into functional group abundances, or the ‘SOC-fga model’. This method uniquely enables carbon speciation to be tracked within a reactive transport framework to partition carbon among different pathways and storage zones within the soil. The third modeling tool builds strongly on the previous approaches and captures the microbial processes driving heterotrophic respiration. This ‘dormancy model’ allows the native soil microbial population to respond transiently to the presence or absence of water in order to catalyze carbon respiration. This approach was also compared to the simpler and more widely used first-order model using two experimental datasets with different temporal and spatial resolutions. Our results illustrate that the simpler first order model provides a robust and efficient representation of deep (>1 meter) soil respiration, but that shallow soils, where most respiration occurs, require explicit representation of moisture-dependent activation and dormancy rates. By assessing soil organic carbon turnover at multiple scales, we see a complex array of controls emerge. Atte scale of a hillslope, spatial heterogeneity in soil respiration rates dominates and is uncorrelated with instantaneous soil moisture and plant community. At the profile scale, the balance between plant inputs and water availability is the dominant control. At the microbial to molecular scale, physiological processes associated with carbon use and carbon speciation are important controls. These scale-dependent controls emphasize the need for new modeling approaches that examine their interactions and hierarchies.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Investigation of Coupled Processes in Fractures and the Bordering Matrix via a Micro‐Continuum Reactive Transport Model

Abstract In multi‐mineral fractured rocks, the altered porous layer on the fracture surface resulting from preferential dissolution of the fast‐reacting minerals can have profound impacts on subsequent chemical‐physical alteration of the fractures. This study adopts the micro‐continuum approach to provide further understanding of reactive transport processes in the altered layer (AL), and mass exchanges with the bordering matrix and fracture. The modeling framework couples the Darcy‐Brinkman‐Stokes (DBS) solver in COMSOL Multiphysics and the geochemical modeling capability of CrunchFlow. Three‐dimensional steady state simulations with systematically varied chemical‐physical parameters of the AL were performed to examine the impacts of individual factors and processes. Our simulation results confirm previous observations that dissolution of the fast‐reacting mineral (i.e., calcite) is largely controlled by diffusion across the AL. We also show that dissolution of the slow‐reacting mineral (i.e., dolomite), which controls AL development and fracture enlargement, increases with surface area and has a complex dependence on different local rate‐limiting processes. In particular, advection can result in evident spatial variations in the local dissolution rates of dolomite, although it does not affect the bulk chemistry significantly. The difference in the spatial patterns between simulations with and without advection in the AL is more noticeable in the locations with smaller apertures, with up to 20% difference in local reaction rates. Therefore, it is important to include a full depiction of advection, diffusion, and reactions for accurately capturing local dynamics that control long‐term fracture evolution.

58 GEOSCIENCES↗