Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “two-phase bubbly flow”

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 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS1979

Marangoni bubble motion in zero gravity

It is shown experimentally that the Marangoni phenomenon is a primary mechanism for the movement of a gas bubble in a nonisothermal liquid in a low-gravity environment. In such two-phase flow systems, local variations in bubble surface tension are caused by a temperature gradient in the liquid. Shearing stresses thus generated at the bubble surface lead to convection in both media, as a result of which the bubble begins to move. A mathematical model consisting of the Navier-Stokes equations and the thermal energy equations, along with the appropriate boundary conditions for both media, is proposed.

Thompson, R. L.↗
DOE OSTI2025

Advanced Computational Modeling of High-Level Waste Vitrification at the Hanford Site

The U.S. Department of Energy (DOE) has selected vitrification for stabilizing legacy tank waste at the Hanford site, where radioactive waste from plutonium production was historically stored in underground tanks. This waste will be separated into low-activity waste (LAW) and high-level waste (HLW) fractions and processed at the Waste Treatment and Immobilization Plant (WTP). At WTP, glass melters are used for the vitrification of radioactive tank waste, transforming it into a stable borosilicate glass form for safe long-term storage. The melter vessel is constructed from highly durable and heat-resistant materials, where the vitrification process occurs. The main regions that are modeled are the melt pool, plenum, cold cap, riser/discharge chamber, and surrounding structure with insulation layers. Forced convection induced by air bubblers at the base of the melter ensure uniform temperature distribution and provide heat to the cold cap layer. The cold cap is a region of reacting batch feed that floats on top of the molten glass and is where the batch-to-glass reactions occur. Joule heating provided by electrodes mounted along the vertical walls of the melter and immersed directly in the glass, generates the necessary heat for the net endothermic conversion processes that occur in the cold cap. The high temperatures, radioactivity, and opaque nature of the glass prevent direct observation inside the melters. Therefore, computational models are essential for providing insight into factors that affect melter throughput. Thermocouples in the plenum provide operators with plenum temperature measurements. Operational adjustments include bubbling rate, voltage supplied to the electrodes, feed adjustments, and glass removal rate. Different computational fluid dynamics (CFD) models have been developed, each serving a specific purpose. There are CFD models of different scale melters, as well as models that capture the two-phase flow interfaces of rising bubbles in the molten glass or models with a simplified molten glass region so that the surrounding structure and plenum can be feasibly incorporated. Pilot-scale melter models have been developed to serve as validation of the methods employed in the simulation of the full-scale WTP melters. Models incorporating resolved bubbling are used to develop momentum source terms to implement into a single phase, multi-region, steady-state flow model that is being validated by measured process parameters such as glass production rate, voltage, input power, plenum temperatures, etc. The resolved bubbling model uses the multiphase volume of fluid approach to model the system with a high-resolution interface capturing scheme to maintain sharp interfaces between the molten glass and the air phase. The suite of CFD models is continually being improved to incorporate more realistic physics and achieve faster turnaround time. For example, an incremental controller is implemented to automatically adjust electrode voltage within the simulation to a molten glass set point temperature of 1150°C. Newer models feature improved meshes to ensure conformal meshes between regions and eliminate unnecessary mesh refinement in areas that are not of interest (such as boundary layers in offgas ports). Instead of explicitly modeling the structural, refractory, and insulation layers of the melter, a thermal resistance approach is used with published correlations used for boundary conditions. The development of robust and efficient CFD models will be instrumental in enabling the WTP to successfully fulfill its mission of safely stabilizing legacy nuclear waste.

12 - MGMT OF RADIOACTIVE AND NON-RADIOACTIVE WASTE↗
NASA NTRS1977

Experimental study of the structure of isotropic turbulence with intermediate range of Reynolds number

The energetic isotropic turbulence generated by a waterfall of low head was found to be developed in part through the unstable two-phase flow of entrained air bubbles. The resulting turbulent field had a turbulent Reynolds number in excess of 20,000 and maintained a self-similar structure throughout the decay period studied. The present study may provide some insight into the structure of turbulence produced by breaking waves over the ocean.

Ling, S. C.↗
DOE OSTI2025

MOSCATO Development and Integration in Fiscal Year 2025: Implementation of Multiphase, Multiphysics Modeling Capabilities for Molten Salt Systems

MOSCATO (Molten Salt Chemistry and Transport) is a multiphysics code that provides high-fidelity, coupled simulations of fluid flow, heat transfer, mass transfer, chemistry, electrochemical phenomena, and alloy corrosion for molten salt systems. In FY25, significant developments were made to the code package, enhancing its capabilities for modeling all relevant phenomena within operating moltens salt reactors (MSRs). The developments and activities in FY25 included: 1. Implementation of Level-Set methods to enable modeling of single-bubble behavior in molten salts. In FY25, the Level-Set two-phase flow modeling implementation was improved to simulate single bubble behavior with molten salt media. The large density and viscosity ratios between typical gases and molten salt liquids present challenges for these types of numerical solvers. With enhancements to the pressure projection method, MOSCATO’s Level-Set solver was able to be successfully validated to experiments related to helium bubble rise in stagnant molten salt. The simulated bubble rising velocity showed reasonable good agreement with experimental measurements. The bubble shape and dynamics were also visually compared with experimental snapshots, demonstrating a good qualitative match. 2. Generation of mass transfer correlations for multiphase flow systems. To enable calculations of the tritium transport across the interface between gas bubbles and salt, we modeled high- Schmidt-number mass transfer around a sphere across a broad range of Reynolds numbers. The mesh near the sphere surface was highly refined to resolve steep concentration gradients caused by the low diffusion coefficient. Literature-based mass transfer correlations were compared with the numerical results, and modifications were proposed to improve agreement, particularly at higher Schmidt numbers. These mass transfer correlations were subsequently provided to other national laboratories to help enable high quality mass transfer simulations using lower-order solvers under development within the NEAMS program. 3. Preliminary implementation of a bubbly flow solver. To model bubbly flow in molten salt, we implemented a bubbly flow solver for void fractions less than 5%. To do so, an algebraic relative velocity model that assumes small bubbles with rapid momentum equilibration was added to MOSCATO to compute bubble velocities. Preliminary comparisons with experimental data showed reasonable agreement, and further improvements are underway. 4. Generation of mass transfer correlations for MSRE subchannel The Molten-Salt Reactor Experiment (MSRE) was a landmark historical project that demonstrated the feasibility of molten-salt reactor technology. The MSRE campaign also generated a significant body of experimental data and reports that continue to support molten-salt–related research. In this report, the MSRE core subchannel was used as the reference geometry for a mass transfer study performed with MOSCATO. The geometry and computational mesh were adapted from a previous study, providing adequate resolution for the relatively low Reynolds number in this case. Additional mesh refinement was applied to reach higher Schmidt numbers, enabling the derivation of a reliable mass-transfer correlation for the present scenario. 5. Simulations of oxygen ingressions into molten salt. In the previous fiscal year, we initiated a study to simulate oxygen ingression in stagnant salt. As oxygen enters the salt through its surface, it reacts with Ce 3+ to form solid CeO 2 and other reaction products. To more fully capture the complex diffusion-convection-reaction mechanisms, capabilities for modeling natural convection in the salt vessel were added. These were needed as the flow of the ingressed gas induced flow in the salt caused by surface shear and non-isothermal effects. With these updated physics in place, we were able to successfully reproduce the experimental results for the rate of change of CeCl 3 concentrations versus time. 6. Flow corrosion model validation. In FY24, MOSCATO’s corrosion model was validated against static corrosion experiments. In FY25, this work was extended to a flow corrosion experiment, where FLiNaK salt was driven by natural convection, with initial salt impurities to initiate corrosion. Despite uncertainties in parameters such as elemental diffusion coefficients in the alloy and unknown H + concentrations, the simulations achieved good agreement with experimental results, especially in predicting sample mass losses.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1975

Liquid jet pumped by rising gas bubbles

A two-phase mathematical model is proposed for calculating the induced turbulent vertical liquid flow. Bubbles provide a large buoyancy force and the associated drag on the liquid moves the liquid upward. The liquid pumped upward consists of the bubble wakes and the liquid brought into the jet region by turbulent entrainment. The expansion of the gas bubbles as they rise through the liquid is taken into account. The continuity and momentum equations are solved numerically for an axisymmetric air jet submerged in water. Water pumping rates are obtained as a function of air flow rate and depth of submergence. Comparisons are made with limited experimental information in the literature.

Hussain, N. A.↗
NASA NTRS1975

Liquid jet pumped by rising gas bubbles

From observations of a stream of gas bubbles rising through a liquid, a two-phase mathematical model is proposed for calculating the induced turbulent vertical liquid flow. The bubbles provide a large buoyancy force and the associated drag on the liquid moves the liquid upward. The liquid pumped upward consists of the bubble wakes and the liquid brought into the jet region by turbulent entrainment. The expansion of the gas bubbles as they arise through the liquid is taken into account. The continuity and momentum equations are solved numerically for an axisymmetric air jet submerged in water. Water pumping rates are obtained as a function of air flow rate and depth of submergence. Comparisons are made with limited experimental information in the literature.

Hussain, N. A.↗
DOE OSTI2025

Using X-ray radiography to study oxygen flow in a proton exchange membrane electrolyzer operating under balanced pressure conditions

Of the various water electrolyzer technologies, the proton exchange membrane electrolyzer (PEMWE) is one of the best solutions for producing clean hydrogen without releasing CO 2 . In order to allow for widespread use of clean hydrogen, it is necessary to decrease its cost, which is intrinsically related to system operation. Current PEMWE plants operate in differential mode, directly pressurizing hydrogen and benefiting from thermodynamic compression, which increases overall system efficiency. However, high differential pressure above 30 bar can cause membrane stress, resulting in membrane creeping and failure. Pressurizing the water and operating at balanced pressure allows hydrogen to be produced at higher pressures while preserving the integrity of the membrane and porous layers. Nevertheless, the impact of pressurizing water on PEMWE performance must be better understood to maximize performance under balanced pressure conditions. Here, this study examined the impact of water pressure on two-phase flow. A high-pressure electrolyzer setup was developed to perform operando X-ray radiography and examine oxygen transport with high temporal resolution. The imaging segmentation process, developed to capture bubble properties in the channel, was applied to a specific experiment. The results clearly showed that as pressure increased up to 30 bars, the initial bubbly flow transitioned to slug flow, which led to channel saturation with oxygen. This work demonstrates that two-phase flow in an electrolyzer can be studied using X-ray radiography, which has the advantages of fast measurements and the ability to probe dense materials, such as those required for pressurized electrolyzers.

Balanced high pressure operation↗
NASA NTRS1986

Estimating Transient Pressure Surges in Cryogenic Systems

Potentially-damaging pressure waves anticipated and, therefore, avoided. Mathematical model developed for prediction of pressure behavior in single-and two-phase cryogenic systems. Transient liquid-flow analysis modified to incorporate behavior of vapor bubbles and used to predict maximum pressure in cryogenic transfer systems consisting of complex pipe and valve arrangements under both steady-state and transient conditions. Simulation compared favorably with data obtained during transfer of liquid oxygen from ground storage tanks to Space Shuttle orbiter external tanks. Program written in FORTRAN 77 for batch execution.

Pfister, P.↗
DOE OSTI2026

Two-phase flow numerical analysis of electrode geometry for alkaline water electrolyzers

Hydrogen is a promising component of a future energy-secure and efficient economy, but its competitiveness depends on reducing production costs. One strategy is to operate alkaline water electrolyzers at higher current densities to increase output. However, this intensifies performance losses due to gas bubble accumulation, which blocks transport pathways and deactivates electrochemically active surfaces. Enhancing bubble evacuation through electrode design is therefore essential. Previous studies have explored various approaches — such as modifying surface morphology, applying sonication or pressure modulation, and introducing surfactants — but these efforts have addressed a limited range of conditions due to the complexity of two-phase flow and electrode geometries. Experiments have also largely been focused on either cell level improvements, which lack the information necessary to isolate each contributing factor, or on modified geometries that are not relevant to practical cell operation. From a modeling perspective, conventional Eulerian multiphase models do not track the complex gas–liquid interfacial dynamics and often neglect surface tension and contact angle effects, reducing their predictive accuracy. To provide insights on the effects of different electrode geometries on the performance of alklaine water electrolyzers this work employs an immersed boundary volume-of-fluid method to simulate bubble behavior in 3D porous electrodes. Multiple base electrode geometries, typically used in practice, with varying porosity are evaluated under a constant surface gas generation rate. Simulation data is analyzed to quantify electrode gas coverage, bubble size dynamics and other relevant metrics. Results show that porosity strongly influences bubble accumulation on electrode surfaces, with higher porosity reducing gas coverage, and its not strictly dependent on the electrode geometry. However, the electrode’s base geometry significantly affects gas accumulation at the separator gap, independent of porosity. A foam electrode geometry resulted in the lowest gas coverage of all electrodes with a median volumetric gas coverage of 11%, but at the cost of a 70% reduction in active area compared with the largest surface area electrode, while gyroid electrodes showed the best trade-off between gas coverage, particularly at the separator surface, and electrochemically active area. In conclusion, the results highlight the need for holistic electrode design strategies.

Alkaline↗
NASA NTRS1983

Analogy between fluid cavitation and fracture mechanics

When the stresses imposed on a fluid are sufficiently large, rupture or cavitation can occur. Such conditions can exist in many two-phase flow applications, such as the choked flows, which can occur in seals and bearings. Nonspherical bubbles with large aspect ratios have been observed in fluids under rapid acceleration and high shear fields. These bubbles are geometrically similar to fracture surface patterns (Griffith crack model) existing in solids. Analogies between crack growth in solid and fluid cavitation are proposed and supported by analysis and observation (photographs). Healing phenomena (void condensation), well accepted in fluid mechanics, have been observed in some polymers and hypothesized in solid mechanics. By drawing on the strengths of the theories of solid mechanics and cavitation, a more complete unified theory can be developed.

Hendricks, R. C.↗
NASA NTRS2015

Numerical Modeling of the Transient Chilldown Process of a Cryogenic Propellant Transfer Line

Before cryogenic fuel depots can be fully realized, efficient methods with which to chill down the spacecraft transfer line and receiver tank are required. This paper presents numerical modeling of the chilldown of a liquid hydrogen tank-to-tank propellant transfer line using the Generalized Fluid System Simulation Program (GFSSP). To compare with data from recently concluded turbulent LH2 chill down experiments, seven different cases were run across a range of inlet liquid temperatures and mass flow rates. Both trickle and pulse chill down methods were simulated. The GFSSP model qualitatively matches external skin mounted temperature readings, but large differences are shown between measured and predicted internal stream temperatures. Discrepancies are attributed to the simplified model correlation used to compute two-phase flow boiling heat transfer. Flow visualization from testing shows that the initial bottoming out of skin mounted sensors corresponds to annular flow, but that considerable time is required for the stream sensor to achieve steady state as the system moves through annular, churn, and bubbly flow. The GFSSP model does adequately well in tracking trends in the data but further work is needed to refine the two-phase flow modeling to better match observed test data.

Chilldown↗
NASA NTRS2023

Cryogenic Flow Boiling in Microgravity: Effects of Reduced Gravity on Two-Phase Fluid Physics and Heat Transfer

With the growing interest in space exploration, cryogenic technologies involving two-phase flow and heat transfer are in high demand to successfully procure advanced space applications such as fuel depots and nuclear thermal propulsion (NTP) systems for deep space missions. However, the unique and extreme thermal properties of cryogenic fluids introduce distinct flow boiling fluid physics and energy transport phenomena, which differ significantly from those observed with conventional fluids. Understanding the unique two-phase physics in cryogenic flow boiling remains an ongoing challenge. Furthermore, the lack of readily available microgravity cryogenic steady-state heat transfer data hinders the assessment of gravitational effects on cryogenic flow boiling. This study aims to elucidate the gravitational effects on two-phase fluid physics and heat transfer by conducting the first-ever experimental measurement of cryogenic flow boiling performance using a steady-state heated method in a reduced gravity environment. Parabolic flight experiments were performed to acquire both heat transfer measurements and high-speed video of interfacial behaviors, under varying gravity levels (microgravity, hypergravity, Lunar gravity, and Martian gravity). The experiments involved flow boiling of liquid nitrogen (LN 2 ) with a near-saturated inlet along a circular heated tube of dimensions 8.5-mm inner diameter and 680-mm heated length. The operating parameters varied are mass velocity of 398.3 - 1342.8 kg/m2s, inlet quality of -0.08 to -0.01, and inlet pressure of 413.68 - 689.48 kPa. Captured microgravity flow patterns range from bubbly to annular, all having vapor structures that are larger than those under higher gravity levels. Under microgravity, absence of buoyancy yields symmetrical vapor structures without flow stratification, laying a physical foundation for the distinct two-phase heat transfer trends during LN 2 flow boiling in microgravity. Transient data collected during the flight parabolas exhibited decreasing heated wall temperature as the aircraft transitioned from hypergravity to microgravity phases. The temperature variation indicated an enhancement in flow boiling heat transfer with decreasing gravity levels and a reduction with increasing gravity levels. The effect of reduced gravity on cryogenic flow boiling heat transfer coefficient (HTC) is discussed based on steady state heat transfer analysis. Seminal HTC correlations are evaluated against the measured microgravity HTC data, of which one is identified for superior accuracy in predicting microgravity data. Finally, a new HTC correlation is proposed to improve accuracy of microgravity predictions, yet there still exists room for further improvement with future terrestrial flow boiling experiments at different flow orientations relative to Earth gravity.

Microgravity↗
NASA NTRS2019

Rewet Temperature Correlations for Liquid-Nitrogen Boiling Pipe Flows Across Varying Flow Conditions and Orientations

In many convective liquid-vapor phase change heat transfer engineering applications, cryogenic fluids are widely used in industrial processes, spacecraft and cryosurgery systems, and so on. For example, cryogens are usually used as liquid fuels such as liquid hydrogen and oxygen in the rocket industry, liquid nitrogen (LN2) and helium are frequently used to cool superconducting magnetic device for medical applications. In these systems, proper transport, handling, and storage of cryogenic fluids are of extreme importance. Among all the cryogenic transport processes performed in room temperatures, quenching, also termed chilldown, is a unavoidable initial, transient phase-change heat transfer process that brings the system down to the cryogenic condition. The Leidenfrost temperature or rewet temperature that signals the end of film boiling is practically considered the completion point of a quenching process. Therefore, rewet temperature has been considered the most important parameter for the engineering design of cryogenic thermal management systems. As most of the previous correlations for predicting the Leidenfrost temperature and the rewet temperature have been basically developed for water, they are shown to disagree with recent liquid nitrogen pipe chilldown experiments in upward and downward flow directions over a wide range of flow rates, pressures, and degrees of inlet subcooling. In addition to a complete review of the literature, two new correlations are presented in this work, one based on bubble growth and another based on the theoretical maximum limit of superheat. Each correlation performs well over the entire data set.

liquid nitrogen↗
DOE OSTI2025

Integrated Molten Salt Reactor Modeling Capabilities in NEAMS Thermal Hydraulics Tools

The DOE neams program supports a full range of computational thermal fluids analysis capabilities and code developments for a broad range of advanced reactor concepts. The research and development approach under the thermal fluids technical area synergistically combines three length and time scales in a hierarchical multi-scale approach. To enable multi-scale thermal fluids capability using these codes, a key joint effort has been underway to develop an integrated system- and engineering-scale thermal fluids analysis capability, through integration of SAM and Pronghorn codes, both based on the MOOSE framework. This report summarizes recent advances in developing an integrated system- and engineering-scale modeling capability for the msr concept, which has gained significant interest in recent years. A consistent framework was established by coupling Pronghorn and SAM through the Saline interface, with thermophysical properties provided by the Molten Salt Thermal Property Database (MSTDB-TP). Further improvements were made to the coupling schemes and domain-overlapping strategies, enhancing the stability and robustness of multi-code simulations. Verification and validation efforts demonstrate the accuracy of this integration across a range of benchmark problems, including one-dimensional heated pipe flows, three-dimensional natural convection loops with evolving isotopic compositions, and \gls{msre} demonstration cases. Within Pronghorn, new capabilities were introduced to model corrosion and noble-metal plating phenomena, supported by an extended thermal-hydraulics framework and refined turbulence treatments. To capture two-phase flow behavior, a multiphase Euler–Euler model was implemented in Pronghorn, including advanced closure relations, high-resolution advection techniques, and capillary force reconstruction. Preliminary verification cases confirm the fidelity of the approach, while planned validation efforts target canonical multiphase benchmarks and application to msr components such as the msre pump bowl. Finally, updates to SAM’s msr mass transfer modeling were extended to consider noble gas migration into porous structures like graphite. The point kinetics model was updated to include reactivity feedback contributions from any defined species, such as xenon. The gas transport model was expanded for applicability to gas mixtures, bubble efflux phenomena, and species transport between liquid and gas phases. A selection of multi-scale Sherwood number correlations from MOSCATO/NekRS and multi-phase correlations from literature have been added for improved accuracy in calculating mass transfer coefficients. A companion effort on developing system-level redox corrosion has also been incorporated into SAM. Collectively, these enhancements strengthen the predictive capability of SAM and Pronghorn for simulating MSR thermal-hydraulics, corrosion, multiphase behavior, and fission-product transport, providing a more complete toolset for design, safety analysis, and licensing support of next-generation \gls{msr}s.

42 - ENGINEERING↗
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)↗
NASA NTRS1975

Non-equilibrium one-dimensional two-phase flow in variable area channels

A one-dimensional nonequilibrium flow analysis has been formulated for a one component two phase flow. The flow is considered homogeneous and essentially isothermal. Phase change is assumed to occur at heterogeneous nucleation sites and the growth of the vapor bubbles is governed by heat conduction from the liquid to the bubble. The analysis adjusted for friction is applied to liquid nitrogen flow in a venturi and comparison is made with the NASA experimental results of Simoneau. Good agreement with the experiments is obtained when one assumes the effective activation energy for nucleus formation to be small but nonzero. The computed pressure distributions deviate from the experimental results in the throat region of the venturi in a manner consistent with centrifugal effects not accounted for in the one-dimensional theory. The results are shown to depend not only on cavitation number but on additional dimensionless parameters governing the nonequilibrium production and subsequent growth of nuclei.

Rohatgi, U. S.↗
NASA NTRS1980

Conceptual design of two-phase fluid mechanics and heat transfer facility for spacelab

Five specific experiments were analyzed to provide definition of experiments designed to evaluate two phase fluid behavior in low gravity. The conceptual design represents a fluid mechanics and heat transfer facility for a double rack in Spacelab. The five experiments are two phase flow patterns and pressure drop, flow boiling, liquid reorientation, and interface bubble dynamics. Hardware was sized, instrumentation and data recording requirements defined, and the five experiments were installed as an integrated experimental package. Applicable available hardware was selected in the experiment design and total experiment program costs were defined.

North, B. F.↗