SEARCH · Search NASA

Results for “ELECTROSTATIC PLASMA”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Nonlinear susceptibilities for weakly turbulent magnetized plasma: Electrostatic approximation

The plasma weak turbulence theory is a perturbative nonlinear theory, which has been proven to be quite valid in a number of applications. However, the standard weak turbulence theory found in the literature is fully developed for highly idealized unmagnetized plasmas. As many plasmas found in nature and laboratory are immersed in a background static magnetic field, it is necessary to extend the existing discussions to include the effects of ambient magnetic field. Such a task is quite formidable, however, which has prevented fundamental and significant progresses in the subject matter. The central difficulty lies in the formulation of the complete nonlinear response functions for magnetized plasmas. The present paper derives the nonlinear susceptibilities for weakly turbulent magnetized plasmas up to the third order nonlinearity, but in doing so, a substantial reduction in mathematical complexity is achieved by the use of Bessel function addition theorem (or sum rule). The present paper also constructs the weak turbulence wave kinetic equation in a formal sense. For the sake of simplicity, however, the present paper assumes the electrostatic interaction among plasma particles. Fully electromagnetic generalization is a subject of a subsequent paper.

Physics↗

DOE OSTI2026

Local reduced-order modeling for electrostatic plasmas by physics-informed solution manifold decomposition

Despite advancements in high-performance computing and modern numerical algorithms, computational cost remains prohibitive for multi-query kinetic plasma simulations. Here, in this work, we develop data-driven reduced-order models (ROMs) for collisionless electrostatic plasma dynamics, based on the kinetic Vlasov-Poisson equation. Our ROM approach projects the equation onto a linear subspace defined by the proper orthogonal decomposition (POD) modes. We introduce an efficient tensorial method to update the nonlinear term using a precomputed third-order tensor. We capture multiscale behavior with a minimal number of POD modes by decomposing the solution manifold into multiple time windows and creating temporally local ROMs. We consider two strategies for decomposition: one based on the physical time and the other based on the electric field energy. Applied to the 1D1V Vlasov–Poisson simulations, that is, prescribed E-field, Landau damping, and two-stream instability, we demonstrate that our ROMs accurately capture the total energy of the system both for parametric and time extrapolation cases. The temporally local ROMs are more efficient and accurate than the single ROM. In addition, in the two-stream instability case, we show that the energy-windowing reduced-order model (EW-ROM) is more efficient and accurate than the time-windowing reduced-order model (TW-ROM). With the tensorial approach, EW-ROM solves the equation approximately 90 times faster than Eulerian simulations while maintaining a maximum relative error of 7.5% for the training data and 11% for the testing data.

Electrostatic plasmas↗

DOE OSTI2024

Formulation of a one-dimensional electrostatic plasma model for testing the validity of kinetic theory

Here, we present a one-dimensional (1-D) model composed of aligned, electrostatically interacting charged disks, conceived to address in a computable model the validity of the Bogoliubov assumption on the decay of particle correlations in the Born–Bogoliubov–Green–Kirkwood–Yvon hierarchy. This assumption is a basic premise of plasma kinetic theory. The disk model exhibits spatially 1-D features at short distances, but retains 3-D features at large distances. Here the collective dynamics of this model plasma is investigated by solving the corresponding Vlasov equation. In addition, the implementation of the model for the numerical validation of the Bogoliubov assumption is formulated.

1-D plasma model↗

DOE OSTI2025

Similarity for downscaled kinetic simulations of electrostatic plasmas: Reconciling the large system size with small Debye length

A simple similarity has been proposed for kinetic (e.g., particle-in-cell) simulations of plasma transport that can effectively address the long-standing challenge of reconciling the tiny Debye length with the vast system size. This applies to both transport in unmagnetized plasma and parallel transport in magnetized plasmas, where the characteristics length scales are given by the Debye length, collisional mean free paths, and the system or gradient lengths. The controlled scaled variables are the configuration space, x/L, and an artificial Coulomb Logarithm, L ln Λ, for collisions, while the scaled time, t/L, and electric field, LE, are automatic outcomes. The similarity properties are examined, demonstrating that the macroscopic transport physics is preserved through a similarity transformation while keeping the microscopic physics at its original scale of Debye length. To showcase the utility of this approach, two examples of 1D plasma transport problems were simulated using the VPIC code: the plasma thermal quench in tokamaks [Li et al., Nuclear Fusion 63, 066030 (2023)] and the plasma sheath in the high-recycling regime [Li et al., Physics of Plasmas 30, 063505 (2023)].

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Kinetic theory of particle-in-cell simulation plasma and the ensemble averaging technique

Abstract We derive the kinetic theory of fluctuations in physically and numerically stable particle-in-cell (PIC) simulations of electrostatic plasmas. The starting point is the single-time correlation at the start of the simulation between the statistical fluctuations of the weighted densities of macroparticle centers in the plasma particle phase-space. The fluctuations are associated with different initial conditions, typically due to the random initial conditions (in velocity space) of the macroparticles/simulation plasma, assigned according to their initial distribution of probability. The single-time correlations at all time steps and in each spatial grid cell are then determined from the Laplace–Fourier transforms of the discretized Klimontovich-like equation for the macroparticles and Maxwell’s equations for the fields, as computed by modern PIC codes. We recover the expressions for the electrostatic field and the plasma particle density fluctuation autocorrelation spectra as well as the kinetic equations describing the average evolution of PIC-simulated plasma particles, first derived by Langdon (1970b Proc. 4th Conf. Numerical Simulation of Plasmas ) using a test macroparticle approach perturbing a discretized Vlasovian plasma and then averaging the obtained physical quantity over the initial macroparticle velocity distribution. We generalize and extend these results to the modern algorithms in PIC codes using arbitrary macroparticle weights. Analytical estimates of statistical fluctuation amplitudes are derived as a function of the plasma simulation parameters, using the central limit theorem in the limit of a large number of macroparticles per cell. The theory is then used to analyze the ensemble averaging technique of PIC simulations where statistical averages are performed over ensembles of PIC simulations, modeling the same plasma physics problem but using different statistical realizations of the initial distribution functions of the macroparticles. This method is illustrated by linear Landau damping uncovering (from noise, which is usually considered numerical) the physical fluctuations driven by a single small amplitude electrostatic wave perturbing a PIC simulation plasma in equilibrium.

fluctuations correlations↗

DOE OSTI2024

Laser-stimulated photodetachment of electrons from the negatively charged dielectric substrates

The photon-stimulated emission of bulk electrons has been extensively studied for various types of materials, while the photodetachment of surplus surface electrons has not been fully explored. The photodetachment barrier energy is commonly defined by the surface electron affinity of material, which is typically less than the work function and more pronounced for non-conducting substrates and in environments with a continuous flux of electrons to the surface, such as in gas discharge plasmas. Herein, it is experimentally shown that the photodetachment yield of surplus electrons created by plasma-induced charging of non-conductive surfaces of dielectric materials depends on the initial surface charge density and do not correlate with the tabulated affinity values of these materials under gas discharge charging conditions. In cocnclusion, this result obtained using laser-stimulated photodetachment for fused silica, boron nitride, and alumina, is critically important for the understanding of charging and discharging dynamics, secondary electron emission, and photo emission effects affecting plasma–wall interactions relevant to surface and capacitively coupled discharges, dusty plasmas, electrostatic probe diagnostics, and applications for plasma processing of materials, plasma propulsion, and gas breakdown.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Dispersion and the speed-limited particle-in-cell algorithm

This paper discusses temporally continuous and discrete forms of the speed-limited particle-in-cell (SLPIC) method first treated by Werner et al. [Phys. Plasmas 25, 123512 (2018)]. The dispersion relation for a 1D1V electrostatic plasma whose fast particles are speed-limited is derived and analyzed. By examining the normal modes of this dispersion relation, we show that the imposed speed-limiting substantially reduces the frequency of fast electron plasma oscillations while preserving the correct physics of lower-frequency plasma dynamics (e.g. ion acoustic wave dispersion and damping). We then demonstrate how the timestep constraints of conventional electrostatic particle-in-cell methods are relaxed by the speed-limiting approach, thus enabling larger timesteps and faster simulations. Here, these results indicate that the SLPIC method is a fast, accurate, and powerful technique for modeling plasmas wherein electron kinetic behavior is nontrivial (such that a fluid/Boltzmann representation for electrons is inadequate) but evolution is on ion timescales.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Close-Out Project Description for Koepke's Dept of Energy grant DE-SC0021405

Objectives: To establish, for lab & space conditions, EM-IEDDI’s (electromagnetic shear-driven instability's) dispersion relation, unstable range, instability threshold, and mode characteristics, we need LAPD’s Alfven-wave-favorable electromagnetic-style conditions, including higher "beta" (0.001 < beta ≤ 0.3) and low-collisionality. Also, we attempted to intentionally launch or spontaneously destabilize compressional and shear Alfven waves in the strong, localized, perpendicular-velocity-shear region at the interface between coaxial plasmas (one plasma cylinder inside an outer, otherwise hollow, tube, each having a different, controllable, value of plasma electrostatic potential, i.e., “space” potential (not to be confused with the temperature-dependent “floating” potential of an object immersed in the plasma). Nonlinear wave-wave interactions between same-family (EM-IEDD or Alfven) and cross-family (EM-IEDD-with-Alfven) fluctuations were targeted for documentation over a range of spectral overlap. Although laboratory experiments were conducted, the following theoretical work was left unfinished: Analytical non-modal prediction Computational non-modal prediction Check to see if Mikhailenko’s theory formulation leads to his published graphs

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

An explicit scheme to enforce charge conservation in transient Particle-in-Cell simulations with Maxwell-Boltzmann electrons

The Maxwell-Boltzmann electron model is a very popular approximation of the electronic behavior in electrostatic plasmas, adopted in all cases where the electron dynamics is governed only by a balance between electric and pressure forces. In such a model, the electron dynamics is reduced to solving for a nonlinear Poisson problem together with an additional expression enforcing charge conservation. In this work we derive an expression of charge conservation which can be conveniently applied to explicit schemes to update the reference Boltzmann density, particularly useful when simulating plasma sheaths. The scheme can be equally applied to both steady-state and transient problems. Two examples are shown, a steady-state plasma sheath, and a radio-frequency magnetic presheath. Our proposed scheme allows to enforce global charge conservation locally in time, and can thus be applied to the simulation of transient phenomena.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Space-Time Structured Plasma Waves

Electrostatic waves play a critical role in nearly every branch of plasma physics from fusion to advanced accelerators, to astro, solar, and ionospheric physics. The properties of planar electrostatic waves are fully determined by the plasma conditions, such as density, temperature, ionization state, or details of the distribution functions. Furthermore we demonstrate that electrostatic wave packets structured with space-time correlations can have properties that are independent of the plasma conditions. For instance, an appropriately structured electrostatic wave packet can travel at any group velocity, even backward with respect to its phase fronts, while maintaining a localized energy density. These linear, propagation-invariant wave packets can be constructed with or without orbital angular momentum by superposing natural modes of the plasma and can be ponderomotively excited by space-time structured laser pulses like the flying focus.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

The integrated miniaturized electrostatic analyzer: A space plasma environment sensor

The integrated Miniaturized Electrostatic Analyzer (iMESA) was a satellite-based ionospheric sensor that operated on NASA’s Space Test Program Satellite (STPSat-3) from December 2013 to July 2019. The instrument’s scientific objective was to (1) measure the plasma density in low Earth orbit, (2) measure the plasma temperature in low Earth orbit, and (3) quantify the spacecraft potential with respect to the ambient plasma potential in the ionosphere. iMESA sampled the ionosphere every 10 s by measuring the ion current density through the ESA as a result of the motion of the spacecraft through the plasma. Current density spectra were transmitted to the ground where they were post-processed into ion density spectra and then analyzed numerically to determine the ion density, ion temperature, and spacecraft potential. This article discusses the instrument design and simulation, the determination of a geometric factor, and data processing procedures and evaluates the final data product with regard to the mission success criteria. Here, the ion density and ion temperature captured by the iMESA instrument are on the same order and range as the values predicted in the literature. The spacecraft potential was also quantified. The conclusion after the evaluation of the instrument’s data product is that the scientific mission is successful on all three points.

47 OTHER INSTRUMENTATION↗

DOE OSTI2024

Kinetic Plasma Simulation in the MOOSE Framework: Verification of Electrostatic Particle In Cell Capabilities

In magnetic confinement nuclear fusion reactors, the interaction between the plasma edge and plasma facing components is extremely important. At the plasma edge, a kinetic representation such as particle-in-cell (rather than a fluid representation) is required to accurately capture the plasma behavior. General purpose particle-in-cell plasma simulation capabilities have been developed in the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework. This new capability is a part of the development of a new MOOSE-based framework for modeling plasma facing components, the Fusion ENergy Integrated multiphys-X (FENIX) framework. In this work, the verification of foundational particle-in-cell capabilities in FENIX is presented. This new plasma simulation capability has three main components: moving particles in discrete steps on the finite element mesh, mapping charge density from the particle's location to the finite element mesh, and solving for the electrostatic potential based on the charge density mapped from particles to the mesh. In this paper, simple verification problems demonstrating each of these new capabilities are presented, and future work includes electromagnetic capabilities and Monte Carlo collisions with neutral gas particles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Eulerian simulations of electrostatic waves in plasmas with a single sign of charge

An Eulerian, numerical simulation is used to model the launching of plasma waves in a non-neutral plasma that is confined in a Penning–Malmberg trap. The waves are launched by applying an oscillating potential to an electrically isolated sector at one end of the conducting cylinder that bounds the confinement region and are received by another electrically isolated sector at the other end of the cylinder. The launching of both Trivelpiece–Gould waves and electron acoustic waves is investigated. Adopting a stratagem, the simulation captures essential features of the finite length plasma, while retaining the numerical advantages of a simulation employing periodic spatial boundary conditions. As a benchmark test of the simulation, the results for launched Trivelpiece–Gould waves of small amplitude are successfully compared to a linearized analytic solution for these fluctuations.

Physics↗

DOE OSTI2020

Momentum-exchange current drive by electrostatic waves in an unmagnetized collisionless plasma

For a planar electrostatic wave interacting with a single species in a collisionless plasma, momentum conservation implies current conservation. However, when multiple species interact with the wave, they can exchange momentum, leading to current drive. A simple, general formula for this driven current is derived. As examples, we show how currents can be driven for Langmuir waves in electron–positron–ion plasmas, and for ion-acoustic waves in electron–ion plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Electrostatic turbulence in EAST plasmas with internal transport barrier

Based on first-principles nonlinear gyrokinetic simulations, the electrostatic turbulence properties in the internal transport barrier (ITB) region of an Experimental Advanced Superconducting Tokamak discharge (#93890) are investigated. Specifically, ITBs with steep density and temperature gradients are located in the weakly negative magnetic shear region at the plasma center. In the linear stage, the growth rate and frequency of the ion temperature gradient (ITG) mode increase significantly due to resonant excitation by trapped electrons. That is, the resonance between trapped electrons and the ITG becomes strong due to the precession drift reversal of trapped electrons by the negative magnetic shear and Shafranov shift. Meanwhile, the trapped electron mode is stable in the ITB region due to only a very small fraction of electrons precessing in the direction of the electron diamagnetic drift. Nonlinear simulations show that, after considering the non-adiabatic effect of trapped electrons, the heat conductivity of ions and the turbulence intensity increase by at least a factor of 7 compared with the results only considering the adiabatic effect of electrons. The zonal charge density of trapped electrons can partially cancel that of ions, which weakens the intensity of the zonal flow, and consequently reduces the zonal flow regulation and enhances the turbulent transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Finite ion temperature effects on electrostatic instabilities in partially magnetized plasmas

In this paper, the effects of warm ions on electrostatic, kinetic instabilities observed in partially magnetized plasmas, such as the electron cyclotron drift instability (ECDI) and modified two-stream instability (MTSI), are investigated. The kinetic dispersion relation that accounts for warm, non-magnetized ions and warm, magnetized electrons is solved using a complex root-finding technique [A. C. Denig and K. Hara, Phys. Plasmas 30, 032108 (2023)]. It is observed that finite ion temperature introduces ion Landau damping effects, decreasing the growth rate of most unstable modes. The high-frequency, short-wavelength ECDI modes are more strongly affected by ion Landau damping than the low-frequency, long-wavelength MTSI modes. The solutions of the kinetic dispersion relation for MTSI are compared with the fluid dispersion relation of MTSI, which was originally proposed by McBride et al. [Phys. Fluids 15 2367–2383 (1972)]. The differences in the dispersion characteristics between the kinetic and fluid MTSI limits are assessed, illustrating the effects of electron and ion Landau damping on kinetic instabilities in partially magnetized plasmas.

Denig, A. C. (ORCID:0000000253339572)↗

DOE OSTI2021

1D EDIPIC

EDIPIC codes were developed for simulation of kinetic effects in low-temperature plasmas. They can be used for simulation of plasmas in electrostatic approximation and include many atomistic and plasma-surface interaction processes needed for simulations of low temperature plasmas. These codes include a detailed model of secondary electron emission induced by electrons and ions. The Monte-Carlo model of electron-neutral collisions and the Langevin model of Coulomb collisions for electrons are also implemented.

Sydorenko, Dmytro↗

DOE OSTI2021

EDIPIC-2D

Sydorenko, Dmytro↗