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 793 records · Page 44

Understanding the surface wave characteristics using 2D particle-in-cell simulation and deep neural network

Here, the characteristics of the surface waves along the interface between a plasma and a dielectric material have been investigated using kinetic particle-in-cell simulations. A microwave source of GHz frequency has been used to trigger the surface wave in the system. The outcome indicates that the surface wave gets excited along the interface of plasma and the dielectric tube and appears as light and dark patterns in the electric field profiles. The dependency of radiation pressure on the dielectric permittivity and supplied input frequency has been investigated. Further, we assessed the capabilities of neural networks to predict the radiation pressure for a given system. The proposed deep neural network model is aimed at developing accurate and efficient data-driven plasma surface wave devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Spatiotemporal evolution of emission and absorption signatures in a laser-produced plasma

Here we report spatiotemporal evolution of emission and absorption signatures of Al species in a nanosecond (ns) laser-produced plasma (LPP). The plasmas were generated from an Inconel target, which contained ~0:4 wt. % Al, using 1064 nm, ≈6 ns full width half maximum pulses from an Nd:YAG laser at an Ar cover gas pressure of ≈34 Torr. The temporal distributions of the Al I (394.4 nm) transition were collected from various spatial points within the plasma employing time-of-flight (TOF) emission and laser absorption spectroscopy, and they provide kinetics of the excited state and ground state population of the selected transition. The emission and absorption signatures showed multiple peaks in their temporal profiles, although they appeared at different spatial locations and times after the plasma onset. The absorption temporal profiles showed an early time signature representing shock wave propagation into the ambient gas. We also used emission and absorption spectral features for measuring various physical properties of the plasma. The absorption spectral profiles are utilized for measuring linewidths, column density, and kinetic temperature, while emission spectra were used to measure excitation temperature. A comparison between excitation and kinetic temperature was made at various spatial points in the plasma. Our results highlight that the TOF measurements provide a resourceful tool for showing the spatiotemporal LPP dynamics with higher spatial and temporal resolution than is possible with spectral measurements but are difficult to interpret without additional information on excitation temperatures and linewidths. The combination of absorption and emission TOF and spectral measurements thus provides a more complete picture of LPP spatiotemporal dynamics than is possible using any one technique alone.

(140.3440) Laser-induced breakdown, (280.4788) Opt↗

Plasma flow generation and particle acceleration from expanding magnetic bubbles

Impulsive plasma dynamics in the laboratory are often driven by rising electric currents, yet their quantitative plasma response has not been well established. By means of fully kinetic particle-in-cell simulations and laser-driven capacitor-coil experiments, we show that a rising current expels plasma, forming an expanding magnetic bubble and accelerating particles. The expansion front velocity scales with the Alfvén speed determined by the magnetic field at its inner edge and the plasma density at its outer edge. This mechanism establishes impulsive current drive as a fundamental way that generates plasma flows and accelerates particles in laboratory plasmas, with potential relevance to astrophysics.

Zhang, Yang [Princeton University, NJ (United Stat↗

High-Energy Two-Stage Pulsed Plasma Thruster

A high-energy (28 kJ per pulse) two-stage pulsed plasma thruster (MSFC PPT-1) has been constructed and tested. The motivation of this project is to develop a high power (approximately 500 kW), high specific impulse (approximately 10000 s), highly efficient (greater than 50%) thruster for use as primary propulsion in a high power nuclear electric propulsion system. PPT-1 was designed to overcome four negative characteristics which have detracted from the utility of pulsed plasma thrusters: poor electrical efficiency, poor propellant utilization efficiency, electrode erosion, and reliability issues associated with the use of high speed gas valves and high current switches. Traditional PPTs have been plagued with poor efficiency because they have not been operated in a plasma regime that fully exploits the potential benefits of pulsed plasma acceleration by electromagnetic forces. PPTs have generally been used to accelerate low-density plasmas with long current pulses. Operation of thrusters in this plasma regime allows for the development of certain undesirable particle-kinetic effects, such as Hall effect-induced current sheet canting. PPT-1 was designed to propel a highly collisional, dense plasma that has more fluid-like properties and, hence, is more effectively pushed by a magnetic field. The high-density plasma loading into the second stage of the accelerator is achieved through the use of a dense plasma injector (first stage). The injector produces a thermal plasma, derived from a molten lithium propellant feed system, which is subsequently accelerated by the second stage using mega-amp level currents, which eject the plasma at a speed on the order of 100 kilometers per second. Traditional PPTs also suffer from dynamic efficiency losses associated with snowplow loading of distributed neutral propellant. The twostage scheme used in PPT-I allows the propellant to be loaded in a manner which more closely approximates the optimal slug loading. Lithium propellant was chosen to test whether or not the reduced electrode erosion found in the Lithium Lorentz Force Accelerator (LiLFA) could also be realized in a pulsed plasma thruster. The use of the molten lithium dense plasma injector also eliminates the need for a gas valve and electrical switch; the injector design fulfills both roles, and uses no moving parts to provide, in principle, a highly reliable propellant feed and electrical switching system. Experimental results reported in this paper include: second-stage current traces, high-speed photographic and holographic imaging of the thruster exit plume, and internal mapping of the discharge chamber magnetic field from B-dot probe data. The magnetic field data is used to create a two-dimensional description of the evolution of the current sheet inside the thruster.

Markusic, Tom↗

Modelling general-relativistic plasmas with collisionless moments and dissipative two-fluid magnetohydrodynamics

ABSTRACT Relativistic plasmas are central to the study of black hole accretion, jet physics, neutron star mergers, and compact object magnetospheres. Despite the need to accurately capture the dynamics of these plasmas and the implications for relativistic transients, their fluid modelling is typically done using a number of (overly) simplifying assumptions, which do not hold in general. This is especially true when the mean free path in the plasma is large compared to the system size, and kinetic effects start to become important. Going beyond common approaches used in the literature, we describe a fully relativistic covariant 14-moment based two-fluid system appropriate for the study of electron–ion or electron–positron plasmas. This generalized Israel-Stewart-like system of equations of motion is obtained directly from the relativistic Boltzmann–Vlasov equation. This new formulation can account for non-ideal effects, such as anisotropic pressures and heat fluxes, not present in previous formulations of two-fluid magnetohydrodynamics. We show that a relativistic two-fluid plasma can be recast as a single fluid coupled to electromagnetic fields with (potentially large) out-of-equilibrium corrections. We keep all electron degrees of freedom, which provide self-consistent evolution equations for electron temperature and momentum. The out-of-equilibrium corrections take the form of a collisional 14-moment closure previously described in the context of viscous single fluids. The equations outlined in this paper are able to capture the full two-fluid character of collisionless plasmas found in black hole accretion and flaring processes around compact objects, as well Braginskii-like two-fluid magnetohydrodynamics applicable to weakly collisional plasmas inside accretion discs.

79 ASTRONOMY AND ASTROPHYSICS↗

CAKE: Consistent Automatic Kinetic Equilibrium reconstruction

Magnetohydrodynamic equilibrium reconstructions with kinetic constraints are an essential input of many tokamak stability and transport analysis work- flows. The reconstructions themselves also contain useful information about pressure and current density distributions within the plasma which can be used to drive other calculations. While standardized magnetic-only, and partially kinetic equilibrium reconstruction is widely available and used, re- constructions for well diagnosed machines such as DIII-D, that incorporates a large and heterogeneous set of internal diagnostic measurements yet robustly reaching high numerical accuracy have been a time intensive, manual process. The Consistent Automatic Kinetic Equilibrium reconstruction tool (CAKE) has been developed and implemented at DIII-D for producing low error, kinetically constrained reconstructions without human intervention, and to identify strategies for doing so that would be broadly applicable. Here, the techniques employed to allow CAKE to handle data of limited quality and availability, as well as the full automation of the workflow may assist future efforts to form kinetic constraint profiles in real time plasma control applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Whistler–Alfvén Turbulence in a Nonneutral Ultrarelativistic Pair Plasma

The large-scale dynamics of most conventional space and astrophysical plasmas are predominantly governed by Alfvén modes, which are low-frequency magnetohydrodynamic modes existing in magnetized media. At scales smaller than the ion gyroscale or frequencies exceeding the ion cyclotron frequency, the Alfvén modes transform into kinetic-Alfvén or whistler modes that significantly contribute to plasma dynamics. However, this scenario reverses in nonneutral pair plasmas, such as those found in the magnetospheres of pulsars and magnetars, around rotating black holes, and in their relativistic jets, as well as in certain laboratory plasmas. In these systems, the large-scale dynamics are governed by hybrid whistler–Alfvén modes, which transform into pure Alfvén modes at smaller scales. We derive the nonlinear equations that describe the dynamics of whistler–Alfvén modes in ultrarelativistic nonneutral magnetically dominated pair plasma and discuss the spectrum of turbulence governed by these equations.

79 ASTRONOMY AND ASTROPHYSICS↗

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

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

Kaltsas, Dimitrios A. (ORCID:0000000300769015)↗

Using Kinetic Entropy to Understand Dissipation in Re-connection and Turbulence

This report summarizes results obtained at Space Science Institute (SSI) under a collaborative project funded by the NSF-DOE partnership in Basic Plasma Physics. The overall goal of the project was to perform a study of kinetic entropy as a diagnostic of dissipation in several archetypal models, such as magnetic reconnection and decaying plasma turbulence. The specific task of the SSI group was to perform collisionless and weakly collisional Particle-in-Cell (PIC) simulations modeling these systems. Overall, the project has demonstrated that, as expected, kinetic entropy carries an imprint of weakly collisional processes commonly thought to be associated with dissipation and thus provides an important diagnostic of dissipation. Specific entropy-related metrics have been proposed that could be used as quantitative measures of dissipation. These metrics have been implemented as diagnostics in PIC simulations and their utility, as well as relation to other, previously proposed, measures of dissipation have been demonstrated. Further, PIC simulations were used to investigate a general question of how energy is dissipated in magnetized weakly collisional turbulence in both non-relativistic and relativistic plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Interpretation of electromagnetic modes in the sub-TAE frequency range in JET plasmas with elevated monotonic q-profiles

Recent JET deuterium experiments with an advanced tokamak scenario using an internal transport barrier (ITB) have been performed to clearly observe destabilized toroidicity-induced Alfvén eigenmodes (TAEs) by fast ions; interestingly, these also exhibit unstable electromagnetic (EM) perturbations in the sub-TAE frequency range. We identify such EM perturbations to be beta-induced ion temperature gradient (BTG) eigenmodes and not beta-induced Alfvén eigenmodes nor beta-induced Alfvén acoustic eigenmodes, which are often unstable in such high-beta plasmas with high-power neutral beam injection (NBI). The BTG modes are the most unstable modes due to the high thermal ion temperature gradient related to the ITB, high thermal ion temperature compared to thermal electron temperature (high T i /T⁠ e ), and a high ion beta regime. BTG mode experimental characteristics match the analytical theory, i.e., location in the vicinity of a rational magnetic surface with a low magnetic shear, mode frequency scaling with the ion diamagnetic frequency (⁠ω$^*_i$⁠), and a coupling among Alfvén and drift waves. We also perform linear gyrokinetic simulations with validated plasma profiles and equilibrium, and find a mode kinetically driven by thermal ions with similar characteristics as the experimental BTG modes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Drift-cyclotron loss-cone instability in 3-D simulations of a sloshing-ion simple mirror

The kinetic stability of collisionless, sloshing beam-ion (45° pitch angle) plasma is studied in a three-dimensional (3-D) simple magnetic mirror, mimicking the Wisconsin high-temperature superconductor axisymmetric mirror experiment. The collisional Fokker–Planck code CQL3D-m provides a slowing-down beam-ion distribution to initialize the kinetic-ion/fluid-electron code Hybrid-VPIC, which then simulates free plasma decay without external heating or fuelling. Over 1 – 10 μs, drift-cyclotron loss-cone (DCLC) modes grow and saturate in amplitude. The DCLC scatters ions to a marginally stable distribution with gas-dynamic rather than classical-mirror confinement. Sloshing ions can trap cool (low-energy) ions in an electrostatic potential well to stabilize DCLC, but DCLC itself does not scatter sloshing beam-ions into the said well. Instead, cool ions must come from external sources such as charge-exchange collisions with a low-density neutral population. Manually adding cool ∼1keV ions improves beam-ion confinement several-fold in Hybrid-VPIC simulations, which qualitatively corroborates prior measurements from real mirror devices with sloshing ions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Computational Fluid Dynamics Modeling of Low Temperature Ignition Processes From a Nanosecond Pulsed Discharge at Quiescent Conditions

Recent interest in nonequilibrium plasma discharges as sources of ignition for the automotive industry has not yet been accompanied by the availability of dedicated models to perform this task in computational fluid dynamics (CFD) engine simulations. The need for a low-temperature plasma (LTP) ignition model has motivated much work in simulating these discharges from first principles. Most ignition models assume that an equilibrium plasma comprises the bulk of discharge kernels. LTP discharges, however, exhibit highly nonequilibrium behavior. In this work, a method to determine a consistent initialization of LTP discharge kernels for use in engine CFD codes like CONVERGE is proposed. The method utilizes first principles discharge simulations. Such an LTP kernel is introduced in a flammable mixture of air and fuel, and the subsequent plasma expansion and ignition simulation is carried out using a reacting flow solver with detailed chemistry. Finally, the proposed numerical approach is shown to produce results that agree with experimental observations regarding the ignitability of methane-air and ethylene-air mixtures by LTP discharges.

33 ADVANCED PROPULSION SYSTEMS↗

On the importance of excited state species in low pressure capacitively coupled plasma argon discharges

In the past three decades, first principles-based fully kinetic particle-in-cell Monte Carlo collision (PIC/MCC) simulations have been proven to be an important tool for the understanding of the physics of low pressure capacitive discharges. However, there is a long-standing issue that the plasma density determined by PIC/MCC simulations shows quantitative deviations from experimental measurements, even in argon discharges, indicating that certain physics may be missing in previous modeling of the low pressure radio frequency (rf) driven capacitive discharges. In this work, we report that the energetic electron-induced secondary electron emission (SEE) and excited state atoms play an important role in low pressure rf capacitive argon plasma discharges. The ion-induced secondary electrons are accelerated by the high sheath field to strike the opposite electrode and produce a considerable number of secondary electrons that lead to additional ionizing impacts and further increase of the plasma density. Importantly, the presence of excited state species even further enhances the plasma density via excited state neutral and resonant state photon-induced SEE on the electrode surface. The PIC/MCC simulation results show good agreement with the recent experimental measurements in the low pressure range (1–10 Pa) that is commonly used for etching in the semiconductor industry. At the highest pressure (20 Pa) and driving voltage amplitudes 250 and 350 V explored here, the plasma densities from PIC/MCC simulations considering excited state neutrals and resonant photon-induced SEE are quantitatively higher than observed in the experiments, requiring further investigation on high pressure discharges.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Analysis of Static Spacecraft Floating Potential at Low Earth Orbit (LEO)

Spacecraft floating potential is the charge on the external surfaces of orbiting spacecraft relative to the space. Charging is caused by unequal negative and positive currents to spacecraft surfaces. The charging process continues until the accelerated particles can be collected rapidly enough to balance the currents at which point the spacecraft has reached its equilibrium or floating potential. In low inclination. Low Earth Orbit (LEO), the collection of positive ion and negative electrons. in a particular direction. are typically not equal. The level of charging required for equilibrium to be established is influenced by the characteristics of the ambient plasma environment. by the spacecraft motion, and by the geometry of the spacecraft. Using the kinetic theory, a statistical approach for studying the interaction is developed. The approach used to study the spacecraft floating potential depends on which phenomena are being applied. and on the properties of the plasma. especially the density and temperature. The results from kinetic theory derivation are applied to determine the charging level and the electric potential distribution at an infinite flat plate perpendicular to a streaming plasma using finite-difference scheme.

Herr, Joel L.↗

Role of perturbed parallel magnetic field effects in predicting turbulent transport in NSTX

This study presents analysis of gyrokinetic simulations on the National Spherical Torus Experiment (NSTX) to investigate the effects of electromagnetic fields on plasma turbulence and transport. The simulations, performed with varying levels of fidelity using the gyrokinetic CGYRO code, include electrostatic (ES), single-field electromagnetic (EM1), and two-field electromagnetic (EM2) models. A detailed comparison across the simulation database reveals that electromagnetic effects increase both predicted growth rates and quasilinear fluxes, with EM2 simulations producing stronger turbulence than ES and EM1 cases. Quasilinear modeling using QLGYRO demonstrates that while the perturbed parallel magnetic field (δB ∥ ) does not drastically affect the total flux at experimental gradients, it leads to a shift in the dominant instability, altering mode structures from microtearing to kinetic ballooning modes (KBMs). The proximity of the plasma profiles to the KBM threshold is explored, with the experimental conditions being near the onset of KBM-driven transport. The KBM, with its large growth rates, is identified as a potential driver of electron temperature flattening, as it can rapidly transport heat across flux surfaces. Performing stability analysis shows core-localized unstable a low- mode that could contribute to the flattening at the early times of the discharge. TGYRO predictive modeling, incorporating both TGLF and QLGYRO, indicates that the inclusion of δB ∥ significantly improves the accuracy of temperature profile predictions in NSTX high-beta plasmas, although challenges remain in modeling the sharp flux discontinuities caused by KBM-driven instabilities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Simple model of plasma acceleration in a magnetic nozzle

A collisionless, steady-state, cold-plasma model is used to calculate the three-dimensional trajectory of a plasma as it is accelerated through a diverging magnetic field. The magnetic field is assumed to be axisymmetric with nonzero radial and axial components and zero azimuthal component. Although random thermal motion of plasma species is neglected in the cold plasma approximation, an important effect of plasma thermal energy is accounted for in the model as the kinetic energy of electrons in their Larmor motion about magnetic field lines. Calculations based on this model confirm previous studies which suggested that plasma separation from the field of a magnetic nozzle can take place even in the absence of collisional diffusion. It is shown that plasma divergence in a magnetic nozzle can be controlled by tailoring the field geometry.

Sercel, Joel C.↗

Diffusion regime of electron–electron collisions in weakly ionized plasmas

We consider weakly ionized plasma where an elastic scattering of electrons on neutrals changes the individual acts and the rate of electron-electron collisions significantly. In this case, the kinetics of electron thermalization is very different from that in fully ionized plasma. Here, the colliding electrons do not move freely. They, instead, diffuse because of fast scattering on neutrals. We demonstrate how a proper account of this diffusion enables one to estimate the characteristic time of electron thermalization. We also present a rigorous derivation of the kinetic equation for electrons by using Bogolyubov’s method based on Liouville equations for multi-particle distribution functions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗