Search NASA⌕ Search

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.

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

At least 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
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 NTRS1994

Characteristics of Vaporizing Cryogenic Sprays for Rocket Combustion Modeling

Experimental measurements of the volume-median drop diameter, Dv.5e, of vaporizing cryogenic sprays were obtained with a drop size measuring instrument developed at NASA Lewis Research Center. To demonstrate the effect of atomizing-gas properties on characteristic drop size, a two-fluid fuel nozzle was used to break up liquid-nitrogen, LN2, jets in high-velocity gasflows of helium argon and gaseous nitrogen, GN2. Also, in order to determine the effect of atomizing-gas temperature on specific surface-areas of LN2 sprays, drop size measurements were made at gas temperatures of 111 and 293 K.

Ingebo, Robert D.↗
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.↗
NASA NTRS1994

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

We present the first results of time-dependent two-fluid cosmic-ray (CR) modified magnetohydrodynamic (MHD) shock simulations. The calculations were carried out with a new numerical code for one-dimensional 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 apppear to have converged to dynamical steady states. We also present results from an initial time-dependent simulation which extends 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 timescales than parallel shocks.

Frank, Adam↗
NASA NTRS1994

Particle/fluid simulations of an eruptive flare: Identifying the field-aligned currents responsible for the hard x-rays

While magnetohydrodynamics (MHD) can provide a reasonable description of the overall magnetic reconnection that is believed to drive flares, additional, and often separate processes have to be envoked to in order to explain the electron acceleration that is responsible for many of the observed flare emissions. A new model that incorporates the dynamic coronal current sheets, the reconnection site, and possible electron acceleration processes is developed through the use of two-dimensional particle and modified two-fluid simulations. The specific example of an eruptive flare driven by the coalescence of flux tubes supported by prescribed photospheric current elements is evaluated. It is shown that the electrons and ions have differential trajectories through the coronal current sheet which leads to the development of additonal plasma currents that flow around the surface of the current sheet. These surface currents are explicitly neglected in MHD but they are vital to the flare dynamics because they divert current from the coronal current sheet into the chromosphere, producing an effective resistivity that aids the development of fast reconnection. Because the surface currents are in the plane of the magnetic field, electrons in them experience strong acceleration and can account for the observed hard X-ray emissions. Model predictions are compared with observed time profiles of hard X-ray emissions and Doppler shifts seen in soft X-ray line emissions and are able to account for such features as (1) the asymmetry in the rise and decay time of the hard X-rays, (2) the apparent delay between the largest Doppler shifts and the hard X-ray peak, and (3) the relatively low intensity of the blue-shifted component. The use of particle and fluid simulations is important because it provides different, but complementary treatments of the electron acceleration, the global magnetic morphology, and the flare current system.

Winglee, R. M.↗
NASA NTRS1994

Time-dependent evolution of cosmic-ray-modified shock structure: Transition to steady state

Steady state solutions to the two-fluid equations of cosmic-ray-modified shock structure were investigated first by Drury and Volk (1981). Their analysis revealed, among other properties, that there exist regions of upstream parameter space where the equations possess three different downstream solutions for a given upstream state. In this paper we investigate whether or not all these solutions can occur as time-asymptotic states in a physically realistic evolution. To do this, we investigate the time-dependent evolution of the two-fluid cosmic-ray equations in going from a specified initial condition to a steady state. Our results indicate that the time-asymptotic solution is strictly single-valued, and it undergoes a transition from weakly to strongly cosmic-ray-modified at a critical value of the upstream cosmic ray energy density. The expansion of supernova remnant shocks is considered as an example, and it is shown that the strong to weak transition is in fact more likely. The third intermediate solution is shown to influence the time-dependent evolution of the shock, but it is not found to be a stable time-asymptotic state. Timescales for convergence to these states and their implications for the efficiency of shock acceleration are considered. We also investigate the effects of a recently introduced model for the injection of seed particles into the shock accelerated cosmic-ray population. The injection is found to result in a more strongly cosmic-ray-dominated shock, which supports our conclusion that for most classes of intermediate and strong cosmic-ray-modified shocks, the downstream cosmic-ray pressure component is at least as large as the thermal gas pressure, independent of the upstream state. As a result, cosmic rays almost always play a significant role in determining the shock structure and dissipation and they cannot be regarded as test particles.

Donohue, D. J.↗
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)↗
NASA NTRS1972

A solar-wind model including proton thermal anisotropy.

A relatively simple model is proposed to study the magnetohydrodynamic expansion of the solar wind. The interplanetary space is divided into an 'inner region' and an 'outer region.' In the inner region the solar wind is one-fluid and thermally isotropic. The solar-wind protons are assumed to be completely collisionless in the outer region. Thus the solar wind is two-fluid in nature, and the protons are thermally anisotropic in the outer region. It is assumed that throughout the interplanetary space electrons are thermally isotropic due to collisions or other processes. The electron heat flux is supposed to obey the well-known relationships that the flux is parallel to the magnetic-field vector and is directly proportional to the spatial derivative of the electron temperature.

Whang, Y. C.↗