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 [1]. 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 r0 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 r0 by the diffusion model simplifies to the so-called Maxwell construction [2, 3]. 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.
Read full abstract