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 901 records · Page 50
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

A reaction mechanism for plasma electrolysis of AgNO 3 forming silver nanoclusters and nanoparticles

In plasma-driven solution electrolysis (PDSE), gas-phase plasma-produced species interact with an electrolytic solution to produce, for example, nanoparticles. An atmospheric pressure plasma jet (APPJ) directed onto a liquid solution containing a metallic salt will promote reduction of metallic ions in solution, generating metallic clusters that nucleate to form nanoparticles. In this article, results from a computational investigation are discussed of a PDSE process in which a radio-frequency APPJ sustained in helium impinges on a silver nitrate solution, resulting in growth of silver nanoparticles. A reaction mechanism was developed and implemented in a global plasma chemistry model to predict nanoparticle growth. To develop the reaction mechanism, density functional theory was used to generate probable silver growth pathways up to Ag 9 . Neutral clusters larger than Ag 9 were classified as nanoparticles. Kinetic reaction rate coefficients for thermodynamically favorable growth pathways were estimated based on an existing, empirically determined base reaction mechanism for smaller Ag particle interactions. These rates were used in conjunction with diffusion-controlled reaction rate coefficients that were calculated for other Ag species. The role of anions in reduction of Agn ions in forming nanoparticles is also discussed. Oxygen containing impurities or admixtures to the helium, air entrainment into the APPJ, and dissociation of saturated water vapor above the solution can produce additional reactive oxygen species in solution, resulting in the production of anions and O 2 – in particular. For a given molarity, delivering a sufficient fluence of reducing species will produce similar nanoparticle densities and sizes for all applied power levels. In conclusion, comparisons are made to alternate models for nanoparticle formation, including charged nanoparticles and use of direct current plasmas.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Towards fully predictive gyrokinetic full- f simulations: validation and triangularity studies in TCV

Designing economical magnetic confinement fusion power plants motivates computational tools that can estimate plasma behavior from engineering parameters without direct reliance on experimental measurement of the plasma profiles. In this work, we present full-f global long-wavelength gyrokinetic simulations of edge and scrape-off layer turbulence in tokamaks that use only magnetic geometry, heating power, and particle inventory as inputs. Unlike many modeling approaches that employ free parameters fitted to experimental data, raising uncertainties when extrapolating to reactor scales. This approach directly simulates turbulence and resulting profiles through gyrokinetics without such empirical adjustments. This is achieved via an adaptive sourcing algorithm in Gkeyll that strictly controls energy injection and emulates particle sourcing due to neutral recycling. We show that the simulated kinetic profiles compare reasonably well with Thomson scattering and Langmuir probe data for Tokamak á Configuration Variable (TCV) discharge #65125, and that the simulations reproduce characteristic features such as blob transport and self-organized electric fields. Applying the same framework to study triangularity effects suggests mechanisms contributing to the improved confinement reported for negative triangularity (NT). Simulations of TCV discharges #65125 and #65130 indicate that NT increases the E x B flow shear (by about 20% in these cases), which correlates with reduced turbulent losses and a modest change in the distribution of power exhaust to the vessel wall. While the physical models contain approximations that can be refined in future work, the predictive capability demonstrated here, evolving multiple profile relaxation times with kinetic electron and ion models in hundreds of GPU hours, indicates the feasibility of using Gkeyll to support design studies of fusion devices.

Hoffmann, Antoine Cyril David [Princeton Plasma Ph↗
DOE OSTI2020

Disruption Prediction and Avoidance in High Beta Long Pulse KSTAR Plasmas

DoE Grant DE-SC0016409 (Disruption Prediction and Avoidance in High Beta Long Pulse KSTAR Plasmas) supported a collaboration among Columbia University, the Princeton Plasma Physics Laboratory (PPPL), and the Plasma Science and Fusion Center (PSFC) at MIT to conduct research on the KSTAR tokamak at the National Fusion Research Institute (NFRI) in Korea to improve predictive capability for disruptions in tokamaks. This initiative included upgrades to kinetic diagnostics at KSTAR whose measurements form the basis for the disruption prediction models. The major thrust of the PPPL and PSFC work was implementing a capability at KSTAR to allow partially-polarized background light seen by the KSTAR Motional Stark Effect (MSE) diagnostic to be compensated by ‘wavelength interpolation’ rather than the traditional ‘beam modulation’ technique. This objective was completed on time and on budget, and the MSE Background Polychromator is currently functioning at KSTAR. Primary responsibility for maintaining and operating the MSE-BP was transferred to Nova Photonics, Inc. in 2019.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Understanding cold electron impact on parallel-propagating whistler chorus waves via moment-based quasilinear theory

Earth's magnetosphere hosts a wide range of collisionless particle populations that interact through various wave-particle processes. Among these, cold electrons, with energies below 100 eV, often dominate the plasma density but remain poorly characterized due to measurement challenges such as spacecraft charging and photoelectron contamination. Understanding the contribution of these cold populations to wave–particle interaction is of significant interest. Recent kinetic simulations identified a secondary drift-driven instability, in which parallel-propagating whistler-mode chorus waves excite oblique electrostatic whistler waves near the resonance cone and Bernstein-mode turbulence. These secondary modes enable a new channel of energy transfer from the parallel-propagating whistler wave to the cold electrons. In this work, we develop a moment-based quasilinear theory of the secondary instabilities to quantify such energy exchange. Our results show that these secondary instabilities persist for a wide range of parameters and, in many cases, lead to nearly complete damping of the primary wave. Such secondary instability might limit the amplitude of parallel-propagating whistler waves in Earth's magnetosphere and might explain why high-amplitude oblique whistler or electron Bernstein waves are rarely observed simultaneously with high-amplitude field-aligned whistler waves in the inner magnetosphere.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Kinetic Modeling and Experimental Studies of Capacitively Coupled Macro- and Microplasmas at Microwave Frequencies

The effect of driving frequency in the range of 13.56 MHz to 73 MHz on electron energy distribution and electron heating modes in a 50 mTorr capacitively coupled argon plasma discharge is studied using 1D-3V particle-in-cell simulations. Calculated electron energy probability functions exhibit three distinct temperatures for low-, mid-, and high-energy electrons at all the studied driving frequencies. When compared to published experimental data, the calculated probability functions show a reasonable agreement for the energy range resolved in the measurements (about 2 eV to 10 eV). Discrepancies due to limitations in experimental energy resolution outside this range lead to differences between computational and experimental values of the electron number density determined from the distribution functions, and the predicted effective electron temperature is within 25 % of experimental values. The impedance of the discharge is interpreted in terms of a homogeneous equivalent circuit model and the driving frequency dependence of the inferred combined sheath thickness is found to obey a known, theoretically-derived, power law. The average power transferred from the field to the electrons (electron heating) is computed, and a region of negative heating near the sheath edge, particularly at higher driving frequencies, is identified. Analysis of the electron momentum equation shows that electron inertia, which on temporal averaging would be zero in a linear regime, is responsible for negative values of power deposition near the sheath edge at high driving frequencies due to the highly nonlinear behavior of the discharge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Entrapment Behavior of Solid Surrogate Fission Products at Engineered UN Nano‐Hetero‐Interfaces Within Metallic Nuclear Fuels

Nanometric hetero-interfaces provide a wealth of scientific and engineering opportunities due to their complex and often misunderstood properties that can differ from their respective bulk constituents. In this work, the ability for engineered nanostructures within a bulk U─Mo alloy to arrest simulant fission products is investigated experimentally and computationally. Nanostructured 90 wt% U/ 10 wt% Mo (U-10Mo) with 7.1 at% Nd is consolidated using spark-plasma- sintering (SPS) techniques and is heat-treated at 500 °C under vacuum for 24, 100, 500, and 1000 h. Analysis on the sintered and heat-treated U-10Mo reveals rapid kinetics in Nd diffusion to nanocluster sites, with evidence of Nd diffusion occurring during sintering and during the following heat-treatment. The segregation behavior of Nd at two different U─Mo/UN interfaces is computationally verified using density functional theory (DFT) to reinforce experimental data.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Single-shot ultrafast dynamics of nanosecond pulsed plasmas: Transition from ps–ns nonequilibrium to near-full ionization

The ultrafast multi-stage evolution of state-defining properties in atmospheric-pressure nanosecond pulsed plasmas is quantified using single-discharge, jitter-free, continuous streak-sweep spectroscopy of N 2 (C → B) molecular and N + /O + ionic species with a single-shot time resolution as short as ∼160 ps. An early period of extreme nonequilibrium is observed with vibrational temperatures [T V (C)] dropping from ∼8000 K at <200 ps after plasma breakdown to <4000 K within ∼1 ns, with near-ambient rotational temperatures [T R (C)] ∼ 300 K due to limited collisional energy transfer. This early-time T V (C) trend is representative of a direct and unquenched look into the high-energy tail of the electron energy distribution function; thus, it tracks real-time changes in the mean electron energy via the N2(C) emission signatures. This is followed by the rapid onset of N + /O + ionic emission after a distinct time delay of ∼14.9 ns due to multistep chemical kinetics. The ionic emission enables determination of electron densities (n e ) >2 × 10 19 cm −3 and electron temperatures (T e ) >36 000 K, indicating the transition to a nearly fully ionized regime. The ps–ns temporal dynamics are also compared between air and N 2 plasmas to assess the influence of collisional partners, as well as across the anode, cathode, and central gap regions to identify spatial variations in plasma behavior. Finally, this work demonstrates versatile continuously probing approach capable of spectro-temporal and spatially resolved single-shot measurements of key state variables in the ps–ns evolution of atmospheric-pressure nanosecond pulsed plasmas.

Electric discharges↗
DOE OSTI2023

Jet thermalization in QCD kinetic theory

We perform numerical studies in the framework of QCD kinetic theory to investigate the energy and angular profiles of a high energy parton — as a proxy for a jet produced in heavy ion collisions — passing through a Quark-Gluon Plasma (QGP). We find that the fast parton loses energy to the plasma mainly via a radiative turbulent quark and gluon cascade that transports energy locally from the jet down to the temperature scale where dissipation takes place. In this first stage of the system time evolution, the angular structure of the turbulent cascade is found to be relatively collimated. However, when the lost energy reaches the plasma temperature it is rapidly transported to large angles w.r.t. the jet axis and thermalizes. We investigate the contribution of the soft jet constituents to the total jet energy. We show that for jet opening angles of about 0.3 rad or smaller, the effect is negligible. Conversely, larger opening angles become more and more sensitive to the thermal component of the jet and thus to medium response. Our result showcases the importance of the jet cone size in mitigating or enhancing the details of dissipation in jet quenching observables.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

How eigenmode self-interaction affects zonal flows and convergence of tokamak core turbulence with toroidal system size

Self-interaction is the process by which a microinstability eigenmode that is extended along the direction parallel to the magnetic field interacts non-linearly with itself. This effect is particularly significant in gyrokinetic simulations accounting for kinetic passing electron dynamics and is known to generate stationary $E\times B$ zonal flow shear layers at radial locations near low-order mode rational surfaces (Weikl et al. Phys. Plasmas, vol. 25, 2018, 072305). Here, we find that self-interaction, in fact, plays a very significant role in also generating fluctuating zonal flows, which is critical to regulating turbulent transport throughout the radial extent. Unlike the usual picture of zonal flow drive in which microinstability eigenmodes coherently amplify the flow via modulational instabilities, the self-interaction drive of zonal flows from these eigenmodes are uncorrelated with each other. It is shown that the associated shearing rate of the fluctuating zonal flows therefore reduces as more toroidal modes are resolved in the simulation. In simulations accounting for the full toroidal domain, such an increase in the density of toroidal modes corresponds to an increase in the toroidal system size, leading to a finite system size effect that is distinct from the well-known profile shearing effect.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Density and Magnetic Field Asymmetric Kelvin‐Helmholtz Instability

Abstract The Kelvin‐Helmholtz (KH) instability can transport mass, momentum, magnetic flux, and energy between the magnetosheath and magnetosphere, which plays an important role in the solar‐wind‐magnetosphere coupling process for different planets. Meanwhile, strong density and magnetic field asymmetry are often present between the magnetosheath (MSH) and magnetosphere (MSP), which could affect the transport processes driven by the KH instability. Our magnetohydrodynamics simulation shows that the KH growth rate is insensitive to the density ratio between the MSP and the MSH in the compressible regime, which is different than the prediction from linear incompressible theory. When the interplanetary magnetic field (IMF) is parallel to the planet's magnetic field, the nonlinear KH instability can drive a double mid‐latitude reconnection (DMLR) process. The total double reconnected flux depends on the KH wavelength and the strength of the lower magnetic field. When the IMF is anti‐parallel to the planet's magnetic field, the nonlinear interaction between magnetic reconnection and the KH instability leads to fast reconnection (i.e., close to Petschek reconnection even without including kinetic physics). However, the peak value of the reconnection rate still follows the asymmetric reconnection scaling laws. We also demonstrate that the DMLR process driven by the KH instability mixes the plasma from different regions and consequently generates different types of velocity distribution functions. We show that the counter‐streaming beams can be simply generated via the change of the flux tube connection and do not require parallel electric fields.

Astronomy & Astrophysics↗
DOE OSTI2024

Understanding Environmental Barrier Coating Lifetimes and Performance for Industrial Gas Turbines

Hydrogen or hydrogen blend fuels are expected to replace natural gas in land-based industrial gas turbines (IGTs) to support a greener power economy. Silicon carbide (SiC) base ceramic matrix composites (CMCs) are considered for replacement of Ni-based superalloys to facilitate future efficiency improvements. SiC CMCs require environmental barrier coatings (EBCs) to mitigate volatilization from high-temperature steam, thus making the EBC lifetime critical information for identifying CMC component lifetimes. Here, the goal of this project is to determine the maximum bond coating temperature underneath the EBC for achieving an IGT component lifetime goal of 25,000 h, which is far greater than current CMC component lifetime requirements for aeroturbine applications. To provide data for the lifetime model, laboratory testing used atmospheric plasma-sprayed rare-earth silicate EBCs on monolithic SiC substrates with an intermediate Si bond coating. Specimens exposed to 1-h thermal cycles in flowing air–steam environments and reaction kinetics were assessed from 700 °C to 1350 °C by measuring the thickness of the thermally grown silica scales. The silica growth and phase transformation appear critical in predicting EBC lifetime and several strategies have been explored to reduce the oxide growth rate and improve EBC durability at elevated temperatures. Advanced characterization using Raman spectroscopy has helped clarify this system.

08 HYDROGEN↗
DOE OSTI2020

Recrystallization suppression through dispersion-strengthening of tungsten

Tungsten is the material of choice for the divertor region of future nuclear fusion reactors, an environment that will expose plasma-facing components (e.g. divertor, etc...) to high temperatures and transient high heat flux events. Under these conditions, recrystallization and grain growth of tungsten can occur, leading to undesirable microstructural and mechanical property changes. Therefore, there is a need to raise the recrystallization temperature of tungsten and limit the kinetics of the recrystallization and grain growth processes. In this work, we examine the effect of different types (TiC vs. TaC vs. ZrC) and different concentrations (1.1 vs. 5 vs. 10 wt.%) of dispersed second phase particles in a tungsten matrix on the high temperature performance. Additionally, the addition of second-phase particles effectively increases the temperature of and time for recrystallization and slow grain growth; however, the addition of a high weight fraction of particles alters the surface chemistry, which may impact subsequent plasma-surface interactions. These results show that the addition of small concentrations of dispersed particles can be effectively employed in tungsten to raise the upper operating temperature limit for tungsten in a fusion reactor.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Reconstruction of current profiles from external magnetic and MSE measurements on DIII-D

This work introduces a new current model into the EFIT code for plasma current profile reconstruction. Using this model, edge-localized current density profiles are reconstructed solely from external magnetic measurements. Quantitative comparisons show that edge current profiles obtained with the new current model align closely with those from kinetic equilibrium reconstruction. The new model enables current density reconstruction with fewer diagnostics than kinetic equilibrium reconstruction (Lao et al 1990 Nucl. Fusion 30 1035–49), demonstrating that edge-localized current density information can be obtained from external magnetic measurements only, without the help of kinetic profiles and internal diagnostics, and without using the bootstrap current model. Furthermore, with the constraints of external magnetic and MSE measurements, current profiles in full radius obtained from the new model are in good agreement with the results of kinetic equilibrium reconstruction, which not only improves the efficiency of current reconstruction but also helps avoid possible errors caused by different kinds of diagnoses. Possible applications are also discussed, and a brief analysis of the evolution of the full-radius current density is performed. Noticeable decreases in both the magnitude and gradient of the edge current around the time of ELM bursts are observed. Moreover, analysis of the full radius current density profiles ( ρ N = 0 − 1.0, ρN denotes to square root of the normalized toroidal flux) reveals that as the edge current changes, the inner current also changes following a sinusoidal-like pattern at almost all time slices during four type-I ELMs in discharge 161493. The current evolution in the inner region ( ρ N = 0.2 − 0.6) is synchronized with the boundary region ( ρ N = 0.89 − 1.0), while the intermediate region ( ρ N = 0.6 − 0.89) shows opposite behavior. This sinusoidal-like variations in full-radius current density profile indicates strong core-pedestal coupling, suggesting ELMs significantly affect core plasmas, though current changes manifest differently across radial regions.

EFIT↗
DOE OSTI2025

A neural master equation framework for multiscale modeling of molecular processes: application to atomic-scale plasma processes

Plasma-surface interactions (PSI) play a crucial role in microelectronics fabrication; however, their multiscale nature and array of complex, often unknown interactions make computational modeling of PSIs extremely difficult. To this end, we propose a general neural master equation (NME) framework that uses master equations to describe the dynamics of a molecular process, wherein neural networks learned from atomistic simulations represent unknown transitions between different system states. By leveraging the physics-based structure of master equations and data-driven state transitions, the NME framework promotes generalizability and physics interpretability, and can bridge disparate length and time scales. The framework is demonstrated for multiscale modeling of Si atomic layer etching and reactive ion etching, where the learned NME-based surface kinetic models exhibit good predictive and extrapolative capabilities for predicting experimentally relevant observables as a function of process parameters. The NME-based surface kinetic models obey physical constraints, which are violated in models based on neural ordinary differential equations. The proposed NME framework for multiscale modeling of molecular processes can pave the way for the discovery of new chemistries and materials in atomic-scale plasma processes.

Chemical engineering↗
DOE OSTI2024

Structural transformation of dusty plasma crystal in dc discharge plasma by changing confinement ring bias

We report an experimental study of the structural transition of a stable complex plasma crystal to a solid-liquid phase coexistence by the controlled adjustment of the confinement potential, while keeping all other parameters constant. The experiments are carried out on a tabletop linear dusty plasma experimental device, which consists of a circular powered electrode and an extended grounded cathode plate. A stationary crystal of melamine formaldehyde particles is formed in a background of argon plasma inside a confining ring that is isolated from the cathode by a ceramic cover. The stable crystal structure breaks in the core region and transitions to a coexistent state by carefully changing the confining potential, thereby modifying the sheath structure. The transition is confirmed by evaluating the variation in different characteristic parameters such as the pair correlation function, local bond order parameter, and dust kinetic temperature as a function of confining bias potential. It is found that melting in the core is due to the onset of dust fluctuations in the layers beneath the topmost layer, which grow in amplitude as the confining bias potential is reduced below a threshold value. The present technique of changing confinement provides a unique feature to study structural transitions of plasma crystals without affecting the overall plasma parameters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Structural transformation of dusty plasma crystal in dc discharge plasma by changing confinement ring bias

We report an experimental study of the structural transition of a stable complex plasma crystal to a solid-liquid phase coexistence by the controlled adjustment of the confinement potential, while keeping all other parameters constant. The experiments are carried out on a tabletop linear dusty plasma experimental device, which consists of a circular powered electrode and an extended grounded cathode plate. A stationary crystal of melamine formaldehyde particles is formed in a background of argon plasma inside a confining ring that is isolated from the cathode by a ceramic cover. The stable crystal structure breaks in the core region and transitions to a coexistent state by carefully changing the confining potential, thereby modifying the sheath structure. The transition is confirmed by evaluating the variation in different characteristic parameters such as the pair correlation function, local bond order parameter, and dust kinetic temperature as a function of confining bias potential. It is found that melting in the core is due to the onset of dust fluctuations in the layers beneath the topmost layer, which grow in amplitude as the confining bias potential is reduced below a threshold value. The present technique of changing confinement provides a unique feature to study structural transitions of plasma crystals without affecting the overall plasma parameters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Implementation of a 9-point stencil in SOLPS-ITER and implications for Alcator C-Mod divertor plasma simulations

The SOLPS-ITER code suite is used worldwide for plasma edge modeling, the interpretation of experiments, as well as for the design of the ITER divertor. The numerical scheme of the plasma solver of the code, B2.5, is based on the assumption of perfectly field-aligned grids, while in practice grids are often strongly distorted to match divertor target shapes. Neglecting these grid distortion leads to qualitatively and quantitatively incorrect results for fluid neutral simulations, and may affect results in cold (detached) divertors even when using kinetic neutral simulations. In this contribution, we present the first results of a newly implemented 9-point stencil in B2.5 to properly handle misaligned grids. The new scheme is then applied to fluid neutral simulations of a well-diagnosed and previously modeled Alcator C-Mod discharge. Results are compared with the original 5-point scheme neglecting grid distortion effects, as well as with simulations including a full kinetic neutral model. We conclude that the 9-point stencil is essential to correctly model the transport of fluid neutrals on distorted grids, and to capture the effects of divertor closure on the fluid neutral behavior.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗