Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reactive flow”

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

Unraveling membrane electrode assembly design for electrochemical conversion of carbon dioxide to formate/formic acid

This work presents a one-dimensional continuum modeling approach to investigate various cell architectures used for electrochemical conversion of CO 2 to formate/formic acid. Ion transport is simulated by a system of generalized modified Poisson–Nernst–Planck (GMPNP) equations that reflect the reactive transport phenomena including steric effects as the electrolyte solutions become concentrated. In the cathode catalyst layer, ionic current contributions from both the supporting electrolyte and solid-state ionomer are considered. Voltage and CO 2 utilization breakdowns are utilized to deconvolute the impacts of the cell architecture. The origins of (bi)carbonate formation in the cathode are explored, as the subsequent decrease in CO 2 availability is a key reason for low faradaic efficiencies to formate/formic acid. In addition, the role of a supporting electrolyte (KOH) is investigated to understand its tradeoffs: while the K + ions can improve both conductivity and electrochemically active surface area in the cathode, the presence of OH − ions raises the pH and leads to deleterious formation of (bi)carbonates. To this end, we also present parametric studies on the concentration and flow rate of supplied KOH to the cell, to establish a path towards eliminating the need for a supporting electrolyte.

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

Non-Boltzmann Effects in Chain Branching and Pathway Branching for Diethyl Ether Oxidation

Low-temperature (LT) engine applications have several potential benefits, including reduced emissions and increased efficiency. Attaining these benefits requires accurate kinetic modeling of LT chain branching, which depends heavily on ketohydroperoxide (KHP) decomposition. For diethyl ether (DEE), a promising biofuel, current estimates of the KHP decomposition rate constant are largely based on empirical fits to data. In this study, we investigate the most important KHP isomer in DEE LT oxidation by applying variable reaction coordinate transition state theory to the main pathway for KHP decomposition: OO bond fission to produce •OH and a keto-alkoxy radical, •OQ'O. We also use ab initio kinetics methods to investigate the decomposition of •OQ'O, where we find dominant branching to acetic acid, with the remaining flux going to CH 3 C(O)OCHO. Additionally, new time-resolved measurements of DEE and acetic acid concentrations during LT (450–600 K) DEE oxidation are obtained in a laser photolysis flow reactor coupled with multiplexed photoionization mass spectrometry. These new experimental data, along with jet-stirred reactor data in the literature, are compared with the predictions of a recent DEE mechanism (Tran et al. Proc. Comb. Inst. 2019, 37, 511-519) that was modified with the newly calculated ab initio rate constants for KHP and •OQ'O decomposition. The predictions of the modified mechanism are quite poor when compared to the experimental data; this is primarily due to the new KHP ⇌ •OQ'O + •OH rate constant, which is 1–2 orders of magnitude slower than empirical values employed in recent mechanisms. To reconcile the new KHP rate constant and the experimental data, we explore and quantify the possible role of non-Boltzmann (nB) reaction sequences. The nB reactions have a substantial effect on both the overall mechanism reactivity and the •OQ'O branching to acetic acid. We also provide guidance on the proper implementation of nB reactions in kinetic mechanisms.

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

Analysis and design of DC-link capacitor in cascaded H-bridge multilevel active front-end converter

Medium-voltage grid-tied systems often use a cascaded H-bridge multi-level active front-end. In this converter, dc link bus capacitors play an important role in stabilizing the converter and enabling both active and reactive power injections. The present work provides analytical expressions for the capacitor current, which are essential for optimizing system design (especially capacitor size vs. lifetime). Then, the expression is incorporated into the grid connected bidirectional power system model. Consequently, this work contributes to the guiding principles to choose accurate dc link capacitor ratings against grid-side power delivery requirements. DC-link capacitors come with an expiration date which is dependent on the operating power and voltage values. Lack of information about the lifetime of the dc-link capacitor creates uncertainty in the duration of online operation of the converter, thus increasing the probability of a contingency outage. The lifetime of the capacitor can range from only a few seconds of online operation of the station due to under-sizing, to a compromised quality of power flowing through other solid state devices which may not have been designed for the transients that are caused by over-sizing of the dc-link capacitor. A better estimate of the lifetime of these capacitors can help maximize usage and maintain scheduled outages more efficiently, without disruption in service.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Factors controlling injection-induced rupture of intersecting faults during geological sequestration of CO 2

This study addresses coupled multiphase fluid flow and geomechanics effects on potential fault activation associated with subsurface CO 2 injection around intersecting faults. An enhanced fault-representation model is used to capture geomechanical responses of two intersecting faults with finite length during CO 2 injection. The faults are embedded in a strike-slip stress regime of a caprock-reservoir-basement system with the faults represented by zero-thickness interfaces with adjacent finite-thickness damage zones. A sensitivity analysis is conducted to study the effect of fault permeability, slip-weakening behavior, well location relative to the orientation of faults, and well placement (the number and location of injection wells). Five metrics (pressure, CO 2 plume, shear state on the fault, as well as shear displacement and stress path at selected fault monitoring points) are selected to assess CO 2 migration and reactivation of intersecting faults. The results show that induced ruptures are favored by low permeability faults due to high pressure buildup and by slip-weakening behavior resulting from fault strength reduction. The location of one injection well relative to fault orientation determines the magnitude of changes in effective normal stress and shear stress, affecting the location of induced ruptures. Well placement (two injection wells used in the paper) dominates pressure diffusion around the intersection and tips of faults. This redistributes changes in effective normal stress caused by each injection well, influencing the spatial distribution of ruptures along faults. A larger injection volume induces far-field ruptures that are controlled by stress transfer within the injection layer. The findings presented here can provide valuable insights into engineering operations for a long-term, safe, and reliable geologic CO 2 storage.

Fault permeability↗
DOE OSTI2023

Model simulations of Plum Island Ecosystems LTER low marsh site using ELM-PFLOTRAN

Model simulations using the E3SM Land Model (ELM) coupled to the PFLOTRAN reactive transport model via the Alquimia interface. The simulations were conducted for a tidal salt marsh at the Plum Island Ecosystems LTER near Rowley, Massachusetts, USA. Model simulations were forced using site-specific tidal cycles and salinity, and the simulations used a biogeochemical reaction network including aerobic decomposition, sulfate reduction, iron reduction, and methanogenesis. Model outputs include simulated carbon stocks, carbon dioxide and methane fluxes, and porewater concentrations of key solutes related to sulfur, iron, and carbon cycling. The model simulations included a saline simulation (with tidal sulfate inputs), a fresh simulation (with low salinity and low sulfate inputs), and a saline simulation with lower vegetation productivity to represent the effect of salinity on vegetation. These simulations were conducted to demonstrate that a new model framework incorporating subsurface redox and biogeochemical interactions into a land surface model could reproduce measured surface greenhouse gas fluxes and biogeochemical dynamics in tidal marsh ecosystems, and to test whether including redox interactions in a land surface model would allow the model to resolve contrasts in biogeochemical cycling and greenhouse gas production between saline and freshwater wetlands.The data package includes gzipped tar archives (which can be expanded using standard tar and gzip utilities) of model outputs from three model configurations: saline subsurface and reduced vegetation productivity related to salinity; saline subsurface with vegetation productivity not reduced; and freshwater. Also included are code for the modified E3SM model, Alquimia interface, and PFLOTRAN reactive transport simulator in gzipped tar format; plain text parameter and configuration files; python code files for visualizing model output and defining model configurations; and model output, tide and salinity forcing, and configuration files in netCDF format. See the README.md file in the data package for a detailed description of all files contained in the package. All files are in netCDF (.nc), gzipped tar archive (.tar.gz or .tgz), or text (all other files).Updated: May 13, 2024. Model output, E3SM code, PFLOTRAN input files, and python codes for visualizing results were updated to reflect changes made for the manuscript revision. The updated archive reflects the code and model output from the final accepted manuscript. Changes included updated reaction parameters reflecting improved parameterization and additional comparisons with field measurements. E3SM code changes included better support for multiple grid cells and improved flow and transport parameterization.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Nafion Inhibits Polysulfide Crossover in Hybrid Nonaqueous Redox Flow Batteries

Here, we report on the effectiveness of sodium-exchanged Nafion membranes for inhibiting polysulfide crossover during redox flow battery operation. The solubility of polysulfides allows them to be used as a high-capacity catholyte for nonaqueous redox flow batteries (NARFB). The NARFB cathode capacity is controlled by the total tank catholyte volume and polysulfide concentration and, as such, is independent of polysulfide adsorption sites on carbon fibers or nanotubes as found within traditional Na or Li–S batteries. However, one of the major barriers to the realization of polysulfide NARFBs is associated with developing a robust membrane that has minimal polysulfide crossover under the operating conditions. We found that sodium-exchanged Nafion membrane separators effectively inhibit polysulfide crossover within a Na–sulfur NARFB. These membranes significantly improve both capacity retention (70%) and Coulombic efficiencies (99%) after 50 cycles. Commercial porous membranes showed a large crossover as detected by UV–vis spectroscopy and resulted in low capacity retention (20%) and Coulombic efficiency (74%) after 45 cycles. Electrochemical impedance spectroscopy (EIS) measurements highlighted the trade-off between innate reactivity and ionic conduction of the membranes. The results show that dense, single ion conducting Nafion enables a long cycle life; however, it reacts with Na metal to form a resistive passivation layer and increases the cell resistance. On the contrary, the open pore structure of Celgard allows for higher current charge/discharge, and its chemical nature is compatible with Na; however, it has a high degree of polysulfide crossover.

25 ENERGY STORAGE↗
DOE OSTI2025

Resonance-stabilized radical clustering bridges the gap between gaseous precursors and soot in the inception stage

Carbonaceous particles are widespread in combustion, atmospheric, extraterrestrial, and nanomaterials environments. Resonance-stabilized radicals (RSRs) are commonly identified in fuel combustion and pyrolysis processes and play an essential role in carbonaceous particle formation. Despite their importance, comprehensive experimental and mechanistic understanding of particle inception through RSR reactions is lacking. This work investigated particle size distribution, chemical composition, and thermal behavior of soot particles generated by the flow reactor pyrolysis reactions of typical RSRs, in particular, 1-indenyl, 1-methylnaphthyl, and 2-methylnaphthyl radicals, and by the pyrolysis of hydrocarbons with a variety of structures. Particle size distributions show soot particles with mobility diameters in an incipient-particle range of 1.3 to 1.6 nm. Laser desorption/ionization mass spectrometry results suggest that soot products consist of much larger covalently bound clusters (CBCs) than those observed in the gas phase. Under our experimental conditions, the CBCs exhibit a phase transition for particles with calculated molecular diameters of around 1.5 nm. Evaporation experiments and thermogravimetric analysis of the soot products reveal distinct thermal characteristics for small and large CBCs. These results implicate CBCs as bridges between gas-phase species and soot particles. The present work provides a soot-inception mechanism called RSR clustering (RSRC) that is characterized by the reactive clustering of RSRs. The RSRC mechanism contrasts with conventional soot formation models that attribute soot inception primarily to the aggregation of large-size polycyclic aromatic hydrocarbons.

carbonaceous particle↗
DOE OSTI2023

Fracture Intensity Impacts on Reaction Front Propagation and Mineral Weathering in Three‐Dimensional Fractured Media

Abstract Studying reaction front propagation in heterogeneous natural settings is challenging, but numerical simulations can provide insight into the varying spatial and temporal scales of reaction front propagation. Here, the impact of increasing fracture intensity on mineral dissolution rates, and the extent of reaction front propagation is investigated using reactive transport simulations in upscaled discrete fracture network domains with varied fracture intensity. Domain‐averaged dissolution rates vary less than 0.5 log units regardless of the fracture intensity, but the spatial distribution of reactions is controlled by the location and number of dead‐end fractures and the number of connected flowpaths through the domain. Higher fracture intensities lead to more weathering in the domain because of more available mineral for water‐rock interactions. We find that reaction fronts propagate through the primary flowpaths in the first 10,000 years of the simulation for a 10‐m length domain, then propagate into secondary flowpaths and dead‐end fractures between 10,000 and 100,000 years, and finally into the matrix over timescales of hundreds of thousands of years. The domain‐averaged reaction rates decrease through time corresponding to a transition from dissolution in advection‐dominated, fast‐flowing pathways, to dissolution in transport‐limited zones of disconnected fractures and matrix. Matrix dissolution, or dissolution under transport‐limited conditions, is the dominant process at late times in these simulations. The results of these simulations recreate the observed paradox found in nature where highly fractured hillslopes tend to be more weathered but have slower weathering rates, while hillslopes with fewer fractures, are less weathered but have higher dissolution rates.

58 GEOSCIENCES↗
DOE OSTI2023

Roles of Mn and Co in the Air Synthesizability of Layered Oxide Cathodes for Lithium-Based Batteries

High-nickel layered oxides (LiNi 1-x-y Mn x Co y O 2 ) are the prevailing cathode materials for high-energy-density lithium-based batteries, but they are plagued with deleterious surface air instabilities stemming from residual lithium formation. These issues severely hinder mass production as cathode calcination is limited to a flowing oxygen atmosphere, which entails high manufacturing costs as opposed to simpler and more economical air calcination. Here, while higher Ni contents are known to worsen air instabilities, the influence of Mn and Co contents on impacting these phenomena are less elucidated. We herein present the synthesis in ambient air and flowing oxygen atmospheres of three cathode variants with the same Ni contents, but varying Mn and Co contents: LiNi 0.7 Mn 0.3 O 2 , LiNi 0.7 Mn 0.15 Co 0.15 O 2 , and LiNi 0.7 Co 0.3 O 2 . It is found that the critical parameter influencing the air stability of the cathodes is the average Ni oxidation state, which is greatly dependent on the Mn and Co contents. Substitution of Mn for Ni drives down the Ni oxidation state as Mn exists as Mn4 + and reduces surface residual lithium formation, which vastly improves the overall air stability and, therefore, the synthesizability in air, but with a penalty of lowered capacity. In contrast, substitution of Co for Ni maintains Ni 3+ as Co exists as Co 3+ , offering increased initial capacity, but worsens the air stability and cyclability as the driving force for residual lithium formation and surface reactivity is increased.

25 ENERGY STORAGE↗
DOE OSTI2024

Simulation of Particulate Transport for Delivery of Solid Amendments into the Subsurface: FY24 Status Report

For particulate-based amendments to be viable for field-scale remediation at the Hanford Site (e.g., 200 DV-1 Operable Unit), particles need to be delivered a sufficient radial distance from an injection well and retained at concentrations high enough for effective treatment. An accurate description of the particle radius of influence (ROI) is critical for developing an overall remediation strategy. However, field-scale particle simulations are currently limited due to insufficient simulation capabilities and a lack of experimental data to validate and parameterize particle transport models. To help build toward field-scale deployment, this fiscal year (FY) we have (1) developed a pre screening tool to estimate particle transport, (2) implemented particle transport models within PFLOTRAN, and (3) conducted preliminary estimations of particle ROI. While field-scale numerical simulations will ultimately be necessary before remedy design and field implementation, we have developed a pre-screening tool that offers valuable estimations of expected particle injectability and ROI in a 1-D system. The advantage of the tool is that it does not require extensive laboratory experiments and instead makes predictions based solely on routine laboratory measurements. This tool can assist in down-selection and decision-making by identifying which particle amendment systems are worth pursuing in future laboratory experiments, such as 1-D column tests and beyond. With any system, scaling up from the lab to the field presents challenges. Currently, there is no field data available for model calibration or validation. However, the theoretical particle models being developed herein are the best tools available to guide progress toward field deployment. To help bridge this gap and verify model predictions, larger-scale lab experiments are being proposed. To advance simulation capabilities, six particle transport models are being integrated into the reactive transport simulator PFLOTRAN. These include colloid filtration theory (CFT) and five additional particle transport models (M1-M5). Each model, from M1 to M5, progressively incorporates additional particle transport and retention processes. Ultimately, the simplest model capable of accurately describing 1-D column data will be selected and parameterized. During FY24, the CFT and M1 model have been fully implemented within PFLOTRAN. Using an existing 1 D column experiment, the two currently implemented particle transport models (CFT and M1), and associated parameters, were fit to this experiment. While simpler model formulations are helpful for estimations, these formulations could not fully describe particle transport and retention behavior in the previous 1-D column experiment. Thus, additional complexities will need to be considered, which will be accounted for in the M2-M5 model formulations. Additionally, because a viscous, shear thinning fluid was required to keep particles in suspension, considerations for flow will also need to also be accounted for. Therefore, a new immiscible two-phase flow mode is currently being implemented in PFLOTRAN. With some modifications, this new flow module could also support simulation of non-Newtonian liquid amendments, foams, and emulsions. We also estimated the expected ROI of solid amendments using 1-D simulations. The average predicted ROI was approximately 15 ft for micron-sized zero valent iron (mZVI) suspended in xanthan gum (XG). Using the pre screening tool and ROI estimates, additional amendment-delivery laboratory characterization and experiments are proposed. The results from additional experiments can be used to validate and parametrize particulate transport model formulations, which will ultimately provide predictive capabilities for field amendment-delivery systems.

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

Predicting Fracture Porosity Evolution in Sandstone (Final Report)

To better understand porosity, strength, chemical reactivity, and patterns of fractures, we developed methods using mineral deposits to unravel how fracture growth and diagenesis interact to create and destroy fracture porosity. Quartz and dolomite cement textures and associated fluid inclusion assemblage sequences and thermal histories provided data to infer the duration and rates at which fractures open. We developed and tested a model that accounts for how fracture porosity and rock properties change as a function of thermal exposure and how fracture size and arrangement evolve. Effects of cement deposits on growth were investigated with new geomechanical models and new methods to quantify fracture spatial arrangement. We tested these concepts in a wide range of rock types using outcrop and core-based datasets. We documented opening histories of fractures and fracture arrays in contrasting tectonic and thermal regimes. Our methods set a new standard for fracture characterization and our diagenetic-geomechanical model accurately predicts fluid flow characteristics of fractured rocks based on coupled effects of diagenesis and deformation. Results represent a significant advance in fundamental understanding of fractures in the Earth and have application to improving subsurface engineering operations in a wide range of contexts including geothermal and sequestration applications. The breakthrough of successfully reconstructing the timing and opening rates of parts of single fractures and portions of fracture arrays to document the growth and linkage of fracture arrays was an accomplishment of this grant that is being followed up in our current BES grant DE-SC0022968, ‘Reconstructing and Predicting Fracture Pattern Evolution‘

58 GEOSCIENCES↗
DOE OSTI2020

Colloid-facilitated transport of 238 Pu, 233 U and 137 Cs through fractured chalk: Laboratory experiments, modelling, and implications for nuclear waste disposal

The influence of montmorillonite colloids on the mobility of 238 Pu, 233 U and 137 Cs through a chalk fracture was investigated to assess the transport potential for radioactive waste. Radioisotopes of each element, along with the conservative tracer tritium, were injected in the presence and absence of montmorillonite colloids into a naturally fractured chalk core. In parallel, batch experiments were conducted to obtain experimental sorption coefficients (K d , mL/g) for both montmorillonite colloids and the chalk fracture material. Breakthrough curves were modelled to determine diffusivity and sorption of each radionuclide to the chalk and the colloids under advective conditions. Uranium sorbed sparingly to chalk (log K d = 0.7 ± 0.2) in batch sorption experiments. 233 U(VI) breakthrough was controlled primarily by the matrix diffusion and sorption to chalk (15 and 25% recovery with and without colloids, respectively). Cesium, in contrast, sorbed strongly to both the montmorillonite colloids and chalk (batch log K d = 3.2 ± 0.01 and 3.9 ± 0.01, respectively). The high affinity to chalk and low colloid concentrations overwhelmed any colloidal Cs transport, resulting in very low 137 Cs breakthrough (1.1–5.5% mass recovery). Batch and fracture transport results, and the associated modelling revealed that Pu migrates both as Pu (IV) sorbed to montmorillonite colloids and as dissolved Pu(V) (7% recovery). Transport experiments revealed differences in Pu(IV) and Pu(V) transport behavior that could not be quantified in simple batch experiments but are critical to effectively predict transport behavior of redox-sensitive radionuclides. Finally, a brackish groundwater solution was injected after completion of the fracture flow experiments and resulted in remobilization and recovery of 2.2% of the total sorbed radionuclides which remained in the core from previous experiments. In general, our study demonstrates consistency in sorption behavior between batch and advective fracture transport. Finally, the results suggest that colloid-facilitated radionuclide transport will enhance radionuclide migration in fractured chalk for those radionuclides with exceedingly high affinity for colloids.

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

Coupled Decay Heat and Thermal Hydraulic Capability for Loss-of-Coolant Accident Simulations

As the nuclear energy industry considers ways to achieve improved economics in the current fleet of light-water reactors (LWRs), one possible approach is to operate each cycle for longer durations. This causes a greater portion of the fuel to be burned and reduces the frequency of outages, which ultimately reduces the cost to operate the reactor. However, this also leads to higher burnup fuels than has traditionally been allowed in these reactors. Thus, there are concerns about integrity of high-burnup (HBu) fuel, especially during accident conditions such as loss-of-coolant accidents (LOCAs), as shown by Capps et al.. To investigate these concerns, advanced modeling and simulation capabilities are being leveraged to determine the susceptibility of HBu fuel to fuel fragmentation, relocation, and dispersion (FFRD). Improvements have previously been made to fuel performance capabilities to more accurately model these phenomena; multiphysics simulations have also been conducted to determine the power and burnup histories of the HBu fuel, which are needed as inputs for the fuel performance calculations. Most recently, new statistical approaches have been developed to identify a subset of fuel rods that have greater FFRD susceptibility, reducing the total number of fuel performance simulations required. Prior LOCA simulations have relied on the TRACE systems code, which can model the core and primary loop during accident conditions. TRACE includes many models for various aspects of the primary loop, but two sets of models are important for this report. First, TRACE uses a lumped-fuel approach for modeling the core. This approximates the ~50,000 fuel rods in the core with a much smaller number of rods. The rods can be lumped in various ways as determined by the user. For example, one lumped rod may be used to represent all rods in an assembly, sometimes with an additional rod representing the hottest fuel rod. However, due to runtime constraints and complexity of modeling, a more common approach is to group several assemblies or larger regions of the core into single lumped rods. These lumping schemes apply not only to fuel rods but to flow channels as well. Second, TRACE has several different models for treating decay heat, ranging from pregenerated decay heat curves based on an ANSI/ANS-5.1 standard (hereinafter abbreviated simply as ANSI) to explicit time-dependent heat inputs from the user. None of these models account for differences in isotopics between different rods, which is an approximation the work in this report seeks to eliminate. This report focuses on the implementation of coupled decay heat capabilities in the Virtual Environment for Reactor Applications (VERA) code suite to address a gap identified in previous LOCA simulations. This constitutes an improvement for both the lumped-fuel and decay heat models in TRACE. VERA has been developed to perform high-fidelity, whole-core multiphysics simulations for LWRs. Previously, during the Consortium for Advanced Simulation of LWRs (CASL) program, the emphasis was on providing accurate steady-state analysis—with a secondary focus on reactivity insertion accident (RIA) analysis—to address operational challenges in the nuclear energy industry. Under the Department of Energy (DOE) Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, these capabilities are being extended to a broader range of transient analyses with the goal of quantifying the risk of fuel failures such as FFRD. To properly model such conditions with VERA, decay heat calculations have been integrated with the multiphysics to enable rod-by-rod thermal hydraulic (TH) conditions to be driven by the decay heat in long-running accidents such as LOCAs.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

2D-to-3D image translation of complex nanoporous volumes using generative networks

Image-based characterization offers a powerful approach to studying geological porous media at the nanoscale and images are critical to understanding reactive transport mechanisms in reservoirs relevant to energy and sustainability technologies such as carbon sequestration, subsurface hydrogen storage, and natural gas recovery. Nanoimaging presents a trade off, however, between higher-contrast sample-destructive and lower-contrast sample-preserving imaging modalities. Furthermore, high-contrast imaging modalities often acquire only 2D images, while 3D volumes are needed to characterize fully a source rock sample. In this work, we present deep learning image translation models to predict high-contrast focused ion beam-scanning electron microscopy (FIB-SEM) image volumes from transmission X-ray microscopy (TXM) images when only 2D paired training data is available. We introduce a regularization method for improving 3D volume generation from 2D-to-2D deep learning image models and apply this approach to translate 3D TXM volumes to FIB-SEM fidelity. We then segment a predicted FIB-SEM volume into a flow simulation domain and calculate the sample apparent permeability using a lattice Boltzmann method (LBM) technique. Results show that our image translation approach produces simulation domains suitable for flow visualization and allows for accurate characterization of petrophysical properties from non-destructive imaging data.

58 GEOSCIENCES↗
DOE OSTI2022

Fast modular reactor nuclear design parameters of fuel cycle and power distributions

This paper describes conceptual design of a 112 MW thermal (50 MW electric) Fast Modular Reactor (FMR) system operating at 7 MPa with inlet/outlet temperatures of 509/800 Celsius degrees. The reactor system includes the reactor core, fuel assemblies, fuel rods, reactor internals, reflector, neutron control system, flow control, and structural components. The nuclear design and analysis were conducted to search for a baseline core with a cycle length greater than 8-year and the power peaking factor less than 1.5 by adjusting the fuel assembly, reactor core, and reflector configurations. The neutron transport calculations of the baseline core showed that a refueling interval of 9-year is achieved with a total peaking factor of 1.47 and a fuel rod average linear power of 3.6 kW/m. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Experimental study of fracture slip due to stress perturbation in fractured geo-resources

Anthropogenic activities such as underground fluid injection/extraction, well stimulation, and dam impoundment will perturb the state-of-stress in the subsurface, which could activate critically stressed pre-existing fractures and faults. Fault and fracture activation can induce seismicity, change flow behavior, and damage asperities on the fracture surface to alter reservoir permeability. Here, in this study, triaxial shear tests were conducted on saw cut and rough fracture Barre granite specimens to better characterize how a decrease in normal stress and increase in shear stress will influence crystalline rock fracture hydromechanical properties. In a triaxial set up, confining pressure (CP) and differential stress (DS) were used to induce normal and shear stress on a 45° oriented fracture. Fracture slip and the corresponding stress relaxation was observed when the confining stress decreased. Fracture permeability increased when the fracture slipped. Moreover, new micro-cracks observed by x-ray images indicated the creation of new flow paths for fluid movement. The experimental outcomes of this study indicate that decrease in normal stress due to stress perturbations along the faults/fractures in the extensional stress regime can induce the slip along the fracture. This could have a significant influence on slope stability in quarries, stimulation and injection strategy in deep geo-resources, and potentially minimize the risk of induced seismicity. The observed stress dependency of hydraulic aperture in rough fracture during pre- and post-slip stages could provide key inputs for field-scale simulation of geo-resources.

58 GEOSCIENCES↗
DOE OSTI2023

Control of Cement Timing, Mineralogy, and Texture on Hydro-chemo-mechanical Coupling from CO2 Injection into Sandstone: A Synthesis

Carbon capture, utilization, and storage (CCUS) has been widely applied to enhance oil recovery (CO2-EOR). A thorough investigation of the impact of injecting CO2 into a heterogeneous reservoir is critical to understanding the overall reservoir robustness and storage performance. We conducted fifteen flow-through tests on Morrow B sandstone that allowed for chemical reactions between a CO2-rich brackish solution and the sandstones, and four creep/flow-through tests that simultaneously allowed for chemical reactions and stress monitoring. From fluid chemistry and X-ray computed tomography, we found that the dissolution of disseminated cements and the precipitation of iron-rich clays did not significantly affect the permeability and geomechanical properties. Minor changes in mechanical properties from Brazilian and creep tests indicated that the matrix structure was well-supported by early diagenetic quartz overgrowth cement and the reservoir’s compaction history at deep burial depths. However, one sample experienced a dissolution of poikilotopic calcite, leading to a permeability increase and significant tensile strength degradation due to pore opening, which overcame the effect of the early diagenetic cements. We concluded that the Morrow B sandstone reservoir is robust for CO2 injection. Most importantly, cement timing, the abundance and texture of reactive minerals, and the reservoir’s burial history are critical in predicting reservoir robustness and storage capacity for CO2 injection.

Wu, Zhidi (ORCID:0000000338869492)↗
DOE OSTI2021

SPRUCE Outflow Chemistry Data for Experimental Plots Beginning in 2016

This data set includes the chemistry of outflow waters in the SPRUCE experiment plots located in the S1 bog of the USDA Forest Service Marcell Experimental Forest (MEF) in northern Minnesota, 40 km north of Grand Rapids. These data are post-treatment data from the warming and elevated carbon dioxide (CO2) treatments associated with the SPRUCE experiment. There are ten experimental plots in SPRUCE: five temperature treatments (+0, +2.25, +4.5, +6.75, +9°C) at ambient CO2, and the same five temperature treatments at elevated CO2 (+500 ppm). Sample collection and analyses started in March of 2016 and will continue for the duration of the experiment. Samples were analyzed for pH, specific conductivity, anion concentrations (chloride, sulfate), cation concentrations (calcium, potassium, magnesium, sodium, aluminum, iron, manganese, silicon, strontium), nutrient concentrations (ammonium, nitrate, soluble reactive phosphorus, total nitrogen, total phosphorus), and total organic carbon (TOC) concentrations. Subsets of the samples have been analyzed for natural-abundance stable isotopes of water (δ18O-H2O and δD-H2O), dissolved carbon dioxide (δ13C-CO2), and dissolved methane (δ13C-CH4), and concentrations of total mercury (Thg), methylmercury (MeHg), ferrous iron (Fe2+), ferric iron (Fe3+), dissolved inorganic carbon (DIC), and dissolved methane (CH4).

54 ENVIRONMENTAL SCIENCES↗