No Arabic abstract
The aim of this study is to analyze the stability of the Alfven eigenmodes (AE) in the Chinese First Quasi-axisymmetric Stellarator (CFQS). The AE stability is calculated using the code FAR3d that solves the reduced MHD equations to describe the linear evolution of the poloidal flux and the toroidal component of the vorticity in a full 3D system, coupled with equations of density and parallel velocity moment for the energetic particles (EP) species including the effect of the helical couplings and acoustic modes. The Landau damping and resonant destabilization effects are added in the model by a given closure relation. The simulation results indicate the destabilization of n = 1 to 4 AEs by EP during the slowing down process, particularly n = 1 and n = 2 toroidal AEs (TAE), n = 3 elliptical AE (EAE) and n = 4 non circular AE (NAE). If the resonance is caused by EPs with an energy above 17 keV (weakly thermalized EP), n = 2 EAEs and n = 3 NAEs are unstable. On the other hand, EPs with an energy below 17 keV (late thermalization stage) lead to the destabilization of n = 3 and n = 4 TAEs. The simulations for an off-axis NBI injection indicate the further destabilization of n = 2 to 4 AEs although the growth rate of the n = 1 AEs slightly decreases, so no clear optimization trend with respect to the NBI deposition region is identified. In addition, n = 2, 4 helical AE (HAE) are unstable above an EP b{eta} threshold. Also, if the thermal b{eta} of the simulation increases (higher thermal plasma density) the AE stability of the plasma improves. The simulations including the effect of the finite Larmor radius and electron-ion Landau damping show the stabilization of the n = 1 to 4 EAE/NAEs as well as a decrease of the growth rate and frequency of the n = 1 to 4 BAE/TAEs.
Alfven Eigenmodes (AE) can be destabilized during ITER discharges driven by neutral beam injection (NBI) energetic particles (EP) and alpha particles. The aim of the present study is to analyze the AE stability of different ITER operation scenarios considering multiple energetic particle species. We use the reduced magneto-hydrodynamic (MHD) equations to describe the linear evolution of the poloidal flux and the toroidal component of the vorticity in a full 3D system, coupled with equations of density and parallel velocity moments for the energetic particles species including the effect of the acoustic modes. The AEs driven by the NBI EP and alpha particles are stable in the configurations analyzed, only MHD-like modes with large toroidal couplings are unstable, although both can be destabilized if the EP beta increases above a threshold. The threshold is two times the model beta value for the NBI EP and alpha particles in the reverse shear case, leading to the destabilization of Beta induced AE (BAE) near the magnetic axis with a frequency of 25-35 kHz and Toroidal or Elliptical AE (TAE/EAE) in the reverse shear region with a frequency of 125-175 kHz, respectively. On the other hand, the hybrid and steady state configurations show a threshold 3 times larger with respect to the model beta for the alpha particle and 40 times for the NBI EP, also destabilizing BAE and TAE between the inner and middle plasma region. In addition, a extended analysis of the reverse shear scenario where the beta of both alpha particles and NBI EP are above the AE threshold, multiple EP damping effects are also identified as well as optimization trends regarding the resonance properties of the alpha particle and NBI EP with the bulk plasma.
The aim of this study is to analyze the destabilization of Alfven Eigenmodes (AE) by multiple energetic particles (EP) species in DIII-D and LHD discharges. We use the reduced MHD equations to describe the linear evolution of the poloidal flux and the toroidal component of the vorticity in a full 3D system, coupled with equations of density and parallel velocity moments for the energetic particles species, including the effect of the acoustic modes, diamagnetic currents and helical couplings. We add the Landau damping and resonant destabilization effects using a closure relation. The simulations with multiple NBI lines show three different regimes: the non damped regime where the multi beam AEs growth rate is larger compared to the growth rate of the AEs destabilized by the individual NBI lines, the interaction regime where the multi beam AEs growth rate is smaller than the single NBI AEs and the damped regime where the AEs are suppressed. Operations in the damped regime requires EP species with different density profile flatness or gradient locations. In addition, the AEs growth rate in the interaction regime is further reduced if the combined NBI lines have similar beam temperatures and the beta of the NBI line with flatter EP density profile increases. Then, optimization trends are identified in DIII-D high poloidal beta and LHD low density / magnetic field discharges with multiple NBI lines as well as the configuration requirements to operate in the damped and interaction regimes. DIII-D simulations show a decrease of the n=2 to 6 AEs growth rate and n=1 AE are stabilized in the LHD case. The helical coupling effects in LHD simulations lead to a transition from the interaction to the damped regime of the n=2,-8,12 helical family.
The aim of this study is to perform a theoretical analysis of the magnetohydrodynamic (MHD) stability and energetic particle effects on a LHD equilibria, calculated during a discharge where energetic-ion-driven resistive interchange mode (EIC) events were triggered. We use the reduced MHD equations to describe the linear evolution of the poloidal flux and the toroidal component of the vorticity in a full 3D system, coupled with equations of density and parallel velocity moments for the energetic particles species. We add the Landau damping and resonant destabilization effects using a closure relation. The simulations suggest that the helically trapped EP driven by the perpendicular neutral beam injector (NBI) further destabilizes the 1/1 MHD-like mode located at the plasma periphery (r/a = 0.88). If the beta of the EP driven by the perpendicular NBI is larger than 0.0025 a 1/1 EIC with a frequency around 3 kHz is destabilized. If the effect of the passing EP driven by the tangential NBI is included on the model, any enhancement of the injection intensity of the tangential NBI below beta=0.025 leads to a decrease of the instability growth rate. The simulations indicate that the perpendicular NBI EP is the main driver of the EIC events, as it was observed in the experiment. If the effect of the helical couplings are added in the model, an 11/13 EIC is destabilized with a frequency around 9 kHz, inward shifted (r/a = 0.81) compared to the 1/1 EIC. Thus, one possible explanation for the EIC frequency chirping down from 9 to 3 kHz is a transition between the 11/13 to the 1/1 EIC due to a weakening of the destabilizing effect of the high n modes, caused by a decrease of the EP drive due to a loss of helically trapped EP or a change in the EP distribution function after the EIC burst.
Recent upgrades in H-1 power supplies have enabled the operation of the H-1 experiment at higher heating powers than previously attainable. A heating power scan in mixed hydrogen/helium plasmas reveals a change in mode activity with increasing heating power. At low power (<50 kW) modes with beta-induced Alfven eigenmode (BAE) frequency scaling are observed. At higher power modes consistent with an analysis of nonconventional Global Alfven Eigenmodes (GAEs) are observed, the subject of this work. We have computed the mode continuum, and identified GAE structures using the ideal MHD solver CKA and the gyrokinetic code EUTERPE. An analytic model for ICRH-heated minority ions is used to estimate the fast ion temperature from the hydrogen species. Linear growth rate scans using a local flux surface stability calculation, LGRO, are performed. These studies demonstrate growth from circulating particles whose speed is significantly less than the Alfven speed, and are resonant with the mode through harmonics of the Fourier decomposition of the strongly-shaped heliac magnetic field. They reveal drive is possible with a small, hot energetic tail of the hydrogen species. Local linear growth rate scans are also complemented with global calculations from CKA and EUTERPE. These qualitatively confirm the findings from the LGRO study, and show that the inclusion of finite Larmor radius effects can reduce the growth rate by a factor of three, but do not affect marginal stability. Finally, a study of damping of the global mode with the thermal plasma is conducted, computing continuum, and the damping arising from parallel electric fields. We find that continuum damping is of order 0.1% for the configuration studied. The inclusion of resistivity lifts the damping to 19%. Such large damping is consistent with experimental observations that in absence of drive the mode decays rapidly (~0.1 ms).
Spontaneous nonlinear excitation of geodesic acoustic mode (GAM) by toroidal Alfven eigenmode (TAE) is investigated using nonlinear gyrokinetic theory. It is found that, the nonlinear decay process depends on thermal ion beta value. Here, beta is the plasma thermal to magnetic pressure ratio. In the low-beta limit, TAE decays into a GAM and a lower TAE sideband in the toroidicity induced shear Alfven wave continuous spectrum gap; while in the high-beta limit, TAE decays into a GAM and a propagating kinetic TAE in the continuum. Both cases are investigated for the spontaneous decay conditions. The nonlinear saturation levels of both GAM and daughter wave are derived. The corresponding power balance and wave particle power transfer to thermal plasma are computed. Implications on thermal plasma heating are also discussed.