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 1,603 records · Page 89
All sourcesNASA NTRSDOE OSTI
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↗
NASA NTRS2010

Compact Rare Earth Emitter Hollow Cathode

A compact, high-current, hollow cathode utilizing a lanthanum hexaboride (LaB6) thermionic electron emitter has been developed for use with high-power Hall thrusters and ion thrusters. LaB6 cathodes are being investigated due to their long life, high current capabilities, and less stringent xenon purity and handling requirements compared to conventional barium oxide (BaO) dispenser cathodes. The new cathode features a much smaller diameter than previously developed versions that permit it to be mounted on axis of a Hall thruster ( internally mounted ), as opposed to the conventional side-mount position external to the outer magnetic circuit ("externally mounted"). The cathode has also been reconfigured to be capable of surviving vibrational loads during launch and is designed to solve the significant heater and materials compatibility problems associated with the use of this emitter material. This has been accomplished in a compact design with the capability of high-emission current (10 to 60 A). The compact, high-current design has a keeper diameter that allows the cathode to be mounted on the centerline of a 6- kW Hall thruster, inside the iron core of the inner electromagnetic coil. Although designed for electric propulsion thrusters in spacecraft station- keeping, orbit transfer, and interplanetary applications, the LaB6 cathodes are applicable to the plasma processing industry in applications such as optical coatings and semiconductor processing where reactive gases are used. Where current electrical propulsion thrusters with BaO emitters have limited life and need extremely clean propellant feed systems at a significant cost, these LaB6 cathodes can run on the crudest-grade xenon propellant available without impact. Moreover, in a laboratory environment, LaB6 cathodes reduce testing costs because they do not require extended conditioning periods under hard vacuum. Alternative rare earth emitters, such as cerium hexaboride (CeB6) can be used in this configuration with possibly an even longer emitter life. This cathode is specifically designed to integrate on the centerline of a high-power Hall thruster, thus eliminating the asymmetries in the plasma discharge common to cathodes previously mounted externally to the thruster s magnetic circuit. An alternative configuration for the cathode uses an external propellant feed. This diverts a fraction of the total cathode flow to an external feed, which can improve the cathode coupling efficiency at lower total mass flow rates. This can improve the overall thruster efficiency, thereby decreasing the required propellant loads for different missions. Depending on the particular mission, reductions in propellant loads can lead to mission enabling capabilities by allowing launch vehicle step-down, greater payload capability, or by extending the life of a spacecraft.

Watkins, Ronald↗
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 OSTI2019

Coupled Geomechanical Modeling to Assess Cap Rock Integrity and Mechanical Fault Stability: Application to Farnsworth Field Unit Project

Abstract The SWP project is located in a mature waterflood undergoing conversion to CO2-WAG operations in the Farnsworth, Texas, USA. Anthropogenic CO2 is sourced from a fertilizer and an ethanol plant. This work utilizes Farnsworth’s full-field, history-matched, compositional hydrodynamic coupled geomechanical model for assessing the impact of stress changes observed through the history matched field life. Production and injection induced stress changes, fault stability and caprock integrity investigations are performed to project the potential for fault reactivation and the loss of caprock integrity under shear failure. A static mechanical earth model (MEM) was constructed for use in transient coupled geomechanical model based on the existing Southwest Regional Partnership (SWP) geological model. The static MEM inherits the stratigraphic and structural features of the geologic model and incorporates additional overburden, underburden, and sideburden formations required to impose mechanical boundary conditions. Mechanical properties were distributed in the 3D MEM through integration of geophysical logs and 3-dimensional seismic elastic inversion properties using a combination of Bayesian and stochastic interpolation methods. These data are further enhanced by lab derived strength and failure criteria for the caprock interval. Additionally, interpreted faults and other geological features were included as part of the static structural framework to fully represent subsurface hydraulic and mechanical systems and appropriately integrate heterogeneity. Two-way coupling of hydrodynamic flow and geomechanical simulations incorporates the Kozeny-Carman relationship for updating permeability and is history matched through primary, secondary (waterflooding) and tertiary (CO2 WAG) recovery periods before performing a twenty year forecast. Two-way simulations are performed to understand the effective stress perturbations imposed by field operations: water injection and the more recent implementation of WAG. Evaluation of Mohr circles with liquid production and field pressure charts as well as slip tendency and distance to failure metrics indicate that neither our faults or cap rock are critically stressed. This paper presents the results of the Farnsworth initial attempts to integrate the seismic driven 3D MEM into coupled hydrodynamic geomechanical history match simulation workflow.

02 PETROLEUM↗
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 OSTI2019

Analysis of the Detonation Wave Structure in a Linearized Rotating Detonation Engine

The rotating detonation engine is increasingly favored as a viable pressure gain combustion technology for both propulsion and power generation applications. Practical designs involve the discrete injection of fuel and air, which then partially mix to produce the reactive mixture that is processed by a continuously moving detonation wave within the detonation chamber. One of the fundamental challenges in the successful operation of rotating detonation engines is the design of a robust fuel-oxidizer injection system. To minimize pressure losses while ensuring efficient mixing that promotes stable detonation, the injection has to satisfy many different constraints. In this work, the impact of discrete injection on wave propagation is studied in an effort to understand the impact of fuel stratification on detonation stability. To this end, the structure of the detonation wave within a linearized rotating detonation engine with fully premixed hydrogen-air and hydrogen-oxygen injection at various equivalence ratios is analyzed. The direct numerical simulation of the linearized model detonation engine is compared with experimental results, and wave behavior is correlated to the wave structure and flow turbulence. As a result, the flow properties and chemical composition of the gases across the detonation wave are studied to examine the characteristics of the reaction zone.

42 ENGINEERING↗
NASA NTRSDate not supplied

Investigating Material Behavior in Atmospheric Entry Conditions: Arc-Jet Testing Insights from Meteorite Ablation to High-Temperature Coatings.

Arc-jet testing is an indispensable tool in elucidating the complex interactions materials undergo when subjected to the extreme thermal and mechanical stresses of atmospheric entry. By replicating these harsh conditions, arc-jets enable researchers to explore critical phenomena such as ablation, melting, and spallation, generating invaluable data that informs the development of sophisticated, physics-based models. These models are pivotal not only in predicting material performance for spacecraft re-entry systems but also in assessing the threat posed by celestial bodies as they encounter Earth’s atmosphere. One of the leading facilities for such investigations is the Hypersonic Materials Environmental Test System (HyMETS) at NASA’s Langley Research Center. A pioneering test campaign conducted at HyMETS focused on unraveling the ablation mechanisms of an ordinary chondrite meteorite (Tamdakht H5) and a terrestrial analog (basalt). These studies unveiled distinct material behaviors: Tamdakht demonstrated a remarkably stable melt flow, with mass loss primarily driven by the volatilization of elements such as iron and sodium, whereas basalt exhibited more aggressive surface degradation due to the rapid decomposition of hydrated minerals, culminating in significant spallation and surface material ejection. Beyond meteorite analysis, HyMETS has also been integral in the exploration of advanced Thermal Protection Systems (TPS) such as Phenolic Impregnated Carbon Ablator (PICA), a material with extensive flight heritage in missions like Stardust, Mars Science Laboratory, Mars 2020, and Osiris-REx. PICA, when coated with NuSil CV-1144-0—a polysiloxane resin designed to prevent particle shedding—undergoes a remarkable transformation upon heating. The resin pyrolyzes to form a thin, oxidation-resistant silicon oxycarbide layer that profoundly influences the material’s thermal response through a sophisticated four-stage process. Initially, the silicon oxycarbide acts as a formidable barrier, impeding reactive interactions with the PICA char and effectively suppressing surface temperature rise and material recession. However, as the heating persists, the protective layer decomposes via carbothermal reduction, exposing the underlying char, which in turn leads to a dramatic spike in surface temperature and accelerated material erosion. In the final stage, equilibrium is reached, with recession rates aligning with those of the virgin material. This intricate understanding of PICA-NuSil behavior under extreme conditions offers valuable insights, enhancing material response models and advancing the development of next-generation TPS for space exploration.

Chemistry↗
NASA NTRS2002

Sympathetic vascular transduction is augmented in young normotensive blacks

The purpose of the present study was to determine sympathetic vascular transduction in young normotensive black and white adults. We hypothesized that blacks would demonstrate augmented transduction of muscle sympathetic nerve activity (MSNA) into vascular resistance. To test this hypothesis, MSNA, forearm blood flow, heart rate, and arterial blood pressure were measured during lower body negative pressure (LBNP). At rest, no differences existed in arterial blood pressure, heart rate, forearm blood flow, and forearm vascular resistance (FVR). Likewise, LBNP elicited comparable responses of these variables for blacks and whites. Baseline MSNA did not differ between blacks and whites, but whites demonstrated greater increases during LBNP (28 +/- 7 vs. 55 +/- 18%, 81 +/- 21 vs. 137 +/- 42%, 174 +/- 81 vs. 556 +/- 98% for -5, -15, and -40 mmHg LBNP, respectively; P < 0.001). Consistent with smaller increases in MSNA but similar FVR responses during LBNP, blacks demonstrated greater sympathetic vascular transduction (%FVR/%MSNA) than whites (0.95 +/- 0.07 vs. 0.82 +/- 0.07 U; 0.82 +/- 0.11 vs. 0.64 +/- 0.09 U; 0.95 +/- 0.37 vs. 0.35 +/- 0.09 U; P < 0.01). In summary, young whites demonstrate greater increases in MSNA during baroreceptor unloading than age-matched normotensive blacks. However, more importantly, for a given increase in MSNA, blacks demonstrate greater forearm vasoconstriction than whites. This finding may contribute to augmented blood pressure reactivity in blacks.

NASA Discipline Neuroscience↗
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 OSTI2019

Plasma kinetic effects on interfacial mix and burn rates in multispatial dimensions

The physics of mixing in plasmas is of great significance to inertial confinement fusion and high energy density laboratory experiments. Two- and three-dimensional (2D and 3D) particle-in-cell simulations with a binary collision model are used to explore kinetic effects arising during the mixing of plasma media. The applicability of the one-dimensional (1D) ambipolarity condition is evaluated in 2D and 3D simulations of a plasma interface with a sinusoidal perturbation. The 1D ambipolarity condition is found to remain valid in 2D and 3D, as electrons and ions flow together required for $J$ = 0. Simulations of perturbed interfaces show that diffusion-induced total pressure imbalance and hydroflows flatten fine interface structures and drive rapid atomic mix. The atomic mix rate from a structured interface is faster than the ~$\sqrt{\bar{t}}$ scaling obtained from 1D theory in the small-Knudsen-number limit. Plasma kinetic effects inhibit the growth of the Rayleigh-Taylor instability at small wavelengths and result in a nonmonotonic growth rate scaling with wavenumber k with a maximum at a low k value, much different from $\sqrt{Agk}$ (where A is the Atwood number and g is the gravitational constant) as expected in the absence of plasma kinetic effects. Simulations under plasma conditions relevant to MARBLE separated-reactant experiments on Omega and the NIF show kinetic modification of DT fusion reaction rates. With non-Maxwellian distributions and relative drifts between D and T ions, DT reactivity is higher than that inferred from rates using stationary Maxwellian distributions. Reactivity is additionally revealed to be reduced in the presence of finite-Knudsen-layer losses.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗