Impurity behavior during sawtooth activity in tokamak plasmas

The transport of impurities by a sawtooth crash is simulated with the XTOR-2F code. Impurities are modeled as passive scalars, evolving in the compressible MHD flow inferred from the main MHD plasma. For a peaked impurity density profile, the non-linear kink flow of the sawtooth crash redistributes the profile efficiently and most of the particles in the peak inside the q = 1 surface are expelled. For an initially hollow impurity density profile, the crash leads to a significant penetration up to the magnetic axis. The results are compared with Kadomtsev's model. Despite essentially different mechanisms, the evolution of the particle content inside the q = 1 surface for Kadomtsev's model and for the non-linear case are virtually identical for the peaked profile, while the model slightly overestimates penetration for the hollow case.