Since 18 of December 2019 uses Nucleus credentials. Visit our help pages for information on how to Register and Sign-in using Nucleus.
May 10 – 15, 2021
Virtual Event
Europe/Vienna timezone
The Conference will be held virtually from 10-15 May 2021

Extended Bounce-Kinetic Model for Trapped Particle Mode Turbulence

May 13, 2021, 8:30 AM
Virtual Event

Virtual Event

Regular Poster Magnetic Fusion Theory and Modelling P5 Posters 5


Taik Soo Hahm (Seoul National University)


Bounce-kinetic model based on the modern nonlinear bounce-kinetic equations[1] has been used for gKPSP[2] gyrokinetic simulations and produced useful and promising results[3]. However, magnetically trapped particles were treated as deeply trapped in that TEM and ITG simulation. This paper reports on an extension including the barely trapped particles. This will allow simulations addressing the precession reversed particles’ effect, reversed shear plasmas, and more precise neoclassical polarization shielding[4]. Modern bounce-kinetic equation advances the distribution function $F(\bar{Y_1}, \bar{Y_2}, \bar{\mu}, \bar{J})$ according to
$$\begin{equation} {{\partial}\over{\partial t}} F + {d \bar{Y_1}\over{dt}} {{\partial F}\over{\partial \bar{Y_1}}}+{d \bar{Y_2}\over{dt}} {{\partial F}\over{\partial \bar{Y_2}}}=0 \tag{1} \end{equation}$$ where $\bar{Y_1}$ and $\bar{Y_2}$ are bounce-averaged magnetic flux coordinates of gyrocenter, $\bar{\mu}$ and $\bar{J}$ are the first and the second adiabatic invariant respectively. With the total bounce-center Hamiltonian including the perturbation $〈H〉$, ${d \bar{Y_1}\over{dt}} = {{c}\over{q}} {{\partial 〈H〉}\over{\partial \bar{Y_2}}}$ and ${d \bar{Y_2}\over{dt}} = - {{c}\over{q}} {{\partial 〈H〉}\over{\partial \bar{Y_1}}}$ describe the motion of bounce-centers. While the expression of $\bar{J}$ in terms of particle’s energy and pitch angle is well-known in terms of elliptic functions[5], their inversion is necessary to express Maxwellian distribution in terms of the action-angle variables. This is straightforward for deeply trapped particles. In this work, we find analytic expressions for barely trapped particles in terms of Lambert function. The associated Poisson equation in terms of F is derived via pull-back transformation from the bounce-center coordinates to gyro-center coordinates [4]. The neoclassical polarization density which quantifies the Rosenbluth-Hinton residual zonal flow level[6] is also calculated. Initial simulation results using this scheme will be reported.

This work was supported by R&D Program of ITER Burning Plasma Research and Development of ITER Plasma Exploitation Plan (Code No. IN1904) through the National Fusion Research Institute of Korea (NFRI) funded by the Government funds.

[1] B.H. Fong and T.S. Hahm, Phys. Plasmas 6, 188 (1999).
[2] J.M. Kwon et al, Nucl. Fusion 52, 013004 (2012).
[3] J.M. Kwon, Lei Qi, S. Yi, and T.S. Hahm, Comp. Phys. Commun. 215, 81 (2017).
[4] Lu Wang and T.S. Hahm, Phys. Plasmas 16, 062309 (2009).
[5] B.B. Kadomtsev and O.P. Pogutse, Soviet Phys. JETP 24, 1172 (1967).
[6] M.N. Rosenbluth and F.L. Hinton, Phys. Rev. Lett. 80, 724 (1998).

Affiliation Seoul National University
Country or International Organization Korea, Republic of

Primary authors

Taik Soo Hahm (Seoul National University) Yong Jik Kim (Seoul National University) Jae-Min Kwon (National Fusion Research Institute) Lei Qi (National Fusion Research Institute of South Korea)

Presentation materials