Film boiling of cryogenic hydrogen during upward and downward flow
Buoyancy effect on cryogenic hydrogen film boiling during upward and downward flows
SEARCH · Search NASA
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.
Buoyancy effect on cryogenic hydrogen film boiling during upward and downward flows
Buoyancy effects on the critical heat flux and general data trends for a liquid nitrogen internal flow system were determined by comparison of upflow and downflow data under identical test conditions. The test section had a 1.28 cm diameter flow passage and a 30.5 cm heated length which was subjected to uniform heat fluxes through resistance heating. Test conditions covered a range of pressures from 3.4 to 10.2 atm, inlet velocities from 0.23 to 3.51 m/sec, with the liquid nitrogen temperature at saturated inlet conditions. Data comparisons showed that the critical heat flux for downflow could be up to 36 percent lower than for upflow. A nonmonotonic relationship between the critical heat flux and velocity was determined for upflow but not for downflow. A limiting inlet velocity of 4.12 m/sec was determined to be the minimum velocity required to completely suppress the influence of buoyancy on the critical heat flux for this saturated inlet flow system. A correlation of this limiting fluid velocity is presented that was developed from previously published subcooled liquid nitrogen data and the saturated data of this investigation.
A meso-scale model developed at MIT [2][11] for boiling processes in water was adapted for cryogens and demonstrated for chilldown of propellant transfer lines in both liquid nitrogen and hydrogen. The sub-grid boiling model accurately captures the contributions to heat transfer from the generation of bubble nuclei, growth, and interaction of the bubbles in the microlayer as well as quenching of the boiling surface following bubble departure. It was adapted for cryogenic fluids using thermodynamic scaling concepts taking into account non-dimensional pressures and temperatures that are scaled by the corresponding critical values for the fluid. The boiling model was demonstrated for line chilldown in liquid nitrogen (Darr et al. [6]). The predicted wall temperature at which quenching occurs was close to the test data while the slope of the temperature curve after quenching is initiated, showing a steeper variation than the test data. Simulations were also performed for liquid hydrogen by simulating experiments of Hartwig et al. [7]. Chilldown times in liquid hydrogen are much more rapid due to higher heat transfer in the film boiling regime and accounting for the higher turbulence levels were found to be important. Furthermore, at the much lower fluid temperatures in liquid hydrogen flows, accounting for the variable thermal properties of the solid material is critical and has a dramatic impact on the quench times. The need for additional studies to better understand the evolution of the quench front in liquid hydrogen was noted.
A meso-scale model developed at MIT [2][11] for boiling processes in water was adapted for cryogens and demonstrated for chilldown of propellant transfer lines in both liquid nitrogen and hydrogen. The sub-grid boiling model accurately captures the contributions to heat transfer from the generation of bubble nuclei, growth, and interaction of the bubbles in the microlayer as well as quenching of the boiling surface following bubble departure. It was adapted for cryogenic fluids using thermodynamic scaling concepts taking into account non-dimensional pressures and temperatures that are scaled by the corresponding critical values for the fluid. The boiling model was demonstrated for line chilldown in liquid nitrogen (Darr et al. [6]). The predicted wall temperature at which quenching occurs was close to the test data while the slope of the temperature curve after quenching is initiated, showing a steeper variation than the test data. Simulations were also performed for liquid hydrogen by simulating experiments of Hartwig et al. [7]. Chilldown times in liquid hydrogen are much more rapid due to higher heat transfer in the film boiling regime and accounting for the higher turbulence levels were found to be important. Furthermore, at the much lower fluid temperatures in liquid hydrogen flows, accounting for the variable thermal properties of the solid material is critical and has a dramatic impact on the quench times. The need for additional studies to better understand the evolution of the quench front in liquid hydrogen was noted.
Forced convection and subcooled boiling heat transfer data for liquid nitrogen and liquid neon were obtained in support of a design study for a 30 tesla cryomagnet cooled by forced convection of liquid neon. This design precludes nucleate boiling in the flow channels as they are too small to handle vapor flow. Consequently, it was necessary to determine boiling incipience under the operating conditions of the magnet system. The cryogen data obtained over a range of system pressures, fluid flow rates, and applied heat fluxes were used to develop correlations for predicting boiling incipience and convective boiling heat transfer coefficients in uniformly heated flow channels. The accuracy of the correlating equations was then evaluated. A technique was also developed to calculate the position of boiling incipience in a uniformly heated flow channel. Comparisons made with the experimental data showed a prediction accuracy of plus or minus 15 percent
Laminar two-phase boundary layer flow in film boiling, obtaining asymptotic solutions for low and high subcooled liquids
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.
One dimensional flow model to correlate subcooled boiling pressure drop in circular tubes with constant heat flux
Tests of cryogenic chilldown in a vertical stainless steel tube are simulated using computational fluid dynamics (CFD). Subcooled liquid nitrogen enters the room-temperature tube at one end, resulting in film boiling and eventual quenching of the tube walls. A CFD-based modeling framework employing the mixture multiphase model has been developed to simulate this chilldown process and applied to validate it against experimental data. Extensive sensitivity studies have been performed to create a consistent modeling framework, and a two test conditions are simulated to test the modeling capability. The importance of the rewetting temperature and its relation to the foam limit temperature is discussed and is considered in the CFD model. The presented results demonstrate a qualitatively consistent representation of the film boiling and quenching regime. However, the quantitative accuracy of the simulations varies with the test conditions. Simulated vapor volume fraction distributions are also presented, although the current lack of experimental volume fraction data does not permit a comparison.
Tests of cryogenic chilldown in a vertical stainless-steel tube are simulated using computational fluid dynamics (CFD). Subcooled liquid nitrogen enters a room-temperature tube at one end, resulting in film boiling and eventual quenching of the tube walls. A CFDbased modeling framework employing the mixture multiphase model was developed to simulate this chilldown process and was applied to validate it against experimental data. Extensive sensitivity studies have been performed to create a consistent modeling framework, and a range of test conditions are simulated to assess the modeling capability. The importance of the rewetting temperature and its relation to the foam limit temperature is discussed and considered in the CFD framework. The presented results demonstrate a qualitatively consistent representation of the film boiling and quenching regime. However, the quantitative accuracy of the simulations varies with the test conditions. In addition, simulated distributions of the vapor volume fraction are presented, but the current lack of experimental volume fraction data prevents comparisons.
Tests of cryogenic chilldown in a vertical stainless-steel tube are simulated using computational fluid dynamics (CFD). Subcooled liquid nitrogen enters a room-temperature tube at one end, resulting in film boiling and eventual quenching of the tube walls. A CFDbased modeling framework employing the mixture multiphase model was developed to simulate this chilldown process and was applied to validate it against experimental data. Extensive sensitivity studies have been performed to create a consistent modeling framework, and a range of test conditions are simulated to assess the modeling capability. The importance of the rewetting temperature and its relation to the foam limit temperature is discussed and considered in the CFD framework. The presented results demonstrate a qualitatively consistent representation of the film boiling and quenching regime. However, the quantitative accuracy of the simulations varies with the test conditions. In addition, simulated distributions of the vapor volume fraction are presented, but the current lack of experimental volume fraction data prevents comparisons.
Closed loop low pressure forced flow convective boiling system analyzed for instabilities due to feed system/boiler interaction
Explore the source record for details and available documents.
A multi-purpose thermal insulation test apparatus is used for testing insulation materials, or other components. The test apparatus is a fluid boil-off calorimeter system for calibrated measurement of the apparent thermal conductivity (k-value) of a specimen material at a fixed vacuum level. The apparatus includes an inner vessel for receiving a fluid with a normal boiling point below ambient temperature, such as liquid nitrogen, enclosed within a vacuum chamber. A cold mass assembly, including the inner vessel and thermal guards, is suspended from the top of the vacuum chamber. Handling tools attach to the cold mass assembly for convenient manipulation of the assembly and for the installation or wrapping of insulation test materials. Liquid nitrogen is typically supplied to the inner vessel using a fill tube with funnel. A single port through the top of the vacuum chamber facilitates both filling and venting. Aerogel composite stacks with reflective films are fastened to the top and the bottom of the inner vessel as thermal guards. The comparative k-value of the insulation material is determined by measuring the boil-off flow rate of gas, the temperature differential across the insulation thickness, and the dimensions (length and diameters) of the test specimen.
Geomorphic evidence suggests that recent gullies on Mars were formed by fluvial activity. The Martian gully features are significant because their existence implies the presence of liquid water near the surface on Mars in geologically recent times. Irrespective of the ultimate source of the fluid carving the gullies, we seek to understand the behavior of this fluid after it reaches the Martian surface. We find that, contrary to popular belief, the fluvially-carved Martian gullies require formation conditions such as now occur on Mars, outside of the temperature-pressure stability regime of liquid water. Mars Global Surveyor observations of gully length and our modeling of water stability are consistent with gully formation from the action of pure liquid water that is simultaneously boiling and freezing.
A test apparatus and method of its use for evaluating various performance aspects of a test specimen is disclosed. A chamber within a housing contains a cold mass tank with a contact surface in contact with a first surface of a test specimen. The first surface of the test specimen is spaced from the second surface of the test specimen by a thickness. The second surface of the test specimen is maintained at a a constant temperature by a liquid disposed within the cold mass tank. A boil-off flow rate of the gas is monitored and provided to a processor along with the temperature of the first and second surfaces of the test specimen. The processor calculates thermal insulation values of the test specimen including comparative values for heat flux and apparent thermal conductivity k-value). The test specimen may be placed in any vacuum pressure level ranging from about 0.01 millitorr to 1,000,000 millitorr with different residual gases as desired. The test specimen may be placed under a mechanical load with the cold mass tank and another factors may be imposed upon the test specimen so as to simulate the actual use conditions.
A test apparatus and method of its use for evaluating various performance aspects of a test specimen is disclosed. A chamber within a housing contains a cold mass tank with a contact surface in contact with a first surface of a test specimen. The first surface of the test specimen is spaced from the second surface of the test specimen by a thickness. The second surface of the test specimen is maintained at a desired warm temperature. The first surface is maintained at a constant temperature by a liquid disposed within the cold mass tank. A boil-off flow rate of the gas is monitored and provided to a processor along with the temperature of the first and second surfaces of the test specimen. The processor calculates thermal insulation values of the test specimen including comparative values for heat flux and apparent thermal conductivity (k-value). The test specimen may be placed in any vacuum pressure level ranging from about 0.01 millitorr to 1,000,000 millitorr with different residual gases as desired. The test specimen may be placed under a mechanical load with the cold mass tank and another factors may be imposed upon the test specimen so as to simulate the actual use conditions.
Two-phase thermal management systems, with both boiling and condensation processes, offer great potential and heat transfer coefficients that are orders of magnitude higher than traditional single-phase systems. However, two-phase flows can suffer from a wide range of interfacial instabilities leading to significant thermal performance degradation. In this study, we aim to detect regime transitions and characterize dominating physical mechanisms of flow condensation, such as turbulent diffusion in annular liquid film and interfacial waves, using an integrated system of acoustic, modal, and optical sensing techniques and thermofluidic characterizations. A wideband acoustic emissions sensor and high-sensitivity accelerometer are utilized to capture acoustic and vibrational signatures that signal the onset of liquid film formation and interfacial waves during flow pattern transitions. Compared to optical imaging, wideband acoustic emission sensing allows for higher sampling rates to capture high-frequency interface oscillations critical to the flow regime transitions and works well even for condensation in opaque tubes. Acoustic features (e.g., amplitude, frequency, energy, duration) are correlated with thermofluidic processes (e.g., capillary flows, turbulent flows, boiling, condensation). By relating thermal performance metrics with these dynamic signatures in acoustic and modal regimes, we explore the ability to probe and monitor critical flow regime transitions and transport efficiency in flow condensation.