SEARCH · Engineering Papers

Results for “collisional plasma”

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 631 records · Page 35

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Catastrophic Cooling Instability in Optically Thin Plasmas

The solar corona is the prototypical example of a low-density environment heated to high temperatures by external sources. The plasma cools radiatively, and because it is optically thin to this radiation, it becomes possible to model the density, velocity, and temperature structure of the system by modifying the MHD equations to include an energy source term that approximates the local heating and cooling rates. The solutions can be highly inhomogeneous and even multiphase because the well-known linear instability associated with this source term, thermal instability, leads to a catastrophic heating and cooling of the plasma in the nonlinear regime. Here we show that there is a separate, much simpler linear instability accompanying this source term that can rival thermal instability in dynamical importance. The stability criterion is the isochoric one identified by Parker (1953), and we demonstrate that cooling functions derived from collisional ionization equilibrium are highly prone to violating this criterion. If catastrophic cooling instability can act locally in global simulations, then it is an alternative mechanism for forming condensations, and due to its nonequilibrium character, it may be relevant to explaining a host of phenomena associated with the production of cooler gas in hot, low density plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Parallel expansion of a fuel pellet plasmoid

The problem of the assimilation of a cryogenic fuel pellet injected into a hot plasma is considered. Due to the transparency to ambient particles of the plasmoid, the localised region of high-density plasma created by ionisation of the ablated pellet material, electrons reach a ‘quasiequilibrium’ (QE) state which is characterised by a steady-state on the fastest collisional time scale. The simplified electron kinetic equation of the QE state is solved. Taking a velocity moment of the higher-order electron kinetic equation, which is valid on the expansion time scale, permits a fluid closure, yielding an evolution equation for the macroscopic parameters describing the QE distribution function. In contrast to the Braginskii equations, the closure does not require that electrons have a short mean free path compared with the size of density perturbations, and permits an anisotropic and highly non-Maxwellian distribution function. As the QE distribution function accounts for both trapped and passing electrons, the self-consistent electric potential that causes the expansion can be properly described, in contrast to earlier models of pellet plasmoid expansion with an unbounded potential. The plasmoid expansion is simulated using both a Vlasov model and a cold-fluid model for the ions. During the expansion plasmoid ions and electrons obtain nearly equal amounts of energy; as hot ambient electrons provide this energy in the form of collisional heating of plasmoid electrons, the expansion of a pellet plasmoid is expected to be a potent mechanism for the transfer of energy from electrons to ions on a time scale shorter than that of ion–electron thermalisation.

Physics↗

DOE OSTI2024

The impact of collisionality on the runaway electron avalanche during a tokamak disruption

The exponential growth (avalanching) of runaway electrons (REs) during a tokamak disruption continues to be a large uncertainty in RE modeling. The present work investigates the impact of tokamak geometry on the efficiency of the avalanche mechanism across a broad range of disruption scenarios. It is found that the parameter ν *,crit ⁠, describing the collisionality at the critical energy to run away, delineates how toroidal geometry impacts RE formation. In particular, utilizing a reduced but self-consistent description of plasma power balance, it is shown that for a high-density deuterium-dominated plasma, ν *,crit is robustly less than one, resulting in a substantial decrease in the efficiency of the RE avalanche compared to predictions from slab geometry. In contrast, for plasmas containing a substantial quantity of neon or argon, ν *,crit ≳1, no reduction of the avalanche is observed due to toroidal geometry. This sharp contrast in the impact of low- vs high-Z material results primarily from the relatively strong radiative cooling from high-Z impurities, enabling the plasma to be radiatively pinned at low temperatures and thus large electric fields, even for modest quantities of high-Z material.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Regulation of Solar Wind Electron Temperature Anisotropy by Collisions and Instabilities

Abstract Typical solar wind electrons are modeled as being composed of a dense but less energetic thermal “core” population plus a tenuous but energetic “halo” population with varying degrees of temperature anisotropies for both species. In this paper, we seek a fundamental explanation of how these solar wind core and halo electron temperature anisotropies are regulated by combined effects of collisions and instability excitations. The observed solar wind core/halo electron data in ( β ∥ , T ⊥ / T ∥ ) phase space show that their respective occurrence distributions are confined within an area enclosed by outer boundaries. Here, T ⊥ / T ∥ is the ratio of perpendicular and parallel temperatures and β ∥ is the ratio of parallel thermal energy to background magnetic field energy. While it is known that the boundary on the high- β ∥ side is constrained by the temperature anisotropy-driven plasma instability threshold conditions, the low- β ∥ boundary remains largely unexplained. The present paper provides a baseline explanation for the low- β ∥ boundary based upon the collisional relaxation process. By combining the instability and collisional dynamics it is shown that the observed distribution of the solar wind electrons in the ( β ∥ , T ⊥ / T ∥ ) phase space is adequately explained, both for the “core” and “halo” components.

Yoon, Peter H. (ORCID:0000000181343790)↗

DOE OSTI2021

RF sheath induced sputtering on Proto-MPEX part 2: Impurity transport modeling and experimental comparison

The Prototype Material Plasma Exposure eXperiment (Proto-MPEX) is a pulsed linear plasma device that uses a radio frequency (RF) helicon antenna and a helicon ceramic window as the plasma ionization source. Modeling performed within this work on the helicon operations predicts hot spots on the helicon ceramic window due to an RF-induced sheath. The RF sheath potential leads to physical sputtering of the ceramic window predominantly where the helicon antenna sits around the helicon window. The role of oxygen in sputtering within Proto-MPEX is also investigated and is needed to match the experimental data. Experiments on Proto-MPEX show reasonable agreement with the modeled erosion/redeposition pattern on the helicon window, and the modeled impurity flux to the target shows radial transport within the device as expected from classical collisional impurity transport. Temperature screening of impurities is shown to be needed for Proto-MPEX to reduce impurity flux to the target. The modeling approach used in this work can be applied to other plasma RF applications for the determination of impurity production and transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Multi-fluid and kinetic models of partially ionized magnetic reconnection

Magnetic reconnection in partially ionized plasmas is a ubiquitous and important phenomenon in both laboratory and astrophysical systems. Here, simulations of partially ionized magnetic reconnection with well-matched initial conditions are performed using both multi-fluid and fully-kinetic approaches. Despite similar initial conditions, the time-dependent evolution differs between the two models. In multi-fluid models, the reconnection rate locally obeys either a decoupled Sweet–Parker scaling, where neutrals are unimportant, or a fully coupled Sweet–Parker scaling, where neutrals and ions are strongly coupled, depending on the resistivity. In contrast, kinetic models show a faster reconnection rate that is proportional to the fully-coupled, bulk Alfvén speed, $v^*_A$. In this work, these differences are interpreted as the result of operating in different collisional regimes. Multi-fluid simulations are found to maintain $ν_{ni}L/v^*_A$ ≳1, where $ν_{ni}$ is the neutral–ion collision frequency and L is the time-dependent current sheet half-length. This strongly couples neutrals to the reconnection outflow, while kinetic simulations evolve to allow $ν_{ni}L/v^*_A$ <1, decoupling neutrals from the reconnection outflow. Differences in the way reconnection is triggered may explain these discrepancies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Effects of midplane density gradient in the propagation of high-harmonic fast waves considering high temperature anisotropy in the scrape-off layer of NSTX-U

High-Harmonic Fast Wave (HHFW) heating experiments in NSTX have shown that up to 60% of the injected power can be lost in the Scrape-Off Layer (SOL) when the density is above the fast wave cutoff density in front of the antenna, for which the fast wave is able to propagate into the plasma. This work models HHFW propagation in the SOL plasmas of NSTX-U using a 2D divertor SOL profile derived from the pressure assumption and the finite element temperature solution, which accounts for the high anisotropy of heat conduction in a 2D axisymmetric geometry. In this work, the two-dimensional axisymmetric SOL temperature profile is first evaluated by solving the steady-state non-linear heat conduction equation, in which thermal conductivity depends on temperature, using a finite element approach in the Petra-M workbench. A 2D density profile is then obtained from a prescribed density profile at the outer midplane assuming pressure is uniform along a flux tube in the SOL. This approach results in density and temperature profiles in which the strong asymmetric nature of heat conduction between the parallel and perpendicular background magnetic field is successfully captured. Furthermore, this work focuses on investigating the effect of the SOL plasma density profile on parasitic HHFW propagation in the SOL. The simulation results show that the radial gradient of the density profile affects the wavefield propagation in the SOL. As the density profile broadens, the wavefield intensity and its poloidal extent are reduced in the SOL, and the core coupling increases. Similarly, as shown by collisional power deposition, a proxy for power absorption, the fraction of power deposited in the SOL decreases with a broadening profile.

HHFW↗

DOE OSTI2023

Hybrid-VPIC: An open-source kinetic/fluid hybrid particle-in-cell code

Hybrid-VPIC is an extension of the open-source high-performance particle-in-cell (PIC) code VPIC incorporating hybrid kinetic ion/fluid electron solvers. This paper describes the models that are available in the code and gives an overview of applications of the code to space and laboratory plasma physics problems. Particular choices in how the hybrid solvers were implemented are documented for reference by users. A few solutions for handling numerical complications particular to hybrid codes are also described. Finally, special emphasis is given to the computationally taxing problem of modeling mix in collisional high-energy-density regimes, for which more accurate electron fluid transport coefficients have been implemented for the first time in a hybrid PIC code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1984

Electron precipitation in solar flares - Collisionless effects

A large fraction of the electrons which are accelerated during the impulsive phase of solar flares stream towards the chromosphere and are unstable to the growth of plasma waves. The linear and nonlinear evolution of plasma waves as a function of time is analyzed with a set of rate equations that follows, in time, the nonlinearly coupled system of plasma waves-ion fluctuations. As an outcome of the fast transfer of wave energy from the beam to the ambient plasma, nonthermal electron tails are formed which can stabilize the anomalous Doppler resonance instability responsible for the pitch angle scattering of the beam electrons. The non-collisional losses of the precipitating electrons are estimated, and the observational implication of these results are discussed.

Vlahos, L.↗

NASA NTRS1998

Modeling the Plasmasphere

The plasmasphere has often been considered one of the more boring regions in the magnetosphere. Its low energy plasma doesn't begin to compete against the free sources of energy available in the ring current, auroral zone, or plasma sheet. Its best known feature is its relatively highly density, archived as a result of prolonged accumulation of ionospheric outflow onto corotating flux tubes. On second look, however, the plasmasphere can be found to exhibit a remarkable influence on its more energetic cousins and display convection behavior indicative of physical processes acting throughout the magnetosphere for which we have no explanation. Plasmaspheric plasma densities and composition of heavy ions are particularly sensitive to heating by processes active in the ionosphere and all along field lines. Wave propagation and instabilities, collisional losses in the ring current, and heat transport from superthermal electrons are all equally sensitive to dense, heavy plasmaspheric densities and density gradients. It is in this context that we seek to characterize plasmaspheric populations using event based, empirical, and physical modeling methods. The modeling approaches, the challenges, and some of the results of these efforts will be presented.

Gallagher, Dennis L↗

NASA NTRS2012

Three-Dimensional Numerical Simulations of Equatorial Spread F: Results and Observations in the Pacific Sector

A three-dimensional numerical simulation of plasma density irregularities in the postsunset equatorial F region ionosphere leading to equatorial spread F (ESF) is described. The simulation evolves under realistic background conditions including bottomside plasma shear flow and vertical current. It also incorporates C/NOFS satellite data which partially specify the forcing. A combination of generalized Rayleigh-Taylor instability (GRT) and collisional shear instability (CSI) produces growing waveforms with key features that agree with C/NOFS satellite and ALTAIR radar observations in the Pacific sector, including features such as gross morphology and rates of development. The transient response of CSI is consistent with the observation of bottomside waves with wavelengths close to 30 km, whereas the steady state behavior of the combined instability can account for the 100+ km wavelength waves that predominate in the F region.

ALTAIR↗

DOE OSTI2021

W7-X and the sawtooth instability: towards realistic simulations of current-driven magnetic reconnection

Magnetic reconnection in W7-X is studied by means of global numerical simulations in a series of models of increasing complexity. The magnetic geometry ranges from that of a cylinder to the full three-dimensional field of W7-X, and the equations solved range from ideal magnetohydrodynamics (MHD) to gyrokinetics. We simulate plasmas from the first operation phase with electron cyclotron current drive (ECCD). These are characterized by an equilibrium magnetic field featuring an ECCD-distorted 'humped' profile of the rotational transform ι, with ι = 1 in two radial locations. Such plasmas generally show sawtooth activity, hence motivating the present study. We pay particular attention to the role of equilibrium current density gradients in the destabilization of reconnecting modes. When the equilibrium temperature and density gradients are artificially suppressed (to eliminate the pressure gradient drive), the perturbed electrostatic potential is radially localized between the locations at which ι = 1. This is shown with a purely collisionless gyrokinetic model, in cylindrical geometry. In the real toroidal geometry of W7-X, for a non-ideal MHD model including a uniform resistivity, electron inertia and (numerical) viscosity, the same qualitative behaviour is observed. In particular, even if a resonant (m, n) = (1, -1) perturbation is initialized, the most unstable mode is the (m, n) = (-4, 4), where m and n are the poloidal and toroidal mode numbers, respectively. Other modes are destabilized due to geometric coupling. The growth rate of this instability scales as η 1/3 , where η is the plasma resistivity, thus suggesting that ECCD drives ideal MHD stable W7-X plasmas towards non-ideal marginality. An ideal magnetohydrodynamic analysis confirms the result. A fluid-kinetic hybrid version of the EUTERPE code shows that gyrokinetic ions have a stabilizing effect on these modes. For W7-X relevant collisionalities, the growth rate scales linearly with the electron skin depth, d e . Implications of our results for sawtoothing W7-X operation are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1983

Ionosphere

The original interest in an ionosphere on Jupiter was generated by the discovery of strong radio-frequency emissions at approximately 20 MHz which were thought to be plasma frequencies associated with Jupiter's ionosphere. The ionosphere of Jupiter provides a means to couple the magnetosphere to the atmosphere by virtue of its high conductivity and collisional interaction with the neutral atmosphere. The Pioneer and Voyager have provided direct measurements of profiles of electron concentration at selected locations on Jupiter. Attention is given to basic principles regarding the characteristics of the Jovian ionosphere, the ionization sources, aspects of ion recombination, ion chemistry, observations of Jupiter's ionosphere, the structure of Jupiter's upper atmosphere, and questions of ionospheric modeling. On the basis of the Pioneer and Voyager observations it appears that Jupiter's ionosphere and thermosphere undergo significant solar cycle changes.

Strobel, D. F.↗

DOE OSTI2022

Cross-beam energy transfer saturation: ion heating and pump depletion

Cross-beam energy transfer (CBET) was measured in two regimes where the energy transfer saturation mechanism was determined by the plasma and laser beam conditions. Here, linear kinetic CBET theory agrees well with the measured energy transfer in all experimental configurations and at all probe beam intensities when accounting for pump depletion and the plasma conditions measured using Thomson-scattering. Simultaneous CBET and Thomson-scattering measurements enabled uncertainties in the plasma conditions to be isolated from CBET theory, which allowed the saturation mechanisms to be identified. In the perpendicular-beam configuration the saturation mode was through ion heating, which resulted from ion trapping in the driven waves and subsequent ion-ion collisional heating. In the co-propagating beam configuration there was minimal ion heating and the saturation mode was through pump depletion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS2012

Electric Field and Plasma Density Observations of Irregularities and Plasma Instabilities in the Low Latitude Ionosphere Gathered by the C/NOFS Satellite

The Vector Electric Field Investigation (VEFI) on the C/NOFS equatorial satellite provides a unique data set which includes detailed measurements of irregularities associated with the equatorial ionosphere and in particular with spread-F depletions. We present vector AC electric field observations gathered on C/NOFS that address a variety of key questions regarding how plasma irregularities, from meter to kilometer scales, are created and evolve. The talk focuses on occasions where the ionosphere F-peak has been elevated above the C/NOFS satellite perigee of 400 km as solar activity has increased. In particular, during the equinox periods of 2011, the satellite consistently journeyed below the F-peak whenever the orbit was in the region of the South Atlantic anomaly after sunset. During these passes, data from the electric field and plasma density probes on the satellite have revealed two types of instabilities which had not previously been observed in the C/NOFS data set: The first is evidence for 400-500km-scale bottomside "undulations" that appear in the density and electric field data. In one case, these large scale waves are associated with a strong shear in the zonal E x B flow, as evidenced by variations in the meridional (outward) electric fields observed above and below the F-peak. These undulations are devoid of smaller scale structures in the early evening, yet appear at later local times along the same orbit associated with fully-developed spread-F with smaller scale structures. This suggests that they may be precursor waves for spread-F, driven by a collisional shear instability, following ideas advanced previously by researchers using data from the Jicamarca radar. A second result is the appearance of km-scale irregularities that are a common feature in the electric field and plasma density data that also appear when the satellite is near or below the F-peak at night. The vector electric field instrument on C/NOFS clearly shows that the electric field component of these waves is strongest in the zonal direction. These waves are strongly correlated with simultaneous observations of plasma density oscillations and appear both with, and without, evidence of larger-scale spread-F depletions. These km-scale, quasi-coherent waves strongly resemble the bottomside, sinusoidal irregularities reported in the Atmosphere Explorer satellite data set by Valladares et al. and are believed to cause scintillations of VHF radiowaves. We interpret these new observations in terms of fundamental plasma instabilities associated with the unstable, nighttime equatorial ionosphere.

Pfaff, Robert F.↗

DOE OSTI2021

An experimental characterization of core turbulence regimes in Wendelstein 7-X

First results from the optimized helias Wendelstein 7-X stellarator (W7-X) have shown that core transport is no longer mostly neoclassical, as is the case in previous kinds of stellarators. Instead, power balance analysis has shown that turbulent transport poses a serious limitation to the global performance of the machine. Furthermore, several studies have found this particularly relevant for ion transport, with core ion temperatures becoming clamped at relatively low values of T i ≃ 1.7 keV, except in the few scenarios in which turbulence can be suppressed. In order to understand the precise turbulent mechanisms at play and thus design improved performance scenarios, it is important to have a clear understanding of the parametric dependencies of turbulent fluctuations, and the relation between them and turbulent transport. As a first step in this direction, in this work we use Doppler reflectometry measurements carried out during a number of relevant operational scenarios to provide a systematic characterization of ion-scale ( k ⊥ ρ i ≃ 1) density fluctuations in the core of W7-X. Then, we study the relation between fluctuation amplitude and plasma profiles and show how distinct regimes can be defined for the former, depending on normalized gradients and . Furthermore, we discuss the importance of other potentially relevant parameters such as T e / T i , E r or collisionality. Comparing the different regimes, we find that turbulence amplitude depends generally on the gradient ratio η i = , as would be expected for ITG modes, with the exception of a range of discharges, for which turbulence suppression may be better explained by an ITG to TEM transition triggered by a drop in collisionality. Finally, we show a number of scenarios under which T i,core > 1.7 keV is achieved and how core fluctuations are suppressed in all of them, thus providing experimental evidence of microturbulence being the main responsible for the limited ion confinement in W7-X.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Verification of a Monte Carlo binary collision model for simulating elastic and inelastic collisions in particle-in-cell simulations

We present the development and verification of a Monte Carlo binary collision model for simulating elastic and inelastic collisions in particle-in-cell simulations. We apply the corrected binary collision model originally developed for charged-particles collisions to all considered scattering channels, including Coulomb collisions, elastic neutral–neutral and charged–neutral collisions, ionization, excitation, and fusion. The model's implementation is described and verified through a series of simulations, including charged- and neutral-particle thermal equilibration, slowing of electrons in warm solid-density aluminum, collisional damping of a Langmuir wave, helium gas breakdown in an applied electric field, and thermonuclear and beam–target fusion. Then, we demonstrate the model within simulations of hydrogen plasma formation in the Princeton Field-Reversed Configuration as well as of the burning of aneutronic fusion fuel p-11B. The latter includes measurement of the fusion power density in a low-density plasma and fusion production due to the stopping of a proton ignitor beam in a compressed boron target.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Close-Out Project Description for Koepke's Dept of Energy grant DE-SC0021405

Objectives: To establish, for lab & space conditions, EM-IEDDI’s (electromagnetic shear-driven instability's) dispersion relation, unstable range, instability threshold, and mode characteristics, we need LAPD’s Alfven-wave-favorable electromagnetic-style conditions, including higher "beta" (0.001 < beta ≤ 0.3) and low-collisionality. Also, we attempted to intentionally launch or spontaneously destabilize compressional and shear Alfven waves in the strong, localized, perpendicular-velocity-shear region at the interface between coaxial plasmas (one plasma cylinder inside an outer, otherwise hollow, tube, each having a different, controllable, value of plasma electrostatic potential, i.e., “space” potential (not to be confused with the temperature-dependent “floating” potential of an object immersed in the plasma). Nonlinear wave-wave interactions between same-family (EM-IEDD or Alfven) and cross-family (EM-IEDD-with-Alfven) fluctuations were targeted for documentation over a range of spectral overlap. Although laboratory experiments were conducted, the following theoretical work was left unfinished: Analytical non-modal prediction Computational non-modal prediction Check to see if Mikhailenko’s theory formulation leads to his published graphs

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗