Assessment of validity of local neoclassical transport theory for studies of electric-field root-transitions in the W7-X stellarator

Abstract The neoclassical ambipolarity condition governing the radial electric field in stellarators can have several solutions, and sudden transitions (in radius) between these can then take place. The radial position and structure of such a transition cannot be determined from local transport theory, and instead a non-rigorous model based on a diffusion equation for the electric field is usually employed for this purpose (Turkin et al 2011 Phys. Plasmas 18 022505). We compare global (full plasma volume) drift-kinetic simulations of neoclassical transport in the Wendelstein 7-X stellarator with this model and find significant discrepancies. The position r 0 of the transition is not predicted correctly by the diffusion model, but the radial structure of the transition layer is in reasonable agreement if the diffusion coefficient is chosen appropriately. In particular, it should depend on the plasma temperature in the same way as the plateau-regime coefficient of neoclassical transport theory or the gyro-Bohm diffusion coefficient. In the small-gyroradius limit, the prediction of r 0 by the diffusion model simplifies to the so-called Maxwell construction (Shaing 1984 Phys. Fluids 27 1567–9; Shaing 1984 Phys. Fluids 27 1924–6). However, this property also emerges from a wide range of other mathematical models in the appropriate limit. The basic assumption underlying these models is that the diffusion, or generalisations thereof, is independent of the radial electric field, which is however unlikely to be the case in practice. Presumably this fact explains the discrepancy between the diffusion model and the drift-kinetic simulations. Finally, it is found that global simulations replicate the phenomenon of spontaneous root transitions driven by variations in the electron-to-ion temperature ratio, as predicted by local theory in the small-gyroradius limit.