Kinetic calculation of neoclassical transport including self-consistent electron and impurity dynamics

Numerical studies of neoclassical transport, beginning with the fundamental drift-kinetic equation (DKE), have been extended to include the self-consistent coupling of electrons and multiple ion species. The code, NEO, provides a first-principles based calculation of the neoclassical transport coefficients directly from solution of the distribution function by solving a hierarchy of equations derived by expanding the DKE in powers of rho*(i), the ratio of the ion gyroradius to system size. This includes the calculation of the first-order electrostatic potential via the Poisson equation, although this potential has exactly no effect on the steady-state transport. Systematic calculations of the second-order particle and energy fluxes and first-order plasma flows and bootstrap current and comparisons with existing theories are given for multi-species plasmas. The ambipolar relation Sigma(a)z(a)Gamma(a) = 0, which can only be maintained with complete cross-species collisional coupling, is confirmed, and finite mass-ratio corrections due to the collisional coupling are identified. The effects of plasma shaping are also explored, including a discussion of how analytic formulae obtained for circular plasmas (i.e. Chang-Hinton) should be applied to shaped cases. Finite-orbit-width effects are studied via solution of the higher-order DKEs and the implications of non-local transport on the validity of the delta f formulation are discussed.