Search NASA⌕ Search

SEARCH · Search NASA

Results for “Two-fluid 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.

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 NTRS1976

A two-region model of the solar wind including azimuthal velocity

The two-region model of the solar wind divides the interplanetary space into two regions: it assumes that the solar wind is one-fluid in an inner region within 0.4 AU and two-fluid in an outer region beyond 0.4 AU. This paper includes the angular motion of the solar wind in the two-region model. The flow in the one-fluid region is governed by the one-fluid magnetohydrodynamic equations. The second and third moment equations of the Vlasov equation together with other conservation equations are used to describe the solar-wind flow in the two-fluid region. The predicted azimuthal velocity at 1 AU is less than 2 km/s. All other macroscopic and microscopic properties from this model are in good agreement with experimental quiet-time observations at 1 AU. The numerical results also confirm that when the azimuthal velocity is included in the analysis, the amount of magnetic-field energy converted into kinetic energy in the solar wind is only a small fraction of the total expansion energy flux and has little effect upon the final radial expansion velocity.

Acuna, M. H.↗
DOE OSTI2025

A Finite Element Method for Compressible and Turbulent Multiphase Flow Instabilities with Heat Transfer

We present a new finite element framework for modeling compressible, turbulent multiphase flows with heat transfer. For two-fluid systems with a free surface, the Volume of Fluid (VOF) method is implemented without the need for interface reconstruction, while turbulence is resolved using a dynamic Vreman large eddy simulation (LES) model. Unlike most two-phase VOF studies, which neglect heat transfer, the present approach incorporates energy transport equations within the VOF formulation to account for heat exchange, an effect particularly important in turbulent flows. Conjugate heat transfer is often challenging in finite volume methods, which require explicit specification of heat fluxes at the solid–fluid interface, limiting accuracy and predictive capability. By contrast, the finite element formulation does not require heat flux inputs, allowing more accurate and robust simulation of heat transfer between solids and fluids. The method is demonstrated through three representative cases. First, a two-fluid instability with a single-mode perturbation is simulated and validated against analytical growth rates. Second, conjugate heat transfer is examined in a high-temperature flow over a cold metal cylinder, with validation performed both quantitatively—via pressure coefficient comparisons with experimental data—and qualitatively using vector field topology. Finally, compressible spray injection and breakup are modeled, demonstrating the ability of the framework to capture interfacial dynamics and atomization under turbulent, high-speed conditions. In the compressible spray injection and breakup case, the results indicate that the finite element formulation achieved higher predictive accuracy and robustness than the finite-volume method. With the same mesh resolution, the FEM reduced the root mean square error (RMSE) and mean absolute percentage error (MAPE) from 6.96 mm and 26.0% (for the FVM) to 4.85 mm and 12.7%, respectively, demonstrating improved accuracy and robustness in capturing interfacial dynamics and heat transfer. The study also introduced vector field topology to visualize and interpret coherent flow structures and instabilities, offering insights beyond conventional scalar-field analyses.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Effects beyond ideal MHD on stability of wide and enhanced pedestal regimes in NSTX

Stability of edge-localized modes (ELMs) in spherical tokamaks is explored using the extended MHD model. Linear NIMROD simulations have been performed for three NSTX discharges 132543, 132588, and 141133, to investigate the role of resistivity, diffusivity, and shear flows on the onset of ELMs. The first discharge represents the wide pedestal regime and the later two discharges represent ELM-free enhanced pedestal H-mode. We first present the effect of toroidal rotation shear and find a flow shear destabilizing effect in these NSTX discharges. Simulations are also extended to include the two-fluid and ion gyroviscosity effects. Simulations show that the flow shear can shift the mode spectrum and alter the critical condition of ELM onset. We also uncover that ELM onset prediction in spherical tokamaks requires effects beyond MHD, in particular gyroviscosity and diamagnetic terms could stabilize Peeling-Ballooning modes consistent with the experimental observation of ELM-free regimes in NSTX. The findings give new insight into the nature of the interplay between resistivity, flows, and diamagnetic stabilization in ELM suppression and have potential applications to ELM control schemes in NSTX-U and next-generation spherical tokamaks. This study identifies the essential physical effects that must be included in future predictive and validation simulations.

NSTX↗
DOE OSTI2026

Modeling of injected current stream-induced 3D perturbations in local helicity injection plasmas

Solenoid-free tokamak startup techniques are essential for spherical tokamaks and offer a pathway to cost reduction and design simplification in fusion energy systems. Local helicity injection (LHI) is one such approach, employing compact edge current sources to drive open field line current that initiates and sustains tokamak plasmas. The recently commissioned Pegasus-III Experiment provides a platform for advancing this and other solenoid-free startup methods. This study investigates the effect of LHI on magnetic topology in Pegasus-III plasmas. A helical filament model represents the injected current, and the linear plasma response to its three-dimensional field is calculated with M3D-C1. Poincaré mapping reveals substantial flux surface degradation in all modeled cases. The onset of overlapping magnetic structures and large-scale surface deformation begins near Ψ N ≈ 0.37, indicating a broad region of perturbed topology extending toward the edge. In rotating plasmas, both single-fluid and two-fluid models exhibit partial screening of the n = 1 perturbation, with two-fluid calculations showing stronger suppression near the edge. In contrast, the absence of rotation leads to strong resonant field amplification in the single-fluid case, while the two-fluid case with zero electron rotation mitigates this amplification and preserves edge screening. Magnetic probe measurements indicate that modeling the stream with spatial spreading—representing distributed current and/or oscillatory motion—better reproduces measured magnetic power profiles than a rigid filament model. The results underscore the role of rotation and two-fluid physics in screening stream perturbations and point to plasma flow measurements and refined stream models as key steps toward improving predictive fidelity.

Schaefer, Carolyn E. [Univ. of Wisconsin, Madison,↗
DOE OSTI2026

A filter-dependent granular temperature model from large-scale CFD-DEM data

The computational study of strongly-coupled, gas–solid flows at scales relevant to most environmental and engineering applications requires the use of ‘coarse-grained’ methodologies such as the two-fluid model, particle-in-cell approach or the multiphase Reynolds Averaged Navier–Stokes equations. While these strategies enable computations at desirable length- and time-scales, they rely heavily on models to capture important flow physics that occur at scales smaller than the mesh. To date, the models that do exist are based on a limited set of flow conditions, such as very dilute particle phase. To this end, we leverage a large-scale repository of CFD-DEM data to develop filter-size dependent models for the mean variance in particle volume fraction, a quantity commonly used to assess the degree of clustering, and the granular temperature, a key quantity for accurately predicting gas–solid flows. In conclusion, because of its filter-size dependence, the granular temperature model can be directly translated to coarse-grained approaches and tied directly to grid size.

AMReX↗
NASA NTRS1968

Two-fluid, impinging-sheet injector

Two-fluid, impinging-sheet propellant injector reduces the severe erosion found to occur when ejector elements are directly exposed during throttling without the benefits of a cooling flow of the propellant liquids. It greatly improves combustion efficiency by venting the secondary stream of combustion gases generated by backspray reaction.

Riebling, R. W.↗
NASA NTRS1972

A solar-wind model including proton thermal anisotropy.

A model of the solar wind is presented which divides the interplanetary space into two regions. It assumes that the solar wind is one-fluid in an 'inner region' of radius less than 0.4 A, and two-fluid and collisionless in an 'outer region' of radius greater than 0.4 A. The second and third moment equations of the Vlasov equation together with conservation equations are used to govern the solar wind flow in the outer region. The model produces solutions for all macroscopic quantities as well as for the microscopic proton distribution function, both as functions of heliocentric distances up to 50 A. Calculated results at 1 A include the solar wind velocity, electron temperature, proton temperature, proton thermal anisotropy and magnetic field magnitude. They are in good agreement with observations.

Whang, Y. C.↗
DOE OSTI2025

Peeling-ballooning modes in spherical tokamaks: Multi-branch instabilities and effects beyond ideal MHD

A number of important physics effects on the stability of relatively high-n (n is the toroidal mode number) peeling-ballooning modes (PBMs) are investigated based on an equilibrium reconstructed from a NSTX discharge, utilizing extended magnetohydrodynamic (MHD) eigenvalue solvers. For a given toroidal mode number n, multiple branches of instabilities are computed, with the total number of unstable branches roughly linearly scaling with n. Most of the unstable branches are located in the plasma core region, but edge-localized branches, i.e., PBMs, are also identified at higher n-numbers. For the single-fluid-wise most unstable PBM with n = 19⁠, stabilizing/destabilizing effects due to various physics beyond ideal MHD are systematically investigated. Plasma toroidal flow is found to be weakly stabilizing. Local flow shear is generally stabilizing as well, with the degree of stabilization depending on the initial growth rate (without flow shear) of the mode. The plasma resistivity can strongly destabilize the PBM within the single-fluid framework. Anisotropic thermal transport, strong parallel sound wave damping, as well as two-fluid effects are all stabilizing to the mode. In particular, diamagnetic stabilization (within the two-fluid model) is found to be very strong for this mode.

Linear stability analysis↗
DOE OSTI2025

Modeling and Simulation of Fuel Dispersal During the Loss-of-Coolant Accident

This document is the compilation of the milestone portion to a larger end of project NEUP report. The executive summary of the modeling portion is provided below: In the event of cladding rupture during a postulated LOCA in a pressurized water reactor, fuel particles, along with fission gases, can be expelled into the reactor core from the fractured fuel rod, a phenomenon referred to as fuel dispersal. The initial stage of fuel dispersal is strongly influenced by the high-pressure ejection of fuel fragments, the size and geometry of the ruptured cladding, and the depressurization history of the fuel rod during the postulated LOCA transient. Depending on the location of the burst orifice relative to the quench front, the dispersal event represents an intricate three-phase flow and heat transfer phenomenon, where high-temperature fuel particles carried by the fission gases interact with the coolant within the narrow subchannels of the fuel assemblies, inducing localized phase change. Given the unique multiphysics nature of this phenomena, the current study develops a dedicated computational framework to predict the mass distribution and cooling of dispersing fuel particles, facilitating post-accident assessment and management of the fuel assemblies. Considering the scale of nuclear reactor applications, a continuum three-fluid model is proposed for simulating the transport of solids within the reactor core. With high-temperature fuel fragments within the liquid media, nucleation sites inducing phase changes are dispersed within the flow domain. Coupled with the fact that the transient dispersal event occurs on different time scales than other three-phase flow applications, this study derives a time-averaged three-fluid flow model without losing generality. The assumptions regarding the continuum treatment of the solid phase and the modeling of fuel dispersal behavior are incorporated to simplify the governing equations and derive applicable closure relations. The computational validation of the model was conducted using adiabatic experimental results obtained from ongoing research at Oregon State University, focusing on characterizing fuel dispersal behavior during simulated LOCA conditions. Settlement characteristics of the solids, quantified by the probability distribution of equivalent particles, closely matched the probability density functions reported in experimental studies. The transport of fuel particles within a scaled 5 × 5 lattice of a pressurized-water reactor rod bundle geometry was modeled through a two-fluid Eulerian framework. The required boundary conditions were evaluated from the fuel performance code BISON in a postulated large-break LOCA scenario. The modeling framework considered solid fuel particles as granular matter, interacting with the gaseous dry steam phase and fission gases through the governing interfacial momentum exchange between the participating fluids. The simulation results provided the volume fraction of the solids obtained at the bottom surface of the enclosing tank geometry. Postulated LOCA leading to fuel dispersal phenomena involves the strong coupling between fuel thermomechanics, cladding deformation, thermal-hydraulics, and fuel particle transport. Incorporation of such a strong coupling in numerical simulation is performed by coupling the multiphysics solvers. In the case of fuel dispersal, a strong coupled simulation can be performed by coupling the BISON code for fuel performance, the TRACE code for system-level thermal hydraulics, and fuel particle transport in Multiphysics Object-Oriented Simulation Environment (MOOSE). For such intricate infrastructure, the MOOSE Framework eases the data transfer between codes. The recent version of MOOSE has incorporated the Navier-Stokes module for the fluid flow. An exploratory exercise was done to gain familiarity with finite volume capabilities in the MOOSE framework to incorporate the Spalart-Allmaras (SA) turbulence model. New finite-volume and auxiliary kernels were introduced to assemble the SA transport equation, compute turbulent viscosity, and evaluate wall distance and diagnostic turbulence terms, fully integrated with existing Navier-Stokes modules. A turbulent lid-driven cavity at a Reynolds number of approximately 10,000 is used for verification. MOOSE shows the robust solver convergence and produces the turbulent features. But it underpredicts the velocity profile and turbulent quantities, emphasizing the need to develop improved SA near-wall treatments (e.g., low-Re corrections or wall functions) as a key direction for future work.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Combined Influence of Rotation and Scrape-Off Layer Drifts on Recycling Asymmetries in Tokamak Plasmas

Coupled 2D fluid-kinetic simulations of a DIII-D high confinement tokamak plasma show that plasma rotation coupled with drift effects near the plasma edge play a significant role in the creation of the observed poloidal distribution of neutrals. It is observed that including either drift or rotation effects enhances particle flux at the inner target in the case of ion 𝐵×∇𝐵 drift toward the 𝑋-point. However, the particle flux asymmetry is significantly higher with the combination of drifts and rotation than either effect alone. The heightened particle flux asymmetry allows for improved simulation of the strong in-out asymmetry of the Lyman-𝛼 brightness profiles measured in the experiment. Enhancement of radial transport of parallel momentum changes the upstream scrape-off layer flow pattern, increasing the fraction of deuterium flux that reaches the inboard divertor entrance while lowering that which arrives at the outboard. In conclusion, this Letter indicates that by combining drifts, rotation, and viscous coupling, existing boundary plasma models can achieve a satisfactory agreement with experimentally measured neutral asymmetries.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Filament extension atomizers: A study of integration potential into spray dryers for high-viscosity and non-Newtonian fluids

Filament extension atomizers (FEAs) are an emerging class of spray nozzles designed to atomize high-viscosity and non-Newtonian fluids that are challenging for conventional pressure or two-fluid nozzles. In this study, we investigate the influence of roller geometry, surface velocity, roller material, and fluid rheology on the atomization performance of FEA systems using concentrated whey protein suspensions (50–70 wt. %). Extensional and shear rheology experiments, along with high-speed imaging and particle image velocimetry, reveal that filament breakup dynamics are governed by competition between inertial, capillary, and viscoelastic stresses. High roller rotational velocity leads to narrower spray cones, contradicting rheology experiments and suggests a significant inertial contribution to filament breakup. Smaller rollers operating at the same rotational velocity led to broader spray cones consistent with expectations. An FEA nozzle was integrated into a conventional dryer producing particles in the 100 μm range. Results suggest that FEA technology enables atomization of highly viscous fluids at industrially relevant spray cone angles similar to those generated by pressure nozzles, offering a pathway to improve energy efficiency in spray drying by enabling higher solids loading feedstocks. Furthermore, these insights provide critical guidance for optimizing FEA nozzle designs and process parameters across a range of applications.

Energy efficiency↗