Search NASA⌕ Search

Engineering topics

Breizman, B. N.

Publications and source records attributed to Breizman, B. N..

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Nonlinear saturation of whistler modes driven by runaway electrons

Runaway electrons exhibit kinetic instabilities with potentially beneficial consequences. The anomalous Doppler resonance between the electrons and whistler modes is a primary underlying mechanism. These instabilities require a first-principle nonlinear theoretical analysis, which would ultimately enable assessment of their impact on disruption mitigation and diagnostics, especially in ITER-relevant conditions. This paper presents recent progress in developing the required theoretical framework for isolated nonlinear resonances. By employing the generic bump-on-tail model, we predict the wave saturation levels in both near-threshold and strongly driven regimes. Remarkably, for the waves of interest, the parallel component of the wave vector dictates the parallel momentum of the resonant particles. This feature, together with the expected wave saturation level, provides a complete description of the resonance impact on the runaway electron distribution function. We show that the parametric dependence of the nonlinearly saturated mode provides a way to recover certain features of the runaway momentum distribution function.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Marginal stability constraint on runaway electron distribution

High-frequency kinetic instabilities of the strongly anisotropic runaway electrons (RE) can enhance the pitch-angle scattering of the runaways significantly. This wave-induced scattering can easily prevail over runaway scattering on high-Z impurities. In a steady state, collisional damping balances the kinetic drive of the unstable waves, keeping the RE distribution function at marginal stability. The marginal stability constraint limits the achievable RE densities and the shape of the RE distribution function. In this study, we consider whistler and compressional Alfvén waves as the primary source of enhanced elastic scattering of the runaways. By balancing the anomalous Doppler resonance drive with the collisional wave damping, we find the RE distribution function in the ultra-relativistic range of the phase space. We also derive an expression for the spectral energy density of the waves. We show that the power needed to compensate for the wave dissipation is negligible compared to the work of the electric field. The latter is in balance with the synchrotron losses of the runaway electrons.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Bifurcation-driven vertical plasma displacement

This paper considers vertical plasma motion resulting from plasma current decay during the disruption event. The presented filament-based model describes the motion in the ideal wall limit as an adiabatically slow evolution of the plasma equilibrium. Here, the equilibrium exhibits a pitchfork bifurcation when the decaying plasma current passes a critical value determined by the external magnetic field. This bifurcation affects the disruption-induced mechanical loads on the first wall.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Simulation of convective transport during frequency chirping of a TAE using the MEGA code

Here, we present a procedure to examine energetic particle phase-space during long range frequency chirping phenomena in tokamak plasmas. To apply the proposed method, we have performed self-consistent simulations using the MEGA code and analyzed the simulation data. We demonstrate a traveling wave in phase-space and that there exist specific slices of phase-space on which the resonant particles lie throughout the wave evolution. For non-linear evolution of an n = 6 toroidicity-induced Alfvén eigenmode (TAE), our results reveal the formation of coherent phase-space structures (holes/clumps) after coarse-graining of the distribution function. These structures cause a convective transport in phase-space which implies a radial drift of the resonant particles. We also demonstrate that the rate of frequency chirping increases with the TAE damping rate. Our observations of the TAE behavior and the corresponding phase-space dynamics are consistent with the Berk–Breizman theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Experimental observation of the geodesic acoustic frequency limit for the NBI-driven Alfvén eigenmodes in TJ-II

We study Alfvén Eigenmodes (AEs) in the TJ-II heliac in hydrogen plasmas heated by hydrogen cofield NBI. Taking advantage of the unique TJ-II flexibility in varying plasma current, we have observed strong variation of the AE frequency from f AE ~ 30 to ~ 220 kHz for selected modes. An advanced heavy ion beam probe diagnostic determines the spatial location and internal amplitudes of the modes. The modes satisfy a local AE dispersion relation including the geodesic acoustic frequency that represents the lowest frequency of the mode. Linear MHD-modeling with STELLGAP and FAR3D codes shows that the calculated temporal evolution of the mode frequency reproduces the observed maxima and minima at the same time intervals with a similar frequency range and the radial profile peaks near the outer edge of the observed one.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Long range frequency chirping of Alfvén eigenmodes

In this work, a theoretical framework has been developed for an NBI scenario to model the hard non-linear evolution of global Alfvén eigenmodes (GAEs) where the adiabatic motion of phase-space structures (holes and clumps), associated with frequency chirping, occurs in generalized phase-space of slowing down energetic particles. The radial profile of the GAE is expanded using finite elements, which allows an update of the mode structure as the mode frequency chirps. Constants of motion are introduced to track the dynamics of energetic particles during frequency chirping by implementing proper action-angle variables and canonical transformations which reduce the dynamics essentially to 1D. Consequently, we specify whether the particles are drifting inward/outward as the frequency deviates from the initial MHD eigenfrequency. Using the principle of least action, we have derived the non-linear equation describing the evolution of the radial profile by varying the total Lagrangian of the system with respect to the weights of the finite elements. For the choice of parameters in this work, it is shown that the peak of the radial profile is shifted and also broadens due to frequency chirping. The time rate of frequency change is also calculated using the energy balance and we show that the adiabatic condition remains valid once it is satisfied. This model clearly illustrates the theoretical treatment to study the long-range adiabatic frequency sweeping events observed for Alfvén gap modes in real experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Modification of the Alfvén wave spectrum by pellet injection

Alfvén eigenmodes driven by energetic particles are routinely observed in tokamak plasmas. These modes consist of poloidal harmonics of shear Alfvén waves coupled by inhomogeneity in the magnetic field. In this work, further coupling is introduced by 3D inhomogeneities in the ion density during the assimilation of injected pellets. This additional coupling modifies the Alfvén continuum and discrete eigenmode spectrum. We find that the frequencies of Alfvén eigenmodes drop dramatically when a pellet is injected in JET. From these observations, information about the changes in the ion density caused by a pellet can be inferred. To use Alfvén eigenmodes for MHD spectroscopy of pellet injected plasmas, the 3D MHD codes Stellgap and AE3D were generalised to incorporate 3D density profiles. Additionally, a model for the expansion of the ionised pellet plasmoid along a magnetic field line was derived from the fluid equations. Thereby, the time evolution of the Alfvén eigenfrequency is reproduced. By comparing the numerical frequency drop of a toroidal Alfvén eigenmode (TAE) to experimental observations, the initial ion density of a cigar-shaped ablation region of length 4 cm is estimated to be n * =6.8×10 22 m -3 at the TAE location (r/a≈0.75). The frequency sweeping of an Alfvén eigenmode ends when the ion density homogenises poloidally. Modelling suggests that the time for poloidal homogenisation of the ion density at the TAE position is τ h =18±4 ms for inboard pellet injection, and τ h =26±2 ms for outboard pellet injection. By reproducing the frequency evolution of the elliptical Alfvén eigenmode (EAE), the initial ion density at the EAE location (r/a≈0.9) can be estimated to be n * =4.8×10 22 m -3 . Poloidal homogenisation of the ion density takes 2.7 times longer at the EAE location than at the TAE location for both inboard and outboard pellet injection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

HALO: A full-orbit model of nonlinear interaction of fast particles with eigenmodes

HALO (HAgis LOcust) solves the initial value Vlasov-Maxwell problem perturbatively for application to certain nonlinear wave-particle problems in tokamak plasmas. It uses the same basic approach as the HAGIS code for wave evolution but is built on the LOCUST-GPU full-orbit code for the solution of the Hamiltonian fast particle motion in cylindrical coordinates. The wave amplitude and particle evolution include all finite Larmour radius effects. We describe and benchmark the currently implemented Alfvén eigenmode workflow, demonstrating correct particle motion, linear and nonlinear power transfer. Here, the formulation and numerical scheme are sufficiently general as to allow easy future implementation of different kinds of eigenmodes, such as modes close to the ion-cyclotron frequency. The code can model multiple eigenmodes and multiple fast ion species simultaneously, and supports the general form of the equilibrium distribution in constants of motion.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2018

Force-free motion of a cold plasma during the current quench

Cold disruptive plasma tends to move during the current quench. Its motion is essentially force-free since the current quench timescale is resistive rather than Alfvénic. In contrast to the hot vertical displacement events, the frozen-in condition is violated in the cold plasma case, and the plasma motion is not governed by magnetic flux conservation but rather by its dissipation. In this paper, we present a numerical model of the cold plasma dynamics. This model predicts electromagnetic loads on the vacuum vessel, the plasma flow and density evolution, and the plasma centroid evolution. Our calculations include poloidal wall currents. Finally, we demonstrate their significant contribution to the force acting on the vacuum vessel.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Axisymmetric global Alfvén eigenmodes within the ellipticity-induced frequency gap in the Joint European Torus

Alfvén eigenmodes (AEs) with toroidal mode number n = 0 (i.e., axisymmetric) have been observed in the ellipticity-induced frequency range in the Joint European Torus. The axisymmetric modes are of interest because they can be used to diagnose fast particle energy distributions at the mode location. The modes were identified here as global Alfvén eigenmodes (GAEs), with the ellipticity of the plasma cross-section preventing strong continuum damping of the modes. The MHD codes CSCAS, MISHKA, and AEGIS were used to compute the n = 0 Alfvén continuum, eigenmode structure, and continuum damping. For zero ellipticity, a single mode exists at a frequency below the Alfvén continuum branch. This mode has two dominant poloidal harmonics with poloidal mode numbers m = ±1 that have the same polarity; therefore, it is an even mode. For finite ellipticity, the continuum branch splits into two branches and the single GAE splits into two modes. An even mode exists below the minimum of the top continuum branch, and the frequency of this mode coincides with the experimentally observed AE frequency. The other mode is found below the lower continuum branch with opposite signs between the two poloidal harmonics (an odd mode structure). This mode was not excited in our experiment. Analytical theory for the n = 0 GAE in an elliptical cylinder shows the n = 0 Alfvén continuum agrees with the numerical modelling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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 OSTI2015

Continuum absorption in the vicinity of the toroidicity-induced Alfvén gap

Excitation of Alfvén modes is commonly viewed as a concern for energetic particle confinement in burning plasmas. The 3.5 MeValpha particles produced by fusion may be affected as well as other fast ions in both present and future devices. Continuum damping of such modes is one of the key factors that determine their excitation thresholds and saturation levels. This work examines the resonant dissipative response of the Alfvén continuum to an oscillating driving current when the driving frequency is slightly outside the edges of the toroidicity-induced spectral gap. The problem is largely motivated by the need to describe the continuum absorption in the frequency sweeping events. Akey element of this problem is the negative interference of the two closely spaced continuum crossing points.Weexplain why the lower and upper edges of the gap can have very different continuum absorption features. Lastly, the difference is associated with an eigenmode whose frequency can be arbitrarily close to the upper edge of the gap whereas the lower edge of the gap is always a finite distance away from the closest eigenmode.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Overview of results from the MST reversed field pinch experiment

An overview of recent results from the MST reversed field pinch programme is presented here. With neutral beam injection, bursty energetic particle (EP) modes are observed. The profiles of the magnetic and density fluctuations associated with these EP modes are measured using a far infrared interferometer–polarimeter. Equilibrium reconstructions of the quasi-single-helicity 3D helical state are provided by the V3FIT code that now incorporates several of MST's advanced diagnostics. The orientation of the helical structure is controlled using a new resonant magnetic perturbation technique. Gyrokinetic simulations based on experimental equilibria predict unstable trapped-electron modes (TEMs), and small-scale density fluctuations are detected in improved-confinement plasmas with TEM-like features. Upgraded pellet injection permits study of density and beta limits over MST's full range of operation, and an MST-record line-average density of 0.9 × 1020 m3 (n/nG = 1.4) has been obtained. Impurity ion temperature measurements reveal a charge-to-mass-ratio dependence in the rapid heating that occurs during a sawtooth crash. Runaway of NBI-born fast ions during the impulsive sawtooth event agrees with test-particle theory. Magnetic self-organization studies include measurements of the dynamo emf with an applied ac inductive electric field using oscillating field current drive.

beta limit↗
DOE OSTI2014

Plasma Energetic Particles Simulation Center (PEPSC)

The main effort of the Texas group was to develop theoretical and simplified numerical models to understand chirping phenomena often seen for Alfven and geodesic acoustic waves in experimental plasmas such as D-III-D, NSTX and JET. Its main numerical effort was to modify the AEGIS code, which was originally developed as an eigenvalue solver. To apply to the chirping problem this code has to be able to treat the linear response to the continuum and the response of the plasma to external drive or to an internal drive that comes from the formation of phase space chirping structures. The theoretical underpinning of this investigation still needed to be more fully developed to understand how to best formulate the theoretical problem. Considerable progress was made on this front by B.N. Breizman and his collaborators and a new reduced model was developed by H. L. Berk and his PhD student, G. Wang which can be uses as simplified model to describe chirping in a large aspect ratio tokamak. This final report will concentrate on these two directions that were developed as well as results that were found in the work with the AEGIS code and in the progress in developing a novel quasi-linear formulation for a description of Alfvenic modes destabilized by energetic particles, such as alpha particles in a burning plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Nonlinearly driven harmonics of Alfvén modes

In order to study the leading order nonlinear magneto-hydrodynamic (MHD) harmonic response of a plasma in realistic geometry, the AEGIS code has been generalized to account for inhomogeneous source terms. These source terms are expressed in terms of the quadratic corrections that depend on the functional form of a linear MHD eigenmode, such as the Toroidal Alfvén Eigenmode. The solution of the resultant equation gives the second order harmonic response. Preliminary results are presented here.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Modeling of long range frequency sweeping for energetic particle modes

Long range frequency sweeping events are simulated numerically within a one-dimensional, electrostatic bump-on-tail model with fast particle sources and collisions. The numerical solution accounts for fast particle trapping and detrapping in an evolving wave field with a fixed wavelength, and it includes three distinct collisions operators: Drag (dynamical friction on the background electrons), Krook-type collisions, and velocity space diffusion. The effects of particle trapping and diffusion on the evolution of holes and clumps are investigated, and the occurrence of non-monotonic (hooked) frequency sweeping and asymptotically steady holes is discussed. The presented solution constitutes a step towards predictive modeling of frequency sweeping events in more realistic geometries.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗