Search NASA⌕ Search

SEARCH · Search NASA

Results for “Two-phase heat transfer coefficient”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS2012

Measurement of Two-Phase Flow and Heat Transfer Parameters using Infrared Thermometry

A novel technique to measure heat transfer and liquid film thickness distributions over relatively large areas for two-phase flow and heat transfer phenomena using infrared (IR)thermometry is described. IR thermometry is an established technology that can be used to measure temperatures when optical access to the surface is available in the wavelengths of interest. In this work, a midwave IR camera (3.6-5.1 microns) is used to determine the temperature distribution within a multilayer consisting of a silicon substrate coated with a thin insulator. Since silicon is largely transparent to IR radiation, the temperature of the inner and outer walls of the multilayer can be measured by coating selected areas with a thin, IR opaque film. If the fluid used is also partially transparent to IR, the flow can be visualized and the liquid film thickness can be measured. The theoretical basis for the technique is given along with a description of the test apparatus and data reduction procedure. The technique is demonstrated by determining the heat transfer coefficient distributions produced by droplet evaporation and flow boiling heat transfer.

Kim, Tae-Hoon↗
NASA NTRSDate not supplied

Analysis of Heat Transfer from Local Heating and Cooling Sources at Cryogenic Temperatures

Understanding the dispersion of heat around a cryogenic fluid tank, specifically the interaction between the cryogenic fluid and the tank wall is critical in the analysis of long duration cryogen storage in microgravity. The heat transfer interaction between a cryogenic storage tank and heat sources from external spacecraft structures is also one of the many factors that determine how much heat enters a tank. Recent flight experiments with two-phase fluids have indicated that local concentrations of heat input (also known as “hot spots”) can cause unwanted affects including local boiling. Computational fluid dynamic (CFD) models can provide a detailed assessment of the heat transfer occurring across a cryogenic storage system. However, CFD modeling takes time to construct and run. A simpler approach that can act as initial guidance for later CFD modeling analyzes external “hot spots” as point or finite heat sources. A radial, finite element network or a local direct solution can effectively estimate the heat spread across a cryogenic storage tank by calculating the temperature and heat load as a function of distance from the heat source. This calculation accounts for the convective heat transfer between the cryogenic fluid and storage tank surface. Similar approaches can be used to determine the effectiveness of cooling from a cryocooler as a finite, local heat sink. This approach allows for quick approximations of the thermal map across a cryogenic tank as well as sensitivity analysis under a wide range of design parameters including gravitational fields as implied through natural convection coefficients.

Cryogenic Fluid Management↗
DOE OSTI2025

MOSCATO Development and Integration in Fiscal Year 2025: Implementation of Multiphase, Multiphysics Modeling Capabilities for Molten Salt Systems

MOSCATO (Molten Salt Chemistry and Transport) is a multiphysics code that provides high-fidelity, coupled simulations of fluid flow, heat transfer, mass transfer, chemistry, electrochemical phenomena, and alloy corrosion for molten salt systems. In FY25, significant developments were made to the code package, enhancing its capabilities for modeling all relevant phenomena within operating moltens salt reactors (MSRs). The developments and activities in FY25 included: 1. Implementation of Level-Set methods to enable modeling of single-bubble behavior in molten salts. In FY25, the Level-Set two-phase flow modeling implementation was improved to simulate single bubble behavior with molten salt media. The large density and viscosity ratios between typical gases and molten salt liquids present challenges for these types of numerical solvers. With enhancements to the pressure projection method, MOSCATO’s Level-Set solver was able to be successfully validated to experiments related to helium bubble rise in stagnant molten salt. The simulated bubble rising velocity showed reasonable good agreement with experimental measurements. The bubble shape and dynamics were also visually compared with experimental snapshots, demonstrating a good qualitative match. 2. Generation of mass transfer correlations for multiphase flow systems. To enable calculations of the tritium transport across the interface between gas bubbles and salt, we modeled high- Schmidt-number mass transfer around a sphere across a broad range of Reynolds numbers. The mesh near the sphere surface was highly refined to resolve steep concentration gradients caused by the low diffusion coefficient. Literature-based mass transfer correlations were compared with the numerical results, and modifications were proposed to improve agreement, particularly at higher Schmidt numbers. These mass transfer correlations were subsequently provided to other national laboratories to help enable high quality mass transfer simulations using lower-order solvers under development within the NEAMS program. 3. Preliminary implementation of a bubbly flow solver. To model bubbly flow in molten salt, we implemented a bubbly flow solver for void fractions less than 5%. To do so, an algebraic relative velocity model that assumes small bubbles with rapid momentum equilibration was added to MOSCATO to compute bubble velocities. Preliminary comparisons with experimental data showed reasonable agreement, and further improvements are underway. 4. Generation of mass transfer correlations for MSRE subchannel The Molten-Salt Reactor Experiment (MSRE) was a landmark historical project that demonstrated the feasibility of molten-salt reactor technology. The MSRE campaign also generated a significant body of experimental data and reports that continue to support molten-salt–related research. In this report, the MSRE core subchannel was used as the reference geometry for a mass transfer study performed with MOSCATO. The geometry and computational mesh were adapted from a previous study, providing adequate resolution for the relatively low Reynolds number in this case. Additional mesh refinement was applied to reach higher Schmidt numbers, enabling the derivation of a reliable mass-transfer correlation for the present scenario. 5. Simulations of oxygen ingressions into molten salt. In the previous fiscal year, we initiated a study to simulate oxygen ingression in stagnant salt. As oxygen enters the salt through its surface, it reacts with Ce 3+ to form solid CeO 2 and other reaction products. To more fully capture the complex diffusion-convection-reaction mechanisms, capabilities for modeling natural convection in the salt vessel were added. These were needed as the flow of the ingressed gas induced flow in the salt caused by surface shear and non-isothermal effects. With these updated physics in place, we were able to successfully reproduce the experimental results for the rate of change of CeCl 3 concentrations versus time. 6. Flow corrosion model validation. In FY24, MOSCATO’s corrosion model was validated against static corrosion experiments. In FY25, this work was extended to a flow corrosion experiment, where FLiNaK salt was driven by natural convection, with initial salt impurities to initiate corrosion. Despite uncertainties in parameters such as elemental diffusion coefficients in the alloy and unknown H + concentrations, the simulations achieved good agreement with experimental results, especially in predicting sample mass losses.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Non-dimensional performance and safety parameters for heat pipes

The use of heat pipes in safety-critical systems such as nuclear microreactors dictates the development of generalized, practical, scalable performance and safety parameters. Traditional dimensional metrics, while informative, lack the universality required for comparative analysis across varying designs and operating regimes. Here, this work introduces a comprehensive set of non-dimensional parameters to characterize heat pipe performance and safety, including capillary performance, effective thermal conductivity, response time, exergetic efficiency, allowable temperature gradients, allowable rate of temperature change, priming coefficients, and factor of safety. A reference heat pipe design representative of microreactor applications was analyzed via the developed parameters using both traditional analytical models and Sockeye simulations under transient and steady-state conditions. Sodium, potassium, and water were evaluated as working fluids to demonstrate the applicability of the framework across a broad temperature range. The proposed non-dimensional parameters effectively captured key thermal-hydraulic behaviors and safety concerns, as was demonstrated via Sockeye simulations. This framework supports the development of design optimization strategies, operational protocols, and safety assurance practices for advanced reactor systems and other high-reliability applications.

42 - ENGINEERING↗
DOE OSTI2025

Integrated Molten Salt Reactor Modeling Capabilities in NEAMS Thermal Hydraulics Tools

The DOE neams program supports a full range of computational thermal fluids analysis capabilities and code developments for a broad range of advanced reactor concepts. The research and development approach under the thermal fluids technical area synergistically combines three length and time scales in a hierarchical multi-scale approach. To enable multi-scale thermal fluids capability using these codes, a key joint effort has been underway to develop an integrated system- and engineering-scale thermal fluids analysis capability, through integration of SAM and Pronghorn codes, both based on the MOOSE framework. This report summarizes recent advances in developing an integrated system- and engineering-scale modeling capability for the msr concept, which has gained significant interest in recent years. A consistent framework was established by coupling Pronghorn and SAM through the Saline interface, with thermophysical properties provided by the Molten Salt Thermal Property Database (MSTDB-TP). Further improvements were made to the coupling schemes and domain-overlapping strategies, enhancing the stability and robustness of multi-code simulations. Verification and validation efforts demonstrate the accuracy of this integration across a range of benchmark problems, including one-dimensional heated pipe flows, three-dimensional natural convection loops with evolving isotopic compositions, and \gls{msre} demonstration cases. Within Pronghorn, new capabilities were introduced to model corrosion and noble-metal plating phenomena, supported by an extended thermal-hydraulics framework and refined turbulence treatments. To capture two-phase flow behavior, a multiphase Euler–Euler model was implemented in Pronghorn, including advanced closure relations, high-resolution advection techniques, and capillary force reconstruction. Preliminary verification cases confirm the fidelity of the approach, while planned validation efforts target canonical multiphase benchmarks and application to msr components such as the msre pump bowl. Finally, updates to SAM’s msr mass transfer modeling were extended to consider noble gas migration into porous structures like graphite. The point kinetics model was updated to include reactivity feedback contributions from any defined species, such as xenon. The gas transport model was expanded for applicability to gas mixtures, bubble efflux phenomena, and species transport between liquid and gas phases. A selection of multi-scale Sherwood number correlations from MOSCATO/NekRS and multi-phase correlations from literature have been added for improved accuracy in calculating mass transfer coefficients. A companion effort on developing system-level redox corrosion has also been incorporated into SAM. Collectively, these enhancements strengthen the predictive capability of SAM and Pronghorn for simulating MSR thermal-hydraulics, corrosion, multiphase behavior, and fission-product transport, providing a more complete toolset for design, safety analysis, and licensing support of next-generation \gls{msr}s.

42 - ENGINEERING↗
DOE OSTI2026

Deflagration to Detonation Transition Update: XDDT Code Modularization

A legacy FORTRAN 77 implementation of the Baer–Nunziato two-phase mixture theory for deflagration-to-detonation transition (DDT) in reactive granular materials—hereafter the XDDT (eXplosive DDT) code—has been modularized to Fortran 90 with modular structure, external input files, and adaptive mesh capability. During validation, two code defects were identified and corrected: an inconsistency in the nodal solid pressure evaluation and a nonphysical burn-front tracking criterion. The ignition criterion was also corrected to use the granular surface temperature from the interface heat transfer model, matching the original Baer implementation. An initial attempt to validate against Figure 3 of the original Baer and Nunziato (1986) paper revealed that the code’s detonation velocity on a 201-node mesh (5.5 km/s) was approximately 21% below the expected Chapman–Jouguet value for 70% TMD HMX (∼7 km/s). Validation was redirected to the piston-driven DDT experiments of McAfee et al. (1989), Shot B-9036, for which well-characterized ionization-pin data are available. With the compaction-burn coefficient calibrated to 𝐶 𝛼 = 75, the XDDT code reproduces the DDT transition time to within 0.4% and produces a steady-state detonation velocity within 4% of the McAfee experimental value of 6.36 km/s. The burn model was generalized to support pressure-dependent exponents, enabling application to nitrocellulose-based ball propellants (TS3659) with a cube-root pressure dependence. Validation against the Sandusky/Baer PDC82 piston-impact experiment yielded a reactive wave velocity of 2.3–2.8 km/s, in good agreement with the experimental value of ∼2.2 km/s, and wave coalescence within 5% of the experimental timing. The mathematical model, input parameter requirements, and a roadmap for extending XDDT to PETN with an autocatalytic burn model are presented.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2025

A Finite Element Method for Compressible and Turbulent Multiphase Flow Instabilities with Heat Transfer

We present a new finite element framework for modeling compressible, turbulent multiphase flows with heat transfer. For two-fluid systems with a free surface, the Volume of Fluid (VOF) method is implemented without the need for interface reconstruction, while turbulence is resolved using a dynamic Vreman large eddy simulation (LES) model. Unlike most two-phase VOF studies, which neglect heat transfer, the present approach incorporates energy transport equations within the VOF formulation to account for heat exchange, an effect particularly important in turbulent flows. Conjugate heat transfer is often challenging in finite volume methods, which require explicit specification of heat fluxes at the solid–fluid interface, limiting accuracy and predictive capability. By contrast, the finite element formulation does not require heat flux inputs, allowing more accurate and robust simulation of heat transfer between solids and fluids. The method is demonstrated through three representative cases. First, a two-fluid instability with a single-mode perturbation is simulated and validated against analytical growth rates. Second, conjugate heat transfer is examined in a high-temperature flow over a cold metal cylinder, with validation performed both quantitatively—via pressure coefficient comparisons with experimental data—and qualitatively using vector field topology. Finally, compressible spray injection and breakup are modeled, demonstrating the ability of the framework to capture interfacial dynamics and atomization under turbulent, high-speed conditions. In the compressible spray injection and breakup case, the results indicate that the finite element formulation achieved higher predictive accuracy and robustness than the finite-volume method. With the same mesh resolution, the FEM reduced the root mean square error (RMSE) and mean absolute percentage error (MAPE) from 6.96 mm and 26.0% (for the FVM) to 4.85 mm and 12.7%, respectively, demonstrating improved accuracy and robustness in capturing interfacial dynamics and heat transfer. The study also introduced vector field topology to visualize and interpret coherent flow structures and instabilities, offering insights beyond conventional scalar-field analyses.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1977

Assessment of analytical techniques for predicting solid propellant exhaust plumes

The calculation of solid propellant exhaust plume flow fields is addressed. Two major areas covered are: (1) the applicability of empirical data currently available to define particle drag coefficients, heat transfer coefficients, mean particle size and particle size distributions, and (2) thermochemical modeling of the gaseous phase of the flow field. Comparisons of experimentally measured and analytically predicted data are made. The experimental data were obtained for subscale solid propellant motors with aluminum loadings of 2, 10 and 15%. Analytical predictions were made using a fully coupled two-phase numerical solution. Data comparisons will be presented for radial distributions at plume axial stations of 5, 12, 16 and 20 diameters.

Tevepaugh, J. A.↗
NASA NTRS1985

A swirl flow evaporative cold plate

A forced flow evaporative cold plate is under development for future application to the thermal bus concept being pursued by NASA for Space Station Thermal Control. The vaporizer is a swirl-flow device employing a spiral tube coil geometry sandwiched between conductive metal plates upon which electric components could be mounted. This concept is based on the inherent phase separation that occurs in a two phase stream in curvilinear flow. This is a zero 'g' design with one 'g' all-attitude capability and is capable of high heat transfer coefficients, good isothermality, and the ability to function at heat fluxes approaching 5w/sq cm on the cold plates (10w/sq cm on the tube wall) with Freon 114. The advantages of this design over other two phase evaporator approaches are high heat flux capability, simplified control requirements, insensitivity to micro-gravity oscillations, and inexpensive manufacturability. The program included design, fabrication, and test of such a cold plate utilizing an existing test stand developed for two-phase thermal management system (TPTMS) testing. Test results analysis and conclusions are included.

Niggemann, R. E.↗
NASA NTRS1990

High heat flux single phase heat exchanger

This paper presents the results obtained to date in a program to develop a high heat flux, single-phase heat exchanger for spacecraft thermal management. The intended application is a net generation interface heat exchanger to couple the crew module water thermal bus to the two-phase ammonia main thermal bus in the Space Station Freedom. The large size of the interface heat exchanger is dictated by the relatively poor water-side heat transfer characteristics. The objective of this program is to develop a single-phase heat transfer approach which can achieve heat fluxes and heat transfer coefficients comparable to those of the evaporation ammonia side. A new heat exchanger concept has been developed to meet these objecties. The main feature of this heat exchanger is that it can achieve very high heat fluxes with a pressure drop one to two orders of magnitude lower than those of previous microchannel or jet impingement high heat flux heat exchangers. This paper describes proof-of-concept experiments performed in air and water and presents analytical model of the heat exchanger.

Valenzuela, Javier A.↗
NASA NTRS2004

Gas-Liquid Flows and Phase Separation

Common issues for space system designers include:Ability to Verify Performance in Normal Gravity prior to Deployment; System Stability; Phase Accumulation & Shedding; Phase Separation; Flow Distribution through Tees & Manifolds Boiling Crisis; Heat Transfer Coefficient; and Pressure Drop.The report concludes:Guidance similar to "A design that operates in a single phase is less complex than a design that has two-phase flow" is not always true considering the amount of effort spent on pressurizing, subcooling and phase separators to ensure single phase operation. While there is still much to learn about two-phase flow in reduced gravity, we have a good start. Focus now needs to be directed more towards system level problems .

McQuillen, John↗
NASA NTRS2006

Miniature Distillation Column for Producing LOX From Air

The figure shows components of a distillation column intended for use as part of a system that produces high-purity liquid oxygen (LOX) from air by distillation. (The column could be easily modified to produce high-purity liquid nitrogen.) Whereas typical industrial distillation columns for producing high-purity liquid oxygen and/or nitrogen are hundreds of feet tall, this distillation column is less than 3 ft (less than about 0.9 m) tall. This column was developed to trickle-charge a LOX-based emergency oxygen system (EOS) for a large commercial aircraft. A description of the industrial production of liquid oxygen and liquid nitrogen by distillation is prerequisite to a meaningful description of the present miniaturized distillation column. Typically, such industrial production takes place in a chemical processing plant in which large quantities of high-pressure air are expanded in a turboexpander to (1) recover a portion of the electrical power required to compress the air and (2) partially liquefy the air. The resulting two-phase flow of air is sent to the middle of a distillation column. The liquid phase is oxygen-rich, and its oxygen purity increases as it flows down the column. The vapor phase is nitrogen-rich and its nitrogen purity increases as it flows up the column. A heater or heat exchanger, commonly denoted a reboiler, is at the bottom of the column. The reboiler is so named because its role is to reboil some of the liquid oxygen collected at the bottom of the column to provide a flow of oxygen-rich vapor. As the oxygen-rich vapor flows up the column, it absorbs the nitrogen in the down-flowing liquid by mass transfer. Once the vapor leaves the lower portion of the column, it interacts with down-flowing nitrogen liquid that has been condensed in a heat exchanger, commonly denoted a condenser, at the top of the column. Liquid oxygen and liquid nitrogen products are obtained by draining some of the purified product at the bottom and top of the column, respectively. Because distillation is a mass-transfer process, the purity of the product(s) can be increased by increasing the effectiveness of the mass-transfer process (increasing the mass-transfer coefficient) and/or by increasing the available surface area for mass transfer through increased column height. The diameter of a distillation column is fixed by pressure-drop and mass-flow requirements. The approach taken in designing the present distillation column to be short yet capable of yielding a product of acceptably high purity was to pay careful attention to design details that affect mass-transfer processes.

Rozzi, Jay C.↗
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 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

CFD Model Development of a Cryogenic Storage Tank Self-Pressurization in Normal Gravity and Validation against SHIIVER Experiment

Two-phase flow and heat transfer simulations with interfacial phase change of the Structural Heat Intercept, Insulation, and Vibration Evaluation Rig (SHIIVER) self-pressurization experiment were conducted using storage tank CFD model in the framework of the ANSYS Fluent CFD code. The simulations were performed for the 70% fill level case with MLI on domes and no vapor cooling. All the phase change calculations in these simulations were generated by in-house Schrage-based evaporation-condensation model. The calculations were performed using both the Volume of Fluid (VOF) and Sharp Interface multiphase 2D axisymmetric models. A number of parametric and sensitivity studies were performed to check the various aspects of the CFD model. These studies helped to understand the effects of varying several parameters on the tank pressure and temperature during self-pressurization. Turbulence modeling; turbulence damping at the interface; and using constant vs. temperature dependent fluid properties were shown to have the most profound influence on predicted tank pressures and temperatures. Current study indicates that including phase change at the interface into the computational model is crucial for accurate prediction of the tank self-pressurization process. The effect of accommodation coefficient was also studied. Tank pressure values predicted by the VOF, and Sharp Interface models are within 4% of the experimental ones.

Computational Fluid Dynamics↗
NASA NTRSDate not supplied

CFD Model Development of a Cryogenic Storage Tank Self-Pressurization in Normal Gravity and Validation against SHIIVER Experiment

Two-Phase flow and heat transfer simulations with interfacial phase change of the Structural Heat Intercept, Insulation, and Vibration Evaluation Rig (SHIIVER) self-pressurization experiment were conducted using storage tank CFD model in the framework of the ANSYS Fluent CFD code. The simulations were performed for the 70% fill level case with MLI on domes and no vapor cooling. All the phase change calculations in these simulations were generated by in-house Schrage-based evaporation-condensation model. The calculations were performed using both the VOF and Sharp Interface multiphase 2Daxisymmetric models. A number of parametric and sensitivity studies were performed to check the various aspects of the CFD model. These studies helped to understand the effects of varying several parameters on the tank pressure and temperature during self-pressurization. These parameters include: the effect of turbulence modeling and turbulence damping at the interface; the effect of constant vs. temperature dependent fluid properties; the effect of accommodation coefficient; the effect of modeling phase change at the interface and the effect of the fill level.

Computational Fluid Dynamics↗