On neoclassical impurity transport in stellarator geometry

The impurity dynamics in stellarators has become an issue of moderate concern due to the inherent tendency of the impurities to accumulate in the core when the neoclassical ambipolar radial electric field points radially inwards (ion root regime). This accumulation can lead to collapse of the plasma due to radiative losses, and thus limit high performance plasma discharges in non-axisymmetric devices.\\ A quantitative description of the neoclassical impurity transport is complicated by the breakdown of the assumption of small $\mathbf{E}\times \mathbf{B}$ drift and trapping due to the electrostatic potential variation on a flux surface $\tilde{\Phi}$ compared to those due to the magnetic field gradient. The present work examines the impact of this potential variation on neoclassical impurity transport in the Large Helical Device (LHD) stellarator. It shows that the neoclassical impurity transport can be strongly affected by $\tilde{\Phi}$. The central numerical tool used is the $\delta f$ particle in cell (PIC) Monte Carlo code EUTERPE. The $\tilde{\Phi}$ used in the calculations is provided by the neoclassical code GSRAKE. The possibility of obtaining a more general $\tilde{\Phi}$ self-consistently with EUTERPE is also addressed and a preliminary calculation is presented.