Account of non-standard orbits in computations of neoclassical toroidal viscous torque in the resonant plateau regime of a tokamak

This article extends theoretical details based on a short paper originally submitted
to the 2022 EPS conference in plasma physics [1]. The quasilinear theory of resonant transport regimes in a tokamak is developed for the general case of orbits forming various classes separated in phase space by homoclinic orbits with infinite bounce time. Beyond standard orbits (banana and passing orbits) also all types of non-standard orbits (e.g. “potato” orbits) are taken into account. In case of a weak radial electric field, such orbits are usually present only near the magnetic axis. If the radial electric field cannot be treated as weak, there can be arbitrary many classes, located elsewhere. The present approach covers all such cases and is demonstrated on a specific example of a radial electric field profile. The resulting quasilinear kinetic equation is applicable to compute neoclassical toroidal viscous (NTV) torque in a tokamak with non-axisymmetric magnetic field perturbations. A fully non-local approach to NTV computation has been realized in the upgraded version of the code NEO-RT. Based on a generalization of magnetic flux surfaces to drift surfaces, the notion of a local thermodynamic equilibrium is extended for our purpose. We obtain an expression for the integral toroidal torque within a chosen flux surface and dicuss means to compute such integrals taking singularities in bounce and precession frequencies into account.