Bubble formation and growth - Study of the boundary conditions at a liquid-vapor interface through irreversible thermodynamics Final report, Feb. 1966 - Feb. 1967
Irreversible thermodynamics of liquid-vapor phase change process
SEARCH · Search NASA
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.
Irreversible thermodynamics of liquid-vapor phase change process
Theoretical and experimental approaches to exhaust plume analysis are compared. A two-phase model is extended to include treatment of reacting gas chemistry, and thermodynamical modeling of the gaseous phase of the flow field is considered. The applicability of empirical data currently available to define particle drag coefficients, heat transfer coefficients, mean particle size, and particle size distributions is investigated. Experimental and analytical comparisons are presented for subscale solid rocket motors operating at three altitudes with attention to pitot total pressure and stagnation point heating rate measurements. The mathematical treatment input requirements are explained. The two-phase flow field solution adequately predicts gasdynamic properties in the inviscid portion of two-phase exhaust plumes. It is found that prediction of exhaust plume gas pressures requires an adequate model of flow field dynamics.
Homogeneous nucleation kinetics are rederived in a manner fundamentally similar to the approach of classical nucleation theory with the following modifications and improvements. First, the cluster is a parent phase cluster and does not require energization to the parent state. Second, the thermodynamic potential used to describe phase stability is a continuous function along the pathway of phase decomposition. Third, the kinetics of clustering corresponds directly to the diffusional flux of monomers through the cluster distribution and are formally similar to classical theory with the resulting kinetic equation modified by two terms in the preexponential factor. These terms correct for the influence of a supersaturation dependent clustering within the parent phase and for the influence of an asymmetrical cluster concentration as a function of cluster size at the critical cluster size. Fourth, the supersaturation dependence of the nucleation rate is of the same form as that given by classical nucleation theory. This supersaturation dependence must however be interpreted in terms of a size dependent surface tension. Finally, there are two scaling laws which describe supersaturation to either constant nucleation rate or to the thermodynamically determined physical spinodal.
Study of thermodynamic phenomena occurring during transfer of cryogenic liquids from dewar to receiver tank reveals that the basic cause of tank implosion is evaporation rate of droplets entering the tank in the early transfer phase. Analyses of the thermodynamics involved and implosion prevention techniques are included.
Computer program subroutine FLUID calculates thermodynamic and transport properties of pure fluids in liquid, gas, or two-phase (liquid/gas) conditions. Program determines thermodynamic state from assigned values for temperature and density, pressure and density, temperature and pressure, pressure and entropy, or pressure and enthalpy.
The primary challenge of the X-33 CAN is to build and test a prototype LO2 and LH2 densification ground support equipment (GSE) unit, and perform tank thermodynamic testing within the 15 month phase 1 period. The LO2 and LH2 propellant densification system will be scaled for the IPTD LO2 and LH2 tank configurations. The IPTD tanks were selected for the propellant technology demonstration because of the potential benefits to the phase 1 plan: tanks will be built in time to support thermodynamic testing; minimum cost; minimum schedule risk; future testing at MSFC will build on phase 1 data base; and densification system will be available to support X-33 and RLV engine test at IPTD. The objective of the task 1 effort is to define the preliminary requirements of the propellant densification GSE and tank recirculation system. The key densification system design parameters to be established in Task 1 are: recirculation flow rate; heat exchanger inlet temperature; heat exchanger outlet temperature; maximum heat rejection rate; vent flow rate (GN2 and GH2); densification time; and tank pressure level.
Rare earth oxide stabilized zirconias find a wide range of applications. An understanding of phase equilibria is essential to all applications. In this study, the available phase boundary data and thermodynamic data is collected and assessed. Calphad-type databases are developed to completely describe the Y2O3-ZrO2, Yb2O3-ZrO2, and Y2O3-Yb2O3 systems. The oxide units are treated as components and regular and subregular solution models are used. The resultant calculated phase diagrams show good agreement with the experimental data. Then the binaries are combined to form the database for the Y2O3-Yb2O3-ZrO2 psuedo-ternary.
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.
We present a formal thermodynamic treatment of superfluid flow in a Josephson junction. We show that the current and the phase difference are thermodynamic conjugate variables. We derive quantitative expressions for the rms fluctuations of these variables. Also, we discuss the thermodynamic stability and the thermal activation to the phase slip region. We apply the developed formalism to show why an array of apertures in He-4 can exhibit the Josephson effect near the Lambda transition despite strong thermal fluctuations.
Constraints on the clustered mass density of the universe derived from the observed population mean intracluster gas fraction of x-ray clusters may be biased by reliance on a single-phase assumption for the thermodynamic structure of the intracluster medium (ICM). We propose a descriptive model for multiphase structure in which a spherically symmetric ICM contains isobaric density perturbations with a radially dependent variance. Fixing the x-ray emission and emission weighted temperature, we explore two independently observable signatures of the model in the parameter space. For bremsstrahlung dominated emission, the central Sunyaev-Zel'dovich (SZ) decrement in the multiphase case is increased over the single-phase case and multiphase x-ray spectra in the range 0.1-20 keV are flatter in the continuum and exhibit stronger low energy emission lines than their single-phase counterpart. We quantify these effects for a fiducial 10e8 K cluster and demonstrate how the combination of SZ and x-ray spectroscopy can be used to identify a preferred location in the plane of the model parameter space. From these parameters the correct value of mean intracluster gas fraction in the multiphase model results, allowing an unbiased estimate of clustered mass density to he recovered.
The steady motion, thermal and free surface behavior of a volatile, wetting liquid in microgravity are studied using scaling and numerical techniques. The objective is to determine whether the thermocapillary and two-phase convection arising from thermodynamic nonequilibrium along the porous surfaces of spacecraft liquid acquisition devices could cause the retention failures observed with liquid hydrogen and heated vapor pressurant. Why these devices seem immune to retention loss when pressurized with heated helium or heated directly through the porous structure was also examined. Results show that highly wetting fluids exhibit large negative and positive dynamic pressure gradients towards the meniscus interline when superheated and subcooled, respectively. With superheating, the pressure variation and recoil force arising from liquid/vapor phase change exert the same influence on surface morphology and promote retention. With subcooling, however, the pressure distribution produces a suction that degrades mechanical equilibrium of the surface. This result indicates that thermocapillary-induced deformation arising from subcooling and condensation is the likely cause for retention loss. In addition, increasing the level of nonequilibrium by reducing accommodation coefficient suppresses deformation and explains why this failure mode does not occur in instances of direct screen heating or pressurization with a heated inert gas.
In systems where the design inlet and outlet pressures P sub amb are maintained above the thermodynamic critical pressure P sub c, it is often assumed that heat and mass transfer are governed by single-phase relations and that two-phase flows cannot occur. This simple rule of thumb is adequate in many low-power designs but is inadequate for high-performance turbomachines, boilers, and other systems where two-phase regions can exist even though P sub amb P sub c. Heat and mass transfer and rotordynamic-fluid-mechanic restoring forces depend on momentum differences, and those for a two-phase zone can differ significantly from those for a single-phase zone. By using a laminar, variable-property bearing code and a rotating boiler code, pressure and temperature surfaces were determined that illustrate nesting of a two-phase region within a supercritical pressure region. The method of corresponding states is applied to bearings with reasonable rapport.
In systems where the design inlet and outlet pressure P sub amb are maintained above the thermodynamic critical pressure P sub c, it is often assumed that heat and mass transfer are governed by single-phase relations and that two-phase flows cannot occur. This simple rule of thumb is adequate in many low-power designs but is inadequate for high-performance turbomachines, boilers, and other systems where two-phase regions can exist even though P sub amb P sub c. Heat and mass transfer and rotordynamic-fluid-mechanic restoring forces depend on momentum differences, and those for a two-phase zone can differ significantly from those for a single-phase zone. By using a laminar, variable-property bearing code and a rotating boiler code, pressure and temperature surfaces were determined that illustrate nesting of a two-phase region within a supercritical pressure region. The method of corresponding states is applied to bearings with reasonable rapport.
Nonlinear equations governing complex chemical and phase equilibria derived from thermodynamic principles and supplemented by solving algorithms and practical applications
High temperature thermodynamics requiring species and phases identification, crystal structures, molecular geometries and vibrational, rotational and electronic energy levels and equilibrium constants
Power supply utilizes small quantities of two-phase fluid of suitable thermodynamic properties for circulation in a capillary-pumped heat transfer loop. Fluid is vaporized in evaporator, passed through fluidic system load, condensed, pumped by multistage capillary pump, and returned to the evaporator.
The sputtering and ion plating technology is reviewed in terms of potential, uses, and performance. Sputtering is not regulated by classical thermodynamics and Gibb's phase rule relationships, thus eliminating the restrictions for materials combinations, and consequently making possible tailoring of coatings in any preferred chemical combination. The future of sputtered and ion plated films for industrial applications is expected to depend primarily on greater comprehension of materials selection, and the utilization of the proper deposition parameters.
The primary phases of Al-rich inclusions in amoeboid olivine aggregates have undergone alteration reactions with the solar nebular gas. The simplest interpretation of the present observations is that melilite was the first primary phase to disappear with falling temperature, and was replaced by grossular + anorthite + feldspathoids, followed by fassaite; spinel was the last phase to be altered. Thermodynamic calculations suggest that Na-rich phlogopite could have formed at about 470 K and chlorite at about 328 K at a water fugacity of 0.000001, which is that of a gas of solar composition in this temperature range. The olivine around Al-rich inclusions is not serpentized, indicating the cessation of gas-solid equilibrium above 274 K.