Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “two-phase flows”

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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
NASA NTRS1998

Selection and Assessment of Predictions of the Mars Pathfinder Landing Site

The successful landing of the Mars Pathfinder spacecraft on Mars allows the review of the process of selecting the landing site and assessing predictions made for the site based on Viking and Earth-based data. Selection of the landing site for Mars Pathfinder was a two-phase process. The first phase took place from October 1993 to June 1994 and involved: initial identification of engineering constraints, definition of environmental conditions at the site for spacecraft design, and evaluation of the scientific potential of different landing sites. This phase culminated with the first "Mars Pathfinder Landing Site Workshop", held at the Lunar and Planetary Institute in Houston, Texas on April 18-19, 1994, in which suggested approaches and landing sites were solicited from the entire scientific community. A preliminary site was selected by the project for design purposes in June 1994. The second phase took place from July 1994 to March 1996 and involved: developing criteria for evaluating site safety using images and remote sensing data, testing of the spacecraft and landing subsystems (with design improvements) to establish quantitative engineering constraints on landing site characteristics, evaluating all potential landing sites on Mars, and certification of the site by the project. This phase included a second open workshop, "Mars Pathfinder Landing Site Workshop II: Characteristics of the Ares Vallis Region and Field Trips in the Channeled Scabland, Washington" held in Spokane and Moses Lake September 24-30, 1995 and formal acceptance of the site by NASA Headquarters. Engineering constraints on Pathfinder landing sites were developed from the initial design of the spacecraft and the entry, descent and landing scenario. The site must be within 5 degrees of the subsolar latitude at the time of landing (15N for maximum solar power and flexible communications with Earth. It also must be below 0 km elevation to enable enough time for the parachute to bring the lander to the proper terminal velocity for landing. The entire landing ellipse, which is 70 km by 200 km due to navigational, ephemeris and atmospheric uncertainties, must be free of steep slopes, scarps and obvious hazards in Viking orbiter images, have acceptable radar reflectivity, moderate rock abundances and have little or no dust. Scientific considerations of the Mars Pathfinder payload and mission indicate that analyses of "grab bag" samples at the mouths of outflow channels can offer a first order assessment of a variety of rock types on Mars. Highland sites offer the advantage of in situ analysis of ancient rocks on Mars that record crustal differentiation and the nature of the early environment. Dark gray sites offer the potential of analyzing unweathered and unoxidized materials. Following a general assessment of the safety of different sites, a preliminary selection of a "grab bag" site was made. This site, Ares Vallis, is near the mouth of an outflow channel that may contain ancient Noachian terrain, Hesperian ridged plains, and reworked channel materials. All potential landing sites on Mars that met basic safety criteria were analyzed in detail. Sites (100 by 200 km target ellipses) were considered safe if they were below 0 km elevation, were free of obvious hazards (high relief surface features) in high-resolution (< 50 m/pixel) Viking orbiter images and had acceptable reflectivity and roughness at radar wavelengths, high thermal inertia, moderate rock abundance, low red to violet ratio, and low albedo. Only 4 sites on Mars met all the above criteria, which included 1995 opposition 3.5 cm delay-Doppler radar data. Complete data were evaluated for 7 sites and the Viking landing sites for comparison for all the above criteria as well as crater abundance, hill and mesa abundance, slopes over meter to kilometer scales, low altitude winds (from global circulation models and slopes), the size-frequency distribution of large rocks, as well as rover trafficability and science potential. Discussion of potential hazards at Ares Vallis using a variety of data sets (including radar) at a second open workshop, indicated this site cannot be shown to be any more hazardous than the Viking landing sites. Field trips to the Channeled Scabland and the Ephrata Fan, analogs for Ares Vallis and the landing site, respectively, provided valuable insight into possible geologic processes and potential surface characteristics. Three sites met all the data requirements and safety criteria for landing Pathfinder. Ares Vallis was selected by the project because it appeared acceptably safe (although it appeared to have greater rock abundances than other sites, its elevation was likely the best known) and offered the prospect of analyzing a variety of rock types expected to be deposited by catastrophic floods, which would enable addressing first-order scientific questions such as differentiation of the crust, the development of weathering products, and the nature of the early martian environment and its subsequent evolution. The selection was reviewed by an external board at a number of meetings and accepted, and the site was approved by NASA Headquarters. Data gathered by the Pathfinder lander' and rover provides the opportunity to test the predictions made for the site in the selection process based on remote observations from Earth, orbit, and the surface. The discussion below is taken from Golombek et al. to which the reader is referred for a more complete discussion and a complete list of references, which are omitted here for brevity. Many characteristics of the landing site are consistent with its being shaped and deposited by the Ares and Tiu catastrophic floods. The rocky surface is consistent a depositional plain comprising semi-rounded pebbles, cobbles and tabular boulders (some of which appear imbricated and/or inclined in the direction of flow) that appear similar to depositional plains in terrestrial catastrophic floods. The Twin Peaks appear to be streamlined hills in lander images, which is consistent with interpretations of larger hills in Viking orbiter images of the region that suggest the lander is on the flank of a broad, gentle ridge trending northeast from Twin Peaks. This ridge, which is the rise to the north of the lander, is aligned in the downstream direction from the Ares and Tiu Valles floods, and may be a debris tail deposited in the wake of the Twin Peaks. Channels visible throughout the scene may be a result of late stage drainage. As predicted by delay-Doppler radar measurements and tracking results, the average elevation of the center of the site was about the same as Viking Lander I relative to the 6.1 mbar geoid. The Doppler tracking and two-way ranging estimate for the elevation of the spacecraft is only 45 in lower than the Viking I Lander and within 100 in of that expected, which is within the uncertainties of the measurements. After landing, surface pressures and winds (5-10 m/s) were found to be similar to expectations based on Viking data, although temperatures were about 10 K warmer. The temperature profile below 50 km was also roughly 20 K warmer. As a result, predicted densities were 5% higher near the surface and up to 40% lower at 50 km but within the entry, descent and landing design margins. The populations of craters and small hills and the slopes of the hills measured in high-resolution (38 m/pixel) Viking orbiter images and the radar derived slopes of the landing site are all consistent with observations of these properties in the lander images. A rocky surface was expected from Viking Infra-Red Thermal Mapper (IRTM) observations and comparisons with the Viking landing sites. The observed cumulative fraction of area covered by rocks with diameters greater than 3 cm and heights greater than 0.5 in (potentially hazardous to landing) at Ares is similar to that predicted by IRTM observations and models of Viking lander and Earth analog rock size-frequency distributions. The IRTM prediction postulated an effective thermal inertia of 30 (10(exp -3) cgs units - cal/cubic cm/s(exp 0.5)/K) for the rock population, but we obtain a slightly different effective thermal inertia for the actual rock population. The validity of interpretations of radar echoes prior to landing are supported by a simple radar echo model, an estimate of the reflectivity of the soil from its bulk density, and the fraction of area covered by rocks. In the calculations, the soil produces the quasi-specular echo and the rocks produce the diffuse echo. The derived quasispecular cross section is comparable to the cross-sections and reflectivities reported for 3.5-cm wavelength observations. The model yields a diffuse echo that is modestly larger than the polarized diffuse echo reported for 3.5-cm wavelength observations. At 12.5-cm wavelength, similar rock populations at Ares and the Viking I site were expected because the diffuse echoes are comparable, but the large normal reflectivities suggests that bulk densities of the soils at depth are greater than those at the surface. We also obtain a fine-component inertia near 8.4 which agrees with the fine-component inertia of 8.7 (in 10(exp -3) cgs units) estimated from thermal observations from orbit by the IRTM; for this estimate, we used a bulk thermal inertia of 10.4 for the landing site, an effective thermal inertia near 40 (10(exp -3) cgs units) for the rock population, and a graphical representation of Kieffer's model. Color and albedo data for Ares suggested surfaces of materials at Ares Vallis would be relatively dust free or unweathered prior to landing compared with the materials at the Viking landing sites. This suggestion is supported by the abundance of relatively dark-gray rocks at Ares and their relative rarity at the Viking landing sites, where rocks are commonly coated with bright red dust. Finally, the 40 km long Ephrata Fan of the Channeled Scabland in Washington state, which was deposited where c

Golombek, M. P.↗
NASA NTRS2000

Computations of Boiling in Microgravity

The absence (or reduction) of gravity, can lead to major changes in boiling heat transfer. On Earth, convection has a major effect on the heat distribution ahead of an evaporation front, and buoyancy determines the motion of the growing bubbles. In microgravity, convection and buoyancy are absent or greatly reduced and the dynamics of the growing vapor bubbles can change in a fundamental way. In particular, the lack of redistribution of heat can lead to a large superheat and explosive growth of bubbles once they form. While considerable efforts have been devoted to examining boiling experimentally, including the effect of microgravity, theoretical and computational work have been limited. Here, the growth of boiling bubbles is studied by direct numerical simulations where the flow field is fully resolved and the effects of inertia, viscosity, surface deformation, heat conduction and convection, as well as the phase change, are fully accounted for. Boiling involves both fluid flow and heat transfer and thus requires the solution of the Navier-Stokes and the energy equations. The numerical method is based on writing one set of governing transport equations which is valid in both the liquid and vapor phases. This local, single-field formulation incorporates the effect of the interface in the governing equations as source terms acting only at the interface. These sources account for surface tension and latent heat in the equations for conservation of momentum and energy as well as mass transfer across the interface due to phase change. The single-field formulation naturally incorporates the correct mass, momentum and energy balances across the interface. Integration of the conservation equations across the interface directly yields the jump conditions derived in the local instant formulation for two-phase systems. In the numerical implementation, the conservation equations for the whole computational domain (both vapor and liquid) are solved using a stationary grid and the phase boundary is followed by a moving unstructured two-dimensional grid. While two-dimensional simulations have been used for preliminary studies and to examine the resolution requirement, the focus is on fully three-dimensional simulations. The numerical methodology, including the parallelization and grid refinement strategy is discussed, and preliminary results shown. For buoyancy driven flow, the heat transfer is in good agreement with experimental correlations. The changes when gravity is turned off and/or fluid shear is added are discussed, as well as the difference between simulations of a layer freely releasing bubbles versus simulations using only one wavelength initial perturbation. Figure 1 shows the early stages of the formation of a three-dimensional bubble from a thin vapor layer. The boundary conditions are periodic in the x and y direction, the bottom is a hot and the top allows a free outflow. The jagged edge of the surface close to the bottom of the computational domain is due to some of the surface elements being on the other side of the domain and some elements not plotted by our plotting routine. In the second figure, we show the temperature distribution through two perpendicular planes.

Tryggvason, G.↗
NASA NTRS2001

Robust Strategy for Rocket Engine Health Monitoring

Monitoring the health of rocket engine systems is essentially a two-phase process. The acquisition phase involves sensing physical conditions at selected locations, converting physical inputs to electrical signals, conditioning the signals as appropriate to establish scale or filter interference, and recording results in a form that is easy to interpret. The inference phase involves analysis of results from the acquisition phase, comparison of analysis results to established health measures, and assessment of health indications. A variety of analytical tools may be employed in the inference phase of health monitoring. These tools can be separated into three broad categories: statistical, rule based, and model based. Statistical methods can provide excellent comparative measures of engine operating health. They require well-characterized data from an ensemble of "typical" engines, or "golden" data from a specific test assumed to define the operating norm in order to establish reliable comparative measures. Statistical methods are generally suitable for real-time health monitoring because they do not deal with the physical complexities of engine operation. The utility of statistical methods in rocket engine health monitoring is hindered by practical limits on the quantity and quality of available data. This is due to the difficulty and high cost of data acquisition, the limited number of available test engines, and the problem of simulating flight conditions in ground test facilities. In addition, statistical methods incur a penalty for disregarding flow complexity and are therefore limited in their ability to define performance shift causality. Rule based methods infer the health state of the engine system based on comparison of individual measurements or combinations of measurements with defined health norms or rules. This does not mean that rule based methods are necessarily simple. Although binary yes-no health assessment can sometimes be established by relatively simple rules, the causality assignment needed for refined health monitoring often requires an exceptionally complex rule base involving complicated logical maps. Structuring the rule system to be clear and unambiguous can be difficult, and the expert input required to maintain a large logic network and associated rule base can be prohibitive.

Santi, L. Michael↗
NASA NTRS2006

Evaporating Spray in Supersonic Streams Including Turbulence Effects

Evaporating spray plays an important role in spray combustion processes. This paper describes the development of a new finite-conductivity evaporation model, based on the two-temperature film theory, for two-phase numerical simulation using Eulerian-Lagrangian method. The model is a natural extension of the T-blob/T-TAB atomization/spray model which supplies the turbulence characteristics for estimating effective thermal diffusivity within the droplet phase. Both one-way and two-way coupled calculations were performed to investigate the performance of this model. Validation results indicate the superiority of the finite-conductivity model in low speed parallel flow evaporating sprays. High speed cross flow spray results indicate the effectiveness of the T-blob/T-TAB model and point to the needed improvements in high speed evaporating spray modeling.

Balasubramanyam, M. S.↗
NASA NTRS2004

Freeze-Tolerant Condensers

Two condensers designed for use in dissipating heat carried by working fluids feature two-phase, self-adjusting configurations such that their working lengths automatically vary to suit their input power levels and/or heat-sink temperatures. A key advantage of these condensers is that they can function even if the temperatures of their heat sinks fall below the freezing temperatures of their working fluids and the fluids freeze. The condensers can even be restarted from the frozen condition. The top part of the figure depicts the layout of the first condenser. A two-phase (liquid and vapor) condenser/vapor tube is thermally connected to a heat sink typically, a radiatively or convectively cooled metal panel. A single-phase (liquid) condensate-return tube (return artery) is also thermally connected to the heat sink. At intervals along their lengths, the condenser/vapor tube and the return artery are interconnected through porous plugs. This condenser configuration affords tolerance of freezing, variable effective thermal conductance (such that the return temperature remains nearly constant, independently of the ultimate sink temperature), and overall pressure drop smaller than it would be without the porous interconnections. An additional benefit of this configuration is that the condenser can be made to recover from the completely frozen condition either without using heaters, or else with the help of heaters much smaller than would otherwise be needed. The second condenser affords the same advantages and is based on a similar principle, but it has a different configuration that affords improved flow of working fluid, simplified construction, reduced weight, and faster recovery from a frozen condition.

Crowley, Christopher J.↗
NASA NTRS2011

Analytical Approach in DeCoM

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

Patel, Deepak↗
NASA NTRS2013

Effect of Residual Noncondensables on Pressurization and Pressure Control of a Zero-Boil-Off Tank in Microgravity

The Zero-Boil-Off Tank (ZBOT) Experiment is a small-scale experiment that uses a transparent ventless Dewar and a transparent simulant phase-change fluid to study sealed tank pressurization and pressure control with applications to on-surface and in-orbit storage of propellant cryogens. The experiment will be carried out under microgravity conditions aboard the International Space Station in the 2014 timeframe. This paper presents preliminary results from ZBOT's ground-based research that focuses on the effects of residual noncondensable gases in the ullage on both pressurization and pressure reduction trends in the sealed Dewar. Tank pressurization is accomplished through heating of the test cell wall in the wetted and un-wetted regions simultaneously or separately. Pressure control is established through mixing and destratification of the bulk liquid using a temperature controlled forced jet flow with different degrees of liquid jet subcooling. A Two-Dimensional axisymmetric two-phase CFD model for tank pressurization and pressure control is also presented. Numerical prediction of the model are compared to experimental 1g results to both validate the model and also indicate the effect of the noncondensable gas on evolution of pressure and temperature distributions in the ullage during pressurization and pressure control. Microgravity simulations case studies are also performed using the validated model to underscore and delineate the profound effect of the noncondensables on condensation rates and interfacial temperature distributions with serious implications for tank pressure control in reduced gravity.

Marangoni and Natural Convection↗
NASA NTRS2014

On the TFNS Subgrid Models for Liquid-Fueled Turbulent Combustion

This paper describes the time-filtered Navier-Stokes (TFNS) approach capable of capturing unsteady flow structures important for turbulent mixing in the combustion chamber and two different subgrid models used to emulate the major processes occurring in the turbulence-chemistry interaction. These two subgrid models are termed as LEM-like model and EUPDF-like model (Eulerian probability density function), respectively. Two-phase turbulent combustion in a single-element lean-direct-injection (LDI) combustor is calculated by employing the TFNS/LEM-like approach as well as the TFNS/EUPDF-like approach. Results obtained from the TFNS approach employing these two different subgrid models are compared with each other, along with the experimental data, followed by more detailed comparison between the results of an updated calculation using the TFNS/LEM-like model and the experimental data.

Liu, Nan-Suey↗
NASA NTRS2015

Investigations of Physical Processes in Microgravity Relevant to Space Electrochemical Power Systems

NASA has performed physical science microgravity flight experiments in the areas of combustion science, fluid physics, material science and fundamental physics research on the International Space Station (ISS) since 2001. The orbital conditions on the ISS provide an environment where gravity driven phenomena, such as buoyant convection, are nearly negligible. Gravity strongly affects fluid behavior by creating forces that drive motion, shape phase boundaries and compress gases. The need for a better understanding of fluid physics has created a vigorous, multidisciplinary research community whose ongoing vitality is marked by the continuous emergence of new fields in both basic and applied science. In particular, the low-gravity environment offers a unique opportunity for the study of fluid physics and transport phenomena that are very relevant to management of fluid - gas separations in fuel cell and electrolysis systems. Experiments conducted in space have yielded rich results. These results provided valuable insights into fundamental fluid and gas phase behavior that apply to space environments and could not be observed in Earth-based labs. As an example, recent capillary flow results have discovered both an unexpected sensitivity to symmetric geometries associated with fluid container shape, and identified key regime maps for design of corner or wedge-shaped passive gas-liquid phase separators. In this presentation we will also briefly review some of physical science related to flight experiments, such as boiling, that have applicability to electrochemical systems, along with ground-based (drop tower, low gravity aircraft) microgravity electrochemical research. These same buoyancy and interfacial phenomena effects will apply to electrochemical power and energy storage systems that perform two-phase separation, such as water-oxygen separation in life support electrolysis, and primary space power generation devices such as passive primary fuel cell.

energy storage↗
NASA NTRS2018

Validation of a CFD Model Predicting the Effect of High Level Lateral Acceleration Sloshing on the Heat Transfer and Pressure Drop in a Small-Scale Tank in Normal Gravity

A two-phase CFD model is developed to study the effects of sloshing with high level lateral acceleration on the heat transfer and pressure drop in a small scale tank. Computational results are compared to the data provided by a non-isothermal sloshing experiment without phase change conducted by T. Himeno et al. at the University of Tokyo and JAXA in 2011. The results of the current model are, also, compared to CFD predictions reported by Himeno et al. A step change in lateral acceleration was applied in the experiment. Different levels of lateral acceleration amplitude, varying between 0G and 0.5G, were considered. CFD results for interface movement and tank pressure are presented and compared in this paper to the experimental data for the case in which the value of lateral acceleration was set to 0.5G. The effects of initial and boundary conditions and turbulence modeling approach on the tank pressure change during sloshing are discussed in detail. The effect of conjugate heat transfer in the tank wall is also studied to show its important role in determining the tank pressure evolution. The results of the Reynolds Averaged Navier Stokes (RANS) models are compared to the results of the Large Eddy Simulation model (LES) to underscore the importance of correctly capturing the effects of turbulence for high fidelity predictions.

multiphase flow↗
NASA NTRS2018

Validation of A CFD Model Predicting the Effect of High Level Lateral Acceleration Sloshing on the Heat Transfer and Pressure Drop in A Small-Scale Tank in Normal Gravity

A two-phase CFD model is developed to study the effects of sloshing with high level lateral acceleration on the heat transfer and pressure drop in a small scale tank. Computational results are compared to the data provided by a non-isothermal sloshing experiment without phase change conducted by T. Himeno et al. at the University of Tokyo and JAXA in 2011 [1]. The results of the current model are, also, compared to CFD predictions reported by Himeno et al. [2]. A step change in lateral acceleration was applied in the experiment. Different levels of lateral acceleration amplitude, varying between 0G and 0.5G, were considered. CFD results for interface movement and tank pressure are presented and compared in this paper to the experimental data for the case in which the value of lateral acceleration was set to 0.5G. The effects of initial and boundary conditions and turbulence modeling approach on the tank pressure change during sloshing are discussed in detail. The effect of conjugate heat transfer in the tank wall is also studied to show its important role in determining the tank pressure evolution. The results of the Reynolds Averaged Navier Stokes (RANS) models are compared to the results of the Large Eddy Simulation model (LES) to underscore the importance of correctly capturing the effects of turbulence for high fidelity predictions.

multiphase flow↗
NASA NTRS2014

A Concept of Operations for Far-Term Surface Trajectory-Based Operations (STBO)

The goal of this far-term STBO (Surface Trajectory-Based Operations) ConOps (Concept of Operations) is to increase the efficiency and predictability of airport surface operations, and reduce the environmental impact, by incorporating a time-based component to surface operations. In the far-term NextGen timeframe, airport surface operations will transition from current-day first-come, first-served operations, to strategically scheduled operations in which pilots are recruited as active participants in meeting the precise time-based goals of NextGen surface operations. The far-term STBO concept includes two-phases. Phase 1 introduces time-based traffic flow constraint points, which divide the taxi route into segments with an assigned Required Time of Arrival (RTA). This Phase 1 approach provides temporal certainty only near the traffic flow constraint points, but not in between. Minimal augmentations to the flight deck are required to support required time of arrival (RTA) management. Phase 2 further increases precision and efficiency by introducing full four-dimensional (4D) trajectories, with an x-y location for all times t. This phase assumes adoption of advanced flight deck equipage enabling higher temporal precision sufficient to support aircraft conformance to 4D trajectories. This allows more precision and less temporal uncertainty at all times along the route.

Surface Trajectory-Based Operations (STBO)↗
NASA NTRS2019

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

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

critical heat flux↗
NASA NTRSDate not supplied

Modeling Autogenous Pressurization and Draining of a Cryogenic Storage Tank in Normal Gravity

A two-phase CFD model for autogenous pressurization and draining of a cryogenic storage tank is presented using both the Sharp Interface and Volume-Of-Fluid (VOF) approaches for capturing the front and the associated interfacial heat, mass and momentum transfer between the liquid and vapor regions. Both models are validated against data provided by the Cryogenic Propellant Storage and Transfer (CPST) Engineering Development Unit (EDU) experiment. The results of the autogenous pressurization are presented first, focusing on the phase change and turbulence effects on the tank pressure and temperature predictions. Both the Sharp Interface (SI-CFD) and VOF (VOF-CFD) multiphase models predict tank pressure during pressurization within 3% of the measured values. The sensitivity of key physical and numerical parameters of the problem are tested using the Sharp Interface model. Effects of the accommodation coefficient (AC), and the computational grid structure are studied. The second part of this paper is devoted to validating the SI-VOF model, with an enhanced capability of moving the liquid-vapor interface, against the draining segment of the EDU experiment. The VOF-CFD model was also used to simulate tank draining and its results were compared with the results of the Sharp Interface model with the moving interface. Both models predict tank pressure during draining within 3.5% of the measured values. Pressure decrease rate is underpredicted by both models during the first 100 seconds of draining but matches the experimental rate for the rest of the simulation.

Computational Fluid Dynamics↗
NASA NTRSDate not supplied

Modeling Autogenous Pressurization and Draining of a Cryogenic Storage Tank in Normal Gravity

A two-phase CFD model for autogenous pressurization and draining of a cryogenic storage tank is presented using both the Sharp Interface and Volume-Of-Fluid (VOF) approaches for capturing the front and the associated interfacial heat, mass and momentum transfer between the liquid and vapor regions. Both models are validated against data provided by the Cryogenic Propellant Storage and Transfer (CPST) Engineering Development Unit (EDU) experiment1. The results of the autogenous pressurization are presented first, focusing on the phase change and turbulence effects on the tank pressure and temperature predictions. Both the Sharp Interface (SI-CFD) and VOF (VOF-CFD) multiphase models predict tank pressure during pressurization within 3% of the measured values. The sensitivity of key physical and numerical parameters of the problem are tested using the Sharp Interface model. Effects of the accommodation coefficient (AC), and the computational grid structure are studied. The second part of this paper is devoted to validating the SI-VOF model, with an enhanced capability of moving the liquid-vapor interface, against the draining segment of the EDU experiment. The VOF-CFD model was also used to simulate tank draining and its results were compared with the results of the Sharp Interface model with the moving interface. Both models predict tank pressure during draining within 3.5% of the measured values. Pressure decrease rate is underpredicted by both models during the first 100 seconds of draining but matches the experimental rate for the rest of the simulation.

Computational Fluid Dynamics↗
NASA NTRSDate not supplied

Development and Validation of Two-Phase CFD Models for Key Elements of Propellant Tank CFM Operations in 1G and Microgravity –An Overview

This paper presents an overview of the state-of-the-art two-phase CFD models that have been developed for various propellant tank storage and transfer operations in 1g, partial gravity, and microgravity. The models are developed in the framework of the industry standard ANSYS/Fluent CFD code, the capabilities of which have been significantly enhanced and customized through the incorporation of unique submodels via User Defined Functions (UDF)s to satisfy the requirements of its intended variable gravity Cryogenic Fluid Management (CFM) applications. These models/submodels have been validated against experimental data that cross spatial scales, fluid types, and gravity levels in order to properly anchor the models’ physical and numerical fidelity. The CFM application/processes that have been modeled include tank self-pressurization, tank autogenous pressurization, tank pressure control using both subcooled jet mixing and droplet spray injection mechanisms, tank chilldown/filling, tank drainage, tank slosh for both volatile and non-volatile fluids. The strength and shortcomings of the two-phase models for each application are highlighted and discussed briefly.

Computational Fluid Dynamics↗
NASA NTRS2023

Universal Two-Phase Convection Heat Transfer Correlations for Cryogenic Pipe Chilldown

This paper presents a set of universal two-phase convection heat transfer correlations for modeling boiling heat transfer during pipe chilldown fit over the widest available range of cryogenic fluids and thermodynamic conditions. The correlations improve upon prior correlations that were developed separately for liquid nitrogen (LN2) and liquid hydrogen (LH2) pipe quenching datasets. The new correlations include equations to calculate the single-phase vapor heat transfer, film boiling heat transfer, transition boiling heat transfer, nucleate boiling heat transfer, single-phase liquid heat transfer, bulk vapor temperature during high quality film boiling, Leidenfrost temperature, critical heat flux, critical heat flux temperature, and the onset of nucleate boiling temperature. The correlations were validated against LH2, LN2, liquid methane, liquid oxygen, and liquid argon pipe quenching datasets. The eight datasets cover the following parameter ranges: pipe lengths of 0.1 to 6.5m; outer pipe diameters of 12.7 to 25.4 mm; pipe wall thicknesses from 0.51 to 1.64 mm; flow directions of upward, downward, and horizontal; gravity levels of 1g and 0g±0.01g. A numerical model assumes homogeneous mixing between the vapor and liquid and implicitly integrates the coupled energy equations for the pipe and fluid, as well as continuity for the fluid. The model was used to estimate the equilibrium quality and fluid mass accumulation along the pipes during chilldown from estimates of the fluid-pipe heat transfer extracted from temperature measurements. The correlations can be implemented in lumped parameter codes such as SINDA/FLUINT and the Generalized Fluid System Simulation Program (GFSSP) to improve the accuracy of chilldown time and chilldown boiloff mass predictions. Such predictions are useful for designing ground or in-space cryogenic liquid transfer systems.

Universal Correlations↗