Investigation of transient flow boiling heat transfer physics and system-level thermal-hydraulic responses during line chilldown
Explore the source record for details and available documents.
Jiayi Zhang↗
SEARCH · Search NASA
Results for “flow boiling”
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.
Comparison of Cryogenic Flow Boiling in Liquid Nitrogen and Liquid Hydrogen Chilldown Experiments
This paper presents a comparison between experimental results from recent liquid hydrogen (LH 2 ) transfer line chilldown experiments at high Reynolds (Re) numbers versus liquid nitrogen (LN 2 ) experiments conducted at low Re numbers. Parasitic heat leak, inner wall temperatures, inner wall heat fluxes, and heat transfer coefficients are computed to compare between the two systems. Analysis of temperature traces and flow visualization indicates that the chilldown process evolves much more rapidly at higher Re numbers due to a quick transition from vapor flow to annular liquid flow and near immediate liquid contact along the pipe walls. The lower kinematic viscosity and surface tension of LH 2 , along with reduced parasitic heat leak and higher Re numbers relative to the LN 2 experiments, causes chilldown to proceed almost immediately into the nucleate boiling regime, in comparison to low Re flows where >75% of the chilldown is spent in vapor film boiling.
Cryogenic Heat Transfer Coefficient↗
Cryogenic Flow Boiling in Microgravity: Effects of Reduced Gravity on Two-Phase Fluid Physics and Heat Transfer
With the growing interest in space exploration, cryogenic technologies involving two-phase flow and heat transfer are in high demand to successfully procure advanced space applications such as fuel depots and nuclear thermal propulsion (NTP) systems for deep space missions. However, the unique and extreme thermal properties of cryogenic fluids introduce distinct flow boiling fluid physics and energy transport phenomena, which differ significantly from those observed with conventional fluids. Understanding the unique two-phase physics in cryogenic flow boiling remains an ongoing challenge. Furthermore, the lack of readily available microgravity cryogenic steady-state heat transfer data hinders the assessment of gravitational effects on cryogenic flow boiling. This study aims to elucidate the gravitational effects on two-phase fluid physics and heat transfer by conducting the first-ever experimental measurement of cryogenic flow boiling performance using a steady-state heated method in a reduced gravity environment. Parabolic flight experiments were performed to acquire both heat transfer measurements and high-speed video of interfacial behaviors, under varying gravity levels (microgravity, hypergravity, Lunar gravity, and Martian gravity). The experiments involved flow boiling of liquid nitrogen (LN 2 ) with a near-saturated inlet along a circular heated tube of dimensions 8.5-mm inner diameter and 680-mm heated length. The operating parameters varied are mass velocity of 398.3 - 1342.8 kg/m2s, inlet quality of -0.08 to -0.01, and inlet pressure of 413.68 - 689.48 kPa. Captured microgravity flow patterns range from bubbly to annular, all having vapor structures that are larger than those under higher gravity levels. Under microgravity, absence of buoyancy yields symmetrical vapor structures without flow stratification, laying a physical foundation for the distinct two-phase heat transfer trends during LN 2 flow boiling in microgravity. Transient data collected during the flight parabolas exhibited decreasing heated wall temperature as the aircraft transitioned from hypergravity to microgravity phases. The temperature variation indicated an enhancement in flow boiling heat transfer with decreasing gravity levels and a reduction with increasing gravity levels. The effect of reduced gravity on cryogenic flow boiling heat transfer coefficient (HTC) is discussed based on steady state heat transfer analysis. Seminal HTC correlations are evaluated against the measured microgravity HTC data, of which one is identified for superior accuracy in predicting microgravity data. Finally, a new HTC correlation is proposed to improve accuracy of microgravity predictions, yet there still exists room for further improvement with future terrestrial flow boiling experiments at different flow orientations relative to Earth gravity.
Microgravity↗
Rewet Temperature Correlations for Liquid-Nitrogen Boiling Pipe Flows Across Varying Flow Conditions and Orientations
In many convective liquid-vapor phase change heat transfer engineering applications, cryogenic fluids are widely used in industrial processes, spacecraft and cryosurgery systems, and so on. For example, cryogens are usually used as liquid fuels such as liquid hydrogen and oxygen in the rocket industry, liquid nitrogen (LN2) and helium are frequently used to cool superconducting magnetic device for medical applications. In these systems, proper transport, handling, and storage of cryogenic fluids are of extreme importance. Among all the cryogenic transport processes performed in room temperatures, quenching, also termed chilldown, is a unavoidable initial, transient phase-change heat transfer process that brings the system down to the cryogenic condition. The Leidenfrost temperature or rewet temperature that signals the end of film boiling is practically considered the completion point of a quenching process. Therefore, rewet temperature has been considered the most important parameter for the engineering design of cryogenic thermal management systems. As most of the previous correlations for predicting the Leidenfrost temperature and the rewet temperature have been basically developed for water, they are shown to disagree with recent liquid nitrogen pipe chilldown experiments in upward and downward flow directions over a wide range of flow rates, pressures, and degrees of inlet subcooling. In addition to a complete review of the literature, two new correlations are presented in this work, one based on bubble growth and another based on the theoretical maximum limit of superheat. Each correlation performs well over the entire data set.
liquid nitrogen↗
Numerical Predictions of the Flow and Heat Transfer Characteristics in the Film Boiling Regime During Tube Quenching
Cryogenic fluid management plays a major role in refueling of spacecrafts while in space for NASA’s future human space exploration missions. Due to the low boiling points of cryogens, storage, transport and handling of these fluids becomes difficult and may result in inefficient operation of the space propulsion systems. For refueling applications in space, the cryogenic fluids have to be transported across different locations and hence, the transfer of cryogenic fluids through pipes become critical. The cryogenic chill-down process is characterized by different regimes of flow boiling, viz., film boiling, transition boiling and nucleate boiling. The prediction of these regimes in a single CFD framework available in the literature is challenging and the present work attempts to address this challenge by initially modeling the film boiling regime accurately and to incorporate an user-defined function for transition and nucleate boiling at a later stage. Hence, the aim of the present work is to numerically model and validate the film boiling regime of the chilldown curve for liquid nitrogen experiments available in the literature. The validations are carried out at different inlet mass fluxes to have a robust simulation methodology. A dispersed mixture model is used to predict the vapor-liquid interface dynamics with the phase change phenomena modeled using the Lee model.
line chilldown↗
Heat Transfer Enhancement in Cryogenic Quenching Process
This paper reports a heat transfer advancement in the cryogenic quenching process. An experiment was performed to evaluate the enhancement of quenching heat transfer by the use of metal tubes with low thermal conductivity coating layers. Four coating thicknesses with various coolant mass flow rates of liquid nitrogen were investigated. The results indicated that the tube inner surface coating greatly enhanced the quenching efficiency. The quenching efficiency was found to increase with increasing number of coating layers, and the efficiency also increased with decreasing mass flow rates. In general, the efficiencies cover a range between 40.6% and 80%. Comparing to the bare surface case, the percentage increase in the quenching efficiency was the minimum at 4.2% for a single coated layer at the highest flow rate and the maximum of 109.1% for four coated layers at the lowest flow rate. The coated tubes could save up to 53% in the amount of cryogen consumption.
two-phase flow↗
Numerical Modeling of the Transient Chilldown Process of a Cryogenic Propellant Transfer Line
Before cryogenic fuel depots can be fully realized, efficient methods with which to chill down the spacecraft transfer line and receiver tank are required. This paper presents numerical modeling of the chilldown of a liquid hydrogen tank-to-tank propellant transfer line using the Generalized Fluid System Simulation Program (GFSSP). To compare with data from recently concluded turbulent LH2 chill down experiments, seven different cases were run across a range of inlet liquid temperatures and mass flow rates. Both trickle and pulse chill down methods were simulated. The GFSSP model qualitatively matches external skin mounted temperature readings, but large differences are shown between measured and predicted internal stream temperatures. Discrepancies are attributed to the simplified model correlation used to compute two-phase flow boiling heat transfer. Flow visualization from testing shows that the initial bottoming out of skin mounted sensors corresponds to annular flow, but that considerable time is required for the stream sensor to achieve steady state as the system moves through annular, churn, and bubbly flow. The GFSSP model does adequately well in tracking trends in the data but further work is needed to refine the two-phase flow modeling to better match observed test data.
Temperature Sensors↗
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↗
Characterization Testing of a 20 Kelvin Cryocooler for Space Applications
Future NASA space exploration missions will require long-duration storage and liquefaction of cryogenic liquids, enabled by active cooling provided by cryocoolers. Recent gap analyses of Lunar and Mars transportation systems have identified 20 K-class cryocoolers as a critical enabling technology for chemical and nuclear thermal propulsion architectures using liquid hydrogen propellant. To address this technology gap, NASA has undergone the development of a high-efficiency, high-capacity 20 K cryocooler via an SBIR partnership with Creare. While the January 2025 testing demonstrated functionality and compliance with contractual requirements, the objective of the NASA-led characterization effort was to generate a dataset for supporting future mission designs across a broader operating envelope, including off-nominal conditions. The cryocooler demonstrated strong performance, achieving a peak coefficient of performance of 17.68% relative to Carnot efficiency and a maximum lift capacity of 24.4 W at 21 K. Overall, the results confirm that the 20 W 20 K cryocooler provides a flexible range of capabilities to enable zero boil-off storage of liquid hydrogen for future Lunar and Mars missions. The data collected provide a strong foundation for model validation and future system design efforts.
Cryocooler↗
Characterization Testing of a 20 Kelvin Cryocooler for Space Applications
Future NASA space exploration missions will require long-duration storage and liquefaction of cryogenic liquids, enabled by active cooling provided by cryocoolers. Recent gap analyses of Lunar and Mars transportation systems have identified 20 K-class cryocoolers as a critical enabling technology for chemical and nuclear thermal propulsion architectures using liquid hydrogen propellant. To address this technology gap, NASA has undergone the development of a high-efficiency, high-capacity 20 K cryocooler via an SBIR partnership with Creare. While the January 2025 testing demonstrated functionality and compliance with contractual requirements, the objective of the NASA-led characterization effort was to generate a dataset for supporting future mission designs across a broader operating envelope, including off-nominal conditions. The cryocooler demonstrated strong performance, achieving a peak coefficient of performance of 17.68% relative to Carnot efficiency and a maximum lift capacity of 24.4 W at 21 K. Overall, the results confirm that the 20 W 20 K cryocooler provides a flexible range of capabilities to enable zero boil-off storage of liquid hydrogen for future Lunar and Mars missions. The data collected provide a strong foundation for model validation and future system design efforts.
Zero boil-off↗
Vaporization Would Cool Primary Battery
Temperature of discharging high-power-density primary battery maintained below specified level by evaporation of suitable liquid from jacket surrounding battery, according to proposal. Pressure-relief valve regulates pressure and boiling temperature of liquid. Less material needed in cooling by vaporization than in cooling by melting. Technique used to cool batteries in situations in which engineering constraints on volume, mass, and location prevent attachment of cooling fins, heat pipes, or like.
Bhandari, Pradeep↗
Dry Pressed, High Areal Loading Electrode Architectures Enabled by Holey Graphene
For future electric aviation, advanced battery cell chemistry beyond lithium ion batteries are required to meet mission requirements. High energy density battery concepts such as lithium-sulfur (Li-S) and lithiumoxygen (Li-O2) chemistries are being intensively investigated to realize their extraordinary theoretical promise in terms of energy density. Most fabrication methods of cathodes for these novel battery chemistries followed a conventional approach. In this approach, the active material is mixed with a polymer binder and a conductive carbon in a high-boiling organic solvent to form a slurry, followed by casting onto a current collector and solvent evaporation. The process is usually lengthy and poses environmental hazards due to the use of organic solvents.
Lin Yi↗
Three-Dimensional, Multi-Phase Unsteady Flow of Water Undergoing Multiple Flow Regimes (Free Surface, Droplet, and Fluid Film)
Upgrades to the Mobile Launcher Rainbird nozzles were implemented for the Artemis II launch to increase water flow rates and improve the mobile launcher deck water coverage while minimizing the water throw distance to avoid impacting the Space Launch System (SLS) rocket nozzles. This analysis describes the use of a three-dimensional, multi-phase, free surface transient computational fluid dynamics (CFD) model to simulate the water flowing out through the Rainbird nozzle and ensure compliance with current Artemis requirements. Flow rates, pressures, velocities, nozzle throw distances and coverage areas were analyzed to ensure that system requirements are met. The completed results were used to make necessary modifications to the Ignition Overpressure Protection and Sound Suppression (IOP/SS) system Rainbird design by reshaping the nozzle used in the system.
CFD↗
Heat Flow and Thermal History of the Moon
An analysis was made of current heat flow data and thermal models of lunar evolution which satisfy the diverse information that has accumulated on internal processes.
Ye A Lyubimova↗
Thermochemically-Closed Sonic-Flow Inversion for Enthalpy and Temperature in Multispecies Arc-Jet Flows
A thermochemically-closed sonic-flow inversion framework (TSIF) is developed to infer bulk enthalpy and total temperature upstream of a choked nozzle in arc-jet flows. The formulation recasts a pressure-rise total enthalpy quantification technique as an inverse problem in characteristic-velocity c * space using measured mass flow rate, upstream total pressure, gas composition, and nozzle throat geometry as inputs. Unlike calorimetric energy-balance approaches or optical diagnostics, the method relies primarily on routinely measured facility quantities combined with explicit thermochemical closure. Thermochemical states are obtained using NASA’s open-source Chemical Equilibrium with Applications (CEA) code, enabling construction of a chemistry-consistent relation between characteristic velocity, total enthalpy, and total temperature under equilibrium or frozen assumptions. A discharge coefficient is self-calibrated using cold-flow (arc-off) operation data and applied to hot-flow (arc-on) measurements, enabling upstream losses to be accounted for without empirical correlations. The framework is applied to air, N 2 , and CO 2 –N 2 arc-jet flows and demonstrates expected trends for the inferred thermochemical states as function of arc power, specific energy input, mass-flow, heater configuration, and test gas. In the air limit, under equilibrium assumptions, the method recovers the classical high-enthalpy asymptotic correlation of Winovich with a mean residual of 4.4%, demonstrating compatibility with established sonic-flow scaling, while extending applicability to arbitrary multi-species mixtures and non-equilibrium chemistry. The framework provides a mixture-flexible methodology for determining bulk thermochemical states in modern arc-jet environments using routine facility pressure, mass-flow, gas-composition, and nozzle-geometry information together with a cold-flow calibration.
stagnation heat flux↗
Thermochemically-Closed Sonic-Flow Inversion for Enthalpy and Temperature in Multispecies Arc-Jet Flows
A thermochemically-closed sonic-flow inversion framework (TSIF) is developed to infer bulk enthalpy and total temperature upstream of a choked nozzle in arc-jet flows. The formulation recasts a pressure-rise total enthalpy quantification technique as an inverse problem in characteristic-velocity c * space using measured mass flow rate, upstream total pressure, gas composition, and nozzle throat geometry as inputs. Unlike calorimetric energy-balance approaches or optical diagnostics, the method relies primarily on routinely measured facility quantities combined with explicit thermochemical closure. Thermochemical states are obtained using NASA’s open-source Chemical Equilibrium with Applications (CEA) code, enabling construction of a chemistry-consistent relation between characteristic velocity, total enthalpy, and total temperature under equilibrium or frozen assumptions. A discharge coefficient is self-calibrated using cold-flow (arc-off) operation data and applied to hot-flow (arc-on) measurements, enabling upstream losses to be accounted for without empirical correlations. The framework is applied to air, N 2 , and CO 2 –N 2 arc-jet flows and demonstrates expected trends for the inferred thermochemical states as function of arc power, specific energy input, mass-flow, heater configuration, and test gas. In the air limit, under equilibrium assumptions, the method recovers the classical high-enthalpy asymptotic correlation of Winovich with a mean residual of 4.4%, demonstrating compatibility with established sonic-flow scaling, while extending applicability to arbitrary multi-species mixtures and non-equilibrium chemistry. The framework provides a mixture-flexible methodology for determining bulk thermochemical states in modern arc-jet environments using routine facility pressure, mass-flow, gas-composition, and nozzle-geometry information together with a cold-flow calibration.
inviscid theory↗
Mitigating Grazing Flow Impedance Eduction Errors in Additively Manufactured Porous Liners
Explore the source record for details and available documents.
liners↗
Turbulance Boundary Conditions for Shear Flow Analysis, Using the DTNS Flow Solver
The effects of different turbulence boundary conditions were examined for two classical flows: a turbulent plane free shear layer and a flat plate turbulent boundary layer with zero pressure gradient. The flow solver used was DTNS, an incompressible Reynolds averaged Navier-Stokes solver with k-epsilon turbulence modeling, developed at the U.S. Navy David Taylor Research Center. Six different combinations of turbulence boundary conditions at the inflow boundary were investigated: In case 1, 'exact' k and epsilon profiles were used; in case 2, the 'exact' k profile was used, and epsilon was extrapolated upstream; in case 3, both k and epsilon were extrapolated; in case 4, the turbulence intensity (I) was 1 percent, and the turbulent viscosity (mu(sub t)) was equal to the laminar viscosity; in case 5, the 'exact' k profile was used and mu(sub t) was equal to the laminar viscosity; in case 6, the I was 1 percent, and epsilon was extrapolated. Comparisons were made with experimental data, direct numerical simulation results, or theoretical predictions as applicable. Results obtained with DTNS showed that turbulence boundary conditions can have significant impacts on the solutions, especially for the free shear layer.
M Mizukami↗