No Arabic abstract
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).
Linear and nonlinear modelling of Alfvenic instabilities, most notably toroidal Alfven eigenmodes (TAEs), obtained by using the global nonlinear electromagnetic gyrokinetic model of the code ORB5 are presented for the 15 MA scenario of the ITER tokamak. Linear simulations show that elliptic Alfven eigenmodes and odd-parity TAEs are only weakly damped but not excited by alpha particles, whose drive favours even-parity TAEs. Low mode number TAEs are found to be global, requiring global treatment. Nonlinearly, even with double the nominal EP density, single mode simulations lead to saturation with negligible EP transport however multi-mode simulations predict that with double the nominal EP density, enhanced saturation and significant EP redistribution will occur.
Nonlinear generation of high frequency mode (HFM) by toroidal Alfven eigenmode (TAE) observed in HL-2A tokamak is analyzed using nonlinear gyrokinetic theory. It is found that, the HFM can be dominated by $|nq-m|=1$ perturbations with predominantly ideal magnetohydrodynamic if the two primary TAEs are co-propagating; while the HFM can be characterized by $nq-m=0$ electrostatic perturbations if the two primary TAEs are counter-propagating. Here, $n$ and $m$ are respectively the toroidal and poloidal mode numbers, and $q$ is the safety factor. The nonlinear process is sensitive to the equilibrium magnetic geometry of the device.
Gyrokinetic theory of nonlinear mode coupling as a mechanism for toroidal Alfven eigenmode (TAE) saturation in the fusion plasma related parameter regime is presented, including 1) para- metric decay of TAE into lower kinetic TAE (LKTAE) and geodesic acoustic mode (GAM), and 2) enhanced TAE coupling to shear Alfven wave (SAW) continuum via ion induced scattering. Our theory shows that, for TAE saturation in the parameter range of practical interest, several processes with comparable scattering cross sections can be equally important.
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.
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.