Numerical Verification of Bounce-Harmonic Resonances in Neoclassical Toroidal Viscosity for Tokamaks

This Letter presents the first numerical verification for the bounce-harmonic (BH) resonance phenomena of the neoclassical transport in a tokamak perturbed by nonaxisymmetric magnetic fields. The BH resonances were predicted by analytic theories of neoclassical toroidal viscosity (NTV), as the parallel and perpendicular drift motions can be resonant and result in a great enhancement of the radial momentum transport. A new drift-kinetic $\delta f$ guiding-center particle code, POCA, clearly verified that the perpendicular drift motions can reduce the transport by phase-mixing, but in the BH resonances the motions can form closed orbits and particles radially drift out fast. The POCA calculations on resulting NTV torque are largely consistent with analytic calculations, and show that the BH resonances can easily dominate the NTV torque when a plasma rotates in the perturbed tokamak and therefore, is a critical physics for predicting the rotation and stability in the International Thermonuclear Experimental Reactor.