Search NASA⌕ Search

SEARCH · Search NASA

Results for “drift-kinetic theory”

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.

7 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Fast Particle-Wave Interactions and Alfvén Eigenmodes in JET Tokamak Plasmas

This document serves as the closeout report of DOE Grant Award No. DE-FG02-99ER54563 with project period 1 April 2015 through 31 March 2020. The project comprised the international collaboration between MIT and EU scientists on the JET facility to improve our understanding of the physics of energetic particle-wave interactions by measuring the damping rates of stable Alfvén Eigenmodes (AEs) and unstable energetic particle driven modes. More specifically, this project involved the continued participation of MIT, the Culham Centre for Fusion Energy, the Swiss Plasma Center, and theorists from various European laboratories and from UC Irvine. Past contributions from the University of São Paulo are gratefully acknowledged. The Alfvén Eigenmode Active Diagnostic (AEAD) on JET discharges was successfully upgraded in this period and has extracted physically useful information, in primarily deuterium plasma discharges. It is expected that such experiments would be continued during the DT campaign in CY 2021 to assess the damping rates of similar modes in the presence of alpha particles, a product of fusing burning plasma. Only JET would carry out such experiments in the near term in the world. It is important to note that this collaboration was continued under DOE Grant Award No. DE-SC0014264 from 1 April 2020 through the present. In this grant period, the AEAD was upgraded with a set of individual amplifiers for each of a set of six antennas in two toroidally opposite locations. These new amplifiers allowed targeted selection of the antennas’ toroidal spectrum for AEs of interest with toroidal mode number | n | ≤ 20. The resonant detection and measurement of the damping rates of AEs was obtained from magnetic probes that could be compared with theory and simulations. This is a key topic of investigation for ITER and all other next-step fusion experiments, where such modes will interact with energetic particles produced by fusion reactions (alpha particles), Neutral Beam Injection (NBI) and Ion Cyclotron Resonance Heating (ICRH). The majority of the observations during this period were of Toroidal Alfvén Eigenmodes (TAEs) as these are most commonly observed on JET; a number of measurements were made during dedicated TAE experiments. However, a new set of lower frequency band filters were procured with the goal of studying Geodesic Acoustic Modes (GAMs), Beta (Acoustic) AEs (BAE/BAAEs), and Reverse Shear AEs (RSAEs), also predicted by theory. Commissioning and optimizations were successfully completed in this period, and the diagnostic has been in operation during the more recent JET campaigns. Initial results have been obtained in dedicated TAE experiments and successfully compared to drift-kinetic theory. A wide range of theoretical studies have been undertaken in support of the upcoming JET campaigns experiments. To supplement the ongoing use of ideal MHD codes, such as MISHKA during the studies of TAEs, gyrokinetic simulations of low frequency AEs were performed in collaboration with UC Irvine. The Gyrokinetic Toroidal Code (GTC) was used to determine the structure, frequencies and stability of AEs in JET plasmas. Thus, a solid scientific foundation has been laid for future DT campaigns in JET in the 2021 operation period when the damping rates of relevant modes in the presence of alpha particles could be assessed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Scattering theory in noncanonical phase space: A Drift-Kinetic collision operator for weakly collisional plasmas

After developing a scattering theory for grazing collisions in general noncanonical phase spaces, we introduce a guiding center collision operator in five-dimensional phase space designed for plasma regimes characterized by long wavelengths (relative to the Larmor radius), low frequencies (relative to the cyclotron frequency), and weak collisionality (where repeated Coulomb collisions induce cumulatively small changes in particle magnetic moment). The collision operator is fully determined by the noncanonical Hamiltonian structure of guiding center dynamics and exhibits a metriplectic structure, ensuring the conservation of particle number, momentum, energy, and interior Casimir invariants. It also satisfies an H-theorem, allowing for deviations from an equilibrium Maxwellian distribution due to the nontrivial kernel of the noncanonical guiding center Poisson tensor, spanned by the magnetic moment. We propose that this collision operator and its underlying mathematical structure may offer valuable insight into the study of turbulence, transport, and self-organizing phenomena in both laboratory and astrophysical plasmas.

Hamiltonian mechanics↗
DOE OSTI2021

Accurate numerical, integral methods for computing drift-kinetic Trubnikov-Rosenbluth potentials

A novel numerical method is employed to compute the integral form of the axi-symmetric Trubnikov-Rosenbluth potentials. Two methods for quadrature in pitch-angle are described and their convergence properties are studied. Careful attention is given to quadrature over a singular Green's function. Here it is shown that an infinite series representation of the Green's function can be used more efficiently than its closed form involving complete elliptic integrals. Then a collocation method in speed, with its associated quadrature scheme, is laid out and its convergence properties are studied. Using the proposed scheme, accurate low-order moments of the field collision operator are obtained using relatively few velocity space degrees of freedom. The scheme is showcased by solving for the equilibrium, axi-symmetric bootstrap current in tokamaks. A C 0 Gauss-Lobatto-Legendre finite element pitch-angle basis with vertex nodes at the trapped/passing boundary is shown, in the context of the integral methods used, to be much more efficient than the more common Legendre polynomial expansion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Dispersion relation for gauge-free electromagnetic drift kinetics

Recently, a new approach to gyrokinetics, invariant under electromagnetic gauge transformations, was developed. The gyrocenter equations of motion are now expressed in terms of the perturbed fields instead of the potentials, in a form suitable for numerical simulations and analytic studies. In this paper, we verify that the long-wavelength limit, i.e., the drift-kinetic limit of the new gyrokinetic theory, is in line with existing work, providing a solid foundation for simulations. Here, we compute the dispersion relation of the new drift-kinetic theory in slab geometry and find agreement with a long-wavelength limit of the full Vlasov–Maxwell model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

MAS: A versatile Landau-fluid eigenvalue code for plasma stability analysis in general geometry

We have developed a new global eigenvalue code, Multiscale Analysis for plasma Stabilities (MAS), for studying plasma problems with wave toroidal mode number (n) and frequency (ω) in a broad range of interest in general tokamak geometry, based on a five-field Landau-fluid description of thermal plasmas. Beyond keeping the necessary plasma fluid response, we further retain the important kinetic effects including diamagnetic drift, ion finite Larmor radius, finite parallel electric field (E||), ion and electron Landau resonances in a self-consistent and non-perturbative manner without sacrificing the attractive efficiency in computation. The physical capabilities of the code are evaluated and examined in the aspects of both theory and simulation. In theory, the comprehensive Landau-fluid model implemented in MAS can be reduced to the well-known ideal MHD model, electrostatic ion-fluid model, and drift-kinetic model in various limits, which clearly delineates the physics validity regime. In simulation, MAS has been well benchmarked with theory and other gyrokinetic and kinetic-MHD hybrid codes in a manner of adopting the unified physical and numerical framework, which covers the kinetic Alfv\'en wave (KAW), ion sound wave (ISW), low-n kink, high-n ion temperature gradient mode (ITG) and kinetic ballooning mode (KBM). Moreover, MAS is successfully applied to model the Alfv\'en eigenmode (AE) activities in DIII-D discharge #159243, which faithfully captures the frequency sweeping of reversed shear Alfv\'en eigenmode (RSAE), the tunneling damping of toroidal Alfv\'en eigenmode (TAE), as well as the polarization characteristics of kinetic beta-induced Alfv\'en eigenmode (KBAE) and beta-induced Alfv\'en-acoustic eigenmode (BAAE) being consistent with former gyrokinetic theory and simulation. With respect to the key progress contributed to the community, MAS has the advantage of combining rich physics ingredients, realistic global geometry and high computation efficiency together for plasma stability analysis in linear regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2017

Kinetic Theory and Fast Wind Observations of the Electron Strahl

We develop a model for the strahl population in the solar wind - a narrow, low-density and high-energy electron beam centred on the magnetic field direction. Our model is based on the solution of the electron drift-kinetic equation at heliospheric distances where the plasma density, temperature and the magnetic field strength decline as power laws of the distance along a magnetic flux tube. Our solution for the strahl depends on a number of parameters that, in the absence of the analytic solution for the full electron velocity distribution function (eVDF), cannot be derived from the theory. We however demonstrate that these parameters can be efficiently found from matching our solution with observations of the eVDF made by the Wind satellite's SWE strahl detector. The model is successful at predicting the angular width (FWHM) of the strahl for the Wind data at 1 au, in particular by predicting how this width scales with particle energy and background density. We find that the strahl distribution is largely determined by the local temperature Knudsen number γ ∼ |T dT/dx|/n, which parametrizes solar wind collisionality. We compute averaged strahl distributions for typical Knudsen numbers observed in the solar wind, and fit our model to these data. The model can be matched quite closely to the eVDFs at 1 au; however, it then overestimates the strahl amplitude at larger heliocentric distances. This indicates that our model may be improved through the inclusion of additional physics, possibly through the introduction of 'anomalous diffusion' of the strahl electrons.

Horaites, Konstantinos↗
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↗