Search NASA⌕ Search

SEARCH · Search NASA

Results for “Critical Heat Flux”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS2026

Decoupling Surface Topography from Gravitational Acceleration in Cryogenic Pool Boiling

Boiling heat transfer is governed by a complex interplay between surface conditions and gravitational acceleration. To isolate the sole effects of gravity, we investigated the pool boiling characteristics of liquid nitrogen on atomically smooth silicon dioxide (SiO 2 ) surfaces under terrestrial (1-g) and reduced gravity (0±0.02 g) conditions achieved via parabolic flight. Our results quantify a drastic reduction in the critical heat flux (CHF) in reduced gravity, decreasing from 16.15 W/cm 2 at 1-g to 5−6 W/cm 2 at μ-g due to the suppression of buoyancy. Conversely, we observed a distinct increase in the heat transfer coefficient (HTC) in the reduced gravity condition prior to CHF. By utilizing a surface with a maximum peak-to-valley height of ≈36.7 nm and low contact angle hysteresis (<10°), we confirm this HTC enhancement is an intrinsic response to the gravitational environment, decoupled from surface-defect-induced nucleation. These findings demonstrate that the influence of surface topography is significantly more prominent in reduced gravity than in terrestrial conditions, providing a critical baseline for rationalizing the design of cryogenic thermal management systems in space and quantum applications.

Liquid nitrogen↗
NASA NTRS2017

Gravity Effects in Microgap Flow Boiling

Increasing integration density of electronic components has exacerbated the thermal management challenges facing electronic system developers. The high power, heat flux, and volumetric heat generation of emerging devices are driving the transition from remote cooling, which relies on conduction and spreading, to embedded cooling, which facilitates direct contact between the heat-generating device and coolant flow. Microgap coolers employ the forced flow of dielectric fluids undergoing phase change in a heated channel between devices. While two phase microcoolers are used routinely in ground-based systems, the lack of acceptable models and correlations for microgravity operation has limited their use for spacecraft thermal management. Previous research has revealed that gravitational acceleration plays a diminishing role as the channel diameter shrinks, but there is considerable variation among the proposed gravity-insensitive channel dimensions and minimal research on rectangular ducts. Reliable criteria for achieving gravity-insensitive flow boiling performance would enable spaceflight systems to exploit this powerful thermal management technique and reduce development time and costs through reliance on ground-based testing. In the present effort, the authors have studied the effect of evaporator orientation on flow boiling performance of HFE7100 in a 218 m tall by 13.0 mm wide microgap cooler. Similar heat transfer coefficients and critical heat flux were achieved across five evaporator orientations, indicating that the effect of gravity was negligible.

microscale↗
NASA NTRS2017

Gravity Effects in Microgap Flow Boiling

Increasing integration density of electronic components has exacerbated the thermal management challenges facing electronic system developers. The high power, heat flux, and volumetric heat generation of emerging devices are driving the transition from remote cooling, which relies on conduction and spreading, to embedded cooling, which facilitates direct contact between the heat-generating device and coolant flow. Microgap coolers employ the forced flow of dielectric fluids undergoing phase change in a heated channel between devices. While two phase microcoolers are used routinely in ground-based systems, the lack of acceptable models and correlations for microgravity operation has limited their use for spacecraft thermal management. Previous research has revealed that gravitational acceleration plays a diminishing role as the channel diameter shrinks, but there is considerable variation among the proposed gravity-insensitive channel dimensions and minimal research on rectangular ducts. Reliable criteria for achieving gravity-insensitive flow boiling performance would enable spaceflight systems to exploit this powerful thermal management technique and reduce development time and costs through reliance on ground-based testing. In the present effort, the authors have studied the effect of evaporator orientation on flow boiling performance of HFE7100 in a 218 m tall by 13.0 mm wide microgap cooler. Similar heat transfer coefficients and critical heat flux were achieved across five evaporator orientations, indicating that the effect of gravity was negligible.

microscale↗
NASA NTRS2002

Length Scale and Gravity Effects on Microgravity Boiling Heat Transfer

Boiling is a complex phenomenon where hydrodynamics, heat transfer, mass transfer, and interfacial phenomena are tightly interwoven. An understanding of boiling and critical heat flux in microgravity environments is of importance to space based hardware and processes such as heat exchange, cryogenic fuel storage and transportation, electronic cooling, and material processing due to the large amounts of heat that can be removed with relatively little increase in temperature. Although research in this area has been performed in the past four decades, the mechanisms by which heat is removed from surfaces in microgravity are still unclear. In earth gravity, buoyancy is an important parameter that affects boiling heat transfer through the rate at which bubbles are removed from the surface. A simple model describing the bubble departure size based on a quasistatic force balance between buoyancy and surface tension is given by the Fritz [I] relation: Bo(exp 1/2) = 0.0208 theta where Bo is the ratio between buoyancy and surface tension forces. For small, rapidly growing bubbles, inertia associated with the induced liquid motion can also cause bubble departure. In microgravity, the magnitude of effects related to natural convection and buoyancy are small and physical mechanisms normally masked by natural convection in earth gravity such as Marangoni convection can substantially influence the boiling and vapor bubble dynamics. CHF (critical heat transfer) is also substantially affected by microgravity. In 1 g environments, Bo has been used as a correlating parameter for CHF. Zuber's CHF model for an infinite horizontal surface assumes that vapor columns formed by the merger of bubbles become unstable due to a Helmholtz instability blocking the supply of liquid to the surface. The jets are spaced lambda(sub D) apart, where lambda(sub D) = 2pi square root of 3[(sigma)/(g(rho(sub l) - rho(sub v)](exp 1/2) = 2pi square root of 3 L Bo(exp -1/2) = square root of 3 lambda(sub c) and is the wavelength that amplifies most rapidly. The critical wavelength, lambda(sub c), is the wavelength below which a vapor layer underneath a liquid layer is stable. For heaters with Bo smaller than about 3 (heaters smaller than lambda(sub D)), the above model is not applicable, and surface tension effects dominate. Bubble coalescence is thought to be the mechanism for CHF under these conditions. Small Bo can result by decreasing the size of a heater in earth gravity, or by operating a large heater in a lower gravity environment. In the microgravity of space, even large heaters can have low Bo, and models based on Helmholtz instability should not be applicable. The macrolayer model of Haramura and Katto is dimensionally equivalent to Zuber's model and has the same dependence on gravity, so it should not be applicable as well. The goal of this work is to determine how boiling heat transfer mechanisms in a low-g environment are altered from those at higher gravity levels. Boiling data using a microheater array was obtained under gravity environments ranging from 1.8 g to 0.02 g with heater sizes ranging from 2.7 mm to 1 mm. The boiling behavior for 2.7 mm at 0.02 g looked quite similar to boiling on the 1 mm heater at 1 g-the formation of a large primary bubble surrounded by smaller satellite bubbles was observed under both conditions. The similarity suggests that for heaters smaller than some fraction of I(sub c), coalescence and surface tension dominate boiling heat transfer. It also suggests that microgravity boiling can be studied by studying boiling on very small heaters.

Kim, Jungho↗
DOE OSTI2025

Assessment of CTF for Steady-state and Transient Post-CHF Conditions in Support of Time-at-Temperature Modeling Applications

The US nuclear industry is exploring options to improve operational economics and uprate the current fleet of light-water reactors by investigating transitioning to cladding performance–based safety criteria as opposed to the current limit, which requires complete avoidance of critical heat flux (CHF)/dryout. Past experience has shown that not all events leading to a dryout are severe enough to cause fuel performance degradation. Allowing temporary dryout of the fuel—that is, using a time-at-temperature (TaT) strategy—could allow for economic improvements via large power uprates and enhanced operational flexibility for current plants without compromising fuel integrity. To support this effort, the US Department of Energy is executing a comprehensive program that includes generating cladding material data under TaT conditions, developing new mechanistic models, and demonstrating modeling and simulation capabilities for transients of interest. This paper presents work performed to assess the CTF thermal-hydraulics subchannel code. CTF is a package used in the VERA core simulator, which will ultimately be used for TaT analysis. CTF will provide the thermal-hydraulic boundary conditions that will be needed for fuel performance analysis in the BISON code. Quantifying both the accuracy and uncertainty of post-CHF models will therefore be necessary. This paper outlines the strategy for the assessment of TaT and presents the results of using the steady-state and transient dryout experiments of the Boiling Fine-mesh Bundle Tests for CTF validation. The results show that the current model tends to overpredict steady-state critical power. This behavior translates to the transient tests, in which CTF is unable to capture transient dryout behavior. Some discussion of sensitivity analysis work being performed is provided to indicate which models must be further analyzed to properly model transient dryout and its uncertainty.

Salko Jr, Robert [ORNL] (ORCID:0000000253566679)↗
DOE OSTI2020

Transition Core Planning and Safety Analyses in Support of LEU Fuel Conversion of the University of Missouri Research Reactor (MURR)

The University of Missouri Research Reactor (MURR®) is one of six U.S. High Performance Research Reactors (USHPRR), including one critical facility, that is working with the National Nuclear Security Administration (NNSA) Office of Material Management and Minimization (M3) Reactor Conversion Program to convert from highly enriched uranium (HEU) to low-enriched uranium (LEU) fuel. The M3 Reactor Conversion USHPRR Project objectives include the development of LEU fuel element designs that will ensure safe reactor operations and to maintain the existing experimental performance of each facility. The work is being conducted through many inter-related activities being completed by four Project Pillars: Fuel Qualification (FQ), Fuel Fabrication (FF), Reactor Conversion (RC), and Cross Cutting (CC). A new type of LEU fuel based on an alloy of uranium-10 wt% molybdenum (U-10Mo) is expected to allow the conversion of those USHPRR, like MURR, requiring higher density fuels. The very-high-density LEU U-10Mo monolithic fuel is currently undergoing irradiation testing and post-irradiation examination under a planned and documented fuel qualification effort. The FQ Pillar will document fuel property and fuel performance data and qualify the fuel for use in these reactors. The FF Pillar is fabricating fuel for ongoing and future irradiation tests, as well as conducting fabrication demonstrations to validate or update preliminary fabrication assumptions. The FF Pillar is also working to develop and install commercial manufacturing capacity with the U-10Mo monolithic fuel to produce prototypic fuel. Working with the RC Pillar at Argonne, MURR has progressed through a preliminary fuel element design using preliminary data for the proposed monolithic alloy of U-10Mo. Analyses were completed in previous work that found for typical equilibrium operations with the preliminary LEU fuel element design, in conjunction with a power uprate to 12 MW and appropriate changes to the MURR Limiting safety system settings (LSSS), MURR will have adequate margins to safety for steady-state operations and postulated transient accidents and will have experimental performance in key locations that meets or exceeds current operations with HEU fuel. The purpose of this work is to develop a sequence of transition cycles that will enable MURR to transition from operation with the reactor core loaded with fresh LEU fuel elements only to typical equilibrium operations with mixed-burnup cores following conversion while meeting operational requirements on safety and experimental performance. It is expected that the use of fresh LEU fuel at conversion and subsequent low burnup of the LEU fuel elements that will initially be available for use following conversion will result in critical control blade positions that will substantially change the axial power distribution in the core and the neutron flux available in key experimental locations relative to equilibrium LEU operations. Given the constraints of MURR safety margins, operational practices, and production and research, a novel method has been developed to identify a transition sequence that minimizes the time MURR operates atypically compared to the current prototypic cycles using HEU fuel. The proposed transition sequence moves quickly to the same sort of equilibrium cycles for the LEU fuel that have already been evaluated in documented preliminary safety analyses. Although shifting the neutron flux peak to the lower half of the core during initial cycles with LEU at 12 MW reduces the experiment performance in some key locations relative to current HEU operations at 10 MW, all LEU cores provide an average performance that meets or exceeds that of HEU. An LEU cycle is reached that meets or exceeds the level of experimental performance predicted for current HEU and equilibrium LEU operations in more than 450 key locations identified by a reactor specialist at MURR by the 23rd cycle following conversion and that afterwards will enable MURR to consistently meet its experimental performance requirements. The proposed transition sequence only requires the fabrication of 34 fresh LEU elements in the first year of operation and does not exceed the anticipated availability of fresh elements that can be produced by the fuel fabricator. By the third year after conversion, 22 fresh LEU elements will be required each year, which is the same as expected for equilibrium LEU operations and the same as current operations with HEU fuel. The proposed transition sequence thus combines a relatively short time period before equilibrium burnup is achieved, a temporary increase of fuel elements needed annually relative to typical operations that are within the production capabilities of the fuel fabricator, and demonstrates comparable experimental performance of the LEU cores relative to current HEU operations. Further measures may be taken to reduce any initial experimental performance penalty even further, where possible, by repositioning certain experiments to leverage the increased performance in the lower axial experimental positions in the initial cycles following conversion or leaving the experiments in the irradiation facilities longer in order to achieve the required neutron fluence. This analysis may require refinement depending on the experimental facilities in use at the time of conversion. Nonetheless, the results presented here, including the experimental performance, core burnup, and critical control blade positions throughout the transition cycles, show that the proposed transition cycle fuel management patterns are consistent with what is expected and desired for MURR operation with LEU U-10Mo fuel. Detailed core power distributions from the neutronics models were also used to evaluate safety margins during steady-state operations for the selected transition cycles and the equilibrium LEU core. It is shown that there are adequate safety margins for both steady-state operations and postulated accident scenarios. For the steady-state operations with the preliminary LEU fuel element design the analysis predicts at least 2.49 MW margin to the onset of flow instability at the LSSS power of 15 MW. Considering the LSSS power is 125% of full license power, the margin to OFI is sufficient. In addition, the critical heat flux ratio at LSSS power is well above the requirement of CHFR > 2.0 from NUREG-1537 for all considered cases. For postulated transient accidents, the minimum margin to the fuel temperature safety limit is at least 109 °C. In summary, the proposed sequence of core loadings for MURR operations following conversion to LEU fuel and a power uprate to 12 MW provides sufficient safety margins for both steady-state operations and postulated transient accidents during a proposed sequence of transition cycles to equilibrium operations. Analysis has shown that there are some local experimental performance penalties during the initial cycles. Although there are local shifts in the experimental performance, on average all LEU cores at 12 MW have equal or higher performance than HEU at 10 MW. Temporary adjustments are being planned that will produce suitable experimental performance during these cycles. The results indicate that for the equilibrium LEU core the experimental performance exceeds that of current HEU operations in all key locations while also demonstrating sufficient safety margins.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Initial development of an RIA envelope for dispersed nuclear fuel

A reactivity insertion accident (RIA) is a design basis accident in which reactivity is rapidly injected as a result of a control rod ejection or blade drop scenario. The resultant power increase can result in fuel rod failure and subsequent release of radioactive material into the reactor’s primary system. For light water reactors, failure under RIAs is typically a result of fuel melt, pellet-cladding interaction, or rod over-pressurization. Fuel melt occurs when the fuel system is unable to transport heat out of the fuel system. Pellet-cladding interaction occurs when an aggressive fuel pellet expansion pushes on an embrittled cladding beyond its strain limit. Rod over-pressurization occurs when a fuel rod exceeds the critical heat flux and departs from nucleate boiling, causing the rod internal pressure to rapidly increase, exceeding the systems pressure, and balloon until it bursts. To prevent failure, regulators have developed safety requirements that limit the injected enthalpy based on the state of the fuel system. However, dispersed fuel could enable the removal or relaxation of regulator-imposed safety criteria. Dispersed fuel embeds fuel particles in a highly conductive metal. Dispersed particles can have a particle radius of 1 µm to greater than 100 µm. Traditional fuel systems are plagued by poor thermal conductivity, especially at higher burnups, thus reducing the ability of the fuel system to respond to an RIA. However, dispersed fuel could improve the ability of the fuel system to transport injected energy away from the fissile material and into the coolant more efficiently. Furthermore, dispersing the fuel particles mitigates hard contact that typically occurs in the traditional zirconium-uranium dioxide fuel system, thus removing pellet-cladding mechanical interaction as a potential failure mechanism. This paper evaluates fuel particles dispersed in a metal matrix using the BISON fuel performance code to develop an initial failure threshold based on melt temperature of the fuel particle and/or cladding material as a function of fission density for beginning-of-life and end-of-life conditions. BISON results show that (1) reducing the size of the particles allows fission density to increase as a function of time, (2) particle-to-particle proximity must be considered to evaluate the limiting conditions, and (3) improving the fuel particle thermal conductivity improves thermal performance.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

A study on the impact of using a subchannel resolution for modeling of large break loss of coolant accidents

The nuclear industry is investigating the feasibility of transitioning from 18- to 24-month fuel cycles because of the positive impact it would have on the operational costs for the current fleet of light-water reactors. A challenge to making this change is the increased risk of fuel fragmentation, relocation, and dispersal (FFRD) due to the known potential for ceramic fuel to pulverize into fine particles at the higher discharge burnups. Previous work has been performed by the Nuclear Energy Advanced Modeling and Simulation program to assess FFRD risk in high-burnup cores using the BISON fuel performance code and a coarse mesh thermal hydraulics (T/H) solution for a loss-of-coolant accident (LOCA) using the TRACE system T/H code. Because of the importance of the T/H solution for FFRD assessment, this study seeks to investigate the impact of using higher-fidelity subchannel techniques for modeling of the LOCA transient. CTF was used to model a subregion of a high-burnup core that was depleted by the Virtual Environment for Reactor Applications (VERA) multiphysics core simulator. Both coarse-mesh and pin-resolved models were created in CTF, and a consistent coarse-mesh TRACE model was also developed to allow for benchmarking the code results. Further, a large-break loss-of-coolant accident (LBLOCA) reflood transient was simulated using these three models, and results were compared. Results showed some consistent differences between the CTF and TRACE coarse models, including a higher peak cladding temperature (PCT) prediction in CTF and later quenching in CTF; however, the transient clad temperature behavior was similar, and these differences are likely due to post-critical heat flux heat transfer modeling differences and minimum film boiling temperature model differences. The pin-resolved results indicate that the PCT in the lumped model is often under-predicted by as much as 70 °C and that PCT occurs at a different location than the high-power pin in the assembly. The lumped model predicts a difference of 10 °C or less between the average and hot pins in the assembly, whereas the pin-resolved model predicts a range of over 100 °C. These results indicate that higher-fidelity T/H results may have an impact on predicted core behavior during LOCA, which may be important to consider when assessing FFRD risk.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

Enhanced flow boiling heat transfer on chromium coated zircaloy-4 using cold spray technique for accident tolerant fuel (ATF) materials

Flow boiling heat transfer tests were conducted to evaluate the Critical Heat Flux (CHF) and Heat Transfer Coefficient (HTC) of conventional and accident tolerant fuel (ATF) cladding materials, i.e., bare Zircaloy-4 (Zr4), and Zircaloy-4 coated with Chromium using physical vapor deposition (PVD) (Zr4-Cr-PVD) and cold spray process (Zr4-Cr-CS), respectively. The tests were performed on a single heater rod with a uniform heat flux profile, at the atmospheric pressure, inlet temperature (24 °C) and mass flow rate (750 kg/m 2 s). HTC’s were increased by 5.2% and 3.3% on Zr4-Cr-CS compared to bare Zr4 and Zr4-Cr-PVD. Improved HTCwas attributed to the increased void fraction and high roughness (Ra = 532 nm) of the Zr4-Cr-CS. In detail, the micro-cavities, which act as bubble seeds, were entirely distributed on the Zr4-Cr-CS, and these micro-cavities generated smaller and faster bubbles, thus the void fraction increased by 10.9% compared to bare Zr4. HTC is improved by these smaller and faster bubbles which can remove heat from the heater rod surface more efficiently. To identify the exact location of CHF, the surface temperature profile was measured using advanced fiber-optic sensors which have high temporal/spatial resolution (distance between point to point: ~2.5 mm, frequency ~100 Hz). CHF occurred at 80–95% along the heated length and showed an 11.6% reduction on Zr4-Cr-CS compared to bare Zr4. We hypothesize the CHF reduction was caused by the lower wettability of the Zr4-Cr-CS which reduces the liquid supply to the surface, and vigorous bubble accumulation owing to high void fraction near the outlet, which results in early formation of the vapor film. During the post-CHF quenching phase, the rougher Zr4-Cr-CS showed larger cooling rates compared to the bare Zr4 and Zr4-Cr-PVD, preventing of oxidation by chromium layer. In conclusion, we demonstrate that the high roughness on the Zr4-Cr-CS increased the number of micro-cavities on the surface, resulting in 5.2% increase in HTC and improvement in the quenching heat transfer performance whereas CHF was 11.6% reduced compared to bare Zr4. Overall statistically there is little effect in the coating with regards to CHF and slight improvement in HTC.

42 ENGINEERING↗
DOE OSTI2022

Wettability and CHF limits of Accident-Tolerant nuclear fuel cladding materials in light water reactor conditions

Here, we present the results of experimental investigations aimed at evaluating the thermal–hydraulic performance of chromium-coated zircaloy, i.e., one of the most promising accident tolerant fuel (ATF) cladding material for light water nuclear reactors. Precisely, we investigate the wettability and critical heat flux (CHF) limits of chromium-coated and conventional zircaloy surfaces in prototypical reactor conditions. For both surface types, we measure the contact angle in a vapor-saturated atmosphere from ambient pressure to the operating pressure of pressurized water reactors (PWRs), i.e., ~15 MPa. We measure the ambient-pressure steady-state flow boiling CHF with a spatially uniform heat flux. We measure the high-pressure steady-state flow boiling CHF with a cosine shape heat flux (up to 20 MPa) and with a uniform heat flux (up to 15 MPa), also exploring the effect of CRUD deposits on the chromium-coated surface. Our results reveal that the chromium surface and the bare zircaloy surface have similar wettability and both become super-hydrophilic in PWR conditions, and that there is practically no difference in the steady-state CHF limits, both at low-pressure and high-pressure conditions, also when the chromium-coated surface is covered by a CRUD deposit. However, while the chromium-coating does not improve the CHF compared to the bare zircaloy surface, it improves the post-CHF behavior. The chromium coating prevents the reaction between zircaloy and steam, which results in the formation of a brittle zirconium oxide through the surface of the cladding. We also measure the transient CHF under exponentially escalating heat flux inputs of a nano-smooth and a rough surface mimicking a chromium-coated zircaloy cladding. Interestingly, the results of the transient heat flux tests suggest that the CHF limit for very short periods (i.e., fast transients) is independent of the surface finish, being the same for a rough chromium surface or a nano-smooth surface.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Development of experimental and computational frameworks to predict subcooled flow boiling in the LANL Isotope Production Facility

Cooling is crucial to maintain the integrity of target systems in isotope production facilities. At Los Alamos National Laboratory (LANL)’s Isotope Production Facility (IPF), multiple encapsulated targets are stacked and irradiated in tandem with a 100 MeV, ~250μA proton beam. To facilitate effective heat removal, these stacked targets are separated and cooled via a series of water channels. At these beam currents, this high-energy proton beam heats the target system, likely initiating subcooled flow boiling in the cooling channels. However, in-beam monitoring of the IPF target system is not possible due to the extreme radiation environment, and the necessarily significant shielding. To better understand high-power target performance, we developed ex-situ experimental and computational frameworks to predict the behavior of subcooled flow boiling at IPF. Subcooled flow boiling experiments on Inconel 625 samples under IPF conditions (2 bar pressure, 10 GPM flow rate (i.e., 2249 kg/m 2 /s), 85 K subcooling) revealed that IPF's average operating power is at the early stage of boiling with a heat transfer coefficient of 48,000 W/m 2 /s. The proposed modeling framework enables us to predict a complete boiling curve, i.e., single-phase heat transfer, onset of nucleate boiling, two-phase heat transfer, and critical heat flux (CHF), with specification of input boiling parameters up to intermediate heat flux levels. The estimated CHF under IPF conditions is 5.2 MW/m 2 . Experimental data under reduced conditions (2 bar pressure, 1.5 GPM flow rate (i.e., 337 kg/m 2 /s), 45 K subcooling) served as validation cases for the computational modeling. This computational model can be further extended to more complicated systems replicating the real IPF configuration, for instance, to study void distribution as a function of the incident proton beam profile and coolant velocity profile of multiple cooling channels. Finally, the proposed experimental and computational frameworks provide a means to better understand cooling systems in the isotope production facilities at different accelerators, where in-beam monitoring of the cooling process is not available.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2021

Improved departure from nucleate boiling prediction in rod bundles using a physics-informed machine learning-aided framework

The critical heat flux (CHF) corresponding to the departure from nucleate boiling (DNB) crisis is a regulatory limit for the licensing of pressurized water reactors (PWRs) worldwide. Despite the abundance of predictive tools available to the reactor thermal-hydraulics community, the path for an accurate CHF model remains elusive. This work approaches the prediction of DNB through a physics-informed machine learning-aided framework (PIMLAF) with the objective of achieving superior predictive capabilities for a rod bundle. In view of the limitations in existing macro-scale physics-driven tools, an improved mechanistic model is first proposed, leveraging key concepts in the liquid sublayer dryout and bubble crowding mechanisms. Furthermore, the proposed mechanistic model is able to predict DNB in different heater geometries for a broad range of flow conditions without the need for recalibration. This model is then incorporated as the physics-informed component of the hybrid framework PIMLAF, which takes advantage of established understanding in the field (i.e., domain knowledge [DK]) and uses machine learning (ML) to capture undiscovered information from the mismatch between the actual and DK-predicted output. Two bundle-related case studies using the PWR subchannel and bundle tests (PSBT) database are carried out to illustrate the PIMLAF’s improved performance over traditional approaches for both interpolation and extrapolation purposes. In light of the PIMLAF’s promising potential to reduce prediction error, reactor vendors are encouraged to leverage their in-house experimental efforts and apply the hybrid framework to potentially achieve margin reductions in the minimum DNB ratio (MDNBR) for the designs of interest.

42 ENGINEERING↗
DOE OSTI2024

Flow Instabilities in boiling channels and their suppression methodologies—A review

Small modular reactors (SMRs) are gaining significant attention as a promising solution for clean and sustainable nuclear-power generation. However, the operation of SMRs is subject to various challenges, including two-phase flow instabilities. Flow instability has the potential to trigger flow-induced vibration and cyclic fluctuations in local thermal stress. These instabilities frequently manifest because of the complex interplay among a multitude of factors, encompassing thermal-hydraulic conditions, the geometric configuration of the steam generator, and operational parameters. These conditions could subsequently lead to premature critical heat flux, equipment malfunctions, and other safety concerns. Here, the endeavor to address steam-generator flow instabilities is of utmost importance in augmenting the sustainability and efficiency of contemporary energy production. This study offers a comprehensive review of instabilities in two-phase flow, with a particular focus on the influential factors impacting the stability of flow boiling. Furthermore, it delves into the processes of identifying, characterizing, and ameliorating these instabilities, emphasizing pivotal findings, methodologies employed, and avenues for prospective research. The primary parameters of concern encompass the efficient transfer of thermal energy, the optimization of mass-flow rates, and the establishment of favorable boundary conditions, all in the context of steam generator design to alleviate instability for water-cooled SMRs. These discernments bear substantial ramifications for enhancing reactor performance and ensuring operational safety.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
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 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 OSTI2020

Investigations on the Thermal-hydraulic Behavior of Accident Tolerant Fuel Cladding Materials

Since the 2011 Fukushima accident, significant research has been devoted to developing accident-tolerant fuel (ATF) cladding materials. These investigations have mostly focused on the ability to resist runaway steam oxidation and retain mechanical strength and structural integrity under thermal shocks. However, much remains unknown about the materials’ thermal-hydraulic behavior, particularly under light-water reactor (LWR) operating conditions. Two phenomena that determine safety margins in both normal and off-normal operating conditions - i.e., the critical heat flux (CHF) and Leidenfrost point (LFP) temperature - have not been investigated thoroughly for ATF materials. Surface properties, such as wettability and roughness, are known to influence pool and flow boiling CHF, as well as the LFP temperature. However, little is known about the surface wettability of ATF materials, particularly at LWR pressure and temperature. Little has been done on flow-boiling CHF and quenching heat transfer, especially for droplet quenching. Nonetheless, it is known that CHF in transient conditions, e.g., an exponentially escalating power transient, can be significantly different from that expected in steady-state operation. Very few studies have investigated transient CHF, either on ATF materials or under LWR pressure and temperature. The thorough understanding of transient CHF under prototypical reactor conditions will benefit not only the deployment of ATF but also the upcoming national efforts to study transient ATF behavior in the Transient Reactor Test Facility (TREAT) at Idaho National Laboratory. In summary, as a part of ATF development efforts, there is an urgent need to understand how and how much these materials may affect two-phase heat transfer phenomena in nuclear reactor conditions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Thermocatalytic Heat Pipes for Geothermal Resource Recovery

Heat pipes are an important technology that allow orders of magnitude faster heat transfer than simple conduction. However, operating principles in heat pipes place fundamental bounds on their performance (critical heat flux and efficiency). Conventional heat pipe functionality is inherently tied to vaporization and condensation of the working fluid charged in the heat pipe. These fluids each have different operating temperature ranges based on the capillary, entrainment, sonic, and boiling limits of the heat pipe design. These limits, typically the capillary limit, dictate the maximum heat flux a heat pipe can carry, and most importantly for geothermal systems, the distance over which the pipes can operate (100 to 200 m maximum under optimum conditions). A thermocatalytic heat pipe breaks the inherent limitations of phase change thermo- and hydrodynamics and can transform heat pipe technology as a potentially more efficient means of extracting heat from a geothermal resource. The thermocatalytic heat pipe uses a working fluid to transport both sensible and chemical heat. An endothermic chemical reaction at depth removes heat from the reservoir and produces reactive intermediates, which are transported to the surface and used to run a reverse exothermic reaction that releases heat for use in power generation or other useful purposes. This technology offers two distinct advantages over conventional geothermal heat recovery technologies: (1) lower heat loss to the rock outside of the geothermal reservoir, and (2) higher heat transfer rates to the well field within the geothermal reservoir. Both advantages offer opportunity to reduce risks and lower costs of geothermal energy recovery. In this report, we discuss an initial effort to assess the efficacy and limitations of this technology for extracting heat from both porous/permeable and nominally impermeable geothermal reservoirs. Numerical simulation capabilities of the STOMP-GT code were enhanced to enable simulations of thermochemical heat pipes traversing geothermal reservoirs. An array of potential thermochemical reaction systems was evaluated and screened. Of these, an ethanol dehydration reaction was most promising in the vapor-liquid reaction set. A solid-phase dehydration reaction (CuSO4·5H2O) showed the highest reaction enthalpy per unit volume but would require development of a nonaqueous carrier fluid to implement it in a heat pipe. Subsurface reservoir simulations predicted long-term performance of the heat pipes for each geothermal reservoir type. The performance of U-shaped wells and coaxial wells was evaluated for a suite of reactions for both hydrothermal and hot dry rock reservoirs and was compared with a baseline case of simply pumping water through the wells. The heat pipe technology was additionally evaluated for an enhanced geothermal system (EGS) with an injection borehole, production borehole, and intervening hydraulically conductive fracture. All reservoir types showed significant improvement in heat recovered over a 20-year operating period ranging from a 1.8X increase for the hot dry rock case to more than 2.5X more energy recovered for the EGS case.

15 GEOTHERMAL ENERGY↗
DOE OSTI2020

The STAT7 Code for Statistical Propagation of Uncertainties In Steady-State Thermal Hydraulics Analysis of Plate-Fueled Reactors

STAT7 was written to automate many of the steady-state thermal hydraulic safety calculations for the MIT research reactor, both for conversion of the reactor from highly enriched uranium fuel to low-enriched uranium fuel and for future fuel re-loads after the conversion. A Monte-Carlo statistical propagation approach is used to treat uncertainties in important parameters in the analysis. These safety calculations are ultimately intended to protect against high fuel plate temperatures due to critical heat flux or departure from nucleate boiling or onset of flow instability; but additional margin is obtained by basing the limiting safety settings on avoiding onset of nucleate boiling. STAT7 can simultaneously analyze all of the axial nodes of all of the fuel plates and all of the coolant channels for one stripe of a fuel element. The stripes run the length of the fuel, from the bottom to the top. Power splits are calculated for each axial node of each plate to determine how much of the power goes out each face of the plate. By running STAT7 multiple times, full core analysis can be performed by analyzing the margin to ONB for each axial node of each stripe of each plate of each element in the core.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗