Search NASA⌕ Search

SEARCH · Search NASA

Results for “two-phase flow modeling”

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

Non-dimensional performance and safety parameters for heat pipes

The use of heat pipes in safety-critical systems such as nuclear microreactors dictates the development of generalized, practical, scalable performance and safety parameters. Traditional dimensional metrics, while informative, lack the universality required for comparative analysis across varying designs and operating regimes. Here, this work introduces a comprehensive set of non-dimensional parameters to characterize heat pipe performance and safety, including capillary performance, effective thermal conductivity, response time, exergetic efficiency, allowable temperature gradients, allowable rate of temperature change, priming coefficients, and factor of safety. A reference heat pipe design representative of microreactor applications was analyzed via the developed parameters using both traditional analytical models and Sockeye simulations under transient and steady-state conditions. Sodium, potassium, and water were evaluated as working fluids to demonstrate the applicability of the framework across a broad temperature range. The proposed non-dimensional parameters effectively captured key thermal-hydraulic behaviors and safety concerns, as was demonstrated via Sockeye simulations. This framework supports the development of design optimization strategies, operational protocols, and safety assurance practices for advanced reactor systems and other high-reliability applications.

42 - ENGINEERING↗

Determining the extent of potential fugitive fluid migration from geologic carbon storage in hydrocarbon-bearing reservoirs: Insights from one-dimensional numerical modeling

Numerical modeling of Geologic Carbon Sequestration in permeable reservoirs initially containing hydrocarbons is conducted using the multi-phase, multi-component thermohydrologic simulator TOGA (TOUGH Oil, Gas, Aqueous; TOUGH stands for Transport Of Unsaturated Groundwater and Heat), to determine how phase and composition of the original fluids influence the extent of the zone where upward fugitive fluid migration could potentially occur, denoted R f . The area within R f comprises regions of substantially elevated pressure and free-phase CO 2 saturation, where a breach in reservoir sealing capacity would lead to upward fugitive fluid migration. The model examines the conditions within the storage reservoir that could lead to fugitive flow, but does not model the fugitive flow itself. A one-dimensional radial model of the storage reservoir is used, and three initial phase conditions are considered: single-phase aqueous, two-phase gas-aqueous, and three-phase oil-gas-aqueous. Components that may be present are H 2 O, CO 2 , CH 4 , C 4 H 10 , and C 10 H 22 . The most important factors controlling Rf are (1) the initial gas-phase saturation within the reservoir, and (2) the lateral extent of multi-phase initial conditions, particularly CO 2 . The composition of liquid and gas phases has a secondary effect. The impact of reservoir depth, thickness, injection rate, and hydrologic properties are also briefly examined, with thickness (or equivalently injection rate) having the biggest effect. These results can help to understand important trends in potential response of CO 2 -EOR fields being considered for dedicated CO 2 storage.

CO₂ plume migration↗

A Kinetic Model of the Long-Term Corrosion of Glass-Ceramic Materials

Multiphase waste forms show promise for increased waste loading and for the ability to dispose of contaminated solid and particulate waste through direct densification. However, achieving predictive capability for long-term durability of multiphase waste forms, and thus assessing their possible deployment, requires expanding the current, limited knowledge base. Here, we describe the development of a corrosion model of a two-phase waste form consisting of crystals of known volume fraction embedded in a glass matrix. This model accounts for the dissolution of both the crystalline and glass phases as well as the hydration of the glass phase through an ion exchange reaction. Because of the large difference in solubility between the two phases, the reactive surface of the crystalline phase is a function of the extent of dissolution of the glass phase in this model. Model parameterization was performed using corrosion data, such as from single-pass flow-through tests, for the individual phases. The parameterized corrosion model was evaluated against static dissolution test data for a glass-ceramic multiphase waste form. This evaluation demonstrated the model’s ability to reproduce the time-dependent release of key tracers of glass and crystalline phase dissolution. Hence, the development of a kinetic model provides a pathway for long-term durability predictions and thus the use of multiphase waste forms in nuclear cleanup missions.

Kerisit, Sebastien N.↗

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)↗

zmodel-kokkos

The code implements the so-called ZModel on an unstructured mesh. This models two-phase (or feasibly multi-phase) incompressible flow by meshing and evolving only the interface between materials, and never directly discretizing the bulk flow.

May, Ian↗

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.↗

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.↗

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.↗

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.↗

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.↗

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.↗

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.↗

A pumped two-phase cooling system for spacecraft

A pumped, two-phase heat-transport system is being developed for possible use for temperature control of scientific instruments on future NASA missions. As compared to a single-phase system, this two-phase system can maintain tighter temperature control with less pumping power. A laboratory model of the system has been built and tested. The measured heat transfer coefficients were approximately the same as in heat pipes, 220 Btu/hr-sq ft-F, as compared to 25 Btu/hr-sq ft-F for single-phase liquid flow. Heat shearing between experiments has been demonstrated wherein vapor generated in the cold plate of an active experiment was condensed in a cold, unheated experiment. System stability has been observed. However, additional development is needed. The use of non-azeotropic mixtures of coolants appears especially promising as a simple way to determine exit quality and thus control the flow rates to prevent dryout.

Ollendorf, S.↗

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.↗

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.↗

SSME fuel preburner phase 1/phase 2+ comparison

An analytical comparison of the predicted combustion flowfields for the Phase 1 and Phase 2+ fuel preburner configurations using an advanced CFD combustion code is currently underway at Rocketdyne. The Phase 2+ injector has a modified element design. The redesigned element provides for improved atomization through a shift in the injector element/face bleed flow split while maintaining an acceptable pressure drop. Cold-flow tests of the new design have been encouraging. As a complement to those tests, the Advanced Rocket Injector Combustor Code (ARICC) was selected to conduct a detailed combustion flowfield simulation. ARICC models two-dimensional (or axisymmetric), transient, turbulent, two-phase mixing and combustion flowfields. It is derived from a member of the Los Alamos ICEd-ALE family of codes. Unique features of ARICC include explicit representation of the coaxial LOX jet and hydrogen gas flows, distributed atomization processes, droplet breakup, and supercritical vaporation processes. Predictions of temperature and OH concentration profiles for the Phase 1 and Phase 2+ injector flowfields indicate a trend toward more uniform temperature distributions and shorter flame lengths with the Phase 2+ design.

Jensen, R. J.↗

Numerical analysis of confined recirculating gas-solid turbulent flows

Particle dispersion in confined recirculating turbulent flows has been investigated numerically. The present computational scheme utilizes Eulerian forms of the governing equations and allows two-way coupling between gas and solid phases. A recently developed two-phase closure model based on the multiple-scale turbulence model was used for the estimation of turbulent viscosities and diffusivities. For the particle size and loading considered in this study, the effect of particulate phase on the conveying gas is small, the nonequilibrium between the two phases is significant. Effects of recirculation, expanded chamber size and secondary annular jet momentum on the particle mixing rate are also investigated. In general, the present numerical results are in reasonably good agreement with the available experimental data.

Chen, C. P.↗

Liquid droplet heat exchanger studies

Recent analytical and experimental investigations of the liquid droplet heat exchanger (LDHX) concept for space power applications are described. The performance of the LDHX is compared to that of a conventional heat exchanger for heat rejection applications in a Brayton cycle, using the mass-specific heat exchanger effectiveness as a figure of merit. It is shown that the LDHX has an order of magnitude advantage over the conventional heat exchanger. Furthermore, significant improvement in cycle efficiency and power to mass ratio is possible. Two-phase flow experiments in a laboratory scale LDHX, using air and water as the two media, show very good agreement with the quasi-one-dimensional model used in the parametric studies.

Bruckner, A. P.↗