Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “interfacial evaporation”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Interfacial thermodynamics of cryogenic fluids: The effect on non-condensable gas on fluid storage

Propellant tanks that contain cryogenic fluids (CFs) in low gravity conditions are usually pressurized with non-condensable (NC) gases for fast extraction. Unfortunately, the presence of NC gases causes CFs to exhibit higher boil-offs compared to pure CF systems. For optimal utilization of the cryogenic fuels, these higher boil-off rates must be minimized. Our goal is to quantify the effects NC gases have on the evaporation and condensation dynamics of CFs to develop strategies that mitigate the higher boil-offs. The hypothesis is that the NC gas accumulates around the CF’s liquid vapor interface in the Knudsen layer, creating a kinetic barrier for both evaporation and condensation. Testing this hypothesis with experimental methods is difficult due to the transient nature of the Knudsen layer, which is only a few nanometers thick. Furthermore, experimental capabilities are limited in low gravity conditions and extremely expensive. Accordingly, we approach this problem from a theoretical perspective. In this study, we employ molecular dynamics (MD) simulations to probe the interfacial mechanisms that affect CF’s evaporation and condensation at varying concentrations of NC gases. Specifically, we use nitrogen (N2) and oxygen (O2) as our CFs and neon (Ne) as our NC gas. Using MD simulations, we show that Ne accumulates at N2 liquid vapor interface across a wide range of Ne concentrations, thereby impeding mass transport of N2. Our simulations allow for direct computation of the molar flux as well as the mass accommodation coefficient (MAC) which can then be used as input parameters to continuum fluid dynamics (CFD) models for optimal storage tank design.

Michael Robert DeLyser↗
NASA NTRS1973

Laminar film boiling on inclined isothermal flat plates.

Laminar film boiling from an inclined flat plate has been investigated analytically. Using the singular perturbation scheme, the complete set of Navier-Stokes equations is solved. The zeroth-order perturbation coinciding with the boundary-layer equations for vertical flat plates governs the problem. The higher-order perturbations become important near the leading edge and for large values of the inclination angle from the vertical. The assumption of zero interfacial velocity shows that, except for fluids having large (rho x mu) ratios, the results can be predicted using the vertical flat plate results by defining a modified Grashof parameter containing a cos phi term. When the interfacial shear is considered, the solutions indicate that for fluids having large (rho x mu) ratios, the heat transfer rates will be larger (approximately 15% maximum) than those predicted by the simplified model using zero interfacial velocity. In general, the inclination decreases the rate of heat transfer as well as the rate of evaporation.

Nagendra, H. R.↗
NASA NTRS1982

Non-isothermal effects on SO2 absorption by water droplets. I - Model development. II - Results and discussion

An analytic model of SO2 absorption in a falling water droplet is developed and a simulation of SO2 washout is performed. Nonisothermic effects on drop growth, droplet physical parameters, reaction rates, and multicomponent diffusion are treated in the model. The gas-liquid interface is assumed to be at equilibrium, and interfacial resistance is negligible. Raindrops are simulated as falling from a 2 km height through an atmospheric region containing SO2. The droplets decrease in size from evaporation and cooling, and their slightly basic pH aids SO2 absorption. The simulation indicates higher SO2 absorption at higher altitudes, and desorption may occur at ground level. Isothermal effects are concluded to be significant, and quantification of effects will depend on further modelling.

Reda, M.↗
NASA NTRS2002

New Method Developed to Measure Contact Angles of a Sessile Drop

The spreading of an evaporating liquid on a solid surface occurs in many practical processes and is of importance in a number of practical situations such as painting, textile dyeing, coating, gluing, and thermal engineering. Typical processes involving heat transfer where the contact angle plays an important role are film cooling, boiling, and the heat transfer through heat pipes. The biological phenomenon of cell spreading also is analogous to a drop spreading (ref. 1). In the study of spreading, the dynamic contact angle describes the interfacial properties on solid substrates and, therefore, has been studied by physicists and fluid mechanics investigators. The dynamic contact angle of a spreading nonvolatile liquid drop provides a simple tool in the study of the free-boundary problem, but the study of the spreading of a volatile liquid drop is of more practical interest because the evaporation of common liquids is inevitable in practical processes. The most common method to measure the contact angle, the contact radius, and the height of a sessile drop on a solid surface is to view the drop from its edge through an optical microscope. However, this method gives only local information in the view direction. Zhang and Yang (ref. 2) developed a laser shadowgraphy method to investigate the evaporation of sessile drop on a glass plate. As described here, Zhang and Chao (refs. 3 and 4) improved the method and suggested a new optical arrangement to measure the dynamic contact angle and the instant evaporation rate of a sessile drop with much higher accuracy (less than 1 percent). With this method, any fluid motion in the evaporating drop can be visualized through shadowgraphy without using a tracer, which often affects the field under investigation.

Chao, David F.↗
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 NTRS1995

Thermocapillary flow with evaporation and condensation at low gravity. Part 2: Deformable surface

The free surface behavior of a volatile wetting liquid at low gravity is studied using scaling and numerical techniques. An open cavity model, which was applied in part 1 to investigate fluid flow and heat transfer in non-deforming pores, is used to evaluate the influence of convection on surface morphology with length scales and subcooling/superheating limits of 1 less than or equal to D less than or equal to 10(exp 2) microns and approximately 1 K, respectively. Results show that the menisci shapes of highly wetting fluids are sensitive to thermocapillary flow and to a lesser extent the recoil force associated with evaporation and condensation. With subcooling, thermocapillarity produces a suction about the pore centerline that promotes loss of mechanical equilibrium, while condensation exerts an opposing force that under some conditions offsets this destabilizing influence. With superheating, thermocapillarity and evaporation act in the same direction and mutually foster surface stability. All of these trends are magnified by high capillary and Biot numbers, and the stronger circulation intensities associated with small contact angles. These phenomena strongly depend on the thermal and interfacial equilibrium between the liquid and vapor, and have important ramifications for systems designed to maintain a pressure differential across a porous surface.

Schmidt, G. R.↗
NASA NTRS2015

Effect of Interfacial Turbulence and Accommodation Coefficient on CFD Predictions of Pressurization and Pressure Control in Cryogenic Storage Tank

Laminar models agree closely with the pressure evolution and vapor phase temperature stratification but under-predict liquid temperatures. Turbulent SST k-w and k-e models under-predict the pressurization rate and extent of stratification in the vapor but represent liquid temperature distributions fairly well. These conclusions seem to equally apply to large cryogenic tank simulations as well as small scale simulant fluid pressurization cases. Appropriate turbulent models that represent both interfacial and bulk vapor phase turbulence with greater fidelity are needed. Application of LES models to the tank pressurization problem can serve as a starting point.

Condensation Evaporation↗
NASA NTRS1972

Review of critical flow rate, propagation of pressure pulse, and sonic velocity in two-phase media

For single-phase media, the critical discharge velocity, the sonic velocity, and the pressure pulse propagation velocity can be expressed in the same form by assuming isentropic, equilibria processes. In two-phase mixtures, the same concept is not valid due to the existence of interfacial transports of momentum, heat, and mass. Thus, the three velocities should be treated differently and separately for each particular condition, taking into account the various transport processes involved under that condition. Various attempts are reviewed to predict the critical discharge rate or the propagation velocities by considering slip ratio (momentum change), evaporation (mass and heat transport), flow pattern, etc. Experimental data were compared with predictions based on various theorems. The importance is stressed of the time required to achieve equilibrium as compared with the time available during the process, for example, of passing a pressure pulse.

Hsu, Y.↗
NASA NTRS1994

Effects of thermocapillarity on an evaporating extended meniscus in microgravity

An analytical investigation of the effects of thermocapillarity on the flow field within and heat transfer from the extended meniscus region of a heated meniscus which is re-supplied by capillarity is presented. Microgravity conditions are considered. The analysis shows that even for extremely small temperature differences between the wall and the vapor (less than 1 mK) thermocapillary stresses at the liquid-vapor interface due to a non-uniform interfacial temperature drastically alters the flow field. At the same time, these stresses were shown to have only a slight effect on the heat transfer from the extended meniscus but increasing with an increasing temperature difference. Additionally, thermocapillary effects were shown to be sensitive to pore size. A criterion was established from a scaling analysis identifying the conditions necessary for thermocapillarity to affect the operation of capillary-pumped heat transport devices in microgravity. A critical Marangoni number and corresponding critical temperature difference between wall and vapor were identified.

Hallinan, Kevin P.↗
NASA NTRS1978

A constant heat flux model of the evaporating interline region

An approximate procedure is developed for determining the heat sink capability of an interline region (junction of vapor, evaporating thin film, and nonevaporating thin film), with consideration of the effect of resistance of the solid to heat transfer. The formulation is based on a simple one-dimensional constant heat flux model. The solution is directed toward a limiting case in which the thermal conductance of the solid is controlling. Since the value of the ideal liquid-vapor interfacial heat transfer coefficient is very large, conductance of the solid is found to dominate the overall heat transfer coefficient in the interline region after a small increase in film thickness.

Wayner, P. C., Jr.↗
NASA NTRS1990

Free convection in the Matian atmosphere

The 'free convective' regime for the Martian atmospheric boundary layer (ABL) was investigated. This state occurs when the mean windspeed at the top of the ABL drops below some critical value U(sub c) and positive buoyant forces are present. Such forces can arise either from vertical temperature or water vapor gradients across the atmospheric surface layer. During free convection, buoyant forces drive narrow plumes that ascend to the inversion height with a return circulation consisting of broad slower-moving downdraughts. Horizontal pressure, temperature, windspeed, and water vapor fluctuations resulting form this circulation pattern can be quite large adjacent to the ground (within the surface layer). The local turbulent fluctuations cause non-zero mean surface stresses, sensible heat fluxes, and latent heat fluxes, even when the mean regional windspeed is zero. Although motions above the surface layer are insensitive to the nature of the surface, the sensible and latent heat fluxes are primarily controlled by processes within the interfacial sublayer immediately adjacent to the ground during free convection. Thus the distinction between aerodynamically smooth and rough airflow within the interfacial sublayer is more important than for the more typical situation where the mean regional windspeed is greater than U(sub c). Buoyant forces associated with water vapor gradients are particularly large on Mars at low pressures and high temperatures when the surface relative humidity is 100 percent, enhancing the likelihood of free convection under these conditions. On this basis, Ingersol postulated the evaporative heat losses from an icy surface on Mars at 237 K and current pressures would exceed the available net radiative flux at the surface, thus prohibiting ice from melting at low atmospheric pressures. Schumann has developed equations describing the horizontal fluctuations and mean vertical gradients occurring during free convection. Schumann's model was generalized to include convection driven by water vapor gradients and to include the effects of circulation above both aerodynamically smooth and rough surfaces.

Clow, G. D.↗
NASA NTRS1995

Influence of Two-Phase Thermocapillary Flow on Liquid Retention in Microscopic Pores

An important feature of screened propellant acquisition devices is the retention capability or maximum maintainable pressure difference across the porous barrier separating the liquid and gas. Previous experiments with liquid hydrogen showed a marked reduction in retention when the tank containing the device was pressurized with hydrogen vapor. These tests, however did not indicate any appreciable degradation in retention with helium pressurization or direct heating through the screen. The objective of this article is to determine if the thermocapillary convection arising from phase change in the microscopic pores of such screens could cause these disparities in performance. A numerical model of flow in a single pore suggests that the thermocapillary-induced gradient in liquid pressure along the surface can strongly affect surface morphology. In an evaporative environment, this gradient exerts a stabilizing influence on surface curvature, and preserves the momentum balance between the liquid and gas. With condensation, it causes a force imbalance and a destabilizing suction in the middle of the pore that reduces retention. Results also indicate that introducing an inert gas, such as helium, suppresses this retention loss mechanism by lowering thermocapillary circulation and its associated interfacial pressure gradient.

Schmidt, George R.↗
NASA NTRSDate not supplied

Mastering the Cryogenic Frontier: Predicting Sloshing in Cryogenic Propellant Tanks

In recent years, NASA has undertaken the development and validation of computation fluid dynamic (CFD) models for predicting cryogenic two-phase sloshing within propellant tanks. These models have been rigorously anchored against cryogenic ground tests that induce sloshing, as well as limited microgravity tests. The primary focus of these models is to accurately predict the two-phase condensation and, at times, evaporation within the cryogenic propellant. This approach is important for estimating the required autogenous and/or helium pressurant for completing propulsion system burns or ensuring adequate pressure for in-space propellant transfers. Unlike in hypergolic propulsion systems, where the propellant is stored at room temperature, cryogenic propellant is highly volatile. Therefore, ensuring the accurate estimation of pressurant necessitates the precise modeling of two-phase interfacial heat and mass transfer processes. This paper aims to provide an overview of NASA’s continuing efforts in developing and validating these CFD models, with particular emphasis on benchmark cases with sloshing that enhances phase change, tracing their historical development, and demonstrating their predictive accuracy when compared to test results.

low-gravity fluid modeling↗
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↗
NASA NTRS2022

Flow Visualization, Heat Transfer, and Critical Heat Flux of Flow Boiling in Earth Gravity with Saturated Liquid-Vapor Mixture Inlet Conditions – In Preparation for Experiments Onboard the International Space Station

This study investigates flow boiling of n-Perfluorohexane with saturated two-phase mixture inlet in a rectangular channel of dimensions 114.6-mm heated length, 2.5-mm width, and 5 mm-height. The experiments were performed as part of the Mission Sequence Testing of the Flow Boiling and Condensation Experiment’s (FBCE) Flow Boiling Module (FBM) in the vertical upflow configuration in Earth gravity using the same experimental system that was launched to the International Space Station (ISS) in August 2021. The operating parameters varied are heating configuration (single- and double-sided), mass velocity (380 – 2400 kg/m2s), inlet quality (0.011 – 0.519), and inlet pressure (120 – 179 kPa). High-speed video photographs are presented to explain the two-phase flow patterns within the channel’s heated length. Flow patterns are constituted by low-density and high-density fronts moving along the channel, with the high-density fronts gradually reducing in length due to evaporation. Heat transfer results in terms of flow boiling curves, streamwise wall temperature profiles, streamwise heat transfer coefficient profiles, and average heat transfer coefficients are presented and trends discussed. CHF data from the present experiments are combined with prior databases to compile a consolidated FBCE-CHF database for saturated inlet to expand the ranges of operating conditions and include other flow orientations in Earth gravity. Experimental CHF trends are also discussed. The interfacial lift-off model shows a good CHF predictive accuracy evidenced by a mean absolute error of 11.97% for this consolidated database after constraining it to mass velocities greater than or equal to 500 kg/m2s. Finally, this study confirmed reliability of the upcoming ISS experiments for saturated inlet conditions and the collected Earth-gravity data will be compared to ISS microgravity data.

Flow boiling↗
NASA NTRS2017

Pressurization of a Flightweight, Liquid Hydrogen Tank: Evaporation & Condensation at a Liquid/Vapor Interface

This paper presents an analysis and simulation of evaporation and condensation at a motionless liquid/vapor interface. A 1-D model equation, emphasizing heat and mass transfer at the interface, is solved in two ways, and incorporated into a subgrid interface model within a CFD simulation. Simulation predictions are compared with experimental data from the CPST Engineering Design Unit tank, a cryogenic fluid management test tank in 1-g. The numerical challenge here is the physics of the liquid/vapor interface; pressurizing the ullage heats it by several degrees, and sets up an interfacial temperature gradient that transfers heat to the liquid phase-the rate limiting step of condensation is heat conducted through the liquid and vapor. This physics occurs in thin thermal layers O(1 mm) on either side of the interface which is resolved by the subgrid interface model. An accommodation coefficient of 1.0 is used in the simulations which is consistent with theory and measurements. This model is predictive of evaporation/condensation rates, that is, there is no parameter tuning.

Stewart, Mark E. M.↗
NASA NTRS2003

Combustion Of Moving Droplets And Of Droplets Suspended Within A Convective Environment: Transient Numerical Results

The problem considered is that of a single-component liquid fuel (n-heptane) droplet undergoing evaporation and combustion in a hot, convective, low pressure, zero-gravity environment of infinite expanse. For a moving droplet, the relative velocity (U(sub infinity)) between the droplet and freestream is subject to change due to the influence of the drag force on the droplet. For a suspended droplet, the relative velocity is kept constant. The governing equations for the gas-phase and the liquid-phase consist of the unsteady, axisymmetric equations of mass, momentum, species (gas-phase only) and energy conservation. Interfacial conservation equations are employed to couple the two phases. Variable properties are used in the gas- and liquid-phase. Multicomponent diffusion in the gas-phase is accounted for by solving the Stefan-Maxwell equations for the species diffusion velocities. A one-step overall reaction is used to model the combustion. The governing equations are discretized using the finite volume and SIMPLEC methods. A colocated grid is adopted. Hyperbolic tangent stretching functions are used to concentrate grid points near the fore and aft lines of symmetry and at the droplet surface in both the gas- and liquid-phase. The discretization equations are solved using the ADI method with the TDMA used on each line of the two alternating directions. Iterations are performed within each time-step until convergence is achieved. The grid spacing, size of the computational domain and time-step were tested to ensure that all solutions are independent of these parameters. A detailed discussion of the numerical model is given.

Gogos, George↗