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 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Multi-Fluid and Kinetic Models of Partially Ionized Magnetic Reconnection

Magnetic reconnection in partially ionized plasmas is a ubiquitous and important phenomena 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\'en speed, $v_A^\star$. These differences are interpreted as the result of operating in different collisional regimes. Multi-fluid simulations are found to maintain $\nu_{ni}L/v_A^\star \gtrsim 1$, where $\nu_{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 $\nu_{ni}L/v_A^\star < 1$, decoupling neutrals from the reconnection outflow. Differences in the way reconnection is triggered may explain these discrepancies.

Magnetic reconnection↗
DOE OSTI2024

Spatio-temporal electric field distributions in an atmospheric plasma jet impinging on a microchannel array surface

The electric field distribution in the ionization waves (IWs) propagating over a microchannel array dielectric surface, with the channels either empty or filled with distilled water, is measured by ps electric field induced second harmonic generation. The surface IW is initiated by the atmospheric pressure N 2 –Ar plasma jet impinging on the surface vertically and powered by ns pulse discharge bursts. The results show that the electric field inside the microchannels, specifically its horizontal component, is enhanced by up to a factor of 2. The field enhancement region is localized within the channels. The vertical electric field inside the channels lags in time compared to the field measured at the ridges, indicating the transient reversal of the IW propagation direction across the channels (toward the jet). This is consistent with the phase-locked plasma emission images and confirmed by the kinetic modeling predictions, which show that the IW ‘jumps’ over the empty channels and propagates into the channels only after the jump between the adjacent ridges. When the channels are filled with water, the wave speed increases by up to 50%, due to the higher effective dielectric constant of the surface. No evidence of a significant electric field enhancement near the dielectric surface (ceramic or water) has been detected, within the spatial resolution of the present diagnostic, ~100 μ m.

plasmas and liquids↗
DOE OSTI2023

Machine learning-based real-time kinetic profile reconstruction in DIII-D

Abstract Kinetic equilibrium reconstruction plays a vital role in the physical analysis of plasma stability and control in fusion tokamaks. However, the traditional approach is subjective and prone to human biases. To address this, the consistent automatic kinetic equilibrium reconstruction (CAKE) method was introduced, providing objective results. Nonetheless, its offline nature limits its application in real-time plasma control systems (PCSs). To address this limitation, we present RTCAKENN, a machine learning model that approximates 7 CAKE-level output profiles, namely pressure, inverse q , toroidal current density, electron temperature and density, carbon ion impurity temperature and rotation profiles, using real-time available inputs. The deep neural network consists of an encoder layer, where the scalars and interdependent inputs such as plasma boundary coordinates and motional Stark effect data are encoded using multi-layer perceptrons (MLPs), while profile inputs are encoded by 1D convolutional layers. The encoded data is passed through a MLP for latent feature extraction, before being decoded in the decoding layers, which consist of upsampling and convolutional layers. RTCAKENN has been implemented in the DIII-D PCS and our model achieves accuracy comparable to CAKE and surpasses existing real-time alternatives. Through clever dropout training, RTCAKENN exhibits robustness and can operate even in the absence of Thomson scattering data or charge exchange recombination data. It executes in under 8 ms in the real-time environment, enabling future application in real-time control and analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Electric field distribution in a non-self-sustained RF discharge with ionization generated by Ns discharge pulses

Abstract Electric field in a capacitively coupled, non-self-sustained RF discharge nitrogen plasma with external ionization generated by high-voltage ns pulses has been measured by ps Electric field induced second harmonic generation (EFISH). The measurements are made both in the bulk of the plasma and in the sheaths, using the absolute calibration by the Laplacian field between two plane electrodes. The results are compared with the kinetic modeling calculations. The RF electric field in the sheaths is significantly higher compared to that in the plasma, due to the displacement of the electrons by the drift oscillations and the resultant plasma self-shielding. However, the kinetic modeling predictions indicate that the electron impact ionization in the sheaths is largely ineffective, due to the low electron density. The reduction of the electric field in the plasma by the self-shielding in the sheaths is moderate, such that the energy is coupled to the plasma by the below-breakdown RF field. The peak RF field in the plasma is in the range of 15–25 Td, indicating the efficient vibrational excitation of N 2 by electron impact. The modeling predictions suggest that the targeted vibrational excitation of molecular species in a non-self-sustained RF discharge with external ionization is scalable to high pressures, electron densities, and discharge powers. The present approach can be extended to the vibrational excitation of other molecular species where the vibrational relaxation is relatively slow, such as CO, CO 2 , and H 2 .

Physics↗
DOE OSTI2021

Suppressing the enhancement of stimulated Raman scattering in inhomogeneous plasmas by tuning the modulation frequency of a broadband laser

The stimulated Raman scattering (SRS) instability can inhibit the performance of laser driven inertial confinement fusion (ICF) implosions by scattering light into unwanted directions or by generating hot electrons that preheat the target fuel. In principle, ICF target designs can avoid parameter regimes conducive to large, linear SRS gains. In practice, kinetic inflation—the nonlinear enhancement of SRS due to electron trapping in the excited plasma wave—makes this difficult. Here we show that laser bandwidth in the form of frequency modulation can either decrease or increase the inflationary SRS (iSRS) threshold in inhomogeneous plasmas depending on the maximum chirp of the laser pulse. The threshold, mapped out by a series of particle-in-cell simulations, exhibits a minimum when the frequency change within the pulse cancels the spatial detuning due to density inhomogeneities along the trajectory of the scattered light. In conclusion, by tuning the pump laser parameters away from this minimum, the iSRS threshold can be larger than at zero bandwidth, providing a path to mitigating kinetic inflation in ignition relevant plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A physics-constrained deep learning treatment of runaway electron dynamics

An adjoint formulation leveraging a physics-informed neural network (PINN) is employed to advance the density moment of a runaway electron (RE) distribution forward in time. A distinguishing feature of this approach is that once the adjoint problem is solved, its solution can be used to project the RE density forward in time for an arbitrary initial momentum space distribution of REs. Furthermore, by employing a PINN, a parametric solution to the adjoint problem can be learned. Thus, once trained, this adjoint-deep learning framework is able to efficiently project the RE density forward in time across various plasma conditions while still including a fully kinetic description of RE dynamics. As an example application, the temporal evolution of the density of primary electrons is studied, with particular emphasis on evaluating the decay of a RE population when below threshold. Predictions from the adjoint-deep learning framework are found to be in good agreement with a traditional relativistic electron Fokker–Planck solver, for several distinct initial conditions, and across an array of physics parameters. Once trained, the PINN thus provides a means of generating RE density time histories with exceptionally low online execution time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1992

Observation and simulation of the radial evolution and stream structure of solar wind turbulence

This paper presents a brief overview of the observed evolution in a variety of quantities describing the turbulent evolution of the interplanetary plasma and describes simulation results consistent with many features of the evolution. The turbulence is manifested through a dissipation at small scales in the inner heliosphere with a corresponding evolution in the breakpoint between a relatively flat and a Kolmogoroff spectrum; an evolution from kinetically to (slightly) magnetically dominated energy of the plasma fluctuations; a general decrease in the cross helicity or 'Alfvenicity'; changes in the anisotropy of the fluctuations; and the increasing predominance of quasi-pressure-balanced structures in the compressive component of the fluctuations. MHD simulations with shear layers either side of a central current sheet show that even in the absence of compressibility the lack of a mean field along the direction of the main flow in the current sheet leads to rapid nonlinear evolution and the observed characteristics of 'Elsasser spectra' of the fields in the inner heliosphere. Adding compressibility to the simulations does not greatly change the 'incompressive' quantities but leads in addition to observed correlations between a measure of compression and other quantities.

Roberts, D. A.↗
DOE OSTI2026

Partially Ionized Plasma Physics and Technological Applications

Partially ionized plasma physics has attracted increased attention recently due to numerous technological applications made possible by the increased sophistication of computer modelling, the depth of the theoretical analysis, and the technological applications to a vast field of manufacturing for computer components. Partially ionized plasma is characterized by a significant presence of neutral particles in contrast to the fully ionized plasma. The theoretical analysis is based upon solutions of the kinetic Boltzmann equation, yielding the non-Maxwellian electron energy distribution function (EEDF), thereby emphasizing the difference with a fully ionized plasma. The impact of the effect on discharges in inert and molecular gases is described in detail, yielding the complex nonlinear phenomena resulting in plasma selforganization. A few examples of such phenomena are given, including the non-monotonic EEDFs in the discharge afterglow in a mixture of argon with the molecular gas NF3; the explosive generation of cold electron populations in capacitive discharges, hysteresis of EEDF in inductively coupled plasmas. Recently, highly advanced computer codes were developed in order to address the outstanding challenges in plasma technology. These developments are briefly described in general terms.

non-Maxwellian electron energy distribution functi↗
NASA NTRS1995

Solar wind theory and modelling

The outflow of coronal plasma into interplanetary space is a consequence of the coronal heating process. Therefore the formation of the corona and the acceleration of the solar wind should be treated as a single problem. Traditionally the mass or particle flux emanating from the extended corona has been thought of as being determined by the coronal temperature or scale height and the coronal (base) density. This argument follows from considerations of the momentum balance of the corona-wind system from which one obtains models of a close to hydrostatic corona out to the critical point where the flow becomes supersonic. With this approach to the acceleration of the wind is has been difficult to reconcile the relatively small variation observed in the proton flux at 1 AU with the predicted exponential dependence of the proton flux on the coronal temperature. In this talk we would like to emphasize another approach in which coronal energetics play the primary role. The deposition of energy into the corona through some 'mechanical' energy flux is balanced by the various energy sinks available to the corona and the sum of these processes determine the coronal structure, i.e. its temperature and density. The corona loses energy through heat conduction into the transition region, through radiative losses, and through the gravitational potential energy and kinetic energy put into the solar wind itself. We will show from a series of models of the chromosphere transition region-corona-solar wind system that most of the energy deposited in a magnetically open region will go into the solar wind, with roughly half going into kinetic energy and half into lifting the plasma out of the solar gravity field. The coronal base density will adjust itself in such a way that the heat conductive flux flowing into the transition region is radiated away in the upper chromosphere. The coronal temperature is set by the requirements that most of the deposited energy goes into accelerating the solar wind; the coronal scale height will adjust itself so that the solar wind energy losses conform to the amplitude of the input energy. These processes are modified by the 'mode' of energy deposition, and we will show the effects on coronal structure of changing the parameters describing coronal heating as well as the effects of including a helium fluid in the models. However, the location, scale height and/or form of the energy deposition (i.e. heating or direct acceleration) are not too important for the solar wind, the coronal density and temperature structure will vary with the 'mode' of energy deposition, but the solar wind mass flux depends mainly on the amplitude of the energy flux.

Hansteen, Viggo H.↗
DOE OSTI2021

Koopman wavefunctions and Clebsch variables in Vlasov–Maxwell kinetic theory

Motivated by recent discussions on the possible role of quantum computation in plasma simulations, here, we present different approaches to Koopman's Hilbert-space formulation of classical mechanics in the context of Vlasov–Maxwell kinetic theory. The celebrated Koopman–von Neumann construction is provided with two different Hamiltonian structures: one is canonical and recovers the usual Clebsch representation of the Vlasov density, the other is non-canonical and appears to overcome certain issues emerging in the canonical formalism. Furthermore, the canonical structure is restored for a variant of the Koopman–von Neumann construction that carries a different phase dynamics. Going back to van Hove's prequantum theory, the corresponding Koopman–van Hove equation provides an alternative Clebsch representation which is then coupled to the electromagnetic fields. Finally, the role of gauge transformations in the new context is discussed in detail.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Experimental and computational study of phase space dynamics in strongly coupled plasmas with steep density gradients

Understanding how plasmas thermalize when density gradients are steep remains a fundamental challenge in plasma physics, with direct implications for fusion experiments and astrophysical phenomena. Standard hydrodynamic models break down in these regimes, and kinetic theories make predictions that have never been directly tested. Here, we present the first detailed phase-space measurements of a strongly coupled plasma as it evolves from sharp density gradients to thermal equilibrium. Using laser-induced fluorescence imaging of an ultracold calcium plasma, we track the complete ion distribution function f(x,v,t). We discover that commonly used kinetic models (Bhatnagar–Gross–Krook and Lenard–Bernstein) overpredict thermalization rates, even while correctly capturing the initial counterstreaming plasma formation. Our measurements reveal that the initial ion acceleration response scales linearly with electron temperature, and that the simulations underpredict the initial ion response. In our geometry we demonstrate the formation of well-controlled counterpropagating plasma beams. This experimental platform enables precision tests of kinetic theories and opens new possibilities for studying plasma stopping power and flow-induced instabilities in strongly coupled systems.

Bergeson, Scott (ORCID:0000000231249226)↗
DOE OSTI2025

Analytical models for the enhancement of fusion reactivity by turbulence

The reactivity of fusion plasma depends not only on its local density and temperature but also, through a recently identified kinetic effect, on the relative velocities of nearby fluid elements. Turbulence on fine spatial scales, therefore, enhances fusion reactivity. The enhancement is quantified here for general subsonic turbulent flows. Leveraging this effect in the design of inertial confinement fusion experiments could enable substantial energy savings.

Fluid dynamics↗
DOE OSTI2020

Continued Validation Studies using the MOOSE Framework for Plasma Simulation with Electromagnetics

Resources and tools for the modeling and simulation of low-temperature plasma (LTP) discharges are increasingly vital to progress in the field in order to properly characterize and study complex source designs and plasma chemistries beyond the scope of traditional diagnostics. Open-source software provides powerful platforms for this work and can enable community-driven LTP R&D. Within the Multiphysics Object-Oriented Simulation Environment (MOOSE) open-source framework [1], capabilities have been demonstrated in the areas of plasma fluids (Zapdos [2]), plasma chemistry (CRANE [3]), and general electromagnetic wave theory (Electromagnetic Library for Kinetics & fluids [ELK] [4]). ELK has since been coupled to Zapdos/CRANE to enable fully coupled electromagnetic plasma simulations. This talk will detail the continued validation efforts and discuss Zapdos-ELK-CRANE code coupling with various-low temperature plasma sources. The impact of fully coupled electromagnetics on process parameters (e.g., temperature and electron/ion energy) versus an electrostatic description will also be discussed. [1] Permann et al., SoftwareX, 11 (2020) 100430. [2] Lindsay et al., J. Phys. D: Appl. Phys. 49 (2016) 235204. [3] Keniley et al., Talk, GEC Session DT3.00002 (2019). [4] Icenhour et al., Talk, GEC Session RR1.00006 (2019). Acknowledgements: INL Graduate Fellowship Program, NSF SI2 Grant 1740300, U.S. Dept. of Energy Office of Science Graduate Student Research Program

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1992

On the consequences of bi-Maxwellian plasma distributions for parallel electric fields

The objective is to use the measurements of the equatorial particle distributions to obtain the parallel electric field structure and the evolution of the plasma distribution function along the field line. Appropriate uses of kinetic theory allows us to use the measured ( and inferred) particle distributions to obtain the electric field, and hence the variation on plasma density along the magnetic field line. The approach, here, is to utilize the adiabatic invariants, and assume the plasma distributions are in equilibrium.

Olsen, Richard C.↗
DOE OSTI2023

MAS: A versatile Landau-fluid eigenvalue code for plasma stability analysis in general geometry

We have developed a new global eigenvalue code, Multiscale Analysis for plasma Stabilities (MAS), for studying plasma problems with wave toroidal mode number (n) and frequency (ω) in a broad range of interest in general tokamak geometry, based on a five-field Landau-fluid description of thermal plasmas. Beyond keeping the necessary plasma fluid response, we further retain the important kinetic effects including diamagnetic drift, ion finite Larmor radius, finite parallel electric field (E||), ion and electron Landau resonances in a self-consistent and non-perturbative manner without sacrificing the attractive efficiency in computation. The physical capabilities of the code are evaluated and examined in the aspects of both theory and simulation. In theory, the comprehensive Landau-fluid model implemented in MAS can be reduced to the well-known ideal MHD model, electrostatic ion-fluid model, and drift-kinetic model in various limits, which clearly delineates the physics validity regime. In simulation, MAS has been well benchmarked with theory and other gyrokinetic and kinetic-MHD hybrid codes in a manner of adopting the unified physical and numerical framework, which covers the kinetic Alfv\'en wave (KAW), ion sound wave (ISW), low-n kink, high-n ion temperature gradient mode (ITG) and kinetic ballooning mode (KBM). Moreover, MAS is successfully applied to model the Alfv\'en eigenmode (AE) activities in DIII-D discharge #159243, which faithfully captures the frequency sweeping of reversed shear Alfv\'en eigenmode (RSAE), the tunneling damping of toroidal Alfv\'en eigenmode (TAE), as well as the polarization characteristics of kinetic beta-induced Alfv\'en eigenmode (KBAE) and beta-induced Alfv\'en-acoustic eigenmode (BAAE) being consistent with former gyrokinetic theory and simulation. With respect to the key progress contributed to the community, MAS has the advantage of combining rich physics ingredients, realistic global geometry and high computation efficiency together for plasma stability analysis in linear regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRSDate not supplied

Plasma Abatement of Volatile Organic Compounds

One difficulty with crewed space operations is the potential for the contamination of the spacecraft habitable volume from volatile organic compounds (VOCs) that may be present. One possible VOC abatement process is to apply a plasma to the contaminated gases. Plasmas provide high energy electrons and ions that are capable of breaking down gaseous organic species generally into smaller compounds. To explore plasma VOC abatement, the authors constructed a plasma gas system that recirculates a particular VOC-laden gas mixture through a plasma torch. This allows for the periodic collection of aliquots of gaseous samples for subsequent analyses via gas chromatography mass spectrometry, designed to quantify the remaining VOC materials. The VOCs evaluated included 100 ppm of acetone, benzene, ethanol, and pentane in a background of carbon dioxide. The plasma conversion of all analytes followed first-order rate kinetics, and VOC elimination was observed after 10 min or less of treatment for all analytes measured. Additionally, it was observed that the plasma caused the conversion of CO 2 into oxygen and CO.

Joel A. Olson↗
NASA NTRS1987

Kinetic aspects of tail dynamics - Theory and simulation

Kinetic theories relevant to the geomagnetic tail are reviewed. The topics discussed include kinetic instabilities, simulations, and current-sheet particle acceleration. Tearing mode and reconnection theories are emphasized. Kinetic treatment is appropriate for these topics since the tail plasma is collisionless. Fluid calculations are appropriate when stochastic processes dominate and for studies where long wavelengths are important. However, fluid treatments of tearing modes and reconnection require a finite resistivity in the diffusion region. Thus, although 'anomalous resistivity' can be guessed or in some cases calculated, ideally the kinetic treatment is often to be preferred. Particle motion and acceleration in the current sheet can give rise to beam-like distributions in the plasma-sheet boundary layer. Studies of current-sheet particle motion have also been used as the basis for 'kinetic' tail equilibrium models. Furthermore, quite recently current-sheet particle motion is used directly in Coroniti's explosive tail reconnection model. The 'inertial conductivity' from the equilibrium models provides the 'dissipation' necessary for reconnection.

Speiser, T. W.↗