How to calculate the neoclassical viscosity, diffusion, and current coefficients in general toroidal plasmas

A novel method to obtain the full neoclassical transport matrix for general toroidal plasmas by using the solution of the linearized drift kinetic equation with the pitch-angle-scattering collision operator is shown. In this method, the neoclassical coefficients for both poloidal and toroidal viscosities in toroidal helical systems can be obtained, and the neoclassical transport coefficients for the radial particle and heat fluxes and the bootstrap current with the nondiagonal coupling between unlike-species particles are derived from combining the viscosity-flow relations, the friction-flow relations, and the parallel momentum balance equations. Since the collisional momentum conservation is properly retained, the well-known intrinsic ambipolar condition of the neoclassical particle fluxes in symmetric systems is recovered. Thus, these resultant neoclassical diffusion and viscosity coefficients are applicable to evaluating accurately how the neoclassical transport in quasi-symmetric toroidal systems deviates from that in exactly symmetric systems

Moment-equation methods for calculating neoclassical transport coefficients in general toroidal plasmas

A detailed comparison is made between moment-equation methods presented by H. Sugama and S. Nishimura [Phys. Plasmas 9, 4637 (2002)] and by M. Taguchi [Phys. Fluids B 4, 3638 (1992)] for calculating neoclassical transport coefficients in general toroidal plasmas including nonsymmetric systems. It is shown that these methods can be derived from the drift kinetic equation with the same collision model used for correctly taking account of collisional momentum conservation. In both methods, the Laguerre polynomials of the energy variable are employed to expand the guiding-center distribution function and to obtain the moment equations, by which the radial neoclassical transport fluxes and the parallel flows are related to the thermodynamic forces. The methods are given here in the forms applicable for an arbitrary truncation number of the Laguerre-polynomial expansion so that their accuracies can be improved by increasing the truncation number. Differences between results from the two methods appear when the Laguerre-polynomial expansion is truncated up to a finite order because different weight functions are used in them to derive the moment equations. At each order of the truncation, the neoclassical transport coefficients obtained from the Sugama?Nishimura method show the Onsager symmetry and satisfy the ambipolar-diffusion condition intrinsically for symmetric systems. Also, numerical examples are given to show how the transport coefficients converge with the truncation number increased for the two methods