Search NASA⌕ Search

SEARCH · Search NASA

Results for “SUBCOOLING”

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 73 records · Page 4

Nodal Numerical Modeling of Submerged Helium Injection in a Cryogenic Propellant Tank

Subcooling of cryogenic propellant by helium injection is one of the most effective methods for suppressing bulk boiling and keeping subcooled propellant conditions for pre-launch, launch, and post-launch pressurization applications. For tank pressurization, submerged helium injection can substantially reduce helium consumption by infusing gaseous propellant into the tank ullage. This paper presents a mathematical model of the helium bubbling process in liquid oxygen to estimate the amount of oxygen vapor absorbed by the rising helium bubbles and the amount of subcooling of liquid oxygen due to evaporative heat and mass transfer. This mathematical model was incorporated in a simulation model of tank pressurization built with Generalized Fluid System Simulation Program (GFSSP), a general-purpose flow network code developed at NASA/Marshall Space Flight Center. The numerical predictions of subcooling have been compared with the experimental data of Cho et al. which investigated the propellant subcooling effect as a function of system pressure, helium injection temperature, and flowrate for a non-drained submerged injection system. The numerical predictions of helium consumption have been compared with the test data from a NASA Centaur test vehicle which included both direct and submerged injection with draining of propellants. The numerical model developed with GFSSP has been validated against two sets of experimental data and has been shown to predict both propellant subcooling and helium consumption to within 30% in most cases. The test data used for the model validation were taken in 1-g, but the mass diffusion model was developed to be applied in both 1-g and micro-g environments.

Cryogenics↗

Booster propulsion/vehicle impact study, 2

This is the final report in a study examining the impact of launch vehicles for various boost propulsion design options. These options included: differing boost phase engines using different combinations of fuels and coolants to include RP-1, methane, propane (subcooled and normal boiling point), and hydrogen; variable and high mixture ratio hydrogen engines; translating nozzles on boost phase engines; and cross feeding propellants from the booster to second stage. Vehicles examined included a fully reusable two stage cargo vehicle and a single stage to orbit vehicle. The use of subcooled propane as a fuel generated vehicles with the lowest total vehicle dry mass. Engines with hydrogen cooling generated only slight mass reductions from the reference, all-hydrogen vehicle. Cross feeding propellants generated the most significant mass reductions from the reference two stage vehicle. The use of high mixture ratio or variable mixture ratio hydrogen engines in the boost phase of flight resulted in vehicles with total dry mass 20 percent greater than the reference hydrogen vehicle. Translating nozzles for boost phase engines generated a heavier vehicle. Also examined were the design impacts on the vehicle and ground support subsystems when subcooled propane is used as a fuel. The most significant cost difference between facilities to handle normal boiling point versus subcooled propane is 5 million dollars. Vehicle cost differences were negligible. A significant technical challenge exists for properly conditioning the vehicle propellant on the ground and in flight when subcooled propane is used as fuel.

Johnson, P.↗

Boiling Heat Transfer Mechanisms in Earth and Low Gravity: Boundary Condition and Heater Aspect Ratio Effects

Boiling is a complex phenomenon where hydrodynamics, heat transfer, mass transfer, and interfacial phenomena are tightly interwoven. An understanding of boiling and critical heat flux in microgravity environments is of importance to space based hardware and processes such as heat exchange, cryogenic fuel storage and transportation, electronic cooling, and material processing due to the large amounts of heat that can be removed with relatively little increase in temperature. Although research in this area has been performed in the past four decades, the mechanisms by which heat is removed from surfaces in microgravity are still unclear. Recently, time and space resolved heat transfer data were obtained in both earth and low gravity environments using an array of microheaters varying in size between 100 microns to 700 microns. These heaters were operated in both constant temperature as well as constant heat flux mode. Heat transfer under nucleating bubbles in earth gravity were directly measured using a microheater array with 100 m resolution operated in constant temperature mode with low and high subcooled bulk liquid along with images from below and from the side. The individual bubble departure diameter and energy transfer were larger with low subcooling but the departure frequency increased at high subcooling, resulting in higher overall heat transfer. The bubble growth for both subcoolings was primarily due to energy transfer from the superheated liquid layer relatively little was due to wall heat transfer during the bubble growth process. Oscillating bubbles and sliding bubbles were also observed in highly subcooled boiling. Transient conduction and/or microconvection was the dominant heat transfer mechanism in the above cases. A transient conduction model was developed and compared with the experimental data with good agreement. Data was also obtained with the heater array operated in a constant heat flux mode and measuring the temperature distribution across the array during boiling. The instantaneous heat transfer into the substrate was numerically determined and subtracted from the supplied heat to obtain the wall to liquid heat flux.

Kim, Jungho↗

Nodal Modeling of Submerged Helium Injection Pressurization of a Cryogenic Propellant Tank

Subcooling of cryogenic propellant by helium injection is one of the most effective methods for suppressing bulk boiling and keeping subcooled propellant conditions for pre-launch, launch, and post-launch pressurization applications. For tank pressurization, submerged helium injection can substantially reduce helium consumption by infusing gaseous propellant into the tank ullage. This paper describes a thermodynamic model of the helium bubbling process in liquid oxygen to estimate the amount of oxygen vapor absorbed by the rising helium bubbles and the amount of subcooling of liquid oxygen due to evaporative heat and mass transfer. The physics of helium dissolution during the pressurization process is also modeled, primarily for liquid hydrogen propellant where the dissolution is more significant. The analyses were performed in a simulation model of tank pressurization built with Generalized Fluid System Simulation Program (GFSSP), a general-purpose flow network code developed at NASA/Marshall Space Flight Center. The numerical predictions of subcooling have been compared with the experimental data of Cho et al. which investigated the propellant subcooling effect as a function of system pressure, helium injection temperature, and flowrate for a non-drained submerged injection system. The numerical predictions of helium consumption have been compared with the test data from a NASA Centaur test vehicle which included both direct and submerged injection with draining of propellants. Hydrogen propellant testing of the Cryogenic Propellant Storage and Transfer Engineering Developmental Unit (CPST EDU) conducted at NASA/Glenn Research Facility was also analyzed. The data used for the model validation were taken in 1-g, but the model was developed to be applicable in both multi-g and micro-g environments.

tank pressurization↗

A Mechanistic Study of Nucleate Boiling Heat Transfer Under Microgravity Conditions

Experimental studies of growth and detachment processes of a single bubble and multiple bubbles formed on a heated surface have been conducted in the parabola flights of KC-135 aircraft. Distilled water and PF5060 were used as the test liquids. A micro-fabricated test surface was designed and built. Artificial cavities of diameters 10 microns, 7 microns and 4 microns were made on a thin polished Silicon wafer that was electrically heated by a number of small heating elements on the back side in order to control the surface superheat. Bubble growth period, bubble size and shape from nucleation to departure were measured under subcooled and saturation conditions. Significantly larger bubble departure diameters and bubble growth periods than those at earth normal gravity were observed. Bubble departure diameters as large as 20 mm for water and 6 mm for PF5060 were observed as opposed to about 3 mm for water and less than 1 mm for PF5060 at earth normal gravity respectively. It is found that the bubble departure diameter can be approximately related to the gravity level through the relation D(sub d) proportional 1/g(exp 1/2). For water,the effect of wall superheat and liquid subcooling on bubble departure diameter is found to be small.The growth periods are found to be very sensitive to liquid subcooling at a given wall superheat. However,the preliminary results of single bubble dynamics using PF5060 showed that the departure diameter increases when wall superheat is elevated at the same gravity and subcooling. Growth period of single bubbles in water has been found to vary as t(sub g) proportional g(exp -.93). For water, when the magnitude of horizontal gravitational components was comparable to that of gravity normal to the surface, single bubbles slid along the heater surface and departed with smaller diameter at the same gravity level in the direction normal to the surface. For PF5060, even a very small horizontal gravitational component caused the sliding of bubble along the surface. The numerical simulation has been carried out by solving under the condition of axisymmetry, the mass, momentum, and energy equations for the vapor and the liquid phases. In the model the contribution of micro-layer has been included and instantaneous shape of the evolving vapor-liquid interface is determined from the analysis. Consistent with the experimental results, it is found that effect of reduced gravity is to stretch the growth period and bubble diameter It is found that effect of reduced gravity is to stretch the growth period and bubble diameter at departure. The numerical simulations are in good agreement with the experimental data for both the departure diameters and the growth periods. In the study on dynamics of multiple bubbles, horizontal merger of 2,3 4,and 5 bubbles was observed. It is found that after merger of 2 and 3 bubbles the equivalent diameter of the detached bubble is smaller than that of a single bubble departing at the same gravity level. During and after bubble merger, liquid still fills the space between the vapor stems so as to form mushroom type bubbles. The experimental and numerical studies conducted so far have brought us a step closer to prediction of nucleate boiling heat fluxes under low gravity conditions. Preparations for a space flight are continuing.

Dhir, V. K.↗

Mini-Channel Flow Boiling in 1g and μg

Subcooled flow boiling combines complex processes governed by a strong coupling of heat and mass transfer, hydrodynamics, and liquid-vapor transitions. The gravitational body force plays an important role in these processes. By exploiting the latent heat of vaporization, subcooled flow boiling can transfer up to several orders-of-magnitude more heat compared to forced liquid flow. It is widely used in Earth-based thermal equipment to reduce the cost, size, and power required to move heat. Subcooled flow boiling is now considered a promising technology for space applications since it offers a means to control the motion of bubbles near the heated surface and thereby enhance nucleate boiling. However, a lack of data and predictive models on the role that gravity plays in these phenomena hinders the use of flow boiling in a microgravity environment. We will present results of experiments on transients in nucleate boiling that were carried out in the Flow Boiling Module (FBM) on the International Space Station (ISS). The unique ISS FBM ability is to couple large heat loads with flow visualization in testing flow boiling under strict control of experimental conditions in long-duration microgravity. Comparison of results of microgravity experiments with data obtained in a similar setup on the NASA GRC demonstrates the role of gravity in nucleate boiling. The work is supported by the NASA Flight Opportunities program (grant 80NSSC21K0501) and NSF/CASIS program (NSF CBET grants 1832260, 2126461, 2126462).

Two-phase flow and phase transition↗

Noise analysis of nucleate boiling

The techniques of noise analysis have been utilized to investigate nucleate pool boiling. A simple experimental setup has been developed for obtaining the power spectrum of a nucleate boiling system. These techniques were first used to study single bubbles, and a method of relating the two-dimensional projected size and the local velocity of the bubbles to the auto-correlation functions is presented. This method is much less time consuming than conventional methods of measurement and has no probes to disturb the system. These techniques can be used to determine the contribution of evaporation to total heat flux in nucleate boiling. Also, these techniques can be used to investigate the effect of various parameters upon the frequency response of nucleate boiling. The predominant frequencies of the power spectrum correspond to the frequencies of bubble generation. The effects of heat input, degree of subcooling, and liquid surface tension upon the power spectra of a boiling system are presented. It was found that the degree of subcooling has a more pronounced effect upon bubble size than does heat flux. Also the effect of lowering surface tension can be sufficient to reduce the effect of the degree of subcooling upon the size of the bubbles.

Mcknight, R. D.↗

Centaur Propellant Thermal Conditioning Study

A wicking investigation revealed that passive thermal conditioning was feasible and provided considerable weight advantage over active systems using throttled vent fluid in a Centaur D-1s launch vehicle. Experimental wicking correlations were obtained using empirical revisions to the analytical flow model. Thermal subcoolers were evaluated parametrically as a function of tank pressure and NPSP. Results showed that the RL10 category I engine was the best candidate for boost pump replacement and the option showing the lowest weight penalty employed passively cooled acquisition devices, thermal subcoolers, dry ducts between burns and pumping of subcooler coolant back into the tank. A mixing correlation was identified for sizing the thermodynamic vent system mixer. Worst case mixing requirements were determined by surveying Centaur D-1T, D-1S, IUS, and space tug vehicles. Vent system sizing was based upon worst case requirements. Thermodynamic vent system/mixer weights were determined for each vehicle.

Blatt, M. H.↗

Flow boiling with enhancement devices for cold plate coolant channel design

The use of flow boiling for thermal energy transport is intended to provide an alternative for accommodating higher heat fluxes in commercial space systems. The objectives are to: (1) examine the variations in both the mean and local (axial and circumferential) heat transfer coefficients for a circular coolant channel with either smooth walls, spiral fins, or both spiral fins and a twisted tape; (2) examine the effects of channel diameter and subcooling; and (3) develop an improved reduction analysis and/or suggest possible heat transfer correlation of the present data. Freon-11 is the working fluid. Two-dimensional (circumferential and axial) wall temperature distributions were measured for coolant channels with the above noted internal geometries. The flow regimes which are being studied are: (1) single phase; (2) subcooled flow boiling; and (3) stratified flow boiling. The inside diameter of all test sections is near 1.0 cm. Cicumferentially averaged heat transfer coefficients at several axial locations were obtained for selected coolant channels for a mass velocity of 210 kg/sq m s, an exit pressure of 0.19 MPa (absolute), and an inlet subcooling of 20.8 C. Overall (averaged over the entire channel) heat transfer coefficients were compared for the above channel geometries. This comparison showed that the channel with large pitch spiral fins had higher heat transfer coefficients at all power levels.

Boyd, Ronald D., Sr.↗

Prediction of self-pressurization rate of cryogenic propellant tankage

The SOLA-ECLIPSE code is being developed to enable prediction of the behavior of cryogenic propellants in spacecraft tankage. A brief description of the formulations used for modeling heat transfer and for determining the thermodynamic state is presented. Code performance is verified through comparison to experimental data for the self-pressurization of scale-model liquid hydrogen tanks. SOLA-ECLIPSE is used to examine the effect of initial subcooling of the liquid phase on the self-pressurization rate of an on-orbit full-scale liquid hydrogen tank typical for a chemical-propulsion orbit transfer vehicle. The computational predictions show that even small amounts of subcooling will significantly decrease the self-pressurization rate. Further, if the cooling is provided by a thermodynamic vent system, it is concluded that small levels of subcooling will maximize propellant conservation.

Hochstein, John I.↗

Transient analysis of a capillary pumped loop heat pipe

A bench-top Capillary Pumped Loop (CPL) test system has been developed and tested to investigate the transient mode operation of this system by applying a step power input to the evaporators. Tests were conducted at several power input and evaporator inlet subcooling combinations. In addition, a lumped-heat-capacity model of the CPL test system has been presented which is used for predicting qualitatively the transient operation characteristics. Good agreement has been obtained between the predicted and the measured temperature variations. A simple evaporator inlet subcooler model has also been developed to study effects of inlet subcooling on the steady-state evaporator wall temperature. Results were compared with the test data collected.

Kiper, A. M.↗

Thermocapillary flow with evaporation and condensation and its effect on liquid retention in low-G fluid acquisition devices

The steady motion, thermal and free surface behavior of a volatile, wetting liquid in microgravity are studied using scaling and numerical techniques. The objective is to determine whether the thermocapillary and two-phase convection arising from thermodynamic nonequilibrium along the porous surfaces of spacecraft liquid acquisition devices could cause the retention failures observed with liquid hydrogen and heated vapor pressurant. Why these devices seem immune to retention loss when pressurized with heated helium or heated directly through the porous structure was also examined. Results show that highly wetting fluids exhibit large negative and positive dynamic pressure gradients towards the meniscus interline when superheated and subcooled, respectively. With superheating, the pressure variation and recoil force arising from liquid/vapor phase change exert the same influence on surface morphology and promote retention. With subcooling, however, the pressure distribution produces a suction that degrades mechanical equilibrium of the surface. This result indicates that thermocapillary-induced deformation arising from subcooling and condensation is the likely cause for retention loss. In addition, increasing the level of nonequilibrium by reducing accommodation coefficient suppresses deformation and explains why this failure mode does not occur in instances of direct screen heating or pressurization with a heated inert gas.

Schmidt, George R.↗

Enhancement of Nucleate Boiling Heat Flux on Macro/Micro-Structured Surfaces Cooled by Multiple Impinging Jets

An experimental investigation of nucleate boiling heat transfer from modified surfaces cooled by multiple in-line impinging circular jets is reported and found to agree with single jet results. A copper block is heated from the back by two electrical arcs, and cooled on the opposite side by three identical liquid jets of distilled water at subcoolings of 25 C 50 C and 77 C and Freon 113 at 24 C subcooling. Liquid flow rates are held constant at 5, 10, and 15 GPH for each of the three jets with jet velocities ranging from 1.4 m/s to 1 1.2 m/s and jet diameters from 0.95 mm to 2.2 mm. To increase the maximum heat flux (CHF) and heat removal rate, the boiling surface was modified by both macro and micro enhancements. Macro modification consists of machined radial grooves in the boiling surface arranged in an optimally designed pattern to allow better liquid distribution along the surface. These grooves also reduce splashing of liquid droplets, and provide 'channels' to sweep away bubbles. Micro modification was achieved by flame spraying metal powder on the boiling surface, creating a porous, sintered surface. With the addition of both micro and macro structured enhancements, maximum heat flux and nucleate boiling can be enhanced by more than 200%. Examination of each surface modification separately and together indicates that at lower superheats, the micro structure provides the enhanced heat transfer by providing more nucleation sites, while for higher superheats the macro structure allows better liquid distribution and bubble removal. A correlation is presented to account for liquid subcoolings and surface enhancements, in addition to the geometrical and fluid properties previously reported in the literature.

Kugler, Scott Lee↗

Investigation of Critical Heat Flux in Reduced Gravity Using Photomicrographic Techniques

Experiments were performed to examine the effects of body force on flow boiling critical heat flux (CHF). FC-72 was boiled along one wall of a transparent rectangular flow channel that permitted photographic study of the vapor-liquid interface just prior to CHF. High-speed video imaging techniques were used to identify dominant CHF mechanisms corresponding to different flow orientations and liquid velocities. Six different CHF regimes were identified: Wavy Vapor Layer, Pool Boiling, Stratification, Vapor Counterflow, Vapor Stagnation, and Separated Concurrent Vapor Flow. CHF showed significant sensitivity to orientation for flow velocities below 0.2 m/s, where extremely low CHF values where measured, especially with downward-facing heated wall and downflow orientations. High flow velocities dampened the effects of orientation considerably. The CHF data were used to assess the suitability of previous CHF models and correlations. It is shown the Interfacial Lift-off Model is very effective at predicting CHF for high velocities at all orientations. The flooding limit, on the other hand, is useful at estimating CHF at low velocities and for downflow orientations. A new method consisting of three dimensionless criteria is developed for determining the minimum flow velocity required to overcome body force effects on near-saturated flow boiling CHF. Vertical upflow boiling experiments were performed in pursuit of identifying the trigger mechanism for subcooled flow boiling CHF. While virtually all prior studies on flow boiling CHF concern the prediction or measurement of conditions that lead to CHF, this study was focused on events that take place during the CHF transient. High-speed video imaging and photomicrographic techniques were used to record the transient behavior of interfacial features from the last steady-state power level before CHF until the moment of power cut-off following CHF. The video records show the development of a wavy vapor layer which propagates along the heated wall, permitting cooling prior to CHF only in wetting fronts corresponding to the wave troughs. Image analysis software was developed to estimate void fraction from the individual video images. The void fraction records for subcooled flow boiling show the CHF transient is accompanied by gradual lift-off of wetting fronts culminating in some maximum vapor layer mean thickness, following which the vapor layer begins to thin down as the transition to film boiling ensues. This study proves the Interfacial Lift-off Model, which has been validated for near-saturated flow boiling CHF, is equally valid for subcooled conditions.

Mudawar, Issam↗

Geoscience Laser Altimetry System (GLAS) On-Orbit Flight Report on the Propylene Loop Heat Pipes (LHPs)

The Geoscience Laser Altimetry System (GLAS) instrument which is the sole instrument on ICESat was launched on January 12, 2003. GLAS utilizes two actively controlled propylene Loop Heat Pipes (LHPs) as the core of its thermal system. The LHPs started quickly when the Dale Ohm starter heaters were powered and have as designed. The low control heater power and on-orbit tight temperature control appear independent of gravity effects when comparing ground testing to flight data. The use of coupling blocks was also unique to these LHPs. Their application reduced control heater power by reducing the subcooling from the radiator. The effectiveness in reducing subcooling of the coupler blocks decreased during flight from ground testing, but internal thermal isolation in the compensation chamber between the subcooled returning liquid increased in flight resulting in no net increase in control heater power versus ground measurements. Overall the application of LHPs in the thermal system for GLAS met instrument requirements and provided flexibility for the overall system as last minute requirements became known.

Baker, Charles L.↗

A Densified Liquid Methane Delivery System for the Altair Ascent Stage

The Altair Lunar Lander is currently carrying options for both cryogenic and hypergolic ascent stage propulsion modules. The cryogenic option uses liquid methane and liquid oxygen to propel Altair from the lunar surface back to rendezvous with the Orion command module. Recent studies have determined that the liquid methane should be densified by subcooling it to 93 K in order to prevent over-pressurization of the propellant tanks during the 210 day stay on the lunar surface. A trade study has been conducted to determine the preferred method of producing; loading, and maintaining the subcooled, densified liquid methane onboard Altair from a ground operations perspective. The trade study took into account the limitations in mass for the launch vehicle and the mobile launch platform as well as the historical reliability of various components and their thermal efficiencies. Several unique problems were encountered, namely delivering a small amount of a cryogenic propellant to a flight tank that is positioned over 350 ft above the launch pad as well as generating the desired delivery temperature of the methane at 93 K which is only 2.3 K above the methane triple point of 90.7 K. Over 20 methods of subcooled liquid methane production and delivery along with the associated system architectures were investigated to determine the best solutions to the problem. The top four cryogenic processing solutions were selected for further evaluation and detailed thermal modeling. This paper describes the results of the preliminary trade analysis of the 20 plus methane densification methods considered. The results of the detailed analysis will be briefed to the Altair Project Office and their propulsion team as well as the Ground Operations Project Office before the down-select is made between cryogenic and hypergolic ascent stages in August 2010.

Tomsik, Thomas M.↗

Modeling and Simulation of Tank Pressure Control using Zero-Boil Off Active Thermal Control for LOXSAT Technology Demonstration Mission

To-date, research and modeling of cryogenic fluid management technologies (CFM) for spaceflight has been limited to ground tests, short-duration zero-g simulations (e.g. drop towers), and small-scale experiments on-orbit. There has not been a large-scale flight demonstration of a flight-like system. As future NASA missions to take humans further from Earth will require large, in-space cryogenic propulsion vehicles, it is imperative to begin collecting flight data for these systems to accurately model and design future vehicles. To meet this goal, NASA awarded tipping point technology demonstration awards to Eta Space, Lockheed Martin, Space Exploration Technologies (SpaceX), and United Launch Alliance (ULA) to demonstrate on-orbit storage and transfer of cryogenic propellant. For its award, Eta Space is developing LOXSAT-1. It is a small satellite that will be launched on a Rocket Lab Electron rocket. The spacecraft consists of a Rocket Lab Photon spacecraft bus with a primary payload of a spherical liquid oxygen (LOX) storage tank with thermodynamic control systems. The primary objective of the mission is to demonstrate zero-boil-off storage of liquid oxygen. To accomplish this objective, the payload is equipped with an active thermal control system fluid loop that consists of propellant management device (PMD), positive-displacement pump, heat exchanger connected to a cryocooler, and a spray bar mixing injector. When the fluid loop is operating, fluid is drawn from the tank by the PMD and pumped through the heat exchanger, lowering the fluid temperature below the fluid temperature in the tank. This subcooled liquid is then injected back into the tank through the spray bar. The subcooled injected liquid has two effects. If it is sprayed into the ullage space, the injected liquid will form into jets or droplets and exchange heat with the ullage gas. This will cool and condense the gas, reducing the pressure in the tank. Additionally, the liquid that is not sprayed through the ullage, as well as any remaining liquid spray from the ullage, will rejoin the liquid mass of the tank, lowering the bulk temperature of the liquid. These combines effects provide for zero-boil off pressure control by lowering the tank pressure and the liquid saturation pressure simultaneously, ensuring the liquid stays subcooled. To model these complex mechanics and predict the performance of the active thermal control system, NASA is providing Eta Space with 3 parallel models of the tank. The approach of providing 3 different models allows for cross-checking and comparisons between the three to better understand how different modeling assumptions and selection semi-empirical factors affects the modeling result. Additionally, developing 3 models provides three different schemes for numerical simulation, providing confidence that results depict real physical phenomenon and not numerical quirks of the program. Within the tank thermodynamics, there are two primary areas of heat transfer we concern ourselves with: the heat transfer between the ullage space and the droplet spray, and between the ullage space and bulk liquid. For the droplet heat transfer, there are multiple sets of assumptions that can be made and correlations that can be used. Currently, two working models - the TankSIM model and Easy5 model – provide for an overview of the different approaches available. The TankSIM model and Easy5 model use two different models for droplet heating and evaporation that illustrate how the models use different types of mechanisms to arrive at the same answer. For the TankSIM model, droplets are treated as spheres of constant radius. Heat is transferred from the ullage to the droplet and warms the droplet until it reaches saturation, then the droplet begins evaporating and reducing its radius and mass. To calculate the heat transfer coefficient between the droplet and ullage, the Ranz-Marshall correlation is used. To determine the number of droplets in the ullage, a resident mass approach is used. This approach averages the number of droplets such that residuals at the start-up and shut-off of the spray bar cancel out. This same approach is used in the Easy5 model. The GFSSP model implements a linked list to track individual droplet “nodes” within the model. For the Easy5 model, the droplet is assumed to have an interface at a temperature equal to the saturation temperature corresponding to the pressure of the gas phase. The heat transfer from the gas to the interface and the interface to the droplet bulk is then calculated, and the net mass transfer between the droplet and interface is determined by performing an energy balance across the interface. For the gas side of the interface, the Ranz-Marshall correlation is used. For the liquid side of the interface, a variety of correlations were tried, including Kronig and Brink (1950) and effective conductivity models. As a result of these assumptions, the Easy5 model currently predicts faster depressurization, as at saturated vapor conditions, the heat transfer coeffect on the liquid side for the Easy5 model is greater than the heat transfer coefficient predicted by Ranz-Marshall used in the TankSIM code. This greater heat flux translates into faster condensation of the saturated ullage gas. At the bulk liquid to ullage interface, the models are in much closer agreement. Both models model the ullage as a sphere centered within the bulk liquid in the tank, and both use the energy-jumping boundary condition to model heat and mass transfer across the interface. There are slight differences in how the interface temperature is calculated, however. The Easy5 model assumes the temperature of the interface is equal to the saturation temperature associated with the pressure of the gas phase. The TankSIM model calculates this temperature with Alabovskii’s equation, which provides a correction factor for interface temperatures. Analysis tasks are focused on determining rates of depressurization within the tank during active cooling operation. To maintain net positive suction head at the pump inlet, the tank pressure cannot fall faster than the saturation pressure associated with the temperature of the bulk liquid. Additionally, there is interest in analyzing the performance of the loop at different pump speeds and cryocooler input powers. Adjusting the flowrate affects both the performance of the heat exchanger between the cryocooler and pumped liquid, and the heat transfer between the droplet spray and the ullage. Ideally, a pump speed and cryocooler power can be selected that will allow the tank to operate in zero-boil-off mode with a very narrow range of storage pressure.

Cameron J. Hines↗

Modeling and Simulation of Tank Pressure Control using Zero-Boil Off Active Thermal Control for LOXSAT Technology Demonstration Mission

To-date, research and modeling of cryogenic fluid management technologies (CFM) for spaceflight has been limited to ground tests, short-duration zero-g simulations (e.g. drop towers), and small-scale experiments on-orbit. There has not been a large-scale flight demonstration of a flight-like system. As future NASA missions to take humans further from Earth will require large, in-space cryogenic propulsion vehicles, it is imperative to begin collecting flight data for these systems to accurately model and design future vehicles. To meet this goal, NASA awarded tipping point technology demonstration awards to Eta Space, Lockheed Martin, Space Exploration Technologies (SpaceX), and United Launch Alliance (ULA) to demonstrate on-orbit storage and transfer of cryogenic propellant. For its award, Eta Space is developing LOXSAT-1. It is a small satellite that will be launched on a Rocket Lab Electron rocket. The spacecraft consists of a Rocket Lab Photon spacecraft bus with a primary payload of a spherical liquid oxygen (LOX) storage tank with thermodynamic control systems. The primary objective of the mission is to demonstrate zero-boil-off storage of liquid oxygen. To accomplish this objective, the payload is equipped with an active thermal control system fluid loop that consists of propellant management device (PMD), positive-displacement pump, heat exchanger connected to a cryocooler, and a spray bar mixing injector. When the fluid loop is operating, fluid is drawn from the tank by the PMD and pumped through the heat exchanger, lowering the fluid temperature below the fluid temperature in the tank. This subcooled liquid is then injected back into the tank through the spray bar. The subcooled injected liquid has two effects. If it is sprayed into the ullage space, the injected liquid will form into jets or droplets and exchange heat with the ullage gas. This will cool and condense the gas, reducing the pressure in the tank. Additionally, the liquid that is not sprayed through the ullage, as well as any remaining liquid spray from the ullage, will rejoin the liquid mass of the tank, lowering the bulk temperature of the liquid. These combines effects provide for zero-boil off pressure control by lowering the tank pressure and the liquid saturation pressure simultaneously, ensuring the liquid stays subcooled. To model these complex mechanics and predict the performance of the active thermal control system, NASA is providing Eta Space with 3 parallel models of the tank. The approach of providing 3 different models allows for cross-checking and comparisons between the three to better understand how different modeling assumptions and selection semi-empirical factors affects the modeling result. Additionally, developing 3 models provides three different schemes for numerical simulation, providing confidence that results depict real physical phenomenon and not numerical quirks of the program. Within the tank thermodynamics, there are two primary areas of heat transfer we concern ourselves with: the heat transfer between the ullage space and the droplet spray, and between the ullage space and bulk liquid. For the droplet heat transfer, there are multiple sets of assumptions that can be made and correlations that can be used. Currently, two working models - the TankSIM model and Easy5 model – provide for an overview of the different approaches available. The TankSIM model and Easy5 model use two different models for droplet heating and evaporation that illustrate how the models use different types of mechanisms to arrive at the same answer. For the TankSIM model, droplets are treated as spheres of constant radius. Heat is transferred from the ullage to the droplet and warms the droplet until it reaches saturation, then the droplet begins evaporating and reducing its radius and mass. To calculate the heat transfer coefficient between the droplet and ullage, the Ranz-Marshall correlation is used. To determine the number of droplets in the ullage, a resident mass approach is used. This approach averages the number of droplets such that residuals at the start-up and shut-off of the spray bar cancel out. This same approach is used in the Easy5 model. The GFSSP model implements a linked list to track individual droplet “nodes” within the model. For the Easy5 model, the droplet is assumed to have an interface at a temperature equal to the saturation temperature corresponding to the pressure of the gas phase. The heat transfer from the gas to the interface and the interface to the droplet bulk is then calculated, and the net mass transfer between the droplet and interface is determined by performing an energy balance across the interface. For the gas side of the interface, the Ranz-Marshall correlation is used. For the liquid side of the interface, a variety of correlations were tried, including Kronig and Brink (1950) and effective conductivity models. As a result of these assumptions, the Easy5 model currently predicts faster depressurization, as at saturated vapor conditions, the heat transfer coeffect on the liquid side for the Easy5 model is greater than the heat transfer coefficient predicted by Ranz-Marshall used in the TankSIM code. This greater heat flux translates into faster condensation of the saturated ullage gas. At the bulk liquid to ullage interface, the models are in much closer agreement. Both models model the ullage as a sphere centered within the bulk liquid in the tank, and both use the energy-jumping boundary condition to model heat and mass transfer across the interface. There are slight differences in how the interface temperature is calculated, however. The Easy5 model assumes the temperature of the interface is equal to the saturation temperature associated with the pressure of the gas phase. The TankSIM model calculates this temperature with Alabovskii’s equation, which provides a correction factor for interface temperatures. Analysis tasks are focused on determining rates of depressurization within the tank during active cooling operation. To maintain net positive suction head at the pump inlet, the tank pressure cannot fall faster than the saturation pressure associated with the temperature of the bulk liquid. Additionally, there is interest in analyzing the performance of the loop at different pump speeds and cryocooler input powers. Adjusting the flowrate affects both the performance of the heat exchanger between the cryocooler and pumped liquid, and the heat transfer between the droplet spray and the ullage. Ideally, a pump speed and cryocooler power can be selected that will allow the tank to operate in zero-boil-off mode with a very narrow range of storage pressure.

zero boil-off↗