Search NASA⌕ Search

SEARCH · Search NASA

Results for “thermohydraulics”

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

Thermohydraulic experiments on a supercritical carbon dioxide–air microtube heat exchanger

Heat exchangers are critical components in supercritical CO 2 Brayton cycles and provide necessary heat input, recovery, and dissipation. Tubular heat exchangers with unconventionally small tube sizes (tube diameters less than 5 mm) are promising components for supercritical CO 2 cycles and potentially provide excellent structural stability with wide scope of application. This paper provides essential design and fabrication guidelines for a compact shell-and-tube heat exchanger with microtubes (with an inner diameter equal to 1.75 mm). A heat exchanger test rig is used to evaluate the thermohydraulic performance of this heat exchanger with supercritical CO 2 and air as working fluids. Thermohydraulic data are reported for more than forty sets of experiments with varying Reynolds numbers for shell and tube flows. Critical performance metrics are calculated from the data and compared with predictions from a previously described numerical model. The average deviations between the experimental and model results fall within 10% for all critical metrics. This excellent agreement validates the numerical model for supercritical CO 2 heat exchanger optimization and scale-up.

42 ENGINEERING↗

A Novel Approach for Bounding the Stress Experienced by the Core of Utility-Scale Printed Circuit Heat Exchangers Under Thermohydraulic Loads

Printed circuit heat exchangers (PCHEs) have rapidly gained popularity since being introduced nearly three decades ago, and they are currently widely deployed in the petrochemical and aviation industry. Their compactness, thermohydraulic efficiency, inherent suitability for high temperature/pressure fluids containment, and demonstrated durability are some of the reasons the nuclear industry is seeking to adopt this technology as well. However, the relatively strict nuclear-related regulatory design code, especially when classified as critical to the safety of the reactors, is posing challenges to adopting the technology. From stress analysis point of view, one undesirable feature of PCHEs is their geometrical complexity, which is implied by their multilength-scale features. As a result, a full-scale model of a utility-scale exchanger cannot simply be solved on a computer because meshing such components results in a vast number of degrees-of-freedom. Here, this work seeks to address the challenge of stress analyses to PCHEs by presenting a method to simplify the geometry of PCHE designs. The models proposed by this work can be practically analyzed on a standard computer and provide a path for implementing ASME design rules. The analyses presented herein are divided into five separate investigations. Each is carried out to incrementally simplify the analyzed model by addressing features such as the shapes of the flow passages, the complex distribution of stress in large components, the three-dimensionality of the stress and strain, the thermal stresses caused by thermohydraulic operation observed experimentally and more.

42 ENGINEERING↗

Physical Modeling of Coupled Thermohydraulic Behavior of Compacted MX80 Bentonite during Heating

Abstract This paper describes a tank-scale test setup and experimental methodology used to investigate coupled heat transfer and water flow processes during heating of compacted MX80 bentonite to high temperatures. Specifically, a temperature of 200°C was maintained by a cylindrical heating element at the center of a compacted bentonite layer containing an array of temperature, dielectric, and relative humidity sensors. In addition to providing an evaluation of the spatiotemporal variations in temperature, relative humidity, degree of saturation, and global volume, the coupled thermohydraulic properties of the bentonite were assessed. A wetting front was initially observed to move away from the central heater, followed by a drying process until reaching thermohydraulic equilibrium. The soil–water retention curve (SWRC) of the bentonite followed a wetting scanning path before following the primary drying path exhibiting a shift in water retention with elevated temperature. Results from the tank-scale test can be used for validation of numerical simulations of drying processes in the engineered barrier system of a high-level radioactive waste geological disposal repository and confirm that a temperature-dependent hysteretic SWRC with scanning paths is required to accurately capture the bentonite response.

Engineering↗

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↗

Thermohydraulic Design Analysis of the Target Assembly in the Material Plasma Exposure Experiment Facility

The Material Plasma Exposure eXperiment (MPEX) project seeks to design a steady-state linear plasma facility at Oak Ridge National Laboratory that will be used to study plasma-material interactions (PMIs) at fusion prototypic levels, supporting the evaluation and development of materials for the next generation of fusion devices. This study is focused on PMI exposure of small-size neutron-irradiated specimens, which are clamped onto an actively cooled component. A thermohydraulic evaluation of a new MPEX target assembly design to assess the appropriate operation during MPEX operation is presented. To further guide the design and assess the structural integrity of the components under expected loads, preliminary thermomechanical stress analyses were also conducted. To ensure good thermal contact between the components, thermal interface materials, such as silver flexible graphite, were used in the assembly. It was found that the maximum target temperatures of 1572, 1463, and 1315 K were obtained for Grafoil thicknesses of 0.61, 0.38, and 0.25 mm, respectively. The distribution of the axial deformation at high heat fluxes showed that there are no gaps between components, indicating good contact at material interfaces. Moreover, the contact pressure between the target and other components indicated that very good contact was established at these interfaces. Furthermore, the stress-strain conditions for the target will be further used to assess the appropriate operation during MPEX experiments and gain insight into materials science phenomena during PMI experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Effect of Inlet Throttling on Thermohydraulic Instability in a Large Scale Waterbased RCCS: A System-Level Analysis with RELAP5-3D

This paper presents results from system -level modeling of a water -based reactor cavity cooling system using RELAP5-3D. The computational model is benchmarked with experimental data from a half -scale RCCS test facility at Argonne National Laboratory. The model prediction is first compared with a two-phase oscillatory baseline experimental case where mixed accuracy is obtained. The model shows reasonable prediction of mass flow rate, pressure, and temperature but significant overprediction of void fraction. The model prediction is then compared with a fault case where the inlet of the risers is gradually reduced using a throttling valve. As the valve is closed, the model is able to predict some major flow phenomena observed in the experiment such as the dampening of oscillations, the reintroduction of oscillations, as well as boiling, flashing, and geysering in the risers. However, the timeline of these events are not well captured by the model. The model is also used to investigate the evolution of flow regime in the chimney. This work highlights that the semiempirical constitutive relations used in RELAP-3D could have a strong influence on the accuracy of the model in two-phase oscillatory flows.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Validation of Thermohydraulic Simulations using RELAP for Critical Dual Purpose Canisters

Simulations using the RELAP5-3D systems analysis code and STAR-CCM+, a high-fidelity CFD code, were performed on a model of a real-scale MPC-32 dual purpose canister undergoing a criticality event. These simulations were performed to provide code-to-code verification of the natural convection capabilities in RELAP5. The STAR-CCM+ simulations revealed various natural convection flow patterns that aided in the cooling of fission heat produced in the fuel rods in the canister. Negligible spatial temperature variations were observed across the canister in the STAR-CCM+ results, supporting the use of RELAP5 in predicting the sub-cooled regime. The surface-averaged rod temperature predicted by RELAP5 compared well with STAR-CCM+, verifying the predictive capability of RELAP5 for sub-cooled conditions.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

A fully coupled variable properties thermohydraulic model for a cryogenic hydrostatic journal bearing

The goal set forth here is to continue the work started by Braun et al. (1984-1985) and present an integrated analysis of the behavior of the two row, 20 staggered pockets, hydrostatic cryogenic bearing used by the turbopumps of the Space Shuttle main engine. The variable properties Reynolds equation is fully coupled with the two-dimensional fluid film energy equation. The three-dimensional equations of the shaft and bushing model the boundary conditions of the fluid film energy equation. The effects of shaft eccentricity, angular velocity, and inertia pressure drops at pocket edge are incorporated in the model. Their effects on the bearing fluid properties, load carrying capacity, mass flow, pressure, velocity, and temperature form the ultimate object of this paper.

Braun, M. J.↗

Transient modeling of the thermohydraulic behavior of high temperature heat pipes for space reactor applications

Many proposed space reactor designs employ heat pipes as a means of conveying heat. Previous researchers have been concerned with steady state operation, but the transient operation is of interest in space reactor applications due to the necessity of remote startup and shutdown. A model is being developed to study the dynamic behavior of high temperature heat pipes during startup, shutdown and normal operation under space environments. Model development and preliminary results for a hypothetical design of the system are presented.

Hall, Michael L.↗

Multi-physics Topology OPtimization and Additive Manufacturing for High-temperature Heat Exchangers

This research significantly advances the understanding of high-temperature heat exchanger design through an integrated approach that combines topology optimization (TO), triply periodic minimal surface (TPMS) structures, additive manufacturing (AM) and thermohydraulic testing. Each of these components contributes uniquely to a unified, high-performance design, fabrication and testing workflow. Topology optimization serves as the foundation of the design methodology by providing a systematic way to determine the most effective material layout for separating hot and cold fluids while maximizing thermal performance. The researchers introduced a novel three-material optimization framework using two density fields to represent hot fluid, cold fluid, and solid domains. This approach enables automated discovery of optimal shapes and flow paths that cannot be intuitively designed, especially under constraints imposed by manufacturing technologies. Furthermore, constraints such as minimal wall thickness and overhang angles were embedded into the optimization process, ensuring that resulting designs are not only thermally efficient but also manufacturable using modern additive techniques. In parallel, the study delves into the use of Gyroid-based TPMS geometries for constructing the core of the heat exchanger. TPMS structures are known for their high surface area, excellent fluid mixing capabilities, and minimal pressure drop characteristics. The researchers applied a data-driven modeling framework using Heteroscedastic Sparse Gaussian Process Regression (HSGPR) combined with genetic algorithms. This allowed for the rapid evaluation and optimization of key geometric parameters such as frequency, iso-value, and phase shift. The result was a set of Gyroid structures tailored for high heat transfer and low flow resistance, demonstrating clear improvements over conventional straight-channel designs. After the designing process, additive manufacturing played a critical role by turning these highly complex, optimized geometries into physical components. Utilizing Laser Powder Bed Fusion (LPBF) with Haynes 282, the study demonstrated the feasibility of fabricating these heat exchangers at high precision. Post-processing methods, including dilation-erosion operations, were applied to ensure local features adhered to self-supporting constraints. The fabricated structures were then subjected to thermohydraulic testing under conditions representative of supercritical CO 2 Brayton cycles, validating the predicted performance and confirming the viability of the full design-to-fabrication pipeline. Finally, thermohydraulic testing across the above studies served as a crucial experimental validation of advanced heat exchanger. Under consistent high-temperature and high-pressure conditions using supercritical CO 2 , the testing demonstrated that both TO and Gyroid-based TPMS designs significantly outperformed conventional straight-channel HXs. The TO design achieved a 115% increase in UA and NTU and a 27.6% boost in gravimetric power density, while the data-driven optimized Gyroid design delivered a 166% increase in UA and NTU and improved effectiveness from 68.7% to 86.1%. These results validate the simulation models, confirm the manufacturability of complex geometries under AM constraints, and provide key insights into design-performance trade-offs, thereby advancing the development of high-efficiency, compact heat exchangers for extreme environments.

36 MATERIALS SCIENCE↗

Modeling of Supercritical CO2 Shell-and-Tube Heat Exchangers Under Extreme Conditions. Part I: Correlation Development

Abstract High-temperature supercritical CO2 Brayton cycles are promising possibilities for future stationary power generation and hybrid electric propulsion applications. Heat exchangers are critical components in supercritical CO2 thermal cycles and require accurate correlations and comprehensive performance modeling under extreme temperatures and pressures. In this paper (Part I), new Colburn and friction factor correlations are developed to quantify shell-side heat transfer and friction characteristics of flow within heat exchangers in the shell-and-tube configuration. Using experimental and computational fluid dynamics (CFD) data sets from existing literature, multivariate regression analysis is conducted to achieve correlations that capture the effect of multiple critical geometric parameters. These correlations offer superior accuracy and versatility as compared to previous studies and predict the thermohydraulic performance of about 90% of the existing experimental and CFD data within ±15%. Supplementary thermohydraulic performance data are acquired from CFD simulations with supercritical CO2 as working fluid to validate the developed correlations and demonstrate its capability to be applied to supercrtical CO2 heat exchangers.

Engineering↗

PIPII cryoplant non thermal cycling updates

The Proton Improvement Plan-II (PIP-II) is a crucial upgrade to the Fermilab accelerator complex, featuring a new 800-MeV Superconducting Radio-Frequency (SRF) linear accelerator (LINAC) with 23 cryomodules operating at 2K. The LINAC thermohydraulic conditions are satisfied by the cryogenic subsystems: Cryogenic Distribution System (CDS), a helium refrigerator cold box (CB), a warm compression station (WCS) and a helium recovery system (RSYS). The accelerator has a strict requirement of non-thermal cycling of the LINAC cryomodules during planned and unplanned subsystem outages. This paper presents an integrated operating modes analysis of the of the LINAC/CDS thermohydraulic loads satisfied by the CB/WCS cooling system supported by the RSYS inventory management system. The study is based on latest cryoplant and CDS engineering deliverables, as well as the recent performance data from single cryomodule qualification tests performed at the Fermilab PIP-II Injector Test test s tand. The study evaluates both normal and abnormal operating modes, with a focus on identifying integrated scenarios that put subsystems components under stress. The conclusions of this study will help build redundancy to reduce the risk of thermal cycles during planned and unplanned subsystem outages.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Native Natural circulation in VERA

This work seeks to introduce a native natural circulation solver to VERA's thermohydraulics capabilities. This work begins with initial investigations into the viability of such a solver by determining the fundamental equations as well as necessary inputs. The method was then implemented in the form of a Python script coupled to CTF which was compared to TRACE. Once viability of the natural circulation solver was established, it was then implemented natively in MPACT for use with previously created simplified thermohydraulics solvers. This new implementation of the native natural circulation solver in MPACT is then tested for single physics problems for both a BWR and PWR single channel model and compared to TRACE where it shows an error of about 30% for BWRs and about 15% for PWRs. The solver is additionally tested for NuScale based models where it demonstrates an error to the reported flow rate of about 10%, resulting in an eigenvalue error of only 138 pcm when comparing fixed flow rate multiphysics runs. Finally, the solver is tested in multiphysics systems for a small PWR assembly, a small BWR assembly, and a full core approximate model of the NuScale design. In these multiphysics runs, it demonstrates no impact on multiphysics convergence stability, as well as negligible impact on runtimes. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Multiphysics modeling of a critical dual-purpose canister in a saturated geological repository

Dual purpose dry cask storage canisters for spent nuclear fuel are designed for storage and transportation, but are not licensed for permanent disposal in a geological repository. If dual purpose canisters were to be used to dispose of spent nuclear fuel in a geological repository, they would be expected to eventually breach and be flooded with groundwater, and it is shown that some fraction of these canisters will achieve criticality. In order to evaluate the consequences of canisters going critical in a repository, an initial capability has been developed for estimating the quasi-static power level of a critical canister using loosely coupled multiphysics simulations. The low power level in a critical canister enables coupling through precomputed physics proxies. This calculated power level is then used to compute the change in the critical canister’s isotopic inventory as a function of time. Three as-loaded canisters are evaluated and two were found to have power levels below 4 kW, with a modest effect on the radiological inventory over time. This effort also shows that although some DPCs will have extremely peaked power shapes, the relatively low power and long-time scales result in relatively homogenized thermohydraulic properties in the water within the DPC.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

srlife : a software tool for estimating the life of high temperature concentrating solar receivers. Part I – metallic receivers

Here, this paper introduces srlife, a tool for estimating the structural service life of concentrating solar power (CSP) receivers operating at high temperatures. Supporting both metallic and ceramic receiver designs, srlife is available as open-source software at https://github.com/applied-material-modeling/srlife and can be installed via the PyPi package manager (https://pypi.org). Given basic receiver geometry and incident heat flux, the tool performs thermohydraulic and structural analysis and estimates the life of a receiver. Designed for easy integration into a software stack, including solar field and levelized cost analysis, the tool can be utilized for optimizing receiver designs to meet service life and economic targets. This paper is Part I in a two-part series. Part I discusses the analysis process used to estimate the life of metallic receivers, along with a description of the required input data. Additionally, several heuristics applied within srlife can reduce analysis time significantly while maintaining accurate life estimations for metallic receivers when compared to full analyses. Several examples demonstrating the utility of srlife in receiver design are also discussed. Part II focuses on the life estimation of ceramic receivers, using time-dependent reliability analysis and various ceramic failure models implemented in srlife.

Creep-fatigue analysis↗

srlife : A software tool for estimating the life of high temperature concentrating solar receivers. Part II – Ceramic receivers

As Concentrating Solar Power (CSP) technologies aim for higher operating temperatures to enhance efficiency and meet industrial process heat demands, high-temperature metallic materials, including nickel-based superalloys, face challenges in maintaining structural integrity. Advanced ceramics offer a promising alternative due to their superior high-temperature strength. However, accurately assessing the performance of ceramic components requires a fundamentally different approach from that used for metallic components. This Part II of a two-part paper describes the integration of ceramic statistical failure models within srlife – an open-source tool for predicting the life of high-temperature CSP receivers. These models account for the inherent variability in ceramic strength, as well as the effects of subcritical crack growth (SCG) under high temperature cyclic loads. Here, the paper includes an example problem that demonstrate the process of evaluating ceramic receivers using srlife. Part I details the life estimation process for metallic receivers (i.e. creep-fatigue life) along with input and output data structure, thermohydraulic analysis, and structural analysis. The complete tool is available as open-source software at https://github.com/srlife-project/srlife and can be installed via the PyPi package manager (https://pypi.org). By supporting both ceramic and metallic receiver analyses, srlife facilitates fair comparisons between competing metallic and ceramic designs, enabling accurate evaluations of plant efficiency and the economic benefits of ceramic solar receivers and other components.

High temperature ceramic receivers↗

A RELAP5-3D Model of the Lobo Lead Loop

The goal of our research is to build upon the capability of RELAP5-3D to model molten lead systems. Molten lead has several potential uses in future advanced reactors like the lead fast reactor or fusion reactors that utilize dual-coolant lead lithium blankets. This potential for use in future generations of reactors highlights the necessity to develop molten lead models to ensure that they can accurately predict the thermohydraulic behavior. We have developed a RELAP5-3D model of the Lobo Lead Loop facility, located at the University of New Mexico, to verify the accuracy of RELAP5-3D via comparison to existing computational fluid dynamics results and analytical calculations. It was found that RELAP5-3D accurately calculated radiative heat transfer (within < 1%) when compared to theoretical calculations. In addition, pressure drop calculations done in RELAP5-3D demonstrated reasonable agreement within 20 kPa, mostly within ~7-15%, when compared to the computational fluid dynamics model of the facility developed by the University of New Mexico, and captured the dependence of pressure drop on flow velocity accurately. Finally, a hypothetical loss of flow transient was imposed on the RELAP5-3D model to determine the feasibility of performing a similar experiment with the Lobo Lead Loop. It was found that such an experiment could be possible as the RELAP5-3D model indicated that the temperatures of the fluid would not exceed limiting temperatures of the structure (1658 K) nor the maximum temperature of the electromagnetic pump inlet (823 K). Although there is not experimental data to begin validation, the model will be readily available for future validation studies when the experimental data is generated, especially as the model continues to evolve over time. Furthermore, the results so far demonstrate a promising first step in the verification/validation of the RELAP5-3D model of the Lobo Lead Loop.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Numerical investigation of two-dimensional buoyancy-driven eddies in liquid metal magnetohydrodynamic flows in breeding blankets *

Abstract Lead–lithium flows are key features in the design of tokamak breeding blanket concepts such as the dual-coolant lead–lithium (DCLL). Since they flow under magnetic fields, they are affected by magnetohydrodynamic (MHD) effects. The neutron flux originating in the tokamak plasma heats the breeding blanket channels in a non-uniform manner, inducing buoyancy forces in the liquid metal. Buoyancy may become a source of quasi-two-dimensional (Q2D) turbulence, and the appearance of eddies may affect the transport of heat and tritium across the blanket. Blankets characterized by high-speed liquid metal flows (such as DCLL) will need ceramic insulating walls to reduce the MHD-related pressure drop in the channels. In our simulations, we have used the Q2D model proposed by Sommeria and Moreau (SM82) which is especially suitable for modeling electrically insulating channel flows. On top of that, we have modelled buoyancy forces in the momentum equation using the Oberbeck–Boussinesq approximation. In this work, we include a validation of the implemented Q2D model in buoyancy-driven cases and identify a thermohydraulic configuration that promotes the generation of eddies and the accumulation of tritium. We also show the results provided by our post-processing tool based on the bi-dimensional fast Fourier transform for eddy detection and characterization. We complete our investigation by performing an initial assessment of how eddies can accumulate tritium in the breeding blanket. We conclude by discussing the relevance of preventing accumulation of tritium within the blanket and suggest a possible solution.

Suarez, Daniel (ORCID:0000000251740100)↗