Extended Numerical Modeling of Impurity Neoclassical Transport in Tokamak Edge Plasmas

Understanding of impurity transport in tokamaks is an important issue in order to reduce the impurity contamination in fusion core plasmas. Recently, a new kinetic numerical scheme of impurity classical/neoclassical transport has been developed. This numerical scheme makes it possible to include classical self-diffusion (CL SD), classical inward pinch (CL IWP), and classical temperature screening effect (CL TSE) of impurity ions. However, impurity neoclassical transport has been modeled only in the case where background plasmas are in the Pfirsch-Schluter (PS) regime. The purpose of this study is to extend our previous model to wider range of collisionality regimes, i.e., not only the PS regime, but also the plateau regime. As in the previous study, a kinetic model with Binary Collision Monte-Carlo Model (BMC) has been adopted. We focus on the modeling of the neoclassical self-diffusion (NC SD) and the neoclassical inward pinch (NC IWP). In order to simulate the neoclassical transport with the BCM, velocity distribution of background plasma ions has been modeled as a deformed Maxwell distribution which includes plasma density gradient. Some test simulations have been done. As for NC SD of impurity ions, our scheme reproduces the dependence on the collisionality parameter in wide range of collisionality regime. As for NC IWP, in cases where test impurity ions and background ions are in the PS and plateau regimes, parameter dependences have been reproduced.