Search NASA⌕ Search

SEARCH · Search NASA

Results for “LIQUID-VAPOR INTERFACE”

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.

191 records · Page 11

Vaporization and combustion of fuel droplets at supercritical conditions

Vaporization and combustion liquid-fuel droplets in both sub- and super-critical environments have been examined. The formulation is based on the complete conservation equations for both gas and liquid phases, and accommodates finite-rate chemical kinetics and a full treatment of liquid-vapor phase equilibrium at the droplet surface. The governing equations and the associated interface boundary conditions are solved numerically using a fully coupled, implicit scheme with the dual time-stepping integration technique. The model is capable of treating the entire droplet history, including the transition from the subcritical to the supercritical state. As a specific example, the combustion of n-pentane fuel droplets in air is studied for pressures of 5-140 atm. In addition, the dynamic responses of droplet vaporization and combustion to ambient-pressure oscillations are investigated. Results indicate that the droplet gasification and burning mechanisms depend greatly on the ambient pressure. In particular, a rapid enlargement of the vaporization and combustion responses occurs when the droplet surface reaches its critical point, mainly due to the strong variations of latent heat of vaporization and thermophysical properties at the critical state.

Yang, Vigor↗

Capillary condenser/evaporator

A heat transfer device is disclosed for transferring heat to or from a fluid that is undergoing a phase change. The heat transfer device includes a liquid-vapor manifold in fluid communication with a capillary structure thermally connected to a heat transfer interface, all of which are disposed in a housing to contain the vapor. The liquid-vapor manifold transports liquid in a first direction and conducts vapor in a second, opposite direction. The manifold provides a distributed supply of fluid (vapor or liquid) over the surface of the capillary structure. In one embodiment, the manifold has a fractal structure including one or more layers, each layer having one or more conduits for transporting liquid and one or more openings for conducting vapor. Adjacent layers have an increasing number of openings with decreasing area, and an increasing number of conduits with decreasing cross-sectional area, moving in a direction toward the capillary structure.

Valenzuela, Javier A.↗

Progress on Acoustic Measurements of the Bulk Viscosity of Near-Critical Xenon (BVX)

We plan to determine the bulk viscosity of xenon 10 times closer [in reduced temperature tau = (T-Tc)/Tc] to its liquid-vapor critical point than ever before. (Tc is the critical temperature.) To do so, we must measure the dispersion and attenuation of sound at frequencies 1/100 of those used previously. In general, sound attenuation has contributions from the bulk viscosity acting throughout the volume of the xenon as well as contributions from the thermal conductivity and the shear viscosity acting within thin thermoacoustic boundary layers at the interface between the xenon and the solid walls of the resonator. Thus, we can determine the bulk viscosity only when the boundary layer attenuation is small and well understood. We present a comparison of calculations and measurements of sound attenuation in the acoustic boundary layer of xenon near its liquid-vapor critical point.

Gillis, Keith A.↗

Modeling of Bulk Evaporation and Condensation

This report describes the modeling and mathematical formulation of the bulk evaporation and condensation involved in liquid-vapor phase change processes. An internal energy formulation, for these phase change processes that occur under the constraint of constant volume, was studied. Compared to the enthalpy formulation, the internal energy formulation has a more concise and compact form. The velocity and time scales of the interface movement were obtained through scaling analysis and verified by performing detailed numerical experiments. The convection effect induced by the density change was analyzed and found to be negligible compared to the conduction effect. Two iterative methods for updating the value of the vapor phase fraction, the energy based (E-based) and temperature based (T-based) methods, were investigated. Numerical experiments revealed that for the evaporation and condensation problems the E-based method is superior to the T-based method in terms of computational efficiency. The internal energy formulation and the E-based method were used to compute the bulk evaporation and condensation processes under different conditions. The evolution of the phase change processes was investigated. This work provided a basis for the modeling of thermal performance of multi-phase nuclear fuel elements under variable gravity conditions, in which the buoyancy convection due to gravity effects and internal heating are involved.

Anghaie, S.↗

Experimental Investigation of Pool Boiling Heat Transfer Enhancement in Microgravity in the Presence of Electric Fields

Boiling is an effective mode of heat transfer since high heat flux levels are possible driven by relatively small temperature differences. The high heat transfer coefficients associated with boiling have made the use of these processes increasingly attractive to aerospace engineering. Applications of this type include compact evaporators in the thermal control of aircraft avionics and spacecraft environments, heat pipes, and use of boiling to cool electronic equipment. In spite of its efficiency, cooling based on liquid-vapor phase change processes has not yet found wide application in aerospace engineering due to specific problems associated with the low gravity environment. After a heated surface has reached the superheat required for the initiation of nucleate boiling, bubbles will start forming at nucleation sites along the solid interface by evaporation of the liquid. Bubbles in contact with the wall will continue growing by this mechanism until they detach. In terrestrial conditions, bubble detachment is determined by the competition between body forces (e.g. buoyancy) and surface tension forces that act to anchor the bubble along the three phase contact line. For a given body force potential and a balance of tensions along the three phase contact line, bubbles must reach a critical size before the body force can cause them to detach from the wall. In a low gravity environment the critical bubble size for detachment is much larger than under terrestrial conditions, since buoyancy is a less effective means of bubble removal. Active techniques of heat transfer enhancement in single phase and phase change processes by utilizing electric fields have been the subject of intensive research during recent years. The field of electrohydrodynamics (EHD) deals with the interactions between electric fields, flow fields and temperature fields. Previous studies indicate that in terrestrial applications nucleate boiling heat transfer can be increased by a factor of 50 as compared to values obtained for the same system without electric fields. Imposing an external electric field holds the promise to improve pool boiling heat transfer in low gravity, since a phase separation force other than gravity is introduced. The goal of our research is to experimentally investigate the potential of EHD and the mechanisms responsible for EHD heat transfer enhancement in boiling in low gravity conditions.

Herman, Cila↗

Numerical Modeling of Liquid-Vapor Phase Change

We implemented a two- and three-dimensional finite difference/front tracking technique to solve liquid-vapor phase change problems. The mathematical and the numerical features of the method were explained in great detail in our previous reports, Briefly, we used a single formula representation which incorporated jump conditions into the governing equations. The interfacial terms were distributed as singular terms using delta functions so that the governing equations would be the same as conventional conservation equations away from the interface and in the vicinity of the interface they would provide correct jump conditions. We used a fixed staggered grid to discretize these equations and an unstructured grid to explicitly track the front. While in two dimensions the front was simply a connection of small line segments, in three dimensions it was represented by a connection of small triangular elements. The equations were written in conservative forms and during the course of computations we used regriding to control the size of the elements of the unstructured grid. Moreover, we implemented a coalescence in two dimensions which allowed the merging of different fronts or two segments of the same front when they were sufficiently close. We used our code to study thermocapillary migration of bubbles, burst of bubbles at a free surface, buoyancy-driven interactions of bubbles, evaporation of drops, rapid evaporation of an interface, planar solidification of an undercooled melt, dendritic solidification, and a host of other problems cited in the reference.

Esmaeeli, Asghar↗

Platinum liquid-vapor phase boundary mapped by fluid flyer experiments

We report a direct measurement of the temperature and density of a metal along its liquid-vapor coexistence (L-V) curve. By shocking platinum to a high-pressure liquid, we imparted sufficient heat for subsequent isentropic release to place it in a state on the boundary between the liquid and vapor phases. Released material in the liquid phase acted as a high velocity flyer pinned to the L-V curve. We measured velocity and radiant emission of the flyer as well as the interface motion and transiting shock states induced in a downstream window material by its impact. We used these measurements to calculate temperature and density of the L-V curve state which we compare to density functional theory predictions.

Equations of state↗

Sensitivity of magnetohydrodynamic simulations of Joule-heated conductors to the vaporization curve in equations of state

Magnetohydrodynamic (MHD) simulations of electrically exploded aluminum and copper rods demonstrate a technique to validate equations of state (EOS) for rapidly Joule-heated conductors. The balance of internal and magnetic forces at the conductor-insulator interface drives the metal there along the vaporization phase boundary. Variations between critical points and vaporization curves in existing models predict differing densities and temperatures in MHD simulations for these models. Here, the inclusion of Maxwell constructs in the liquid-vapor biphase region of the EOS caused the rod surface to vaporize earlier in time than unmodified tables with van der Waals loops. Velocimetry of recent experiments is used to validate the location of the vaporization curve in existing EOS models and differentiate between the vapor dome treatments. Dielectric coatings applied to the metal surface restricted the conductor’s expansion and diverted the metal into the warm dense matter regime.

36 MATERIALS SCIENCE↗

Modeling of zero gravity venting: Studies of two-phase heat transfer under reduced gravity

The objective is to predict the pressure response of a saturated liquid-vapor system when undergoing a venting or depressurization process in zero gravity at low vent rates. An experimental investigation of the venting of cylindrical containers partially filled with initially saturated liquids was previously conducted under zero-gravity conditions and compared with an analytical model which incorporated the effect of interfacial mass transfer on the ullage pressure response during venting. A new model is presented to improve the estimation of the interfacial mass transfer. Duhammel's superposition integral is incorporated to approximate the transient temperature response of the interface, treating the liquid as a semi-infinite solid with conduction heat transfer. Account is also taken of the condensation taking place within the bulk of a saturated vapor as isentropic expansion takes place. Computational results are presented for the venting of R-11 from a given vessel and initial state for five different venting rates over a period of three seconds, and compared to prior NASA experiments. An improvement in the prediction of the final pressure takes place, but is still considerably below the measurements.

Merte, H., Jr.↗

Finite-Size Effects on the Behavior of the Susceptibility in van der Waals Films Bounded by Strongly Absorbing Substrates

We study critical point finite-size effects in the case of the susceptibility of a film in which interactions are characterized by a van der Waals-type power law tail. The geometry is appropriate to a slab-like system with two bounding surfaces. Boundary conditions are consistent with surfaces that both prefer the same phase in the low temperature, or broken symmetry, state. We take into account both interactions within the system and interactions between the constituents of the system and the material surrounding it. Specific predictions are made with respect to the behavior of 3He and 4He films in the vicinity of their respective liquid-vapor critical points.

phase-transition↗

Use of Generalized Fluid System Simulation Program (GFSSP) for Teaching and Performing Senior Design Projects at the Educational Institutions

This paper describes the experience of the authors in using the Generalized Fluid System Simulation Program (GFSSP) in teaching Design of Thermal Systems class at University of Alabama in Huntsville. GFSSP is a finite volume based thermo-fluid system network analysis code, developed at NASA/Marshall Space Flight Center, and is extensively used in NASA, Department of Defense, and aerospace industries for propulsion system design, analysis, and performance evaluation. The educational version of GFSSP is freely available to all US higher education institutions. The main purpose of the paper is to illustrate the utilization of this user-friendly code for the thermal systems design and fluid engineering courses and to encourage the instructors to utilize the code for the class assignments as well as senior design projects. The need for a generalized computer program for thermofluid analysis in a flow network has been felt for a long time in aerospace industries. Designers of thermofluid systems often need to know pressures, temperatures, flow rates, concentrations, and heat transfer rates at different parts of a flow circuit for steady state or transient conditions. Such applications occur in propulsion systems for tank pressurization, internal flow analysis of rocket engine turbopumps, chilldown of cryogenic tanks and transfer lines, and many other applications of gas-liquid systems involving fluid transients and conjugate heat and mass transfer. Computer resource requirements to perform time-dependent, three-dimensional Navier-Stokes computational fluid dynamic (CFD) analysis of such systems are prohibitive and therefore are not practical. Available commercial codes are generally suitable for steady state, single-phase incompressible flow. Because of the proprietary nature of such codes, it is not possible to extend their capability to satisfy the above-mentioned needs. Therefore, the Generalized Fluid System Simulation Program (GFSSP1) has been developed at NASA Marshall Space Flight Center (MSFC) as a general fluid flow system solver capable of handling phase changes, compressibility, mixture thermodynamics and transient operations. It also includes the capability to model external body forces such as gravity and centrifugal effects in a complex flow network. The objectives of GFSSP development are: a) to develop a robust and efficient numerical algorithm to solve a system of equations describing a flow network containing phase changes, mixing, and rotation; and b) to implement the algorithm in a structured, easy-to-use computer program. The analysis of thermofluid dynamics in a complex network requires resolution of the system into fluid nodes and branches, and solid nodes and conductors as shown in Figure 1. Figure 1 shows a schematic and GFSSP flow circuit of a counter-flow heat exchanger. Hot nitrogen gas is flowing through a pipe, colder nitrogen is flowing counter to the hot stream in the annulus pipe and heat transfer occurs through metal tubes. The problem considered is to calculate flowrates and temperature distributions in both streams. GFSSP has a unique data structure, as shown in Figure 2, that allows constructing all possible arrangements of a flow network with no limit on the number of elements. The elements of a flow network are boundary nodes where pressure and temperature are specified, internal nodes where pressure and temperature are calculated, and branches where flowrates are calculated. For conjugate heat transfer problems, there are three additional elements: solid node, ambient node, and conductor. The solid and fluid nodes are connected with solid-fluid conductors. GFSSP solves the conservation equations of mass and energy, and equation of state in internal nodes to calculate pressure, temperature and resident mass. The momentum conservation equation is solved in branches to calculate flowrate. It also solves for energy conservation equations to calculate temperatures of solid nodes. The equations are coupled and nonlinear; therefore, they are solved by an iterative numerical scheme. GFSSP employs a unique numerical scheme known as simultaneous adjustment with successive substitution (SASS), which is a combination of Newton-Raphson and successive substitution methods. The mass and momentum conservation equations and the equation of state are solved by the Newton-Raphson method while the conservation of energy and species are solved by the successive substitution method. GFSSP is linked with two thermodynamic property programs, GASP2 and WASP3 and GASPAK4, that provide thermodynamic and thermophysical properties of selected fluids. Both programs cover a range of pressure and temperature that allows fluid properties to be evaluated for liquid, liquid-vapor (saturation), and vapor region. GASP and WASP provide properties of 12 fluids. GASPAK includes a library of 36 fluids. GFSSP has three major parts. The first part is the graphical user interface (GUI), visual thermofluid analyzer of systems and components (VTASC). VTASC allows users to create a flow circuit by a 'point and click' paradigm. It creates the GFSSP input file after the completion of the model building process. GFSSP's GUI provides the users a platform to build and run their models. It also allows post-processing of results. The network flow circuit is first built using three basic elements: boundary node, internal node, and branch.

Majumdar, A. K.↗