Search NASA⌕ Search

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.

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

At least 91 records · Page 5

Collisionless ion-electron energy exchange in magnetized shocks

Energy partition between ions and electrons in collisionless shocks has been a long-standing unsolved fundamental physical question. Here, we show that kinetic simulations of moderate Alfv´enic Mach number, magnetized, collisionless shocks reveal rapid, faster-than-Coulomb, energy exchange between ions and electrons when the plasma is sufficiently magnetized. Using kinetic and multi-fluid models with counter-streaming ions, we identify resonances between electron whistler and ion magnetohydrodynamic waves that account for this rapid energy exchange.

High-energy-density plasmas↗

The kinetic analog of the pressure–strain interaction

Energy transport in weakly collisional plasma systems is often studied with fluid models and diagnostics. However, the applicability of fluid models is limited when collisions are weak or absent, and using a fluid approach can obscure kinetic processes that provide key insights into the physics of energy transport. Kinetic diagnostics retain all of the information in 3D-3V phase space and thereby reach beyond the insights of fluid models to elucidate the mechanisms responsible for collisionless energy transport. In this work, we derive the Kinetic Pressure–Strain (KPS): a kinetic analog of the pressure–strain interaction, which is the channel between flow energy density and internal energy density in fluid models. Through two case studies of electron Landau damping, we demonstrate that the KPS diagnostic can elucidate kinetic mechanisms that are responsible for energy transport in this channel, just as the related field–particle correlation is known to identify kinetic mechanisms of transport between electromagnetic field energy density and kinetic energy density in particle flows. In addition, we show that resonant electrons play a major role in transferring energy between fluid flows and internal energy during the process of Landau damping.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A drift kinetic model for the expander region of a magnetic mirror

Here, we present a drift kinetic model for the free expansion of a thermal plasma out of a magnetic nozzle. This problem relates to plasma space propulsion systems, natural environments such as the solar wind, and end losses from the expander region of mirror magnetically confined fusion concepts such as the gas dynamic trap. The model incorporates trapped and passing orbit types encountered in the mirror expander geometry and maps to an upstream thermal distribution. This boundary condition and quasineutrality require the generation of an ambipolar potential drop of ~5T e /e, forming a thermal barrier for the electrons. The model for the electron and ion velocity distributions and fluid moments is confirmed with data from a fully kinetic simulation. Finally, the model is extended to account for a population of fast sloshing ions arising from neutral beam heating within a magnetic mirror, again resulting in good agreement with a corresponding kinetic simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Three dimensional simulation of the operation of a hollow cathode electron emitter on the Shuttle orbiter

Several researchers have suggested using hollow cathodes as plasma contactors for electrodynamic tethers, particularly to prevent the shuttle orbiter from charging to large negative potentials. Previous studies have shown that fluid models with anomalous scattering can describe the electron transport in hollow cathode generated plasmas. An improved theory of the hollow cathode plasmas is developed and computational results using the theory are compared with laboratory experiment. Numerical predictions for a hollow cathode plasma source of the type considered for use on the shuttle are presented as are three-dimensional NASCAP/LEO calculations of the emitted ion trajectories and the resulting potentials in the vicinity of the orbiter. The computer calculations show that the hollow cathode plasma source makes vastly superior contact with the ionospheric plasma compared with either an electron gun or passive ion collection by the orbiter.

Davis, V. A.↗

Implications of a two-fluid heat diffusion model for RF condensation effects in magnetic islands

Theoretical calculations have predicted an “RF condensation” effect through which the temperature perturbation and current drive may be nonlinearly enhanced in electron cyclotron (EC) heated magnetic islands. Previously, two-fluid effects had been neglected in the OCCAMI code used to model RF condensation. For realistic experimental conditions, we show that two-fluid effects in an EC heated plasma may have a major impact on RF condensation. Less collisional plasmas, where electron and ion temperatures are decoupled, exhibit an enhanced linear temperature perturbation in the island, and an even stronger nonlinear enhancement of the perturbed temperature. This can occur even when nonlinear effects are small in the single-fluid limit. We also calculate corresponding hysteresis effects. Finally, we demonstrate that other physical effects, including the rotation of the magnetic island and the Gaussian spread of the beam from a gyrotron, can significantly decrease the predicted observable RF condensation effect.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Electron heating in superhigh Mach number shocks

Fluid and MHD models, as well as direct extrapolation of the earth's bow shock measurements in the high Mach number (HMN) range to the superhigh Mach number (SHMN) range predict that the downstream electron pressure is only a negligible fraction of the Rankine-Hugoniot downstream pressure. Following Alfven, plasma physics experimental-theoretical data combined with magnetospheric observations were used to probe the physics of the SHMN shocks. It is shown that inclusion of proper plasma physics considerations in the interaction of the reflected and transmitted ions and the electrons at the 'foot' of the shock leads to the surprising result that electron heating can dominate in the SHNM range. A stationary model of the shock structure is derived and shown to be the result of extrapolation of the high Mach number shock physics wiht incorporation of collective interactions at the foot.

Papadopoulos, K.↗

Particle and fluid simulations of resistive current-driven electrostatic ion cyclotron waves

The results from 1-D numerical simulations of electrostatic ion cyclotron waves (EIC) are presented for a model in which the electrons are a resistive (collisional) fluid. Simulations of both the kinetic and fluid descriptions are performed and compared in order to assess the fundamental limitations of fluid theory for EIC waves. The effect of ion-neutral collisions is also included using a simple Monte Carlo technique. It is found that a small ion-neutral collision frequency destroys the frequency harmonic coupling of kinetic EIC waves and tends to validate the fluid description. The saturation amplitude of the current-driven EIC instability is in agreement with recent laboratory experiments. The coherent nature (extremely narrow spectral width) and phase velocity agree with ground based (coherent backscatter radars) and in situ observations of current-driven EIC waves in the high latitude ionosphere.

Seyler, Charles E.↗

Understanding the depletion of electrons in dusty plasmas at atmospheric pressure

The nucleation and growth of NPs in the gas phase using atmospheric-pressure plasma systems is an important approach to synthesizing novel dimensionally-controlled materials. In this paper, we investigated the effect of the NPs on a typical type of continuous-flow, substrate-free plasma at atmospheric pressure to understand their potential contribution to electron density changes. A tandem plasma system was set up consisting of one plasma reactor that produced carbonaceous NPs from mixtures of argon and hexane, and another identical plasma reactor where the as-grown particles were injected and non-intrusive electrical and optical measurements were performed. The electron densities obtained from conductivity measurements and a plasma fluid model were found to decrease in the presence of NPs. However, control experiments revealed that the main source of the electron depletion was residual vapor or reaction byproducts in the form of molecular species or nanoclusters and not the particles themselves. These results were validated by constant number Monte Carlo simulations which showed that at the experimentally-measured conditions, the NPs were not of sufficiently high enough concentration to reduce the electron density; however, if molecules or clusters are ionizable, they remain in sufficient concentration to deplete electron densities. Our study shows that at atmospheric pressure, because of their typically larger electron density values, particle-producing plasmas are distinct from those at low pressure, and nanoparticle formation does not have the same impact while molecular-scale species may be a more important consideration.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

On rotational discontinuities in both two-fluid and hybrid models

Rotational discontinuities are studied in a two-fluid model that includes finite ion inertia dispersion and in a hybrid model in which the full ion dynamics is retained while the electrons are treated as a massless fluid. It is shown that as in previous dissipative MHD studies, a rotational discontinuity is unstable in both models and evolves to a 2-3 intermediate shock, a slow rarefaction wave, and other waves. In addition, it is shown that the so-called Walen relation, which holds exactly for rotational discontinuities, can also be well satisfied by tran-Alfvenic intermediate shocks. Thus intermediate shocks can be candidates for those observed structures that satisfy the Walen relation.

Wu, C. C.↗

Role of the edge stochastic layer in density pump-out by resonant magnetic perturbations

The role of the edge stochastic layer on particle transport is addressed in DIII-D plasmas with applied resonant magnetic perturbations (RMPs) causing density pump-out (a term used to denote the density reduction at the top of the pedestal due to the applied RMPs). Using an analytical model that adds the stochastic parallel transport of electrons (Harvey et al 1981 Phys. Rev. Lett. 47 102) to the fluid equations, the ambipolar radial electric field and particle flux are calculated simultaneously. In this model the nonambipolar electron particle flux, driven by the stochastic magnetic field, is predominantly balanced by the nonambipolar perpendicular ion flux, driven by toroidal viscosity, across a narrow stochastic layer of order two percent of the plasma radius (0.98 < ψn < 1). The model reproduces the level of density pump-out and its dependence on RMP amplitude near the separatrix for three plasma discharges for which the pedestal foot is at medium to high collisionality $(2\leqslant {\nu }_{\text{sep}}^{\ast }\leqslant 30)$. In this work, the experimentally observed and numerically tested inverse density dependence of density pump-out (Hu et al 2020 Nucl. Fusion 60 076001) is accurately captured by the model: an increase in collisionality with density in the pedestal foot results in a decrease in stochastic diffusivity, and hence a decrease in the level of pump-out.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Establishing criteria for the transition from kinetic to fluid modeling in hollow cathode analysis

Hollow cathodes for plasma switch applications are investigated via 2D3V particle-in-cell simulations of the channel and plume region. The kinetic nature of the plasma within the channel is dependent on the thermalization rate of electrons, emitted from the insert. When Coulomb collisions occur at a much greater rate than ionization or excitation collisions, the electron energy distribution function rapidly relaxes to a Maxwellian and the plasma within the channel can be described accurately via a fluid model. In contrast, if inelastic processes are much faster than Coulomb collisions, then the electron energy distribution function in the channel exhibits a notable high-energy tail, and a kinetic treatment is required. This criterion is applied to hollow cathodes from the literature, revealing that a fluid approach is suitable for most electric propulsion applications, whereas a kinetic treatment can be more critical to accurate modeling of plasma switches.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Magnetofluid Simulations of the Global Solar Wind Including Pickup Ions and Turbulence Modeling

I will describe a three-dimensional magnetohydrodynamic model of the solar wind that takes into account turbulent heating of the wind by velocity and magnetic fluctuations as well as a variety of effects produced by interstellar pickup protons. In this report, the interstellar pickup protons are treated as one fluid and the protons and electrons are treated together as a second fluid. The model equations include a Reynolds decomposition of the plasma velocity and magnetic field into mean and fluctuating quantities, as well as energy transfer from interstellar pickup protons to solar wind protons that results in the deceleration of the solar wind. The model is used to simulate the global steady-state structure of the solar wind in the region from 0.3 to 100 AU. Where possible, the model is compared with Voyager data. Initial results from generalization to a three-fluid model is described elsewhere in this session.

Goldstein, Melvyn L.↗

Verification of gyrokinetic particle simulations of neoclassical tearing modes in fusion plasmas

Abstract The ability to simulate neoclassical tearing modes (NTMs) in the gyrokinetic toroidal code (GTC) has been developed and verified, in which ions are treated with a gyrokinetic model and electrons are treated as a resistive fluid. The simulation results demonstrate that the neoclassical bootstrap current effect can destabilize an otherwise stable classical tearing mode. In the cylindrical geometry, GTC simulations in the magnetohydrodynamic limit show quantitative agreement with the modified Rutherford theory, both in terms of the scaling law in the small island limit and in terms of the saturation level and pressure flattening effect in the large island limit. The toroidal effects are slightly destabilizing for the NTM, while the kinetic effects of thermal ions are stabilizing for the NTM and increase its excitation threshold.

Physics↗

Propagation of acoustic modes in the transitional ionosphere

Continuum theory is unsatisfactory for describing wave propagation in the transitional ionosphere, since the mean free path is of the same order as the characteristic length (i.e., Knudsen number in this region is of the order of unity). Therefore a transition model in which isotropic electrons behave as a fluid while the ions are governed by the kinetic equation modified by ion-electron collision effects is proposed to study acoustic wave propagation in the upper atmosphere in the altitude region from 500 to 2000 km. The results show that the dissipation of acoustic waves by electrons is through the viscous and thermal conduction effects, and the dissipation of acoustic waves by ions is through Landau damping. Comparisons of the characteristics of acoustic mode propagation and wave dissipation mechanisms in the collision-dominated, transitional from collisional to collisionless, and collisionless media are also discussed.

Hung, R. J.↗

A mathematical model of the structure and evolution of small scale discrete auroral arcs

A three dimensional fluid model which includes the dispersive effect of electron inertia is used to study the nonlinear macroscopic plasma dynamics of small scale discrete auroral arcs within the auroral acceleration zone and ionosphere. The motion of the Alfven wave source relative to the magnetospheric and ionospheric plasma forms an oblique Alfven wave which is reflected from the topside ionosphere by the negative density gradient. The superposition of the incident and reflected wave can be described by a steady state analytical solution of the model equations with the appropriate boundary conditions. This two dimensional discrete auroral arc equilibrium provides a simple explanation of auroral acceleration associated with the parallel electric field. Three dimensional fully nonlinear numerical simulations indicate that the equilibrium arc configuration evolves three dimensionally through collisionless tearing and reconnection of the current layer. The interaction of the perturbed flow and the transverse magnetic field produces complex transverse structure that may be the origin of the folds and curls observed to be associated with small scale discrete arcs.

Seyler, C. E.↗

Simulations of a molecular plasma in collisional-radiative nonequilibrium

A code for the simulation of nonequilibrium plasmas is being developed, with the capability to couple the plasma fluid-dynamics for a single fluid with a collisional-radiative model, where electronic states are treated as separate species. The model allows for non-Boltzmann distribution of the electronic states. Deviations from the Boltzmann distributions are expected to occur in the rapidly ionizing regime behind a strong shock or in the recombining regime during a fast expansion. This additional step in modeling complexity is expected to yield more accurate predictions of the nonequilibrium state and the radiation spectrum and intensity. An attempt at extending the code to molecular plasma flows is presented. The numerical techniques used, the thermochemical model, and the results of some numerical tests are described.

Cambier, Jean-Luc↗

Fluid-kinetic modeling of a high power density radio frequency inductively coupled positive hydrogen ion source

High power density radio-frequency (RF) inductively coupled positive ion sources are attractive candidates for next-generation neutral beam injection (NBI) systems, where higher injected power and longer pulse lengths are desired without sacrificing source reliability. Operating at absorbed power densities of order $\gt 1~\mathrm{W\,cm}^{-3}$ places these sources in a regime with stronger gas heating, higher dissociation, and non-Maxwellian electron energy distributions. The Large Uniform Plasma for Ionizing Neutrals (LUPIN) is an RF inductively coupled plasma source designed to explore this high power density regime and to provide guidance for a positive ion source upgrade for the DIII-D NBI system. LUPIN is designed to operate at up to 20 kW of RF power at 2 MHz, coupling energy through a cylindrical quartz vessel to achieve target ion current densities of $2100\,\mathrm{A\,m}^{-2}$ . This paper presents fluid-kinetic modeling of LUPIN using the hybrid plasma equipment model where electrons are treated kinetically, and the simulations reveal that electron energy distribution function transitions from nearly Maxwellian in the core to bi-Maxwellian towards the edge. Parametric simulations investigate the effects of RF power, gas pressure, and frequency on plasma density, ion flux, and uniformity. Parametric sweeps reveal that increasing power shifts the primary ionization channel from molecular to atomic with diminishing flux gains due to skin-depth contraction and gas rarefaction. Higher frequency localizes heating and increases $\mathrm{H}_2^+$ and $\mathrm{H}_3^+$ delivery to the grid, while elevated pressure boosts ionization yet hinders ion transport due to increase in collisionality.

inductively coupled plasma↗

Investigation of the Polytropic Relationship Between Density and Temperature Within Interplanetary Coronal Mass Ejections Using Numerical Simulations

Single-point spacecraft measurements within coronal mass ejections (CMEs) often exhibit a negative correlation between electron density and temperature. At least two opposing interpretations have been suggested for this relationship. If, on one hand, these single spacecraft observations provide direct measures of the polytropic properties of the plasma, then they imply that the polytropic index for the electrons gamma(sub e) is often < 1. Moreover, since the electrons carry the bulk of the pressure (via their significantly higher temperature), this further implies that the dynamics of CME evolution are dominated by an effective polytropic index gamma(sub e)ff < 1. On the other hand, gamma < 1 implies that as the ejecta propagate away from the Sun and expand, they also heat up; a result clearly at odds with in situ observations. In contrast to these CME intervals, many studies have shown that the quiescent solar wind exhibits a positive correlation between electron density and temperature, suggesting that gamma(sub e) > 1. In this study we simulate the evolution of a variety of CME-like disturbances in the solar wind using a one-dimensional, single-fluid model, to address the interpretation of the relationship between electron density and temperature within CMEs at fixed locations in space. Although we strictly impose a polytropic relationship (with gamma = constant) throughout our simulations, we demonstrate that a variety of correlations can exist between density and temperature at fixed points. Furthermore, we demonstrate that the presence of only local uncorrelated random fluctuations in density and temperature can produce a negative correlation. Consequently, we conclude that these single-point observations of negative correlations between electron density and temperature cannot be used to infer the value of gamma(sub e). Instead, we suggest that entropy variations, together with the plasma's tendency to achieve pressure balance with its surroundings, are responsible for the observed profiles.

Riley, Pete↗