Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Vapor liquid equilibrium”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
NASA NTRS2008

The Stationary Condensation and Radial Outflow of a Liquid Film on a Horizontal Disk

The application of capillary screen liquid acquisition devices to space-based cryogenic propulsion systems is expected to necessitate thermodynamic conditioning in order to stabilize surface tension retention characteristics. The present results have been obtained in the framework of the research of low gravity condensation-flow processes for conditioning cryogenic liquid acquisition devices. The following system is studied: On the top of a subcooled horizontal disk, a liquid film condenses from the ambient saturated vapor. The liquid is forcedly removed at the disk edge, and there is an outward radial flow of the film. Stationary regimes of the flow are uncovered such that (i) the gravity is negligible, being eclipsed by the capillary forces; (ii) the film thickness is everywhere much smaller than the disk radius; and (iii) the slow-flow lubrication approximation is valid. A nonlinear differential equation for the film thickness as a function of the radial coordinate is obtained. The (two-dimensional) fields of velocities, temperature and pressure in the film are explicitly determined by the radial profile of its thickness. The equilibrium is controlled by two parameters: (i) the vapor-disk difference of temperatures and (ii) the liquid exhaust rate. For the flow regimes with a nearly uniform film thickness, the governing equation linearizes, and the film interface is analytically predicted to have a concave-up quartic parabola profile. Thus, perhaps counter-intuitively, the liquid film is thicker at the edge and thinner at the center of the disk.

Bolshinskiy, Leonid↗
DOE OSTI2021

Influence of Non-Equilibrium Conditions on Liquid Hydrogen Storage Tank Behavior [Slides]

Understanding liquid hydrogen tank fluid dynamics is key for modeling liquid hydrogen systems. The tank is the source for nearly all liquid hydrogen systems. Accurate flow modeling out of the tank is needed to predict flows through downstream components. Tank contains liquid and gas that may not be at equilibrium. Questions to be addressed are: Does heat and mass transfer between liquid and vapor affect the flow rate? Is boiling an important consideration? For what conditions is a pressure relief valve (PRV) sufficient to relieve pressure and when is the burst disc needed?

08 HYDROGEN↗
NASA NTRS1973

Cavitation in liquid cryogens. 2: Hydrofoil

Boundary layer principles, along with two-phase concepts, are used to improve existing correlative theory for developed cavity data. Details concerning cavity instrumentation, data analysis, correlative techniques, and experimental and theoretical aspects of a cavitating hydrofoil are given. Both desinent and thermodynamic data, using liquid hydrogen and liquid nitrogen, are reported. The thermodynamic data indicated that stable thermodynamic equilibrium exists throughout the vaporous cryogen cavities. The improved correlative formulas were used to evaluate these data. A new correlating parameter based on consideration of mass limiting two-phase flow flux across the cavity interface, is proposed. This correlating parameter appears attractive for future correlative and predictive applications. Agreement between theory and experiment is discussed, and directions for future analysis are suggested. The front half of the cavities, developed on the hydrofoil, may be considered as parabolically shaped.

Hord, J.↗
NASA NTRS2008

Carbonate in Comets: A Comparison of Comets 1P/Halley, 9P/Temple 1, and 81P/Wild 2

Comets are generally believed to have formed in a cold region, trapping in the cometary ices the original low-temperature condensate grains of our Solar System. These grains would have been preserved in cold-storage, at a temperature below the freezing point of CO2, for the last 4.5+ billion years. Carbonates are common in hydrous meteorites and hydrous interplanetary dust particles (IDPs), where they are believed to have formed by parent-body aqueous processing. Since simple models of cometary evolution involve no aqueous processing, carbonates were generally presumed not to occur in comets. However, Toppani et al. [1] have performed experiments that indicate carbonate can be formed by non-equilibrium condensation in circumstellar environments where water is present as a vapor, not as a liquid. This suggests carbonate might have condensed in cold regions of the Solar Nebula, and might be present in comets.

Flynn, G. J.↗
NASA NTRS1994

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

Previous experiments indicate that the bubble point pressure of spacecraft liquid hydrogen acquisition devices is reduced substantially when the ullage is pressurized with heated hydrogen vapor. The objective is to determine whether the two-phase thermocapillary convection arising from thermodynamic non-equilibrium along the porous surfaces of such devices could lead to this observed degradation in retention performance. We also examine why retention capability appears to be unaffected by pressurization with heated helium or direct heating through the porous structure. Computational assessments based on coupled solution of the flowfield and liquid free surface indicate that for highly wetting fluids in small pores, dynamic pressure and vapor recoil dictate surface morphology and drive meniscus deformation. With superheating, the two terms exert the same influence on curvature and promote mechanical equilibrium, but with subcooling, the pressure distribution produces a suction about the pore center-line that degrades retention. This result points to thermocapillary-induced deformation arising from condensation as the cause for retention loss. It also indicates that increasing the level of non-equilibrium by reducing accommodation coefficient restricts deformation and explains why retention failure does not occur with direct screen heating or helium pressurization.

Schmidt, G. R.↗
DOE OSTI2024

Generalized Brunauer–Emmett–Teller Isotherm for Mixed-Gas Multilayer Adsorption Equilibria

Here, this work presents a rigorous thermodynamic framework to calculate mixed-gas multilayer adsorption equilibria from single-component isotherms and vapor–liquid equilibria. Named the generalized Brunauer–Emmett–Teller (gBET) isotherm, the newly formulated isotherm considers adsorbent surface heterogeneity, competitive adsorption on the monolayer, condensation–evaporation on the subsequent layers, and adsorbed phase nonideality for the monolayer and the subsequent layers. The monolayer adsorbed phase nonideality is tracked by using an area-based adsorption nonrandom two-liquid activity coefficient model. The adsorbed phase composition and corresponding nonideality in the subsequent layers are calculated at the dew point condition of the mixed gas with either an equation of state or an activity coefficient model for the vapor–liquid equilibria. The proposed model is validated with six single and three binary multilayer adsorption equilibrium systems, and the model results are compared against those from the classical BET isotherm for single-component adsorption and ideal adsorbed solution theory for mixed-gas adsorption equilibria.

09 BIOMASS FUELS↗
NASA NTRSDate not supplied

Parametric Analysis of the Charge-Hold-Vent Method for Cryogenic Propellant Tank Chill Down

In the absence of external heat exchangers, the on-orbit transfer of cryogenic propellants requires the receiver tank to first be quenched to a sufficiently low energy state to allow for a continuous no-vent fill to avoid unnecessary venting of liquid. One proposed method for tank chilldown that minimizes the potential for venting liquid is the charge hold vent (CHV) method. CHV follows a cyclic process that gradually removes thermal energy from the receiver tank by injecting liquid with the vent valve closed and allowing the fluid and wall to reach near-thermal equilibrium before venting the superheated vapor. However, the CHV method must be optimized to minimize complexity, mass, and time. This paper presents a modular CHV analytical model used to quantify the number of cycles and propellant mass consumed based on first principles. The model is used to examine the effect of eight parameters: receiver tank material, volume, mass, maximum expected operating pressure, and initial pressure, liquid injection pressure and temperature, and the target temperature. The model is validated against the only two available CHV datasets. Based on results, the tank mass-to-volume ratio is the most important factor in determining the number of CHV cycles and thus degree of difficulty in tank chilldown. The model can easily be used for early-stage design, sizing, and analysis of cryogenic propellant transfer systems.

Tank Chilldown↗
NASA NTRSDate not supplied

Parametric Analysis of the Charge-Hold-Vent Method for Cryogenic Propellant Tank Chill Down

Abstract. In the absence of external heat exchangers, the on-orbit transfer of cryogenic propellants requires the receiver tank to first be quenched to a sufficiently low energy state to allow for a continuous no-vent fill to avoid unnecessary venting of liquid. One proposed method for tank chilldown that minimizes the potential for venting liquid is the charge hold vent (CHV) method. CHV follows a cyclic process that gradually removes thermal energy from the receiver tank by injecting liquid with the vent valve closed and allowing the fluid and wall to reach near-thermal equilibrium before venting the superheated vapor. However, the CHV method must be optimized to minimize complexity, mass, and time. This paper presents a modular CHV analytical model used to quantify the number of cycles and propellant mass consumed based on first principles. The model is used to examine the effect of eight parameters: receiver tank material, volume, mass, maximum expected operating pressure, and initial pressure, liquid injection pressure and temperature, and the target temperature. The model is validated against the only two available CHV datasets. Based on results, the tank mass-to-volume ratio is the most important factor in determining the number of CHV cycles and thus degree of difficulty in tank chilldown. The model can easily be used for early-stage design, sizing, and analysis of cryogenic propellant transfer systems. Keywords: tank chilldown, charge-hold-vent, no-vent fill

Tank Chilldown↗
NASA NTRS1986

A search for the prewetting line

This paper describes efforts to locate the prewetting line in a binary liquid system (isopropanol-perfluoromethylcyclohexane) at the vapor-liquid interface. Tight upper bounds were placed on the temperature separation (0.2 K) between the prewetting line and the line of bulk liquid phase separation. The prewetting line in systems at equilibrium was not detected. Experimental signatures indicative of the prewetting line occurred only in nonequilibrium situations. Several theories predict that the adsorption of one of the components (the fluorocarbon, in this case) at the liquid-vapor interface should increase abruptly, at a temperature sightly above the temperature at which the mixture separates into two liquid phases. A regular solution calculation indicates that this prewetting line should have been easily detectable with the instruments used in this experiment. Significant features of the experiment are: (1) low-gradient thermostatting, (2) in situ stirring, (3) precision ellipsometry from the vapor-liquid interface, (4) high resolution differential index of refraction measurements using a novel cell design, and (5) computer control.

Schmidt, J. W.↗
DOE OSTI2020

RELAP-7 Application and Enhancement for FLEX Strategies and ATF Behavior under Extended Loss of AC Power Conditions

This report summarizes the results of a three-year research project sponsored by the U.S. Department of Energy (DOE) Nuclear Energy University Program (NEUP) to enhance and apply the RELAP-7 code by adding and improving several important components (e.g., a mechanistic Reactor Core Isolation Cooling (RCIC) system model) for thermal hydraulic studies of LWRs under ELAP conditions and evaluating the time available for transition to portable FLEX equipment. The project team included University of Massachusetts–Lowell, The Ohio State University, Texas A&M University, Idaho National Laboratory and Oak Ridge National Laboratory. In the Fukushima accident, it was found that the RCIC system played a crucial role in delaying core meltdown by almost three days in Fukushima Daiichi Unit 2, because of self-regulated operation of the steam driven RCIC turbine-pump injection system. Steam flow in the convergent-divergent nozzles of the RCIC Terry turbine is two-phase non-equilibrium transonic flow with homogenous nucleation condensation. To more accurately predict the dynamic process and behavior of the transonic compressible steam flow, a one-dimensional transient two-phase analytical model is presented. A simplified four-fluid model was employed in the present work with the consideration of four separate fluid fields: vapor, liquid film, entrained droplets and condensed droplets. The mass, momentum and energy interactions between the fluids were considered and modeled. An extended seven-equation non-equilibrium critical flow model was developed to obtain the critical pressure and velocities of each phase at the nozzle throat. To predict the wetness in the divergent section, a mechanistic nucleation condensation model was integrated in the nozzle analysis model, considering the generation and consequent growth of droplets. The governing differential equations on a staggered grid were discretized using the second-order Lax-Wendroff scheme with a flux limiter, and the Semi-Implicit Method for Pressure-Linked Equation (SIMPLE) algorithm was employed to solve the discrete linear system. To demonstrate the predictability and reliability of the physical models and the numerical method proposed in the present work, three representative nozzles were modeled and simulated. The results show good agreement with the available experimental data, even for condensation shock. Then, the 1D nozzle model was employed to obtain nozzle flow tables of the Terry turbine nozzle for different working pressures which can cover the operation pressure range of the RCIC system. A mechanistic RCIC turbine-pump system model was developed and implemented in the system code TRACE to simulate dynamic responses of the RCIC system under Beyond Design Basis Accident (BDBA) conditions. The turbine-pump governing equations are based on the control volume approach of the angular momentum balance. The physics based mechanistic RCIC model was developed using the TRACE control system components (i.e., signal variables, control blocks, and tables), and incorporated into a TRACE boiling water reactor (BWR) model. The TRACE model in this report has a detailed nodalization of the reactor pressure vessel (RPV), and all of the major flow paths and system components, including the safety relief valves (SRVs) and the containment suppression pool and drywell. Based on the nozzle flow tables generated from the 1D nozzle model developed, the turbine drive torque can be calculated from table lookup. Since the detailed specifications of the RCIC pump are unavailable, the homologous curves for a Bingham pump were used in the current pump component. A station black-out (SBO) accident test problem was selected to demonstrate the TRACE RCIC model. The short-term SBO simulations were performed for two cladding materials: Zircaloy and FeCrAl, to demonstrate the effect of the accident tolerant fuel cladding on fuel heat-up under BDBA accident conditions. The wetwell plays a vital safety role in SBO and other BWR accident scenarios in that it can reduce containment pressure and supply additional core make-up water. The suppression pool temperature distribution has a very large impact on both RPV and containment pressure. Thus, another novel contribution of the project comes mainly from an improved, systems-level wetwell model which can capture buoyancy-induced thermal stratification effects due to steam injection and condensation. A two-zone stratified wetwell model has been implemented in RELAP-7 and some results from that model are presented. This wetwell model is capable of simulating thermal stratification due to a low steam mass injection rate. With a low mass flow rate, the model assumes that all the steam condenses within the pipe and the resulting plume can be approximated with a purely buoyant, heat-source driven model. The wetwell model developed with these assumptions is adequate to simulate slow transients such as extended SBO transients.

42 ENGINEERING↗
DOE OSTI2026

Intense cavitation-assisted electric discharge as a promising tool for water treatment

This study investigates interrelations between one-electrode Cavitation-Assisted Electric Discharge (CAED), two-electrode CAED, and recently discovered Intense CAED (I-CAED). The one-electrode CAED is a self-triggered nanosecond discharge with pulse energy in the micro-Joule range, which can be generated even by a DC voltage. I-CAED consists of a non-equilibrium part within a low-pressure cavitating region and a micro-spark traversing a liquid film. We hypothesize that CAED propagates from the high-voltage electrode as an ionization wave through bubbles of saturated vapor. Subsequently, the streamer-like discharges in the bubbles may form a continuous plasma channel. Inside the cavitating region, the plasma is strongly non-equilibrium, providing an ideal environment for generating chemically unstable species such as hydrogen peroxide (H 2 O 2 ). Plasma of I-CAED spark is characterized by high electron density and near-thermal equilibrium, emitting a continuous ultraviolet spectrum. The combination of these different discharge parts makes I-CAED in water a highly effective tool for the Advanced Oxidation Process, particularly in water disinfection. Experimentally demonstrated Electric Energy per Order value for disinfection of E. coli-contaminated water is as low as 0.135 ± 0.035 kWh/m 3 /order. Estimates show that the implementation of “dry electrodes” configuration reduces the erosion rate of the electrode material by at least one order of magnitude. Spectral analysis reveals that the continuum emission generated by I-CAED in proximity to metal electrodes deviates from the spectra of discharges spatially decoupled from the electrodes. We assume that this spectral divergence is attributable to blackbody-like emission originating from metallic nanoparticles form during the electrode's erosion process.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

How does humidity data impact land surface modeling of hydrothermal regimes at a permafrost site in Utqia$\mathrm{\dot{g}}$vik, Alaska?

Humidity is a basic and crucial meteorological indicator commonly measured in several forms, including specific humidity, relative humidity, and absolute humidity. These different forms can be inter-derived based on the saturation vapor pressure (SVP). In past decades, dozens of formulae have been developed to calculate the SVP with respect to, and in equilibrium with, liquid water and solid ice surfaces, but many prior studies use a single function for all temperature ranges, without considering the distinction between over the liquid water and ice surfaces. These different approaches can result in humidity estimates that may impact our understanding of surface-subsurface thermal-hydrological dynamics in cold regions. In this study, we compared the relative humidity (RH) downloaded and calculated from four data sources in Alaska based on five commonly used SVP formulas. These RHs, along with other meteorological indicators, were then used to drive physics-rich land surface models at a permafrost-affected site. Here, we found that higher values of RH (up to 40%) were obtained if the SVP was calculated with the over-ice formulation when air temperatures were below freezing, which could lead to a 30% maximum difference in snow depths. The choice of whether to separately calculate the SVP over an ice surface in winter also produced a significant range (up to 0.2 m) in simulated annual maximum thaw depths. The sensitivity of seasonal thaw depth to the formulation of SVP increases with the rainfall rate and the height of above-ground ponded water, while it diminishes with warmer air temperatures. These results show that RH variations based on the calculation of SVP with or without over-ice calculation meaningfully impact physically-based predictions of snow depth, sublimation, soil temperature, and active layer thickness. Under particular conditions, when severe flooding (inundation) and cool air temperatures are present, care should be taken to evaluate how humidity data is estimated for land surface and earth system modeling.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1981

Defining critical point experiments for a space laboratory

Three representative low-gravity experiments for a fluid near its liquid-vapor critical point are being defined. Two of these experiments require very careful measurements of properties of the fluid in thermodynamic equilibrium, while the third experiment is a series of optical observations of the phenomena which occur as a fluid is changed from one phase to two phases, either by cooling through the critical point, or by adiabatic expansion. There is concern that residual spacecraft motions may complicate the interpretation of the data from the proposed experiments. It is possible that the Spacelab environment will render certain desirable experiments impractical.

Moldover, M. R.↗
NASA NTRS1996

Analysis of Two-Phase Flow in Damper Seals for Cryogenic Turbopumps

Cryogenic damper seals operating close to the liquid-vapor region (near the critical point or slightly su-cooled) are likely to present two-phase flow conditions. Under single phase flow conditions the mechanical energy conveyed to the fluid increases its temperature and causes a phase change when the fluid temperature reaches the saturation value. A bulk-flow analysis for the prediction of the dynamic force response of damper seals operating under two-phase conditions is presented as: all-liquid, liquid-vapor, and all-vapor, i.e. a 'continuous vaporization' model. The two phase region is considered as a homogeneous saturated mixture in thermodynamic equilibrium. Th flow in each region is described by continuity, momentum and energy transport equations. The interdependency of fluid temperatures and pressure in the two-phase region (saturated mixture) does not allow the use of an energy equation in terms of fluid temperature. Instead, the energy transport is expressed in terms of fluid enthalpy. Temperature in the single phase regions, or mixture composition in the two phase region are determined based on the fluid enthalpy. The flow is also regarded as adiabatic since the large axial velocities typical of the seal application determine small levels of heat conduction to the walls as compared to the heat carried by fluid advection. Static and dynamic force characteristics for the seal are obtained from a perturbation analysis of the governing equations. The solution expressed in terms of zeroth and first order fields provide the static (leakage, torque, velocity, pressure, temperature, and mixture composition fields) and dynamic (rotordynamic force coefficients) seal parameters. Theoretical predictions show good agreement with experimental leakage pressure profiles, available from a Nitrogen at cryogenic temperatures. Force coefficient predictions for two phase flow conditions show significant fluid compressibility effects, particularly for mixtures with low mass content of vapor. Under these conditions, an increase on direct stiffness and reduction of whirl frequency ratio are shown to occur. Prediction of such important effects will motivate experimental studies as well as a more judicious selection of the operating conditions for seals used in cryogenic turbomachinery.

Arauz, Grigory L.↗
NASA NTRS2021

Experimental and Computational Investigation of Flow Boiling in Microgravity

This study explores use of Computational Fluid Dynamics (CFD) to predict near-saturated flow boiling of FC-72 in microgravity. The computational method employs transient analysis to predict detailed interfacial behavior and heat transfer characteristics along a rectangular channel heated along two opposite walls. Predicted results are validated against experimental temperature measurements and high-speed video images captured during a series of parabolic aircraft maneuvers for three sets of operating conditions which include variations of both mass velocity and wall heat flux. The computational method is based on the multi-phase volume of fluid (VOF) model, which is combined with appropriate phase change and turbulence models, and accounts for both shear-lift force on bubbles and conjugate heat transfer along the heating walls. A key advantage of the CFD method is ability to capture details that are very difficult to measure experimentally, including detailed spatial variations of bubble shape, void fraction, mixture fluid temperature, liquid velocity, and vapor velocity, results for which are presented for each of the three test cases. Different flow regimes predicted along the heated length exhibit a number of dominant mechanisms including bubble nucleation, bubble growth, coalescence, vapor blankets, interfacial waviness, and residual liquid sub-layer, all of which agree well with experiment. Vapor velocity is shown to increase appreciably along the heated length because of increased void fraction, while liquid velocity experiences large fluctuations. Non-equilibrium effects are accentuated with increasing mass velocity, contributing minor deviations of fluid temperature from simulations compared to those predicted by the analytical method. Predicted wall temperature is fairly uniform in the middle of the heated length but increases in the entrance region, due to sensible heat transfer in the subcooled liquid, and decreases toward the exit, mostly because of flow acceleration resulting from increased void fraction.

two-phase cooling; flow boiling; CFD; microgravity↗
NASA NTRS1990

An analysis of the vapor flow and the heat conduction through the liquid-wick and pipe wall in a heat pipe with single or multiple heat sources

A numerical analysis is presented for the overall performance of heat pipes with single or multiple heat sources. The analysis includes the heat conduction in the wall and liquid-wick regions as well as the compressibility effect of the vapor inside the heat pipe. The two-dimensional elliptic governing equations in conjunction with the thermodynamic equilibrium relation and appropriate boundary conditions are solved numerically. The solutions are in agreement with existing experimental data for the vapor and wall temperatures at both low and high operating temperatures.

Chen, Ming-Ming↗
NASA NTRS1999

Determination of the Accomodation Coefficient Using Vapor/Gas Bubble Dynamics in an Acoustic Field

Non-equilibrium liquid/vapor phase transformations can occur in superheated or subcooled liquids in fast processes such as in evaporation in a vacuum, in processing of molten metals, and in vapor explosions. The rate at which such a phase transformation occurs, Xi, can be described by the Hertz-Knudsen-Langmuir formula. More than one century of the history of the accommodation coefficient measurements shows many problems with its determination. This coefficient depends on the temperature, is sensitive to the conditions at the interface, and is influenced by small amounts of impurities. Even recent measurements of the accommodation coefficient for water (Hagen et al, 1989) showed a huge variation in Beta from 1 for 1 micron droplets to 0.006 for 15 micron droplets. Moreover, existing measurement techniques for the accommodation coefficient are complex and expensive. Thus development of a relatively inexpensive and reliable technique for measurement of the accommodation coefficient for a wide range of substances and temperatures is of great practical importance.

Gumerov, Nail A.↗
NASA NTRS1992

A comparison between protein crystals grown with vapor diffusion methods in microgravity and protein crystals using a gel liquid-liquid diffusion ground-based method

Crystals of human serum albumin have been successfully grown in a variety of gels using crystallization conditions otherwise equivalent to those utilized in the popular hanging-drop vapor-equilibrium method. Preliminary comparisons of gel grown crystals with crystals grown by the vapor diffusion method via both ground-based and microgravity methods indicate that crystals superior in size and quality may be grown by limiting solutal convection. Preliminary X-ray diffraction statistics are presented.

Miller, Teresa Y.↗