Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Model calibration”

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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Informing electrification strategies of residential neighborhoods with urban building energy modeling

Electrifying end uses is a key strategy to reducing GHG emissions in buildings. However, it may increase peak electricity demand that triggers the need to upgrade the existing power distribution system, leading to delays in electrification and needs of significant investment. There is also concern that building electrification may cause an increase of energy costs, leading to further energy burden for low-income communities. This study uses the urban scale building modeling tool CityBES to assess the electrification impacts of more than 43,000 residential buildings in a neighborhood of Portland, Oregon, USA. Energy efficiency upgrades were investigated on their potential to mitigate the increase of peak electricity demand and energy burden. Simulation results from the calibrated EnergyPlus models show that electrification with heat pumps for space heating and cooling as well as for domestic water heating can reduce CO2e emissions by 38%, but increase peak electricity demand by about 9% from the baseline building stock. Combining electrification measures and energy efficiency upgrades can reduce CO2e emissions by 48% while reducing peak electricity demand by 6% and saving the median household energy costs by 28%. City and utility decision makers should consider integrating energy efficiency upgrades with electrification measures as an effective residential building electrification strategy, which significantly reduces carbon emissions, caps or even decreases peak demand while reducing energy burden of residents.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

Dynamic flow stress of pure polycrystalline aluminum: Pressure-shear plate impact experiments and extension of dislocation-based modeling to large strains

The dynamic thermo-mechanical behavior of pure aluminum has attracted renewed interest lately due to experimental observations of an anomalous increase in Hugoniot Elastic Limit (HEL) at incipient plasticity and elevated temperatures in polycrystalline pure metals. In context of current dislocation-mediated plasticity models for metals, this increase in dynamic strength is indicative of a transition in the rate controlling mechanism for dislocation glide from being thermally assisted to phonon-drag restricted due to increase in phonon viscosity at elevated temperatures. Though these studies have helped to shed light on these important mechanisms operative in FCC metals at incipient plasticity, the extent to which they contribute to flow stress, particularly at larger strains, remains unclear. Here, we address these questions through a combined experimental and modeling effort focused on investigating the evolution of dynamic flow stress in polycrystalline aluminum using experimental data gathered from a series of combined pressure-and-shear plate impact (PSPI) experiments designed to reveal the flow stress of pure aluminum at strain rates ~ 10 5 /s, plastic strains of up to 40% and temperatures ranging from room to 866 K. In all cases, the flow stress of aluminum, as inferred from the measured transverse particle velocity histories at the free surface of an elastic tungsten carbide target plate, reveals saturation with increasing plastic strains at stress levels that decrease with increasing test temperatures. Numerical simulations are performed to correlate the experimentally observed temperature and strain rate dependence of flow stress at small and large plastic strains using the Austin-McDowell dislocation-mediated plasticity model parametrized to normal plate impact experiments conducted in an earlier study by Zaretsky and Kanel (2012). Extensions to the model are made to better represent dynamic behavior of pure aluminum at larger plastic strains as observed in elevated temperature split Hopkinson Pressure bar (SPHB) experiments of Samanta (1971) and Lindholm and Yeakly (1965), and the combined pressure-and-shear plate-impact experiments conducted in the present study. The main theoretical extension to the Austin-McDowell model made in this paper is the introduction of a new rate- and temperature- dependent dynamic recovery function, which can potentially allow for an effective reduction in the rate of accumulation of dislocations at large plastic strains. The numerical predictions of the revised and re-calibrated plasticity model are brought into agreement with the experimental observations, correlating sufficiently well with dynamic yield stress at incipient plasticity and the flow stress levels at larger plastic strains, plastic strain rates in the range 10 3 – 10 6 /s, and elevated temperatures up to near melt. The model, in concert with the experimental measurements, suggests that at incipient plasticity the rate governing mechanism for plastic flow is phonon-drag restricted dislocation glide, whereas, at higher magnitudes of plastic strain, it transitions to stress-assisted thermally activated glide. The transition between these two rate governing mechanisms is controlled by the evolution of dislocations throughout the deformation process.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Modeling Assessment of Tidal Energy Extraction in the Western Passage

Numerical models have been widely used for the resource characterization and assessment of tidal instream energy. The accurate assessment of tidal stream energy resources at a feasibility or project-design scale requires detailed hydrodynamic model simulations or high-quality field measurements. This study applied a three-dimensional finite-volume community ocean model (FVCOM) to simulate the tidal hydrodynamics in the Passamaquoddy–Cobscook Bay archipelago, with a focus on the Western Passage, to assist tidal energy resource assessment. IEC Technical specifications were considered in the model configurations and simulations. The model was calibrated and validated with field measurements. Energy fluxes and power densities along selected cross sections were calculated to evaluate the feasibility of the tidal energy development at several hotspots that feature strong currents. When taking both the high current speed and water depth into account, the model results showed that the Western Passage has great potential for the deployment of tidal energy farms. The maximum extractable power in the Western Passage was estimated using the Garrett and Cummins method. Different criteria and methods recommended by the IEC for resource characterization were evaluated and discussed using a sensitivity analysis of energy extraction for a hypothetical tidal turbine farm in the Western Passage.

16 TIDAL AND WAVE POWER↗
DOE OSTI2023

Battery state-of-health diagnostics during fast cycling using physics-informed deep-learning

Rapid, in-situ Li-ion battery state-of-health (SOH) quantification is challenging. Li-ion battery aging can vary significantly with chemistry, operating conditions, cycling demands, electrode design, and operation history. As a cell ages, optimal and safe operating conditions need to be adapted to account for battery degradation by tracking critical aging modes such as loss-of-lithium-inventory (LLI), loss-of-active-material (LAM) in either electrode, and/or impedance rise. This manuscript describes a framework for identifying battery aging modes in-operando using fast-rate voltage charge/discharge responses. The framework uses a physically based Li-ion battery model to produce synthetic high-rate responses at aged states. The aging model is calibrated against experimental data from cells with different electrode loadings and cycled under a variety of fast-charging conditions (1 h, 15 min, 10 min, and 7 min charging). The synthetically generated high-rate responses at aged states are then used to train a deep-learning model to identify real cell state-of-health from fast charge/discharge battery voltage responses. The synthetically trained deep-learning model performance is validated by comparing to standard incremental capacity analysis and half-cell measurements. Finally, the framework demonstrates the benefits of using high-rate physics-based models to generate synthetic data for training deep-learning models.

25 ENERGY STORAGE↗
DOE OSTI2022

Modeling deformation and failure in AlSi-polyester abradable sealcoating material using microstructure-based finite element simulation

A plasma-sprayed aluminum-silicon (AlSi)/polyester coating is applied in modern gas turbine engines as an abradable sealcoating to maintain tight clearances between the rotating blades and the static casing. While running the engine, the rotating blades “rub” with the abradable coating, which results in extreme strain rate (up to 10 6 ) dynamics and a high-temperature environment. Due to the difficulty of collecting direct measurements, predictive computational models are important for analyzing the deformation and failure of the abradable material, to help meet the design target of avoiding damage to the blade tip and maintaining high fuel efficiency. In this research, a microstructure-based finite element (FE) computational model was developed to capture the complex mechanical behavior of the AlSi/polyester microstructure. The model is based on a virtual representative-volume-element (RVE) of a metal-polymer microstructure, reconstructed from x-ray computed tomography. It models the plastic deformation of, and damage to, each AlSi and polyester constituents, as well as the failure at their interface. The model was calibrated and validated with uniaxial tension and compression experiments, conducted at two temperatures (298 K and 533 K) at an applied strain rate of 10 3 - 10 4 s -1 . The material exhibited strongly asymmetric tension-compression behavior and a sensitivity to temperature, which was well captured by the model. The model was further applied to investigate changes in mechanical behavior due to variations in constituents’ volume fractions, which provides guidance to the microstructural design of AlSi/polyester abradable materials. The model is expected to facilitate the development of improved abradable materials by bypassing the conventional trial-and-error approach and extensive testing requirements.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Three-Dimensional Evaluation of Sand Particle Fracture Using Discrete-Element Method and Synchrotron Microcomputed Tomography Images

Recent research showed that fracture of sand particles plays a significant role in determining the plastic bulk volumetric changes of granular materials under different loading conditions. One of the major tools used to better understand the influence of particle fracture on the behavior of granular materials is discrete-element modeling (DEM). This paper employed the bonded block model (BBM) to simulate the fracture behavior of sand. Each sand particle is modeled as an agglomerate of rigid blocks bonded at their contacts using the linear-parallel contact model, which can transmit both moment and force. DEM simulated particles closely matched the actual three-dimensional (3D) shape of sand particles acquired using high-resolution 3D synchrotron microcomputed tomography (SMT). Results from unconfined one-dimensional (1D) compression of a single synthetic silica cube were used to calibrate the model parameters. Particle fracture was investigated for specimens composed of three sand particles that were loaded under confined 1D compression. Breakage energy measured from DEM models matched well with that measured experimentally. The paper studied the effects of contact loading condition and particle interaction on the fracture mode of particles using BBM that can closely capture the 3D shape of real sand particles.

58 GEOSCIENCES↗
DOE OSTI2025

A comparative analysis of numerical approaches for the description of gas flow in clay-based repository systems: From a laboratory to a large-scale gas injection test

There is nowadays a consensus among many countries that geological disposal is a favourable solution for the long-term management. Although different host formations and different barrier systems are under consideration around the world, clay-based materials form an important component for waste isolation in most national programmes. Hence, a good comprehension of the effect of gas flow on the hydro-mechanical behaviour of clay-based soils is essential, both at laboratory and field scale. Task B under the international cooperative project DECOVALEX-2023 has recently shown that, after some enhancement, models can be employed to reproduce laboratory scale tests, even with different sample geometries 37 . However, further work is required to understand whether they can be applied to simulate a large-scale experiment. Up-scaling of models for the advective transport of gas through clay-based low permeable material presents a number of problems related to the difficulty in obtaining consistent hydrogeological parameters and constitutive relationships at both laboratory and field scale. Based on a unique dataset from a large-scale gas injection test (Lasgit) performed at the Äspö Hard Rock Laboratory (Sweden), Task B within DECOVALEX-2023 has explored the refinement of these numerical strategies applied to the simulation of gas flow. Work performed within the task reveals that codes do not need to be substantially modified from the laboratory models to reproduce full-scale tests: indeed, model parameters calibrated and validated at laboratory scale have been applied to predict field scale gas flow at Lasgit, including peak gas pressure and injected cumulative gas volume. By means of (1) the introduction of interfaces between blocks to reflect the experimental configuration and the (2) adjustment of some parameters (e.g., higher permeability), the updated models are able to represent most of the key features observed in the experimental data, even at a large scale.

Tamayo-Mas, E↗
DOE OSTI2025

Complex Dependence of Calcite Crack Kinetics on Salinity: The Role of DLVO and Hydration Forces

Abstract Subcritical crack growth (SCG) plays an important role in many geological processes such as delayed earth rupture and rock weathering. The complex dependency of SCG on the in‐crack fluid chemistry, however, is still poorly understood. In this study, we utilize the newly developed surface force‐based fracture theory (SFFT) to elucidate the relative contributions of surface forces and solute transport to the crack growth kinetics of calcite in NaCl solutions. Expanding on Barenblatt's cohesive crack model, SFFT introduces an effective stress intensity at the crack tip that encompasses all the relevant intermolecular forces across the crack in addition to the external far‐field stresses. The nonlinear system of equations portraying the crack opening profile, the solute distribution in a propagating crack, and the crack growth velocity are numerically solved via an implicit scheme. After carefully calibrating the model for calcite‐water systems, the SFFT is used to predict the SCG response of calcite at different NaCl concentrations, based on various hypotheses. These predictions are then compared to existing SCG data from the literature. We demonstrate that the experimentally observed variation of SCG rate with NaCl concentration cannot be explained solely by DLVO forces (electrostatic and Van der Waals interactions). This can be remediated by introducing an exponentially decaying hydration force with a nonlinear, nonmonotonic dependence on NaCl concentration. Furthermore, we demonstrate that accounting for both diffusive and advective transport of ions is important in explaining the absence of a stage‐II SCG response for calcite in electrolyte solutions. Plain Language Summary Subcritical crack growth (SCG) refers to the slow propagation of cracks in materials under a stress below the threshold for catastrophic failure. SCG is a key process in many geological events, for example, delayed earth ruptures and rock weathering. New initiatives such as underground CO 2 and H 2 storage in carbonate reservoirs further call for better understanding of SCG in carbonate minerals subjected to varying fluid chemistry. This study examines the SCG of calcite, a key mineral found in carbonate rocks, intergranular cement in sandstones, and filling material in mineral veins and faults, determining their deformation and strength. A mathematical model is developed to describe how the crack opens and propagates, how solutes (like salts) distribute within the crack, and how the crack surfaces interact with each other. We used the model to predict calcite SCG in water at different salt concentrations and compared it with experimental data. Our results revealed that the hydration force is the dominating factor in determining the complex, non‐linear dependency of SCG on salinity. We also found that both the movement of ions by diffusion and by bulk water flow are crucial for explaining the SCG rates, especially when the cracks grow quickly. Key Points Surface Force‐Based Fracture Theory predicts the complex subcritical crack growth patterns of calcite crystals immersed in NaCl solutions Results highlight the dominant role of hydration forces in altering the fracture behavior of calcite compared to VdW and electric double‐layer forces Advective solute transport explains the absence of stages‐II and ‐III subcritical crack growth responses in solid‐liquid systems

DLVO↗
DOE OSTI2025

Heat Transfer in Void Generating Foam Decomposition: Further Development

Continued development of the additive conductivity material model, used to simulate changes in heat transfer that occurs in void generating foam decomposition, has resulted in an improved model and new features. The previous version of the model was calibrated against the Aria Bulk Fluid Element (BFE) solution and proposed a third-order polynomial correction term best captured the increased heat transfer due to voids in the foam. An investigation of the Fuego Conjugate Heat Transfer (CHT) and Aria BFE solutions at several geometries revealed the CHT solution and BFE solution had differing behavior across length scales, especially at smaller scales. Five calibration studies, using the Fuego CHT as the calibration data, were carried out with polynomial functions of 4-th, 3-rd, 2-nd, 1-st and 0-th orders to determine the best correction function that generalized well across length scales. Each polynomial function was calibrated/trained on six different sized geometries and then tested on three uniquely sized geometries. This study revealed that the 1-st order additive conductivity model performed the best. A new feature of void formation scaling was implemented to more realistically capture the heat transfer as voids are created. A scaling term was added to the model to activate the conductivity correction as decomposition progresses.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Ectomycorrhizal effects on decomposition are highly dependent on fungal traits, climate, and litter properties: A model-based assessment. Dataset.

To simulate the effects of mycorrhizal fungi on soil organic matter cycling, we incorporated mycorrhizal processes into the Carbon, Organisms, Rhizosphere, and Protection in the Soil Environment (CORPSE) model to develop a new soil model Myco-CORPSE. The new model was calibrated and evaluated against soil measurements taken at temperate forests in New Hampshire (NH) and Georgia (GA). A series of scenario analysis were also conducted to explore the conditions under which ectomycorrhizal (ECM) N acquisition processes can induce different soil C accumulation in ECM systems compared to arbuscular (AM) systems.In this data package, we included:-The Python codes of the standard Myco-CORPSE model we developed: "Standard Myco_CORPSE python codes.zip". The main program is the "gradient_sim.py" which calculates the bulk soil microbes and CN content along a user defined gradient of clay, soil temperature, soil moisture and mycorrhizal dominance, and relies on two subprograms "CORPSE_deriv.py" and "CORPSE_integrate.py". "CORPSE_deriv.py" calculated the changes in all simulated soil stock within every time step and "CORPSE_integrate.py" integrate the changes in all simulated soil stock within simulated time period. The program "Plot.py" is used to plot the major outputs produced by the main program "gradient_sim.py".-The modified Python codes of Myco-CORPSE models with site-level environmental inputs (in NH and GA) used to conduct simulations in NH and GA sites: "NH_GA model simulations.zip". -The Python codes used to evaluate the Myco-CORPSE simulation outputs in NH and GA sites against site-level measurements: "Plot NH_GA simulation against measurements.zip". It includes both the evaluation Python code, the model outputs on NH and GA sites, and the measured soil properties in both sites.-The modified Python codes of Myco-CORPSE models "Scenario analysis_model simulations.zip" that is used to conduct scenario analysis of how different litter properties, mycorrhizal fungal traits, climate, and seasonal variation in temperature and vegetation phenology impact the mycorrhizal effects on soil CN properties. The sub file folder "Scenario analysis_litter traits" contains the codes for scenario analysis of different litter properties; The sub file folder "Scenario analysis_ECM types" contains the codes for scenario analysis of different ECM fungal traits; The sub file folder "Scenario analysis_climate&seasonality" contains the codes for scenario analysis of different climate and seasonalities;-"Scenario analysis_model results and plotting codes.zip" contains all the output files from the the scenario analysis of Myco-CORPSE model as described above and the plotting codes used to the generate the heatmaps and scatterplots shown in the manuscript "Ectomycorrhizal effects on decomposition are highly dependent on fungal traits, climate, and litter properties: A model-based assessment"The majority of the model outputs did not have specific geographic information or temporal coverage because the analysis we conducted are mainly hypothetical model simulations. We only provided geographic description, coordinates and temporal coverage for those soil measurements which we used for model evaluations (included in the "Plot NH_GA simulation against measurements.zip").

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Improvements to CTF Closure Models for Modelling of Two-Phase Flow

This report documents the efforts to improve the CTF prediction of void fraction and two-phase pressure drop by improving two-phase closure models. Previous validation activities have revealed that CTF tends to overpredict void fraction and two-phase pressure drop. In response to this, two approaches were taken to improve CTF’s predictive capabilities. Based on findings that the interfacial drag and subcooled boiling models significantly impact void prediction, alternative closure models for these physical effects were found and implemented into the code. An extensive assessment was performed by using the existing and newly added two-phase experimental data, which showed that void prediction, wall temperature, and two-phase pressure-drop results are improved by using the newly implemented models. A second approach for improving CTF modeling accuracy involved the use of a Bayesian calibration process that uses CTF validation data to optimize selected modeling coefficients and closure model multipliers to achieve a more accurate prediction of experimental results. A similar assessment was performed with calibrated models that improved the void and pressure-drop prediction for test cases that were both included and not included in the calibration dataset. Results of this study will be used to change the models used in CTF to achieve more accurate BWR analyses moving forward. This study also identified opportunities for improving additional closure models and other opportunities to use calibration techniques to improve CTF.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Identification of Faults Susceptible to Induced Seismicity (Final Report)

Central to the work documented in this report is the capability of geocellular models to represent the geologic conceptual model updated with fault identification from machine learning and joint inversion modeling of microseismic data measured and recorded as a consequence of CO 2 injection at a field demonstration site: the Illinois Basin - Decatur Project (IBDP). This work required seven unique geocellular models with 100s of simulated variations to gain a very high degree of confidence in the identification of geologic features present that contributed to induced microseismicity at IBDP. All forward modeling: pressure modeling, stress modeling, and seismic modeling used the same geologic conceptual model and representations of that model at different scales. The pressure modeling and poroelastic modeling created “snapshots” of pore pressure and stress field changes at different times during CO 2 injection, in which microseismic events were clustered (in time). These pressure and stress snapshots, within the framework and architecture of the geologic conceptual model via the geocellular model, informed the single fault and fault network models to ascertain the likelihood of fault movement (seismic or aseismic). The outcomes of the pressure, stress, and fault/fault network (seismic) modeling confirmed that the faults in the geologic conceptual model in Task 2 were likely the source of microseismic events measured at IBDP and acted as conduits for pressure to be transmitted from the injection interval into the Precambrian crystalline basement rock. This closely coordinated and integrated unique modeling approach was conducted to prove the viability of our proposed workflow 1) to better resolve crystalline basement faults, 2) detect subseismic faults that could be activated by injection, 3) increase the certainty in fault detection and their susceptibility to release seismic energy, and 4) understand transmission of pressure vertically from the well to the underlying fractured crystalline basement. The proposed methodology was effective in guiding an iterative process of calibrating forward modeling results based on similar geocellular models while honoring the geologic conceptual model (i.e., characterization data and knowledge of regional geology); this led to higher level of certainty in the identification of fault/faults zones to control seismicity and transmission of pressure to the regions of recorded and located injection induced seismicity.

58 GEOSCIENCES↗
DOE OSTI2024

FORGE STRESS annual report

The project's goal is to combine high-fidelity numerical models and true-triaxial block fracturing tests at high temperatures to understand the relationship between in situ stress, thermal effects, wellbore orientations and hydraulic fracture patterns. The numerical models are calibrated against field data, such as well pressures and microseismic data, and employed to estimate the in situ stress at the FORGE site. Laboratory experiments investigate the complex physics driving hydraulic fracture nucleation in EGS, enhancing understanding of the role of parameters like temperature, well orientation, and stress. Additionally, they are employed to validate some of the numerical tools used in the project. The project will have a significant impact by: (1) improving the characterization of the in-situ stress field at FORGE; (2) demonstrating the use of high-fidelity modeling tools for EGS; (3) providing a unique set of high-temperature hydraulic fracturing results to identify key components for in-situ stress estimation, validate current theories, and propose new ones; (4) offering a validated set of numerical tools within an open-source simulation framework, GEOS, that will be available to any future user.

15 GEOTHERMAL ENERGY↗
DOE OSTI2026

Towards modeling phage therapy

Patients infected with life-threatening multi-drug resistant (MDR) bacteria have been treated with cocktails of bacteriophages. This is a complicated form of personalized medicine as the phages given to a patient have to be selected beforehand on the basis of their lytic capacity of the infecting bacteria. Because bacteria rapidly become resistant, the evolution of resistance to a diverse cocktail of phages is a complicated dynamical process, during which competing bacterial strains replace one another by accumulating several resistance mechanisms, each of which may involve a fitness cost. As a consequence, it is typically not known why a particular phage therapy succeeded or failed, and how one can optimize the composition of the cocktails to maximize the rate of success. To improve upon this, we extend an existing in vivo -calibrated mouse model into a novel mathematical model for the human situation, and include multiple phages infecting multiple bacterial strains, differing in their resistance to each of the phages. We adjust several parameter estimates of the bacterial model to the human situation, and use the model to describe a successful case of phage therapy involving several cocktails, each containing several phages. In the model, treatment success crucially depended on pretreatment resistance levels, and on the diversity and the timing of the cocktails. Once an appropriate cocktail is found, it is less important to further optimize the infection rates of the phages. Resistant bacterial strains expand rapidly when sensitive strains decline, and the higher the infectivity of the phages, the faster resistant strains expand. Because resistance evolves rapidly, it is best to provide a diverse set of phages right from the start of therapy, i.e., to hit hard and early, and create a high genetic barrier to bacterial resistance.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Training calibration-based counterfactual explainers for deep learning models in medical image analysis

The rapid adoption of artificial intelligence methods in healthcare is coupled with the critical need for techniques to rigorously introspect models and thereby ensure that they behave reliably. This has led to the design of explainable AI techniques that uncover the relationships between discernible data signatures and model predictions. In this context, counterfactual explanations that synthesize small, interpretable changes to a given query while producing desired changes in model predictions have become popular. This under-constrained, inverse problem is vulnerable to introducing irrelevant feature manipulations, particularly when the model’s predictions are not well-calibrated. Hence, in this paper, we propose the TraCE (training calibration-based explainers) technique, which utilizes a novel uncertainty-based interval calibration strategy for reliably synthesizing counterfactuals. Given the wide-spread adoption of machine-learned solutions in radiology, our study focuses on deep models used for identifying anomalies in chest X-ray images. Using rigorous empirical studies, we demonstrate the superiority of TraCE explanations over several state-of-the-art baseline approaches, in terms of several widely adopted evaluation metrics. Our findings show that TraCE can be used to obtain a holistic understanding of deep models by enabling progressive exploration of decision boundaries, to detect shortcuts, and to infer relationships between patient attributes and disease severity.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
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

Coupled Thermo–Hydro–Mechanical Modeling of Hydro-Shearing Stimulation in an Enhanced Geothermal System in the Raft River Geothermal Field, USA

The Raft River Enhanced Geothermal System (EGS) demonstration project aims to improve the geothermal production by enhancing the reservoir permeability via shear stimulation in well RRG-9. In this paper, we performed a series of 3D thermo– hydro–mechanical (THM) simulations to investigate the influence of hydraulic and thermal effects on the development of this EGS. The model includes synthetic fracture populations based on borehole televiewer images and in-situ stress measurements from well RRG-9. Fracture permeability evolution is determined using an empirical permeability law developed from laboratory experiments. The model was calibrated by comparing the hydraulic response to field observations, including wellhead pressure, injection rate, and well injectivity. Particularly, we analyzed the enhancement of reservoir permeability and the spatial extent of the stimulation zone for the given injection schedule. Our results indicate that the permeability enhancement of fractured geothermal reservoir is caused by the combined effects of injection-induced cooling and fluid pressure increase. The decrease of temperature plays a dominant role in reactivation of natural fractures under the hydro-shearing mechanism, while the higher injection pressure promotes shear failure and enlarges the stimulation zone. Finally, for the specified extensional stress state, the model favors greatest permeability enhancement along the maximum horizontal principal stress, a moderate vertical enhancement, and a smaller gain along the minimum horizontal principal stress.

58 GEOSCIENCES↗