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

ZrN Phase Formation, Hardening and Nitrogen Diffusion Kinetics in Plasma Nitrided Zircaloy-4

Plasma nitridation was conducted to modify the surfaces of Zircaloy-4. Scanning electron microscopy (SEM), transmission electron microscopy (TEM), and Raman analysis were used to characterize microstructures and phases. Surface indentation and cross-sectional indentation were performed to evaluate mechanical property changes. Nitridation forms a thin layer of ZrN phase, followed by a much deeper layer affected by nitrogen diffusion. The ZrN phase is confirmed by both TEM and Raman characterization. The Raman peaks of ZrN phase show a temperature dependence. The intensity increases with increasing nitridation temperatures, reaches a maximum at 700 °C, and then decreases at higher temperatures. The ZrN layer appears as continuous small columnar grains. The surface polycrystalline ZrN phase is harder than the bulk by a factor of ~8, and the nitrogen diffusion layer is harder by a factor of ~2–5. The activation energy of nitrogen diffusion was measured to be 2.88 eV. The thickness of the nitrogen-hardened layer is controllable by changing the nitridation temperature and duration.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A deterministic verification strategy for electrostatic particle-in-cell algorithms in arbitrary spatial dimensions using the method of manufactured solutions

As simulations of kinetic plasmas continue to increase in scope and complexity, a rigorous and straightforward method for verifying particle-in-cell (PIC) implementations is necessary to ensure their correctness. Here, in this paper, we present a deterministic method for the rigorous verification of multidimensional, multispecies, electrostatic particle-in-cell codes based on the method of manufactured solutions. Specifically, we prove that rigorous verification is possible through the exclusive examination of errors of grid quantities (i.e., moments and/or fields), allowing for a very light-weight and non-intrusive implementation in existing PIC codes. This is a marked improvement over earlier PIC verification approaches (only demonstrated with one species in 1D-1V), which rely on the comparison of cumulative distribution functions, and are invasive on the code base, introduce additional stochastic noise, are significantly more computationally expensive, and lack rigorous convergence properties. Interestingly, we show that different grid quantities feature different rates of convergence with the number of particles and mesh size, impacting the verification process. These theoretical results are confirmed numerically with a multi-species 2D-2V particle-in-cell code, including the ability of the method to detect order reduction due to an incorrect implementation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Multi-Scale Modeling for Plasma-Enhanced Ammonia Decomposition over Carbides and Nitrides

Ammonia is a carbon-free hydrogen carrier, but its decomposition typically requires high temperatures over costly Ru-based catalysts due to the large barrier for N≡N bond formation. We develop a multiscale framework combining density functional theory, zero-dimensional plasma kinetics, and microkinetic modeling to elucidate how non-thermal plasma (NTP) enables low-temperature NH 3 decomposition over Co-based carbides and nitrides, benchmarked against Ru and Co. Under thermal conditions, all catalysts are limited by N≡N bond formation, with Co 3 C(001) most active owing to its negatively charged surface, strong N* binding, and low activation barriers of N≡N bond formation. Plasma-induced vibrational excitation of NH 3 and its reactive radicals promotes a radical-driven •NH 2 –N* coupling pathway that dominates on Co 3 C(001) and Co 3 N(001), shifting the rate-limiting step to NH 3 (v1) dissociation, increasing turnover frequencies by up to 6 orders of magnitude, and reducing the temperature needed to reach a turnover frequency of 5 s –1 from >680 °C (Ru and Co under thermal condition) to 267 °C (Co 3 C) and 415 °C (Co 3 N). These results identify Co-based carbides and nitrides as promising plasma-active catalysts for energy-efficient hydrogen production from ammonia.

ammonia decomposition↗
NASA NTRS1966

Chemical reactions in electrical plasmas.

Chemical kinetics of electron plasma reactions, discussing energy states, ion-molecule reactions, charge transfer, transport properties, phase interactions, etc

ENERGY LEVEL↗
DOE OSTI2023

Plasma-based global pathway analysis to understand the chemical kinetics of plasma-assisted combustion and fuel reforming

The Global Pathway Analysis (GPA) algorithm helps analyze the chemical kinetics of complex combustion systems by identifying important global reaction pathways connecting a source species to a sink species through various important intermediate species (i.e., hub species). Here, the present work aims to extend GPA algorithm to plasma-assisted combustion and fuel reforming systems to identify the dominant global pathways in such systems at various conditions. In addition, the present study extends the ability of GPA algorithm to identify reaction cycles involving the excitation of high-concentration species (e.g., O 2 , N 2 , and fuel) to their vibrational and electronic states and the subsequent de excitation to their ground state, based on their significance on the reactivity of plasma-assisted systems in terms of gas heating and radical production. Provisions are made in the GPA algorithm to evaluate the reactivity of identified re action pathways and cycles based on the element-flux transfer (i.e., dominance), heat release, and radical production rate. The newly developed Plasma-based Global Pathway Analysis (PGPA) algorithm is then used to analyze the plasma assisted combustion of ammonia and reforming of methane. The PGPA analyses elucidated the significance of vibrational-translational cycles on the reactivity of NH 3 /air mixtures. Further, analyses on the production of NO ascribed the early reforming of NH 3 to N 2 and H 2 in impeding the production of NO during plasma-assisted NH 3 ignition. Lastly, the enhanced reforming of CH 4 /N 2 mixtures using plasma has been attributed to electron impact dissociation of CH 4 when compared to thermal reforming. In contrast, conventional path-Flux analysis (PFA) was found to require significant manual effort and pre-analysis intuitions from expert knowledge, making it arduous to provide valuable in sights into plasma chemistry. The user-friendly and automated nature of PGPA thus provides a valuable tool for assessing the kinetics of plasma-assisted systems helpful in analyzing and, further, a foundation in reducing plasma-assisted chemistry, without the needs of expert knowledge.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2024

A Study on Particle Trajectory Error in Finite-Element Particle-in-Cell Algorithms

Particle-in-cell (PIC) algorithms are widely used for the simulation of kinetic plasmas. PIC algorithms account for the interaction between charged particles in a plasma and the electromagnetic field in ambient space, including self-field interactions. The objective of this article is to study the error in charged particle trajectories present in finite-element (FE)-based PIC algorithms on unstructured meshes. We study how the trajectory error behaves according to the FE mesh resolution and the matrix solver employed in the FE algorithm. The study is performed by considering a trajectory established by a parabolic electric potential and an axial magnetic force acting on the charged particle. Under a proper choice combination of electric and magnetic field parameters, the 2-D particle trajectories comprise closed orbits. Numerical errors cause small perturbations on the trajectories, with cumulative effects. As a result, the resulting orbital trajectories exhibit distortions including spurious apsidal precession. These distortions provide a clear imprint of the numerical error. Here, we also study the numerical error in a quantitative fashion by computing the distance norm between the trajectories obtained by the exact fields and numerical fields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

An unstructured body-of-revolution electromagnetic particle-in-cell algorithm with radial perfectly matched layers and dual polarizations

A novel electromagnetic particle-in-cell algorithm has been developed for fully kinetic plasma simulations on unstructured (irregular) meshes in complex body-of-revolution geometries. The algorithm, implemented in the BORPIC++ code, utilizes a set of field scalings and a coordinate mapping, reducing the Maxwell field problem in a cylindrical system to a Cartesian finite element Maxwell solver in the meridian plane. The latter obviates the cylindrical coordinate singularity in the symmetry axis. The choice of an unstructured finite element discretization enhances the geometrical flexibility of the BORPIC++ solver compared to the more traditional finite difference solvers. Symmetries in Maxwell’s equations are explored to decompose the problem into two dual polarization states with isomorphic representations that enable code reuse. The particle-in-cell scatter and gather steps preserve charge conservation at the discrete level. Our previous algorithm (BORPIC+) discretized the E and B field components of TE Φ and TM Φ polarizations on the finite element (primal) mesh. Here, we employ a new field-update scheme. Using the same finite element (primal) mesh, this scheme advances two sets of field components independently: (1) E and B of TE Φ polarized fields, (E z , E ρ , B Φ ) and (2) D and H of TM Φ polarized fields, (D Φ , H z , H ρ ). Since these field updates are not explicitly coupled, the new field solver obviates the coordinate singularity, which otherwise arises at the cylindrical symmetric axis, ρ = 0 when defining the discrete Hodge matrices (generalized finite element mass matrices). Here, a cylindrical perfectly matched layer is implemented as a boundary condition in the radial direction to simulate open space problems, with periodic boundary conditions in the axial direction. We investigate effects of charged particles moving next to the cylindrical perfectly matched layer. We model azimuthal currents arising from rotational motion of charged rings, which produce TMΦ polarized fields. Several numerical examples are provided to illustrate the first application of the algorithm.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Temporally-consistent koopman autoencoders for forecasting dynamical systems

Absence of sufficiently high-quality data often poses a key challenge in data-driven modeling of high-dimensional spatio-temporal dynamical systems. Koopman Autoencoders (KAEs) harness the expressivity of deep neural networks (DNNs), the dimension reduction capabilities of autoencoders, and the spectral properties of the Koopman operator to learn a reduced-order feature space with simpler, linear dynamics. However, the effectiveness of KAEs is hindered by limited and noisy training datasets, leading to poor generalizability. To address this, we introduce the Temporally-Consistent Koopman Autoencoder (tcKAE), designed to generate accurate long-term predictions even with limited and noisy training data. This is achieved through a consistency regularization term that enforces prediction coherence across different time steps, thus enhancing the robustness and generalizability of tcKAE over existing models. We provide analytical justification for this approach based on Koopman spectral theory and empirically demonstrate tcKAE’s superior performance over state-of-the-art KAE models across a variety of test cases, including simple pendulum oscillations, kinetic plasma, and fluid flow data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Electron-scale temperature gradients in kinetic equilibrium: MMS observations and Vlasov–Maxwell solutions

For an ideal gas, the notion of thermodynamic equilibrium requires that the system’s temperature be spatially uniform. We demonstrate that for a collisionless plasma, kinetic equilibrium can be achieved even in the presence of a temperature gradient. Motivated by recent MMS magnetopause observations of thin current layers exhibiting electron-scale temperature gradients, here we present an exact solution to the Vlasov-Maxwell system which accounts for the observed electron temperature gradient and reproduces much of the plasma’s spatial variation along one dimension. Furthermore, the current layer in the Vlasov-Maxwell model is self-consistently supported by a prominent crescent shaped electron velocity distribution that naturally arises due to the temperature transition, which agrees with the electron distribution structures measured by MMS during two asymmetric current layers encountered at the magnetopause.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Reconstruction and interpretation of ionization asymmetry in magnetic confinement via synthetic diagnostics

Abstract Strong poloidal refueling asymmetry in the DIII-D tokamak is inferred from line radiation measurements. Synthetic diagnostics in neutral transport modeling coupled to gyrokinetic simulations illuminate implications for the plasma flow profile in the scrape-off layer of single-null beam-driven discharges. Recycling occurs primarily either on the inner or outer divertor legs, depending on the toroidal magnetic field direction. By reversing the toroidal magnetic field, the observed line radiation asymmetry is nearly eliminated or reversed. It is determined that, while relatively simple physics can describe the observed ionization asymmetry, predicting the overall brightness of the hydrogenic Lyman- α signal requires detailed simulation of the plasma and resulting turbulence. To this end, kinetic plasma simulations fully coupled to comprehensive neutral transport calculations—a novel capability—provide first-principles reproduction of Lyman- α observations on DIII-D.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Closure models for the feedback of energetic particles on plasma turbulence

Energetic particles interact with the plasma surrounding them, resonating with certain types of plasma waves to stabilize them while destabilizing others, and changing the character of the background turbulence in ways that have not been fully quantified or understood. Interaction with the turbulent background plasma is key to the acceleration of many types of energetic particles including high-energy cosmic rays, solar energetic particles, and pick-up ions. The acceleration of particles is a process that would ideally be described by a kinetic model, a type of model that follows a probability distribution function (PDF) for all particles in 7-dimensional (x, y, z, v x , v y , v z , t) space. Because of the high dimensionality of a kinetic model, simulations that solve kinetic equations use the largest computational resources currently available, and are yet unable to simulate a realistic number of particles, reach the large scales necessary for astrophysical problems, and use high-precision numerical methods. Two available alternatives to kinetic plasma models have been explored for this problem, with limited success. One is a multi-fluid model produced by a cumulant discarding closure, which evolves coupled equations for the velocity, magnetic field, and internal energy for both the background plasma and the fluid of energetic particles. However, simulations that solve multi-fluid magnetohydrodynamic (MHD) equations are able to include the interaction with energetic particles only in crude ways, typically as an add-on pressure term. The second alternative is to use a hybrid method to couple a fluid description of the background plasma to a kinetic model or a Fokker–Planck model for the energetic particles. These methods are hampered by the physical modeling of the coupling. In this work, we develop a new model, which follows the PDF for all particles; this can be viewed as a step toward physical realism above a multi-fluid MHD model, while also being more computationally efficient than a kinetic model. The equations we develop model both the background plasma and the energetic particles self-consistently. Over the last decade, similar PDF methods have been developed to a high level of sophistication to model reactive flows and turbulent combustion for engineering applications. For treatment of the feedback of the energetic particles on a background plasma, a PDF closure approach should evaluate the mean characteristics, including the density, with better statistical quality than will particle-sampling procedures.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1982

Ion distribution effects of turbulence on a kinetic auroral arc model

An inverted-V auroral arc structure plasma-kinetic model is extended to phenomenologically include the effects of electrostatic turbulence, with k-parallel/k-perpendicular being much less than unity. It is shown that, unless plasma sheet ions are very much more energetic than the electrons, anomalous resistivity is not a large contributor to parallel electrostatic potential drops, since the support of the observed potential drop requires a greater dissipation of energy than can be provided by the plasma sheet. Wave turbulence can, however, be present, with the ion cyclotron turbulence levels suggested by the ion resonance broadening saturation mechanism of Dum and Dupree (1970) being comparable to those observed on auroral field lines. The diffusion coefficient and net growth rate are much smaller than estimates based solely on local plasma properties.

Cornwall, J. M.↗