SEARCH · Search NASA
Results for “Spray Boiling”
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Evaluation of Spray and Combustion Models for Simulating Dilute Combustion in a Direct-Injection Spark-Ignition Engine
Dilute combustion in spark-ignition engines has the potential to improve thermal efficiency by mitigating knock and by reducing throttling and wall heat losses. However, ignition and combustion processes can become unstable for dilute operation due to a lowered laminar flame speed, resulting in excessive cycle-to-cycle variability (CCV) of the combustion process. To compensate for the slower combustion in less reactive mixtures, a modified intake port geometry can be employed to generate a strong tumble flow in the cylinder and elevate turbulence levels around the spark plug, thereby promoting a faster transition to turbulent deflagration. Consequently, optimizing combustion chamber geometry and operating strategy is crucial to maximizing the benefits of using dilute combustion with enhanced in-cylinder turbulence across a wide range of operating conditions. Computational fluid dynamics (CFD) simulations can be utilized for virtual engine optimization tasks, but this would require the models to be truly predictive regarding the impact of changes to the engine design and operational parameters.In this study, multicycle large-eddy simulations (LES) are performed for a direct-injection spark-ignition engine to investigate the model performance in predicting engine combustion characteristics with respect to changes in the intake configuration. A tumble plate that blocks the lower part of the intake port inlet is used to vary the tumble. A set of CFD models that have been recently developed are employed, which takes into account the drag of nonspherical droplets, flash-boiling behavior of liquid sprays, spray-wall interaction, surrogate formulation of a research-grade E10 gasoline, and fast chemical kinetic solvers. Simulation results are compared to experimental engine data in terms of cylinder pressure, apparent heat release rate, mass fraction burned timing, and flame images. It is found that LES employing the state-of-the-art CFD models are capable of properly predicting the spray processes and reproducing the measured mean cylinder pressure for the case with the tumble plate. On the other hand, the LES over-predicts the combustion rate during the early combustion stage and under-estimates the CCV, and these discrepancies become larger when the tumble plate is removed.
Numerical Simulation on Liquid Hydrogen Chill-Down Process of Vertical Pipeline
In order to improve the cryogenic propellant management technologies for a liquid hydrogen rocket with high specific impulse, JAXA, the University of Tokyo, and the NASA Glenn Research Center have jointly organized a multi-agency model validation collaboration project. As part of this project, JAXA's boiling simulation was validated with NASA's experimental data on vertical pipeline chill-down. Simulation results were in good agreement with the experimental data obtained using an improved boiling model to reproduce the spray flow. This activity achieved liquid hydrogen turbo-pump simulation at JAXA for grasping the boiling flow phenomenon from engine cut-off to re-ignition. This joint research resulted in an international cooperative relationship for discussing the cryogenic propellant management technologies necessary to develop next-generation liquid rockets.
A Phenomenological Thermal Spray Wall Interaction Modeling Framework Applied to a High Temperature Ignition Assistant Device
Airborne compression ignition engines must operate with reliable ignition systems to achieve proper ignition at every cycle, particularly at high altitudes. Glow-plug-based ignition-assistant (IA) devices can provide the necessary energy to preheat the fuel and ensure ignitability of the fuel-air mixture. Ignitability of liquid sprays can be facilitated via direct impingement onto the hot IA surface, however this comes with adverse effects on the IA durability. Therefore, optimizing an IA’s design requires detailed understanding of the physics of fuel spray impingement of superheated surfaces. While spray impingement on relatively low wall temperatures has been extensively studied and appropriate numerical models have been proposed through the years, fundamental understanding of high-speed liquid spray impingement on superheated walls is still elusive. This work aims to formulate a phenomenological thermal spray-wall interaction framework for modeling the film-boiling-induced heat transfer, atomization, and dispersion of fuel spray droplets impinging onto a superheated IA device. A qualitative comparison of the new phenomenological model is performed against optical experiments from the literature of an F-24 fuel spray injected onto an IA device located 12 mm away from the injector tip. The temperature of the IA was set at 1400 K. The fuel injection pressure was 400 bar, while the ambient gas pressure and temperature were 30 bar and 800 K, respectively. The performance of the phenomenological model is evaluated in comparison with two other state-of-art models from the literature. A qualitative analysis of the different spray and fuel-air mixture characteristics is performed to outline the differences in the predictions offered by the new phenomenological model and the two state-of-art spray-wall interaction models.
Assessment of an E10 gasoline surrogate: Qualitative and quantitative comparisons of in-cylinder spray morphology
A minimum-component gasoline fuel surrogate that captures both chemical and physical behaviors of a full-distillate fuel is needed for high-fidelity CFD simulations. This study evaluates gasoline spray characteristics in a direct-injection spark-ignition engine under motored operation. Two fuels are compared; PACE-20, which is a 9-component surrogate formulation of RD5-87, is compared with its target fuel RD5-87, which is a full-boiling range research grade E10 gasoline. The spray morphologies of both fuels are recorded for a centrally-located direct-injection 8-hole spray subject to intake air cross-flow during the early part of the intake stroke. High-speed imaging recorded scattered light of the side and axial projections of the liquid spray. Quantitative metrics were developed and employed to facilitate comparison of spray morphologies as well as to identify the transition in spray morphology due to flash boiling. This paper builds on a previous study of RD5-87 where coolant temperature (20°C–100°C), in-cylinder pressure (40–110 kPa), engine speed (650–1950 rpm), and injection pressure (60–180 bar) were systematically changed to span operating conditions with and without flash boiling. Images of the PACE-20 morphology are selected for a sub-set of operating conditions from the previous study where distinctive morphology changes occurred. Visual inspection of the images and quantitative metrics demonstrate that the PACE-20 spray morphology is equivalent to that of the RD5-87 in most cases. The exception was for changes in the ambient-gas pressure where the flash-boiling transition occurred at ∼5 kPa higher in-cylinder pressure for PACE-20. Three empirical metrics, Merging Index, Asymmetry, and Flash Index are proposed here and they were found to be useful both as quantitative comparisons of the fuel morphologies, and for identifying the transition in spray morphology due to flash boiling.
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.
A Summary of Recent SBIR and STTR Projects in Support of Cryogenic Fluid Management Modeling
Over the last five years, NASA has invested in development of sub-models to improve the current capabilities of Computational Fluid Dynamics (CFD) codes utilized for cryogenic fluid management (CFM) predictions. The goal is for the submodels to be modular so that they can easily be implemented into NASA’s CFD codes such as ANSYS Fluent, Siemens STAR-CCM+, and Streamline Numerics Loci-Stream. These technologies are being developed through Small Business Innovative Research (SBIR) and Small Business Technology Transfer (STTR) contracts under the CFM subtopic solicitation and funded by the CFM Portfolio Project. All contracts presented include cryogenic experiments to enable model validation. A sub-grid boiling model has been developed by Combustion Research and Flow Technology Inc. and Massachusetts Institute of Technology for cryogenic line chilldown predictions and has completed Phase III. A spray chill-and-fill model for cryogenic tankage is being developed by Combustion Research and Flow Technology Inc. and University of Connecticut and is currently in Phase II. Last, film condensation modeling to support cryogenic liquefaction was progressed through a Phase I contract with Combustion Research and Flow Technology Inc. and University of Connecticut. A summary of accomplishments for each contract is presented.
System Study: High-Pressure Core Spray 1998-2020
This report presents an unreliability evaluation of the high-pressure core spray (HPCS) at eight U.S. commercial boiling water reactors. Demand, run hour, and failure data from 1998–2020 for selected components were obtained from the Institute of Nuclear Power Operations (INPO) Industry Reporting and Information System (IRIS), formerly the INPO Consolidated Events Database (ICES). The unreliability results are trended for the most recent 10 year period while yearly estimates for system unreliability are provided for the entire active period. No statistically significant increasing or decreasing trends were identified in the HPCS results.
Stress corrosion cracking mechanism of cold spray coating on a galvanically similar substrate
Here, the chloride-induced stress corrosion cracking (CISCC) mechanism of cold spray (CS) coating on a galvanically similar substrate is investigated. Arc welded 304L stainless steel (SS) specimens are loaded into four-point bend fixtures, cold sprayed with 304L SS, then immersed in boiling MgCl 2 . Interconnected porosity forms through crevice corrosion along CS splat boundaries, allowing corrosive species to penetrate through the CS layer. Nevertheless, the substrate is resistant to CISCC likely because of compressive stress introduced by peening during CS particle impacts. These findings underscore the importance of residual stress in the environmental degradation of CS coatings or repairs of engineering structures.
System Study: High-Pressure Core Spray 1998-2022
This report presents an unreliability evaluation of the high-pressure core spray (HPCS) at eight U.S. commercial operating boiling water reactors. New Standardized Plant Analysis Risk (SPAR) models with the most recent SPAR parameter update results were used in this report. Demand, run hour, and failure data from 1998–2022 for selected components were obtained from the Institute of Nuclear Power Operations Industry Reporting and Information System. The unreliability results are trended for the most recent 10-year period while yearly estimates for system unreliability are provided for the entire active period. Statistically significant increasing trends were identified in both the HPCS system start-only unreliability and 8-hour mission unreliability.
System Study: High-Pressure Core Spray 1998-2024
This report presents an unreliability evaluation of the high-pressure core spray (HPCS) at eight U.S. commercial operating boiling water reactors. New Standardized Plant Analysis Risk (SPAR) models with the most recent SPAR parameter update results were used in this report. Demand, run hour, and failure data from 1998 to 2024 for selected components were obtained from the Institute of Nuclear Power Operations Industry Reporting and Information System. The unreliability results are trended for the most recent 10-year period and yearly estimates for system unreliability are provided for the entire active period. No statistically significant increasing or decreasing trends were identified in the HPCS results.
Testing of a Spray-Bar Zero Gravity Cryogenic Vent System for Upper Stages
The capability to vent in zero gravity without resettling is a fundamental technology need that involves practically all uses of subcritical cryogenics in space. Venting without resettling would extend cryogenic orbital transfer vehicle capabilities. However, the lack of definition regarding liquid/ullage orientation coupled with the somewhat random nature of the thermal stratification and resulting pressure rise rates, lead to significant technical challenges. Typically a zero gravity vent concept, termed a thermodynamic vent system (TVS), consists of a tank mixer to destratify the propellant, combined with a Joule- Thomson (J-T) valve to extract then-nal energy from the propellant. In a cooperative effort, Marshall Space Flight Center's (MSFC's) Multipurpose Hydrogen Test Bed (N4HTB) was used to test a unique "spray bar" TVS system developed by Boeing. A schematic of this system is included in Figure 1. The system consists of a recirculation pump, a parallel flow concentric tube, heat exchanger, and a spray bar positioned close to the longitudinal axis of the tank. In the mixing mode, the recirculation pump withdraws liquid from the tank and sprays it radially into the tank liquid, ullage, and exposed tank surfaces. When energy extraction is required, a small portion of the recirculated liquid is passed sequentially through the J-T expansion valve, the spray bar heat exchanger element, and is vented overboard. The vented vapor cools the circulated bulk fluid, thereby removing thermal energy and reducing tank pressure. Figure 2 is a plot of ullage pressure (P4) and liquid vapor pressure (PSAI) versus time. The pump operates alone, cycling on and off, to destratify the tank liquid and ullage until the liquid vapor pressure reaches the lower set point. At that point, the J-T valve begins to cycle on and off with the pump. Thus, for short duration missions, only the mixer may operate, thus minimizing or even eliminating boil-off losses. The primary advantage of the spray bar configuration is that pressure reduction is achieved independent of liquid and vapor location, thereby enhancing the applicability of normal gravity test data to zero gravity conditions. The in-tank components are minimized with the proposed TVS design. Because the recirculation pump is external to the tank, no electrical power penetration of the tank is required for pump or valve operation. This is especially desirable for L02 tanks since the presence of an electrical ignition source in oxygen represents a critical failure mode. Also, since the critical components (pump, motor, valve, orifice) are external to the tank, system checkout and ground servicing/replacement are easier. For zero-g operation, component replacement external to the tank may be a significant benefit. In addition to satisfying the zero g TVS design objectives, the TVS concept tested offers additional benefits to the integrated subcritical cryogenic storage and launch system.
Axisymmetric Two-Dimensional Modeling of No Vent Filling of a Cryogenic Tank using Generalized Fluid System Simulation Program
This paper presents preliminary analysis using a finite volume based computational method for network flow analysis for the simulation of chill and fill of a cryogenic tank. An axisymmetric two-dimensional model of a spherical tank was developed using a nodal approach. No Vent Filling (NVF) was accomplished using a Thermodynamic Vent System assisted injector where the injector was cooled by the cold gas generated by Joule-Thompson expansion of the liquid propellant. The numerical model accounts for a) different regimes of pool boiling heat transfer, b) condensation of vapor around spray droplets, and c) condensation of vapor at the interface of the cooled injector and ullage. The numerical predictions of tank pressure, filling rate, and wall temperatures were compared with the test data and one-dimensional (1D) model predictions. It was demonstrated that the axisymmetric 2D model provided a more realistic prediction than the 1D model.
Two-Dimensional Network Flow Modeling of No-Vent Tank Filling of A Cryogenic Tank With Thermo-Dynamic Vent System Assisted Injector
A finite volume-based network flow modeling tool, Generalized Fluid System Simulation Program (GFSSP), was used to simulate the no-vent filling of a cryogenic tank with a thermodynamic vent system assisted injector. The tank was discretized into nodes in the radial and axial directions to calculate two-dimensional axisymmetric flow with heat transfer between the solid wall and fluid. The nodes are connected by branches where momentum equations are solved to calculate the flowrates. The mass and energy conservation equations and the equation of state are solved to calculate pressure, temperature, and resident mass at the nodes. The energy conservation equations at the solid nodes were solved to calculate the temperature of the solid wall. The system of equations was solved by a combination of successive substitution and simultaneous Newton-Raphson method. The numerical model accounts for a) different regimes of pool boiling heat transfer, b) condensation of vapor around the spray droplets, c) condensation of vapor at the interface of the cooled injector and ullage, and d) condensation at the liquid-vapor interphase. The numerical predictions of tank pressure, filling rate, and wall temperatures were compared with the test data. The paper also presents the effect of nodal discretization by comparing the results of the 1D and 2D models and examines the effect of the droplet diameter on tank filling.
Two-Dimensional Network Flow Modeling of No-Vent Tank Filling of A Cryogenic Tank With Thermo-Dynamic Vent System Assisted Injector
A finite volume-based network flow modeling tool, Generalized Fluid System Simulation Program (GFSSP), was used to simulate the no-vent filling of a cryogenic tank with a thermodynamic vent system assisted injector. The tank was discretized into nodes in the radial and axial directions to calculate two-dimensional axisymmetric flow with heat transfer between the solid wall and fluid. The nodes are connected by branches where momentum equations are solved to calculate the flowrates. The mass and energy conservation equations and the equation of state are solved to calculate pressure, temperature, and resident mass at the nodes. The energy conservation equations at the solid nodes were solved to calculate the temperature of the solid wall. The system of equations was solved by a combination of successive substitution and simultaneous Newton-Raphson method. The numerical model accounts for a) different regimes of pool boiling heat transfer, b) condensation of vapor around the spray droplets, c) condensation of vapor at the interface of the cooled injector and ullage, and d) condensation at the liquid-vapor interphase. The numerical predictions of tank pressure, filling rate, and wall temperatures were compared with the test data. The paper also presents the effect of nodal discretization by comparing the results of the 1D and 2D models and examines the effect of the droplet diameter on tank filling.
Collaborative investigation of the internal flow and near-nozzle flow of an eight-hole gasoline injector (Engine Combustion Network Spray G)
The internal details of fuel injectors have a profound impact on the emissions from gasoline direct injection engines. However, the impact of injector design features is not currently understood, due to the difficulty in observing and modeling internal injector flows. Gasoline direct injection flows involve moving geometry, flash boiling, and high levels of turbulent two-phase mixing. In order to better simulate these injectors, five different modeling approaches have been employed to study the engine combustion network Spray G injector. Here these simulation results have been compared to experimental measurements obtained, among other techniques, with X-ray diagnostics, allowing the predictions to be evaluated and critiqued. The ability of the models to predict mass flow rate through the injector is confirmed, but other features of the predictions vary in their accuracy. The prediction of plume width and fuel mass distribution varies widely, with volume-of-fluid tending to overly concentrate the fuel. All the simulations, however, seem to struggle with predicting fuel dispersion and by inference, jet velocity. This shortcoming of the predictions suggests a need to improve Eulerian modeling of dense fuel jets.
Modeling and Simulation of Tank Pressure Control using Zero-Boiloff 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.
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.