The effect of a saturated kink on the dynamics of tungsten impurities in the plasma core

Abstract

We studied the effect of saturated kink-like instabilities on the redistribution of W impurities. A numerical code that follows the exact trajectories of the impurity ions in the total fields and includes collisions was employed. The code is written in CUDA C and runs in graphical processing units, allowing simulations with a large number of particles with modest resources. The total electric and magnetic fields are calculated as the sum of a simple analytical equilibrium (large aspect ratio, circular cross section) plus the fields of the kink mode. Our simulations show that when the W ions have a thermal velocity distribution, with zero mean velocity, the kink has no effect on the W density. However, when we include the plasma rotation and the drag on the impurities, the kink can affect the W density. When the average passing frequency of the W ions is similar to the frequency of the kink mode, the expulsion of W ions from the plasma core is maximized and the W density profile shows a hollow structure. This could have implications for the mitigation of W accumulation.