SEARCH · Search NASA

Results for “TWO-FLUID MODEL”

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 109 records · Page 6

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Linear simulations of a wide pedestal quiescent H-mode plasma with the extended-MHD code NIMROD

We present linear simulations of a Wide Pedestal QH (WPQH)-mode DIII-D plasma with the NIMROD code to address the role of ion multispecies collisionality on WPQH-modes. We use a full extended-MHD model that captures Alfvénic and magnetoacoustic waves, including two-fluid, ion gyroviscosity and cross heat flux. We show that two-fluid and ion gyroviscous effects destabilize electron-directed peeling–ballooning modes in the pedestal. The linear growth rates of the modes decrease when carbon (C) is added to a deuterium plasma because the multispecies collisional effects increase the plasma resistivity. When replacing C with tungsten (W), the multispecies collisional effects increase the resistivity further, as well as the stabilizing effect on the instabilities. The resistivity impact on mode stability is reversed when both two-fluid and ion gyroviscous effects are ignored, which is consistent with previous works, proving that multispecies collisional effects, together with two-fluid and ion gyroviscosity, should be considered to address the stability of plasmas with significant impurity content.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1982

Solar transition region response to variations in the heating rate

The response of a numerical model for the upper chromosphere, transition region, and corona to variations in the energy input has been examined. The numerical model solves the set of one-dimensional two-fluid hydrodynamic equations in a simple vertical magnetic flux tube. The atmosphere responds to both the increase and decrease in energy deposition by smoothly readjusting the temperature gradient and the amount of material in the region of peak radiating efficiency to radiate away energy being deposited. At no time during this readjustment is a departure from a thin laminar transition region structure seen. In addition, a time-dependent description of the nonequilibrium ionization of all of the ionization stages of oxygen has been included. This calculation is coupled with the self-consistent calculations of the dynamical variables. It is found that the nonequilibrium ionization balance calculations for both heating and cooling small loops in the quiet sun predict relative ionic abundances which differ substantially from those which would be predicted by an equilibrium calculation

Mariska, J. T.↗

NASA NTRS1990

A self-consistent two-dimensional resistive fluid theory of field-aligned potential structures including charge separation and magnetic and velocity shear

A self-consistent two-fluid theory that includes the magnetic field and shear patterns is developed to model stationary electrostatic structures with field-aligned potential drops. Shear flow is also included in the theory since this seems to be a prominent feature of the structures of interest. In addition, Ohmic dissipation, a Hall term, and pressure gradients in a generalized Ohm's law, modified for cases without quasi-neutrality, are included. In the analytic theory, the electrostatic force is balanced by field-aligned pressure gradients (i.e., thermal effects in the direction of the magnetic field) and by pressure gradients and magnetic stresses in the perpendicular direction. Within this theory, simple examples of applications are presented to demonstrate the kind of solutions resulting from the model. The results show how the effects of charge separation and shear in the magnetic field and the velocity can be combined to form self-consistent structures such as are found to exist above the aurora, suggested also in association with solar flares.

Hesse, Michael↗

NASA NTRS1989

Scaling and modeling of turbulent suspension flows

Scaling factors determining various aspects of particle-fluid interactions and the development of physical models to predict gas-solid turbulent suspension flow fields are discussed based on two-fluid, continua formulation. The modes of particle-fluid interactions are discussed based on the length and time scale ratio, which depends on the properties of the particles and the characteristics of the flow turbulence. For particle size smaller than or comparable with the Kolmogorov length scale and concentration low enough for neglecting direct particle-particle interaction, scaling rules can be established in various parameter ranges. The various particle-fluid interactions give rise to additional mechanisms which affect the fluid mechanics of the conveying gas phase. These extra mechanisms are incorporated into a turbulence modeling method based on the scaling rules. A multiple-scale two-phase turbulence model is developed, which gives reasonable predictions for dilute suspension flow. Much work still needs to be done to account for the poly-dispersed effects and the extension to dense suspension flows.

Chen, C. P.↗

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 OSTI2024

Interpreting Mass and Radius Measurements of Neutron Stars with Dark Matter Halos

Abstract The high densities of neutron stars (NSs) could provide astrophysical locations for dark matter (DM) to accumulate. Depending on the DM model, these DM admixed NSs (DANSs) could have significantly different properties than pure baryonic NSs, accessible through X-ray observations of rotation-powered pulsars. We adopt the two-fluid formalism in general relativity to numerically simulate stable configurations of DANSs, assuming a fermionic equation of state (EOS) for the DM with repulsive self-interaction. The distribution of DM in the DANS as a halo affects the path of X-rays emitted from hot spots on the visible baryonic surface, causing notable changes in the pulse profile observed by telescopes such as NICER, compared to pure baryonic NSs. We explore how various DM models affect the DM mass distribution, leading to different types of dark halos. We quantify the deviation in observed X-ray flux from stars with each of these halos. We identify the pitfalls in interpreting mass and radius measurements of NSs inferred from electromagnetic radiation and constraining the baryonic matter EOS if these dark halos exist.

Shawqi, Shafayat (ORCID:0000000210956183)↗

NASA NTRS1974

The electromagnetic interchange mode in a partially ionized collisional plasma

A collisional electromagnetic dispersion relation is derived from two-fluid theory for the interchange mode coupled to the Alfven, acoustic, drift and entropy modes in a partially ionized plasma. The fundamental electromagnetic nature of the interchange model is noted; coupling to the intermediate Alfven mode is strongly stabilizing for finite k sub z. Both ion viscous and ion-neutral stabilization are included, and it was found that collisions destroy the ion finite Larmor radius cutoff at short perpendicular wavelengths.

Hudson, M. K.↗

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↗

NASA NTRS1992

Time-dependent cosmic ray shocks with injection - A progress report

The study presents preliminary results of a two-fluid hydrodynamical simulation of cosmic-ray-modified shock structure by extending the steady-state calculations of Zank et al. (1991) to the time-dependent case. The present simulations model diffusive shock acceleration of the energetic particle component at planar shocks. A new model for the injection of cosmic ray particles is presented. The quadratic artificial viscosity algorithm is found to be numerically unstable and therefore unusable for strong cosmic ray shocks. There is visible ringing and overshoot in the various profiles. A large enhancement in the fluid density appears downstream of the subshock. This 'overcompression' results from the upstream fluid being compressed in the precursor, then traversing a strong shock with compression ratio near the adiabatic limit of four.

Donohue, D. J.↗

DOE OSTI2025

Modeling Volumetric Energy Deposition-driven Rayleigh-Taylor Instability

The goal of this work is to simulate a volumetric energy deposition (VED) driven Rayleigh-Taylor instability (RTI). Efforts will focus on simulating effects observed in two-fluid VED-driven RTI experiments. The VED triggers the RTI by changing the density of the fluids through time, hence a variable Atwood number is achieved. While initially stable (heavier fluid on the bottom), the microwave heating creates a density inversion and leads to a RTI. We hope to explore the mechanisms in VED-driven RTI and see how initial perturbations change the dynamics of growth. This work will primarily focus on early-time behavior but will also include steps towards simulating late-time nonlinear behavior.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1986

Thermofluid analysis of the SSME preburner using a gas-gas diffusion model for oxygen and hydrogen combustion at supercritical pressures

The paper discusses the thermofluid analysis of the Space Shuttle Main Engine (SSME) fuelside preburner. The governing equations have been solved numerically to predict flow, heat transfer, mixing, and combustion. A two-fluid approach is adopted in which oxygen is regarded as one fluid and hydrogen is regarded as the other fluid. The chemical kinetics is assumed to be very fast so that combustion is primarily controlled by the rate of mixing between oxygen and hydrogen. The preburner pressure is much greater than the critical pressures of oxygen and hydrogen; hence, a gas-gas diffusion model (rather than an evaporation model) has been developed to compute the rate of interphase mixing. Empirical correlations have been incorporated to account for the effect of slip on the interphase exchange. A sensitivity study has been performed with various model parameters. It is observed that the model can predict possibility of incomplete combustion and local regions of high temperatures under steady operating conditions. Some of these anomalies have been observed in actual tests, and the numerical model is useful for understanding possible causes and remedies. At least some measurements are needed for quantitative verification of the model.

Prakash, C.↗

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

Nonequilibrium effects in high-gain inertial confinement fusion

Recent experimental demonstrations of ignition and target gain in inertial confinement fusion (ICF) have stimulated interest in exploring the fundamental physics of violent deuterium-tritium (DT) burn in high-gain ICF targets. A significant DT-burn fraction is a necessary condition for high energy gain and large neutron yields (>100MJ). Using classical molecular-dynamics (MD) simulations and a hybrid fluid-kinetic model, we examine how a large fraction of low-energy 𝛼 particles can kick D and T ions out of equilibrium in high-gain ICF targets. The MD results suggest that (1) temperatures of 𝑇 𝐷 and 𝑇 𝑇 can differ by as much as ∼20% of their mean temperature and (2) the deviation of the DT energy distribution from the Maxwell-Boltzmann function can exceed ∼30%. Some of these MD observations, such as the preferential heating of D ions by low-energy 𝛼 particles and the temperature separation, can be explained by a proposed hybrid fluid-kinetic model. Furthermore, the implication of such nonequilibrium effects on the DT reactivity is also discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Planetary waves can activate resonant drag instabilities in 3D dusty gaseous discs

Resonant drag instabilities (RDIs) in protoplanetary discs are driven by the aerodynamic back-reaction of dust on gas and occur when the relative dust-gas motion resonate with a wave mode intrinsic to the gas fluid. Axisymmetric models indicate that the RDI generates filamentary perturbations, leading to grain clumping and planetesimal formation. Motivated by these findings, we investigate the dust–gas interaction in a non-axisymmetric inviscid protoplanetary disc with an embedded low-mass planet ($M_{\mathrm{p}}\in [0.3, 3] \, \rm M_{\oplus }$, here $\, \rm M_{\oplus }$ is the Earth mass). We conduct global 3D high-resolution two-fluid simulations, with the dust being parametrized by the Stokes number $\mathrm{St}\in [0.01,0.5]$. We find that planetary waves (PWs; also known as Rossby waves), which propagate along the downstream separatrices of the horseshoe region, resonate with the streaming motion and trigger the RDI. The consequent development of a global-scale filamentary dust distribution does not sensitively depend on the Stokes number, nor does it depend on the fast dust settling that takes place in an inviscid disc. The rapid onset of this instability, which is comparable to the dynamical orbital time-scale, suppresses the formation of asymmetric structures in the dust in the vicinity of the planet (such as dust voids and filaments). Additionally, we find that the dust feedback enables buoyancy resonances in an otherwise non-buoyant (globally isothermal) disc. Therefore, our results provide the first numerical evidence of RDIs generation driven by planetary waves.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1991

Oblique MHD cosmic-ray modified shocks: Two-fluid numerical simulations

We present the first results of time dependent, two-fluid, cosmic-ray (CR) modified, MHD shock simulations. The calculations were carried out with a new numerical code for 1-D ideal MHD. By coupling this code with the CR energy transport equation we can simulate the time-dependent evolution of MHD shocks including the acceleration of the CR and their feedback on the shock structures. We report tests of the combined numerical method including comparisons with analytical steady state results published earlier by Webb, as well as internal consistency checks for more general MHD CR shock structures after they appear to have converged to dynamical steady states. We also present results from an initial time dependent simulation which extend the parameter space domain of previous analytical models. These new results support Webb's suggestion that equilibrium oblique shocks are less effective than parallel shocks in the acceleration of CR. However, for realistic models of anisotropic CR diffusion, oblique shocks may achieve dynamical equilibrium on shorter timescale than parallel shocks.

Frank, Adam↗

NASA NTRS1984

The collisionless macroscopic Kelvin-Helmholtz instabiity. I Transverse electrostatic mode

The electrostatic limit of the Kelvin-Helmholtz instability in a collisionless plasma is investigated for a magnetic field perpendicular to the velocity flow. A two-fluid eigenvalue analysis of the transverse electrostatic Kelvin-Helmholtz instability is performed using the Chew-Goldberger-Low hydromagnetic equations with first-order finite Larmer radius corrections, establishing that the transverse electrostatic mode reduces to the incompressible transverse MHD mode in the limit of zero ion gyroradius. The linear growth rates are determined and a strong stabilization of the instability for (k-perpendicular) (rho sub i) of 0.2 or greater is observed. The nonlinear stage of the instability displays large vortices whose size is many times the initial velocity shear length. At saturation the simulation is dominated by the longest-wavelength mode permitted in the system. The simulation results can be understood in terms of a simple model which assumes that saturation occurs when the angular frequency of the vortex motion reaches the linear growth rate of the instability.

Prichett, P. L.↗

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 OSTI2025

Topological Mixed Valence Model for Twisted Bilayer Graphene

Song and Bernevig (SB) have recently proposed a topological heavy-fermion description of the physics of magic angle twisted bilayer graphene (MATBG), involving the hybridization of flat-band electrons with a relativistic conduction sea. Here, we explore the consequences of this model, seeking a synthesis of understanding drawn from heavy-fermion physics and MATBG experiments. Our work identifies a key discrepancy between measured and calculated on-site Coulomb interactions, implicating renormalization effects that are not contained in the current model. With these considerations in mind, we consider a SB model with a single, renormalized on-site interaction between the f electrons, containing a phenomenological heavy-fermion binding potential on the moiré A A sites. This feature allows the simplified model to capture the periodic reset of the chemical potential with filling and the observed stability of local moment behavior. We argue that a two-stage Kondo effect will develop in MATBG as a consequence of the relativistic conduction band: Kondo I occurs at high temperatures, establishing a coherent hybridization at the Γ points and a non-Fermi liquid of incoherent fermions at the moiré K points; at much lower temperatures, Kondo II leads to a Fermi liquid in the flat band. Utilizing an auxiliary-rotor approach, we formulate a mean-field treatment of MATBG that captures this physics, describing the evolution of the normal state across a full range of filling factors. By contrasting the relative timescales of phonons and valence fluctuations in bulk heavy-fermion materials with that of MATBG, we are led to propose a valley-polaron origin to the Coulomb renormalization and the heavy-fermion binding potential identified from experiment. We also discuss the possibility that the two-fluid, non-Fermi liquid physics of the relativistic Kondo lattice is responsible for the strange-metal physics observed in MATBG. Published by the American Physical Society 2025

Lau, Liam L. H. (ORCID:0000000166039088)↗