Numerical modelling of sawtooth crash using two-fluid equations

The nonlinear growth of internal kink modes is studied numerically using two-fluid equations. As already reported earlier (Günter et al 2015 Plasma Phys. Control. Fusion 57 014017), short sawtooth crash times (<100 μ s) are found for typical ASDEX Upgrade parameters, in agreement with experimental observations. These fast sawtooth crashes are associated with large parallel electric field perturbations, giving rise to the generation of supra-thermal electrons. Slow sawtooth crashes (~ms) are obtained only for a sufficiently small ( ) and a large local electron diamagnetic drift frequency at the surface, where q 0 is the value of the safety factor q at the magnetic axis. Kink modes are shown to drive sheared plasma rotation, propagating from the surface towards the magnetic axis during the nonlinear phase. After the sawtooth crash, the driven plasma rotation is in the co- (counter-) current direction inside (outside) the surface, as observed in experiments.