Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reactive transport modeling”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2012

Coupling in silico microbial models with reactive transport models to predict the fate of contaminants in the subsurface (Final Report)

This project successfully accomplished its goal of coupling genome-scale metabolic models with hydrological and geochemical models to predict the activity of subsurface microorganisms during uranium bioremediation. Furthermore, it was demonstrated how this modeling approach can be used to develop new strategies to optimize bioremediation. The approach of coupling genome-scale metabolic models with reactive transport modeling is now well enough established that it has been adopted by other DOE investigators studying uranium bioremediation. Furthermore, the basic principles developed during our studies will be applicable to much broader investigations of microbial activities, not only for other types of bioremediation, but microbial metabolism in diversity of environments. This approach has the potential to make an important contribution to predicting the impact of environmental perturbations on the cycling of carbon and other biogeochemical cycles.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Reactive transport modeling of microbial dynamics

Reactive transport modeling of microbially mediated processes has contributed significantly to an improved understanding of elemental cycling in Earth near-surface environments. We describe key characteristics of microbial reactive transport models (RTMs), recent advances in modeling approaches, and their application to terrestrial and marine environmental problems. We introduce relevant case studies and discuss ways to integrate -omics data that can inform and validate microbial RTMs, improving our ability to address grand challenges in a changing world.

Meile, Christof↗
DOE OSTI2015

Maximum likelihood Bayesian model averaging and its predictive analysis for groundwater reactive transport models

While Bayesian model averaging (BMA) has been widely used in groundwater modeling, it is infrequently applied to groundwater reactive transport modeling because of multiple sources of uncertainty in the coupled hydrogeochemical processes and because of the long execution time of each model run. To resolve these problems, this study analyzed different levels of uncertainty in a hierarchical way, and used the maximum likelihood version of BMA, i.e., MLBMA, to improve the computational efficiency. Our study demonstrates the applicability of MLBMA to groundwater reactive transport modeling in a synthetic case in which twenty-seven reactive transport models were designed to predict the reactive transport of hexavalent uranium (U(VI)) based on observations at a former uranium mill site near Naturita, CO. Moreover, these reactive transport models contain three uncertain model components, i.e., parameterization of hydraulic conductivity, configuration of model boundary, and surface complexation reactions that simulate U(VI) adsorption. These uncertain model components were aggregated into the alternative models by integrating a hierarchical structure into MLBMA. The modeling results of the individual models and MLBMA were analyzed to investigate their predictive performance. The predictive logscore results show that MLBMA generally outperforms the best model, suggesting that using MLBMA is a sound strategy to achieve more robust model predictions relative to a single model. MLBMA works best when the alternative models are structurally distinct and have diverse model predictions. When correlation in model structure exists, two strategies were used to improve predictive performance by retaining structurally distinct models or assigning smaller prior model probabilities to correlated models. Since the synthetic models were designed using data from the Naturita site, the results of this study are expected to provide guidance for real-world modeling. Finally, limitations of applying MLBMA to the synthetic study and future real-world modeling are discussed.

58 GEOSCIENCES↗
DOE OSTI2016

Reactive Transport Models with Geomechanics to Mitigate Risks of CO2 Utilization and Storage

Reactivity of carbon dioxide (CO 2 ), rocks and brine is important in a number of practical situations in carbon dioxide sequestration. Injectivity of CO 2 will be affected by near wellbore dissolution or precipitation. Natural fractures or faults containing specific minerals may reactivate leading to induced seismicity. In this project, we first examined if the reactions between CO 2 , brine and rocks affect the nature of the porous medium and properties including petrophysical properties in the timeframe of the injection operations. This was done by carrying out experiments at sequestration conditions (2000 psi for corefloods and 2400 psi for batch experiments, and 600°C) with three different types of rocks – sandstone, limestone and dolomite. Experiments were performed in batch mode and corefloods were conducted over a two-week period. Batch experiments were performed with samples of differing surface area to understand the impact of surface area on overall reaction rates. Toughreact, a reactive transport model was used to interpret and understand the experimental results. The role of iron in dissolution and precipitation reactions was observed to be significant. Iron containing minerals – siderite and ankerite dissolved resulting in changes in porosity and permeability. Corefloods and batch experiments revealed similar patterns. With the right cationic balance, there is a possibility of precipitation of iron bearing carbonates. The results indicate that during injection operations mineralogical changes may lead to injectivity enhancements near the wellbore and petrophysical changes elsewhere in the system. Limestone and dolomite cores showed consistent dissolution at the entrance of the core. The dissolution led to formation of wormholes and interconnected dissolution zones. Results indicate that near wellbore dissolution in these rock-types may lead to rock failure. Micro-CT images of the cores before and after the experiments revealed that an initial high-permeability pathway facilitated the formation of wormholes. The peak cation concentrations and general trends were matched using Toughreact. Batch reactor modeling showed that the geometric factors obtained using powder data that related effective surface area to the BET surface area had to be reduced for fractured samples and cores. This indicates that the available surface area in consolidated samples is lower than that deduced from powder experiments. Field-scale modeling of reactive transport and geomechanics was developed in parallel at Idaho National Laboratory. The model is able to take into account complex chemistry, and consider interactions of natural fractures and faults. Poroelastic geomechanical considerations are also included in the model.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2012

Reactive Transport Modeling

This special section in the Vadose Zone Journal focusing on reactive transport modeling was developed from a special symposium jointly sponsored by the Soil Physics and Soil Chemistry Divisions of the Soil Science Society of America at the 2010 annual meetings held in Long Beach, CA. It contains eight contributions reflecting a diverse range of reactive transport modeling applications that address the fate and transport of both reactive solutes and colloidal material (i.e., viruses). Further, the eight studies range in scale from laboratory batch and column experiments to the regional scale. Three of the submissions reflect efforts to model experimental data that were already reported in the literature, providing additional insight concerning underlying physical and chemical processes. In addition to the obvious predictive applications, the eight contributions to this special section illustrate the additional advantages to including reactive transport modeling as an active component within a diverse research agenda. In addition to scenario testing, reactive transport modeling provides a framework for synthesizing diverse experimental data, evaluating the relative importance of various competing and interrelated biogeochemical processes, and identifying important data gaps and research priorities under more natural conditions at multiple temporal and spatial scales.

Seaman, J. C.↗
DOE OSTI2022

Integration of Omics into a New Comprehensive Rate Law for Competitive Terminal Electron-Accepting Processes in Reactive Transport Models: Application to N, Fe, S, and Contaminant Transformations in Stream and Wetland Sediments

Surface waters represent important sources of alternate energy and drinking water in the United States, and characterizing the biogeochemical processes that affect surface water quality is relevant to the DOE-BER mission. Sediment biogeochemical processes regulate the release of carbon (C), nutrients, and contaminants to surface waters and thus influence water quality. Sediment biogeochemical processes are dynamic and affected by the deposition and remobilization of solid material and changes in environmental conditions driven by water discharge variations. Wetlands are important natural filters of surface waters which may either trap, metabolize, or mobilize nutrients and contaminants. Despite their importance, biogeochemical processes regulating nutrient and contaminant release and C transformation in stream and wetland sediments cannot be predicted accurately by current mathematical models. These reactive transport models largely rely on detectable changes in geochemical conditions to activate metabolic processes, do not accurately account for the competition between microbial processes, and poorly constrain effects of hydrological perturbations on biogeochemical processes. In this BER-SBR exploratory project, metagenomic and geochemical signatures were combined to identify microbially-mediated redox processes in anaerobic stream and wetland sediments from the Savannah River Site (SRS, ANL SFA) and East Fork Poplar Creek (EFPC, ORNL SFA) that play important roles in C, uranium (U), and mercury (Hg) transformations. In addition, sediment incubations were conducted to examine the competition between anaerobic respiration processes Finally, new rate laws were developed for reactive transport models that rely on complementary metagenomic and geochemical signatures to identify the underlying anaerobic microbial processes in stream and wetland sediments, describe the competition between the dominant metabolic processes involved in nutrient release and U and Hg mobilization, and more accurately quantify carbon transformation and the response of microbial processes to changes in redox conditions associated with hydrological forcing. These rate laws were optimized in batch reactors with SRS wetland sediments, where iron and sulfate reduction dominate. Anaerobic carbon remineralization processes followed the expected thermodynamic sequence of microbial respiration with depth in both sediments, except that geochemical signals indicated that sulfate reduction was inactive in EFPC sediments and moderate in SRS wetland sediments. Estimates indicated that microbial iron reduction contributed to at least half of the production of reduced iron in these sediments. Incubations demonstrated that nitrate reduction, denitrification, and dissimilatory nitrate reduction to ammonium were active in the natural EFPC sediment and activated upon nitrate amendment in these nitrate-rich sediments. In turn, these processes were outcompeted by the addition of either iron oxides or sulfate as alternative terminal electron acceptors. Although geochemical products of sulfate reduction were not detected in the incubations, the abundance of sulfate reduction genes increased with depth in the sediment and was equally more pronounced in treatments amended with either iron oxides or sulfate. Simultaneously, anaerobic sulfide oxidizing bacteria coupling sulfide oxidation to DNRA (and not conventional denitrification) were apparently enriched over time, regardless of the treatments. These findings indicate that sulfate reduction is important in freshwater stream sediments and probably catalyzed by a cryptic sulfur cycle involving nitrogen species, in which the sulfur products from sulfate reduction are immediately removed by side reactions and not detectable by geochemical measurements alone. Similar experiments in SRS sediments, however, demonstrated little interaction between nitrogen and sulfur cycling microorganisms. Sulfate reduction was impacted by the addition of more thermodynamically favorable electron acceptors, suggesting either that iron-reducing microorganisms outcompeted sulfate-reducing microorganisms for organic substrate, depleting the stock of electron donor available for sulfate reduction, or that the cryptic sulfur cycle was shunted by the precipitation of FeS generated as a result of the abiotic reduction of iron oxides by dissolved sulfide. A diagnostic modeling exercise was conducted to further investigate the competition between terminal electron accepting processes. As conventional kinetic models typically do not account for cryptic cycles and used inaccurate formulations to describe competition between microbial communities, new metabolic rate laws were developed that explicitly express the electron acceptor-specific enzyme of each energetically favorable metabolic process based on gene abundance detected in the incubations. The model was tested with the sediment slurry incubation data to determine whether substrate competition could explain the decrease in sulfate reduction observed in the presence of iron oxide competitor. The model was able to reproduce geochemical concentrations really well in each treatment once the model was calibrated with the unamended control, suggesting that microbial competition was indeed driven by thermodynamic considerations. Overall, carbon remineralization processes and rates will be reproduced much more realistically with the new metabolic rate laws.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

High-Performance Simulation of Dynamic Hydrologic Exchange and Implications for Surrogate Flow and Reactive Transport Modeling in a Large River Corridor

Hydrologic exchange flows (HEFs) have environmental significance in riverine ecosystems. Key river channel factors that influence the spatial and temporal variations of HEFs include river stage, riverbed morphology, and riverbed hydraulic conductivity. However, their impacts on HEFs were often evaluated independently or on small scales. In this study, we numerically evaluated the combined interactions of these factors on HEFs using a high-performance simulator, PFLOTRAN, for subsurface flow and transport. The model covers 51 square kilometers of a selected river corridor with large sinuosity along the Hanford Reach of the Columbia River in Washington, US. Three years of spatially distributed hourly river stages were applied to the riverbed. Compared to the simulation when riverbed heterogeneity is not ignored, the simulation using homogeneous riverbed conductivity underestimated HEFs, especially upwelling from lateral features, and overestimated the mean residence times derived from particle tracking. To derive a surrogate model for the river corridor, we amended the widely used transient storage model (TSM) for riverine solute study at reach scale with reactions. By treating the whole river corridor as a batch reactor, the temporal changes in the exchange rate coefficient for the TSM were derived from the dynamic residence time estimated from the hourly PFLOTRAN results. The TSM results were evaluated against the effective concentrations in the hyporheic zone calculated from the PFLOTRAN simulations. Our results show that there is potential to parameterize surrogate models such as TSM amended with biogeochemical reactions while incorporating small-scale process understandings and the signature of time-varying streamflow to advance the mechanistic understanding of river corridor processes at reach to watershed scales. However, the assumption of a well-mixed storage zone for TSM should be revisited when redox-sensitive reactions in the storage zones play important roles in river corridor functioning.

Fang, Yilin↗
DOE OSTI2021

Local and Global Sensitivity Analysis of a Reactive Transport Model Simulating Floodplain Redox Cycling

Reactive transport models (RTMs) are essential tools that simulate the coupling of advective, diffusive, and reactive processes in the subsurface, but their complexity makes them difficult to understand, develop and improve without accompanying statistical analyses. Although global sensitivity analysis (SA) can address these issues, the computational cost associated with most global SA techniques limits their use with RTMs. In this study, we apply distance-based generalized sensitivity analysis (DGSA), a novel and computationally efficient method of global SA, to a floodplain-scale RTM and compare DGSA results to those from local SA. Our test case focuses on the impact of 17 uncertain environmental parameters on spatially and temporally variable redox conditions within a floodplain aquifer. The input parameters considered include flow and diffusion rates, geochemical reaction rates, and the spatial distribution of sediment facies. Sensitivity was evaluated for three distinct components of the model response, encompassing both multidimensional and categorical output. Parameter rankings differ between local SA and DGSA, due to nonlinear effects of individual parameters and interaction effects between parameters. DGSA results show that fluid residence time, which is controlled by aquifer permeability, generally exerts a stronger control on redox conditions than do geochemical reaction rates. Sensitivity indices also demonstrate that sulfate reduction is key for establishing and maintaining reducing conditions throughout the aquifer. Furthermore, these results provide insights into the key drivers of heterogeneous redox processes within floodplain aquifers, as well as the main sources of uncertainty when modeling complex subsurface systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

BioRT‐HBV 1.0: A Biogeochemical Reactive Transport Model at the Watershed Scale

Abstract Reactive Transport Models (RTMs) are essential tools for understanding and predicting intertwined ecohydrological and biogeochemical processes on land and in rivers. While traditional RTMs have focused primarily on subsurface processes, recent watershed‐scale RTMs have integrated ecohydrological and biogeochemical interactions between surface and subsurface. These emergent, watershed‐scale RTMs are often spatially explicit and require extensive data, computational power, and computational expertise. There is however a pressing need to create parsimonious models that require minimal data and are accessible to scientists with limited computational background. To that end, we have developed BioRT‐HBV 1.0, a watershed‐scale, hydro‐biogeochemical RTM that builds upon the widely used, bucket‐type HBV model known for its simplicity and minimal data requirements. BioRT‐HBV uses the conceptual structure and hydrology output of HBV to simulate processes including advective solute transport and biogeochemical reactions that depend on reaction thermodynamics and kinetics. These reactions include, for example, chemical weathering, soil respiration, and nutrient transformation. The model uses time series of weather (air temperature, precipitation, and potential evapotranspiration) and initial biogeochemical conditions of subsurface water, soils, and rocks as input, and output times series of reaction rates and solute concentrations in subsurface waters and rivers. This paper presents the model structure and governing equations and demonstrates its utility with examples simulating carbon and nitrogen processes in a headwater catchment. As shown in the examples, BioRT‐HBV can be used to illuminate the dynamics of biogeochemical reactions in the invisible, arduous‐to‐measure subsurface, and their influence on the observed stream or river chemistry and solute export. With its parsimonious structure and easy‐to‐use graphical user interface, BioRT‐HBV can be a useful research tool for users without in‐depth computational training. It can additionally serve as an educational tool that promotes pollination of ideas across disciplines and foster a diverse, equal, and inclusive user community.

Sadayappan, Kayalvizhi↗
DOE OSTI2019

Tracking Diverse Minerals, Hungry Organisms, and Dangerous Contaminants Using Reactive Transport Models

Beneath our feet is a fascinating world of flowing water, cosmopolitan microbes, and complex mineral assemblages. Yet we see none of it from above. Our quest to investigate these complex subsurface interactions has led to the development of reactive transport models. These are computer algorithms that allow us to explore, in a virtual way, the natural dynamics of Earth's systems and our anthropogenic impact on those systems. Here, we explain the concepts behind reactive transport models—which include the transport of aqueous species and the descriptions of biogeochemical reactions involving solutes, surfaces and microorganisms—and introduce to reactive transport applications in terrestrial and marine environments.

58 GEOSCIENCES↗
DOE OSTI2013

A Dual Regime Reactive Transport Model for Simulation of High Level Waste Tank Closure Scenarios - 13375

A numerical simulation framework is presented in this paper for estimating evolution of pH and release of major species from grout within high-level waste tanks after closure. This model was developed as part of the Cementitious Barriers Partnership. The reactive transport model consists of two parts - (1) transport of species, and (2) chemical reactions. The closure grout can be assumed to have varying extents of cracking and composition for performance assessment purposes. The partially or completely degraded grouted tank is idealized as a dual regime system comprising of a mobile region having solid materials with cracks and macro-pores, and an immobile/stagnant region having solid matrix with micropores. The transport profiles of the species are calculated by incorporating advection of species through the mobile region, diffusion of species through the immobile/stagnant region, and exchange of species between the mobile and immobile regions. A geochemical speciation code in conjunction with the pH dependent test data for a grout material is used to obtain a mineral set that best describes the trends in the test data of the major species. The dual regime reactive transport model predictions are compared with the release data from an up-flow column percolation test. The coupled model is then used to assess effects of crack state of the structure, rate and composition of the infiltrating water on the pH evolution at the grout-waste interface. The coupled reactive transport model developed in this work can be used as part of the performance assessment process for evaluating potential risks from leaching of a cracked tank containing elements of human health and environmental concern. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2016

Impacts of elevated dissolved CO 2 on a shallow groundwater system: Reactive transport modeling of a controlled-release field test

We report one of the risks that CO 2 geological sequestration imposes on the environment is the impact of potential CO 2 /brine leakage on shallow groundwater. The reliability of reactive transport models predicting the response of groundwater to CO 2 leakage depends on a thorough understanding of the relevant chemical processes and key parameters affecting dissolved CO 2 transport and reaction. Such understanding can be provided by targeted field tests integrated with reactive transport modeling. A controlled-release field experiment was conducted in Mississippi to study the CO 2 -induced geochemical changes in a shallow sandy aquifer at about 50 m depth. The field test involved a dipole system in which the groundwater was pumped from one well, saturated with CO 2 at the pressure corresponding to the hydraulic pressure of the aquifer, and then re-injected into the same aquifer using a second well. Groundwater samples were collected for chemical analyses from four monitoring wells before, during and after the dissolved CO 2 was injected. In this paper, we present reactive transport models used to interpret the observed changes in metal concentrations in these groundwater samples. A reasonable agreement between simulated and measured concentrations indicates that the chemical response in the aquifer can be interpreted using a conceptual model that encompasses two main features: (a) a fast-reacting but limited pool of reactive minerals that responds quickly to changes in pH and causes a pulse-like concentration change, and (b) a slow-reacting but essentially unlimited mineral pool that yields rising metal concentrations upon decreased groundwater velocities after pumping and injection stopped. During the injection, calcite dissolution and Ca-driven cation exchange reactions contribute to a sharp pulse in concentrations of Ca, Ba, Mg, Mn, K, Li, Na and Sr, whereas desorption reactions control a similar increase in Fe concentrations. Finally, after the injection and pumping stops and the groundwater flow rate decreases, the dissolution of relatively slow reacting minerals such as plagioclase drives the rising concentrations of alkali and alkaline earth metals observed at later stages of the test, whereas the dissolution of amorphous iron sulfide causes slowly increasing Fe concentrations.

58 GEOSCIENCES↗
DOE OSTI2026

Hyporheic‐Zone Processes and Stream Oxygen Dynamics: Insights From a Multiscale Reactive Transport Model

Aquatic ecosystem metabolism encapsulates the daily fixation (gross primary production, GPP d ) and mineralization (ecosystem respiration, ER d ) of organic carbon. In fluvial systems, these are commonly estimated by inverse solutions to field observations using a model that describes oxygen concentrations varying in the water column in response to metabolic fluxes and air‐water gas exchange controlled by a rate coefficient (K 600 ). The most common conceptual model is the single‐station metabolism (SSM) model. The simplicity and flexibility of this conceptualization make it attractive; however, it implicitly assumes that all the processes that consume oxygen in fluvial systems can be lumped into a bulk estimate of respiration with poorly understood consequences for estimates of GPP d , ER d , and K 600 . Here, we focus on the implications of using SSM conceptualization when estimating metabolic fluxes from oxygen dynamics in channels where hyporheic exchange occurs. We use a new multiscale numerical model for reactive transport in streams that represents hyporheic exchange and streambed heterotrophic respiration. Nondimensionalization of this model reveals dimensionless groups that collectively control oxygen dynamics. Numerical experiments offer a mechanistic understanding of the impacts of hyporheic exchange on diel oxygen dynamics revealing that potential biases arise from neglecting mass transfer limitations. Specifically, we found that hyporheic exchange significantly affects diel oxygen dynamics, even for nonreactive streambed sediments. Moreover, while the SSM performs well in many situations, we find conditions where significant bias is produced by hyporheic exchange, even when oxygen data are well‐fitted. These situations pose a major challenge in the interpretation of metabolism assessment estimates.

Gomez‐Velez, Jesus D. [Oak Ridge National Laborato↗