-
Effect of the neutral beam injector operational regime on the Alfven eigenmode saturation phase in DIII-D plasma
Authors:
J Varela,
D A Spong,
L Garcia,
Y Ghai,
D Zarzoso,
D del-Castillo-Negrete,
H Betar,
J Ortiz,
D C Pace,
M A Van Zeeland,
X Du,
R Sanchez,
V Tribaldos,
J M Reynolds-Barredo
Abstract:
The aim of this study is to analyze the effect of the neutral beam injector (NBI) operation regime on the saturation phase of the Alfven Eigenmodes (AEs) in DIII-D plasma. The analysis is done using the linear and nonlinear versions of the gyro-fluid code FAR3d. A set of parametric analyses are performed modifying the nonlinear simulation EP \b{eta} (NBI injection power), EP energy (NBI voltage) a…
▽ More
The aim of this study is to analyze the effect of the neutral beam injector (NBI) operation regime on the saturation phase of the Alfven Eigenmodes (AEs) in DIII-D plasma. The analysis is done using the linear and nonlinear versions of the gyro-fluid code FAR3d. A set of parametric analyses are performed modifying the nonlinear simulation EP \b{eta} (NBI injection power), EP energy (NBI voltage) and the radial location of the EP density profile gradient (NBI radial deposition). The analysis indicates a transition from the soft (local plasma relaxation) to the hard MHD (global plasma relaxation) limit if the simulation EP \b{eta} >= 0.02, leading to bursting MHD activity caused by radial AEs overlapping. MHD bursts cause an enhancement of the EP transport showing ballistic-like features as avalanche-like events. Simulations in the soft MHD limit show an increment of the EP density gradient as the EP \b{eta} increases. On the other hand, there is a gradient upper limit in the hard MHD limit, consistent with the critical-gradient behavior. AEs induce shear flows and zonal current leading to the deformation of the flux surfaces and the safety factor profile, respectively, particularly strong for the simulation in the hard MHD limit. Simulations in the hard MHD regime show a decrease of the AE frequency in the saturation phase; this is caused by the destabilization of a transitional mode between a 9/3-10/3 TAE and a 9/3 RSAE that may explain the AE frequency down-sweeping observed in some DIII-D discharges. Reducing the EP energy in the nonlinear simulations leads to a weakening of the plasma perturbation. On the other hand, increasing the EP energy causes the opposite effect. Nonlinear simulations of off-axis NBI profiles indicate a lower plasma perturbation as the EP density gradient is located further away from the magnetic axis.
△ Less
Submitted 7 November, 2023;
originally announced November 2023.
-
Analysis of the ECH effect on the EPM/AEs stability in Heliotron J plasma using a Landau closure model
Authors:
J. Varela,
K. Nagasaki,
S. Kobayashi,
K. Nagaoka,
P. Adulsiriswad,
A. Cappa,
S. Yamamoto,
K. Y. Watanabe,
D. A. Spong,
L. Garcia,
Y. Ghai,
J. Ortiz
Abstract:
The aim of the present study is to analyze the effect of the electron cyclotron heating (ECH) on the linear stability of Alfven Eigenmodes (AE) and energetic particle modes (EPM) triggered by energetic ions in Heliotron J plasma. The analysis is performed using the FAR3d code that solves a reduced MHD model to describe the thermal plasma coupled with a gyrofluid model for the energetic particles (…
▽ More
The aim of the present study is to analyze the effect of the electron cyclotron heating (ECH) on the linear stability of Alfven Eigenmodes (AE) and energetic particle modes (EPM) triggered by energetic ions in Heliotron J plasma. The analysis is performed using the FAR3d code that solves a reduced MHD model to describe the thermal plasma coupled with a gyrofluid model for the energetic particles (EP) species. The simulations reproduce the AE/EPM stability trends observed in the experiments as the electron temperature (Te) increases, modifying the thermal plasma beta, EP beta and EP slowing down time. Particularly, the n/m=1/2 EPM and 2/4 Global AE (GAE) are stabilized in the low bumpiness (LB) configuration due to an enhancement of the continuum, Finite Larmor radius (FLR) and e-i Landau damping effects as the thermal beta increases. On the other hand, a larger ECH injection power cannot stabilize the AE/EPM in Medium (MB) and High bumpiness (HB) configurations because the damping effects are weaker compared to the LB case, unable to balance the further destabilization induced by an enhanced EP resonance as the EP slowing down time and EP beta increases with Te
△ Less
Submitted 18 April, 2023;
originally announced April 2023.
-
Simulation of the TAEs saturation phase in Large Helical Device device: MHD burst
Authors:
J. Varela,
D. A. Spong,
Y. Todo,
L. Garcia,
Y. Ghai,
J. Ortiz,
R. Seki
Abstract:
The aim of the present study is to analyze the saturation regime of the Toroidal Alfven Eigenmodes (TAE) in the LHD plasma, particularly the MHD burst. The linear and nonlinear evolution of the TAEs are simulated by the FAR3d code that uses a reduced MHD model for the thermal plasma coupled with a gyrofluid model for the energetic particles (EP) species. The linear simulations indicate the overlap…
▽ More
The aim of the present study is to analyze the saturation regime of the Toroidal Alfven Eigenmodes (TAE) in the LHD plasma, particularly the MHD burst. The linear and nonlinear evolution of the TAEs are simulated by the FAR3d code that uses a reduced MHD model for the thermal plasma coupled with a gyrofluid model for the energetic particles (EP) species. The linear simulations indicate the overlapping of 1/2-1/1, 2/3-2/4 and 3/5-3/6 TAEs in the inner-middle plasma region and frequency range of 45-75 kHz, triggered by EPs with an energy of Tf = 45 keV and EP beta = 0.022. The nonlinear simulations show that 2/3-2/4 and 3/4-3/5 TAEs are further destabilized due to the energy transfer from 1/1-1/2 TAE, leading to a broad TAEs radial overlapping and the MHD burst triggering. The energy of 1/1-1/2 TAE is also nonlinearly transferred to the thermal plasma destabilizing the 0/0 and 0/1 modes, inducing the generation of shear flows and zonal currents as well as large deformations in the thermal pressure and EP density radial profiles. The nonlinear simulation reproduces the same succession of instabilities and the same frequency range with respect to the experiment. The instability propagates outward during the bursting phase, showing a large decrease of the EP density profile between the middle-outer plasma, pointing out the loss of part of the EP population that explains the decrease of the plasma heating efficiency observed during the MHD burst.
△ Less
Submitted 17 April, 2023;
originally announced April 2023.
-
Theoretical study of the Alfven Eigenmode stability in CFETR steady state discharges
Authors:
J. Varela,
J. Huang,
D. A. Spong,
J. Chen,
V. Chan,
L. Garcia,
A. Wingen,
Y. Ghai,
Y. Zou
Abstract:
The aim of this study is to analyze the stability of Alfven Eigenmodes (AE) in the China Fusion Engineering Test Reactor (CFETR) plasma for steady state operations. The analysis is done using the gyro-fluid code FAR3d including the effect of the acoustic modes, EP Finite Larmor radius damping effects and multiple energetic particle populations. Two high poloidal beta scenarios are studied with res…
▽ More
The aim of this study is to analyze the stability of Alfven Eigenmodes (AE) in the China Fusion Engineering Test Reactor (CFETR) plasma for steady state operations. The analysis is done using the gyro-fluid code FAR3d including the effect of the acoustic modes, EP Finite Larmor radius damping effects and multiple energetic particle populations. Two high poloidal beta scenarios are studied with respect to the location of the internal transport barrier (ITB) at r/a = 0.45 (case A) and r/a = 0.6 (case B). Both operation scenarios show a narrow TAE gap between the inner-middle plasma region and a wide EAE gap all along the plasma radius. The AE stability of CFETR plasmas improves if the ITB is located inwards, case A, showing AEs with lower growth rates with respect to the case B. The AEs growth rate is smaller in the case A because the modes are located in the inner-middle plasma region where the stabilizing effect of the magnetic shear is stronger with respect to the case B. Multiple EP populations effects (NBI driven EP + alpha articles) are negligible for the case A, although the simulations for the case B show a stabilizing effect of the NBI EP on the n=1 BAE caused by alpha particles during the thermalization process. If the FLR damping effects are included in the simulations, the growth rate of the EAE/NAE decreases up to 70 %, particularly for n > 3 toroidal families. Low n AEs (n<6) show the largest growth rates. On the other hand, high $n$ modes (n=6 to 15) are triggered in the frequency range of the NAE, strongly damped by the FLR effects.
△ Less
Submitted 17 April, 2023;
originally announced April 2023.
-
Theoretical analysis of the saturation phase of the $1/1$ energetic-ion-driven resistive interchange mode
Authors:
J. Varela,
D. A. Spong,
L. Garcia,
S. Ohdachi,
K. Y. Watanabe,
R. Seki,
Y. Ghai
Abstract:
The aim of the present study is to analyze the saturation regime of the energetic-ion-driven resistive interchange mode (EIC) in the LHD plasma. A set of non linear simulations are performed by the FAR3d code that uses a reduced MHD model for the thermal plasma coupled with a gyrofluid model for the energetic particles (EP) species. The hellically trapped EP component is introduced through a modif…
▽ More
The aim of the present study is to analyze the saturation regime of the energetic-ion-driven resistive interchange mode (EIC) in the LHD plasma. A set of non linear simulations are performed by the FAR3d code that uses a reduced MHD model for the thermal plasma coupled with a gyrofluid model for the energetic particles (EP) species. The hellically trapped EP component is introduced through a modification of the averaged drift velocity operator to include their precessional drift. The non linear simulation results show similar 1/1 EIC saturation phases with respect to the experimental observations, reproducing the enhancement of the n/m = 1/1 resistive interchange modes (RIC) amplitude and width as the EP $β$ increases, the EP beta threshold for the 1/1 EIC excitation, the further destabilization of the 1/1 EIC as the population of the helically trapped EP increases and the triggering of burst events. The frequency of the 1/1 EIC calculated during the burst event is 9.4 kHz and the 2/2 and 3/3 overtones are destabilized, consistent with the frequency range and the complex mode structure measured in the experiment. In addition, the simulation shows the inward propagation of the 1/1 EIC due to the non linear destabilization of the 3/4 and 2/3 EPMs, leading to the partial overlapping between resonances during the burst event. Finally, the analysis of the 1/1 EIC stabilization phase shows the excitation of the 1/1 RIC as soon as the flattening induced by the 1/1 EIC in the pressure profile vanishes, leading to the retrieval of the pressure gradient at the plasma periphery and the overcoming of the RIC stability limit.
△ Less
Submitted 3 January, 2022;
originally announced January 2022.
-
Study of the Alfven Eigenmodes stability in CFQS plasma using a Landau closure model
Authors:
J. Varela,
A. Shimizu,
D. A. Spong,
L. Garcia,
Y. Ghai
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 velo…
▽ More
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.
△ Less
Submitted 16 April, 2021;
originally announced April 2021.