Modeling of tungsten impurity transport and distribution in EAST based on multi-fluid and kinetic Monte Carlo simulations

Abstract

In recent years, the materials of plasma facing components, such as divertor target plates, domes, and outer walls of tokamaks, such as ASDEX Upgrade, WEST, JET, EAST, and ITER, have been changed from carbon to tungsten because of its lower erosion and tritium retention rates. Impurities are produced by interactions between the plasma and the first wall. This study provides an investigation to simulate the transport and distribution of tungsten impurities in the edge plasma on EAST. The 2D multi-fluid edge plasma transport code SOLPS-ITER and 2D kinetic Monte Carlo impurity transport code DIVIMP were used in the simulations. The multi-fluid model in SOLPS-ITER and the kinetic Monte Carlo model in DIVIMP were employed to treat tungsten impurity ions. The 2D density contour distributions in the computational region and the 1D density radial profiles at the inner and outer midplanes of tungsten impurity particles with ionization states (W0–W+74) and the total tungsten particles with all charge states were obtained. With the heating power 1.5 MW and the line-averaged plasma density 2 × 1019 m−3, the results from SOLPS-ITER and DIVIMP show that the maximum density of tungsten ion with single ionization state is about 1014 m−3 and the total density of tungsten impurities with all charge states is about 1015 m−3 at the core boundary. To the best of our knowledge, the simulation results from SOLPS-ITER and DIVIMP are compared for the first time to benchmark SOLPS-ITER with the multi-fluid mode and DIVIMP with the kinetic model for tungsten impurity transport. The density distributions of tungsten impurities with different ionization states from SOLPS-ITER and DIVIMP are highly similar, and good agreement can be found under similar conditions involved in the calculation. From the comparison benchmark between SOLPS-ITER and DIVIMP for tungsten impurity transport, it can be concluded that the impurity transport approximation used by DIVIMP is good.