Search NASA⌕ Search

SEARCH · Search NASA

Results for “Two-phase heat transfer coefficient”

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

Porous mesh manifold for enhanced boiling performance

High-performance electronics are continuously demanding cooling of higher heat fluxes. Phase-change cooling, including pool boiling, is a useful approach to address this challenge; however, competition between liquid and vapor flows generally limit the heat fluxes that can be dissipated. A range of strategies to control these flows have been investigated previously, including capillary guides. Here a manifold structure formed from a metallic mesh is investigated to control the disposition of liquid and vapor phases above a pool fed boiling surface enhanced with porous structures. Copper mesh forms defined liquid flow paths, using capillary action to guide and distribute liquid evenly over the heated surface, along with open channels to facilitate vapor escape. The mesh provides a novel structure for liquid guidance that imposes low resistance to liquid flow while occluding a minimal area of heated surface underneath. The manifold performance is characterized in boiling fed by a pool of water above a laser-textured aluminum nitride heat dissipation surface with pin–fin structures having heights of 110 µm and spacing of 30 µm with a heated area of 5 mm x 5 mm. A maximum heat flux of 490 W/cm 2 is reached with the manifold in the pool fed configuration, representing an increase of more than 65% over the porous pin fin surface alone. The maximum stable superheat observed for the manifold of 36K is 14K higher than that for the porous surface without the manifold. The factors limiting performance of the manifold are analyzed. High superheat is attributed to partial flooding of the boiling surface as suggested by the reduction in superheat using external suction. Similar systems and structures for enhanced two-phase cooling are compared.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1987

Test results of a pumped two-phase mounting plate with ammonia

The design, fabrication, and testing of full-scale prototype units of a two-phase mounting plate (TPMP), which will be used in a two-phase ammonia-based thermal control system for a large spacecraft, are described. The mounting plate uses an evaporator design in which liquid is mechanically pumped through porous feed tubes within the plate. The prototype TPMPs were tested with ammonia at heat loads over 3000 W (3.2 W/sq cm) and local heat fluxes of up to 4 W/sq cm. Calculated total heat transfer coefficients from these tests were between 0.8 and 1.0 W/sq cm per C. This represents a better than twenty-fold improvement over comparable single-phase heat transfer coefficients. Design diagrams are included.

Swanson, T. D.↗
NASA NTRS2008

Spectral Attenuation of Sound in Dilute Suspensions with Nonlinear Particle Relaxation

Theoretical studies on the dissipation and dispersion of sound in two-phase suspensions have been briefly reviewed. Previous studies on the sound attenuation in particle-laden flows under Stokesian drag and conduction-controlled heat transfer have been extended to accommodate the nonlinear drag and heat transfer. It has been shown that for large particle-to-fluid density ratio, the particle Reynolds number bears a cubic relationship with Omega Tau(sub d) (where Omega is the circular frequency and Tau(sub d) the Stokesian particle relaxation time). This dependence leads to the existence of a peak value in the linear absorption coefficient occurring at a finite value Omega Tau (sub d). Comparison of the predictions with the test data for the spectral attenuation of sound with water injection in a perfectly expanded supersonic air jet shows a satisfactory trend of the theory accounting for nonlinear particle relaxation processes.

Kandula, M.↗
NASA NTRS2007

Transient Numerical Modeling of Catalytic Channels

This paper presents a transient model of catalytic combustion suitable for isolated channels and monolith reactors. The model is a lumped two-phase (gas and solid) model where the gas phase is quasi-steady relative to the transient solid. Axial diffusion is neglected in the gas phase; lateral diffusion, however, is accounted for using transfer coefficients. The solid phase includes axial heat conduction and external heat loss due to convection and radiation. The combustion process utilizes detailed gas and surface reaction models. The gas-phase model becomes a system of stiff ordinary differential equations while the solid phase reduces, after discretization, into a system of stiff ordinary differential-algebraic equations. The time evolution of the system came from alternating integrations of the quasi-steady gas and transient solid. This work outlines the numerical model and presents some sensitivity studies on important parameters including internal transfer coefficients, catalytic surface site density, and external heat-loss (if applicable). The model is compared to two experiments using CO fuel: (1) steady-state conversion through an isothermal platinum (Pt) tube and (2) transient propagation of a catalytic reaction inside a small Pt tube. The model requires internal mass-transfer resistance to match the experiments at lower residence times. Under mass-transport limited conditions, the model reasonably predicted exit conversion using global mass-transfer coefficients. Near light-off, the model results did not match the experiment precisely even after adjustment of mass-transfer coefficients. Agreement improved for the first case after adjusting the surface kinetics such that the net rate of CO adsorption increased compared to O2. The CO / O2 surface mechanism came from a sub-set of reactions in a popular CH4 / O2 mechanism. For the second case, predictions improved for lean conditions with increased external heat loss or adjustment of the kinetics as in the first case. Finally, the results show that different initial surface-species distribution leads to different steady-states under certain conditions. These results demonstrate the utility of a lumped two-phase model of a transient catalytic combustor with detailed chemistry.

Struk, Peter M.↗
NASA NTRS1992

A prototype heat pipe heat exchanger for the capillary pumped loop flight experiment

A Capillary Pumped Two-Phase Heat Transport Loop (CAPL) Flight Experiment, currently planned for 1993, will provide microgravity verification of the prototype capillary pumped loop (CPL) thermal control system for EOS. CAPL employs a heat pipe heat exchanger (HPHX) to couple the condenser section of the CPL to the radiator assembly. A prototype HPHX consisting of a heat exchanger (HX), a header heat pipe (HHP), a spreader heat pipe (SHP), and a flow regulator has been designed and tested. The HX transmits heat from the CPL condenser to the HHP, while the HHP and SHP transport heat to the radiator assembly. The flow regulator controls flow distribution among multiple parallel HPHX's. Test results indicated that the prototype HPHX could transport up to 800 watts with an overall heat transfer coefficient of more than 6000 watts/sq m-deg C. Flow regulation among parallel HPHX's was also demonstrated.

Ku, Jentung↗
NASA NTRS2012

Measurement of Two-Phase Flow and Heat Transfer Parameters using Infrared Thermometry

A novel technique to measure heat transfer and liquid film thickness distributions over relatively large areas for two-phase flow and heat transfer phenomena using infrared (IR)thermometry is described. IR thermometry is an established technology that can be used to measure temperatures when optical access to the surface is available in the wavelengths of interest. In this work, a midwave IR camera (3.6-5.1 microns) is used to determine the temperature distribution within a multilayer consisting of a silicon substrate coated with a thin insulator. Since silicon is largely transparent to IR radiation, the temperature of the inner and outer walls of the multilayer can be measured by coating selected areas with a thin, IR opaque film. If the fluid used is also partially transparent to IR, the flow can be visualized and the liquid film thickness can be measured. The theoretical basis for the technique is given along with a description of the test apparatus and data reduction procedure. The technique is demonstrated by determining the heat transfer coefficient distributions produced by droplet evaporation and flow boiling heat transfer.

Kim, Tae-Hoon↗
DOE OSTI2025

MOSCATO Development and Integration in Fiscal Year 2025: Implementation of Multiphase, Multiphysics Modeling Capabilities for Molten Salt Systems

MOSCATO (Molten Salt Chemistry and Transport) is a multiphysics code that provides high-fidelity, coupled simulations of fluid flow, heat transfer, mass transfer, chemistry, electrochemical phenomena, and alloy corrosion for molten salt systems. In FY25, significant developments were made to the code package, enhancing its capabilities for modeling all relevant phenomena within operating moltens salt reactors (MSRs). The developments and activities in FY25 included: 1. Implementation of Level-Set methods to enable modeling of single-bubble behavior in molten salts. In FY25, the Level-Set two-phase flow modeling implementation was improved to simulate single bubble behavior with molten salt media. The large density and viscosity ratios between typical gases and molten salt liquids present challenges for these types of numerical solvers. With enhancements to the pressure projection method, MOSCATO’s Level-Set solver was able to be successfully validated to experiments related to helium bubble rise in stagnant molten salt. The simulated bubble rising velocity showed reasonable good agreement with experimental measurements. The bubble shape and dynamics were also visually compared with experimental snapshots, demonstrating a good qualitative match. 2. Generation of mass transfer correlations for multiphase flow systems. To enable calculations of the tritium transport across the interface between gas bubbles and salt, we modeled high- Schmidt-number mass transfer around a sphere across a broad range of Reynolds numbers. The mesh near the sphere surface was highly refined to resolve steep concentration gradients caused by the low diffusion coefficient. Literature-based mass transfer correlations were compared with the numerical results, and modifications were proposed to improve agreement, particularly at higher Schmidt numbers. These mass transfer correlations were subsequently provided to other national laboratories to help enable high quality mass transfer simulations using lower-order solvers under development within the NEAMS program. 3. Preliminary implementation of a bubbly flow solver. To model bubbly flow in molten salt, we implemented a bubbly flow solver for void fractions less than 5%. To do so, an algebraic relative velocity model that assumes small bubbles with rapid momentum equilibration was added to MOSCATO to compute bubble velocities. Preliminary comparisons with experimental data showed reasonable agreement, and further improvements are underway. 4. Generation of mass transfer correlations for MSRE subchannel The Molten-Salt Reactor Experiment (MSRE) was a landmark historical project that demonstrated the feasibility of molten-salt reactor technology. The MSRE campaign also generated a significant body of experimental data and reports that continue to support molten-salt–related research. In this report, the MSRE core subchannel was used as the reference geometry for a mass transfer study performed with MOSCATO. The geometry and computational mesh were adapted from a previous study, providing adequate resolution for the relatively low Reynolds number in this case. Additional mesh refinement was applied to reach higher Schmidt numbers, enabling the derivation of a reliable mass-transfer correlation for the present scenario. 5. Simulations of oxygen ingressions into molten salt. In the previous fiscal year, we initiated a study to simulate oxygen ingression in stagnant salt. As oxygen enters the salt through its surface, it reacts with Ce 3+ to form solid CeO 2 and other reaction products. To more fully capture the complex diffusion-convection-reaction mechanisms, capabilities for modeling natural convection in the salt vessel were added. These were needed as the flow of the ingressed gas induced flow in the salt caused by surface shear and non-isothermal effects. With these updated physics in place, we were able to successfully reproduce the experimental results for the rate of change of CeCl 3 concentrations versus time. 6. Flow corrosion model validation. In FY24, MOSCATO’s corrosion model was validated against static corrosion experiments. In FY25, this work was extended to a flow corrosion experiment, where FLiNaK salt was driven by natural convection, with initial salt impurities to initiate corrosion. Despite uncertainties in parameters such as elemental diffusion coefficients in the alloy and unknown H + concentrations, the simulations achieved good agreement with experimental results, especially in predicting sample mass losses.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Non-dimensional performance and safety parameters for heat pipes

The use of heat pipes in safety-critical systems such as nuclear microreactors dictates the development of generalized, practical, scalable performance and safety parameters. Traditional dimensional metrics, while informative, lack the universality required for comparative analysis across varying designs and operating regimes. Here, this work introduces a comprehensive set of non-dimensional parameters to characterize heat pipe performance and safety, including capillary performance, effective thermal conductivity, response time, exergetic efficiency, allowable temperature gradients, allowable rate of temperature change, priming coefficients, and factor of safety. A reference heat pipe design representative of microreactor applications was analyzed via the developed parameters using both traditional analytical models and Sockeye simulations under transient and steady-state conditions. Sodium, potassium, and water were evaluated as working fluids to demonstrate the applicability of the framework across a broad temperature range. The proposed non-dimensional parameters effectively captured key thermal-hydraulic behaviors and safety concerns, as was demonstrated via Sockeye simulations. This framework supports the development of design optimization strategies, operational protocols, and safety assurance practices for advanced reactor systems and other high-reliability applications.

42 - ENGINEERING↗
DOE OSTI2025

Integrated Molten Salt Reactor Modeling Capabilities in NEAMS Thermal Hydraulics Tools

The DOE neams program supports a full range of computational thermal fluids analysis capabilities and code developments for a broad range of advanced reactor concepts. The research and development approach under the thermal fluids technical area synergistically combines three length and time scales in a hierarchical multi-scale approach. To enable multi-scale thermal fluids capability using these codes, a key joint effort has been underway to develop an integrated system- and engineering-scale thermal fluids analysis capability, through integration of SAM and Pronghorn codes, both based on the MOOSE framework. This report summarizes recent advances in developing an integrated system- and engineering-scale modeling capability for the msr concept, which has gained significant interest in recent years. A consistent framework was established by coupling Pronghorn and SAM through the Saline interface, with thermophysical properties provided by the Molten Salt Thermal Property Database (MSTDB-TP). Further improvements were made to the coupling schemes and domain-overlapping strategies, enhancing the stability and robustness of multi-code simulations. Verification and validation efforts demonstrate the accuracy of this integration across a range of benchmark problems, including one-dimensional heated pipe flows, three-dimensional natural convection loops with evolving isotopic compositions, and \gls{msre} demonstration cases. Within Pronghorn, new capabilities were introduced to model corrosion and noble-metal plating phenomena, supported by an extended thermal-hydraulics framework and refined turbulence treatments. To capture two-phase flow behavior, a multiphase Euler–Euler model was implemented in Pronghorn, including advanced closure relations, high-resolution advection techniques, and capillary force reconstruction. Preliminary verification cases confirm the fidelity of the approach, while planned validation efforts target canonical multiphase benchmarks and application to msr components such as the msre pump bowl. Finally, updates to SAM’s msr mass transfer modeling were extended to consider noble gas migration into porous structures like graphite. The point kinetics model was updated to include reactivity feedback contributions from any defined species, such as xenon. The gas transport model was expanded for applicability to gas mixtures, bubble efflux phenomena, and species transport between liquid and gas phases. A selection of multi-scale Sherwood number correlations from MOSCATO/NekRS and multi-phase correlations from literature have been added for improved accuracy in calculating mass transfer coefficients. A companion effort on developing system-level redox corrosion has also been incorporated into SAM. Collectively, these enhancements strengthen the predictive capability of SAM and Pronghorn for simulating MSR thermal-hydraulics, corrosion, multiphase behavior, and fission-product transport, providing a more complete toolset for design, safety analysis, and licensing support of next-generation \gls{msr}s.

42 - ENGINEERING↗
DOE OSTI2026

Deflagration to Detonation Transition Update: XDDT Code Modularization

A legacy FORTRAN 77 implementation of the Baer–Nunziato two-phase mixture theory for deflagration-to-detonation transition (DDT) in reactive granular materials—hereafter the XDDT (eXplosive DDT) code—has been modularized to Fortran 90 with modular structure, external input files, and adaptive mesh capability. During validation, two code defects were identified and corrected: an inconsistency in the nodal solid pressure evaluation and a nonphysical burn-front tracking criterion. The ignition criterion was also corrected to use the granular surface temperature from the interface heat transfer model, matching the original Baer implementation. An initial attempt to validate against Figure 3 of the original Baer and Nunziato (1986) paper revealed that the code’s detonation velocity on a 201-node mesh (5.5 km/s) was approximately 21% below the expected Chapman–Jouguet value for 70% TMD HMX (∼7 km/s). Validation was redirected to the piston-driven DDT experiments of McAfee et al. (1989), Shot B-9036, for which well-characterized ionization-pin data are available. With the compaction-burn coefficient calibrated to 𝐶 𝛼 = 75, the XDDT code reproduces the DDT transition time to within 0.4% and produces a steady-state detonation velocity within 4% of the McAfee experimental value of 6.36 km/s. The burn model was generalized to support pressure-dependent exponents, enabling application to nitrocellulose-based ball propellants (TS3659) with a cube-root pressure dependence. Validation against the Sandusky/Baer PDC82 piston-impact experiment yielded a reactive wave velocity of 2.3–2.8 km/s, in good agreement with the experimental value of ∼2.2 km/s, and wave coalescence within 5% of the experimental timing. The mathematical model, input parameter requirements, and a roadmap for extending XDDT to PETN with an autocatalytic burn model are presented.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2025

A Finite Element Method for Compressible and Turbulent Multiphase Flow Instabilities with Heat Transfer

We present a new finite element framework for modeling compressible, turbulent multiphase flows with heat transfer. For two-fluid systems with a free surface, the Volume of Fluid (VOF) method is implemented without the need for interface reconstruction, while turbulence is resolved using a dynamic Vreman large eddy simulation (LES) model. Unlike most two-phase VOF studies, which neglect heat transfer, the present approach incorporates energy transport equations within the VOF formulation to account for heat exchange, an effect particularly important in turbulent flows. Conjugate heat transfer is often challenging in finite volume methods, which require explicit specification of heat fluxes at the solid–fluid interface, limiting accuracy and predictive capability. By contrast, the finite element formulation does not require heat flux inputs, allowing more accurate and robust simulation of heat transfer between solids and fluids. The method is demonstrated through three representative cases. First, a two-fluid instability with a single-mode perturbation is simulated and validated against analytical growth rates. Second, conjugate heat transfer is examined in a high-temperature flow over a cold metal cylinder, with validation performed both quantitatively—via pressure coefficient comparisons with experimental data—and qualitatively using vector field topology. Finally, compressible spray injection and breakup are modeled, demonstrating the ability of the framework to capture interfacial dynamics and atomization under turbulent, high-speed conditions. In the compressible spray injection and breakup case, the results indicate that the finite element formulation achieved higher predictive accuracy and robustness than the finite-volume method. With the same mesh resolution, the FEM reduced the root mean square error (RMSE) and mean absolute percentage error (MAPE) from 6.96 mm and 26.0% (for the FVM) to 4.85 mm and 12.7%, respectively, demonstrating improved accuracy and robustness in capturing interfacial dynamics and heat transfer. The study also introduced vector field topology to visualize and interpret coherent flow structures and instabilities, offering insights beyond conventional scalar-field analyses.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1977

Assessment of analytical techniques for predicting solid propellant exhaust plumes

The calculation of solid propellant exhaust plume flow fields is addressed. Two major areas covered are: (1) the applicability of empirical data currently available to define particle drag coefficients, heat transfer coefficients, mean particle size and particle size distributions, and (2) thermochemical modeling of the gaseous phase of the flow field. Comparisons of experimentally measured and analytically predicted data are made. The experimental data were obtained for subscale solid propellant motors with aluminum loadings of 2, 10 and 15%. Analytical predictions were made using a fully coupled two-phase numerical solution. Data comparisons will be presented for radial distributions at plume axial stations of 5, 12, 16 and 20 diameters.

Tevepaugh, J. A.↗
NASA NTRS1985

A swirl flow evaporative cold plate

A forced flow evaporative cold plate is under development for future application to the thermal bus concept being pursued by NASA for Space Station Thermal Control. The vaporizer is a swirl-flow device employing a spiral tube coil geometry sandwiched between conductive metal plates upon which electric components could be mounted. This concept is based on the inherent phase separation that occurs in a two phase stream in curvilinear flow. This is a zero 'g' design with one 'g' all-attitude capability and is capable of high heat transfer coefficients, good isothermality, and the ability to function at heat fluxes approaching 5w/sq cm on the cold plates (10w/sq cm on the tube wall) with Freon 114. The advantages of this design over other two phase evaporator approaches are high heat flux capability, simplified control requirements, insensitivity to micro-gravity oscillations, and inexpensive manufacturability. The program included design, fabrication, and test of such a cold plate utilizing an existing test stand developed for two-phase thermal management system (TPTMS) testing. Test results analysis and conclusions are included.

Niggemann, R. E.↗
NASA NTRS1990

High heat flux single phase heat exchanger

This paper presents the results obtained to date in a program to develop a high heat flux, single-phase heat exchanger for spacecraft thermal management. The intended application is a net generation interface heat exchanger to couple the crew module water thermal bus to the two-phase ammonia main thermal bus in the Space Station Freedom. The large size of the interface heat exchanger is dictated by the relatively poor water-side heat transfer characteristics. The objective of this program is to develop a single-phase heat transfer approach which can achieve heat fluxes and heat transfer coefficients comparable to those of the evaporation ammonia side. A new heat exchanger concept has been developed to meet these objecties. The main feature of this heat exchanger is that it can achieve very high heat fluxes with a pressure drop one to two orders of magnitude lower than those of previous microchannel or jet impingement high heat flux heat exchangers. This paper describes proof-of-concept experiments performed in air and water and presents analytical model of the heat exchanger.

Valenzuela, Javier A.↗
NASA NTRS2004

Gas-Liquid Flows and Phase Separation

Common issues for space system designers include:Ability to Verify Performance in Normal Gravity prior to Deployment; System Stability; Phase Accumulation & Shedding; Phase Separation; Flow Distribution through Tees & Manifolds Boiling Crisis; Heat Transfer Coefficient; and Pressure Drop.The report concludes:Guidance similar to "A design that operates in a single phase is less complex than a design that has two-phase flow" is not always true considering the amount of effort spent on pressurizing, subcooling and phase separators to ensure single phase operation. While there is still much to learn about two-phase flow in reduced gravity, we have a good start. Focus now needs to be directed more towards system level problems .

McQuillen, John↗
NASA NTRS2006

Miniature Distillation Column for Producing LOX From Air

The figure shows components of a distillation column intended for use as part of a system that produces high-purity liquid oxygen (LOX) from air by distillation. (The column could be easily modified to produce high-purity liquid nitrogen.) Whereas typical industrial distillation columns for producing high-purity liquid oxygen and/or nitrogen are hundreds of feet tall, this distillation column is less than 3 ft (less than about 0.9 m) tall. This column was developed to trickle-charge a LOX-based emergency oxygen system (EOS) for a large commercial aircraft. A description of the industrial production of liquid oxygen and liquid nitrogen by distillation is prerequisite to a meaningful description of the present miniaturized distillation column. Typically, such industrial production takes place in a chemical processing plant in which large quantities of high-pressure air are expanded in a turboexpander to (1) recover a portion of the electrical power required to compress the air and (2) partially liquefy the air. The resulting two-phase flow of air is sent to the middle of a distillation column. The liquid phase is oxygen-rich, and its oxygen purity increases as it flows down the column. The vapor phase is nitrogen-rich and its nitrogen purity increases as it flows up the column. A heater or heat exchanger, commonly denoted a reboiler, is at the bottom of the column. The reboiler is so named because its role is to reboil some of the liquid oxygen collected at the bottom of the column to provide a flow of oxygen-rich vapor. As the oxygen-rich vapor flows up the column, it absorbs the nitrogen in the down-flowing liquid by mass transfer. Once the vapor leaves the lower portion of the column, it interacts with down-flowing nitrogen liquid that has been condensed in a heat exchanger, commonly denoted a condenser, at the top of the column. Liquid oxygen and liquid nitrogen products are obtained by draining some of the purified product at the bottom and top of the column, respectively. Because distillation is a mass-transfer process, the purity of the product(s) can be increased by increasing the effectiveness of the mass-transfer process (increasing the mass-transfer coefficient) and/or by increasing the available surface area for mass transfer through increased column height. The diameter of a distillation column is fixed by pressure-drop and mass-flow requirements. The approach taken in designing the present distillation column to be short yet capable of yielding a product of acceptably high purity was to pay careful attention to design details that affect mass-transfer processes.

Rozzi, Jay C.↗
NASA NTRS2011

Analytical Approach in DeCoM

There are many papers on describing a LHP as an overall system, but few detail on the condenser section of a loop heat pipe. The DeCoM (Deepak Condenser Model) method utilizes user set initial parameters in-order to simulate a condenser by calculating the interactions between the fluid and the wall. Equations are derived for two sections of the condenser: a two-phase section and a subcooled (liquid) section. All Equations are based upon the conservation of energy theory, from which fluid temperature, and fluid quality values are solved. In order to solve for the heat transfer value, between fluid and the wall in two phase section, the Lockhart-Martinelli correlation method was implemented as a solution approach. For Liquid phase, the Reynolds number was used in-order to differentiate the flow state, from either turbulent or laminar, and Nusselt number was used to solve for the film coefficient. To represent these calculations for both sections a flow chart is presented in order to display the execution process of DeCoM. The benefit of DeCoM is that it is capable of performing preliminary analysis without requiring a license and without much of users knowledge on condensers.

Patel, Deepak↗
NASA NTRS2019

Two-Phase Convection Heat Transfer Correlations for Liquid Hydrogen Pipe Chilldown

Recently, heat transfer correlations based on liquid nitrogen (LN2) and liquid hydrogen (LH2) pipe quenching data were developed to improve the predictive accuracy of lumped node codes like SINDA/FLUINT and the Generalized Fluid System Simulation Program (GFSSP). After implementing these correlations into both programs, updated model runs showed strong improvement in LN2 pipe chilldown modeling but only modest improvement in LH2 modeling. Due to large differences in thermal and fluid properties between the two fluids, results indicated a need to develop a separate set of LH2-only correlations to improve the accuracy of the simulations. This paper presents a new set of two-phase convection heat transfer correlations based on LH2 pipe quenching data. A correlation to predict the bulk vapor temperature was developed after analysis showed that high amounts of thermal nonequilibrium of the liquid and vapor phases occurred during film boiling of LH2. Implemented in a numerical model, the new correlations achieve a mean absolute error of 19.5 K in the predicted wall temperature when compared to recent LH2 pipe chilldown data, an improvement of 40% over recent GFSSP predictions. This correlation set can be implemented in simulations of the transient LH2 chilldown process. Such simulations are useful for predicting the chilldown time and boil-off mass of LH2 for applications such as the transfer of LH2 from a ground storage tank to the rocket vehicle propellant tank, or through a rocket engine feedline during engine startup.

critical heat flux↗