Toroidal modeling of thermal particle drift kinetic effects and sub-sonic plasma flow on internal kink mode

The stability of the n = 1 ( n is the toroidal mode number) internal kink mode in a tokamak plasma is numerically investigated, utilizing the full toroidal, ideal magnetohydrodynamic (MHD) code MARS-F [Y. Q. Liu et al., Phys. Plasmas 7, 3681 (2000)] and the nonperturbative MHD-kinetic hybrid code MARS-K [Y. Liu et al., Phys. Plasmas 15, 112503 (2008)]. This study focuses on two physics effects: the sheared toroidal flow of plasma and the drift kinetic effects from thermal particles (ions and electrons) on the internal kink instability. Within the subsonic flow assumption, the flow and flow shear effects on the internal kink are generally weak, but the kinetic effects can be strong. A significant reduction in the mode growth rate is predicted by the nonperturbative MHD-kinetic hybrid computations, when the precessional drift resonance contributions from both particle species are included and when the mode is not too strongly unstable according to the fluid theory. The stabilization/destabilization of sheared toroidal flow depends on the radial location of the local flow shear. Strong destabilization (stabilization) occurs when a negative (positive) flow shear is located near the q = 1 rational surface.The stability of the n = 1 ( n is the toroidal mode number) internal kink mode in a tokamak plasma is numerically investigated, utilizing the full toroidal, ideal magnetohydrodynamic (MHD) code MARS-F [Y. Q. Liu et al., Phys. Plasmas 7, 3681 (2000)] and the nonperturbative MHD-kinetic hybrid code MARS-K [Y. Liu et al., Phys. Plasmas 15, 112503 (2008)]. This study focuses on two physics effects: the sheared toroidal flow of plasma and the drift kinetic effects from thermal particles (ions and electrons) on the internal kink instability. Within the subsonic flow assumption, the flow and flow shear effects on the internal kink are generally weak, but the kinetic effects can be strong. A significant reduction in the mode growth rate is predicted by the nonperturbative MHD-kinetic hybrid computations, when the precessional drift resonance contributions from both particle species are included and when the mode is not too strongly unstable according to the fluid theory. The stabilization/destabilization of s...