Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Kinetic plasmas”

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 307 records · Page 17

Kinetic Modeling Analysis of Ar Addition to Atmospheric Pressure N 2 –H 2 Plasma for Plasma-Assisted Catalytic Synthesis of NH 3

Zero-dimensional kinetic modeling of atmospheric pressure Ar–N 2 –H 2 nonthermal plasma was carried out to gain mechanistic insights into plasma-assisted catalytic synthesis of ammonia. Ar dilution is a common technique for tailoring plasma discharge properties and has been shown to enhance NH 3 formation when added to N 2 –H 2 plasma. The kinetic model was developed for a coaxial dielectric barrier discharge quartz wool-packed bed reactor operating at near room temperature using a kHz-frequency plasma source. With 30% Ar mixed in a 1:1 N 2 –H 2 plasma at 760 Torr, we find that NH 3 production is dominated by Eley–Rideal (E-R) surface reactions, which heavily involve surface NH x species derived from N and H radicals in the gas phase, while the influence of excited N 2 molecules is negligible. This is contrary to the commonly proposed mechanism that excited N 2 molecules created by Penning excitation of N 2 by Ar(4s) and Ar(4p) play a significant role in assisting NH 3 formation. Our model shows that the enhanced NH 3 formation upon Ar dilution is unlikely due to the interactions between Ar and H species, as excited Ar atoms have a weak effect on H radical formation through H 2 dissociation compared to electrons. We find that excited Ar atoms contribute to 28% of the N radical production in the gas phase via N 2 dissociation, while the rest are dominated by electron-impact dissociation. Furthermore, Ar species play a negligible role in the product NH 3 dissociation. N 2 conversion sensitivity analyses were carried out for electron number density (ne) and reduced electric field (E/N), and contributions from Ar to gas-phase N radical production were quantified. Finally, the model can provide guidance on potential reasons for observing enhanced NH 3 formation upon Ar dilution in N 2 –H 2 plasma beyond changes in the discharge characteristics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Axisymmetric gyrokinetic simulation of ASDEX-Upgrade scrape-off layer using a conservative implicit BGK collision operator

Collisions play an important role in turbulence and transport of fusion plasmas. For kinetic simulations, as the collisionality increases in the domain of interest, the size of the time step to resolve the collisional physics can become overly restrictive in an explicit time integration scheme, leading to high computational cost. With the aim of overcoming such restriction, we have implemented an implicit Bhatnagar–Gross–Krook (BGK) collision operator for use in the discontinuous Galerkin full-f gyrokinetic solver within the Gkeyll framework, which, when combined with Gkeyll's traditional explicit time integrator for collisionless advection, can significantly increase the time step in gyrokinetic simulations of highly collisional regimes. To ensure conservation of density, momentum, and energy, we utilize an iterative scheme to correct the discretized approximation to the equilibrium Maxwellian distribution to which the BGK collision operator relaxes. We have further generalized the BGK infrastructure, both the implicit scheme and the correction routine, to handle cross-species collisions. This improved implicit and conservative BGK operator is benchmarked against the more accurate but more computationally expensive Lenard–Bernstein–Dougherty (LBD) operator, which has been utilized in prior studies with Gkeyll. The implicit BGK operator enables 2D axisymmetric simulations of the ASDEX-Upgrade scrape-off layer to run 56 times faster to completion than the simulations with the LBD operator, because the BGK operator is more robust and converges at a lower resolution than is required by the LBD operator. Additionally, in this more collisional limit, we demonstrate that the results of our simulations utilizing the implicit BGK operator agreed well with simulations utilizing the more computationally expensive LBD operator.

Gyrokinetic simulations↗

Linear embedding of nonlinear dynamical systems and prospects for efficient quantum algorithms

The simulation of large nonlinear dynamical systems, including systems generated by discretization of hyperbolic partial differential equations, can be computationally demanding. Such systems are important in both fluid and kinetic computational plasma physics. This motivates exploring whether a future error-corrected quantum computer could perform these simulations more efficiently than any classical computer. In this work, we describe a method for mapping any finite nonlinear dynamical system to an infinite linear dynamical system (embedding) and detail three specific cases of this method that correspond to previously studied mappings. Then we explore an approach for approximating the resulting infinite linear system with finite linear systems (truncation). Using a number of qubits only logarithmic in the number of variables of the nonlinear system, a quantum computer could simulate truncated systems to approximate output quantities if the nonlinearity is sufficiently weak. Other aspects of the computational efficiency of the three detailed embedding strategies are also discussed.

97 MATHEMATICS AND COMPUTING↗

Interpreting ion-energy distributions using charge exchange emitted from deeply kinetic field-reversed-configuration plasmas

Energy resolving the atomic hydrogen outflux created within plasmas by charge exchange (CX) of hot plasma ions with cooler hydrogen neutrals is used to infer the ion energy distribution within the plasma, IED p . In high- β plasma with field nulls and ion gyro-radii comparable to the plasma size, the measured ion energy distribution (IED m ) of the CX outflux will depend on the viewing angle and position of the detector. Here, we describe the physics for this, results from a synthetic diagnostic code that contrasts the IED m to the IED p within relatively small and hot field-reversed-configuration plasmas, and how these data can show the presence of magnetic nulls and different orbit classes in the plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Global Theory of Microtearing Modes in the Tokamak Pedestal

The pedestal of H-mode tokamaks displays strong magnetic fluctuations correlated with the evolution of the electron temperature. The microtearing mode (MTM), a temperature-gradient-driven instability that alters magnetic topology, is compatible with these observations. Here we extend the conventional theory of the MTM to include the global variation of the temperature and density profiles. The offset between the rational surface and the location of the pressure gradient maximum (μ) emerges as a crucial parameter for MTM stability. In conclusion, the extended theory matches observations on the Joint European Torus tokamak.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Piecewise Omnigenous Stellarators

In omnigenous magnetic fields, charged particles are perfectly confined in the absence of collisions and turbulence. For this reason, the magnetic configuration is optimized to be close to omnigenity in any candidate for a stellarator fusion reactor. However, approaching omnigenity imposes severe constraints on the spatial variation of the magnetic field. In particular, the topology of the contours of constant magnetic field strength on each magnetic surface must be such that there are no particles transitioning between different types of wells. This, in turn, usually leads to complicated plasma shapes and coils. This Letter presents a new family of optimized fields that display tokamak-like collisional energy transport while having transitioning particles. This result radically broadens the space of accessible reactor-relevant configurations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A heat-driven monochromatic light source

The efficiency with which heat may be converted into resonance radiation in a cesium thermionic diode is investigated theoretically. An analytical model of a thermionic converter is used which combines the coupled effects of line radiation transport, excited-state kinetics, and plasma diffusion. Operating regimes are established for various degrees of optical density in the plasma. The results indicate that monochromatic radiation can be produced with efficiencies on the order of 30 percent, provided that there is an adequate voltage drop across the plasma. A drop of 1 V was used since it can be maintained without any electrical power input to the device. It is found that high efficiencies are ude to the higher interelectrode distances which the solutions accommodate, and that radiation can be generated efficiently, even with optically dense gases.

Stefani, Francis↗

Kinetic closures for unmagnetized and magnetized plasmas

Parallel and perpendicular closures with cyclotron resonance effects retained for the five-moment (density, temperature, and flow velocity) fluid equations are derived by solving the kinetic equation with the Bhatnagar–Gross–Krook operator in Fourier space. For parallel propagation, the parallel closures are reduced to those of Ji et al. [Phys. Plasmas 20, 082121 (2013)]. The closures when combined to the fluid equations reproduce the fully kinetic dispersion relation that can be directly derived from the kinetic equation. The closures for the five-moment fluid system can be utilized to derive closures for the extended fluid system, which is demonstrated by deriving closures for the ten-moment system consisting of density, flow velocity, temperature, and viscosity tensor equations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Momentum and Heat Flux Measurements in the Exhaust of VASIMR using Helium Propellant

Interplanetary travel requires propulsion systems that can provide high specific impulse (Isp), while also having sufficient thrust to rapidly accelerate large payloads. One such propulsion system is the Variable Specific Impulse Magneto-plasma Rocket (VASIMR), which creates, heats, and ejects plasma to provide variable thrust and Isp, designed to optimally meet the mission requirements. The fraction of the total energy invested in creating the plasma, as compared to the plasma's total kinetic energy, is an important factor in determining the overall system efficiency. In VASIMR, this 'frozen flow loss' is appreciable when at high thrust, but negligible at high Isp. The loss applies to other electric thrusters as well. If some of this energy could be recovered through recombination processes, and reinjected as neutral kinetic energy, the efficiency of VASIMR, in its low Isp/high thrust mode may be improved. An experiment is being conducted to investigate the possibility of recovering some of the energy used to create the plasma by studying the flow characteristics of the charged and neutral particles in the exhaust of the thruster. This paper will cover the measurements of momentum flux and heat flux in the exhaust of the VASIMR test facility using helium as the propellant where the heat flux is comprised of both kinetic and plasma recombination energy. The flux measurements also assist in diagnosing and verifying the plasma conditions in the existing experiment.

Chavers, D. Gregory↗

Comparison of excitation and kinetic temperatures in a laser-produced plasma using absorption spectroscopy

High-resolution tunable laser absorption spectroscopy is used to measure time-resolved absorption spectra for six neutral uranium transitions in a laser-produced plasma. Here analysis of the spectra shows that kinetic temperatures are similar for all six transitions, but excitation temperatures are higher than kinetic temperatures from 10–100 μs, indicating departures from local thermodynamic equilibrium.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Ground state rotational and kinetic temperatures of C 2 molecules in a laser-produced plasma

This study presents an experimental comparison between rotational and kinetic temperatures of molecular species in a laser-produced plasma, utilizing tunable laser absorption spectroscopy. We analyzed high-resolution and time-resolved absorption spectra of 31 rotational transitions within the Phillips A 1 II u –X 1 Σ g 3-0 band of the C 2 molecule, measured from a laser-produced graphite plasma in ~2 Torr argon, within the 771–781 nm spectral range. Spectral fitting provided quantitative state column densities of C 2 rotational energy levels in the ground electronic and vibrational states. Boltzmann analysis revealed thermalized rotational distributions showing rotational temperatures decreasing from 3000 to 600 K within 5–500 μs following the onset of the laser plasma. The measured kinetic and rotational temperatures of C 2 molecules are comparable, suggesting that collisional processes in the laser-produced plasma are sufficient to maintain thermal equilibrium between the C 2 molecules and the Ar background gas.

Absorption spectroscopy↗

Electron Density and Two-Channel Neutron Emission Measurements in Steady-State Spherical Inertial-Electrostatically Confined Plasmas, with Review of the 1-D Kinetic Model

Electron density measurements have been made in steady-state plasmas in a spherical inertial electrostatic confinement (IEC) discharge using microwave interferometry. Plasma cores interior to two cathodes, having diameters of 15 and 23 cm, respectively, were probed over a transverse range of 10 cm with a spatial resolution of about 1.4 cm for buffer gas pressures from 0.2 to 6 Pa in argon and deuterium. The transverse profiles are generally flat, in some cases with eccentric symmetric minima, and give mean densities of from approx. = 0.4 to 7x 10(exp 10)/cu cm, the density generally increasing with the neutral gas pressure. Numerical solutions of the 1-D Poisson equation for EC plasmas are reviewed and energy distribution functions are identified which give flat transverse profiles. These functions are used with the plasma approximation to obtain solutions which also give densities consistent with the measurements, and a double potential well solution is obtained which has minima qualitatively similar to those observed. Explicit consideration is given to the compatibility of the solutions interior and exterior to the cathode, and to grid transparency. Deuterium fusion neutron emission rates were also measured and found to be isotropic, to within the measurement error, over two simultaneous directions. Anisotropy was observed in residual emissions during operation with non-fusing hydrogen-1. The deuterium rates are consistent with predictions from the model.

Dobson, Chris C.↗

Excitation of whistler and slow-X waves by runaway electrons in a collisional plasma

Runaway electrons are known to provide robust ideal or collisionless kinetic drive for plasma wave instabilities in both the whistler and slow-X branches, via the anomalous Doppler-shifted cyclotron resonances. In a cold and dense post-thermal-quench plasma, collisional damping of the plasma waves can compete with the collisionless drive. Previous studies have found that, due to their higher wavelength and frequency, slow-X waves suffer stronger collisional damping than the whistlers, while the ideal growth rate of slow-X modes is higher. Here, we study runaway avalanche distributions that maintain the same eigen distribution and increase only in magnitude over time. The distributions are computed from the relativistic Fokker–Planck–Boltzmann solver, upon which a linear dispersion analysis is performed to search for the most unstable or least damped slow-X and whistler modes. Taking into account the effect of plasma density, plasma temperature, and effective charge number, we find that the slow-X modes tend to be excited before the whistlers in a runaway current ramp-up. Furthermore, even when the runaway current density is sufficiently high that both branches are excited, the most unstable slow-X mode has a much higher growth rate than the most unstable whistler mode. The qualitative and quantitative trends uncovered in the current study indicate that even though past experiments and modeling efforts have concentrated on whistler modes, there is a compelling case that slow-X modes should also be a key area of focus in the runaway self-mediation through wave instabilities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

High-fidelity kinetic modeling of instabilities and gyromotion physics in nonuniform low-beta plasmas

A fourth-order accurate continuum kinetic Vlasov solver and a systematic method for constructing customizable kinetic equilibria are demonstrated to be powerful tools for the study of nonuniform collisionless low-beta plasmas. The noise-free methodology is applied to investigate two gradient-driven instabilities in 4D (x,y,vx,vy) phase space: the Kelvin–Helmholtz instability and the lower hybrid drift instability. Nonuniform two-species configurations where ion gyroradii are comparable to gradient scale lengths are explored. The approach sheds light on the evolution of the pressure tensor in Kelvin–Helmholtz instabilities and demonstrates that the associated stress tensor deviates significantly from the gyroviscous stress tensor. Even at high magnetization, first-order approximations to finite-gyromotion physics are shown to be inadequate for the Kelvin–Helmholtz instability, as shear scales evolve to become on par with gyromotion scales. The methodology facilitates exploring transport and energy partitioning properties associated with lower hybrid drift instabilities in low-beta plasma configurations. Distribution function features are captured in detail, including the formation of local extrema in the vicinity of particle-wave resonances. Finally, the approach enables detailed targeted investigations and advances kinetic simulation capability for plasmas in which gyromotion plays an important role.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗