Search NASA⌕ Search

SEARCH · Search NASA

Results for “kinetic 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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Two-fluid and kinetic transport physics of Kelvin–Helmholtz instabilities in nonuniform low-beta plasmas

Hall-magnetohydrodynamic (Hall-MHD) theory, two-fluid simulations, and kinetic simulations are used in this work to investigate the cross-field transport properties of Kelvin–Helmholtz instabilities in nonuniform low-beta collisionless plasmas. Hall-MHD analysis shows how the linear properties of the instability are modified by density gradients and magnetization. High-order accurate two-fluid and kinetic simulations, with complete dynamics of finite-mass electrons and ions, are applied to a suite of parameter cases to systematically assess the effects of diamagnetic drift, magnetization, charge separation, and finite Larmor motion. Initialization of exact two-species kinetic equilibria facilitates the study of isolated physical effects and enables detailed cross-comparisons between two-fluid and kinetic simulations, including for cases where ion gyroradii are comparable to gradient scale lengths. For nonuniform plasmas with significant space charge, the results of two-fluid and kinetic simulations are found to disagree with Hall-MHD predictions. Kelvin–Helmholtz instability growth rates, per unit shear, are shown to be smaller when ion diamagnetic drift and $E \times B$ drift are parallel and larger when the two drifts are antiparallel. The effect is attributed to polarization drift in the shear layer, which leads to redistribution of charge, alters the electric field that drives plasma advection, and consequently modifies growth rates. Instability-induced mass transport for different parameters is characterized in terms of the flux across the shear layer and a simplified diffusion model. Distribution functions from kinetic simulations are shown to deviate substantially from Maxwellian reconstructions, indicating the importance of kinetic physics during the nonlinear phase of the instability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1994

Loop heating by D.C. electric current and electromagnetic wave emissions simulated by 3-D EM particle zone

We have shown that a current-carrying plasma loop can be heated by magnetic pinch driven by the pressure imbalance between inside and outside the loop, using a 3-dimensional electromagnetic (EM) particle code. Both electrons and ions in the loop can be heated in the direction perpendicular to the ambient magnetic field, therefore the perpendicular temperature can be increased about 10 times compared with the parallel temperature. This temperature anisotropy produced by the magnetic pinch heating can induce a plasma instability, by which high-frequency electromagnetic waves can be excited. The plasma current which is enhanced by the magnetic pinch can also excite a kinetic kink instability, which can heat ions perpendicular to the magnetic field. The heating mechanism of ions as well as the electromagnetic emission could be important for an understanding of the coronal loop heating and the electromagnetic wave emissions from active coronal regions.

Sakai, J. I.↗
NASA NTRS1992

Ion heating by kinetic cross-field streaming instability due to reflected ions at a quasiperpendicular shock

The present paper shows that the reflected ion at a supercritical quasi-perpendicular shock wave can excite a purely growing mode propagating parallel to the ambient magnetic field. To discuss the ion heating by such an unstable mode, the self-consistent quasi-linear kinetic equation is solved with the assumption that the present purely growing mode is the dominant unstable mode in the system. In the quasi-linear analysis of the instability, two particular cases are considered: the case of low initial ion beta and that of a high initial ion beta.

Yoon, Peter H.↗
DOE OSTI2024

Onset of kinetic effects on Rayleigh–Taylor instability: Advective–diffusive asymmetry

In nature and engineering applications, the Rayleigh–Taylor instability (RTI) occurs over a wide range of Atwood, Reynolds, Mach, and Knudsen numbers. At low Atwood, Mach, and Knudsen numbers, the classic advective instability causes quasi-symmetric bubble and spike growth on the two sides of the interface. However, recent findings suggest that at high degrees of rarefaction, advective effects are suppressed and molecular diffusion leads to planar growth of the density fronts on either side of the interface. This study aims to investigate the flow physics of the transition from advective to diffusive behavior, focusing on the onset of the kinetic effects. Using the gas kinetic methodology, RTI is simulated over a range of Knudsen and Mach numbers in the transition regime. The simulation results reveal the various stages of transformation from advective instability to diffusive transport. For the first time, the study demonstrates the existence of a Knudsen–Mach parameter regime where the bubble side exhibits advective instability, while the other side shows a planar density front due to molecular diffusion, rather than the canonical advective spike shape. The dominance of different mechanisms on the two sides of the interface leads to the advective–diffusive asymmetry. In conclusion, the findings of this study can lead to a more comprehensive understanding of RTI over a wide range of Mach and Knudsen numbers.

42 ENGINEERING↗
DOE OSTI2022

Verification and validation of linear gyrokinetic and kinetic-MHD simulations for internal kink instability in DIII-D tokamak

Here, verification and linear validation of the internal kink instability in tokamak have been performed for both gyrokinetic (GTC) and kinetic-MHD codes (GAM-solver, M3D-C1-K, NOVA, XTOR-K). Using realistic magnetic geometry and plasma profiles from the same equilibrium reconstruction of the DIII-D shot #141216, these codes exhibit excellent agreement for the growth rate and mode structure of the internal kink mode when all kinetic effects are suppressed. The simulated radial mode structures, obtained from linear simulations, are in reasonable agreement with the normalised electron cyclotron emission measurement after adjusting, within the experimental uncertainty, the safety factor q = 1 flux-surface location in the equilibrium reconstruction. Compressible magnetic perturbations strongly destabilize the kink, while poloidal variations of the equilibrium current density reduce the growth rate of the kink. Furthermore, kinetic effects of thermal ions are found to decrease the kink growth rate in kinetic-MHD simulations, but increase the kink growth rate in gyrokinetic simulations, due to the additional drive of the ion temperature gradient and parallel electric field. Kinetic thermal electrons are found to have negligible effects on the internal kink instability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

An improved ten-moment closure for reconnection and instabilities

The integration of kinetic effects in fluid models is important for global simulations of Earth's magnetosphere. The use of the two-fluid ten-moment model, which includes the pressure tensor for both species, has had some success in simulating Ganymede and Mercury with a simple closure model. Here, we discuss a heat flux closure which accounts for some limitations of the earlier work while remaining computationally tractable. Comparisons with kinetic simulations for magnetic reconnection and lower-hybrid drift instabilities show good agreement with kinetic results and improvements on previous closure models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Parametric study of Alfvénic instabilities driven by runaway electrons during the current quench in DIII-D

Abstract To avoid or mitigate runaway electron (RE) beams in ITER, RE-driven instabilities are actively studied as a complimentary technique to massive material injection. In this work we report experimental dependencies of Alfvénic instabilities driven by REs during the current quench in DIII-D on plasma and material injection parameters. These instabilities, observed in the frequency range of 0.1–3 MHz, correlate with increased RE loss and thus may play a role in non-sustained RE beams. It was found that as the toroidal magnetic field ( B T ) decreases, the RE population becomes more energetic, the energy of instabilities increases, and no RE beam is observed when the maximum energy of REs exceeds 15 MeV (or when B T is below 1.8 T). Analysis of disruptions at plasma core temperature ( T e ) of 1 keV and 8 keV shows that the RE population is much less energetic (with the maximum energy of only about 3 MeV) when T e is high, and no instabilities are observed in this case. Besides disruptions above caused by Ar injection, cases with Ne and D 2 injections were also studied. Both Ne and D 2 injections cause no sustained RE beams, however, for different reasons. Measurements of the instability polarization indicate that it is of predominantly compressional nature at the edge, which is consistent with modeling suggesting excitation of compressional Alfvén eigenmodes. However, drive of global Alfvén eigenmodes is also possible at low frequencies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Gravitational wave modes in matter

Here, a general linear gauge-invariant equation for dispersive gravitational waves (GWs) propagating in matter is derived. This equation describes, on the same footing, both the usual tensor modes and the gravitational modes strongly coupled with matter. It is shown that the effect of matter on the former is comparable to diffraction and therefore negligible within the geometrical-optics approximation. However, this approximation is applicable to modes strongly coupled with matter due to their large refractive index. GWs in ideal gas are studied using the kinetic average-Lagrangian approach and the gravitational polarizability of matter that we have introduced earlier. In particular, we show that this formulation subsumes the kinetic Jeans instability as a collective GW mode with a peculiar polarization, which is derived from the dispersion matrix rather than assumed a priori. This forms a foundation for systematically extending GW theory to GW interactions with plasmas, where symmetry considerations alone are insufficient to predict the wave polarization.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Spherical tokamak physics research in preparation for the operation of NSTX-U

The National Spherical Torus Experiment Upgrade (NSTX-U) is preparing to resume operation, representing a crucial step toward realizing compact, cost-effective fusion pilot plants. In advance of this, extensive modeling and data analysis have been conducted to advance the physics basis for low-aspect-ratio, high-performance plasma regimes, focusing on three core objectives: confinement and stability, power and particle handling, and steady-state operation. Significant progress has been made in understanding the electron temperature flattening in high-β plasmas, which is shown to be driven by a complex interplay of magnetohydrodynamic instabilities (e.g. non-resonant infernal modes), fast-ion-driven Alfvén eigenmodes, and electron and ion-scale micro-instabilities, particularly Kinetic Ballooning Modes (KBMs), whose destabilization is strongly dependent on parallel magnetic field fluctuations (δB ∥ ). Furthermore, a new gyrokinetic critical pedestal model was developed, accurately predicting pedestal structure by identifying KBMs as the primary stability limit, offering a critical constraint for future high-confinement scenarios. To address the challenge of high heat flux, novel liquid lithium plasma-facing components were modeled. The analysis confirmed that lithium vapor shielding is a self-regulating mechanism for heat mitigation, while also emphasizing that strong main ion parallel flow is essential to minimize core lithium contamination. Finally, progress toward steady-state operation was anchored by developing the required physics basis and control tools. This includes predictive modeling for reversed magnetic shear sustainment, demonstrating that magnetic island-induced bootstrap current reduction is negligible in STs, and advancing real-time control and disruption avoidance capabilities. The development of high-speed surrogate models (e.g. MMMNet) provides computationally efficient tools vital for non-inductive scenario optimization and integrated, low-disruptivity operations planned for NSTX-U.

NSTX-U↗
DOE OSTI2024

High-temperature carbon dioxide capture in a porous material with terminal zinc hydride sites

Carbon capture can mitigate point-source carbon dioxide (CO 2 ) emissions, but hurdles remain that impede the widespread adoption of amine-based technologies. Capturing CO 2 at temperatures closer to those of many industrial exhaust streams (>200°C) is of interest, although metal oxide absorbents that operate at these temperatures typically exhibit sluggish CO 2 absorption kinetics and instability to cycling. Here, in this study, we report a porous metal–organic framework featuring terminal zinc hydride sites that reversibly bind CO 2 at temperatures above 200°C—conditions that are unprecedented for intrinsically porous materials. Gas adsorption, structural, spectroscopic, and computational analyses elucidate the rapid, reversible nature of this transformation. Extended cycling and breakthrough analyses reveal that the material is capable of deep carbon capture at low CO 2 concentrations and high temperatures relevant to postcombustion capture.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1987

Electromagnetic beam modes driven by anisotropic electron streams

A theory is presented which provides an explanation for recent laboratory observations of a 'new electromagnetic' mode in a non-Maxwellian high-beta plasma' and suggests a likely relevance of these modes to space plasmas. Unstable parallel-propagating waves, which can be beam modes or whistlers, depending mainly on the degree of velocity anisotropy in a weak electron beam, are found. In the beam-mode limit relevant to the experiment, this is a kinetic Weibel instability, enhanced by a parallel magnetic field. The frequency of the mode can be above or below the electron cyclotron frequency.

Goldman, Martin V.↗
NASA NTRS2022

Revolutionizing Our Understanding of Particle Energization in Space Plasmas Using on-Board Wave-Particle Correlator Instrumentation

A leap forward in our understanding of particle energization in plasmas throughout the heliosphere is essential to answer longstanding questions in heliophysics, including the heating of the solar corona, acceleration of the solar wind, and energization of particles that lead to observable phenomena, such as the Earth’s aurora. The low densities and high temperatures of typical heliospheric environments lead to weakly collisional plasma conditions. Under these conditions, the energization of particles occurs primarily through collisionless interactions between the electromagnetic fields and the individual plasma particles with energies characteristic of a particular interaction. To understand how the plasma heating and particle acceleration impacts the macroscopic evolution of the heliosphere, impacting phenomena such as extreme space weather, it is critical to understand these collisionless wave-particle interactions on the characteristic ion and electron kinetic timescales. Such understanding requires high-cadence measurements of both the electromagnetic fields and the three-dimensional particle velocity distributions. Although existing instrument technology enables these measurements, a major challenge to maximize the scientific return from these measurements is the limited amount of data that can be transmitted to the ground due to telemetry constraints. A valuable, but underutilized, approach to overcome this limitation is to compute on-board correlations of the maximum-cadence field and particle measurements to improve the sampling time by several orders of magnitude. Here we review the fundamentals of the innovative field-particle correlation technique, present a formulation of the technique that can be implemented as an on-board wave-particle correlator, and estimate results that can be achieved with existing instrumental capabilities for particle velocity distribution measurements.

plasma heating↗
DOE OSTI2025

Measurement of post-disruption runaway electron kinetic energy and pitch angle during final loss instability in DIII-D

Post-disruption runaway electron (RE) kinetic energy K and pitch angle sin$\vartheta$ are critical parameters for determining resulting first wall material damage during wall strikes, but are very challenging to measure experimentally. During the final loss instability, confined RE K and sin$\vartheta$ are reconstructed during center-post wall strikes for both high impurity (high-Z) and low impurity (low-Z) plasmas by combining soft x-ray, hard x-ray, synchrotron emission, and total radiated power measurements. Deconfined (wall impacting) RE sin$\vartheta$ is then reconstructed for these shots by using time-decay analysis of infra-red imaging. Additionally, deconfined RE K and sin$\vartheta$ are reconstructed for a low-Z downward loss shot by analyzing resulting damage to a sacrificial graphite dome limiter. The damage analysis uses multi-step modeling simulating plasma instability, RE loss orbits, energy deposition, and finally material expansion (MARS-F, KORC, GEANT-4, and finally COMSOL). Overall, mean kinetic energies are found to be in the range $\langle$$K$$\rangle$ $≈ 3-4$ MeV for confined REs. KORC simulations indicate that the final loss instability process does not change individual RE kinetic energy K. In conclusion, confined RE pitch angles are found to be fairly low initially pre-instability, $\langle$sin$\vartheta$$\rangle$ $≈ 0.1 – 0.2$, but appear to increase roughly 2x, to $\langle$sin$\vartheta$$\rangle$ $≈ 0.3 – 0.4$ for both confined and deconfined REs during instability onset in the low-Z case; this increase is not observed in the high-Z case.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Nonlinear kinetic simulation study of the ion–ion streaming instability in single- and multi-ion species plasmas

The nonlinear evolution of the ion–ion streaming instability (IISI) is studied using numerical techniques novel to this problem that afford direct insight into the evolution of the particle distributions of each species. Here, during the linear phase of the instability, we demonstrate quantitative agreement with linear kinetic theory. Subsequently, the electrostatic field generated by the IISI causes ring-like velocity distributions of ions to form that are both heated and slowed to varying degrees relative to their initial flows. Due to variation in the trapping conditions for ion species of differing charge-to-mass ratio, when flows of multiple species interact, the nonlinear evolution of each species can be starkly different: we show a case where a lighter ion species is completely stopped by a heavier ion species via the IISI alone (i.e., without collisions) and, for the first time, demonstrate how the IISI can introduce a relative flow between ion species that initially have the same flow velocities, thereby separating them.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1992

Nonlocal analysis of finite-beam-driven instabilities

The fully kinetic integral eigenmode equation in wave-number space is used to describe the nonlocal behavior of electrostatic waves in an electron-beam plasma, which are studied in the low-temperature-beam regime and the warm-beam regime. The case of strongly magnetized electrons and unmagnetized ions, which corresponds to the waves in a frequency range from the lower-hybrid to the electron plasma frequency, is examined. Three wave modes are found. The first group consists of modes that have dispersive properties similar to the uniform, infinite beam-plasma system. Depending on the beam width, the growth rates are strongly reduced. The second group, surface modes, are localized at the periphery of the beam region and are less unstable than the unstable modes of the first group. The third group represents natural oscillations of the background plasma. These modes are virtually unaffected by the beam.

Serizawa, Y.↗
NASA NTRS1993

Kinetic simulations of the Kelvin-Helmholtz instability at the magnetopause

Two-dimensional hybrid simulations with particle ions and fluid electrons are used to calculate the kinetic evolution of the Kelvin-Helmholtz instability for a magnetopauselike configuration. The unidirectional magnetic field is essentially transverse to the plasma flow velocity, which is the most unstable case according to linear theory and models the flow dynamics in the subsolar region of the magnetopause for northward interplanetary magnetic field. We recover effects analogous to those found in MHD simulations, including a mode cascade to longer wavelengths. The boundary layer consists of coherent structure and is not well described by a diffusive process. Isolated structures on the order of the ion gyroradius are formed which can cross the boundary in either direction. We describe how the time evolution of these structures represents transport across boundary layers, and we consider the possible connection of these entities to flux transfer events and other structure seen in the low-latitude boundary layer at the Earth's magnetopause as well as to flux ropes commonly observed near the ionopause of Venus. We also discuss the relation of the hybrid calculations to previous MHD simulations and to observations.

Thomas, V. A.↗
DOE OSTI2021

Alfvén eigenmode classification based on ECE diagnostics at DIII-D using deep recurrent neural networks

Modern tokamaks have achieved significant fusion production, but further progress towards steady-state operation has been stymied by a host of kinetic and MHD instabilities. Control and identification of these instabilities is often complicated, warranting the application of data-driven methods to complement and improve physical understanding. In particular, Alfvén eigenmodes are a class of ubiquitous mixed kinetic and MHD instabilities that are important to identify and control because they can lead to loss of confinement and potential damage to the walls of a plasma device. In the present work, we use reservoir computing networks to classify Alfvén eigenmodes in a large labeled database of DIII-D discharges, covering a broad range of operational parameter space. However, despite the large parameter space, we show excellent classification and prediction performance, with an average hit rate of 91% and false alarm ratio of 7%, indicating promise for future implementation with additional diagnostic data and consolidation into a real-time control strategy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗