Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “thermal stratification”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Microbial community dynamics and cycling of plutonium and iron in a seasonally stratified and radiologically contaminated pond

Plutonium (Pu) cycling and mobility in the environment can be impacted by the iron cycle and microbial community dynamics. We investigated the spatial and temporal changes of the microbiome in an iron (Fe)-rich, plutonium-contaminated, monomictic reservoir (Pond B, Savannah River Site, South Carolina, USA). The microbial community composition varied with depth during seasonal thermal stratification and was strongly correlated with redox. During stratification, Fe(II) oxidizers (e.g., Ferrovum, Rhodoferax, Chlorobium) were most abundant in the hypoxic/anoxic zones, while Fe(III) reducers (e.g., Geothrix, Geobacter) dominated the deep, anoxic zone. Sulfate reducers and methanogens were present in the anoxic layer, likely contributing to iron and plutonium cycling. Multinomial regression of predicted functions/pathways identified metabolisms highly associated with stratification (within the top 5%), including iron reduction, methanogenesis, C1 compound utilization, fermentation, and aromatic compound degradation. Two sediment cores collected at the Inlet and Outlet of the pond were dominated by putative fermenters and organic matter (OM) degraders. Overall, microbiome analyses revealed the potential for three microbial impacts on the plutonium and iron biogeochemical cycles: (1) plutonium bioaccumulation throughout the water column, (2) Pu–Fe-OM-aggregate formation by Fe(II) oxidizers under microaerophilic/aerobic conditions, and (3) Pu–Fe-OM-aggregate or sediment reductive dissolution and organic matter degradation in the deep, anoxic waters.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Performance of high-order Godunov-type methods in simulations of astrophysical low Mach number flows

High-order Godunov methods for gas dynamics have become a standard tool for simulating different classes of astrophysical flows. Their accuracy is mostly determined by the spatial interpolant used to reconstruct the pair of Riemann states at cell interfaces and by the Riemann solver that computes the interface fluxes. In most Godunov-type methods, these two steps can be treated independently, so that many different schemes can in principle be built from the same numerical framework. Because astrophysical simulations often test out the limits of what is feasible with the computational resources available, it is essential to find the scheme that produces the numerical solution with the desired accuracy at the lowest computational cost. However, establishing the best combination of numerical options in a Godunov-type method to be used for simulating a complex hydrodynamic problem is a nontrivial task. In fact, formally more accurate schemes do not always outperform simpler and more diffusive methods, especially if sharp gradients are present in the flow. For this work, we used our fully compressible Seven-League Hydro (SLH) code to test the accuracy of six reconstruction methods and three approximate Riemann solvers on two- and three-dimensional (2D and 3D) problems involving subsonic flows only. We considered Mach numbers in the range from 10 −3 to 10 −1 , which are characteristic of many stellar and geophysical flows. In particular, we considered a well-posed, 2D, Kelvin–Helmholtz instability problem and a 3D turbulent convection zone that excites internal gravity waves in an overlying stable layer. Although the different combinations of numerical methods converge to the same solution with increasing grid resolution for most of the quantities analyzed here, we find that (i) there is a spread of almost four orders of magnitude in computational cost per fixed accuracy between the methods tested in this study, with the most performant method being a combination of a low-dissipation Riemann solver and a sextic reconstruction scheme; (ii) the low-dissipation solver always outperforms conventional Riemann solvers on a fixed grid when the reconstruction scheme is kept the same; (iii) in simulations of turbulent flows, increasing the order of spatial reconstruction reduces the characteristic dissipation length scale achieved on a given grid even if the overall scheme is only second order accurate; (iv) reconstruction methods based on slope-limiting techniques tend to generate artificial, high-frequency acoustic waves during the evolution of the flow; and (v) unlimited reconstruction methods introduce oscillations in the thermal stratification near the convective boundary, where the entropy gradient is steep.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Modeling Offshore Wind Farm Performance in Coastal Low-Level Jets Using Coupled Mesoscale-Microscale Large Eddy Simulations

Accurately predicting wind farm reliability under complex offshore atmospheric conditions remains a key challenge, particularly during noncanonical meteorological events such as coastal low-level jets (LLJs). LLJs, characterized by strong nonmonotonic vertical shear and directional veer, depart significantly from the simplified inflow assumptions embedded in conventional design standards, low-fidelity engineering models, and microscale large eddy simulations of the atmospheric boundary layer. In this work, we use the virtual wind farm framework—an exascale, graphics processing unit–accelerated large eddy simulation platform coupled with high-fidelity aeroservoelastic turbine models and advanced mesoscale-microscale coupling via the ExaWind software stack—to investigate turbine responses under realistic LLJ forcing. Simulations are performed over the U.S. North Atlantic offshore domain with the use of meteorological inputs from New York State Energy Research and Development Authority buoy data, focusing on a representative LLJ case impacting the International Energy Agency 15 MW reference turbine. Our results show that LLJs can cause up to 50% power deficits in downstream turbine rows and significantly amplify low-speed shaft and tower loads through nonlinear coupling between complex inflow characteristics and turbine structural dynamics. Two primary mechanisms drive these load amplifications: (1) unique LLJ inflow features—including veer and vertical/lateral shear—and (2) the downstream evolution of the flow under stable thermal stratification, which suppresses turbulence mixing and alters wake recovery. These mechanisms produce streamwise variations in turbine loading not captured by standard hub height–based metrics or existing design load case (DLC) definitions. This study highlights the critical role of rotor-scale flow gradients in driving fatigue and system-level aeroelastic responses, challenging current DLC and control strategies. We advocate the integration of full-flow field, environment-aware wind inputs into load modeling and control algorithms. By leveraging exascale computing to resolve mesoscale-microscale coupling, this work lays the groundwork for next-generation offshore wind turbine design and operation in meteorologically complex marine environments.

17 WIND ENERGY↗
DOE OSTI2021

Measured influence of overhead HVAC on exposure to airborne contaminants from simulated speaking in a meeting and a classroom

Tracer gas experiments were conducted in a 158 m 3 room with overhead supply diffusers to study dispersion of contaminants from simulated speaking in physically distanced meeting and classroom configurations. The room was contained within a 237 m 3 cell with open plenum return to the HVAC system. Heated manikins at desks and a researcher operating the tracer release apparatus presented 8-9 thermal plumes. Experiments were conducted under conditions of no forced air and neutral, cooled, or heated air supplied at 980-1100 cmh, and with/out 20% outdoor air. CO 2 was released at the head of one manikin in each experiment to simulate small (< 5 µm diameter) respiratory aerosols. The metric of exposure relative to perfectly mixed (ERM) is introduced to quantify impacts, based on measurements at manikin heads and at three heights in the center and corners of the room. Chilled or neutral supply air provided good mixing with ERMs close to one. Thermal stratification during heating produced higher ERMs at most manikins: 25% were ≥ 2.5 and the highest were > 5× perfectly mixed conditions. Operation of two within-zone air cleaners together moving ≥400 cmh vertically in the room provided enough mixing to mitigate elevated exposure variations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Assimilation of a Coordinated Fleet of Uncrewed Aircraft System Observations in Complex Terrain: Observing System Experiments

Uncrewed aircraft system (UAS) observations from the Lower Atmospheric Profiling Studies at Elevation–A Remotely-Piloted Aircraft Team Experiment (LAPSE-RATE) field campaign were assimilated into a high-resolution configuration of the Weather Research and Forecasting (WRF) Model. The impact of assimilating targeted UAS observations in addition to surface observations was compared to that obtained when assimilating surface observations alone using observing system experiments (OSEs) for a terrain-driven flow case and a convection initiation (CI) case observed within Colorado’s San Luis Valley (SLV). The assimilation of UAS observations in addition to surface observations results in a clear increase in skill for both flow regimes over that obtained when assimilating surface observations alone. For the terrain-driven flow case, the UAS observations improved the representation of thermal stratification across the northern SLV, which produced stronger upvalley flow over the eastern half of the SLV that better matched the observations. For the CI case, the UAS observations improved the representation of the pre-convective environment by reducing dry biases across the SLV and over the surrounding terrain. This led to earlier CI and more organized convection over the foothills that spilled outflows into the SLV, ultimately helping to increase low-level convergence and CI there. In addition, the importance of UAS capturing an outflow that originated over the Sangre de Cristo Mountains and triggered CI is discussed. These outflows and subsequent CI were not well captured in the simulation that assimilated surface observations alone. We report that observations obtained with a fleet of UAS are shown to notably improve high-resolution analyses and short-term predictions of two very different mesogamma-scale weather events.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

CFD modeling of pre-spark heat release in a boosted direct-injection spark-ignition engine

Accurate predictions of low-temperature heat release (LTHR) are critical for modeling auto-ignition processes in internal combustion engines. While LTHR is typically obscured by deflagration, extremely late ignition phasing can lead to LTHR prior to the spark, a behavior known as pre-spark heat release (PSHR). In this research, PSHR in a boosted direct-injection spark-ignition engine was studied using 3-D computational fluid dynamics (CFD) and detailed chemical kinetics. The turbulent combustion was modeled via a hybrid approach that incorporates the G-equation model for tracking the turbulent flame front, and the well-stirred reactor model with detailed chemistry for assessing the low-temperature reactions in unburnt gas. Simulations were conducted using Co-Optima alkylate and E30 fuels at operating conditions characterized by different PSHR intensities. The predicted in-cylinder pressure and heat release rate were found to agree well with experiments. It was found the estimate of previous-cycle trapped residuals is of utmost importance for capturing PSHR correctly. A simulation best practice was developed which keeps the detailed chemistry solver active throughout the entire simulation, allowing to track the evolution of intermediate species from one cycle to the next. Following the validation, the dynamics of PSHR were analyzed in detail employing the pressure-temperature (P-T) trajectory framework. It was shown that PSHR correlated with the first-stage ignition delay of the fuel, hence showing close relation to the in-cylinder P-T trajectory and the chemical kinetics. Besides, it was indicated that LTHR is a self-limiting process that has the effect of attenuating the thermal stratification in the combustion chamber. Furthermore, it was observed the occurrence of PSHR caused the P-T trajectory of end-gas to overlap with the negative temperature coefficient region of the fuel’s ignition-delay maps. This effect was more significant in the fuel-rich regions where engine knock tendency would be generally higher, with potential implications on knock control and mitigation.

42 ENGINEERING↗
DOE OSTI2023

Fuel-air mixing in motored CFR engine at research octane number (RON) relevant condition

This paper presents a three-dimensional (3-D) computational fluid dynamics (CFD) study of a motored cooperative fuel research (CFR) engine at research octane number (RON) relevant condition. The boundary conditions for 3-D simulations were generated with a one-dimensional GT-Power model. For the first time in literature, a carburetor was added to a virtual CFR engine model with 3-D CFD. Therefore, the proposed setup can simulate the fuel and thermal stratifications inside the engine cylinder with realistic detail. The transient simulations in this work were performed within the Reynolds-averaged Navier-Stokes (RANS) framework with a Realizable k-ε turbulence model. Major conclusions from the present work are: (1) The in-cylinder flow of the CFR engine is swirl-dominated due to the existence of the intake valve shroud. (2) There is a significant amount of liquid droplets entering the cylinder during the intake stroke. The maximum instantaneous amount of liquid for 50% PRF 87 (containing 87% iso-octane and 13% n-heptane (v/v)) and 50% ethanol mixture is indicated to be around 26% of total injected fuel mass. (3) The heat of vaporization (HoV) of the fuel is responsible for creating both temperature and charge stratification inside the cylinder.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2024

System Level Analysis Software Validation with Argonne's THETA Experimental Facility (Final CRADA Report)

Under this CRADA, the Contractor will be generating and collecting sodium-based experimental data with modern instrumentation. The Contractor will work symbiotically with the Participant to generate SAS4A/SASSYS-1 and SAM computational models of the THETA experimental facility for validation of the system-level analysis codes utilized by Oklo, Inc. This effort will be specifically focused on the design basis event space scoped by Chapter 15 of the NRC’s Standard Review Plan (NUREG-0800) with an emphasis on natural circulation and thermal stratification phenomena in a prototypic sodium-cooled fast reactor (SFR).

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Environmentally Assisted Fatigue in Light Water Reactor Environment

This report summarizes the Environmentally Assisted Fatigue (EAF) research conducted at ANL under the US DOE Light Water Reactor Sustainability (LWRS) program. Starting from a rich background in theoretical and experimental EAF, ANL previously developed an approach to evaluate fatigue performance of reactor materials in light water reactor environments with the correction factor F en . The approach was based on a large body of experimental work performed at ANL and elsewhere, and was consistent with American Society of Mechanical Engineers (ASME)’s methodology governing the design and construction of reactor components. In recent years, the program was focused on component fatigue prediction and made several major and fundamental contributions in this area. These accomplishments help meet the needs identified by the industry concerning component level fatigue predictions in complex, transient conditions. The main contribution of the ANL program involved the development of a system-level model for estimating residual strain and life of nuclear reactor coolant system components under connected-system-thermal-mechanical boundary conditions. The goal was to predict the stress hotspots, strain residuals, strain amplitudes and the resulting fatigue lives. Thermal-mechanical stress analysis was performed considering thermal stratification and a design-basis reactor loading cycle. Based on the finite element (FE) model results, the strain residuals, strain amplitudes and resulting fatigue lives of reactor coolant system (RCS) components were predicted. The results show that some of the RCS components can have significantly different strain amplitudes, residual strain, and fatigue lives, despite having similar geometry and material. In addition, the simulated component-level strain profile can guide the selection of appropriate test inputs for conducting laboratory-scale EAF tests. Building upon the system-level model, ANL developed a digital twin (DT) framework to predict the structural states and associated fatigue life of components in real-time. This framework is a comprehensive system designed to predict the structural states and fatigue lives of reactor components. It includes multiple models and integrates artificial intelligence (AI), machine learning (ML), and FE based modeling tools to evaluate the structural states and fatigue lives.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Marine Hydrogen Demonstration

This report summarizes Phase 1 of a project involving the design of a Floating Hydrogen Production and Dispensing Barge destined for the Port of San Francisco (SF). The H 2 Barge is designed to produce renewable H 2 at the rate of ~ 530 kg/day, storing 512 kg of hydrogen at 517-bar, allowing fast refueling of hydrogen fuel cell vessels and land-side hydrogen delivery trailers for distribution into the nascent SF hydrogen ecosystem. The broader considerations that impacted the H2 Barge design are also described. An account is given of a new review process formulated by the United States Coast Guard (USCG) to review this first-of-its-kind maritime implementation of hydrogen technology. The immediate goals of the H 2 Barge Project are to 1) demonstrate the feasibility, viability and methods of hydrogen production, storage and fueling in a maritime context, 2) help shape (where needed) and navigate the required local, state and federal regulatory gauntlet and 3) catalyze a “green hydrogen ecosystem” (both marine and landside) with locally produced renewable hydrogen at the San Francisco waterfront. A summary is also given of the modeling and experimental activity of Phase 1 directed to the development of science-based refueling protocols for large marine Type IV 250-bar hydrogen tanks. Combined modeling and experimental studies are reported of the filling of large (28 kg) 250-bar Type IV hydrogen tanks of the type being deployed on early hydrogen ferries, such as the MV Sea Change. The primary question was to determine how such tanks can be successfully filled (state of charge greater than 97%) within 45 minutes without exceeding the 82 °C temperature limit historically set for such tanks. The studies show that a gas injector is needed avoid thermal stratification during filling which can result in potential hot spots. Pre-cooling of the hydrogen was found to be essential in most cases, as ambient conditions greatly affect the need for a pre-cooling to achieve the 45-minute fill time. Pre-cooling cannot be supplied by nearby water, such as that found in nature (bays, lakes, rivers, etc.) because pre-cooling cooling below 0 °C was found to be necessary to avoid excessive compression heating. The experimental results afforded a calibration of the engineering model (SOFIL) for these large 250-bar tanks, which now enables using SOFIL to predict volume-averaged hydrogen filling temperatures to an accuracy of +/- 2.7°C for these tanks. The model can therefore be used to evaluate potential scenarios for development of a standardized fueling methodology for ferries utilizing large Type-IV tanks of the type examined here.

08 HYDROGEN↗
DOE OSTI2023

Modeling of Liquid Hydrogen Tank Cooled with Para-Orthohydrogen Conversion

With the accelerating development of liquid-hydrogen storage facilities, the problem of boil-off hydrogen losses becomes very important. A promising method to reduce these losses is to utilize the endothermic para-orthohydrogen conversion of vented hydrogen, which can effectively decrease heat loads on a hydrogen tank. To model such a process, a hybrid computational model has been developed, based on the application of computational fluid dynamics for an ullage space, where knowledge of thermal stratification is important, and reduced-order models for other system elements. The simulation results for a spheroidal tank in selected conditions indicated a 10–25% reduction of boil-off losses due to cooling produced by vented hydrogen undergoing para-ortho-conversion.

08 HYDROGEN↗
DOE OSTI2025

Evaluating Chemical Kinetics Predictions for Propane Using 3-D and 0-D Models in a Boosted Spark-Ignited Engine

Propane has been shown to be a promising alternative fuel to reduce emissions while simultaneously achieving high efficiencies in medium- and heavy-duty engines. These high-power density applications require boosted engines which, combined with high compression ratio, can lead to auto-ignition and knock. While three-dimensional (3-D) computational fluid dynamics (CFD) models are often used for resolving the complex fluid flow in engines, these models can become computationally expensive when simulating detailed chemical kinetics. Likewise, zero-dimensional (0-D) models are computationally concise enough for kinetics development, but lack any flow-field information which governs the flame propagation processes in spark ignition (SI) engines. This work presents a comprehensive comparison between 3-D and 0-D closed cycle simulations at knocking conditions in a high compression ratio high stroke-to-bore ratio propane engine. In order to initialize the flow-field for the 3-D closed cycle (intake valve closing, (IVC) to exhaust valve opening, (EVO)) simulation, a motored multi-cycle 3-D model was run using Converge to create a map at IVC, reducing the computational time. The map allowed a non-homogeneous 3-D closed cycle simulation to be satisfactorily validated against experiments, while a homogeneous case using only the turbulence field mapping was also simulated, mimicking 0-D modeling. The 3-D simulations were used to prescribe the initial conditions (e.g., IVC thermodynamics, speciation, burn-rate profile) for a 2-zone 0-D SI engine model in Chemkin Pro for both cases. It was found that 2-zone 0-D modeling underpredicted the knock onset timing, likely due to the lack of thermal stratification in the unburned gas region. Future work will carry multi-zone 0-D modeling to capture the fuel auto-ignition in the unburned region.

Douvry-Rabjeau, Julien [Oakland University, Roches↗
DOE OSTI2024

A chemical kinetic analysis of knock propensity of methanol-to-gasoline fuel

Production of low carbon gasoline-like fuels such as methanol-to-gasoline (MTG) is a promising approach to achieve rapid greenhouse gas emission reduction of the transportation sector. Despite the fact that gasoline that meets the ASTM D4814 standard for automotive spark-ignition engine fuel can be readily produced from these processes, it is unclear how the composition of MTG may affect engine performance and emissions. Here, in this paper, a surrogate for an MTG is used to numerically study the effects of gasoline composition on knock propensity and on the sensitivity of knock to thermal and fuel stratification, to oxygen dilution and to nitric oxide from exhaust gas recirculation of residual gases. Simulations were performed in ANSYS CHEMKIN-PRO using a comprehensive chemical kinetic mechanism for gasoline surrogates, and results of the MTG surrogate were compared against those of a petroleum-based regular E10 gasoline, termed PACE-20. A premium-grade MTG fuel was also formulated by adding ethanol to the MTG surrogate, and results were compared against those of four premium-grade, gasoline-like fuels representative of future alternative gasoline formulations. Surrogates and mechanism were evaluated by comparison against experimental engine data, and the model showed high accuracy at stoichiometric conditions (mean absolute error of ignition timing equal to 1.46 crank angle degrees) but larger deviations at lean conditions (mean absolute error of ignition timing equal to 5.52 crank angle degrees). Despite the fact that the MTG surrogate has a RON 1.1 units higher than that of PACE-20, it may show higher knock propensity at medium temperature conditions due to a less intense NTC behavior. MTG autoignition was more temperature- and equivalence ratio-sensitive than that of PACE20, suggesting that MTG can benefit more from naturally-occurring thermal stratification or from induced fuel stratification of the end gas to mitigate knock intensity. The sensitivity of autoignition reactivity to oxygen dilution and to NO concentration was higher for MTG than for regular gasoline at medium loads, but the opposite trend was observed at high loads due to the effect of pressure on the low-temperature chemistry of regular gasoline. Approximately 14 % vol ethanol content was required to upgrade the octane rating of MTG from regular grade to premium grade. Adding 13.6 % vol ethanol made the fuel autoignition less sensitive to both oxygen dilution and NO content (ignition time varies approx. 17 % and 50 % less with oxygen dilution and NO addition, respectively, when adding ethanol at high engine loads).

02 PETROLEUM↗
DOE OSTI2024

Numerical Simulation and Experimental Comparison of System Analysis Module 1D Mixing Model for Cold Shock Transients in the Gallium Thermal-Hydraulic Mixing Facility

Abstract Liquid metals are being investigated as coolants in many advanced reactor designs because of their high thermal conductivity and effectiveness at high temperatures. However, they often pose challenges to reactor operation and safety because of the complex thermal mixing and stratification in the plenum of pool-type reactor designs. The advanced system analysis code System Analysis Module (SAM) currently under development at Argonne National Laboratory aims to develop and implement thermal mixing models to accurately capture these complex thermal fluid behaviors. In this study, the SAM thermal mixing model was compared against experimental data from the Gallium Thermal-Hydraulic Experiment facility, a scaled liquid metal test facility that uses gallium as a surrogate fluid to investigate the stratification and thermal mixing of low-Prandtl-number fluids in the upper plenum of a liquid metal-cooled reactor. Two cold shock transient cases were used: one with stable stratified flow (Ri = 32) and one with stronger thermal mixing (Ri = 0.5). The resultant temperatures were then compared with the experimental temperatures over the entire plenum to assess the ability of the mixing models to capture the thermal behavior and to better correspond mixing parameters to various flow scenarios. Generally, the zero-dimensional mixing model was more capable of capturing the bulk temperature of the component modeled assuming that an accurate mass flow rate was provided, but it was inherently unable to capture thermal gradients in space. The one-dimensional mixing model was capable of capturing that the thermal gradients provided accurate selection of the mixing coefficients. Further, the temperature at the outlet junction was compared over time for each of the mixing models with the recorded experimental temperature. The implemented mixing models demonstrated the ability to effectively capture the overall thermal behavior for stronger mixing scenarios but struggled with more stably stratified flows. It was found that a system analysis code's covering of the entire range of different operating conditions still remains a challenging task, and it is suggested that further model and closure improvements are necessary to accurately capture complex thermal mixing and stratification phenomena.

stratification↗
DOE OSTI2024

SAM Code Enhancements for Modeling of Liquid Metal-Cooled Fast Reactor Concepts

The SAM code is under development and supported by DOE-NE’s Nuclear Energy Advanced Modeling and Simulation program as a modern system-level modeling and simulation tool for advanced non-light water reactor safety analyses. These advanced reactor concepts incorporate novel and improved approaches to achieve safety and economic feasibility. This report summarizes two major efforts in addressing the modeling gaps in SAM for liquid-metal-cooled fast reactor (LMFR) applications, i.e. thermal mixing and stratification phenomena in large pools and corrosion-oxidation of components in flowing lead. A new one-dimensional model for thermal mixing and stratification effects in large pools and enclosures is developed and implemented. Thermal mixing and stratification occur when fluid enters a pool at a temperature different than the bulk fluid itself, a scenario often encountered during transients in pool-type systems. These phenomena are critical for the safety of reactors, impacting phenomena like natural circulation, which is essential for passive cooling. The improved model in SAM addresses limitations of state-of-the-art approaches by combining one-dimensional (1D) channels, representing the coolant jet flow, with lumped-parameter zero-dimensional (0D) pools, representing the rest of coolant in the tank. Energy exchange between the 1D jet and the 0D pools is based on heat transfer correlations calibrated against 3D simulations. It is verified that this model can handle various flow configurations, including hot jets in colder pools, cold jets in hotter pools, and the presence of features like ceilings, free surfaces, and obstacles. Additionally, validation against experimental data demonstrates the ability of the model to capture mixing and stratification effects in a wide range of conditions. The flexibility and improved accuracy of the new model make it a valuable tool for reactor safety analysis, allowing for the simulation of different geometries encountered in advanced reactors. A system-level corrosion modeling capability is developed and implemented in SAM to support Lead Fast Reactor (LFR) development. Although the initial focus of this capability will be on LFR application, this can later be expanded to include other liquid metals such as Lead-Bismuth Eutectic (LBE) and PbLi. This report summarizes the common corrosion mitigation strategies and outlines the progress on implementing and validating a corrosion-oxidation model in SAM. Verification and validation of the corrosion-oxidation portion of the model was performed using analytical solution and measured data from samples tested in the non-isothermal pumped lead loop at IPPE Obninsk. The iron transport and corrosion/precipitation portion of the model was assessed using an analytical model and measured corrosion depths from a natural convection lead loop experiment performed at CEA. It is demonstrated that the model implemented in SAM performed well in these assessments.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Machine Learning Assisted Safety Modeling and Analysis of Advanced Reactors

With the advances in computational power and numerical methods, analysts can now rely on first-principle simulations to predict ultra-fine details in a variety of applications. Advances in machine learning (ML) have produced algorithms that can now learn high-level abstractions via hierarchical models. This project aims to leverage advances in ML techniques and the available high-resolution simulation data to develop a novel modeling and simulation (M\&S) methodology for reactor safety analysis. While application-agnostic ML techniques are available, complex physics constraints need to be incorporated into ML techniques to build ML-based closures for computationally efficient predictive simulations. This project intends to develop a physics-guided data-driven multi-scale methodology for M\&S of advanced reactors. The project focuses on thermal fluid (T/F) phenomena, which play major roles in advanced reactor safety. Specifically, we propose a data-driven coarse-mesh turbulence model based on local flow features for the transient analysis of thermal mixing and stratification in a sodium-cooled fast reactor (SFR). The model has a coarse-mesh setup to ensure computational efficiency, while it is trained by fine-mesh computational fluid dynamics (CFD) data with Reynolds-averaged Navier-Stokes (RANS) turbulence model to ensure accuracy. Three different neural networks are developed and tested for loss-of-flow transients in the hot pool of SFR, i.e. the densely connected convolutional neural network (DCNN), long-short-term-memory network based on proper orthogonal decomposition (POD-LSTM), and the DCNN informed by LSTM (DCNN-LSTM). The performances of these three neural networks are evaluated based on baseline models. The DCNN-LSTM model has been chosen for further hyperparameter optimization. Furthermore, based on a simplified two-dimensional case, uncertainty quantification (UQ) of the developed ML-based closure are investigated with three methods, i.e. Monte Carlo dropout, deep ensemble, and Bayesian neural network. The developed ML-based turbulent viscosity closure relation based on deep ensemble is then integrated into the system analysis module SAM and serves as a term in the conservation equations. Such a SAM-ML based procedure guarantees that the obtained results are consistent with the physical constraints of the thermal-fluid system. The SAM-ML simulation on the same loss-of-flow transient showed comparable accuracy with the CFD simulation but with a much coarser mesh setup. Last but not least, the ML-based closure improvement with the support of higher-fidelity data from large eddy simulation (LES) is discussed. As a first step towards this direction, a baseline LES simulation is performed to obtain comparable data with RANS results. Based on the early results, future investigation on further improving the ML-based closure is discussed. We believe the developed approach that combines scientific machine learning with nuclear system analysis code can benefit the advanced reactor community as more accurate safety analyses will better characterize reactor safety margins and reduce licensing efforts.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

A liquid stratification model to predict failure in thermally damaged EBW detonators

In previous work, commercially available downward facing exploding bridgewire detonators (EBWs) were exposed to elevated temperatures. These detonators were then initiated using a firing set which discharged a high amplitude short duration electrical pulse into a thin gold bridgewire. Responses of the detonators were measured using photonic doppler velocimetry (PDV) and high-speed photography. A time delay of 4 μs between EBW initiation and first movement of an output flyer separated operable detonators from inoperable detonators or duds. Here, we propose a simple method to determine detonator operability from the calculated state of the detonator at the time the firing set is initiated. The failure criterion is based on the gap distance between the exploding bridgewire (EBW) and the adjacent initiating explosive within the detonator which is low-density pentaerythritol tetranitrate (PETN) that melts between 413-415 K (140-142 ºC). The gap forms as PETN melts and flows to the bottom of the input pellet. Melting of PETN is modeled thermodynamically as an energy sink using a normal distribution spread over a temperature range between the onset temperature of 413 K and the ending temperature of 415 K. The extent of the melt is determined from the average temperature of the PETN. The PETN liquid is assumed to occupy the interstitial gas volume in the lower part of the input pellet. The vacated volume from the relocated liquid forms the gap between the EBW and the PETN. The remaining sandwiched layer consists of solid PETN particles and gas filling interstitial volume. We predict that a threshold gap between 17-27 μm separates properly functioning detonators from duds.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2025

Entrainment, Detrainment, and Dilution of Dry and Moist Atmospheric Thermals

Here this study examines the entrainment, detrainment, and dilution of dry and moist (cloud) atmospheric thermals in large-eddy simulations. In a neutrally stable environment (with respect to dry dynamics), moist thermals have an increase in radius R with thermal height z t (α ≡ dR/dz t ) about 4 times smaller compared to dry thermals when density stratification is considered and ~2.4 times smaller without density stratification (i.e., applying the Boussinesq approximation). An analytic expression relating α to several dimensionless parameters is derived from the thermal impulse–circulation relation to clarify the factors impacting α. This expression shows that the difference in buoyancy structure between moist and dry thermals, with buoyancy concentrated in the central cores of moist thermals owing to latent heating, explains their smaller spreading rates. Individual contributions of entrainment and detrainment are analyzed using a direct parcel-based approach in the simulations. Moist thermals have similar fractional detrainment but much smaller fractional entrainment rates compared to dry thermals, consistent with the differences in α. Despite having smaller α, moist thermals are similarly dilute (quantified by a passive tracer) as dry thermals because of their greater mixing efficiency with the environment. Thus, moist thermals are substantially dilute but expand much less in size/volume as they rise compared to dry thermals. The α values for moist thermals in (dry) neutral and statically stable environments are similar, but fractional entrainment and especially detrainment rates are greater in the stable environment. Large detrainment rates are associated with a breakdown of the broader thermal vortex ring structure, especially with low environmental relative humidity, attributed in part to evaporation and buoyancy reversal.

54 ENVIRONMENTAL SCIENCES↗