Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “CFD Modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

The molten salt tritium transport experiment: A pumped fluoride salt loop for hydrogen isotope experimentation

Molten salt reactors (MSRs) and fusion reactors propose to use molten salts as coolants and breeder blanket materials, respectively. Tritium, however, poses safety concerns in both reactor types due to its ability to permeate through reactor materials and potential for environmental release. This manuscript addresses the tritium transport phenomena in molten salts and presents the design and analysis of the Molten Salt Tritium Transport Experiment (MSTTE). MSTTE is a forced-convection fluoride salt loop intended to measure hydrogen isotope permeation through structural materials in a flowing salt system. In the first phase, MSTTE will use FLiNaK salt and deuterium as surrogates for FLiBe and tritium, with future plans to utilize tritium and FLiBe. MSTTE couples a Copenhagen Atomics pumped salt loop with an external test section that introduces hydrogen isotopes into the loop and measures transport phenomena. The Hydrogen Injection System (HIS) controls hydrogen isotope introduction into the molten salt loop. Here, the permeation test section measures the permeation rate through stainless steel tubing in contact with flowing salt. Computational fluid dynamics (CFD) analysis ensures fully developed salt flow in the permeation test section. MSTTE is modeled with MELCOR-TMAP to predict the permeation rate as a function of experimental variables such as source term, salt flow rate, and salt temperature. Results indicate that the source term is the only parameter with a significant effect on the permeation rate. Pressure drop analysis suggests that the loop should operate below 200 LPM to maintain a pressure drop below 200 kPa. Additionally, finite-element analysis assesses thermal stress during loop operation to ensure the experiment's safe design. MSTTE will provide semi-integral data on tritium transport phenomena in molten salts and serve as a testbed for advancing molten salt technology.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Effect of Thermal Expansion Coefficient, Viscosity and Melting Range in Simulation of PCM Embedded Heat Exchangers With and Without Fins

Phase change material heat exchangers (PCMHX) have an important role in integrating renewable energy systems. PCMHX can offer high storage density at various temperatures of interest and can be used for grid load shifting purposes. Numerical models enable engineers to estimate PCMHX performance for different design parameters and operating conditions. Modeling phase change phenomena is challenging due to the complex time-dependent nature of the process. The accuracy of models is highly sensitive to PCM thermo-physical properties. Thermal expansion coefficient (β), viscosity (μ) and melting temperature range (MR) of a PCM are important properties, especially when natural convection is not negligible. In PCMHX modeling, using less than accurate values for these properties can have significant impact on the simulation outcomes. These properties and discussions thereof are not readily available in the literature. This paper presents a brief review of the literature and a numerical study investigating the model sensitivity to the above-mentioned properties for PCMHX with and without fins. CFD is used to evaluate the charging (melting) phenomena. The study quantifies the impact of uncertainty in these properties on the melting rate and temperature profile. Constant wall temperature was considered as heat source with no heat loss to ambient. Here, the results show that β and μ has significant effect on the melting rate and evolution of the melting front. For a non-finned domain when comparing results for different published values of β and μ, the deviation in melting time can be up to 12.9% and 57.6% respectively. For high wall temperatures, change in melting range did not impact melting time. But when the wall temperature is reduced, up to 9.8% deviation in melting time is observed.

CFD↗
DOE OSTI2024

Development of prechamber enabled mixing-controlled combustion strategy for ultra-low methane emissions from lean burn natural gas engines

This numerical study explores the optimization of Prechamber Enabled Mixing-Controlled Combustion (PC-MCC) using natural gas in heavy-duty engines, aiming to enhance combustion efficiency to minimize methane slip and NOx emissions. The approach involves a prechamber ignition system, distinct from conventional spark ignition (SI) systems, to initiate combustion of direct injected natural gas. By leveraging the robust ignition characteristics of the prechamber, the PC-MCC method demonstrates significant potential in achieving efficient combustion akin to diesel engines but with lower greenhouse gas emissions. The research evaluates the effects of various geometric and operational parameters on the combustion process and emissions, including prechamber volume, nozzle diameter, direct injector (DI) geometry, and engine operating strategies. Computational Fluid Dynamics (CFD) simulations are utilized, focusing on a heavy-duty, single-cylinder engine modeled after the Caterpillar C9.3B engine. Key findings indicate that a prechamber volume of 3 cc, coupled with a nozzle diameter of 2.75 mm for two prechamber holes, strikes an optimal balance between combustion efficiency and emissions reduction. This configuration ensures robust combustion across a range of operating conditions while maintaining methane slip within targeted limits. Further investigation into DI geometry shows the significance of the injector umbrella angle and nozzle diameter in shaping the fuel-air mixing and combustion dynamics. An umbrella angle of 130° and a nozzle diameter of 300 microns are identified as optimal, promoting rapid and efficient combustion with minimized methane and NOx emissions. The study also investigates the impact of injection timing and pressure, highlighting their roles in controlling combustion timing and influencing emissions levels. Advanced injection timing is found to be crucial in achieving the desired low methane slip, whereas retarded injection timing assists to reduce NOx emissions while having a slight increase in methane emissions. Operating strategies incorporating various levels of Exhaust Gas Recirculation (EGR) are assessed for their effectiveness in further reducing emissions. The research demonstrates that a judicious combination of internal hot EGR and careful calibration of DI pressure and SOI timing can achieve significant reductions in NOx emissions while keeping methane slip under control. Specifically, an internal EGR level of 15%–25%, combined with DI pressures of 200–300 bar and injection timings at or after top dead center, is recommended. These findings contribute valuable insights into the development of advanced combustion techniques for natural gas engines, offering a viable pathway to reduce methane slip without compromising engine efficiency or performance. The PC-MCC system presents a promising solution for the future of heavy-duty natural gas engine technology to reduce methane emissions.

Nsaif, Osama↗
DOE OSTI2021

Enabling Load Following Capability in the Transatomic Power MSR

This project is dedicated towards designing a fuel processing system that enables liquid-fueled molten salt reactors (MSR) to load follow by removing the dissolved xenon in the fuel salt. As one of the Gen-IV nuclear reactor concepts, the molten salt reactor receives increasing development interests in the recent years. One distinguishing feature of the liquid-fueled molten salt reactor is its improved ability to operate in a load-following mode by including the unique online fission product removal system. Load-following means that the reactor changes its power output based on the demand on the grid. Most of the current operating nuclear reactors have limited load-following ability and operate as the base load on the grid. Due to the rapid increase of solar energy, the requirement on load-following capacity is significantly increased because of the varying power output of the solar panels, yet the traditional load-following capacity is expected to decrease as the decarbonization of the grid continues. Therefore, the ability to perform load-following operation for the nuclear reactors will greatly enhance the resilience of the grid and make nuclear energy more economically competitive. This load-following feature is included in many commercial molten salt reactor designs, such as the designs by Transatomic Power, Terrestrial Energy, and Flibe Energy. Unfortunately, detailed analysis of the fuel processing system for commercial scale MSRs is still lacking, as well as how the fuel processing quantitively impacts the load-following operation. Moreover, experimental data for many of the underlying physics of fuel processing is limited. This project aims to pave the way for the fuel processing technology to advance to the commercial stage by performing combined experimental and simulation research. During the project period, four interconnected aspects of the development of the fuel processing system in liquid-fueled molten salt reactors are investigated. These aspects are the simulation and analysis of the fission product removal system, the fuel cycle simulation, the coupled reactor neutronics and thermal hydraulics transient simulation, and the gaseous fission product removal experiment. Multiphase CFD simulations are performed for components of the processing systems, and simplified air-water experiments are carried out to provide validation data. It is concluded that the CFD simulation can satisfactorily predict the system level performance of the components, and engineering models are constructed based on this success. Fuel cycle analysis is performed for two representative MSR design, the MSBR and the Transatomic Power MSR. Open-source code SaltProc is developed to incorporate the unique fuel processing system of the MSRs. It is concluded that the removal of xenon is essential for load-following operation in thermal spectrum MSR and Molten Salt Breeder Reactor. For the Transatomic Power MSR, the xenon poisoning effect is negligible due to its relatively fast neutron spectrum, though the overall fuel cycle economics still benefits from the removal of xenon. Coupled reactor neutronics and thermal hydraulics transient simulation is performed specifically for the Transatomic Power MSR. It is concluded that the reactor core design could perform power ramping fast enough to satisfy load-following operation. Combining the findings from each aspect, it is concluded that the load-following operation of a thermal neutron MSR is dependent upon the removal of xenon, which could be achieved for a commercial sized reactor using continuous inert gas sparging in a separate system with reasonable dimensions.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Summary of DASSH-CFD Inter-Assembly Heat Transfer Comparison

Fast reactor core design requires integration of detailed power distributions obtained from neutronics models and thermal fluids calculations to determine temperature distributions and flow conditions. During the design phase, calculations may be repeated as specifications change and computational expense must be balanced with solution accuracy, creating an opportunity for the use of intermediate-fidelity codes. To that end, the Ducted Assembly Steady-State Heat Transfer Software (DASSH) has been developed at Argonne National Laboratory to perform full-core, steady state coolant and fuel pin temperature calculations in reactors consisting of hexagonal, ducted assemblies. DASSH is a first-order, finite volume subchannel code based on methodologies initially introduced in the legacy thermal hydraulics codes SUPERENERGY-2 and SE2-ANL. It can work coupled with the Argonne Reactor Computational (ARC) code suite to obtain neutron and gamma power distributions based on neutron transport solutions and features many improvements and new capabilities compared to its predecessors. DASSH calculates temperatures by balancing energy between subchannels and relies on correlations to account for flow distribution and coolant mixing.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

Thermal Modeling of Advanced Test Reactor Fuel in a Generalized Dry Storage System Under Hypothetical Accident Conditions

Work was performed to investigate the thermal behavior of a sealed aluminum-clad spent nuclear fuel (ASNF) dry storage configuration under hypothetical accident conditions during long-term storage. The considered system consists of a concrete dry storage overpack; a welded over canister, backfilled with argon; and nine Department of Energy (DOE) standardized canisters (DOESCs) loaded with ASNF, backfilled with helium gas. One technical concern associated with the long-term storage of ASNF includes radiolytic hydrogen generation. The yield of large quantities of hydrogen could lead to DOESC over pressurization or a flammable internal atmosphere. While flammability concerns can be resolved by preventing the ingress of oxygen, the canister pressure is partially controlled by the atmosphere temperature. This memorandum summarizes the results of modeling and simulation (M&S) work completed with Star CCM+ (a computational fluid dynamics [CFD] software) considering the thermal effects of a fully engulfing flame and another scenario with the vent ports of the concrete dry storage overpack completely blocked. The goal was to determine the bounding maximum temperature of the DOESC internal atmosphere to evaluate a worst-case pressure scenario, as well as the critical vent port blockage durations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

SHOCK TUBE IGNITION STUDIES OF RENEWABLE DIESEL FUELS FOR MEDIUM AND HEAVY-DUTY TRANSPORTATION

Currently extensive research on alternative fuels is being conducted due to their increasing demand to reduce greenhouse emissions. One renewable fuel studied in this work is dimethyl ether (DME) blended with propane(C3H8) as a potential mixture for heavy-duty engines used in semi-trucks. The blend has the potential to drastically reduce particulate and greenhouse gas emissions compared to a conventional diesel engine operating under similar conditions. To develop the use of mixture, one must conduct detailed conceptual and simulation studies before progressing to detail studies in CFD, engine modifications, and live testing. For simulations, accurate high-fidelity chemical kinetic models are necessary. However, the validity of the chemical kinetic mechanism for operating conditions of a heavy-duty mixing-controlled compression (MCCI) engine was widely unknown until recent work presented here and published. In this work, we studied the ignition of DME and propane blends in a shock tube under MCCI engine conditions. Ignition delay time (IDT) gathered behind the reflected shock for DME-propane mixtures for heavy-duty compression ignition (CI) engine parameters. Testing was conducted for undiluted varieties spanning from temperatures of 700 to 1100 K at pressures ranging from 55 to 84 bar for various blends (100% CH3OCH3, 100% C3H8, 60% CH3OCH3/ 40% C3H8) of DME and propane were combusted in synthetic air (21% O2/ 79% N2). Several experiments were conducted at higher pressures (90-120 bar) to improve the model performance and accuracy. The ignition delay times (IDTs) were compared to recent mechanisms, including Aramco3.0, NUIG, and Dames et al. A common trend among the mechanisms was overpredicted experimental IDTs. Further studies were conducted by a sensitivity analysis using the Dames et al. model, and critical reactions sensitive to IDTs of DME-propane mixture near 60 bar are outlined. Chemical analysis was conducted on the NTC region to explain chemical kinetics which is critical for developing MCCI heavy duty engines.

Mohammed, Zuhayr Pasha↗
DOE OSTI2026

Two-Way Coupled CFD/DEM Simulation of Electromagnetic Heating in a Fluidized Bed Reactor - Extended Abstract

A set of tools has been developed to couple a multiphase reacting computational fluid dynamics solver (MFIX) with two computational electromagnetic solvers (COMSOL and Elmer). In contrast to several recent investigations, which calculate a single static EM field prior to the CFD calculation, in the present work the electromagnetic field is periodically updated during the simulation. This enables the EM field and associated EM heating predictions to reflect the effects of evolving material properties (temperature, composition, size etc.) and the evolving distribution of material in the system. The present work uses effective medium theory (EMT) for the electromagnetic field calculations since the goal is to apply the tool for industrial scale reacting systems. Predictions using MFIX/COMSOL and MFIX/Elmer are compared to in bed thermocouple measurements from a fluidized bed experiment and recent simulation predictions and recent experiments performed at NETL. The NETL experiments also include packed bed operating conditions. Coupled MFIX/COMSOL simulations at these conditions are compared with similar simulations using only COMSOL to measure errors due to the solver coupling approach.

CFD↗
DOE OSTI2026

Development of a coupled experimental–computational approach for engineering optimization of spout-fluidized bed particle coating systems

The design of spout-fluidized bed (SFB) coating systems for nuclear particle fuels typically relies on trial-and-error processes, comprising iterative and time-consuming coating deposition experiments and post-deposition characterization. At an engineering scale, this approach to guided SFB system design is inefficient, highlighting the need for streamlined experimental methodologies which can correlate fluidization conditions to downstream coating outcomes. In this study, we combine time-resolved particle image velocimetry (PIV) with CFD–DEM simulations to benchmark hydrodynamic behavior in a 3D spout-fluidized bed. By exploiting easily accessible optical measurements of particle motion at the bed wall and within the spouting region, we obtain quantitative velocity fields that can be directly compared with model predictions of the occluded bed region, without resorting to complex imaging and characterization techniques such as X-ray or magnetic resonance tomography. Experimental benchmarking reveals strong agreement between CFD–DEM and PIV in the spout and annulus regions, while discrepancies near the wall highlight areas for future model development. Here, the proposed integrated experimental–numerical framework will enable a direct connection between measured variables and numerically predicted fluidization performance of dense, surrogate nuclear particle fuel feedstock such that experimental SFB component design can be rapidly evaluated, informing design decisions for nozzle geometry and operating conditions. Future work will extend this framework by correlating quantified fluidization metrics across nozzle geometries and operating conditions with the resulting coating morphology, microstructure, and uniformity. Establishing these correlations will enable predictive links between hydrodynamic performance and coating quality, providing a rational, scalable basis for optimizing SFB design prior to coating deposition.

CFD/DEM↗
DOE OSTI2025

Videos, photos, and AI-derived grain size data associated with “High-throughput AI Video Surveys Enable Reproducible Multiscale Sediment Size Mapping, with Implications for Hydrobiogeochemical Parameterization”

NOTE: The manuscript associated with this data package is currently in review. The data may be revised based on reviewer feedback. Upon manuscript acceptance, this data package will be updated with the final dataset and additional metadata. This data package is associated with the manuscript “High-throughput AI Video Surveys Enable Reproducible Multiscale Sediment Size Mapping, with Implications for Hydrobiogeochemical Parameterization” under review. This data package includes five data types: 1) raw photos and videos from drone survey and walking smartphone surveys; 2) images derived from raw videos; 3) manual labeling of reference scales; 4) metadata for all images and photo resolution derived from artificial intelligence (AI) models or manual labels, 5) grain size data obtained from AI models for all photos, 6) metadata and grain size data after quality control, 7) summaries of sample efficiency for all data, and 8) computational fluid dynamics (CFD) data used to support hydro-biogeochemical (HBGC) parameter estimation. Such data is used to 1) demonstrate significant improvements in accuracy, efficiency, and quality control for grain size data collection with the help of AI models, 2) study the spatial heterogeneity of grain size and observation reproducibility based on tens of thousands of data points generated by the AI models, and 3) evaluate the impacts of grain size heterogeneity on key HBGC parameters across sediment-to-reach and hourly-to-yearly scales. In particular, the data package contains 116 folders and 179696 files. The files include 41 videos in .mov format, 64047 photos in .jpg format, 13541 video-derived photos in .png format, 12747 segmentation mask data in .tif format, 12747 segmentation data in .json format, 24771 .csv files that with metadata and grain size for each individual photo as well as water depth and velocity data from CFD and observation, 51791 .txt files of raw AI predicted labels, and 11 flight record data in .srt format. The summary for all metadata and grain size statistics information is included in “Scales_V3_NG.csv” and “Statistics_V3_NG.csv”. The summary for data that pass data quality control (QC) level 0-2 is included in “QCStatistics_V3_NG.csv”. The QC level 0 represents photos whose photo resolution is positive, excluding photos that miss reference scale. The QC level 1 means reference scale circularity uncertainty is less than 5% for smartphone images while representing photo resolution is larger than 0.44 mm/pixel for drone images. The QC level 2 means excluding photos whose grain number is less than 100, a minimum number of grains recommended by classic literature. The summary for each video’s name, length, frame rates, survey area, grain number, survey efficiency, etc. can be found in “QCSummary_V3_NG.csv”. The summary for site name, GPS coordinates, and number of images at each site can be found in “SitesSummary_V3_*.csv” files. Overall computational efficiency summary is reported in Table 4 of accompanying manuscript. Additionally, the nitrate concentration data used in this work was downloaded from an existing dataset published on ESS-DIVE (Boat-Dragged Sensor Hanford Reach.csv; Conner A. et al., 2020). We thank the United States Forest Service, Washington Department of Fish and Wildlife, Washington Department of Natural Resources, Cowiche Canyon Conservatory, Port of Benton, and the Confederated Tribes and Bands of the Yakama Nation for access to field locations where the data were collected. We also thank the Yakama Nation Tribal Council and Yakama Nation Fisheries for working with us to facilitate data collection and optimization of data usage according to their values and worldview.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Temperature Homogeneity under Selective and Localized Microwave Heating in Structured Flow Reactors

Selective heating of different phases of multiphase systems via microwaves can result in energy savings and suppression of side reactions. However, materials properties and operating conditions that maximize temperature gradients are poorly understood. Here we utilize computational fluid dynamics (CFD) computations and temperature measurements in structured flow reactors (monoliths) in a monomodal microwave cavity to assess the temperature difference between the walls and the fluid and develop a simple lumped model to estimate when temperature gradients exist. We also explore the material’s thermal and electrical properties of structured reactors for isothermal catalyst conditions. We propose that CFD simulations can be used as a nonintrusive, predictive tool of temperature homogeneity. Importantly, we demonstrate that localized heating in the bed under several conditions rather than selective heating is responsible for the selectivity enhancement. Furthermore, our results indicate that structured beds made of high thermal conductivity materials avoid arcing and enable temperature homogeneity and low electrical conductivity materials allow microwaves to penetrate the domain.

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

Impacts of Biomass Feedstock Pre-Processing on Heat and Mass Transfer During Pyrolysis Using X-Ray Computed Tomography and Multiscale Modeling

Knowledge of the transport properties of biomass particles such as porosity, tortuosity, and permeability is paramount for high-fidelity modeling of biomass pyrolysis due to the heat and mass transfer limitations imposed by particle microstructure. X-ray computed tomography (XCT) is a non-destructive imaging method that enables full 3D reconstructions of the biomass particle microstructure with high resolution, permitting direct calculation of porosity, tortuosity, and permeability from real particle geometries. In this study, XCT imaging revealed the 3D microstructures of particles and chars from pyrolytic conversion of cylindrically cut or milled/pelletized loblolly pine samples. The porosity, tortuosity, and permeability were calculated directly from the XCT geometries via open-source microstructural analysis tool MATBOX+TauFactor (https://github.com/NREL/MATBOX_Microstructure_analysis_toolbox) and computational fluid dynamics (CFD) simulations using our solver, Mesoflow (https://github.com/NREL/mesoflow). These properties were used in a reactor scale model developed in COMSOL of the single particle reactor at NREL to investigate the impact of feedstock pre-processing on biomass conversion during pyrolysis with rigorous experimental validation.

biomass↗
DOE OSTI2021

CFD analysis of a solid pin-fueled small modular fluoride salt-cooled reactor

A computational fluid dynamics (CFD) study is performed for the pre-conceptual solid pin-fueled small modular fluoride salt-cooled reactor developed by the Oak Ridge National Laboratory (ORNL). The solid and fluid regions in a 1/12th section of a fuel assembly are modeled. The solid region includes fuel and non-fuel pins, with heat generation in the active region of the fuel pins. Two different power profiles, uniform and center-peaked, are considered. The FLiBe coolant flows from the bottom to the top of the core, parallel to the bank of fuel and non-fuel pins. The k-ω shear-stress transport (SST) model is chosen as the baseline turbulence model. The effects of grid refinement, inlet turbulence specification, and turbulence models on the temperature and pressure drop predictions are studied. Turbulence model sensitivity is investigated by comparing the results from the k-ω SST model with other two-equation models (k-ω baseline or BSL and k-ε realizable) as well as with anisotropic Reynolds stress transport models (linear pressure-strain and stress-BSL). Here the results show that the choice of turbulence model has a significant impact on the pin temperature and bundle pressure drop predictions. The results from the baseline turbulence model show good agreement between the fuel pin temperatures and bundle pressure drop values predicted using a subchannel model previously developed by the authors of the present study.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Computational Fluid Dynamics (CFD) Simulations of Taylor Bubbles in Vertical and Inclined Pipes with Upward and Downward Liquid Flow

Summary Two-phase flow is a common occurrence in pipes of oil and gas developments. Current predictive tools are based on the mechanistic two-fluid model, which requires the use of closure relations to predict integral flow parameters such as liquid holdup (or void fraction) and pressure gradient. However, these closure relations carry the highest uncertainties in the model. In particular, significant discrepancies have been found between experimental data and closure relations for the Taylor bubble velocity in slug flow, which has been determined to strongly affect the mechanistic model predictions (Lizarraga-García 2016). In this work, we study the behavior of Taylor bubbles in vertical and inclined pipes with upward and downward flow using a validated 3D computational fluid dynamics (CFD) approach with level set method implemented in a commercial code. A total of 56 cases are simulated, covering a wide range of fluid properties, pipe diameters, and inclination angles: Eo ∈ [10, 700]; Mo ∈ [1×10–6, 5×103]; ReSL ∈ [–40, 10]; θ ∈ [5°, 90°]. For bubbles in vertical upward flows, the simulated distribution parameter, C0, is successfully compared with an existing model. However, the C0 values of downward and inclined slug flows where the bubble becomes asymmetric are shown to be significantly different from their respective vertical upward flow values, and no current model exists for the fluids simulated here. The main contributions of this work are (1) the relatively large 3D numerical database generated for this type of flow, (2) the study of the asymmetric nature of inclined and some vertical downward slug flows, and (3) the analysis of its impact on the distribution parameter, C0.

Engineering↗
DOE OSTI2020

A status review on the thermal stratification modeling methods for Sodium-cooled Fast Reactors

The thermal stratification phenomenon plays a crucial role in the safety of various nuclear systems, including the Gen-III + Light Water Reactors (LWR) and the Gen-IV reactors. The phenomenon is of particular importance for the pool-type Sodium-cooled Fast Reactors (SFRs) because it may cause neutronic and thermal-hydraulic in- stabilities in the reactor core, or lead to damages of both the reactor vessel and in-vessel components due to the growth of thermal fatigue cracking. More significantly, thermal stratification could impede the establishment of the natural circulation during accidental scenarios and introduce uncertainties to the core safety of SFRs. Efforts for modeling of the thermal stratification in SFRs have been made for decades to prevent or mitigate the damage caused by the phenomenon. This paper gives a review of the advances that have been made in recent 10 years on the computational modeling methods for thermal stratification phenomenon in SFRs. These methods can be generally drawn into two categories. The first one is the system-level methods which provide fast-running but approximate calculations, and the second one is the CFD methods which provide high-resolution calculations at high computational expense. Finally, after introducing the efforts that have been made to improve the one- dimensional (1-D) models, the paper envisioned the possible research directions that could be pursued to enhance the modeling of thermal stratification in the near future.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Experimental measurements of gas-particle flows in large-scale strippers

Measurements of gas–solid particle flow in a large-scale stripper unit are reported. The experiments target industrial (pilot) scale measurements of gas-particle flow that can be used to validate numerical methods for modeling multiphase flows, such as computational fluid dynamics coupled to the discrete element method (CFD-DEM). Specifically, experiments were performed with Geldart Group B (533 µm) glass beads in a stripper unit. In general, previous CFD-DEM validation datasets were performed in benchtop experimental systems and limited to O(10 5 ) particles. Here, in this work, the number of particles in the isolated section of the stripper is O(10 7 ) for all tests. Measurements are reported for a ~2 m tall section of the ~1 m diameter stripper. Radial mass flux profiles measured at five axial positions and three axial pressure drop measurements are reported for six operating conditions.

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

Parametric Evaluation of RDE Inlet Performance using a Reduced-Order Model

This paper presents the results of a parametric study on the effects of inlet restriction, along with a number of other controlling parameters, using a previously developed reduced-order RDE inlet model. The results of this study will be used to inform subsequent high-fidelity CFD simulations and/or selection of experimental operating conditions which may be conducive to maximizing performance.

Bedick, Clinton↗
DOE OSTI2025

An adaptive, data-driven multiscale approach for dense granular flows

The accuracy of coarse-grained continuum models of dense granular flows is limited by the lack of high-fidelity closure models for granular rheology. One approach to addressing this issue, referred to as the hierarchical multiscale method, is to use a high-fidelity fine-grained model to compute the closure terms needed by the coarse-grained model. The difficulty with this approach is that the overall model can become computationally intractable due to the high computational cost of the high-fidelity model. In this work, we describe a multiscale modeling approach for dense granular flows that utilizes neural networks trained using high-fidelity discrete element method (DEM) simulations to approximate the constitutive granular rheology for a continuum incompressible flow model. Our approach leverages an ensemble of neural networks to estimate predictive uncertainty that allows us to determine whether the rheology at a given point is accurately represented by the neural network model. Additional DEM simulations are only performed when needed, minimizing the number of additional DEM simulations required when updating the rheology. This adaptive coupling significantly reduces the overall computational cost of the approach while controlling the error. In addition, the neural networks are customized to learn regularized rheological behavior to ensure well-posedness of the continuum solution. We first validate the approach using two-dimensional steady-state and decelerating inclined flows. We then demonstrate the efficiency of our approach by modeling three-dimensional sub-aerial granular column collapse for varying initial column aspect ratios, where our multiscale method compares well with the computationally expensive computational fluid dynamics (CFD)-DEM simulation.

Dense granular flows↗