Search NASA⌕ Search

SEARCH · Search NASA

Results for “natural circulation”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

FY24 Report on Water NSTF Testing: Lower Tank Inlet Piping Configuration

The following report serves as a summary of the accomplishments and testing results by the Natural convection Shutdown heat removal Test Facility (NSTF) experimental program over the past 12-month period. A major activity included reconfiguration of the chimney piping geometry which altered the discharge position into the tank from the original 50% tank height to a new lower position at a 10.9% tank height. This modification increases the volume of available coolant thus extending long-term operating capacity, however also results in a decreased liquid driving head which has the potential to reduce the natural circulation efficiency and decrease overall performance. Examination of the tradeoffs for these two configurations is an area of interest for RCCS designers and drove planned test activities this year. Ten matrix tests were performed in FY24, totaling 195 hours of heated operations and 12.9 MWh of electrical heating throughout the year, with eight classified as Accepted per NQA-1 and two classified as Trending. Testing prior to the facility reconfiguration examined the facility response to throttling at the tank inlet, demonstrating increased sensitivity of flow instabilities to throttling within the two-phase region when compared to throttling at the lower sensitivity single-phase inlet region. Testing after the chimney reconfiguration began with a baseline test at conditions of 80% inventory fill and prototypic thermal power input of 2.1 MW t . In addition to establishing a reference for nominal system behavior and performance, repeat testing with multiple subsequent runs demonstrated strong repeatability between tests performed at the same conditions in this new configuration. These tests also examined if influences would occur to system behavior with installation of new higher-resolution instrumentation within the upper chimney. In this critical region where boiling and flashing phenomena dominate, the test results provided confidence that the new instrumentation does not uncharacteristically influence the observed behaviors. An inventory parametric series was then initiated to examine system behavior with the new chimney configuration at six varying initial inventory levels, ranging from high 80% to low 20% fill. Two-phase oscillations with similar peak and mean flow rates were observed when comparing to the previous mid tank configuration. Generally, similar system response trends with inventory were also observed, such as two phase oscillations suppressing as inventories were lowered. However, in one absolute inventory comparison at the highest fill of 80%, flow oscillations saw a gradual growth over the 4 hours of two-phase operation in the new lower tank inlet configuration, opposite of the gradual dampening observed in the mid tank inlet configuration. This can be attributed to, in part, a greater hydrostatic head pressure above the two-phase discharge region where the boiling front is developed. Furthermore, the change resulted in greater sensitivity to liquid degassing phenomena during single-phase heating, causing loop instabilities to form which trigged moderate flow degradation during the period approaching saturation and boiling conditions. This behavior had been observed previously under some conditions but was a common occurrence in recent testing with the newer lower tank inlet configuration. Initial observations from these first data sets suggest an overall larger window of stability for the mid tank inlet configuration compared to the lower tank inlet. Lastly, the lowest inventory fill test was repeated over an extended testing window to examine depletion behavior. Natural circulation flow and effective heat removal performance were observed during most of the testing period; only after the tank became fully drained (0% fill) did flow stagnate and violent geysering events occur. This early observation confirms one relative advantage over the mid tank inlet configuration, which stagnated under comparable conditions at ~20% inventory remaining in the tank.

42 ENGINEERING↗
DOE OSTI2022

Analysis of Postulated Accident Scenarios for the Transformational Challenge Reactor

The Transformational Challenge Reactor (TCR) is a helium-cooled, yttrium-hydride-moderated reactor that was designed for the U.S. Department of Energy Office of Nuclear Energy. A key objective of the TCR was to employ advanced manufacturing techniques in a nuclear system and demonstrate their potential for revolutionizing the nuclear reactor design process. One purpose of the present work is to demonstrate the safety of the TCR under postulated accidents. Based on RELAP5-3D and COMSOL analyses, the TCR remained below all current safety limits and far below the expected failure limits for the core materials. Another purpose of this work is to provide useful insights and recommendations regarding the application of RELAP5-3D to gas-cooled or other advanced reactors. A novel approach was implemented for simultaneously modeling conduction and radiation in RELAP5-3D, which was found to provide reasonable predictions of radial core, vessel, and ex-vessel heat transfer during postulated events. A multicode approach was also applied, in which high-fidelity COMSOL calculations were used to tune the radial heat transfer parameters in RELAP5-3D. The tuned RELAP5-3D model demonstrated comparable peak temperature predictions as COMSOL, despite a coarse treatment of the core in RELAP5-3D consisting of only two lumped heat structures. This high-fidelity tuning approach enabled enhanced accuracy as well as minimal complexity within the RELAP5-3D model, even for complex fuel geometric designs as in the TCR. Finally, investigations were made into the potential for flow reversal during a pressurized loss-of-forced-flow event in the TCR. The TCR is designed with downward helium flow through the core during normal operation. The RELAP5-3D model predicted that this downward flow would persist, without flow reversal, up to several days after the circulator trip. This was attributed to natural circulation hysteresis effects as have been noted in similar thermofluidic systems. Although flow stagnation and eventual reversal did not lead to unsafe TCR conditions, interesting spatial effects were observed which may have safety relevance for other reactor system designs and coolant types that are designed for downward core flow during normal operation, warranting closer investigation of the flow reversal phenomenon.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

3D CFD Model Validation Using Benchmark Data of 1/16th Scaled VHTR Upper Plenum and Development of Wall Heat-Transfer Correlation For Laminar Flow

With support from the U.S. Department of Energy-Office of Nuclear Energy’s (DOE-NE’s) Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, an effort has been pursued to support high-temperature gas-cooled reactor (HTGR) technology development and its modeling and simulation needs. There is a particular need for advanced modeling and simulation tools to predict thermal-fluid behavior in the nuclear reactor primary system, especially in the core and the lower and upper plena, during safety-related transients. In this report, two main such activities are presented relevant to the HTGRs: (1) three-dimensional (3D) computational fluid dynamics (CFD) validation using benchmark data from the upper plenum of Texas A&M University’s 1/16th scaled very-high-temperature gas-cooled reactor (VHTR), and (2) development of wall heat-transfer correlation for laminar flow in a wall-heated pipe. The CFD tool validation exercises can be helpful to choose the models and CFD tools to simulate and design specific components of the HTRGs such as upper plenum where jet mixing is a complex phenomenon. In a loss of forced circulation event, the laminar flow can be observed during the development of natural circulation flow. This work includes the development and validation of heat transfer correlations for laminar flow using the Nek5000 CFD code due to limited available experimental data for laminar flow conditions to guide low-order models (1D). In this report, the flow characteristics of a single isothermal jet discharging into the upper plenum was investigated using the Nek5000 Large-Eddy Simulation (LES) CFD tool. Several numerical simulations were performed for various jet-discharged Reynolds numbers ranging from 3,413 to 12,819. A grid-independent study was performed. The numerical results of mean velocity, root-mean-square fluctuating velocity, and Reynolds stress were compared against the benchmark data. Good agreement was obtained between simulated and measured data for axial mean velocities, except near the upper plenum hemisphere. The maximum predicted errors for axial mean velocities at various normalized coolant channel diameter heights of 1, 5, and 10 are 1.56%, 1.88%, and 3.82%, respectively. In addition, the predicted root-mean-square fluctuating velocity and Reynolds stress are qualitatively in agreement with the experimental data. The Nek5000 code was used to develop wall-heat transfer correlation for laminar flow in a cylindrical tube. Several simulations were performed for various Reynolds flow and wall-heat fluxes. A new heat transfer correlation was developed using data from Nek5000 simulation results and regression functions in Matlab. The developed heat transfer correlation is valid for various Reynolds flows from 200 to 2000. The predicted R² value for model fit was 0.875, which ensures that 87.5% of the model data lies on the Nek5000 data. Moreover, a machine learning (ML) tool was used to train and test the Nek5000 data. A good fit of the ML-based model was observed with the test data.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
NASA NTRS2008

Gravity Scaling of a Power Reactor Water Shield

Water based reactor shielding is being considered as an affordable option for use on initial lunar surface power systems. Heat dissipation in the shield from nuclear sources must be rejected by an auxiliary thermal hydraulic cooling system. The mechanism for transferring heat through the shield is natural convection between the core surface and an array of thermosyphon radiator elements. Natural convection in a 100 kWt lunar surface reactor shield design has been previously evaluated at lower power levels (Pearson, 2007). The current baseline assumes that 5.5 kW are dissipated in the water shield, the preponderance on the core surface, but with some volumetric heating in the naturally circulating water as well. This power is rejected by a radiator located above the shield with a surface temperature of 370 K. A similarity analysis on a water-based reactor shield is presented examining the effect of gravity on free convection between a radiation shield inner vessel and a radiation shield outer vessel boundaries. Two approaches established similarity: 1) direct scaling of Rayleigh number equates gravity-surface heat flux products, 2) temperature difference between the wall and thermal boundary layer held constant on Earth and the Moon. Nussult number for natural convection (laminar and turbulent) is assumed of form Nu = CRa(sup n). These combined results estimate similarity conditions under Earth and Lunar gravities. The influence of reduced gravity on the performance of thermosyphon heat pipes is also examined.

Reid, Robert S.↗
DOE OSTI2020

RELAP5-3D Modeling of High Temperature Test Facility (HTTF) Test PG-26

The High Temperature Test Facility (HTTF) at Oregon State University (OSU) is a scaled integral effects experiment designed to investigate transient behavior in high-temperature gas-cooled nuclear reactors with prismatic fuel and reflector blocks. Several tests have been completed, and more are still planned to at the HTTF, including depressurized conduction cooldown (DCC) and pressurized conduction cooldown (PCC) transients. This report analyses test PG-26, a progression of the Double Ended Inlet-Outlet Crossover Duct Break transient that is referred to as a DCC. PG-26 has been performed at the HTTF between May 30 and June 30, 2019. Core initial conditions (i.e., before the DCC started) have been met using low power (<100 kW) and two of ten available electric heaters. The DCC transient was initiated during the 50 th hour of the test. The break valves were opened, and hot helium from the core and cold helium from the reactor cavity simulation tank (RCST) started mixing. The gases flowed in a countercurrent fashion, where the top half of the hot duct contained hot helium that flowed in one direction and cold helium that flowed in the other direction in the bottom half of the duct. After the pressure and density reached equilibrium, the event entered a diffusion mode. The onset of a reverse natural circulation was not observed during the DCC period of the test. Version 4.4.2ie of the RELAP5-3D computer code has been used to model the HTTF PG-26 test, and results have been compared to available high-quality measured data. The model used in this study is the quality-controlled HTTF RELAP5-3D model (HTTF base 2018-04-19 QA), originally developed by P. Bayless. The report includes RELAP5-3D results of the “base calculations” as well as some sensitivities to important uncertain model inputs, such as primary helium mass flow rate, core ceramic thermal properties, as well as heat evacuation and loop friction models. Using the base RELAP5-3D model predicts a countercurrent helium flow in the hot duct observed at the beginning of the DCC, but instead of going into a molecular diffusion mode, the model predicts the onset of natural convection. Increasing friction in the core and hot duct prevents the natural convection from happening in some of the simulations. Although some temperatures are well predicted (and even overpredicted), the general tendency is to underpredict the ceramic and helium temperatures and heat removal rates during the DCC, resulting in many of the assessment findings being in minimal or insufficient agreement with the data. It is worth noting that the described discrepancies between measured data and RELAP5-3D predictions are not RELAP5-3D code limitations. More so, they reflect limitations in boundary condition and thermal property knowledge. While the RELAP5-3D calculations of the test provide some insights into what happens during the transient, and point to missing or potentially uncertain data to which the experimenters can direct their attention, the principal conclusion is that the PG-26 test data are insufficient for a system code assessment.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
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 OSTI2023

Effect of inlet throttling on thermohydraulic instability in a large scale water-based RCCS: An experimental study

The objective of the present experimental study is to investigate the effect of inlet throttling on the thermohydraulic stability of a large scale water-based Reactor Cavity Cooling System (RCCS). The test was performed using the water-based Natural convection Shutdown heat removal Test Facility (NSTF) at Argonne, which represented a ½ axial scale and 12.5° sector slice of the full scale Framatome 625 MW t SC-HTGR RCCS concept. A two-phase steady state was first established through direct condensate refill, followed by increased inlet throttling over 10 stages, corresponding to a loss coefficient K over the range of 0.05–653. With the inlet throttling gradually increased, the system experienced a unique transition process between stabilization and destabilization. Through a stability analysis, three instability mechanisms were identified in the present test, including a compound mechanism due to both natural circulation oscillations (NCOs) and density wave oscillations (DWOs), Type-II DWOs, and geysering.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

System code evaluation of near-term accident tolerant claddings during pressurized water reactor station blackout accidents

Following the Fukushima Daiichi nuclear accident in 2011, researches on Accident-tolerant fuels (ATFs) are currently of high interest in not only the nuclear industry but also governmental and international organizations. In this work, a quantitative evaluation of the performance of monolithic FeCrAl cladding and Cr-coated Zircaloy cladding has been performed for Pressurized Water Reactor (PWR) Station Blackout (SBO) accidents. A generic PWR model has been built in system thermal-hydraulics code TRACE based on the Surry Nuclear Power Station with counter-current natural circulation modelling capability for hotleg and steam generator U-tube components during the accidents. The base model results are then compared to MELCOR and RELAP simulations to verify the system component implementation in TRACE. Two PWR SBO scenarios were investigated, including: short-term SBO and long-term SBO with early reactor coolant pump (RCP) seal failure. These scenarios are defined to be very similar to the accidents studied in the State-of-the-Art Reactor Consequence Analysis (SOARCA) project. TRACE code is modified to reflect the oxidation kinetics of FeCrAl and Cr-coating. Larson-Miller creep rupture model is also implemented in TARCE using its built-in control systems to simulate the creep rupture of hotlegs. Additionally, the comparison between the TRACE models with and without the counter-current flow modeling resulted in significant difference when comparing ATF cladding to Zircaloy for short term SBO, while it marginal impacted the performance during long term SBO with RCP seal failure. For short term SBO, both ATF cladding underwent hot leg creep rupture ~20 min after Zircaloy cladding. While Zircaloy and Cr-coated cladding had generated significant amount of hydrogen gas (>10 kg) before the creep rupture event, FeCrAl cladding had only generated <0.5 kg of hydrogen gas. For long term SBO with RCP seal failure, significant hydrogen generation and fuel melting was predicted before hot leg creep rupture for the ATF cladding while providing only 10–20 min additional coping time compared to Zircaloy.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Transformational challenge reactor analysis to inform preconceptual core design decisions: Sensitivity study of transient analysis in a hydride-moderated microreactor

The Transformational Challenge Reactor (TCR) program aims to demonstrate a revolutionary design approach enabled by advanced manufacturing and data analytics in the design of nuclear reactors. This article discusses scoping analyses of preconceptual designs to inform TCR design decisions and the evaluation of sensitivities and uncertainties on postulated transient scenarios. The applicability of the systems codes TRACE and RELAP5-3D to TCR transient analysis are examined, and RELAP5-3D models are used to examine the transient response of two candidate core designs at multiple power levels. Then, the uncertainty quantification code RAVEN is used to quantify the effect of several design parameters on reactivity-initiated accident (RIA) progression at hot zero power (HZP) and hot full power (HFP) as well as to assess the impact of uncertainties in transient parameters for the pressurized loss of forced cooling (PLOFC).When results were compared, TRACE and RELAP5-3D showed good agreement in their ability to predict system behavior, but RELAP5-3D calculations were closer to analytical predictions for the RIA. Furthermore, models for a PLOFC accident in two designs (a UO2 and TRISO core) at multiple power levels showed greater temperature margins for the TRISO core at all power levels. Using this information, along with other scoping analyses and constraints, the TCR design team selected a power level of 3 MWth and a TRISO-based core design. For this design RAVEN was applied to vary RIA parameters in RELAP5-3D models at HZP and HFP to understand the effect on figures of merit such as peak power, fuel and coolant temperature, and energy deposition. This sensitivity study found that the inserted reactivity worth was the most important parameter controlling all figures of merit, but for insertion up to 1.5$ no failure of TCR fuel is anticipated. For constant reactivity insertions, the magnitude of the fuel temperature coefficient was found to have the greatest effect on all figures of merit under most circumstances. These results not only demonstrate the anticipated robust safety of the proposed TCR fuel form but also provide a reference for future metal-hydride moderated systems to understand RIA behavior. In the PLOFC, the impact of heat transfer enhancement due to wavy flow channel effects dominated the variance in peak temperatures, and variations in heat exchanger elevation provided the greatest control on natural circulation flow rate. Additionally, no fuel particle failure is anticipated in the PLOFC.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Flow reversal benchmark of a one-sided heated narrow rectangular channel with CATHARE and RELAP5

Flow reversal in narrow coolant channels can be a crucial phenomenon for the safety of research reactors with a downward nominal flow direction. During a loss of forced flow accident, the downward flow stagnates briefly before transitioning into an upward natural circulation flow. The fuel may be damaged if dryout occurs and threshold fuel and/or cladding temperatures are exceeded. A comprehensive study is provided for flow reversal in narrow rectangular channels by examining experimental data and conducting software model analyses. The literature on flow reversal was reviewed, and selected experimental datasets were used to benchmark against CATHARE and RELAP5 models and also compare the code calculations with each other. The experimental data comes from flow reversal tests conducted with a narrow rectangular channel with one-sided heating. The results were compared with experimental data for successful flow reversal tests and predicted dryout power for dryout conditions. Also, the study examined the effects of the pump coastdown period, inlet liquid temperature, system pressure, and localized pressure drops. The experimental results showed that shorter coastdown periods, reduced pressure drops, and lower coolant inlet temperatures increased the dryout power. However, the system pressure did not noticeably affect the results. The simulation results showed that both CATHARE and RELAP5 agreed with experimental data, capturing the trends of the experimental results. Slight differences between each code calculation, as well as the predicted and measured dryout powers, were attributed to experimental uncertainties and the modeling of physical phenomena such as wall nucleation, interfacial heat transfer, drag coefficients, and critical heat flux. Overall, this study provides an understanding of flow reversal and the prediction capabilities of thermal-hydraulics software models. In conclusion, a future study of the flow reversal benchmark of a narrow rectangular channel with two-sided heating may provide additional valuable insights.

CATHARE↗
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 OSTI2024

Radionuclide tracing based in situ corrosion and mass transport monitoring of 316L stainless steel in a molten salt closed loop

In the study, we report an in situ corrosion and mass transport monitoring method developed using a radionuclide tracing technique for the corrosion study of 316L stainless steel (316L SS) in a NaCl–MgCl 2 eutectic molten salt natural circulation loop. This method involves cyclotron irradiation of a small tube section with 16 MeV protons, later welds at the hot leg of the molten salt flow loop, generating radionuclides 51 Cr, 52 Mn, and 56 Co at the salt–alloy interface. By measuring the activity variations of these radionuclides at different sections along the loop, both the in situ monitoring of the corrosion attack depth of 316L SS and corrosion product transport and its precipitation in flowing NaCl–MgCl 2 molten salt are achieved. While 316L SS is the focus of this study, the technique reported herein can be extended to other structural materials being used in a wide range of industrial applications.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Pronghorn: A Multidimensional Coarse Mesh Application for Advanced Reactor Thermal-Hydraulics

This paper presents an overview of Pronghorn, a multiscale thermal-hydraulic (T/H) application developed by Idaho National Laboratory and the University of California, Berkeley. Pronghorn, built on the open-source finite element Multiphysics Object-Oriented Simulation Environment (MOOSE), leverages state-of-the-art physical models, numerical methods, and nonlinear solvers to deliver fast-running advanced reactor T/H simulation capabilities within a modern software engineering environment. This work summarizes the physical models, multiphysics and multiscale coupling, and numerical discretization in Pronghorn with emphasis on our initial target application to pebble bed reactors (PBRs). A diverse set of applications are shown to depressurized natural circulation in the SANA experiments, forced convection in the Pebble Bed Modular Reactor, three-dimensional (3-D)/one-dimensional coupling of Pronghorn and RELAP-7 systems T/H for loop analysis in the High Temperature Reactor Power Module, and forced convection in the Mark-1 Pebble Bed Fluoride-Salt-Cooled High-Temperature Reactor. A multiphysics coupling of Pronghorn, RELAP-7, and Griffin deterministic neutronics for a gas-cooled PBR demonstrates the capability of the MOOSE framework for reactor design calculations. These applications highlight the verification and validation underlying Pronghorn’s software development while emphasizing features that improve upon capabilities offered by legacy tools in areas such as 3-D unstructured meshing, physics modeling, and multiphysics coupling.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Progress in Fast Modular Reactor Conceptual Design

The Fast Modular Reactor (FMR) is a 100-MW(thermal) gas-cooled fast reactor being developed by General Atomics Electromagnetic System with the goal of developing a FMR for flexible and dispatchable power to the U.S. electricity market in the mid-2030s. The conceptual design aims to develop and verify simplified design features. These include an inert helium gas coolant, pellet-loaded fuel rods, installations with air cooling as ultimate heat sink, and small and passive heat removal systems. The goal is to ensure the development of a safe, maintainable, cost-effective, and distributed nuclear energy-generating station. The baseline technologies selected to achieve this goal are a helium coolant that is an inert gas with no chemical reaction with structural components, not activated, single phase, enabling high-temperature operation and a high thermal efficiency Brayton cycle; conventional uranium dioxide (UO 2 ) fuel, which is the most widely used and well-known fuel material, capable of high burnup (100 MWd/kg) and a long fuel life; and silicon carbide composite (SiGA®) cladding and internal structures that are chemically inert in the helium environment, exceptionally radiation tolerant, and being derisked by accident tolerant fuel technology development. Further, the reactor was specifically designed with passive safety features, including high-temperature in-core materials and a reactor vessel cooling system consisting of cooling panels of naturally circulating water. The passive safety of the core was confirmed for the depressurized loss-of–forced cooling accident, which showed the peak cladding temperature at ~1600°C during the transient, which is below the current design limit of 1800°C. The conceptual design of the FMR has been conducted for the reactor system, vessel system, generator and turbomachine, instrumentation and control, residual heat removal system, plant service system, and containment, as well as pre-application licensing documents.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

Validation of a Hybrid Domain Overlapping Coupling Between SAM and CFD Against the TALL-3D Transients

The System Thermal Hydraulics (STH) code SAM has been coupled to the Computational Fluid Dynamics (CFD) code Simcenter STAR-CCM+ utilizing a hybrid domain overlapping method with an explicit coupling in time. The coupling aims to extend the STH code’s applicability to scenarios where local momentum and energy transfers are important yet difficult for STH codes to capture, such as three-dimensional (3D) mixing. The coupling method’s numerical stability was verified in the past against two closed-loop configurations, and it was validated against a double T-junction experiment with 3D scalar mixing. In the present work, the coupling method is validated against the TALL-3D STH/CFD coupling benchmark facility. TALL-3D is a three-legged, liquid-metal facility with a large, pool-type enclosure (test section) that exhibits 3D flow effects to be modeled by a CFD code. The rest of the system exhibits approximately 1D behavior well-predicted by an STH code. First, the present STAR-CCM+ CFD model of the 3D test section is validated against experimental data. Then, the SAM-STARCCM+ coupled model is validated against six different TALL-3D steady states, including SAM standalone model results for comparison. Lastly, the SAM-STARCCM+ coupled model is validated against two TALL-3D transients, one exhibiting flow reversal in the test section and one exhibiting nonlinear, Limit Cycle Oscillations (LCO). For the first transient, the SAM-STARCCM+ coupled model properly predicts an increase in the test section’s inlet temperature during flow reversal, and this is not predicted by the SAM standalone model. Following flow reversal, the SAM-STARCCM+ coupled model better-predicts the initial flow recovery and following oscillations as the system approaches a final natural circulation state. For the second transient, no true final steady state is observed due to LCO. Neither the SAM-STARCCM+ coupled model nor the SAM standalone model can perfectly capture the experiment’s changing oscillation frequency during the transient. However, the SAM-STARCCM+ coupled model does reproduce the oscillatory feedback observed in the system. This is a significant achievement as the SAM-STARCCM+ coupled model only uses an explicit coupling in time, as opposed to a semi-implicit coupling. In comparison, previous STH/CFD coupling efforts of the TALL-3D facility required semi-implicit coupling to obtain similar results.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Transient Modeling and Simulation of a Generic Stable Salt Reactor

A SAM system-level model of a generic stable salt reactor has been developed to investigate thermal-hydraulic behavior and safety performance under steady and transient conditions. The model integrates information generated from a reactor physics analysis using PROTEUS and PERSENT, and a computation fluid dynamics (CFD) analysis using STAR-CCM+. A loose, iterative coupling scheme between PROTEUS and SAM is implemented to calculate the equilibrium power and temperature distributions in the steady-state critical core condition. The converged steady-state model is then used in PERSENT to calculate the four reactivity feedback temperature coefficients (Doppler, fuel density, coolant density, and core radial expansion) and kinetic parameters that are needed in SAM to model the temperature feedback effects in transient simulations. Within the fully enclosed liquid fuel pins, natural convection is the dominant heat transfer mechanism. The STAR-CCM+ model of the fuel pin considers conjugate heat transfer from the liquid fuel salt to the pin cladding and external reactor coolant. The CFD results of the axial and radial temperature profiles are used to empirically determine an effective fuel salt thermal conductivity in the SAM fuel pin model so that the temperatures predicted by the SAM model match as closely as possible the CFD results. In the central region of the fuel pin, the effective thermal conductivity is as high as similar to 60 times the physical fuel salt thermal conductivity. The whole-plant SAM model is then used to simulate an unprotected station blackout transient. The results of this simulation showed that the large negative fuel axial expansion reactivity feedback reduces fission power to similar to 2.4% nominal power. The core is cooled by natural circulation, which removes heat in the core to the emergency heat removal system, and ultimately, to the ambient. However, peak fuel salt and cladding temperatures can potentially reach as high as 1500 K, albeit briefly, if the shutdown mechanism fails to operate.

stable salt reactor; transient simulations; system↗
DOE OSTI2026

A New 1D Model for Thermal Mixing and Stratification in Advanced Reactor Transients

Thermal mixing and stratification in large pools and enclosures play a critical role in the safety and performance of pool-type nuclear reactors, particularly during transient scenarios involving significant temperature differences between incoming and bulk coolant. Accurate modeling of these phenomena is essential for predicting system behavior and supporting passive safety features such as natural circulation. Here, this paper presents a new 1D model for thermal mixing and stratification, developed and implemented in the SAM code. The model represents a large pool as 1D coolant jet channels and zero-dimensional bulk pool volumes, enabling the simulation of a wide range of flow configurations, including hot and cold jet interactions, stratified layers, and the influence of complex geometries such as ceilings, free surfaces, and internal obstacles. Heat exchange between jet and pool regions is governed by closure relations calibrated against 3D computational fluid dynamics (CFD) simulations. The model improves upon earlier approaches by incorporating time-dependent jet characteristics and capturing the associated delay effects more accurately. Code-to-code comparisons and validation against experimental data from the Thermal Stratification Test Facility demonstrate the model’s accuracy and flexibility. This work offers two key contributions: (1) an efficient and robust method for simulating thermal mixing and stratification at the system level, eliminating the need for external coupling between system analysis codes and CFD, and (2) a significant enhancement of SAM’s capabilities to analyze thermal stratification phenomena in advanced reactor systems.

SAM↗
DOE OSTI2026

Large Eddy Simulation of Low-Reynolds-Number Turbulent Flow of Low-Prandtl-Number Fluid in a Tight Lattice Bundle for Assessment of Reynolds-Averaged Navier-Stokes Turbulence Model

The MARVEL (Microreactor Applications Research Validation and Evaluation) microreactor utilizes natural circulation as core cooling mechanism and liquid metal as a primary coolant. Moreover, the reactor core has a pitch-to-diameter ratio of 1.056, which is considered a tight lattice configuration. Numerous studies have widely reported that Reynolds-Averaged Navier-Stokes (RANS) turbulence models inaccurately predict heat transfer in liquid metals and fail to capture flow pulsations that can occur within tight lattices, leading to further inaccuracies in simulation results. Therefore, evaluating the accuracy of RANS turbulence models in the thermal-hydraulic analysis of the MARVEL microreactor core is crucial for assessing reactor safety. In this study, a Large Eddy Simulation (LES) of the MARVEL microreactor core subchannel was conducted and compared with a RANS simulation to evaluate the accuracies and conservatism of the RANS model. The flow pulsation in a tight lattice predicted by LES enhanced the heat transfer, whereas the RANS model underpredicted it. Consequently, the RANS model predicted the peak cladding temperature higher than the LES model, but the discrepancy between the two approaches was not significant due to the good thermal characteristics of the liquid metal. It can be concluded that the steady-state RANS model is effective for the thermal analysis of liquid-metal-cooled MARVEL microreactor core and can provide conservative predictions from a safety analysis perspective.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗