Engineering Papers⌕ Search

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.

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

At least 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Breaking Beams: A Kinetic Treatment of the Filamentation and Forward Scatter of Laser Beams in Plasma

Precise and predictable propagation and deposition of laser energy in high energy density (HED) plasma is threatened by the breakup of the beam into filaments. Such filaments can be deflected and sprayed but also locally intensify the light, increasing the likelihood of other problematic interactions with the plasma. Adequately modeling, anticipating, and avoiding such processes in HED experiments requires a thorough understanding of the thresholds and nonlinear saturation mechanisms, particularly for the weakly collisional, high temperature plasmas achieved at the National Ignition Facility (NIF).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

The constraint of plasma power balance on runaway avoidance

Abstract In a post-thermal-quench plasma, mitigated or unmitigated, the plasma power balance is mostly between collisional or Ohmic heating and plasma radiative cooling. In a plasma of atomic mixture { n α } with α labeling the atomic species, the power balance sets the plasma temperature, ion charge state distribution { n α i } with i the charge number, and through the electron temperature T e and ion charge state distribution { n α i } , the parallel electric field E ∥ . Since the threshold electric field for runaway avalanche growth E a v is also set by the atomic mixture, ion charge state distribution and its derived quantity, the electron density n e , the plasma power balance between Ohmic heating and radiative cooling imposes a stringent constraint on the plasma regime for avoiding and minimizing runaways when a fusion-grade tokamak plasma is rapidly terminated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Numerical minimization of neoclassical poloidal viscosity for supersonic equilibria in tokamak geometry

Extensive experimental evidence has shown that the presence of poloidal flow in tokamaks can dramatically improve transport properties. However, theory indicates that poloidal flows are damped by poloidal viscosity, thus necessitating external drivers, such as neutral beam injection or radio frequency heating. In this work, ideal magnetohydrodynamic equilibria are calculated via the FORTRAN code FLOW [Guazzotto et al., Phys. Plasmas 11, 604 (2004)] and a postprocessor is used to estimate the neoclassical poloidal viscosity. The equilibrium inputs, which correspond to intuitive physical quantities, are then numerically optimized to reduce a viscosity figure of merit. We present supersonic equilibria in tokamak geometry with minimized neoclassical poloidal viscosities for various velocity free function inputs, plasma aspect ratios, and collisionality regimes. Benchmarks are made against an analytic theory as well as a classical expression of poloidal viscosity. Numerical confirmation of the analytic theory is obtained in the high aspect ratio and high collisionality limit. Good agreement is also seen near the plasma core and edge, with discrepancies arising in the intermediate region. Outside of these limits, rotation input function profiles are found that provide ∼order of magnitude improvements over the analytic theory, with additional progress being made toward predictions for tokamak-relevant equilibria.

Physics↗
NASA NTRS1974

The collisional drift mode in a partially ionized plasma

The structure of the drift instability was examined in several density regimes. Let sub e be the total electron mean free path, k sub z the wave-vector component along the magnetic field, and the ratio of perpendicular ion diffusion to parallel electron streaming rates. At low densities (k sub z lambda 1) the drift mode is isothermal and should be treated kineticly. In the finite heat conduction regime square root of m/M k sub z Lambda sub 1) the drift instability threshold is reduced at low densities and increased at high densities as compared to the isothermal threshold. Finally, in the energy transfer limit (k sub z kambda sub e square root of m/M) the drift instability behaves adiabatically in a fully ionized plasma and isothermally in a partially ionized plasma for an ion-neutral to Coulomb collision frequency ratio.

Hudson, M. K.↗
DOE OSTI2023

Density wavenumber spectrum measurements, synthetic diagnostic development, and tests of quasilinear turbulence modeling in the core of electron-heated DIII-D H-mode plasmas

Abstract Measurements of the turbulent density wavenumber spectrum, δ n ˆ e ( k ⊥ ) , using the Doppler Back-Scattering (DBS) diagnostic are reported from DIII-D H-mode plasmas with electron cyclotron heating as the only auxiliary heating method. These electron-heated plasmas have low collisionality, ν e ∗ < 1 , T e / T i > 1 , and zero injected torque—a regime expected to be relevant for future fusion devices. We probe density fluctuations in the core ( ρ ≈ 0.7) over a broad wavenumber range, 0.5 ⩽ k ⊥ ⩽ 16 cm −1 ( 0.1 ⩽ k ⊥ ρ s ⩽ 5 ), to characterize plasma instabilities and compare with theoretical predictions. We present a novel synthetic DBS diagnostic to relate the back-scattered power spectrum, P s ( k ⊥ ) —which is directly measured by DBS—to the underlying electron density fluctuation spectrum, δ n ˆ e ( k ⊥ ) . The synthetic DBS P s ( k ⊥ ) spectrum is calculated by combining the SCOTTY beam-tracing code with a model δ n ˆ e ( k ⊥ ) predicted either analytically or numerically. In this work we use the quasi-linear code Trapped Gyro-Landau Fluid (TGLF) to approximate the δ n ˆ e ( k ⊥ ) spectrum. We find that TGLF, using the experimental profiles, is capable of closely reproducing the DBS measurements. Both the DBS measurements and the TGLF-DBS synthetic diagnostic show a wavenumber spectrum with variable decay. The measurements show weak decay ( k −0.6 ) for k < 3.5 cm −1 , with k −2.6 at intermediate- k ( 3.5 ⩽ k ⩽ 8.5 cm −1 ), and rapid decay ( k −9.4 ) for k > 8.5 cm −1 . Scans of physics parameters using TGLF suggest that the normalized ∇ T e scale-length, R / L T e , is an important factor for distinguishing microturbulence regimes in these plasmas. A combination of DBS observations and TGLF simulations indicate that fluctuations remain peaked at ITG-scales (low k ) while R / L T e -driven TEM/ETG-type modes (intermediate/high k ) are marginally sub-dominant.

synthetic diagnostics↗
NASA NTRS2005

Probing the Reverse Shock in an Oxygen-Rich Supernovae Remnant

The purpose of this investigation is to examine the O VI emission from the X-ray bright ring of the supernova remnant 1 E0102.2-729 in the small Magellanic cloud. Three pointings were positioned tangent to the ring, north (N), northeast (NE), and southeast (SE), to examine a range of X-ray emitting regions overlapping a range of optical [O III] nebulosity and to examine the velocity structure. One background pointing was also obtained, but it was contaminated by a star. The background levels in the pointings on the remnant were low enough that the background pointing was not required for the remaining analysis. The SE pointing was reobserved in August, 2004, in order to bring the total exposure up to the originally requested 15 ks. The archive notified us of the data's availability in mid September. Significant broad O VI 1032 and O VI 1038 emission was found, brightest in the NE and SE pointings. In the NE and SE pointings, the FWHM of the broad O VI component is approx.800-1000 km/s, while in the N pointing, the line is approx.1500 km/s wide. The O VI is redshifted in the N (approx.380 km/s) and NE (approx.60 km/s) but is blueshifted in the SE (approx. -160 km/s). These FUSE O VI velocity dispersions can be compared to the X-ray gas velocities inferred from Doppler distortions in the Chandra X-ray data as reported by Flanagan et al. 2004 (ApJ 605, 230). The bulk velocities in the X-ray bright ring of order +/- 1000 km/s, comparable to the velocity dispersion seen in the FUSE data. However, the X-ray data indicates a redshift of approx.1000 km/s in the SE, while the FUSE data show a blueshift of approx.160 km/s, underscoring the complex velocity structure in this remnant. The O VI fluxes estimated from the fits to the FUSE data were combined with X-ray (XMM- Newton) O VI1 and O VIII fluxes and compared with predictions from a plasma nonequilibrium ionization model in a "line-based" analysis. We found that the plasma departs significantly from collisional ionization equilibrium, particularly in the SE, and that the plasma excitation conditions vary among the pointings: the O plasma sees different conditions than the Ne and Mg plasmas.

Gaetz, Terrance↗
DOE OSTI2025

General kinetic ion-induced electron emission model for metallic walls applied to biased Z-pinch electrodes

A kinetic ion-induced electron emission (IIEE) model for general applications is developed to obtain the emitted electron energy spectrum for a distribution of ion impacts on a metallic surface. We assume an ionization cascade mechanism and use empirical models for the ion and electron stopping powers. The emission spectrum and the secondary electron yield (SEY) are validated for a variety of materials. The IIEE model is used to study the effect of IIEE on the plasma-material interactions of Z-pinch electrodes. Un-magnetized Boltzmann-Poisson simulations are performed for a Z-pinch plasma doubly bounded by two biased copper electrodes with and without IIEE at bias potentials from 0 to 9 kV. At the anode, the SEY decreases from 0 to 1 kV, but then increases at higher bias potentials. At the cathode, the SEY is much larger due to higher energy ion bombardment and grows with bias potential. As the bias potential increases, the emitted cathode electrons are accelerated to higher energies into the domain, collisionally heating the plasma. Above 1 kV, the heating is strong enough to increase the plasma potential. Despite SEY greater than 1, only a classical sheath forms as opposed to a space-charge limited or inverse sheath due to the emitted electron flux not reaching the space charge current saturation limits. Furthermore, the current in the emissionless cases saturates to a value lower than experiment. With IIEE, the current does not saturate and continues to increase with the 4 kV case, matching most closely with the experiment.

Carbon based materials↗
DOE OSTI2024

NSTX-U research advancing the physics of spherical tokamaks

Abstract The objectives of NSTX-U research are to reinforce the advantages of STs while addressing the challenges. To extend confinement physics of low- A , high beta plasmas to lower collisionality levels, understanding of the transport mechanisms that set confinement performance and pedestal profiles is being advanced through gyrokinetic simulations, reduced model development, and comparison to NSTX experiment, as well as improved simulation of RF heating. To develop stable non-inductive scenarios needed for steady-state operation, various performance-limiting modes of instability were studied, including MHD, tearing modes, and energetic particle instabilities. Predictive tools were developed, covering disruptions, runaway electrons, equilibrium reconstruction, and control tools. To develop power and particle handling techniques to optimize plasma exhaust in high performance scenarios, innovative lithium-based solutions are being developed to handle the very high heat flux levels that the increased heating power and compact geometry of NSTX-U will produce, and will be seen in future STs. Predictive capabilities accounting for plasma phenomena, like edge harmonic oscillations, ELMs, and blobs, are being tested and improved. In these ways, NSTX-U researchers are advancing the physics understanding of ST plasmas to maximize the benefit that will be gained from further NSTX-U experiments and to increase confidence in projections to future devices.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS1987

MPD thruster plasma instability studies

Results of an ongoing experimental/theoretical effort to understand the influence of plasma waves and instabilities upon coaxial plasma thruster performance and longevity are discussed. An experiment which permits measurement of the plasma dispersion relation plus a suite of plasma parameters for each one millisecond quasi-steady discharge is described. This investigation is guided by a linearized Boltzmann-Poisson kinetic numerical model which can accommodate the effects of collisions, inhomogeneities and the magnetic field. Preliminary tests reveal the presence of the low-frequency branch of an ion cyclotron wave. Although the stability of these waves still has to be established through spatial growth rate measurements, it is well known that they are prone to 'current-driven' instabilities. 'Current-driven' instabilities in the collisional tensor conducting plasma of the MPD thruster are analytically discussed and their relevance to the MPD plasma is illustrated by estimating a critical total current at which their onset occurs.

Kelly, A. J.↗
DOE OSTI2023

Effects of ion trapping and fluctuations of electron temperature and plasma flow on cross-beam energy transfer

The influences of ion trapping and fluctuations of electron temperature and plasma flow on cross-beam energy transfer (CBET) are examined using two- and three-dimensional particle-in-cell simulations in parameter regimes relevant to recent CBET experiments at the OMEGA laser facility. In mid-Z plasma irradiated by an intense pump beam and weaker probe beam, ion trapping, collisional de-trapping, and plasma flow induced by thermal effects are shown to affect the CBET gain. Ion trapping can enhance or detune the CBET resonance [Nguyen et al., Phys. Plasmas 28, 082705 (2021)]. Collisional de-trapping can affect the CBET gain at low seed beam intensity near the onset threshold for ion trapping. Thermal-effects-induced flow can also detune the CBET resonance at a level comparable to that from trapping at low seed beam intensity. As a consequence, the CBET gain is sensitive to collisions and dimensionality at low seed beam intensity where ion trapping is weak but is insensitive to collisions and dimensionality at high seed beam intensity where ion trapping is strong.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1994

A time-dependent gyro-kinetic model of thermal ion upflows in the high-latitude F region

Ample evidence supports the significance of the high-latitude ionospheric contribution to magnetospheric plasma. Assuming flux conservation along a flux tube, the upward field-aligned ion flows observed in the magnetosphere require high-latitude ionospheric field-aligned ion upflows of the order of 10(exp 8) to 10(exp 9)/sq cm/s. Since radar and satellite observations of high-latitude F region flows at times exceed this flux requirement by an order of magnitude, the thermal ionospheric upflows are not simply the ionospheric response to a magnetospheric flux requirement. Several ionospheric ion upflow mechanisms have been proposed, but simulations based on fluid theory do not reproduce all the observed features of ionospheric ion upflows. Certain asymmetries in the statistical morphology of high-latitude F region ion upflows suggest that the ion upflows may be generated by ion-neutral frictional heating. We developed a single-component (O(+)), time-dependent gyro-kinetic model of the high-latitude F region response to frictional heating in which the neutral exobase is a discontinuous boundary between fully collisional and collisionless plasmas. The concept of a discontinuous neutreal exobase and the assumption of a constant and uniform polarization electric field reduce the ion velocity distribution function, from which we can compute the ion density, parallel velocity, parallel and perpendicular temperature, and parallel flux. Using our model, we simulated the response of a convecting flux tube between 500 km and 2500 km to various frictional heating inputs; the results were both qualitatively and quantitatively different from fluid model results, which may indicate an inadequacy of the fluid theory approach. The gyro-kinetic frictional heating model responses to the various simulations were qualitatively similar: (1) initial perturbations of all the modeled parameters propagated rapidly up the flux tube, (2) transient values of the ion parallel velocity, temperature, and flux exceeded 3 km/s, 2 x 10(exp 4) K, and 10(exp 9)/sq cm/s, respectively, (3) a second transient regime developed wherein the parallel temperature drops to very low values (a few hundred Kelvins), and (4) well after heating ceased, large parallel temperatures and large downward parallel velocities and fluxes developed as the flux tube slowly returned to diffusive equilibrium. The ion velocity distributions during the simulation are often non-Maxwellian and are sometimes composed of two distinct ion populations.

Loranc, M.↗
DOE OSTI2021

Total X-Ray Emission from the LMC Observed with HaloSat

We constrain the global parameters of the Large Magellanic Cloud (LMC) by modeling the first soft X-ray (0.4–6.8 keV) observations of the entire LMC produced in a single pointing with moderate spectral resolution by the HaloSat CubeSat. These data are best fit with the sum of two thermal plasma components in collisional ionization equilibrium and a power law. We find cool (0.210 ± 0.014 keV) and hot (0.89 ± 0.14 keV) components. The total X-ray luminosity of the LMC is (1.08 ± 0.14) × 10{sup 39} erg s{sup −1}. X-ray binaries make up a large fraction of the emission with a luminosity of (6.0 ± 0.8) × 10{sup 38} erg s{sup −1}, followed by cool gas from superbubbles, supernovae, and diffuse emission with a luminosity of (3.0 ± 0.3) × 10{sup 38} erg s{sup −1}. The hot gas from star formation contributes the smallest fraction, with a luminosity of (1.9 ± 0.5) × 10{sup 38} erg s{sup −1}. We estimate the total volume of the cool gas to be between (0.2–1.2) × 10{sup 10} pc{sup 3} and the hot gas to be between (1.0–5) × 10{sup 7} pc{sup 3} for filling factors of f = 1 and 0.2. These volumes result in a total thermal energy for the cool gas between (1.4–3) × 10{sup 54} erg for electron densities of 0.017–0.04 cm{sup −3}, and a thermal energy for the hot gas between (1.7–4) × 10{sup 53} erg for electron densities of 0.12–0.3 cm{sup −3}. This yields cooling timescales for the cool and hot gas of (1.5–3) × 10{sup 8} yr and (1.8–6) × 10{sup 7} yr, respectively.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Pressure–strain interaction. III. Particle-in-cell simulations of magnetic reconnection

How energy is converted into thermal energy in weakly collisional and collisionless plasma processes, such as magnetic reconnection and plasma turbulence, has recently been the subject of intense scrutiny. The pressure–strain interaction has emerged as an important piece, as it describes the rate of conversion between bulk flow and thermal energy density. In two companion studies, we presented an alternate decomposition of the pressure–strain interaction to isolate the effects of converging/diverging flow and flow shear instead of compressible and incompressible flow, and we derived the pressure–strain interaction in magnetic field-aligned coordinates. Here, we use these results to study pressure–strain interaction during two-dimensional anti-parallel magnetic reconnection. We perform particle-in-cell simulations and plot the decompositions in both Cartesian and magnetic field-aligned coordinates. We identify the mechanisms contributing to positive and negative pressure–strain interaction during reconnection. Furthermore, this study provides a roadmap for interpreting numerical and observational data of the pressure–strain interaction, which should be important for studies of reconnection, turbulence, and collisionless shocks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1972

Rocket spectroscopy of zeta Orionis.

Analysis of a spectrum of zeta Ori extending from 922 to 1453 A with approximately 0.8 A resolution recorded at rocket altitudes. All lines used in existing models of stellar atmospheres appear in the recorded spectrum with the exception of those masked by telluric N2 or strong P Cygni-type profiles and by an O V line at 1371.29 A. Fifteen multiplets of subordinate lines have been reliably identified, indicating an approximate range of excitation from 0 to 50 eV. Transitions in C III (1176 A), N III (991 A), N V (1239, 1243 A), O VI (1032, 1038 A), Si IV (1394, 1403 A), S IV (1063, 1074 A), and S VI (933, 944 A) have been observed as P Cygni-type profiles presumably arising in a circumstellar envelope. The degree of ionization, transitions present, and mean radial velocities are all consistent with viewing the envelope as a hot (about 100,000 K), rarefied plasma in which collisional ionization is important. Interstellar lines in C I (1277, 1280 A), C II (1036, 1334 A), N I (1134-1135 A), N I (1200-1201 A), N II (1084-1086 A), O I (1302, 1305 A), Si II (1190, 1193 A), Si II (1260 A), and Si II (1304 A) have been definitely identified. Other transitions in Ar II, S I, C I, and Fe II are tentatively identified. The equivalent width of the L alpha line is found to be 10.4 plus or minus 1.6 A, corresponding to a columnar density of 2.0 plus or minus 0.7 x 10 to the 20th per cu cm.

Smith, A. M.↗
NASA NTRS1977

Current-limited electron beam injection

The injection of an electron beam into a weakly collisional, magnetized background plasma was investigated experimentally. The injected beam was energetic and cold, the background plasma was initially isothermal. Beam and plasma dimensions were so large that the system was considered unbounded. The temporal and spatial evolution of the beam-plasma system was dominated by collective effects. High-frequency electrostatic instabilities rapidly thermalized the beam and heated the background electrons. The injected beam current was balanced by a return current consisting of background electrons drifting toward the beam source. The drift between electrons and ions gave rise to an ion acoustic instability which developed into strong three-dimensional turbulence. It was shown that the injected beam current was limited by the return current which is approximately given by the electron saturation current. Non-Maxwellian electron distribution functions were observed.

Stenzel, R. L.↗
NASA NTRS1978

Experiments on current-driven three-dimensional ion sound turbulence. I - Return-current limited electron beam injection. II - Wave dynamics

Pulsed electron beam injection into a weakly collisional magnetized background plasma is investigated experimentally; properties of the electron beam and background plasma, as well as the low-frequency instabilities and wave dynamics, are discussed. The current of the injected beam closes via a field-aligned return current of background electrons. Through study of the frequency and wavenumber distribution, together with the electron distribution function, the low-frequency instabilities associated with the pulsed injection are identified as ion acoustic waves driven unstable by the return current. The frequency cut-off of the instabilities predicted from renormalized plasma turbulence theory, has been verified experimentally.

Stenzel, R. L.↗
NASA NTRS1981

Heavy ions in Jupiter's environment

The extended atmosphere of the Jupiter system consists of heavy-element atoms and ions. This material originates on the satellite Io. Energy is lost from the thermal plasma in collisionally-excited optical and ultraviolet emission. The juxtaposition of earth and spacecraft measurements provide insight concerning the underlying processes of particle transport and energy supply.

Brown, R. A.↗
NASA NTRS1984

Nonlinear mode coupling theory of the lower-hybrid-drift instability

A nonlinear mode coupling theory of the lower-hybrid-drift instability is presented. A two-dimensional nonlinear wave equation is derived which describes lower-hybrid drift wave turbulence in the plane transverse to B (k.B = 0), and which is valid for finite beta, collisional and collisionless plasmas. The instability saturates by transferring energy from growing, long wavelength modes to damped, short wavelength modes. Detailed numerical results are presented which compare favorably to both recent computer simulations and experimental observations. Applications of this theory to space plasmas, the earth's magnetotail and the equatorial F region ionosphere, are discussed. Previously announced in STAR as N84-17734

Drake, J. F.↗