Since 18 of December 2019 conferences.iaea.org uses Nucleus credentials. Visit our help pages for information on how to Register and Sign-in using Nucleus.

# 28th IAEA Fusion Energy Conference (FEC 2020)

May 10 – 15, 2021
Virtual Event
Europe/Vienna timezone
The Conference will be held virtually from 10-15 May 2021

## Eﬀects of core plasma on the low frequency Alfven and Acoustic eigenmodes

May 13, 2021, 8:30 AM
4h
Virtual Event

#### Virtual Event

Regular Poster Magnetic Fusion Theory and Modelling

### Speaker

Ilija Chavdarovski (National Fusion Research Institute)

### Description

Low frequency modes belonging to the Alfvenic Kinetic Ballooning mode (KBM)[1] and the mixed polarization Beta-induced Alfven Acoustic Eigenmode (BAAE)[2] branches, have been recently observed experimentally [3] and confirmed numerically [4]. Due to the low frequency range of the core ion bounce motion, kinetic treatment of both circulating[5,6] and trapped particles[7,8] is required in order to properly describe these fluctuations. We give the dispersion relation of the aforementioned modes in the framework of the Generalized fishbone-like dispersion relation (GFLDR)[1,9]

$i \Lambda(\omega) = \delta \bar{W}_f + \delta \bar{W}_k \;\; ,$

which is a unifying theoretical description for various Alfvenic fluctuations, as well as Energetic particle continuum modes (EPMs). Here, $\Lambda$ is the generalized inertia representing the physics on short radial scales, while $\delta \bar{W}_f$ and $\delta \bar{W}_k$ are the ideal region fluid and kinetic contribution, respectively.
We assume low $\beta=O(\epsilon^{2})$ axisymmetric $(s,\alpha)$ plasma equilibrium and solve the vorticity equation:

$B \textbf{b} \cdot \nabla \left[{1 \over B } {{k_\perp} ^2\over {k_\vartheta} ^2} \textbf{b}\cdot \nabla \delta\psi \right] + {\omega^2 \over v_A^2} \left(1- {\omega_{\ast pi} \over \omega }\right) \frac{k_\perp^2}{k_\vartheta^2} \delta\phi + {\alpha \over q^2 R^2} g (\theta)\delta \psi = \nonumber \\ \Bigl \langle \frac{4\pi e}{k_\vartheta^2 c^2} \omega \omega_{di} \delta K_i \Bigr \rangle \, ,$

and quasi-neutrality condition, in the long wavelength limit [1,5,7]:

$\left( 1+\frac{T_i}{T_e}\right)(\delta\phi - \delta\psi) = \frac{T_i}{ne}\langle \delta K_i - \delta K_e \rangle \, .$

Here $\langle (...)\rangle=\int d\textbf{v}(...)$, $\delta K_{e,i}$ are the particle non-adiabatic distribution functions, $g(\theta)=\cos \theta +[s\theta-\alpha \sin \theta \,] \sin \theta$ and $n_e = n_i = n$. Expanding the fields $\delta \Phi=(k_\perp/k_\theta) \delta \phi$ and $\delta \Psi=(k_\perp/k_\theta) \delta \psi$ in asymptotic series in powers of $\beta^{1/2}$[5] we obtain, in zeroth order of the quasineutrality condition, $\delta \Phi^{(0)}=I_\Phi(\omega,\overline{\omega}_{Di,e},\omega_{*i,e}) \delta \Psi^{(0)}$. Here, $I_\Phi\simeq1$ for most of the Alfvenic spectrum (except frequencies near precession resonance $\overline{\omega}_{Di,e}$), consistently with the ideal MHD limit $\delta E_\parallel=0$. Further expansion of the quasineutrality condition gives $\delta \Phi_s= S(\omega,\omega_{Bi},\overline{\omega}_{Di,e},\omega_{Ti}) \xi \delta \Phi^{(0)}$, where $\xi\simeq k_\perp/k_\theta$ and $\delta \Phi_s \sim O(\beta^{1/2})$ is the $\sim \sin \theta$ modulation of the potential along magnetic field line, which makes the function $|S|$ a measure of how much the mode polarization deviates from pure Alfv\'enic due to the parallel a.c. electric field.

The generalized inertia is obtained from the vorticity equation expanded to $O(\beta)$[7]:

$\Lambda^2/I_\Phi = \frac{\omega^2}{\omega_A^2} \left ( 1- {\omega_{\ast pi} \over \omega} \right)+ \Lambda^2_{cir} +\Lambda^2_{tra} \; ,$

where $\Lambda^2_{cir}$[5] and $\Lambda^2_{tra}$[7] are the circulating and trapped particles contributions, respectively. Even though this expression of $\Lambda$ is a mixture of deeply trapped and well circulating particle responses, it has been shown[7] that this reduced model recovers well the low and high frequency limits of the continuous spectrum, and further gives good insights into experimental results. The term $I_\phi$ acts as an additional inertia enhancement due to the opposite precessional motion of trapped ions and electrons around the torus. In Ref.[8], the equation $\Lambda=0$ was solved for BAE/KBM/BAAE branches at the accumulation points of each of them, without the presence of EPs. The frequency of KBMs is found to be close to $\omega_{*pi}=(T_ic/e_iB)(\textbf{k} \times \textbf{b}) \cdot \nabla p_i /p_i\,$, although this can change due to the coupling with the BAAE mode. $\omega_{*pi}$ affects the polarization $|S|$ and reduces the otherwise large damping rate of BAAEs 6,8. In our model, modes with Alfvenic polarization have $\delta \Phi^{(0)}=\delta \Psi^{(0)}$ and $|S| \sim \beta^{1/2}$, while modes with significant acoustic component, such as BAAEs are identified by $|S|\gg \beta^{1/2}$.

When RHS of GFLDR is taken into account the potential $\delta \bar{W}_f (s,\alpha)$[10] determines the MHD stability of the mode, while the energetic particle term $\delta \bar{W}_k$ the EP drive [9]. Due to the low frequency of the modes, here we focus on the resonant interaction with the thermal particle precessional motion, which for deeply trapped ions can be described by the EP term [1] in the small FLR/FOW limit:

$\delta \bar{W}_{kt}=\frac{2\pi^2 e^2}{mc^2 |s|}q R_0 B_0 \int d\varepsilon \int d\mu \, (\overline{\omega}_d/k_\theta )^2 \tau_b \, QF_0 \frac{1}{\overline{\omega}_d-\omega} \, .$

Here, $\tau_b$ is the thermal ion bounce period and $QF_0=(\omega \partial_\varepsilon+(\textbf{k} \times \textbf{b})/\omega_c \cdot \nabla \,)F_0$.

In the case of BAAEs, solutions of the GFLDR show, consistently with the experiment, that the core plasma effects ($\omega_*$ and $\overline{\omega}_{Di}$ resonance) play a crucial role in the excitation of the mode, much more than the energetic particles. The low frequency Alfvenic mode, which we identify as KBM, is easier to excite and closely related to the ion diamagnetic frequency. Our work shows that GFLDR is a general and comprehensive tool illuminating the nature of the fluctuations, appropriate for understanding the effects of EPs and core plasma on the low frequency modes, as well as for explaining the experimental observations.

[1] Tsai S T and Chen L 1993 Phys. Fluids B 5 3284
[2] Gorelenkov N N, Berk H L, Fredrickson E and Sharapov S E 2007 Phys. Lett. A 370 70
[3] Heidbrink B 2020 in private communication
[4] H. S. Zhang, Y. Q. Liu, Z. Lin and W. L. Zhang 2016 PHYSICS OF PLASMAS 23 042510
[5] Zonca F, Chen L and Santoro R A 1996 Plasma Phys. Control. Fusion 38 2011
[6] Zonca F, Biancalani A, Chavdarovski I, Chen L, Di Troia C and Wang X 2010 Journal of Physics: Conference Series 260 012022
[7] Chavdarovski I and Zonca F 2009 Plasma Phys. Control. Fusion 51, 115001 (22pp)
[8] Chavdarovski I and Zonca F 2014, Plasma Phys 21, 052506
[9] Fulvio Zonca and Liu Chen 2014, Plasma Phys 21, 072121
[10] Ruirui Ma, Ilija Chavdarovski, Gaoxiang Ye and Xin Wang, Physics of Plasmas 21, 062120 (2014)

Affiliation National Fusion Research Institute Korea, Republic of

### Primary author

Ilija Chavdarovski (National Fusion Research Institute)

### Co-authors

Fulvio Zonca (Associazione EURATOM-ENEA sulla Fusione) Liu Chen (Zhejiang University)