Speaker
Description
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.
Acknowledgements:
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.
References:
[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 |
