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 91 records · Page 5

Progress Toward Simulating Departure from Nucleate Boiling at High-Pressure Applications with Selected Wall Boiling Closures

Recently, a Eulerian-based two-fluid computational fluid dynamics (CFD) framework with a wall heat flux partitioning approach has been intensively investigated for departure from nucleate boiling (DNB) simulation under the U.S. Department of Energy–funded Consortium for Advanced Simulation of Light Water Reactors (CASL) program. Understanding of the DNB characteristics over a range of pressurized water reactor–like operating conditions and accurate prediction of boiling crisis in the nuclear power system have been grand challenges because of the large impact of DNB on reactor safety and operational economics. The ultimate goal of this task in the CASL program is to introduce a robust multiphase CFD–based DNB modeling framework that is capable of characterizing an entire boiling history in which the wall boiling mode experiences the following through multiple stages of heat transfer mode: (1) single-phase convective heat transfer, (2) nucleate boiling heat transfer, and (3) identification of the departure of nucleate boiling. To validate the CASL boiling model, we have benchmarked simulated DNB over three different flow channel configurations (pipe flow, 5 × 5 fuel bundle with mixing vane tests, and 5 × 5 fuel bundle without mixing vane tests) against experimental measurements, and the validation result with open literature is reported. The DNB detection criteria in the simulation are checked by monitoring the peak wall temperature, wall dryout factor, and net energy balance. In addition to the DNB performance test, some preliminary sensitivity results on closure model selection are reported to address the prediction capability of local void profile against measurements. The boiling simulation tested in this study exhibits a maximum deviation of 24% from the measured DNB value in a high-pressure (i.e., 138 bars) subcooled pipe flow test. The ranges of operating conditions are as follows: 1650 to 2650 kg/m 2 ·s for mass flux and 8.5 to 96 K for subcooled inlet temperature. The deviation is even reduced to 7% when the subcooled temperature is less than 40 K. Besides accuracy, base practice guidelines for DNB detection criteria are tested by monitoring three simulation variables: (1) maximum wall temperature, (2) wall dryout factor (i.e., K-value), and (3) energy balance. Numerical robustness of DNB simulation is largely achieved in most of the validation test except for a few high subcooled test cases.

42 ENGINEERING↗

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.↗

Field Demonstration and Adoption Impacts of an Affordable Thermoelectric Heat Pump

This study presents the design, field validation, and utility impacts of an air source heat pump that uses thermoelectric (TE) technology to provide supplemental heating at outside air temperatures below 17°F. Conventional heat pumps experience reduced heating capacity and COP at low outdoor temperatures due to the temperature difference between the source and demand sides of the vapor compression cycle. Consequently, during colder months, supplemental heating such as electric resistance or natural gas furnace is often needed to provide additional warmth, which increases the overall equipment and energy costs. To address this challenge, a TE subcooler (an array of solid-state, thermoelectric modules that provides heat pumping using non-vapor compression technique) is fabricated and integrated with the vapor compression cycle of a conventional air-source heat pump. In colder ambient conditions, TE subcooler exchanges heat with the refrigerant to provide additional subcooling, thereby enhancing the capacity of the heat pump. The operation of thermoelectric-integrated heat pump was demonstrated in an occupied single-family home in Nashville, TN. The TE subcooler increased the system’s heating capacity at 15°F outdoor temperature by 25–30% and reduced electricity demand by 10%. It also met the heating load in the house without supplemental electric resistance down to 5°F outdoors, delivering the additional heat at nearly 80% greater efficiency compared to electric resistance heating. With a modest first cost increment of 10–15% over conventional single-speed heat pumps, this approach offers a promising, energy-efficient solution to accelerate residential heat pump adoption in the U.S.

KRISHNAMOORTHY, Sreenidhi [Electric Power Research↗

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↗

Virtual refrigerant charge sensor for variable-speed heat pumps based on feature selection

The refrigerant charge level in heat pump systems significantly impacts their energy efficiency. Virtual refrigerant charge (VRC) sensing technology has been comprehensively investigated and well-established due to its lower cost compared to physical sensors. However, the previous VRC research often relied on expert judgment and physical reasoning for their variable selection, which can potentially select redundant (or highly correlated) or insignificant features, and it is also primarily focused on single-speed systems. To address these challenges, this study proposes a VRC algorithm for variable-speed heat pumps that selects features through a rigorous feature selection method in combination with physical insights. We also propose a piecewise linear model structure segmented by subcooling temperature to accurately predict charge levels, particularly when subcooling temperatures are substantially low. The proposed algorithm was evaluated using experimental data of a residential R410A heat pump, and the performance was compared with two baseline VRC algorithms. The results are: (1) The proposed algorithm outperforms for the case with subcooling temperature less than 1 °C. (2) The proposed algorithm achieves a tested mean absolute percentage error (MAPE) of 4.23%, and improves the overall accuracy for cooling conditions by approximately 60%, compared with the two baseline algorithms. (3) The proposed algorithm uses two fewer features and improves the accuracy for undercharge cooling conditions by 68.0%, compared with baseline algorithm 2. These improvements enhance prediction accuracy and prevent overfitting, providing a more reliable refrigerant charge level prediction and helping improve the heat pump energy efficiency.

Liang, Chenjiyu↗

An experimentally validated evaporative phase change heat transfer model for low mass flux applications using R134a in plate heat exchangers

Plate frame heat exchangers are common in liquid-coupled vapor compression systems due to their compactness and ease of maintenance. In some of these systems, the refrigerant within the evaporator can enter the heat exchanger as a subcooled liquid. As a result, the refrigerant passes through three different phases: single-phase liquid, two-phase fluid, and a superheated vapor. Unfortunately, there have been limited prior studies that have identified the best method to predict performance under these circumstances. Furthermore, previous investigations have not evaluated refrigerant evaporation in large industrial sized plate frame heat exchangers under conditions experienced in the present study experimental test facility, especially at mass fluxes below 7 kg m –2 s –1 . In the present investigation, a model was developed to predict the performance of plate and frame heat exchangers when the fluid enters as a subcooled liquid and exits as a superheated vapor. The model used full-sized plate heat exchanger geometry that was discretized into 20 sections to accurately capture local heat transfer and pressure drop effects. The model was validated using an R134a counter flow heat exchanger used in a turbo-compression cooling system test facility that had subcooled liquid entering the evaporator at low mass fluxes. A variety of empirical correlations were evaluated to determine what combination yielded the best predictive capability over the following range of conditions: 5.8 < G < 6.8 kg m –2 s –1 , 2 < q” < 2.8 kW m –2 , and 10 °C < T sat < 15 °C. When the correlations used a fixed control volume length in the heat flux calculation, the most accurate combination resulted in a mean absolute error of 5.5%. Finally, future studies can use the approach described here to optimize heat exchanger size and performance.

42 ENGINEERING↗

Preliminary Thermal-Hydraulic Analysis of the UTA–2 Subcritical Linear Accelerator–Driven System

The National Nuclear Security Administration’s mission includes establishing a reliable supply of 99 Mo without highly enriched uranium. Oak Ridge National Laboratory (ORNL) supports this objective through collaborative research and development with industrial partners. Niowave Inc., a current partner, is currently designing a subcritical linear accelerator-driven system (ADS) and preparing for the US Nuclear Regulatory Commission’s licensing process. Niowave’s system has the potential to efficiently supply medical radioisotopes. The technology includes a superconducting electron accelerator and a pile of both natural and low-enriched uranium (LEU) targets. The process fissions uranium and many valuable isotopes can then be extracted from the targets. Niowave is currently iterating through conceptual and detailed design processes for several system sizes. This report discusses UTA–2, which is at the demonstration stage. UTA–2 will validate numerical modeling results with experimental measurements before progressing to the detailed design of UTA–3, the commercial-sized ADS. The thermal-hydraulic behavior of the UTA–2 core design was numerically investigated using STAR-CCM+, as described in this report. STAR-CCM+, a state-of-the-art computational fluid dynamics (CFD) software that was commercially developed by Siemens, has an extensive user base and a set of validation studies. It is also compliant with the American Society of Mechanical Engineers’ Nuclear Quality Assurance 1 standard. Three cases are investigated in this report. Case 1 quantified the temperature field within the UTA–2 assembly and water tank using only conduction as the method of thermal energy transport. This simplified approach was overly conservative and yielded wetted cladding temperatures above the coolant saturation temperature. In this case, the maximum temperature of the wetted cladding surface of the highest power rod exceeded the saturation temperature of water by 63.8°C. Because of the overly conservative approach taken in Case 1 and its negative subcooled margin, buoyancy-driven natural circulation flow physics were implemented in Case 2. Adding coolant motion significantly distributed the thermal energy of the system through convective heat transfer. This relatively small amount of convective heat transfer significantly reduced system temperatures and increased the subcooled margin from –63.8 to 61.3°C. This margin confirmed that no boiling was expected during normal operation of UTA–2 at 230 W. Case 3 had no additional physics models but considered an overpower event during which the power of each LEU and natural uranium rod was at its respective peak values. This resulted in a study with the total assembly power equal to 176% of the nominal power of 230 W considered in Cases 1 and 2. The resulting natural circulation flows were slightly enhanced. The subcooled margin decreased slightly to 46.3°C, which is still a significant margin to local boiling of the water within the tank. This margin confirmed that no boiling was expected during an abnormal operation of UTA–2 at 406 W.

07 ISOTOPE AND RADIATION SOURCES↗

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↗