Search NASA⌕ Search

SEARCH · Search NASA

Results for “PLASMA DYNAMICS”

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 937 records · Page 52
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Simulations of tokamak edge plasma turbulent fluctuations based on a minimal 3D model

In this study, a new simulation model for tokamak boundary plasma, SOLT3D, is implemented in the BOUT++ framework. The simulation model includes a set of dynamic equations describing collisional boundary plasma and neutral gas in the tokamak scrape-off layer and divertor region. The model is verified against standard linear plasma instabilities and available nonlinear results. For L-mode like conditions, SOLT3D reproduces characteristics of boundary plasma turbulent fluctuations that are similar to published experimental data, in terms of the amplitude and spatial dependence of the fluctuations. It also reproduces realistic plasma fluxes on material surfaces and expected Bohm-like effective radial transport. Plasma fluctuations inferred from the simulations imply inevitably a significant level of intrinsic 'noise' for axisymmetric tokamak plasma transport modeling, introducing errors on the order of unity. In particular, the toroidally averaged atomic rates below 5–10 eV are strongly modified by turbulent plasma fluctuations, which should significantly affect the standard axisymmetric modeling of the tokamak edge plasma and divertor.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Identification of the weak-to-strong transition in Alfvénic turbulence from space plasma

Plasma turbulence is a ubiquitous dynamical process that transfers energy across many spatial and temporal scales in astrophysical and space plasma systems. Although the theory of anisotropic magnetohydrodynamic (MHD) turbulence has successfully described natural phenomena, its core prediction of an Alfvénic transition from weak to strong MHD turbulence when energy cascades from large to small scales has not been observationally confirmed. Here we report evidence for the Alfvénic weak-to-strong transition in small-amplitude, turbulent MHD fluctuations in Earth’s magnetosheath using data from four Cluster spacecraft. Our observations demonstrate the universal existence of strong turbulence accompanied by weak turbulent fluctuations on large scales. Moreover, we find that the nonlinear interactions of MHD turbulence are crucial to the energy cascade, as they broaden the cascade directions and fluctuating frequencies. The observed connection between weak and strong MHD turbulence systems may be present in many astrophysical environments, such as star formation, energetic particle transport, turbulent dynamos, and solar corona or solar wind heating.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

The fusion code XGC: Enabling kinetic study of multi-scale edge turbulent transport in ITER [Book Chapter]

The XGC fusion gyrokinetic code combines state-of-the-art, portable computational and algorithmic technologies to enable complicated multiscale simulations of turbulence and transport dynamics in ITER edge plasma on the largest US open-science computer, the CRAY XK7 Titan, at its maximal heterogeneous capability, which have not been possible before due to a factor of over 10 shortage in the time-to-solution for less than 5 days of wall-clock time for one physics case. Frontier techniques such as nested OpenMP parallelism, adaptive parallel I/O, staging I/O and data reduction using dynamic and asynchronous applications interactions, dynamic repartitioning for balancing computational work in pushing particles and in grid related work, scalable and accurate discretization algorithms for non-linear Coulomb collisions, and communication-avoiding subcycling technology for pushing particles on both CPUs and GPUs are also utilized to dramatically improve the scalability and time-to-solution, hence enabling the difficult kinetic ITER edge simulation on a present-day leadership class computer.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

PLASMA INSTABILITIES IN THE CONTEXT OF CURRENT HELIUM SEDIMENTATION MODELS: DYNAMICAL IMPLICATIONS FOR THE ICM IN GALAXY CLUSTERS

Understanding whether Helium can sediment to the core of galaxy clusters is important for a number of problems in cosmology and astrophysics. All current models addressing this question are one-dimensional and do not account for the fact that magnetic fields can effectively channel ions and electrons, leading to anisotropic transport of momentum, heat, and particle diffusion in the weakly collisional intracluster medium (ICM). This anisotropy can lead to a wide variety of instabilities, which could be relevant for understanding the dynamics of heterogeneous media. In this paper, we consider the radial temperature and composition profiles as obtained from a state-of-the-art Helium sedimentation model and analyze its stability properties. We find that the associated radial profiles are unstable to different kinds of instabilities depending on the magnetic field orientation at all radii. The fastest growing modes are usually related to generalizations of the magnetothermal instability (MTI) and the heat-flux-driven buoyancy instability which operate in heterogeneous media. We find that the effect of sedimentation is to increase (decrease) the predicted growth rates in the inner (outer) cluster region. The unstable modes grow quickly compared to the sedimentation timescale. This suggests that the composition gradients as inferred from sedimentation models, which do not fully account for the anisotropic character of the weakly collisional environment, might not be very robust. Our results emphasize the subtleties involved in understanding the gas dynamics of the ICM and argue for the need of a comprehensive approach to address the issue of Helium sedimentation beyond current models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Probing the plasma near high power wave launchers in fusion devices for static and dynamic electric fields (invited)

An exploratory study was carried out in the long-pulse tokamak Tore Supra, to determine if electric fields in the plasma around high-power, RF wave launchers could be measured with non-intrusive, passive, optical emission spectroscopy. The focus was in particular on the use of the external electric field Stark effect. The feasibility was found to be strongly dependent on the spatial extent of the electric fields and overlap between regions of strong (>∼1 kV/cm) electric fields and regions of plasma particle recycling and plasma-induced, spectral line emission. Most amenable to the measurement was the RF electric field in edge plasma, in front of a lower hybrid heating and current drive launcher. Electric field strengths and direction, derived from fitting the acquired spectra to a model including time-dependent Stark effect and the tokamak-range magnetic field Zeeman-effect, were found to be in good agreement with full-wave modeling of the observed launcher.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2016

A novel technique for real-time estimation of edge pedestal density gradients via reflectometer time delay data

A new model-based technique for fast estimation of the pedestal electron density gradient has been developed. The technique uses ordinary mode polarization profile reflectometer time delay data and does not require direct profile inversion. Because of its simple data processing, the technique can be readily implemented via a Field-Programmable Gate Array, so as to provide a real-time density gradient estimate, suitable for use in plasma control systems such as envisioned for ITER, and possibly for DIII-D and Experimental Advanced Superconducting Tokamak. The method is based on a simple edge plasma model with a linear pedestal density gradient and low scrape-off-layer density. By measuring reflectometer time delays for three adjacent frequencies, the pedestal density gradient can be estimated analytically via the new approach. Using existing DIII-D profile reflectometer data, the estimated density gradients obtained from the new technique are found to be in good agreement with the actual density gradients for a number of dynamic DIII-D plasma conditions.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2018

New helium bubble growth mode at a symmetric grain-boundary in tungsten: accelerated molecular dynamics study

This work, with an emphasis on helium irradiation rates appropriate for fusion-plasma conditions, advances the understanding of helium evolution at grain boundaries in W, an important consideration in the understanding of W as a plasma-facing component. Using accelerated molecular dynamics, helium bubble nucleation and growth at a symmetric Σ5[100](310) tilt grain-boundary in W is studied. The simulations reveal that the growth mode associated with bubble growth at the grain-boundary leads to a suppression of the helium supply to the bubble and hence to arrested growth. Such an unconventional bubble growth mode may dominate in materials with a high density of sinks.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Threshold for plasma phase transition of aluminum single crystal induced by hypervelocity impact

Molecular dynamics method is used to study the threshold for plasma phase transition of aluminum single crystal induced by hypervelocity impact. Two effective simulation methods, piston-driven method and multi-scale shock technique, are used to simulate the shock wave. The simulation results from the two methods agree well with the experimental data, indicating that the shock wave velocity is linearly dependent on the particle velocity. The atom is considered to be ionized if the increase of its internal energy is larger than the first ionization energy. The critical impact velocity for plasma phase transition is about 13.0 km/s, corresponding to the threshold of pressure and temperature which is about 220 GPa and 11.0 × 10{sup 3 }K on the shock Hugoniot, respectively.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Customizable two-species kinetic equilibria for nonuniform low-beta plasmas

Here, two methods are developed for constructing self-consistent two-species kinetic equilibria for low-beta electrostatic plasmas, in which the magnetic field is uniform and fixed in time. The first method uses Taylor series approximations to construct distribution functions that can be specified analytically. The second method relies on numerically solving a nonlinear ordinary differential equation and produces exact—to numerical precision—equilibria. In both methods, the equilibrium distribution functions for ions and electrons are expressed in terms of constants of motion and satisfy the steady-state Vlasov-Poisson equation system. Provided that the ion drift speed does not exceed the ion thermal speed, the equilibria can be specified with customizable density and electrostatic potential profiles. The methods can thereby be tailored to different applications and are successfully applied to construct kinetic equilibria for cross-field plasmas with sheared flows, large density variations, and different levels of magnetization. The equilibria are used to initialize fourth-order finite-volume Vlasov-Poisson simulations in (x, v x , v y ) coordinates and the associated temporal evolution is used to assess the accuracy of each method. Lastly, the low-amplitude deviations observed in these simulations demonstrate that the kinetic equilibria are robust and that they provide a valuable means of studying the dynamics of nonuniform magnetized plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Detecting shearless phase-space transport barriers in global gyrokinetic turbulence simulations with test particle map models

In magnetically confined fusion plasmas, the role played by zonal 𝐸 x 𝐵 flow shear layers in the suppression of turbulent transport is relatively well understood. However, less is understood about the role played by the weak shear regions that arise in the non-monotonic radial electric field profiles often associated with these shear layers. In electrostatic simulations from the global total-𝑓 gyrokinetic particle-in-cell code XGC, we demonstrate how shearless regions with non-zero flow curvature form zonal ‘jets’ that, in conjunction with neighbouring regions of shear, can act as robust barriers to particle transport and turbulence spreading. By isolating quasi-coherent fluctuations radially localised to the zonal jets, we construct a map model for the Lagrangian dynamics of gyrokinetic test particles in the presence of drift waves. We identify the presence of shearless invariant tori in this model and verify that these tori act as partial phase-space transport barriers in the simulations. We also demonstrate how avalanches impinging on these shearless tori cause eddy detachment events that form ‘cold/warm core ring’ structures analogous to those found in oceanic jets, facilitating transport across the barriers without destroying them completely. We discuss how shearless tori may generically arise from tertiary instabilities or other types of discrete eigenmodes, suggesting their potential relevance to broader classes of turbulent fluctuations.

fusion plasma↗
DOE OSTI2017

Thermodynamic and Kinetic Properties of Shocks in Two-Dimensional Yukawa Systems [Thermodynamic and Kinetic Properties of Shocks in 2D Yukawa Systems]

Particle-level simulations of shocked plasmas are carried out to examine kinetic properties not captured by hydrodynamic models. In particular, molecular dynamics simulations of 2D Yukawa plasmas with variable couplings and screening lengths are used to examine shock features unique to plasmas, including the presence of dispersive shock structures for weak shocks. A phase-space analysis reveals several kinetic properties, including anisotropic velocity distributions, non-Maxwellian tails, and the presence of fast particles ahead of the shock, even for moderately low Mach numbers. As a result, we also examine the thermodynamics (Rankine-Hugoniot relations) of recent experiments and find no anomalies in their equations of state.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Ion fluxes and memory effects in an Ar–O 2 modulated radiofrequency-driven atmospheric pressure plasma jet

In this work we used molecular beam mass spectrometry to characterize ion fluxes produced by a modulated radiofrequency (RF)-driven atmospheric pressure plasma jet operating in a homogeneous gas environment (Ar + 1% O 2 ). The influence of the RF modulation frequency (100 Hz to 20 kHz) on the ion fluxes was investigated by time-resolved measurements, and the lifetimes of the dominant ions, O 2 + , NO + , O − , O 2 − and O 3 − , were found to be 28 ± 2, 117 ± 8, 7.3 ± 0.4, 17 ± 1 and 23 ± 2 μs, respectively. The absolute ion densities in the near afterglow region were found to be of the order of 1011 cm−3. Significant differences in the dynamics of the positive and negative ions were found; these are explained by the presence of electrons at similar densities in the afterglow produced by electron detachment reactions from negative ions due to the large concentrations of atomic oxygen and singlet delta oxygen. Transitions in ion flux dynamics for different modulation frequencies and at the startup of the plasma were analyzed together with intensified charge coupled device images recording the plasma propagation in order to assess the dynamics of plasma plume propagation and how it is impacted by ‘memory effects’. Quantitative measurements of the ion densities enabling these memory effects are reported. The results in this paper highlight the tremendous impact of memory effects on plasma propagation and the corresponding pre-ionization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Global flow structure and exact formal transseries of the Gubser flow in kinetic theory

In this work we introduce the generic conditions for the existence of a non-equilibrium attractor that is an invariant manifold determined by the long-wavelength modes of the physical system. We investigate the topological properties of the global flow structure of the Gubser flow for the Israel-Stewart theory and a kinetic model for the Boltzmann equation by employing Morse-Smale theory. We present a complete classification of the invariant submanifolds of the flow and determine all the possible flow lines connecting any pair of UV/IR fixed points. The formal transseries solutions to the Gubser dynamical system around the early-time (UV) and late-time (IR) fixed points are constructed and analyzed. It is proven that these solutions are purely perturbative (or power-law asymptotic) series with a finite radius of convergence. Based on these analyses, we find that Gubser-like expanding kinetic systems do not hydrodynamize owing to the failure of the hydrodynamization process which heavily relies on the classification of (non)hydrodynamic modes in the IR regime. This is in contrast to longitudinal boost-invariant plasmas where the asymptotic dynamics is described by a few terms of the hydrodynamic gradient expansion. We finally compare our results for both Bjorken and Gubser conformal kinetic models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Recent studies on nanosecond-timescale pressurized gas discharges

The results of recent experimental and numerical studies of nanosecond high-voltage discharges in pressurized gases are reviewed. The discharges were ignited in a diode filled by different gases within a wide range of pressures by an applied pulsed voltage or by a laser pulse in the gas-filled charged resonant microwave cavity. Fast-framing imaging of light emission, optical emission spectroscopy, X-ray foil spectrometry and coherent anti-Stokes Raman scattering were used to study temporal and spatial evolution of the discharge plasma density and temperature, energy distribution function of runaway electrons and dynamics of the electric field in the plasma channel. The results obtained allow a deeper understanding of discharge dynamical properties in the nanosecond timescale, which is important for various applications of these types of discharges in pressurized gases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

The Dynamics of Very High Alfvén Mach Number Shocks in Space Plasmas

Astrophysical shocks, such as planetary bow shocks or supernova remnant shocks, are often in the high or very-high Mach number regime, and the structure of such shocks is crucial for understanding particle acceleration and plasma heating, as well inherently interesting. Recent magnetic field observations at Saturn’s bow shock, for Alfvén Mach numbers greater than about 25, have provided evidence for periodic non-stationarity, although the details of the ion- and electron-scale processes remain unclear due to limited plasma data. High-resolution, multi-spacecraft data are available for the terrestrial bow shock, but here the very high Mach number regime is only attained on extremely rare occasions. Here we present magnetic field and particle data from three such quasi-perpendicular shock crossings observed by the four-spacecraft Cluster mission. Although both ion reflection and the shock profile are modulated at the upstream ion gyroperiod timescale, the dominant wave growth in the foot takes place at sub-proton length scales and is consistent with being driven by the ion Weibel instability. The observed large-scale behavior depends strongly on cross-scale coupling between ion and electron processes, with ion reflection never fully suppressed, and this suggests a model of the shock dynamics that is in conflict with previous models of non-stationarity. Thus, the observations offer insight into the conditions prevalent in many inaccessible astrophysical environments, and provide important constraints for acceleration processes at such shocks.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

When should PIC simulations be applied to atmospheric pressure plasmas? Impact of correlation heating

Abstract Molecular dynamics simulations are used to test when the particle-in-cell (PIC) method applies to atmospheric pressure plasmas. It is found that PIC applies only when the plasma density and macroparticle weight are sufficiently small because of two effects associated with correlation heating. The first is the physical effect of disorder-induced heating (DIH). This occurs if the plasma density is large enough that a species (typically ions) is strongly correlated in the sense that the Coulomb coupling parameter exceeds one. In this situation, DIH causes ions to rapidly heat following ionization. PIC is not well suited to capture DIH because doing so requires using a macroparticle weight of one and a grid that well resolves the physical interparticle spacing. These criteria render PIC intractable for macroscale domains. The second effect is a numerical error due to artificial correlation heating (ACH). ACH is like DIH in that it is caused by the Coulomb repulsion between particles, but differs in that it is a numerical effect caused by a macroparticle weight larger than one. Like DIH, it is associated with strong correlations. However, here the macroparticle coupling strength is found to scale as Γ w 2 / 3 , where Γ is the physical coupling strength and w is the macroparticle weight. So even if the physical coupling strength of a species is small, as is expected for electrons in atmospheric pressure plasmas, a sufficiently large macroparticle weight can cause the macroparticles to be strongly coupled and therefore heat due to ACH. Furthermore, it is shown that simulations in reduced dimensions exacerbate these issues.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Anomalous transport of multi-species edge plasma with the generalized Hasegawa–Wakatani model and the FLR effects

Anomalous transport of multi-species plasma is considered with the generalized Hasegawa–Wakatani model [A. R. Knyazev and S. I. Krasheninnikov, Phys. Plasmas 31, 012502 (2024)] further extended to incorporate the Finite Larmor Radius (FLR) effects. By introducing the “associated” enstrophy, it is shown that with no FLR effects (where anomalous transport of all ion species is described as a transport of passive scalars in the turbulent fields of the electrostatic potential and electron density fluctuations) the fluctuating densities of ion species converge to the state where they are linearly proportional to electron density and vorticity fluctuations, which confirm previous numerical findings of [A. R. Knyazev and S. I. Krasheninnikov, Phys. Plasmas 31, 012502 (2024)]. However, in contrast to the “cold” ion approximation, with the FLR effects included, both the plasma turbulence and the dynamics of all ion species become interconnected. Therefore, for simplicity, the FLR effects in this work were considered only for a small “trace” impurity fraction. It is found that for light (neon) “trace” impurity, the FLR effects reduce both anomalous flux and density fluctuations. However, for heavy (tungsten) “trace” impurity, the FLR effects exhibit non-monotonic impact on anomalous transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Avalanches driven by pressure gradients in a magnetized plasma

The results are presented for a basic heat transport experiment involving an off-axis heat source in which avalanche events occur. The configuration consists of a long, hollow, cylindrical region of elevated electron temperature embedded in a colder plasma, and far from the device walls [Van Compernolle et al. Phys. Rev. E 91, 031102(R) (2015)]. The avalanche events are identified as sudden rearrangements of the pressure profile following the growth of fluctuations from ambient noise. The intermittent collapses of the plasma pressure profile are associated with unstable drift-Alfvén waves and exhibit both radial and poloidal dynamics. After each collapse, the plasma enters a quiescent phase in which the pressure profile slowly recovers and steepens until a threshold is exceeded, and the process repeats. The use of reference probes as time markers allows for the visualization of the 2D spatio-temporal evolution of the avalanche events. Avalanches are observed only for a limited combination of heating powers and magnetic fields. At higher heating powers, the system transits from the avalanche regime into a regime dominated by sustained drift-Alfvén wave activity. Here, this manuscript focuses on new results that illustrate the individual contributions to the avalanche process from density and temperature gradients in the presence of zero-order, sheared flows.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗