Modeling of multi-dimensional impurity transport in a realistic tokamak geometry

A 3D Monte Carlo transport code of heavy metal impurities is developed. The code includes most of important processes of heavy metal impurities, such as Larmor gyration, friction force, Coulomb collision, multi-step ionization and recombination process. The code outputs the 2D density profiles of tungsten impurity on realistic tokamak geometry. Calculations are made for given background plasma profiles for attached and detached plasma, which are the typical conditions in front of the divertor plate. In the attached plasma state, impurity tungsten particles are ionized to higher charge states near the divertor plate due to high background plasma electron temperature. On the other hand, in the detached plasma state, ions in higher charge states exist in the upstream of the divertor region mainly due to low background plasma electron temperature in front of the target plate. Although the code is still under development, it well describes the qualitative feature of impurity transport in the realistic tokamak geometry.