Search NASA⌕ Search

SEARCH · Search NASA

Results for “PLASMA SIMULATION”

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 1,513 records · Page 84
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Scaling of Magnetic Reconnection in Relativistic Collisionless Pair Plasmas

Using fully kinetic simulations, we study the scaling of the inflow speed of collisionless magnetic reconnection from the non-relativistic to ultra-relativistic limit. In the anti-parallel configuration, the inflow speed increases with the upstream magnetization parameter σ and approaches the light speed when σ>O(100), leading to an enhanced reconnection rate. In all regimes, the divergence of pressure tensor is the dominant term responsible for breaking the frozen-in condition at the x-line. The observed scaling agrees well with a simple model that accounts for the Lorentz contraction of the plasma passing through the diffusion region. The results demonstrate that the aspect ratio of the diffusion region remains ~0.1 in both the non-relativistic and relativistic limits.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Faraday effect in collisional magnetized plasmas

Faraday rotation is a valuable diagnostic tool for investigating laboratory and astrophysical plasmas, but it has almost exclusively been treated in the collisionless limit (despite the non-negligible role of Coulomb collisions in many laboratory plasmas). Here, we show that Faraday rotation can occur in collisional plasmas and that the usual effect is curiously altered by collisions. Namely, an initially linearly polarized light wave—propagating parallel to a uniform magnetic field—becomes elliptically polarized in the collisional plasma, and collisional absorption of the wave also occurs. Moreover, the ellipticity and rotation angle are quantifiably sensitive to the fidelity of the collisional transport coefficients. As we will demonstrate with particle-in-cell EPOCH simulations, these effects offer a stringent diagnostic tool for benchmarking multi-scale plasma simulation codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2007

Space Plasma Testing of High-Voltage Thin-Film Solar Arrays with Protective Coatings

This paper gives an overview of the space plasma test program for thin-film photovoltaics (TFPV) technologies developed at the Air Force Research Laboratory (AFRL). The main objective of this program is to simulate the effects of space plasma characteristic of LEO and MEO environments on TFPV. Two types of TFPV, amorphous silicon (a-Si) and copper-indium-gallium-diselenide (CIGS), coated with two types of thin-film, multifunctional coatings were used for these studies. This paper reports the results of the first phase of this program, namely the results of preliminary electrostatic charging, arcing, dielectric breakdown, and collection current measurements carried out with a series of TFPV exposed to simulated space plasma at the NASA Glenn Plasma Interaction Facility. The experimental data demonstrate that multifunctional coatings developed for this program provide effective protection against the plasma environment while minimizing impact on power generation performance. This effort is part of an ongoing program led by the Space Vehicles Directorate at the AFRL devoted to the development and space qualification of TFPV and their protective coatings.

Tlomak, Pawel↗
DOE OSTI2016

Kinetic modeling of E-to-H mode transition in inductively coupled hydrogen plasmas

Radio Frequency (RF) Inductively Coupled Plasmas (ICPs) are widely known for their two discharge modes, i.e., H-mode and E-mode, where the dynamics of the plasmas are completely different from each other. We have performed a kinetic simulation of a hydrogen plasma discharge in order to clarify the discharge mechanism and the E-to-H transition of the RF ICPs. The numerical simulation results, such as the time variations of spatial distribution of electron density and the power dissipated in the plasma, show the characteristic changes of the plasma dynamics due to E-to-H mode transition. Especially, the drastic change during the mode transition has been observed in the time evolution of the electron energy distribution function (EEDF). The EEDF deviates from a Maxwellian distribution before/after the transition and the deviation is more significant in the E-mode phase. These results indicate the importance of kinetic modeling for the physical understanding of E-to-H transition.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A moment-conserving discontinuous Galerkin representation of the relativistic Maxwellian distribution

Kinetic simulations of relativistic gases and plasmas are critical for understanding diverse astrophysical and terrestrial systems, but the accurate construction of the relativistic Maxwellian, the Maxwell–Jüttner distribution, on a discrete simulation grid is challenging. Difficulties arise from the finite velocity bounds of the domain, which may not capture the entire distribution function, as well as errors introduced by projecting the function onto a discrete grid. Here, we present a novel scheme for iteratively correcting the moments of the projected distribution applicable to all grid-based discretizations of the relativistic kinetic equation. In addition, we describe how to compute the needed nonlinear quantities, such as Lorentz boost factors, in a discontinuous Galerkin scheme through a combination of numerical quadrature and weak operations. The resulting method accurately captures the distribution function and ensures that the moments match the desired values to machine precision.

astrophysical plasmas↗
DOE OSTI2022

Three-dimensional inhomogeneity of electron-temperature-gradient turbulence in the edge of tokamak plasmas

Nonlinear multiscale gyrokinetic simulations of a Joint European Torus edge pedestal are used to show that electron-temperature-gradient (ETG) turbulence has a rich three-dimensional structure, varying strongly according to the local magnetic-field configuration. In the plane normal to the magnetic field, the steep pedestal electron temperature gradient gives rise to anisotropic turbulence with a radial (normal) wavelength much shorter than in the binormal direction. In the parallel direction, the location and parallel extent of the turbulence are determined by the variation in the magnetic drifts and finite-Larmor-radius (FLR) effects. The magnetic drift and FLR topographies have a perpendicular-wavelength dependence, which permits turbulence intensity maxima near the flux-surface top and bottom at longer binormal scales, but constrains turbulence to the outboard midplane at shorter electron-gyroradius binormal scales. Here our simulations show that long-wavelength ETG turbulence does not transport heat efficiently, and significantly decreases overall ETG transport—in our case by ~40%—through multiscale interactions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Simulation study of high-frequency energetic particle driven geodesic acoustic mode

High-frequency energetic particle driven geodesic acoustic modes (EGAM) observed in the large helical device plasmas are investigated using a hybrid simulation code for energetic particles and magnetohydrodynamics (MHD). Energetic particle inertia is incorporated in the MHD momentum equation for the simulation where the beam ion density is comparable to the bulk plasma density. Bump-on-tail type beam ion velocity distribution created by slowing down and charge exchange is considered. It is demonstrated that EGAMs have frequencies higher than the geodesic acoustic modes and the dependence on bulk plasma temperature is weak if (1) energetic particle density is comparable to the bulk plasma density and (2) charge exchange time (τ{sub cx}) is sufficiently shorter than the slowing down time (τ{sub s}) to create a bump-on-tail type distribution. The frequency of high-frequency EGAM rises as the energetic particle pressure increases under the condition of high energetic particle pressure. The frequency also increases as the energetic particle pitch angle distribution shifts to higher transit frequency. It is found that there are two kinds of particles resonant with EGAM: (1) trapped particles and (2) passing particles with transit frequency close to the mode frequency. The EGAMs investigated in this work are destabilized primarily by the passing particles whose transit frequencies are close to the EGAM frequency.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Effects of artificial collisions, filtering, and nonlocal closure approaches on Hermite-based Vlasov–Poisson simulations

Kinetic simulations of collisionless plasmas are computationally challenging due to phase-space mixing and filamentation, resulting in fine-scale velocity structures. This study compares three methods developed to reduce artifacts related to limited velocity resolution in Hermite-based Vlasov–Poisson simulations: artificial collisions, filtering, and nonlocal closure approaches. We evaluate each method's performance in approximating the linear kinetic response function and suppressing recurrence in linear and nonlinear regimes. Numerical simulations of Landau damping demonstrate that artificial collisions, particularly higher orders of the Lenard-Bernstein collisional operator, most effectively recover the correct damping rate across a range of wavenumbers. Moreover, Hou-Li filtering and nonlocal closures underdamp high wavenumber modes in linear simulations, and the Lenard-Bernstein collisional operator overdamps low wavenumber modes in both linear and nonlinear simulations. This study demonstrates that hypercollisions offer a robust approach to kinetic simulations, accurately capturing collisionless dynamics with limited velocity resolution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1984

High-temperature erosion of plasma-sprayed, yttria-stabilized zirconia in a simulated turbine environment

A series of rig calibration and high temperature tests simulating gas path seal erosion in turbine engines were performed at three impingement angles and at three downstream locations. Plasma sprayed, yttria stablized zirconia specimens were tested. Steady state erosion curves presented for 19 test specimens indicate a brittle type of material erosion despite scanning electron microscopy evidence of plastic deformation. Steady state erosion results were not sensitive to downstream location but were sensitive to impingement angle. At difference downstream locations specimen surface temperature varied from 1250 to 1600 C (2280 to 2900 F) and particle velocity varied from 260 to 320 m/s (850 to 1050 ft/s). The mass ratio of combustion products to erosive grit material was typically 240.

Hanschuh, R. F.↗
NASA NTRS1985

High-temperature erosion of plasma-sprayed, yttria-stabilized zirconia in a simulated turbine environment

A series of rig calibration and high temperature tests simulating gas path seal erosion in turbine engines were performed at three impingement angles and at three downstream locations. Plasma sprayed, yttria stabilized zirconia specimens were tested. Steady state erosion curves presented for 19 test specimens indicate a brittle type of material erosion despite scanning electron microscopy evidence of plastic deformation. Steady state erosion results were not sensitive to downstream location but were sensitive to impingement angle. At different downstream locations specimen surface temperature varied from 1250 to 1600 C (2280 to 2900 F) and particle velocity varied from 260 to 320 m/s (850 to 1050 ft/s). The mass ratio of combustion products to erosive grit material was typically 240.

Handschuh, R. F.↗
DOE OSTI2017

Stability of divertor detachment

The 2D simulations of edge plasma transport show that unlike some earlier publications, the impurity radiation loss per se does not cause the bifurcation-like transition to detached divertor regime. However, for the case where anomalous plasma transport is increasing with advancement to detachment, like it was recently observed experimentally, the transition to detachment exhibit the bifurcation-like character. Some other plausible reasons for similar bifurcation-like evolution to detachment are discussed. It is demonstrated that the current convective instability can be triggered in detached inner divertor plasma for the condition when outer divertor is still attached. This can explain the fluctuations of radiation loss observed recently experimentally for similar conditions. The self-sustained oscillations observed recently in numerical simulations and related to the interplay of the thermal force effects in impurity transport and impurity radiation loss are further investigated. It is shown that for some conditions these oscillations are ubiquitous, since no stable solutions possible.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

High-fidelity kinetic modeling of instabilities and gyromotion physics in nonuniform low-beta plasmas

A fourth-order accurate continuum kinetic Vlasov solver and a systematic method for constructing customizable kinetic equilibria are demonstrated to be powerful tools for the study of nonuniform collisionless low-beta plasmas. The noise-free methodology is applied to investigate two gradient-driven instabilities in 4D (x,y,vx,vy) phase space: the Kelvin–Helmholtz instability and the lower hybrid drift instability. Nonuniform two-species configurations where ion gyroradii are comparable to gradient scale lengths are explored. The approach sheds light on the evolution of the pressure tensor in Kelvin–Helmholtz instabilities and demonstrates that the associated stress tensor deviates significantly from the gyroviscous stress tensor. Even at high magnetization, first-order approximations to finite-gyromotion physics are shown to be inadequate for the Kelvin–Helmholtz instability, as shear scales evolve to become on par with gyromotion scales. The methodology facilitates exploring transport and energy partitioning properties associated with lower hybrid drift instabilities in low-beta plasma configurations. Distribution function features are captured in detail, including the formation of local extrema in the vicinity of particle-wave resonances. Finally, the approach enables detailed targeted investigations and advances kinetic simulation capability for plasmas in which gyromotion plays an important role.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
NASA NTRS1987

A comparison of the Mercury and earth magnetospheres - Electron measurements and substorm time scales

The present search for similarities between earth and Mercury plasma electron distribution and large-scale dynamics notes that both spectral shapes are similar to a kappa-distribution. A model distribution of this type which incorporates convective flow is used to simulate the observed plasma electron spectral variations near the Mariner 10-Mercury 1 A event; convection appears to be stronger before, rather than during, the A event, in contradiction to the Baker (1986) convective injection model for Mercury's two relativistic electron flux enhancements. Mercury's postmidnight energetic electron B and B-prime events seem to be multiple onsets in the course of a substorm.

Christon, S. P.↗
NASA NTRSDate not supplied

AI-Enhanced Computational Tools for Entry Systems Modeling

To advance the understanding of complex atmospheric entry phenomena, NASA’s Entry Systems Modeling (ESM) team [1] has developed high-fidelity computational tools addressing multiscale challenges, from material microstructures to full-scale heatshield response. This abstract highlights a subset of ESM tools, focusing on AI integration to enhance workflows and predictive modeling. - PuMA [2] computes effective material properties from high-resolution micro-CT scans, supporting TPS analysis for NASA missions. - TomoSAM [3] automates 3D tomography dataset segmentation for PuMA using the Segment Anything Model, reducing manual effort and improving accuracy. - PATO [4] models porous reactive materials under extreme conditions, with advancements such as unified solvers, mechanical erosion, and TPS coatings for NASA missions. - arcjetCV [5] employs deep learning to analyze arc jet test footage, measuring recession rates, shape changes, and shock standoff distances, bridging simulations, and experiments to reveal TPS ablation behavior. - ARCHeS [6] simulates arc heater plasma flows, modeling turbulence, radiation, and electromagnetic interactions to optimize arc heater performance, validate TPS under extreme conditions, and serve as a foundation for developing digital twins of arc heater facilities. - SPARTA [7] simulates rarefied hypersonic flows and gas-surface interactions for planetary entry missions, leveraging GPU architectures for scalable and efficient aerothermal and ablation analyses. AI-driven solutions, such as deep learning segmentation, have streamlined workflows in ESM tools and still hold significant potential to further accelerate processes and enhance automation in entry systems modeling. [1] Haskins, J.B. (2023), [2] Ferguson, J.C. (2018), [3] Meurisse, J.B.E. (2018), [4] Semeraro, F. (2023), [5] Quintart, A. (2024) [6] Meurisse, J.B.E. (2022), [7] Plimpton, S.J. (2019)

Predictive Modeling↗
DOE OSTI2014

How the “main condition” of phase stability can explain the effect of the velocity deviation of secondary electrons in DC-biased single-sided multipactors

In this work, a “main condition” for phase stability has been employed to investigate the effects of the velocity deviation of the electrons in DC-biased single-sided multipactors (MPs). In a previous study [M. Mostajeran, Phys. Plasmas 21, 053108 (2014)], a stability equation was derived, where the secondary electron was assumed to have zero initial velocity and the phase deviation from the resonant phase was considered. In this work, both deviations in phase and velocity from the resonant condition are taken into account, assuming nonzero initial velocity for the secondary electrons. Using the main condition for stability, it is shown that MP discharge can rise in situations, where large velocity deviations from initial velocity and large phase deviations from resonant phase exist. This is contrary to what can be predicted on the basis of the “simple stability condition.” This result is further confirmed by numerical simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1993

A Search of Large Amplitude Waves in Slow Shock in the Distant Tail and an Examination of Wave Turbulence in the Plasma Sheet Boundary Layer

Two fluid simulations predict that switch-off slow shocks should have large amplitude left-hand wave trains downstream of the shock. Recent simulations show that this feature depends on specific plasma parameters, such as Much number and ratio of electron-to-ion temperature. In this study we carefully examine all distant tail crossings using ISEE-3 field and plasma data. We seek large amplitude rotational waves and compare to shock parameters. We also investigate right-hand waves associated with ion beams in the plasma sheet boundary layer. We examine the role of those back-streaming ions from slow shocks in wave generation. Wave mode cascade and parametric decay in these regions also will be studied.

Ho, C. M.↗