Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Validation of Pronghorn for Natural-Circulation Molten Salt Loops

This paper presents the development and validation of a high-fidelity thermal-hydraulic model of a molten salt natural circulation flow loop, designed for integration within a digital twin framework. The study evaluates the performance of Idaho National Laboratory’s Pronghorn against experimental data from Texas A&M University Molten Salt Flow Loop (MSFL) four Hitec-salt test benchmark data. Natural circulation of high-Prandtl-number fluids exhibits complex, counter-intuitive flow patterns that make pointwise thermocouple readings unreliable. Experimental work at TAMU’s MSFL provides benchmark data, including flow visualization at a test-section and centerline steady-state temperature measurements along the loop. Validation includes four single-phase natural circulation test cases with Hitec salt. Key metrics include flow profile agreement and steady-state temperature accuracy. Pronghorn results for two-dimensional single-phase agree qualitatively with the experimental flow profile. This paper illustrates the importance of Computational Fluid Dynamics (CFD) in elucidating the behavior of high-Prandtl-number thermal-hydraulics, along with how misleading centerline temperature measurements can be. Pronghorn reproduces the axial and radial stratification that makes single thermocouple readings unreliable. Future research will focus on reduced-order modeling techniques to enable rapid simulation suitable for real-time digital twin applications. The validated cases provide a basis for developing reduced-order surrogates aimed at real-time digital-twin applications.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Intracore Natural Circulation Study in the High Temperature Test Facility

The development of the Modular High-Temperature Gas-Cooled Reactor is a significant milestone in advanced nuclear reactor technology. One of the concerns for the reactor’s safe operation is the effects of a loss-of-flow accident (LOFA) where the coolant circulators are tripped, and forced coolant flow through the core is lost. Depending on the steam generator placement, loop or intracore natural circulation develops to help transfer heat from the core to the reactor cavity, cooling system. This paper investigates the fundamental physical phenomena associated with intracore coolant natural circulation flow in a one-sixth Computational Fluid Dynamics (CFD) model of the Oregon State University High Temperature Test Facility (OSU HTTF) following a loss-of-flow accident transient. This study employs conjugate heat transfer and steady-state flow along with an SST k-ω turbulence model to characterize the phenomenon of core channel-to-channel natural convection. Previous studies have revealed the importance of complex flow distribution in the inlet and outlet plenums with the potential to generate hot coolant jets. For this reason, complete upper and lower plenum volumes are included in the analyzed computational domain. CFD results also include parametric studies performed for a mesh sensitivity analysis, generated using the STAR-CCM+ software. The resulting channel axial velocities and flow directions support the test facility scaling analysis and similarity group distortions calculation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1986

Laser anemometry measurements of natural circulation flow in a scale model PWR reactor system

The natural circulation of a single phase fluid in a scale model of a pressurized water reactor system during a postulated grade core accident is analyzed. The fluids utilized were water and SF6. The design of the reactor model and the similitude requirements are described. Four LDA tests were conducted: water with 28 kW of heat in the simulated core, with and without the participation of simulated steam generators; water with 28 kW of heat in the simulated core, with the participation of simulated steam generators and with cold upflow of 12 lbm/min from the lower plenum; and SF6 with 0.9 kW of heat in the simulated core and without the participation of the simulated steam generators. For the water tests, the velocity of the water in the center of the core increases with vertical height and continues to increase in the upper plenum. For SF6, it is observed that the velocities are an order of magnitude higher than those of water; however, the velocity patterns are similar.

Kadambi, J. R.↗
NASA NTRS1987

Measurements of natural circulation flow in a scale model PWR reactor system during postulated degraded core accidents using laser anemometry

The natural circulation of a single-phase fluid in a scale-model pressurized water reactor system was studied during a postulated degraded core accident. A half section of a one-seventh scale model with a plexiglass adiabatic window was employed. Water and SF6 were used as the fluid. LDA was used to perform velocity measurements along the center plane of the model at five elevations. It was found that the recirculation flow patterns are nearly symmetric except near the hot legs and in the upper head and that the fluid in the upper plenum is well mixed.

Kadambi, J. R.↗
DOE OSTI2025

Corrosion of 316 L stainless steel under the natural circulation of molten NaCl-MgCl 2 salt

The corrosion behavior of 316 L stainless steel (SS) was studied via the natural circulation of molten eutectic NaCl-MgCl 2 salt through a microloop. The post-corrosion tested 316 L SS microloop sections were characterized with microscopy techniques to determine the microstructural and microchemical changes that occurred at the alloy/salt interface. It was found that 316 L SS showed heterogeneous dissolution at the hot-leg, whereas deposition of corrosion products occurred at the cold-leg. For the first time, experimentally obtained molten salt flow-induced corrosion of 316 L SS results were combined with computational thermodynamic-kinetic models to validate the dissolution and deposition in terms of elemental compositional changes at the alloy/salt interface. The thermodynamic-kinetic modeling predicted that the heterogeneous dissolution of Cr and Fe from the hot-leg section of 316 L SS persisted throughout the salt circulation. The model also estimated that, despite Cr deposition starting earlier than Fe, the total redeposition of Fe is expected to be significantly greater than that of Cr over the circulation of salt. Furthermore, the modeling accurately predicted the subsurface enrichment of Mo which is attributed to the reduced Cr activity and the relatively higher diffusion rate of Mo within the alloy matrix. Here, the agreement between modeling and experimental results confirms that Fe chlorides dissolve at the hot-leg and subsequently deposit at the cold-leg due to activity changes driven by the thermal gradient. By contrast, Cr was not detected in the cold-leg deposits, which is attributed to its weaker temperature dependence on activity, limiting its redeposition under these conditions.

36 - MATERIALS SCIENCE↗
NASA NTRS1992

Natural circulation decay heat removal from an SP-100, 550 kWe power system for a lunar outpost

This research investigated the decay heat removal from the SP-100 reactor core of a 550-kWe power system for a lunar outpost by natural circulation of lithium coolant. A transient model that simulates the decay heat removal loop (DHRL) of the power system was developed and used to assess the system's decay heat removal capability. The effects of the surface area of the decay heat rejection radiator, the dimensions of the decay heat exchanger (DHE) flow duct, the elevation of the DHE, and the diameter of the rise and down pipes in the DHRL on the decay heat removal capability were examined. Also, to determine the applicability of test results at earth gravity to actual system performance on the lunar surface, the effect of the gravity constant (1 g and 1/6 g) on the thermal behavior of the system after shutdown was investigated.

El-Genk, Mohamed S.↗
DOE OSTI2025

CFD modeling of natural circulation in LiCl-KCl molten salt closed loop

Characterizing flow within a molten salt closed-loop system is crucial for assessing system requirements, evaluating performance, and identifying potential flaws. Direct flow measurement using instrumentation is challenging due to extreme environmental conditions and the limitations associated with measuring molten salt flow under natural convection. Here, this study aims to provide comprehensive insights into the thermal-hydraulic behavior of a closed loop, with a particular focus on temperature distribution and velocity prediction. The Computational Fluid Dynamics (CFD) model demonstrated the capability to effectively simulate and predict both temperature distributions and flow velocities within the molten salt loop. The CFD model's predictive capability was validated by its ability to replicate temperature measurements under varying boundary conditions. The analysis revealed that the CFD model tends to underpredict temperatures in the cold leg and overpredict them in the hot leg, highlighting the need for continuous model refinement and acknowledging the limitations of using a steady-state approach. Furthermore, the potential of using external temperature measurements to estimate internal molten salt temperatures and predict flow velocity was explored, revealing that this approach could introduce up to a 5.5% error in flow velocity calculations. Line probes mapping temperature distributions across the tube's cross-section and molten salt provided valuable insights into temperature gradients, emphasizing the need for a thermal conductivity equation for molten salt with lower uncertainty to achieve more accurate temperature predictions of the system.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Fundamental flow characterization for LiCl–KCl molten salt natural circulation closed loop

Molten salt flow characterization presents numerous challenges due to the hostile environments that cause corrosion and erosion of instrumentation and structural components. To support future molten salt flow characterization under natural convection conditions, a LiCl–KCl eutectic molten salt loop was designed and constructed. The primary objective of this loop is to facilitate natural convection as the principal cooling mechanism, enabling insights into the physics behind natural convection in closed loops and serving as a baseline for loop dimensional and non-dimensional analyses. Here, this study summarizes the fundamental functions employed for flow characterization and introduces a model to estimate flow rate when intrusive instrumentation is unavailable. The approach relies solely on structural surface measurements to avoid flow disturbances. Furthermore, a comprehensive thermal–hydraulic analysis is included to deepen understanding of the system’s behavior.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Fuel Behavior Implications of Reactor Design Choices in Pressurized Water SMRs

Small pressurized water reactors (PWRs) can feature boron free operation, natural circulation mode, reduced height assemblies and/or long refueling cycles. This paper attempts to explore core design optimization for each of these design evolutions. In consequence, five core design layouts are developed incorporating boron free operation with continuous control rods insertion, natural circulation with low burnup/low power density design, natural circulation with high burnup/low power density design, forced circulation with standard core power density design, and forced circulation with high power density design. These cores’ performance is compared to a standard 4-loop PWR. The design process aims to improve the fuel cycle cost under safety constraints through core design optimization using CASMO4E/SIMULATE3 reactor physics codes and FRAPCON4.1 fuel performance assessment tool. Core modeling assumes standard 17x17 PWR fuel assemblies loaded with low enriched uranium (LEU) up to 5wt% or LEU+ (i.e., below 10wt% enrichment) pellets with gadolinium oxide (Gd2O3) as the burnable poison. Satisfactory core and fuel performances are obtained for all the designed cores under steady state and considered overpower transients. For low power density operation, long cycle lengths are achieved reaching a 2.5- and a 5-year cycles and peak rod-average burnup is pushed to 83 MWd/kgU. Other cycle lengths are maintained at 18 months. Boron free operation exhibits the ability to achieve longer cycle lengths at the cost of higher peaking factors leading to high local power and fuel temperatures which prevents sizable power uprates and is deemed uneconomical. Fuel assembly height reduction allows coolant velocity retrofit which enables higher core power density without violating structural integrity of the fuel assembly. As a result, a core power density of 123 kW/l is reached where total cladding hoop strain becomes the limiting parameter.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2026

Thermal-Hydraulic-Mechanical Modeling and Simulation of Sodium-Potassium–Cooled MARVEL Microreactor Core

The U.S. Department of Energy's Microreactor Program, with Idaho National Laboratory's development of a nuclear microreactor applications test bed named MARVEL, aims to support R&D for the deployment of small, transportable reactors across civilian, industrial, and defense sectors. The MARVEL microreactor, an 85-kWth thermal fission reactor, incorporates TRIGA nuclear fuel and a sodium-potassium eutectic as its primary coolant, designed for safety and efficiency, with natural circulation eliminating the risk of critical heat flux conditions. The reliance on natural circulation for primary cooling means the reactor avoids using fuel spacers to minimize core pressure drop, which could disrupt the primary coolant's natural flow. However, the reactor core’s tight P/D ratio of 1.05, in the absence of fuel spacers, could pose a risk of fuel rod contact and increased peak cladding temperatures. To ensure the reactor safety, this study conducted computational modeling and simulations to investigate the reactor's thermal-hydraulic-mechanical characteristics, including the reactor core heat transfer coefficients, the potential for rod-to-rod contact, and assessed its impact on peak cladding temperature and overall reactor safety. The computational analyses of the MARVEL microreactor core revealed that the thermal deformation of fuel rods under worst-case scenario may lead to the fuel rod contact, but the peak cladding temperatures will remain significantly lower than the safety criteria, ensuring the safety operation reactor without fuel spacers under normal operating conditions.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗
DOE OSTI2025

Device for steam cladding oxidation testing at TREAT

To compare the chemical degradation of conventional zirconium alloy (Zry) cladding to advance silicon carbide (SiC) cladding in a post loss of coolant accident (LOCA) environment, new nuclear testing capabilities are necessary. The Transient Reactor Test (TREAT) Facility at Idaho National Laboratory (INL) has matured its transient fuel testing capabilities since its 2017 restart. The most recent experiment architecture is the Transient Water Irradiation System in TREAT (TWIST), which is designed to support qualification of accident tolerant fuels in light water reactors. INL has designed and analyzed a natural circulation steam flow modification for TWIST to produce prototypic conditions of cladding oxidation. The in-situ device will be electrically heated to drive natural circulation. Moreover, the SiC cladding requires heating above 1700 °C to observe failure, thus internal prototypic nuclear heating with radiation effects will be used. Thermal hydraulic analysis with RELAP5-3D (Reactor Excursion and Leak Analysis Program) estimated steam fluxes greater than 50 mg cm −2 s −1 can be achieved. These fluxes are adequate to test Zry cladding according to draft regulatory guides and to test SiC cladding according to past experiments.

Oxidation↗
DOE OSTI2026

Development of SAM Code Capabilities for Safety Analysis of GCR Air-ingress Events

Air-ingress following a depressurized loss-of-forced-cooling (DLOFC) event is a challenging, multiphysics safety scenario for High-Temperature Gas-Cooled Reactors (HTGRs), involving coupled gas composition transport, buoyancy-driven flow redistribution, graphite oxidation, and structural heat-up. Despite its importance — air ingress is a key scenario identified in the PIRT process for the HTGRs — existing system-level safety codes have lacked the integrated capability to simulate the complete event sequence with high confidence. This report documents the development, validation, and demonstration of three new capabilities in the SAM code to address this gap: (1) a multi-component gas mixture flow model with binary diffusion to track the helium-air composition and its effect on system density and flow; (2) a 0-D graphite oxidation model based on the Roes correlation, including oxygen consumption and exothermic heat release; and (3) an isentropic critical flow model for accurate representation of primary system depressurization through a break. These capabilities are validated against two benchmark experiments. The NSTF heavy-gas ingress experiment validates the multi-component flow model: SAM correctly reproduces the rapid buoyancydriven flow stagnation and subsequent natural circulation recovery driven by composition-dependent density changes. The NACOK graphite oxidation experiment validates the oxidation model: SAM predicts a bottom-level graphite weight loss of 25%, in close agreement with the measured 24%, and reproduces the strong axial nonuniformity and block-geometry dependence of oxidation, at a level comparable to the SPECTRA and TINTE codes. The validated capabilities are then exercised together in an integrated, reactor-scale simulation of a DLOFC air-ingress transient in a simplified HTR-PM pebble-bed reactor. In a single calculation spanning approximately 8 days, SAM reproduces the complete accident sequence: rapid depressurization, densityand diffusion-driven air ingress over ˜15 hours, onset of buoyancy-driven natural circulation, exothermic graphite oxidation with a peak fuel temperature at ˜62 hours, and eventual passive cooldown. These results demonstrate that SAM now provides the nuclear community with a preliminarily validated, modern systemlevel tool for HTGR air-ingress safety analysis, filling a recognized capability gap. Future extensions to broaden species tracking, improve oxidation chemistry, and refine the reactor model are discussed.

Yang, Gang↗
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.↗