Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Plasma Instabilities”

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

First observations of edge instabilities in strongly shaped negative triangularity plasmas on DIII-D

Throughout the 2023 DIII-D negative triangularity (NT) campaign, unexpected magnetic fluctuations indicative of an edge localized instability were observed in the majority of strongly shaped (–0.55 ≤ δ ≤ –0.45) discharges. Moreover, these same fluctuations were found to be absent in weaker NT plasma shapes. Based on a database study, three separate mode behaviors are observed based primarily on the edge safety factor and normalized plasma beta β N . The main regime is the bursty state, in which short (~5 ms) magnetic bursts occur every 5–10 ms. Each burst drives a small reduction in the edge electron temperature (T e ) as measured by electron cyclotron emissions (ECE), which then rapidly recovers between the bursts. The quiescent regime appears to have governing mechanics similar to the bursty regime, but with multiple unique mode signatures observed in the magnetics. These additional modes greatly reduce the bursty nature of the plasma. Finally, a saturated state is observed in a number of discharges, particularly at early time slices. The saturated regime is often short-lived (maximum 500 ms), but can have a large impact on the dynamics of the plasma edge region, causing significant drops in electron temperature as observed in the ECE diagnostic. 2D ECE imaging (ECE-I) is utilized, showing large kink-like perturbations localized near the last closed flux surface that extend into the SOL indicative of an external kink/peeling behavior. Additionally, evidence of phase mixing and observed tearing-parity perturbations localized near the edge rational surfaces is observed in the quiescent regime plasmas. As a result, the various modes do not appear to have a large impact on global confinement but may play a role in regulating the transport of the NT edge and the observed edge gradients.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Energetic particle physics: Chapter 7 of the special issue: on the path to tokamak burning plasma operation

We review the physics of energetic particles (EPs) in magnetically confined burning fusion plasmas with focus on advances since the last update of the ITER Physics Basis (Fasoli et al 2007 Nucl. Fusion 47 S264). Topics include basic EP physics, EP generation, diagnostics of EPs and instabilities, the interaction of EPs and thermal plasma instabilities, EP-driven instabilities, energetic particle modes (EPMs), and turbulence, linear and nonlinear stability and simulation of EP-driven instabilities and EPMs, 3D effects, scenario optimization strategies based on EP phase-space control, EPs in reduced field scenarios in ITER before DT, and the physics of runaway electrons. We describe the simulation and modeling of EPs in fusion plasmas, including instability drive and damping as well as EP transport, with a range of approaches from first-principles to reduced models, including gyrokinetic simulations, kinetic-MHD models, gyrofluid models, reduced models, and semi-analytical approaches.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Theory of gradient drift instabilities in low-temperature, partially magnetised plasmas

A fluid dispersion theory in partially magnetised plasmas is analysed to examine the conditions under which large-wavelength modes develop in Penning-type configurations, that is, where an electric field is imposed perpendicular to a homogeneous magnetic field. The fluid dispersion relation assuming a slab geometry shows that two types of low-frequency, gradient drift instabilities occur in the direction of the E×B and diamagnetic drifts. One type of instability, observed when the equilibrium electric field and plasma density gradient are in the same direction, is similar to the classic modified Simon–Hoh instability. A second instability is found for conditions in which (i) the diamagnetic drift is in the direction opposite to the E×B drift and (ii) the magnitude of the diamagnetic drift is sufficiently larger than the electron thermal speed. The present fluid dispersion theory suggests that the rotating spokes driven by such fluid instabilities propagate in the same direction as the diamagnetic drift, which can be in the same direction as or opposite to the E×B drift, depending on the plasma conditions. This finding may account for the observation, in some plasma devices, of the rotation of large-scale structures in both the E×B and -E×B directions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

On parallel laser beam merger in plasmas

Self-focusing instability is a well-known phenomenon of nonlinear optics, which is of great importance in the field of laser–plasma interactions. Self-focusing instability leads to beam focusing and, consequently, breakup into multiple laser filaments. The majority of applications tend to avoid a laser filamentation regime due to its detrimental role on laser spot profile and peak intensity. In our work, using nonlinear Schrödinger equation solver and particle-in-cell simulations, we address the problem of interaction of multiple parallel beams in plasmas. We consider both non-relativistic and moderately relativistic regimes and demonstrate how the physics of parallel beam interaction transitions from the familiar self- and mutual-focusing instabilities in the non-relativistic regime to a moderately relativistic regime, where an analytical description of filament interaction is not available.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

A two-dimensional numerical study of ion-acoustic turbulence

We investigate the linear and nonlinear evolution of the current-driven ion-acoustic instability in a collisionless plasma via two-dimensional (2-D) Vlasov–Poisson numerical simulations. We initialise the system in a stable state and gradually drive it towards instability with an imposed, weak external electric field, thus avoiding physically unrealisable super-critical initial conditions. A comprehensive analysis of the nonlinear evolution of ion-acoustic turbulence (IAT) is presented, including the detailed characteristics of the evolution of the particles’ distribution functions, (2-D) wave spectrum and the resulting anomalous resistivity. Our findings reveal the dominance of 2-D quasi-linear effects around saturation, with nonlinear effects, such as particle trapping and nonlinear frequency shifts, becoming pronounced during the later stages of the system's nonlinear evolution. Remarkably, the Kadomtsev–Petviashvili (KP) spectrum is observed immediately after the saturation of the instability. Another crucial and noteworthy result is that no steady saturated nonlinear state is ever reached: strong ion heating suppresses the instability, which implies that the anomalous resistivity associated with IAT is transient and short-lived, challenging earlier theoretical results. Towards the conclusion of the simulation, electron-acoustic waves are triggered by the formation of a double layer and strong modifications to the particle distribution induced by IAT.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Coherent High-Frequency Axial Oscillations in a Partially Magnetized Direct Current Magnetron Discharge

High-frequency oscillations are observed in a neon plasma of a direct current magnetron discharge. At low discharge currents, we see highly coherent 60 MHz fluctuations. Above a distinct current threshold, secondary 5–10 MHz fluctuations emerge in addition to turbulent fluctuations in the 60–100 MHz range. The oscillations in the total discharge current suggest axial wave propagation. A lower-hybrid wave theory is invoked to model the high-frequency oscillations. Here, we attribute the low-frequency modes to a turbulence-driven inverse cascade process, as suggested by recent simulations.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2021

Dissipation in the magnetic turbulence of reversed field pinch plasmas

Reversed field pinch plasmas are subject to tearing instability that creates a broad spectrum of magnetic fluctuations. The dominant fluctuations have poloidal and toroidal mode numbers (m,n) = (1,6 - 10) and can grow to 2 - 3% of the mean magnetic field. Through non-linear coupling this growth culminates in a strong reconnection event and broadening of the magnetic spectrum extending to the ion gyroradius scale. Multiple developments occur during the reconnection stage: ions and electrons are energized, magnetic fluctuation amplitudes increase, plasma flow is halted, and the toroidal magnetic flux increases in a sawtooth-like fashion as the RFP dynamo becomes stronger. Magnetic fluctuations are measured in the plasma edge at multiple radial locations from r/a = 0.75 to 0.96 to assess and characterize the magnetic turbulence. The measured spectrum perpendicular to the mean field, S(k ⊥ ), can be fit to a model spectrum consisting of power law and exponential component with one free parameter that characterizes dissipation. The measured dissipation is much larger than estimated from classical viscous or resistive dissipation, but it is consistent with a flow damping measurement of anomalous viscosity. The measurements show an evolution of the spectrum during which fluctuation power builds up in the smallest wavenumbers and cascades to the larger wavenumber due to the non-linear coupling between the linear (m,n) = (1, > 6) and the non-linear (m,n) = (0, 1) tearing mo

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Electron Weibel instability induced magnetic fields in optical-field ionized plasmas

Generation and amplification of magnetic fields in plasmas is a long-standing topic that is of great interest to both plasma and space physics. The electron Weibel instability is a well-known mechanism responsible for self-generating magnetic fields in plasmas with temperature anisotropy and has been extensively investigated in both theory and simulations, yet experimental verification of this instability has been challenging. Recently, we demonstrated a new experimental platform that enables controlled initialization of highly nonthermal and/or anisotropic plasma electron velocity distributions via optical-field ionization. Using an external electron probe bunch from a linear accelerator, the onset, saturation, and decay of the self-generated magnetic fields due to electron Weibel instability were measured for the first time to our knowledge. In this paper, we will first present experimental results on time-resolved measurements of the Weibel magnetic fields in non-relativistic plasmas produced by Ti:Sapphire laser pulses (0.8 μm) and then discuss the feasibility of extending the study to a quasi-relativistic regime by using intense CO 2 (e.g., 9.2 μm) lasers to produce much hotter plasmas.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Collaborative Research: Parametric Instabilities of Alfven Waves in Low‐beta Plasmas

Alfven waves are of fundamental importance in magnetized plasmas. This project aims to advance our understanding of the basic physics of the parametric instabilities of Alfven waves using advanced computer simulations. The results will have implications for important processes in the heliosphere involving Alfven waves, such as coronal heating and solar wind heating.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Electron Modulational Instability in the Strong Turbulent Regime for an Electron Beam Propagating in a Background Plasma

We study collective processes for an electron beam propagating through a background plasma using simulations and analytical theory. A new regime where the instability of a Langmuir wave packet can grow locally much faster than ion frequency is clearly identified. The key feature of this new regime is an electron modulational instability that rapidly creates a local Langmuir wave packet, which in its turn produces local charge separation and strong ion density perturbations because of the action of the ponderomotive force, such that the beam-plasma wave interaction stops being resonant. Three evolution stages of the process and observed periodic burst features are discussed. In conclusion, different physical regimes in the plasma and beam parameter space are demonstrated for the first time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

On the study of hydrodynamic instabilities in the presence of background magnetic fields in high-energy-density plasmas

Blast-wave-driven hydrodynamic instabilities are studied in the presence of a background B-field through experiments and simulations in the high-energy-density (HED) physics regime. In experiments conducted at the Laboratoire pour l’utilisation des lasers intenses (LULI), a laser-driven shock-tube platform was used to generate a hydrodynamically unstable interface with a prescribed sinusoidal surface perturbation, and short-pulse x-ray radiography was used to characterize the instability growth with and without a 10-T B-field. The LULI experiments were modeled in FLASH using resistive and ideal magnetohydrodynamics (MHD), and comparing the experiments and simulations suggests that the Spitzer model implemented in FLASH is necessary and sufficient for modeling these planar systems. These results suggest insufficient amplification of the seed B-field, due to resistive diffusion, to alter the hydrodynamic behavior. Although the ideal-MHD simulations did not represent the experiments accurately, they suggest that similar HED systems with dynamic plasma-β (=2μ0ρv2/B2) values of less than ∼100 can reduce the growth of blast-wave-driven Rayleigh–Taylor instabilities. These findings validate the resistive-MHD FLASH modeling that is being used to design future experiments for studying B-field effects in HED plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Verification of electromagnetic simulation capabilities in global gyrokinetic particle-in-cell code GTS

Recently, the numerical scheme presented by Mishchenko et al. enabled explicit gyrokinetic simulations of low-frequency electromagnetic instabilities in tokamaks at experimentally relevant values of plasma β⁠. This scheme resolved the long-standing cancellation problem that previously hindered gyrokinetic particle-in-cell code simulations of magnetohydrodynamic phenomena with inherently small parallel electric fields. Moreover, the scheme did not employ approximations that eliminate critical tearing-type instabilities. Here, we report on the implementation of this numerical scheme in the global gyrokinetic particle-in-cell code GTS. This implementation allows for a more complete and accurate picture of interaction between small scale turbulence and MHD modes in tokamaks. Additionally, we present a comprehensive set of verification simulations of numerous electromagnetic instabilities relevant to present-day tokamaks. These simulations encompass the kinetic ballooning mode, the internal kink mode, the tearing mode, the micro-tearing mode, and the toroidal Alfven eigenmode destabilized by energetic ions, which are all instrumental in understanding tokamak physics. We will also showcase the preliminary nonlinear simulations of kinetic ballooning instabilities and (2,1) island formation due to tearing mode instability. These simulations validate the accuracy of the scheme implementation and pave the way for studying how these instabilities affect plasma confinement and performance.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Transition from ITG to MTM linear instabilities near pedestals of high density plasmas

Here, investigation of linear gyrokinetic ion-scale modes ( k θ ρ s = 0.3) finds that a transition from ion temperature gradient to microtearing mode (MTM) dominance occurs as the density is increased near the pedestal region of a parameterized DIII-D sized tokamak. H-modes profile densities, temperatures, and equilibria are parameterized utilizing the OMFIT PRO_create module. With these profiles, linear gyrokinetic ion-scale instabilities are predicted with CGYRO. This transition (nMTM) has a weak dependence on radial location in the region near the top of the pedestal ( ρ = 0.7 – 0.9), which allows simulating single radii to examine the approximate scaling of nMTM with global parameters. The critical nMTM is found to scale with plasma current. Additionally, increasing the minor radius by decreasing the aspect ratio and increasing the major radius are found to reduce nMTM. However, any relationship between nMTM and density limit physics remains unclear as nMTM increases relative to the Greenwald density with larger minor radius and with larger magnetic field, suggesting that the transport due to MTM may be less important for a reactor. Additionally, nMTM is sensitive to the pedestal temperature, the local electron and ion gradients, the ratio of ion to electron temperature T i / T e, and the current profile. MTMs are predicted to be the dominant instability in the core at similar Greenwald fractions for DIII-D, NSTX, and NSTX-U H-mode experiments, supporting the results of the parameterized study. Additionally, MTMs continue to be the dominant linear instability in a DIII-D L-mode after an H–L transition as the plasma approaches a density limit disruption despite the large change in plasma profiles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Transition from retrograde to prograde drift instabilities in a magnetron microdischarge

Rotating plasma structures, or “spokes,” in magnetized discharges characterized by perpendicular electric and magnetic fields have been seen in a growing number of studies and are believed to be the result of gradient-driven drift instabilities. Under certain conditions, we have shown the spoke’s rotation to be opposite to the E B direction, i.e., retrograde in its expected direction. Recently, we have linked such counter-intuitive rotation to a local inversion of the electric field. Here, we give further experimental evidence for this inversion and attempt to provide an explanation for a relatively distinct transition seen between retrograde and positive (prograde) drift. In the experiments, a partially magnetized plasma forms inside a magnetron device of mm radius operated with argon. Discharge current–voltage measurements are acquired for a range of argon fill pressure and inter-electrode spacing. We find two branches of operation—a low current branch of negative resistance, coinciding with the retrograde spoke rotation, and a higher current branch of positive resistance, coincident with prograde spoke rotation. We postulate that at low discharge currents, high magnetic field confinement leads to a large density gradient, causing more electron transport to the anode than that demanded by the external circuit. At higher currents, anomalous axial electron transport (across the magnetic field lines) becomes dominant, eliminating the conditions for field inversion. Here, the current thresholds for the field inversion are found to be sensitive to the argon pressure and inter-electrode spacing.

42 ENGINEERING↗
DOE OSTI2021

Nonlinear and 3D plasmas. Final report

The project objectives are to characterize the nonlinear properties of edge instabilities in fusion plasmas from the standpoint of MHD and extended MHD models of plasmas and to develop, as much as possible, simplified models that may be useful for prediction and control.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Bootstrap current modeling in M3D-C1

Bootstrap current plays a crucial role in the equilibrium of magnetically confined plasmas, particularly in quasi-symmetric stellarators and in tokamaks, where it can represent bulk of the electric current density. Accurate modeling of this current is essential for understanding the magnetohydrodynamic (MHD) equilibrium and stability of these configurations. This study expands the modeling capabilities of M3D-C1, an extended-MHD code, by implementing self-consistent physics models for bootstrap current. It employs two analytical frameworks: a generalized Sauter model (Sauter et al. 1999 Phys. Plasmas vol. 6, no. 7, pp. 2834–2839), and a revised Sauter-like model (Redl et al. 2021 Phys. Plasmas vol. 28, no. 2, pp. 022502). The isomorphism described by Landreman et al. (2022 Phys. Rev. Lett. vol. 128, pp. 035001) is employed to apply these models to quasi-symmetric stellarators. The implementation in M3D-C1 is benchmarked against neoclassical codes, including NEO, XGCa and SFINCS, showing excellent agreement. These improvements allow M3D-C1 to self-consistently calculate the neoclassical contributions to plasma current in axisymmetric and quasi-symmetric configurations, providing a more accurate representation of the plasma behavior in these configurations. A workflow for evaluating the neoclassical transport using SFINCS with arbitrary toroidal equilibria calculated using M3D-C1 is also presented. This workflow enables a quantitative evaluation of the error in the Sauter-like model in cases that deviate from axi- or quasi-symmetry (e.g. through the development of an MHD instability).

fusion plasma↗
DOE OSTI2023

The theory of kinetic effects on resistive wall mode stability in tokamaks

Tokamak fusion plasmas benefit from high pressures but are then susceptible to modes of instability. These magnetohydrodynamic (MHD) modes are macroscopic distortions of the plasma, but certain collective motions of individual particles can provide stabilizing effects opposing them. The presence of a resistive wall slows the mode growth, converting a kink to a resistive wall mode (RWM). A kinetic MHD model includes Maxwell's equations, ideal MHD constraints, and kinetic effects included through the pressure tensor, calculated with the perturbed drift-kinetic distribution function of the particles. The kinetic stabilizing effects on the RWM arise through resonances between the plasma rotation and particle drift motions: precession, bounce, and transit. A match between particle motions and the mode allows efficient transfer of energy that would otherwise drive the growth of the mode, thus damping the growth. The first approach to calculating RWM stability is to write a set of equations for the complex mode frequency in terms of known quantities and then to solve the system. The “energy principle” approach, which has the advantage of clarity in distinguishing the various stabilizing and destabilizing effects, is to change the force balance equation into an equation in terms of changes of kinetic and potential energies, and then to write a dispersion relation for the mode frequency in terms of those quantities. These methods have been used in various benchmarked codes to calculate kinetic effects on RWM stability. Importantly, the theory has illuminated the important roles of plasma rotation, energetic particles, and collisions in RWM stability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Mitigation of space charge effects using electron column at IOTA ring

We investigate a novel method to mitigate space charge effects of high intensity proton beams propagating in circular accelerators by means of trapping and controlling electrons generated from beam-induced residual gas ionization. This compensation method uses Coulomb repulsion force between a proton beam and electrons to mitigate self-space charge effects of the beam if it passes through a plasma column. The transverse electron-proton (e-p) instability in the plasma column is well controlled by the longitudinal magnetic field of a solenoid magnet and the bias voltages on electrodes. In this report, we will show simulation results how to control distributions of electrons and ions as well as that of the proton beam inside the column.

43 PARTICLE ACCELERATORS↗