SEARCH · Search NASA

Results for “Electron 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 37 records · Page 2

All sources NASA NTRS DOE OSTI

NASA NTRS1984

Nonlinear evolution of the lower-hybrid drift instability

The results of simulations of the lower-hybrid drift instability in a neutral sheet configuration are described. The simulations use an implict formulation to relax the usual time step limitations and thus extend previous explicit calculations to weaker gradients, larger mass ratios, and long times compared with the linear growth time. The numerical results give the scaling of the saturation level, heating rates, resistivity, and cross-field diffusion and a demonstration by comparison with a fluid electron model that dissipation in the lower-hybrid drift instability is caused by electron kinetic effects.

Brackbill, J. U.↗

DOE OSTI2021

Self-consistent simulation of resistive kink instabilities with runaway electrons

A new fluid model for runaway electron simulation based on fluid description is introduced and implemented in the magnetohydrodynamics code M3D-C1, which includes self-consistent interactions between plasma and runaway electrons. The model utilizes the method of characteristics to solve the continuity equation for the runaway electron density with large convection speed, and uses a modified Boris algorithm for pseudo particle pushing. The model was employed to simulate magnetohydrodynamics instabilities happening in a runaway electron final loss event in the DIII-D tokamak. Nonlinear simulation reveals that a large fraction of runaway electrons get lost to the wall when kink instabilities are excited and form stochastic field lines in the outer region of the plasma. Plasma current converts from runaway electron current to Ohmic current, and get pinched at the magnetic axis. Here, given the good agreement with experiment, the simulation model provides a reliable tool to study macroscopic plasma instabilities in existence of runaway electron current, and can be used to support future studies of runaway electron mitigation strategies in ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1984

Numerical investigation of fluid models with full electron and proton thermal conduction equations for the quiet solar wind

The effect of using full equations (rather than Fourier law expressions) for the electron and proton thermal conductivities in conjunction with community, momentum, and temperature equations for the description of spherically symmetric quiet solar wind states is investigated numerically. For this purpose a time-dependent method is used, and steady state solutions in both subsonic and supersonic regimes between the sun and 5 AU are obtained and presented. The consequences for one-fluid and two-fluid models for the solar wind are summarized and discussed.

Cuperman, S.↗

NASA NTRS2019

Understanding the Role of 𝛼 Particles in Oblique Heliospheric Shock Oscillations

Recent observations by DSCOVR provide high temporal resolution (50 samples per second) magnetic vector field data that allows investigating the details of oblique heliospheric shock oscillations. It was found that some of these shocks exhibit magnetic oscillations, both downstream and upstream of the shock front. The DSCOVR/MAG magnetic field data are supplemented by an extensive database of low Mach number (M < 3) low‐β (<1) shock data observed by Wind albeit with lower temporal resolution. Motivated by the observations, we use the 2.5D hybrid model of the oblique shocks with α particles in addition to kinetic protons and electron fluid. We model the properties of the oblique shocks for a number of typical parameters found in observations and study the effects of the shock parameters and the relative α particle abundances on the properties of the shock magnetic field, density, and velocity oscillations. We find the α particles “surf” on the shock front and produce a wake of density oscillations. We examine the details of the phase space of the ions as well as the ion velocity distribution functions in various parts of the shock and study their nonthermal properties. We determine the effects of the α particle kinetic properties and abundances on the structure and dynamics of the shock downstream oscillations for a range of parameters relevant to low Mach number low‐β heliospheric shocks.

Ofman, L.↗

DOE OSTI2024

1D modeling of plasma streamers at ammonia-air flame conditions

Abstract Self-consistent 1D modeling of streamers in ammonia-oxygen-nitrogen-water mixtures has been performed in this work. A fluid model that includes species transport, electrostatic potential, and detailed chemistry was developed and verified. This model is then used to simulate the avalanche, streamer formation and propagation phases, driven by a nanosecond voltage pulse, at different thermochemical conditions derived from a 1D laminar premixed ammonia-air flame. The applicability of the Meek’s criterion in predicting the streamer inception location was successfully confirmed. Streamer formation and propagation duration were found to vary significantly with different thermochemical conditions, due to the difference in ionization rates. The thermochemical state also affected the breakdown characteristics which was tested by maintaining the background reduced electric field constant. Detailed kinetic analyses revealed the importance of O ( 1 D ) in the production of key radicals, such as O, OH, and NH 2 . Furthermore, the contributions of the dissociative electronic excitation of NH 3 towards the production of H and NH 2 radicals have also been reported. Spatial and temporal evolution of the electron energy loss fractions for various inelastic collision processes at different thermochemical states uncovered the input plasma energy spent of fuel dissociation and the large variability in the dominant processes during the avalanche and streamer propagation phases. The methodology and analyses reported in this work are key towards developing effective strategies for controlled nanosecond-pulsed non-equilibrium plasma sources used for ammonia ignition and flame stabilization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRSDate not supplied

Current Capabilities of AFRL’s Spacecraft Simulation Tool

Assessment of spacecraft integration issues is typically accomplished with a combination of numerical tools developed to simulate different regions with different key physics. For instance, a detailed study of spacecraft with an electric propulsion device requires a device model, a plume model, and a spacecraft charging model. AFRL’s in-house spacecraft simulation tool, TURF, is now capable of performing all of these calculations in one simulation. In addition, the plume simulation capability has been expanded significantly to improve speed and accuracy. Other upgrades in TURF include adaptive mesh refinement (AMR) and dynamic load balancing. All of these upgrades are covered in this paper by providing detailed descriptions or pointing to relevant references. This paper also presents three example simulations demonstrating its multi-physics multi-scale capability.

spacecraft simulation↗

NASA NTRS1968

Two-fluid model of the solar wind.

Solar wind studied with two fluid model, using proton and electron fluid type heat equations, equation of motion and continuity equation

Hartle, R. E.↗

NASA NTRS1972

Solar wind for a magnetized plasma with tensor plasma pressure.

A solar wind model for a magnetized solar wind is presented using one-fluid hydromagnetic equations with generalized polytrope equations of state for the two tensor components of the plasma pressure. Fluid and magnetic field variables are calculated at the Earth using certain boundary conditions at the Sun. The azimuthal velocity agrees with observed values and explains the considerable loss of angular momentum from the Sun by the solar wind. The results are good for most variables but suggest that a two-fluid model with electrons at higher temperatures and smaller temperature anisotropy ratios than the ions would give improved agreement for some quantities.

Tan, M.↗

NASA NTRS1981

Fluid model of plasma outside a hollow cathode neutralizer

The present study analyzes the capability of a fluid model of electron transport to explain observed properties of the external plasma of a hollow cathode neutralizer used to neutralize beams emerging from ion thrusters. Calculations reported here show that when the effective collision frequency in such a model is near the plasma frequency, the resulting electric potential and electron temperature variations are in qualitative agreement with values measured in the plume mode of the hollow cathode. Both theory and experiment show strong variations of temperature and potential within a few centimeters of the cathode orifice.

Parks, D. E.↗

NASA NTRS1993

A radiating one-dimensional current sheet configuration

The structure of the x-independent (one-dimensional) forced current sheet including a self consistent By component is investigated for the case of small normal field component, Bz/B0 much less than 1. A hybrid (kinetic ions, massless fluid electrons) simulation model is used to demonstrate that such a current sheet has a time-dependent structure which radiates incompressible Alfven waves with amplitude of the order of the asymptotic (lobe) field strength B0. The central density enhancement acts as the source of a propagating wavetrain in which Bx rotates into By and back again. One of the characteristic signatures of the radiating current sheet is the presence of a reversal in Bx (or By) without a corresponding increase in density.

Pritchett, P. L.↗

NASA NTRS2021

Oblique High Mach Number Heliospheric Shocks: The Role of α Particles

Spacecraft observations of heliospheric shocks often find oscillations in the magnetic field and density both, upstream and downstream. The downstream magnetic oscillations of oblique collisionless shocks were detected by Wind with 10.9 samples/s and DSCOVR spacecraft with high temporal resolution of 50 samples/s. The density oscillations associated with the shocks are also evident in proton and α particle density by Wind (with much lower temporal resolution). Recently, we have investigated low Mach number low-β oblique shock oscillations using satellite data and 2.5D hybrid modeling with electrons modeled as fluid and ions modeled as particles and found that α particles—an important component of heliospheric plasma—may affect considerably the downstream oscillations and the shock structure. The objective of the present study is to investigate the effects of α particles on high Mach number heliospheric shocks dynamics, oscillations, nonstationarity, and shock-front rippling. We extend the study to high Mach number shocks (M > 3), investigate several α particle typical densities, and compare the results for the various shock parameters. We model the effects of α particles on the shock ramp, wake, and downstream oscillation structure and the kinetic properties of proton and α particle velocity distributions at various locations downstream of the shocks. Using the 2.5D hybrid model we found that the modeled high Mach number quasi-perpendicular shock magnetic and density structures are significantly affected by α particles with typical solar wind relative abundances, suggesting that the observed high Mach number shocks are similarly affected by α particles.

Leon Ofman↗

DOE OSTI2025

Ponderomotive electron physics captured in a single-fluid extended magnetohydrodynamics model

The ponderomotive force, arising from the interaction between electromagnetic waves and plasma, plays a critical role in laser fusion, astrophysical plasmas, and laser diagnostics. Traditionally, modeling this force requires multi-fluid or particle-in-cell simulations due to its strong coupling to electron-scale dynamics. In this work, we demonstrate that a one-fluid, two-temperature extended magnetohydrodynamics (XMHD) model—augmented with a generalized Ohm's law (GOL) including electron inertia—can accurately reproduce key ponderomotive effects. We derive the ponderomotive force within this framework using a phasor-based approach and then validate its nonlinear manifestations through direct numerical simulations in the PERSEUS code, where steepening and density modulation phenomena typically associated with kinetic-scale models are reproduced. These results and prior work establish XMHD as a robust and efficient alternative for modeling nonlinear laser–plasma dynamics, bridging the gap between ideal MHD and fully kinetic approaches.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Fluid and hybrid simulations of the ionization instabilities in Hall thruster

Low-frequency axial oscillations in the range of 5–50 kHz stand out as a pervasive feature observed in many types of Hall thrusters. While it is widely recognized that the ionization effects play the central role in this mode, as manifested via the large-scale oscillations of neutral and plasma density, the exact mechanism(s) of the instabilities remain unclear. To gain further insight into the physics of the breathing mode and evaluate the role of kinetic effects, a one-dimensional time-dependent full nonlinear low-frequency model describing neutral atoms, ions, and electrons is developed in full fluid formulation and compared to the hybrid model in which the ions and neutrals are kinetic. Both models are quasi-neutral and share the same electron fluid equations that include the electron diffusion, mobility across the magnetic field, and the electron energy evolution. The ionization models are also similar in both approaches. Further, the predictions of fluid and hybrid simulations are compared for different test cases. Two main regimes are identified in both models: one with pure low-frequency behavior and the other one, where the low-frequency oscillations coexist with high-frequency oscillations in the range of 100–200 kHz, with the characteristic time scale of the ion channel fly-by time, 100–200 kHz. The other test case demonstrates the effect of a finite temperature of injected neutral atoms, which has a substantial suppression effect on the oscillation amplitude.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1991

Simulation of high-speed (orbital) releases of electron attachment materials in the ionosphere

The dynamics of ionospheric plasma irregularities produced by the release of electron attachment materials at orbital velocities across the geomagnetic field is studied. A two-dimensional electrostatic fluid model which includes electron attachment and mutual neutralization chemistry, self-consistent electric fields, and three-species transport is developed. Numerical simulations are performed to study the behavior at early and at late times after the release. At early times, of the order of or less than the attachment material neutral collision time, the negative ion cloud produced by the release may structure owing to the shear in the E x B velocity within the cloud. At high altitudes the cloud may bifurcate and form vortices on the back. At lower altitudes where ion-neutral collisional effects dominate, this structuring is suppressed. At late times, after a plasma depletion has formed due to neutralization chemistry, the cloud structures by the E x B interchange instability. Depending on the release altitude, the depletion structures by the collisional or inertial limit of this instability.

Scales, W. A.↗

DOE OSTI2021

L-H transition trigger physics in ITER-similar plasmas with applied n = 3 magnetic perturbations

The L–H transition power threshold PLH is observed to increase with applied n = 3 resonant magnetic perturbations (RMP) in ITER-similar-shape plasmas with balanced neutral beam torque injection in DIII-D. The increase is most pronounced with added electron–cyclotron heating: PLH increases with decreasing edge plasma collisionality as PLH/PLH-08 ~ (ν*)−0.5, where PLH-08 is the 2008 ITPA multi-machine power threshold scaling. This result raises concerns for H-mode access at low edge collisionality in ITER, where RMP may have to be applied before the L–H transition to safely suppress the first edge-localized mode. Non-axisymmetric modifications with RMP include a simultaneous reduction of the radial electric field (Er) well depth and E × B shear. This can be attributed to increasing edge toroidal co-current rotation, and is consistent with substantially increased local long-wavelength turbulence (measured via beam emission spectroscopy). At high RMP perturbation strength the edge electric field Er reverses sign locally (becomes positive), with changes in dominant turbulence modes. Edge magnetic stochasticity provides an attractive explanation of the observed modifications, and the observed changes in toroidal rotation and Er are consistent with a simple fluid model describing radial electron current flow along stochastic fieldlines. The observed collisionality dependence of the L-mode edge electric field with applied RMP is also qualitatively consistent with this model. Reflectometry data indicate a significant reduction of the normalized L-mode radial density gradient a/Ln at high RMP field with simultaneous increase in radial particle flux and electron thermal flux from power balance analysis. We conjecture that the increase of PLH with RMP results from the combined effects of reduced E × B flow shear (increasing turbulent transport levels) and toroidal/poloidal flow modulation due to edge stochasticity. Initial experiments indicate that non-resonant n = 3 magnetic perturbations lead only to relatively small changes in Er, E × B shear and fluctuation characteristics, and have less impact on the L–H transition power threshold. This motivates further exploration of the RMP spectrum dependence of PLH for possible mitigation of the observed threshold increase.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Electron cyclotron drift instability and anomalous transport: two-fluid moment theory and modeling

In the presence of a strong electric field perpendicular to the magnetic field, the electron cross-field (E × B) flow relative to the unmagnetized ions can cause the so-called electron cyclotron drift instability (ECDI) due to resonances of the ion acoustic mode and the electron cyclotron harmonics. This occurs in, for example, collisionless shock ramps in space, and in E × B discharge devices such as Hall thrusters. A prominent feature of ECDI is its capability to induce an electron flow parallel to the background >E field at a speed greatly exceeding predictions by classical collision theory. Such anomalous transport is important due to its role in particle thermalization at space shocks, and in causing plasma flows towards the walls of E × B devices, leading to unfavorable erosion and performance degradation, etc. The development of ECDI and anomalous transport is often considered requiring a fully kinetic treatment. In this work, however, we demonstrate that a reduced variant of this instability, and more importantly, the associated anomalous transport, can be treated self-consistently in a collisionless two-fluid framework without any adjustable collision parameter. By treating both electron and ion species on an equal footing, the free energy due to the inter-species velocity shear allows the growth of an anomalous electron flow parallel to the background E field. We will first present linear analyses of the instability in the two-fluid five- and ten-moment models, and compare them against the fully-kinetic theory. At low temperatures, the two-fluid models predict the fastest-growing mode in good agreement with the kinetic result. Also, by including more (> = 10) moments, secondary (and possibly higher) unstable branches can be recovered. The dependence of the instability on ion-to-electron mass ratio, plasma temperature, and background B field strength is also thoroughly explored. We then carry out direct numerical simulations of the cross-field setup using the five-moment model. The development of the instability, as well as the anomalous transport, is confirmed and in excellent agreement with theoretical predictions. The force balance properties are also studied using the five-moment simulation data. Here this work casts new insights into the nature of ECDI and the associated anomalous transport and demonstrates the potential of the two-fluid moment model in efficient modeling of E × B plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Chaotic and integrable magnetic fields in one-dimensional hybrid Vlasov–Maxwell equilibria

The construction of kinetic equilibrium states is important for studying stability and wave propagation in collisionless plasmas. Thus, many studies over the past decades have been focused on calculating Vlasov–Maxwell equilibria using analytical and numerical methods. However, the problem of kinetic equilibrium of hybrid models is less studied, and self-consistent treatments often adopt restrictive assumptions ruling out cases with irregular and chaotic behaviour, although such behaviour is observed in spacecraft observations of space plasmas. In this paper, we develop a one-dimensional (1-D), quasineutral, hybrid Vlasov–Maxwell equilibrium model with kinetic ions and massless fluid electrons and derive associated solutions. The model allows for an electrostatic potential that is expressed in terms of the vector potential components through the quasineutrality condition. The equilibrium states are calculated upon solving an inhomogeneous Beltrami equation that determines the magnetic field, where the inhomogeneous term is the current density of the kinetic ions and the homogeneous term represents the electron current density. We show that the corresponding 1-D system is Hamiltonian, with position playing the role of time, and its trajectories have a regular, periodic behaviour for ion distribution functions that are symmetric in the two conserved particle canonical momenta. For asymmetric distribution functions, the system is nonintegrable, resulting in irregular and chaotic behaviour of the fields. The electron current density can modify the magnetic field phase space structure, inducing orbit trapping and the organization of orbits into large islands of stability. Thus, the electron contribution can be responsible for the emergence of localized electric field structures that induce ion trapping. We also provide a paradigm for the analytical construction of hybrid equilibria using a rotating two-dimensional harmonic oscillator Hamiltonian, enabling the calculation of analytic magnetic fields and the construction of the corresponding distribution functions in terms of Hermite polynomials.

Kaltsas, Dimitrios A. (ORCID:0000000300769015)↗

DOE OSTI2020

A Plasma Modeling Hierarchy and Verification Approach

This report reviews a hierarchy of formal mathematical models for describing plasma phenomena. Starting with the Boltzmann equation, a sequence of approximations and modeling assumptions can be made that progressively reduce to the equations for magnetohydrodynamics. Understanding the assumptions behind each of these models and their mathematical form is essential to appropriate use of each level of the hierarchy. A sequence of moment models of the Boltzmann equation are presented, then focused into a generalized three-fluid model for neutral species, electrons, and ions. This model is then further reduced to a two-fluid model, for which Braginskii described a useful closure. Further reduction of the two-fluid model yields a Generalized Ohm's Law model, which provides a connection to magnetohydrodynamic approaches. A verification approach based on linear plasma waves is presented alongside the model hierarchy, which is intended as an initial and necessary but not sufficient step for verification of plasma models within this hierarchy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗