ORBIT modeling of fast particle redistribution induced by sawtooth instability

Initial tests on the National Spherical Torus Experiment Upgrade (NSTX-U) (Menard et al 2012 Nucl. Fusion 083015) device suggest that introducing energy selectivity for sawtooth induced fast ion redistribution is required to improve the agreement between experimental and simulated quantities such as neutron rate and Fast-Ion D-Alpha profiles. The aim of this work is to assess the requirements to properly describe the behaviour of fast ions during a sawtooth crash for predictive sawtooth simulations. As the first step, in this work, we use the particle-following ORBIT code to characterize the redistribution of fast particles. In order for a sawtooth crash to be simulated, a spatial and temporal displacement is implemented into the ORBIT code. The perturbation amplitude is determined by comparison with experimental measurement of the neutron rate drop. The characteristics of fast ions with different orbit types are investigated in phase and real space. Due to a sawtooth crash, fast ion energy and angular momentum are modified resulting in the redistribution in phase space and orbit type change. The redistribution of fast ions in real space shows that the sawtooth instability brings different effect on fast particles with different orbit types as observed in experiments. The initial interpretative TRANSP simulation using the so-called kick model based on the ORBIT modeling result shows an improvement of fast ion redistribution before and after a sawtooth crash but the neutron rate still has discrepancy compared to the experimental measurement.