Search NASA⌕ Search

SEARCH · Search NASA

Results for “two-phase models”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Combining Distributed Electrochemical Measurements with a Semi-Empirical Model to Identify Local Variations in Overpotentials for PEM Water Electrolyzers

Abstract A segmented cell and PCB approach was used to demonstrate an advanced technique for measuring and understanding the electrochemical behavior of PEM water electrolyzers. By employing a semi-empirical model, discrete modes of overpotential and their spatial distributions were identified in an operating cell. Distributed area-specific resistance (DASR measurements indicate that local dehydration of the polymer electrolyte membrane is the primary driver of increased polarization, behaving similarly to mass transport-limited currents at low flow rates due to ineffective gas-phase bubble removal. To improve quantification, the overall ohmic overpotential is divided into background (η Ohm,b ) and mass transport-coupled (η Ohm,mt ) components. Two types of anode diffusion media, a baseline PTL and a thin foil liquid/gas diffusion layer (LGDL), were examined. The results show significant variations in η Ohm,b and highlight that ineffective bubble removal leads to pronounced local dehydration. The effect of cell temperature on various overpotentials was also examined. Negligible differences in diffusion overpotential were observed when comparing PEWE performance using the PTL and LGDL, despite major differences in two-phase fluid transport within each diffusion media. This underscores the critical role of the catalyst layer's micro-porous structure in diffusion processes. The study highlights the complex nature and intricacies of mass transport limitations within PEWEs.

Roy, Anirban (ORCID:0000000163803386)↗
DOE OSTI2026

Deflagration to Detonation Transition Update: XDDT Code Modularization

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

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
NASA NTRS1972

Theoretical analysis of a rotating two-phase detonation in liquid rocket motors.

A nonlinear analysis to study tangential mode shock instabilities in a thin annular chamber is carried out by employing a one dimensional two phase detonation wave as a reaction model. It is assumed that phase change and reaction take place only within the wave, which is treated as a discontinuity. The annulus is unrolled and the flow is considered as two dimensional with the coordinate system fixed on the wave front. Between waves, the flow is assumed to be isentropic with no interaction between droplets and burned gases. Jump conditions across the wave are solved for two phase flow. The average pressure along the injection plate is related to the design chamber pressure by the use of overall conservation equations. The wave strength is written in terms of the design parameters of the chamber. The results compare favorably with existing experiments. Finally, the effects of drop size on the wave strength are discussed and a simple criterion which sets the lower limit of validity for this strong wave analysis, is presented.

Shen, P. I.-W.↗
NASA NTRS1976

A theoretical investigation of thermodynamic effects on developed cavitation

The results of a theoretical investigation of thermodynamic effects on developed cavitation are presented. An approximate solution to the conservation equations for a two-phase laminar boundary layer is obtained. This analysis produces an expression for the temperature difference between the liquid and vapor phases which can be applied to developed cavity flows. Experimental data of cavity temperature depressions are correlated using this result. In addition, a theoretical estimate of the Nusselt number for the cavity is made using a turbulent boundary layer cavity model proposed by Brennen. The result agrees in part with empirically determined expressions for the cavity Nusselt number.

Weir, D. S.↗
NASA NTRS1976

Non-equilibrium radiation from viscous chemically reacting two-phase exhaust plumes

A knowledge of the structure of the rocket exhaust plumes is necessary to solve problems involving plume signatures, base heating, plume/surface interactions, etc. An algorithm is presented which treats the viscous flow of multiphase chemically reacting fluids in a two-dimensional or axisymmetric supersonic flow field. The gas-particle flow solution is fully coupled with the chemical kinetics calculated using an implicit scheme to calculate chemical production rates. Viscous effects include chemical species diffusion with the viscosity coefficient calculated using a two-equation turbulent kinetic energy model.

Penny, M. M.↗
NASA NTRS1976

Diffusion analysis for two-phase metal-matrix composite

Diffusion controlled filament matrix interaction in a metal matrix composite, where the filaments and matrix comprise a two phase binary alloy system, was mathematically modeled. The problem of a diffusion controlled, two phase moving interface by means of a one dimensional, variable grid, finite difference technique was analyzed. Concentration dependent diffusion coefficients and equilibrium solubility limits were used, and the change in filament diameter and compositional changes in the matrix were calculated as a function of exposure time at elevated temperatures. With the tungsten nickel (W-Ni) system as a model composite system, unidirectional composites containing from 0.06 to 0.44 initial filament volume fraction were modeled. Compositional changes in the matrix were calculated by superposition of the contributions from neighboring filaments. Alternate methods for determining compositional changes between first and second nearest neighbor filaments were also considered. The results show the relative importance of filament volume fraction, filament diameter, exposure temperature, and exposure time as they affect the rate and extent of filament matrix interaction.

Tenney, D. R.↗
NASA NTRS1977

FLUID: A numerical interpolation procedure for obtaining thermodynamic and transport properties of fluids

A computer program subroutine, FLUID, was developed to calculate thermodynamic and transport properties of pure fluid substances. It provides for determining the thermodynamic state from assigned values for temperature-density, pressure-density, temperature-pressure, pressure-entropy, or pressure-enthalpy. Liquid or two-phase (liquid-gas) conditions are considered as well as the gas phase. A van der Waals model is used to obtain approximate state values; these values are then corrected for real gas effects by model-correction factors obtained from tables based on experimental data. Saturation conditions, specific heat, entropy, and enthalpy data are included in the tables for each gas. Since these tables are external to the FLUID subroutine itself, FLUID can implement any gas for which a set of tables has been generated. (A setup phase is used to establish pointers dynamically to the tables for a specific gas.) Data-table preparation is described. FLUID is available in both SFTRAN and FORTRAN

Fessler, T. E.↗
NASA NTRS1977

An information system design for watershed-wide modeling of water loss to the atmosphere using remote sensing techniques

Results are presented of a study intended to develop a general location-specific remote-sensing procedure for watershed-wide estimation of water loss to the atmosphere by evaporation and transpiration. The general approach involves a stepwise sequence of required information definition (input data), appropriate sample design, mathematical modeling, and evaluation of results. More specifically, the remote sensing-aided system developed to evaluate evapotranspiration employs a basic two-stage two-phase sample of three information resolution levels. Based on the discussed design, documentation, and feasibility analysis to yield timely, relatively accurate, and cost-effective evapotranspiration estimates on a watershed or subwatershed basis, work is now proceeding to implement this remote sensing-aided system.

Khorram, S.↗
NASA NTRS1977

Assessment of analytical techniques for predicting solid propellant exhaust plumes

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

Tevepaugh, J. A.↗
NASA NTRS1979

Two-phase choked flow of cryogenic fluids in converging-diverging nozzles

Data are presented for the two phase choked flow of three cryogenic fluids - nitrogen, methane, and hydrogen - in four converging-diverging nozzles. The data cover a range of inlet stagnation conditions, all single phase, from well below to well above the thermodynamic critical conditions. In almost all cases the nozzle throat conditions were two phase. The results indicate that the choked flow rates were not very sensitive to nozzle geometry. However, the axial pressure profiles, especially the throat pressure and the point of vaporization, were very sensitive to both nozzle geometry and operating conditions. A modified Henry-Fauske model correlated all the choked flow rate data to within + or - 10 percent. Neither the equilibrium model nor the Henry-Fauske model predicted throat pressures well over the whole range of data. Above the thermodynamic critical temperature the homogeneous equilibrium model was preferred for both flow rate and pressure ratio. The data of the three fluids could be normalized by the principle of corresponding states.

Simoneau, R. J.↗
NASA NTRS1982

The photochemistry and dynamics of a dusty cometary atmosphere

The solving of a simultaneous set of differential equations representing conservation of number density, momentum, and energy together with solar radiation transfer in the streams which result in photolytic processes and the heating of the nucleus yields a self-consistent solution of the dynamical and thermal structure of an H2O-dominated two-phase dusty gas cometary atmosphere. Two models are considered for the transfer of solar radiation through the circumnuclear dust halo. The first considers only the direct extinction by the dust, and in the second the diffuse radiation field due to multiple scattering by the dust halo, which compensates for radiation removed by direct absorption when the optical depth is near unity, is approximated by neglecting the attenuation of radiation given off by the dust. It is shown that while dust attenuation has a strong effect on the H2O production rate, it also increases the electron density in the inner coma over the unattenuated case.

Marconi, M. L.↗
NASA NTRS1982

A theory of electrophoresis of emulsion drops in aqueous two-phase polymer systems

An electrophoresis study has been carried out in an emulsion formed from an electrically neutral aqueous mixture of dextran and polyethylene glycol equilibrated at sufficient concentrations in the presence of electrolytes. Electrophoresis of a drop of one phase suspended in the other is observed, and the direction of the drop's motion is reversed when the disperse phase and the continuous phase are interchanged. In the presence of sulfate, phosphate, or citrate ions, an electrostatic potential difference of the order of a few mV exists between the two phases. The potential implied by the direction of the electrophoretic motion is opposite to the Donnan potential observed between the two phases. The mobility of an emulsion drop increases with the drop radius and depends on ion concentration. These results are explained in terms of a model postulating an electric dipole layer associated with a mixture of oriented polymer molecules at the surface of a drop, with a potential difference between the interiors of the two phases resulting from the unequal ion distribution.

Levine, S.↗
NASA NTRS1983

Modeling of transient two-component flow using a four-point implicit method

The four-point, centered implicit scheme that is extensively used in open channel flow simulation is shown to be applicable to rapid and slow pressure transient problems in conduits with nearly single phase and two-phase flows. It is only necessary to choose the proper weighting factor value, theta, of the Courant number. For rapid pressure transients such as waterhammer, the implicit method can yield reasonable results with limited numerical dispersion and attenuation if theta is only slightly greater than the critical value of 0.5. For slower pressure gradients in single and two-phase flows, reasonable numerical solutions may be achieved for Courant number values as high as 20.

Wiggert, D. C.↗
NASA NTRS1984

Dynamics of two-phase face seals

An analytic study is presented of the effects of phase change on load support for parallel and tapered face seals. Consideration is given to an adiabatic model for low Reynolds number flow. Numerical integration is carried out of the descriptive fluid equations, giving the opening force due to fluid film pressure. The loci of steady-state solutions are then plotted for water to provide curves of load support as a function of film thickness. For axial excursions of the seal rings, a quasi-steady transient analysis is made. It is found that the load support generated by fluid pressure can be multivalued for a given film thickness. Another finding is that axial disturbances of the seal rings may lead to sudden drops in load support generated by fluid pressure with three possible results. The first is that sufficient damping may permit the seal to return to the previous equilibrium operating position. The second is that the seal may collapse to an equilibrium position of smaller film thickness where face contact is more likely and a significantly higher operating temperature is assured. The third is that a limit cycle of self-sustained oscillation in the axial direction may occur if damping is sufficiently low.

Beeler, R. M.↗
NASA NTRS1985

Comprehensive modeling of a liquid rocket combustion chamber

An analytical model for the simulation of detailed three-phase combustion flows inside a liquid rocket combustion chamber is presented. The three phases involved are: a multispecies gaseous phase, an incompressible liquid phase, and a particulate droplet phase. The gas and liquid phases are continuum described in an Eulerian fashion. A two-phase solution capability for these continuum media is obtained through a marriage of the Implicit Continuous Eulerian (ICE) technique and the fractional Volume of Fluid (VOF) free surface description method. On the other hand, the particulate phase is given a discrete treatment and described in a Lagrangian fashion. All three phases are hence treated rigorously. Semi-empirical physical models are used to describe all interphase coupling terms as well as the chemistry among gaseous components. Sample calculations using the model are given. The results show promising application to truly comprehensive modeling of complex liquid-fueled engine systems.

Liang, P.-Y.↗
NASA NTRS1985

Two-phase ultraviolet spectrophotometry of the pulsating white dwarf ZZ Piscium

Spectra of the pulsating white dwarf ZZ Psc (= G29-38) were obtained using the International Ultraviolet Explorer. By using a multiple-exposure technique in conjunction with simultaneous ground-based exposure-metering photometry, it was possible to obtain mean on-pulse and off-pulse spectra in the 1950-1310 A wavelength range. The ratio of the time-averaged on-pulse to off-pulse spectra is best fitted by a temperature variation that is in phase with the optical light variation. This result is consistent with the hypothesis that the observed variation is due to a high-order nonradial pulsation. Conventional ultraviolet spectra of ZZ Psc showed broad absorption features at 1390 and 1600 A. These features are also found in the spectra of the cool DA-type white dwarfs G226-29 and G67-23, and appear to increase in strength with decreasing temperature. A possible explanation for the 1600 A feature is absorption by the satellite band of resonance-broadened hydrogen Ly-alpha. Such absorption would also help explain a discrepancy between the observed pulsation amplitude shortward of 1650 A and the predicted amplitudes based on model atmospheres.

Bond, H. E.↗
NASA NTRS1985

High thermal capacity cold plate/hot plate

The results of an analytical study to determine the feasibility of a novel two-phase cold plate/hot plate (CPHP) are presented. A key feature of the CPHP is the use of capillary forces to separate the liquid and vapor phases and distribute the liquid over the evaporating/condensing surface. The liquid phase is held by capillary forces in a reservoir and is carried to the evaporating surface by a wick. The reservoir is replenished at intervals from a valved external liquid supply line. In the hot plate mode, liquid accumulates in the reservoir and is removed by an external condensate line. Performance requirements for the device were capability of handling a power density of 4 watts/sq cm, an outlet quality (percentage of vapor flow to total flow) greater than 90 percent, and operation in a 0-g environment. An analytical model of CPHP operation was developed which concentrated on the liquid and vapor flows. It was found that the liquid carrying capillary grooves of rectangular cross-section gave significantly better predicted performance than V-shaped grooves. Using the analytical model, capillary groove width and length, vapor channel dimensions, and other parameters were selected to meet the performance requirements with a 100-percent margin.

Hwangbo, H.↗