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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Global Gyrokinetic Simulations of Isotope Effects under Ambipolar Electric Fields and Advances Toward Whole-Volume Modeling

We review global gyrokinetic simulation studies on plasma transport in the Large Helical Device using XGC-S. XGC-S is an extended version of X-point Gyrokinetic Code for stellarators and has been progressively verified throughout the code development process. Verification tests of neoclassical transport successfully demonstrate the generation of an ambipolar electric field due to ripple-trapped particles. We perform quasi-linear analyses of the ion temperature gradient mode under the influence of the ambipolar electric field. The results reveal that the ambipolar electric field and the heavy hydrogen component in mixed isotope plasmas can lead to the favorable isotope effect observed in recent deuterium experiments. We also present recent efforts in code development toward whole-volume simulations, including the helical divertor region. A mesh generation scheme based on field-line tracing and the construction of curved surfaces perpendicular to the magnetic field would be promising for global field calculations in the whole-volume simulations.

Basic Plasma Phenomena and Gas Discharges↗
DOE OSTI2017

The accurate particle tracer code

The Accurate Particle Tracer (APT) code is designed for systematic large-scale applications of geometric algorithms for particle dynamical simulations. Based on a large variety of advanced geometric algorithms, APT possesses long-term numerical accuracy and stability, which are critical for solving multi-scale and nonlinear problems. To provide a flexible and convenient I/O interface, the libraries of Lua and Hdf5 are used. Following a three-step procedure, users can efficiently extend the libraries of electromagnetic configurations, external non-electromagnetic forces, particle pushers, and initialization approaches by use of the extendible module. APT has been used in simulations of key physical problems, such as runaway electrons in tokamaks and energetic particles in Van Allen belt. As an important realization, the APT-SW version has been successfully distributed on the world’s fastest computer, the Sunway TaihuLight supercomputer, by supporting master–slave architecture of Sunway many-core processors. Here, based on large-scale simulations of a runaway beam under parameters of the ITER tokamak, it is revealed that the magnetic ripple field can disperse the pitch-angle distribution significantly and improve the confinement of energetic runaway beam on the same time.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Early- and late-time dynamics of laser-produced plasmas by combining emission and absorption spectroscopy

Optical emission spectroscopy and laser absorption spectroscopy are two highly complementary techniques for diagnosing plasmas, specifically in the determination of important parameters such as temperature and density. In this work, we characterized a laser-produced plasma by combining emission and tunable absorption spectroscopy. Plasmas are produced by focusing pulses from a frequency quadrupled Nd:YAG laser onto a glass target containing a minor concentration of aluminum in an ambient atmosphere of ~100 Torr N 2 . Here, emission spectroscopy is used to investigate the line behavior of Al and Ca, electron density, and excitation temperature at times = 5 μs, while absorption spectroscopy is used to investigate the line behavior of Al and Ca, species density, and kinetic temperature at times = 5-100 μs. A comparison is made between the excitation temperature measured from emission using the Boltzmann plot method and the kinetic temperature measured from absorption; the implications of this comparison on local thermodynamic equilibrium and line-of-sight averaging are discussed. By combining the two methods of optical emission and laser absorption spectroscopy, a comprehensive characterization of the plasma over the entire lifetime of its evolution is possible, leading to a greater understanding of the plasma behavior including particle condensation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Electromagnetic zonal flow residual responses

We report the collisionless axisymmetric zonal flow residual calculation for a tokamak plasma is generalized to include electromagnetic perturbations. We formulate and solve the complete initial value zonal flow problem by retaining the fully self-consistent axisymmetric spatial perturbations in the electric and magnetic fields. Simple expressions for the electrostatic, shear and compressional magnetic residual responses are derived that provide a fully electromagnetic test of the zonal flow residual in gyrokinetic codes. Unlike the electrostatic potential, the parallel vector potential and the parallel magnetic field perturbations need not relax to flux functions for all possible initial conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Practical gyrokinetics

We report thousands of gyrokinetic papers have been published since the introduction of the gyrokinetic change of variables. The intent here is not to review the field, but rather the goal is to present a tutorial providing insight into why gyrokinetics is an appropriate description of turbulent transport. The focus is on turbulent transport in axisymmetric tokamaks, but many of the ideas and techniques are applicable to stellarators and other magnetic fusion devices with nested surfaces of constant pressure. Besides the origins of gyrokinetics and recent insights, gyrokinetic orderings and gyrokinetic variables are summarized. Then a compact, but careful, derivation of the simplest electrostatic gyrokinetic equation for tokamaks is presented, along with a brief mention of gyrokinetics for stellarators. The advantages of assuming scale separation between the finer spatial scales and faster time variation of the turbulence and the global behaviour and slow temporal evolution of the background are stressed. Moreover, the procedure for removing the adiabatic or Maxwell–Boltzmann response is emphasized. Scale separation allows the near Maxwellian behaviour of the background and the rapid variation of the turbulence to be described by separate, local gyrokinetic equations. The turbulent fluctuations are found by solving the nonlinear gyrokinetic equation for the non-adiabatic portion of the fluctuating distribution function. Also, generalizations of the gyrokinetic equation so that it retains electromagnetic and flow modifications are considered in detail along with some symmetry properties. In addition, ambipolarity, particle transport and heat transport are addressed, along with a discussion of the complications associated with describing momentum transport and profile evolution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Improved accuracy in degenerate variational integrators for guiding centre and magnetic field line flow

First-order-accurate degenerate variational integration (DVI) was introduced in Ellison et al. ( Phys. Plasmas , vol. 25, 2018, 052502) for systems with a degenerate Lagrangian, i.e. one in which the velocity-space Hessian is singular. In this paper we introduce second-order-accurate DVI schemes, both with and without non-uniform time stepping. We show that it is not in general possible to construct a second-order scheme with a preserved two-form by composing a first-order scheme with its adjoint, and discuss the conditions under which such a composition is possible. We build two classes of second-order-accurate DVI schemes. We test these second-order schemes numerically on two systems having non-canonical variables, namely the magnetic field line and guiding centre systems. Variational integration for Hamiltonian systems with non-uniform time steps, in terms of an extended phase space Hamiltonian, is generalized to non-canonical variables. It is shown that preservation of proper degeneracy leads to single-step (one-step) methods without parasitic modes, i.e. to non-uniform time step DVIs. This extension applies to second-order-accurate as well as first-order schemes, and can be applied to adapt the time stepping to an error estimate.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Evolution of an arched magnetized laboratory plasma in a sheared magnetic field

Arched magnetized structures are a common occurrence in space and laboratory plasmas. Results from a laboratory experiment on spatio-temporal evolution of an arched magnetized plasma ( $\beta \approx 10^{-3}$ , Lundquist number $\approx 10^{4}$ , plasma radius/ion gyroradius $\approx 20$ ) in a sheared magnetic configuration are presented. The experiment is designed to model conditions relevant to the formation and destabilization of similar structures in the solar atmosphere. The magnitude of a nearly horizontal overlying magnetic field was varied to study its effects on the writhe and twist of the arched plasma. In addition, the direction of the guiding magnetic field along the arch was varied to investigate its role in the formation of either forward- or reverse-S shaped plasma structures. The electrical current in the arched plasma was well below the current required to make it kink unstable. A significant increase in the writhe of the arched plasma was observed with larger magnitudes of overlying magnetic field. A forward-S shaped arched plasma was observed for a guiding magnetic field oriented nearly antiparallel to the initial arched plasma current, while the parallel orientation yielded the reverse-S shaped arched plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Neural network reconstruction of the DIII-D tokamak plasma boundary using a reduced set of diagnostics

This study investigates the feasibility of reconstructing the last closed flux surface in the DIII-D tokamak using neural network models trained on reduced input feature sets, addressing an ill-posed task. Two models are compared: one trained solely on coil currents and another incorporating coil currents, plasma current and loop voltage. The model trained exclusively on coil currents achieved a mean point displacement of $0.04$ m on a held-out test set, while the inclusion of plasma current and loop voltage reduced the error to $0.03$ m. This comparison highlights the trade-offs between input feature complexity and reconstruction accuracy, demonstrating the potential of machine learning algorithms to perform effectively in data-limited environments, such as those expected in fusion power plants due to diagnostic constraints imposed by the presence of blankets and shielding.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Full-polarization millimeter wavelength variability of Sagittarius A * during the 2018 EHT campaign

Context. Sagittarius A* (Sgr A*), the supermassive black hole at the center of the Milky Way, provides a unique laboratory to study accretion dynamics and plasma processes near the event horizon. Aims. We investigated the variability and polarization properties of Sgr A* using ALMA observations during the 2018 Event Horizon Telescope campaign. Methods. We analyzed high-cadence full-polarization light curves from ALMA at millimeter wavelengths, performed time-series analysis, and investigated the temporal behavior during an X-ray flare observed by Chandra on 2018 April 24. The variability characteristics are compared with expectations from standard accretion flow models. Results. We find low variability in total intensity (σ/μ < 10%), but significantly higher variability in linear and circular polarization (∼30% and ∼50%, respectively). A time-series analysis reveals red-noise variability, with power spectral densities between −2 and −3 across all Stokes parameters. Polarized intensity shows stable intra-day timescales, while total intensity exhibits more variable timescales, suggesting distinct emission regions, with polarization likely arising from a coherent structure. On April 24, a statistically significant inter-band delay in polarized intensity coincides with a near-simultaneous X-ray and millimeter peak that deviates from the typical delayed flare scenario. This event also features enhanced millimeter variability and coherent polarization loop evolution. The observed simultaneity challenges standard models of transient synchrotron emission with cooling delays, favoring instead a scenario of continuous energy injection in an optically thin region. Conclusions. Our results offer new constraints on the physical mechanisms driving variability in Sgr A*, and provide key observational input for refining theoretical models of accretion and plasma behavior in the vicinity of supermassive black holes.

Galaxy: center↗
DOE OSTI2016

Influence of field emission on the propagation of cylindrical fast ionization wave in atmospheric-pressure nitrogen

The influence of field emission of electrons from surfaces on the fast ionization wave (FIW) propagation in high-voltage nanosecond pulse discharge in the atmospheric-pressure nitrogen is studied by a one-dimensional Particle-in-Cell Monte Carlo Collisions model. A strong influence of field emission on the FIW dynamics and plasma parameters is obtained. Namely, the accounting for the field emission makes possible the bridging of the cathode–anode gap by rather dense plasma (∼10{sup 13 }cm{sup −3}) in less than 1 ns. This is explained by the generation of runaway electrons from the field emitted electrons. These electrons are able to cross the entire gap pre-ionizing it and promoting the ionization wave propagation. We have found that the propagation of runaway electrons through the gap cannot be accompanied by the streamer propagation, because the runaway electrons align the plasma density gradients. In addition, we have obtained that the field enhancement factor allows controlling the speed of ionization wave propagation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Model of vertical plasma motion during the current quench

Tokamak disruptions impair plasma position control, which allows the plasma column to move and hit the wall. These detrimental events enhance thermal and mechanical loads due to halo currents and runaway electron losses. Their fundamental understanding and prevention is one of the high-priority items for ITER. As commonly observed in experiments (see e.g., Refs. 1 and 2), the disruptive plasma tends to move vertically, and the timescale of this motion is rather resistive than Alfvenic. These observations suggest that the plasma column is nearly force-free during its vertical motion. In fact, the force-free constraint is already implemented in disruption simulators (see Ref. 3). In this work, we consider a geometrically simple system that mimics the tokamak plasma surrounded by the conducting structures. Using this model, we highlight the underlying mechanism of the vertical displacement events during the current quench phase of plasma disruption. The same geometrical simplification was recently employed in Ref. 8 but without taking advantage of the force-free approximation. We find that the force-free ansatz makes the description more transparent by immediate elimination of the fast Alfvenic time-scale. This helps to reveal possible violation of ideal MHD stability during slow resistive motion of the plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

A fast low-to-high confinement mode bifurcation dynamics in the boundary-plasma gyrokinetic code XGC1

Here, a fast edge turbulence suppression event has been simulated in the electrostatic version of the gyrokinetic particle-in-cell code XGC1 in a realistic diverted tokamak edge geometry under neutral particle recycling. The results show that the sequence of turbulent Reynolds stress followed by neoclassical ion orbit-loss driven together conspire to form the sustaining radial electric field shear and to quench turbulent transport just inside the last closed magnetic flux surface. As a result, the main suppression action is located in a thin radial layer around ψ N ≃0.96–0.98, where ψ N is the normalized poloidal flux, with the time scale ~0.1 ms.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Late-time mixing and turbulent behavior in high-energy-density shear experiments at high Atwood numbers

The LANL Shear Campaign uses millimeter-scale initially solid shock tubes on the National Ignition Facility to conduct high-energy-density hydrodynamic plasma experiments, capable of reaching energy densities exceeding 100 kJ/cm3. These shock-tube experiments have for the first time reproduced spontaneously emergent coherent structures due to shear-based fluid instabilities [i.e., Kelvin-Helmholtz (KH)], demonstrating hydrodynamic scaling over 8 orders of magnitude in time and velocity. The KH vortices, referred to as “rollers,” and the secondary instabilities, referred to as “ribs,” are used to understand the turbulent kinetic energy contained in the system. Their evolution is used to understand the transition to turbulence and that transition's dependence on initial conditions. Experimental results from these studies are well modeled by the RAGE (Radiation Adaptive Grid Eulerian) hydro-code using the Besnard-Harlow-Rauenzahn turbulent mix model. Information inferred from both the experimental data and the mix model allows us to demonstrate that the specific Turbulent Kinetic Energy (sTKE) in the layer, as calculated from the plan-view structure data, is consistent with the mixing width growth and the RAGE simulations of sTKE.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Collisional transport coefficients of dense high-temperature plasmas within the quantum Landau-Fokker-Planck framework

In this paper, we extend the long-established formulas for the transport coefficients of classical plasmas inside the dense plasma regime for temperatures and densities where the classical Landau equation breaks down but its quantum extension that includes quantum degeneracy effects is valid. To this end, the quantum Landau kinetic equation is solved by the Chapman-Enskog method. The mathematical derivation is done in full generality, i.e., for multicomponent systems and to all orders of the polynomials expansion used to approximate the distribution functions. We apply the general results to two important examples, the electron gas model and an electron-ion plasma model consisting of one type of ions of any charge. We discuss the combined effects of the Pauli exclusion principle, of the electron-electron, and of the electron-ion collisions on the transport coefficients and on the convergence of the Chapman-Enskog method. For the electron gas model, the effect of the Pauli exclusion principle on the transport coefficients rapidly becomes non-negligible outside the domain of validity of the classical Landau equation. For the electron-ion plasmas, the effect of the Pauli exclusion principle depends sensitively on the ion charge Z and varies non-monotonically with Θ. For instance, for ion charge Z = 1, the electrical conductivity is increased by up to ~30% compared to its classical value over the range of degeneracy parameters studied, the thermal conductivity is reduced by up to ~9%, and the shear viscosity coefficient is increased by up to ~13%. Finally, in the Lorentz gas ( Z → ∞ ) limit, the electrical conductivity is reduced by up to ~14% compared to its classical value over the range of degeneracy parameters studied, the thermal conductivity is reduced by up to ~39%, and the shear viscosity coefficient is not affected.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Kelvin-Helmholtz instability is the result of parity-time symmetry breaking

Parity-Time (PT)-symmetry is being actively investigated as a fundamental property of observables in quantum physics. We show that the governing equations of the classical two-fluid interaction and the incompressible fluid system are PT-symmetric, and the well-known Kelvin-Helmholtz instability is the result of spontaneous PT-symmetry breaking. It is expected that all classical conservative systems governed by Newton's law admit PT-symmetry, and the spontaneous breaking thereof is a generic mechanism for classical instabilities. Here discovering the PT-symmetry of systems in fluid dynamics and plasma physics and identifying the PT-symmetry breaking responsible for instabilities enable new techniques to classical physics and enrich the physics of PT-symmetry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Current channel evolution in ideal Z pinch for general velocity profiles

Recent diagnostic advances in gas-puff Z pinches at the Weizmann Institute for the first time allow the reconstruction of the current flow as a function of time and radius. These experiments show an unexpected radially-outward motion of the current channel, as the plasma moves radially-inward. In this paper, a mechanism that could explain this current evolution is described. Here, we examine the impact of advection on the distribution of current in a cylindrically symmetric plasma. In the case of metric compression, |v r | ∝ r, the current enclosed between each plasma fluid element and the axis is conserved, and so the current profile maintains its shape. We show that for more general velocity profiles, this simple behavior quickly breaks down, allowing for non-conservation of current in a compressing conductor, rapid redistribution of the current density, and even for the formation of reverse currents. In particular, a specific inward radial velocity profile is shown to result in radially-outward motion of the current channel, recovering the surprising current evolution discovered at the Weizmann Institute.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Observation of weakly damped modes using high resolution measurement of turbulence in a dipole confined plasma

When plasma is confined by a levitated superconducting current ring, the turbulent radial particle flux changes the direction depending upon the gradient of the particle number within tubes of equal magnetic flux [Garnier et al., Phys. Plasmas 24, 012506 (2017)]. When the density profile is highly peaked, the particle flux is outward. When the density profile is less peaked, turbulence drives an inward “particle pinch.” This change in the direction of particle flux coincides with a reversal of the mean toroidal propagation direction of turbulent fluctuations. Previously, only the mean wavenumber of the fluctuations was reported. The mean wavenumber was computed from the ensemble correlation between two probes. Here, the full frequency-wavenumber power spectrum of the saturated turbulence is calculated by applying Capon’s “maximum likelihood method” to floating potential measurements from an array of probes. The power spectrum shows that the highest intensity modes are due to unstable entropy modes, and the less intense modes coincide with weakly damped convective cells, which rotate toroidally in the direction opposite to the unstable modes. Because the weakly damped modes make up a significant fraction of the fluctuating power, these measurements imply that saturated turbulence in a dipole-confined plasma involves the excitation of weakly damped rotating convective cells.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Computational analysis of ion orbital loss in diverted positive- and negative-triangularity tokamaks

The mechanism of ion orbital loss in diverted negative-triangularity (NT) tokamaks is significantly different from that of positive-triangularity (PT) tokamaks. The orbital loss from the edge of NT tokamaks can be quite large, because the separatrix X-point is located on the outboard side, meaning that trapped ions crossing the separatrix are lost to the inboard divertor plate after a mirror reflection on the side of the strong magnetic field. NT tokamaks have difficulties in building up steep pressure gradients in the edge regions. However, the majority of trapped ions crossing the separatrix in PT tokamaks survive (through complete banana motions) without reaching the X-point region. For NT tokamaks, topological changes in the trapped ion trajectories can be achieved in the presence of the radial electric field, thus reducing the orbital loss, whereas for PT tokamaks, the orbital loss is relatively insensitive to the radial electric field.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗