Nonlinear simulation of multiple toroidal Alfvén eigenmodes in tokamak plasmas

Nonlinear evolution of multiple toroidal Alfvén eigenmodes(TAEs)driven by fast ions is self-consistently investigated by kinetic simulations in toroidal plasmas.To clearly identify the effect of nonlinear coupling on the beam ion loss,simulations over single-n modes are also carried out and compared with those over multiple-n modes,and the wave-particle resonance and particle trajectory of lost ions in phase space are analyzed in detail.It is found that in the multiple-n case,the resonance overlap occurs so that the fast ion loss level is rather higher than the sum loss level that represents the summation of loss over all single-n modes in the single-n case.Moreover,increasing fast ion betaβh can not only significantly increase the loss level in the multiple-n case but also significantly increase the loss level increment between the single-n and multiple-n cases.For example,the loss level in the multiple-n case forβh=6.0%can even reach 13%of the beam ions and is 44%higher than the sum loss level calculated from all individual single-n modes in the single-n case.On the other hand,when the closely spaced resonance overlap occurs in the multiple-n case,the release of mode energy is increased so that the widely spaced resonances can also take place.In addition,phase space characterization is obtained in both single-n and multiple-n cases.

Numerical simulation of the multiple reversed shear Alfvén eigenmodes associated with the triangularity Alfvén gap

It was found that there are multiplicity of low shear toroidicity-induced Alfv′en eigenmodes in a zero beta limit if the inverse aspect ratio is larger than the magnetic shear at the mode location(Candy 1996 Phys. Lett. A 215 299). Because the reversed shear Alfv′en eigenmode(RSAE) and even the RSAE associated with the non-circular triangularity-induced Alfv′en eigenmode(NAE) gap(NAE–RSAE) usually reside near the shear-reversal point, the condition that the inverse aspect ratio is larger than the magnetic shear is naturally satisfied. For this reason, we numerically investigate the existence of multiplicity of core-localized NAE–RSAEs and mode characteristics in the present work. We firstly verify the existence of the multiplicity for zero beta plasma by using a D-shaped equilibrium. It is pointed out that, for a given toroidal mode number, the Alfv′en cascade spectrum accommodates down-sweeping and up-sweeping modes above and below the NAE range of frequencies. An analytical model for the existence of multiple RSAE modes is in good agreement with the simulation results. One notices that the triangularity has a greater effect on the odd-type modes than that on the even-type modes: the odd-type modes come into existence because of the plasma triangularity.

Alfvén eigenmode stability analysis and energetic particle transport prediction for CFETR hybrid scenario

The hybrid scenario is a projection for CFETR operation with high plasma current and density.Therefore, the energetic particles(EPs) generated by fusion reactions can destabilize Alfvén eigenmodes(AEs), which could result in significant EPs loss and redistribution. Both the eigenvalue code NOVA-K and the wrapped local stability code TGLFEP are used to analyze AE stability. The simulation indicates the beta-induced Alfvén eigenmodes with n?>?5 in the core region are the most unstable. The NOVA-K code is used to benchmark the critical density gradient calculated by TGLFEP. In addition, the EPtran code is employed to predict EP transport induced by unstable AEs and turbulence, which reduce EP density in the core and drive approximately 30% EP transport from the core to the edge, thus the EP density profile flattens and EPs with lower energy deposit near the edge.

Transition from fishbone mode toβ-induced Alfvén eigenmode on HL-2A tokamak

In the presence of energetic particles(EPs)from auxiliary heating and burning plasmas,fishbone instability and Alfvén modes can be excited and their transition can take place in certain overlapping regimes.Using the hybrid kinetic-magnetohydrodynamic model in the NIMROD code,we have identified such a transition between the fishbone instability and theβ-induced Alfvén eigenmode(BAE)for the NBI heated plasmas on HL-2 A.When the safety factor at magnetic axis is well below one,typical kink-fishbone transition occurs as the EP fraction increases.When q0 is raised to approaching one,the fishbone mode is replaced with BAE for sufficient amount of EPs.When q0 is slightly above one,the toroidicity-induced Alfvén eigenmode dominates at lower EP pressure,whereas BAE dominates at higher EP pressure.

Preliminary computation of the gap eigenmode of shear Alfvén waves on the LAPD

Characterizing the gap eigenmode of shear Alfv′en waves(SAWs) and its interaction with energetic ions is important to the success of magnetically confined fusion. Previous studies have reported an experimental observation of the spectral gap of SAW on the on Large Plasma Device(LAPD)(Zhang et al. 2008 Phys. Plasmas 15 012103), a linear large plasma device(Gekelman et al. 1991 Rev. Sci. Instrum. 62 2875) possessing easier diagnostic access and lower cost compared with traditional fusion devices, and analytical theory and numerical gap eigenmode using ideal conditions(Chang 2014 Ph.D Thesis at Australian National University). To guide experimental implementation, the present work models the gap eigenmode of SAWs using exact LAPD parameters. A full picture of the wave field for previous experiment reveals that the previously observed spectral gap is not global but an axially local result. To form a global spectral gap, the number of magnetic mirrors has to be increased and stronger static magnetic field makes it clearer. Such a spectral gap is obtained for the magnetic field of B_0(z) = 1.2 + 0.6 cos[2π(z-33.68)/3.63] with 7.74-m magnetic beach. By introducing two types of local defects(corresponding to E_θ(z_0) = 0 and E'_θ(z_0) = 0 respectively), odd-parity and even-parity discrete eigenmodes are formed clearly inside the gap. The strength of these gap eigenmodes decreases significantly with collision frequency, which is consistent with previous studies. Parameter scans show that these gap eigenmodes can be even formed successfully for the field strength of B_0(z) = 0.2 + 0.1 cos[2π(z-33.68)/3.63] and with only four magnetic mirrors, which are achievable by the LAPD at its present status. This work can serve as a strong motivation and direct reference for the experimental implementation of the gap eigenmode of SAWs on the LAPD and other linear plasma devices.

Mode structure symmetry breaking of reversed shear Alfvén eigenmodes and its impact on the generation of parallel velocity asymmetries in energetic particle distribution

In this work,the gyrokinetic eigenvalue code LIGKA,the drift-kinetic/MHD hybrid code HMGC and the gyrokinetic full-f code TRIMEG-GKX are employed to study the mode structure details of reversed shear Alfvén eigenmodes(RSAEs).Using the parameters from an ASDEXUpgrade plasma,a benchmark with the three different physical models for RSAE without and with energetic particles(EPs)is carried out.Reasonable agreement has been found for the mode frequency and the growth rate.Mode structure symmetry breaking(MSSB)is observed when EPs are included,due to the EPs’non-perturbative effects.It is found that the MSSB properties are featured by a finite radial wave phase velocity,and the linear mode structure can be well described by an analytical complex Gaussian expressionФ(s)=e^(-σ(s-s_(0))^(2))with complex parametersσand s_(0),where s is the normalized radial coordinate.The mode structure is distorted in opposite manners when the EP drive shifted from one side of qminto the other side,and specifically,a non-zero average radial wave number with opposite signs is generated.The initial EP density profiles and the corresponding mode structures have been used as the input of HAGIS code to study the EP transport.The parallel velocity of EPs is generated in opposite directions,due to different values of the average radial wave number,corresponding to different initial EP density profiles with EP drive shifted away from the qmin.

Linear gyrokinetic simulations of reversed shear Alfvén eigenmodes and ion temperature gradient modes in DIII-D tokamak

Global linear gyrokinetic simulations using realistic DIII-D tokamak geometry and plasma profiles find co-existence of unstable reversed shear Alfvén eigenmodes(RSAE)with low toroidal mode number n and electromagnetic ion temperature gradient(ITG)instabilities with higher toroidal mode number n.For intermediate n?=?[10,12],RSAE and ITG co-exist and overlap weakly in the radial domain with similar growth rates but different real frequencies.Both RSAE and ITG growth rates decrease less than 5%when compressible magnetic perturbations are neglected in the simulations.The ITG growth rates increase less than 7%when fast ions are not included in the simulations.Finally,the effects of trapped electrons on the RSAE are negligible.

Nonlinear excitation of a geodesic acoustic mode by reversed shear Alfvén eignemodes

The parametric decay process of a reversed shear Alfvén eigenmeode(RSAE)into a geodesic acoustic mode and a kinetic RSAE is investigated using nonlinear gyrokinetic theory.The excitation conditions mainly require the pump RSAE amplitude to exceed a certain threshold,which could be readily satisfied in burning plasmas operated in steady-state advanced scenario.This decay process can contribute to thermal plasma heating and confinement improvement.

Runaway electron dynamics in Experimental Advanced Superconducting Tokamak helium plasmas

The generation of runaway electrons(REs)is observed during the low-density helium ohmic plasma discharge in the Experimental Advanced Superconducting Tokamak(EAST).The growth rate of hard x-ray(HXR)is inversely proportional to the line-average density.Besides,the RE generation in helium plasma is higher than that in deuterium plasma at the same density,which is obtained by comparing the growth rate of HXR with the same discharge conditions.The potential reason is the higher electron temperature of helium plasma in the same current and electron density plateau.Furthermore,two Alfvén eigenmodes driven by REs have been observed.The frequency evolution of the mode is not fully satisfied with the Alfvén scaling and when extension of the Alfvén frequency is towards 0,the high frequency branch is~50 kHz.The different spatial position of the two modes and the evolution of the helium concentration could be used to understand deviation between theoretical and experimental observation.

Formation of the mass density peak at the magnetospheric equator triggered by EMIC waves

We report a simultaneous observation of two band electromagnetic ion cyclotron(EMIC)waves and toroidal Alfvén waves by the Van Allen Probe mission.Through wave frequency analyses,the mass densityρis found to be locally peaked at the magnetic equator.Perpendicular fluxes of ions(<100 eV)increase simultaneously with the appearances of EMIC waves,indicating a heating of these ions by EMIC waves.In addition,the measured ion distributions also support the equatorial peak formation,which accords with the result of the frequency analyses.The formation of local mass density peaks at the equator should be due to enhancements of equatorial ion concentrations,which are triggered by EMIC waves’perpendicular heating on low energy ions.