SEARCH · Search NASA

Results for “Interfacial flows”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 217 records · Page 12

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Development and Validation of Two-Phase Flow Models in MOOSE and Application to Molten Salt Reactors

Two-phase flow in Molten Salt Reactors (MSRs) is important as it impacts reactivity evolution, reactor transient response, and the removal of species dissolved in the molten salt through gas phase transfer. Therefore, accurately predicting the gas distribution and the associated liquid-gas interface area in MSRs is essential for their design and operation. Recently, we integrated a new two-phase model into Idaho National Laboratory (INL)’s Multiphysics Object-Oriented Simulation Environment (MOOSE): a multi-D generalization of a mixture drift-flux model. It provides greater computational efficiency, which is typically preferred for modeling reactor transients. However, the mixture model’s accuracy in capturing void distribution and interfacial area in MSRs still needs to be assessed. This article begins with a description of the mathematical framework for the two-phase model implemented in MOOSE. It then presents validation of these models against relevant experimental data. Finally, the model is applied to the Molten Salt Reactor Experiment case study, analyzing various operational conditions such as different rates of fission product volatilization and diverse cover gas entrainment scenarios at the reactor pump. The article concludes by assessing the suitability of the mixture drift-flux model for capturing the two-phase flow dynamics critical to MSR operations

42 - ENGINEERING↗

NASA NTRS1999

Numerical Modeling of Three-Dimensional Fluid Flow with Phase Change

We present a numerical method to compute phase change dynamics of three-dimensional deformable bubbles. The full Navier-Stokes and energy equations are solved for both phases by a front tracking/finite difference technique. The fluid boundary is explicitly tracked by discrete points that are connected by triangular elements to form a front that is used to keep the stratification of material properties sharp and to calculate the interfacial source terms. Two simulations are presented to show robustness of the method in handling complex phase boundaries. In the first case, growth of a vapor bubble in zero gravity is studied where large volume increase of the bubble is managed by adaptively increasing the front resolution. In the second case, growth of a bubble under high gravity is studied where indentation at the rear of the bubble results in a region of large curvature which challenges the front tracking in three dimensions.

Esmaeeli, Asghar↗

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 NTRS1985

Surface-Tension Driven Convection

The present experiments are directed to examine surface-tension gradients. The surface-tension induced convection under reduced gravity it's comparison with normal gravity convection with all other identical conditions are studied. The intensities of such flows and their penetration depth below the free surface are also determined. In this way, the dominant role of a free surface in a reduced-gravity environment is indicated. The interface between two fluid phases can influence the motion of fluids when either the interface has finite curvature or when the interfacial tension varies from point to point. In both cases, forces appear in the interfacial region that can affect or generate fluid motions.

Ostrach, S.↗

NASA NTRS2009

Barium Depletion in Hollow Cathode Emitters

The effect of tungsten erosion, transport and redeposition on the operation of dispenser hollow cathodes was investigated in detailed examinations of the discharge cathode inserts from an 8200 hour and a 30,352 hour ion engine wear test. Erosion and subsequent re-deposition of tungsten in the electron emission zone at the downstream end of the insert reduces the porosity of the tungsten matrix, preventing the ow of barium from the interior. This inhibits the interfacial reactions of the barium-calcium-aluminate impregnant with the tungsten in the pores. A numerical model of barium transport in the internal xenon discharge plasma shows that the barium required to reduce the work function in the emission zone can be supplied from upstream through the gas phase. Barium that flows out of the pores of the tungsten insert is rapidly ionized in the xenon discharge and pushed back to the emitter surface by the electric field and drag from the xenon ion flow. This barium ion flux is sufficient to maintain a barium surface coverage at the downstream end greater than 0.6, even if local barium production at that point is inhibited by tungsten deposits. The model also shows that the neutral barium pressure exceeds the equilibrium vapor pressure of the impregnant decomposition reaction over much of the insert length, so the reactions are suppressed. Only a small region upstream of the zone blocked by tungsten deposits is active and supplies the required barium. These results indicate that hollow cathode failure models based on barium depletion rates in vacuum dispenser cathodes are very conservative.

plasma discharges↗

NASA NTRS2015

Modeling of Non-Isothermal Cryogenic Fluid Sloshing

A computational fluid dynamic model was used to simulate the thermal destratification in an upright self-pressurized cryostat approximately half-filled with liquid nitrogen and subjected to forced sinusoidal lateral shaking. A full three-dimensional computational grid was used to model the tank dynamics, fluid flow and thermodynamics using the ANSYS Fluent code. A non-inertial grid was used which required the addition of momentum and energy source terms to account for the inertial forces, energy transfer and wall reaction forces produced by the shaken tank. The kinetics-based Schrage mass transfer model provided the interfacial mass transfer due to evaporation and condensation at the sloshing interface. The dynamic behavior of the sloshing interface, its amplitude and transition to different wave modes, provided insight into the fluid process at the interface. The tank pressure evolution and temperature profiles compared relatively well with the shaken cryostat experimental test data provided by the Centre National D'Etudes Spatiales.

Cryogenics↗

DOE OSTI2025

Development and Validation of Two-Phase Flow Models in MOOSE for Molten Salt Reactor Application

Two-phase flow in Molten Salt Reactors (MSRs) is important as it impacts reactivity evolution, reactor transient response, and the removal of species dissolved in the molten salt through gas phase transfer. Therefore, accurately predicting the gas distribution and the associated liquid-gas interface area in MSRs is essential for their design and operation. Recently, we integrated two new models into Idaho National Laboratory (INL)’s Multiphysics Object-Oriented Simulation Environment (MOOSE): a multi-D generalization of a mixture drift-flux model and a Euler-Euler model. The Euler-Euler model offers higher fidelity, while the mixture drift-flux model provides greater computational efficiency, which is typically preferred for modeling reactor transients. However, the mixture model's accuracy in capturing void distribution and interfacial area in MSRs still needs to be assessed. This article begins with a description of the mathematical framework for the two-phase models implemented in MOOSE. It then presents validation of these models against relevant experimental data. Finally, both models are applied to the Molten Salt Reactor Experiment case study, analyzing various operational conditions such as different rates of fission product volatilization and diverse cover gas entrainment scenarios at the reactor pump. The article concludes by assessing the suitability of both models for capturing the two-phase flow dynamics critical to MSR operations.

42 - ENGINEERING↗

DOE OSTI2025

Development and Validation of Two-Phase Flow Models in MOOSE and Application to Molten Salt Reactors

Two-phase flow in Molten Salt Reactors (MSRs) is important as it impacts reactivity evolution, reactor transient response, and the removal of species dissolved in the molten salt through gas phase transfer. Therefore, accurately predicting the gas distribution and the associated liquid-gas interface area in MSRs is essential for their design and operation. Recently, we integrated two new models into Idaho National Laboratory (INL)’s Multiphysics Object-Oriented Simulation Environment (MOOSE): a multi-D generalization of a mixture drift-flux model and a Euler-Euler model. The Euler-Euler model offers higher fidelity, while the mixture drift-flux model provides greater computational efficiency, which is typically preferred for modeling reactor transients. However, the mixture model's accuracy in capturing void distribution and interfacial area in MSRs still needs to be assessed. This article begins with a description of the mathematical framework for the two-phase models implemented in MOOSE. It then presents validation of these models against relevant experimental data. Finally, both models are applied to the Molten Salt Reactor Experiment case study, analyzing various operational conditions such as different rates of fission product volatilization and diverse cover gas entrainment scenarios at the reactor pump. The article concludes by assessing the suitability of both models for capturing the two-phase flow dynamics critical to MSR operations.

42 - ENGINEERING↗

NASA NTRS1982

An operational large-scale marine planetary boundary layer model

A marine planetary boundary layer (PBL) model is presented and compared with data from sea-based experiments. The PBL model comprises two layers, the outer an Ekman-Taylor layer with stratification-dependent secondary flow, and the logarithmic surface layer corrected for stratification and humidity effects and variable surface roughness. Corrections are noted for air much warmer than water in stable conditions and for low wind speeds. The layers are analytically defined along with similarity relations and a resistance law for inclusion in a program. An additional interfacial layer correction is developed and shown to be significant for heat flux calculations. Experimental data from GOASEX were used to predict the windfield in the Gulf of Alaska, and JASIN data was used for windfields SE of Iceland. The JASIN-derived wind field predictions were accurate to within 1 m/sec and 10 deg in a 200 km triangle.

Brown, R. A.↗

NASA NTRS2023

Hydrodynamics of bubble flow through a porous medium with applications to packed bed reactors

Gas-liquid flows through packed bed reactors (PBRs) are challenging to predict due to the tortuous flow paths that fluid interfaces must traverse. Experiments at the International Space Station showed that bubble and pulse flows are predominately observed under microgravity conditions, while the trickle and spray flows observed under terrestrial conditions are not present in microgravity. To understand the physics behind the former experiments, we simulate bubble flow through a PBR for different packing-particle-diameter-based Weber numbers and under different gravity conditions. We demonstrate different pore-scale mechanisms, such as capillary entrapment, buoyancy entrapment, and inertia-induced bubble displacement. Then, we perform a quantitative analysis by introducing new dynamic scales, dependent upon the evolving gas-liquid interfacial area, to understand the dynamic trade-offs between the inertia, capillary, and buoyancy forces on a bubble passing through a PBR. This analysis leads us to define new dimensionless Weber-like numbers that delineate bubble entrapment from bubble displacement.

Pranay P Nagrani↗

DOE OSTI2024

An integrated approach to derive relative permeability from capillary pressure

Surface tension affects all aspects of fluid flow in porous media. Through measurements of surface tension interaction under multiphase conditions, a relative permeability curve can be determined. Relative permeability is a numerical description of the interaction between two or more fluids and the porous media. It is a critical parameter for various tools that characterize subsurface multiphase flow systems, such as numerical simulation for carbon sequestration, oil and gas development, and groundwater contamination remediation. Therefore, it is critical to get a good statistical distribution of relative permeability in the porous media under study. Empirical formula for determining relative permeability from capillary pressure are already well established but do not provide the needed flexibility that is required to match laboratory-derived relative permeability curves. By expanding the existing methods for calculating relative permeability from capillary pressure data, it is possible to create both two and three-phase relative permeability curves. Mercury intrusion capillary pressure (MICP) data from the Morrow 'B' Sandstone coupled with interfacial tension and contact angle measurements were used to create a suite of relative permeability curves. Furthermore, these curves were then calibrated to a small sample of existing laboratory curves to elucidate common fitting parameters for the formation that were then used to create relative permeability curves from MICP data that does not have an associated laboratory-measured relative permeability curve.

58 GEOSCIENCES↗

NASA NTRS2003

Transient Interfacial Phenomena in Miscible Polymer Systems (TIPMPS)

Almost one hundred years ago Korteweg published a theory of how stresses could be induced in miscible fluids by concentration gradients, causing phenomena that would appear to be the same as with immiscible fluids. Miscible fluids could manifest a transient or effective interfacial tension (EIT). To this day, there has been no definitive experiment to confirm Korteweg's model but numerous fascinating and suggestive experiments have been reported. The goal of TIPMPS is to answer the question: Can concentration and temperature gradients in miscible materials induce stresses that cause convection? Many polymer processes involving miscible monomer and polymer systems could be affected by fluid flow and so this work could help understand miscible polymer processing, not only in microgravity, but also on earth. Demonstrating the existence of this phenomenon in miscible fluids will open up a new area of study for materials science. The science objectives of TIPMPS are: (1) Determine if convection can be induced by variation of the width of a miscible interface; (2) Determine if convection can be induced by variation of temperature along a miscible interface; (3) Determine if convection can be induced by variation of conversion along a miscible interface An interface between two miscible fluids can best be created via a spatially-selective photopolymerization of dodecyl acrylate with a photoinitiator, which allows the creation of precise and accurate concentration gradients between polymer and monomer. Optical techniques will be used to measure the refractive index variation caused by the resultant temperature and concentration fields. The viscosity of the polymer will be measured from the increase in the fluorescence of pyrene. Because the large concentration and temperature gradients cause buoyancy-driven convection that prevents the observation of the predicted flows, the experiment must be done in microgravity. In this report, we will consider our efforts to estimate the square gradient parameter, k, and our use of the estimates in modeling of the planned TIPMPS experiments. We developed a model consisting of the heat and diffusion equations with convective terms and of the Navier-Stokes equations with an additional volume force written in the form of the Korteweg stresses arising from nonlocal interaction in the fluid. The fluid's viscosity dependence on polymer conversion and temperature was taken from measurements of poly(dodecyl acrylate). Numerical modeling demonstrated that significant flows would arise for conditions corresponding to the planned experiments.

Pojman, John A.↗

NASA NTRS1980

Improved refractory coatings

The adhesion, friction and wear properties of sputtered refractory coatings on substrates of materials that form stable nitrides are enhanced by placing each substrate directly below a titanium carbide target of a commercial radiofrequency diode apparatus in a vacuum chamber. Nitrogen is bled into the system through a nozzle resulting in a small partial pressure of about 0.5% to 2.5% during the first two minutes of deposition. The flow of nitrogen is then stopped, and the sputtering ambient is reduced to pure argon through a nozzle without interrupting the sputtering process. When nitrogen is deliberately introduced during the crucial interface formation, some of the titanium at the interface reacts to form titanium nitride while the metal of the substrate also forms the nitride. These two nitrides atomically mixed together in the interfacial region act to more strongly bond the growing titanium carbide coating as it forms on the substrate.

Brainard, W. A.↗

NASA NTRS1982

Refractory coatings and method of producing the same

The adhesion, friction, and wear properties of sputtered refractory coatings on substrates of materials that form stable nitrides is improved by placing each substrate directly below a titanium carbide target of a commercial radiofrequency diode apparatus in a vacuum chamber. Nitrogen is bled into the system through a nozzle resulting in a small partial pressure of about 0.5% to 2.5% during the first two minutes of deposition. The flow of nitrogen is then stopped, and the sputtering ambient is reduced to pure argon through a nozzle without interrupting the sputtering process. When nitrogen is deliberately introduced during the crucial interface formation, some of the titanium at the interface reacts to form titanium nitride while the metal of the substrate also forms the nitride. These two nitrides atomically mixed together in the interfacial region act to more strongly bond the growing titanium carbide coating as it forms on the substrate.

Brainard, W. A.↗

NASA NTRS1992

Adaptive multigrid domain decomposition solutions for viscous interacting flows

Several viscous incompressible flows with strong pressure interaction and/or axial flow reversal are considered with an adaptive multigrid domain decomposition procedure. Specific examples include the triple deck structure surrounding the trailing edge of a flat plate, the flow recirculation in a trough geometry, and the flow in a rearward facing step channel. For the latter case, there are multiple recirculation zones, of different character, for laminar and turbulent flow conditions. A pressure-based form of flux-vector splitting is applied to the Navier-Stokes equations, which are represented by an implicit lowest-order reduced Navier-Stokes (RNS) system and a purely diffusive, higher-order, deferred-corrector. A trapezoidal or box-like form of discretization insures that all mass conservation properties are satisfied at interfacial and outflow boundaries, even for this primitive-variable, non-staggered grid computation.

Rubin, Stanley G.↗

DOE OSTI2024

Effects of particle size and AQDS on the flow of electron equivalents between magnetite and aqueous Fe2+

Magnetite can occur naturally in nano- to micro-size regimes and widely coexists with aqueous Fe2+ (Fe2+ (aq)) in natural environments. However, the effects of magnetite particle size on its interaction with Fe2+ (aq) in anoxic subsurface environments, particularly with redox-active organics, remain unclear. In this study, the interactions of Fe2+ (aq) with magnetite particles of 12 nm versus 109 nm (Mag-12 vs. Mag-109), with/without anthraquinone- 2,6-disulfonate (AQDS), were studied based on equilibrium Fe2+ (aq) concentrations, kinetics of AQDS reduction, and structural versus surface-localized Fe(II)/Fe(III) ratios (xstru and xsurf) of magnetite. In the absence of AQDS, Mag-12 tends to release Fe2+ (aq) at pH 7 but sorb Fe2+ (aq) at pH 8, while Fe2+ (aq) uptake by Mag-109 is observed at both pH 7 and 8. The amounts of Fe2+ (aq) adsorbed per unit area of Mag-109 is higher than that of Mag-12, due to the higher electron-accepting capacity of Mag-109 that facilitates interfacial electron transfer (IET) from surfaceassociated Fe(II) to structural Fe(III). The increases of xstru and xsurf in Mag-109 after reaction with Fe2+ (aq) at pH 7 and 8 suggest Fe2+ (aq) incorporation or electron injection into the structure of Mag-109. The presence of AQDS promotes Fe2+ (aq) uptake by both Mag-12 and Mag-109. However, AQDS reduction by Fe2+-amended Mag-12 results in the decrease of xstru and inhibits Fe2+ (aq) incorporation or electron injection into the structure. On the contrary, the increase of xstru observed in Fe2+-amended Mag-109 after reaction with AQDS suggests that Fe2+ (aq) incorporation or electron injection into the surface structure and then consequently into the interiors is more favorable for magnetite with larger particle sizes. The different flow directions of electron equivalents across the solid-solution interfaces can be attributed to the relatively higher electron-accepting capacity, i.e. redox potential, of Mag-109 than Mag-12; larger particle sizes facilitate IET from surface-associated Fe(II) to structural Fe(III) and promotes further Fe2+ (aq) uptake, culminating in the pronounced changes of redox potentials in magnetitebearing solutions. The results demonstrate that particle size and redox-active organics are important factors to affect reductive activity of Fe2+-magnetite system in redox-oscillating environments.

Peng, Huan↗

DOE OSTI2024

Evaluating Pressure‐dependent Discharge Behavior of Foil Versus In situ Plated Lithium Metal Anodes in Solid‐State Batteries

Abstract Anode‐free manufacturing of solid‐state batteries (SSBs) shows promise to maximize energy density by eliminating excess lithium (Li) and simplifying battery production. However, high reversibility during discharge (stripping of Li) is necessary for long‐lifetime SSBs with a limited Li reservoir. Further, the plastic flow of Li changes depending on the Li thickness, leading to possible differences in discharge performance under stack pressure. This work investigates the pressure‐dependent discharge performance of anode‐free manufactured SSBs with in situ plated Li and compares the performance to that of conventional thick Li foil cells. Distinct stripping behavior is observed at low pressures (0–1 MPa), where Li diffusivity and initial interfacial contact may control accessible capacity, compared to high pressures (3–10 MPa) where mechanical deformation of Li likely governs stripping behavior. Analysis of impedance spectra collected during stripping shows that additional stack pressure delays the formation of deep, as opposed to lateral, voids in the Li anode. These results provide insights to guide the transition from thick Li foil anodes to anode‐free manufactured SSBs.

Haslam, Catherine G.↗

NASA NTRS1990

Kinetics and transport at AMTEC electrodes. I - The interfacial impedance model

Mixed mass-transport and kinetic control of sodium ion reduction at porous inert electrodes on sodium beta-double-prime alumina solid electrolyte (BASE) ceramic in a high-temperature electrochemical cell has been observed and modeled. The high ionic conductivity of BASE and the reversibility of the liquid sodium/BASE anodic half-cell led to assignment of potential-dependent (nonohmic) resistances to kinetic and mass-transport processes associated with the porous electrode. The morphology of these electrodes and typical sodium gas pressures are consistent with Knudsen, or free-molecular, flow through the electrode.

Williams, R. M.↗