A kinetic model for neoclassical transport in ELMO Bumpy Torus including the effects of a self-consistently calculated enhanced high energy tail on the ion distribution function

A drift kinetic equation is bounce-averaged over the rapid motion of particles parallel to the magnetic field, leaving an equation involving only the relatively slow drift motion of particles across the magnetic field. The drift velocity for particles in EBT is derived, and a Fokker-Planck collision operator assuming isotropic velocity distributions for both electrons and ions is used. Fourier expansion in the poloidal angle and a point model approximation for the radial dimension are used to further simplify the problem. The resulting three coupled nonlinear integro-differential equations are solved numerically, showing the formation of an enhanced high energy tail on the ion distribution function. Neoclassical loss rates are calculated with and without including the non-Maxwellian nature of the ion distribution function. The results show that the presence of the high energy tail vastly influences neoclassical losses in EBT.