Modeling of beam ions loss and slowing down with Coulomb collisions in EAST

This paper uses the implicit Monte-Carlo full-orbit-following parallel program ISSDE to calculate the prompt loss and slowing down process of neutral beam injection (NBI)-generated fast ions due to Coulomb collisions in the equilibrium configuration of Experimental Advanced Superconducting Tokamak (EAST). This program is based on the weak equivalence of the Fokker-Planck equation under Rosenbluth MacDonald Judd (RMJ) potential and Stratonovich stochastic differential equation (SDE). The prompt loss with the LCFS boundary and the first wall (FW) boundary of the two co-current neutral injection beams are studied. Simulation results indicate that the loss behavior of fast ions using the FW boundary is very different from that of the LCFS boundary, especially for fast ions with a large gyration radius. According to our calculations, about 5.11% of fast ions generated by perpendicular injection drift out of the LCFS and then return inside the LCFS to be captured by the magnetic field. The prompt loss ratio of fast ions and the ratio of orbital types depend on the initial distribution of fast ions in the P_ζ-$\varLambda$ space. Under the effect of Coulomb collisions, the pitch-angle scattering and stochastic diffusion happens, which will cause more fast ion loss. For short time scales, among the particles lost due to collisions, the fraction of banana ions reaches 92.31% in the perpendicular beam and 58.65% in the tangential beam when the fraction of banana ions in the tangential beam is 3.4% of the total ions, which means that the effect of Coulomb collisions on banana fast ions is more significant. For long time scales, the additional fast ion loss caused by Coulomb collisions of tangential and perpendicular beams accounted for 16.21% and 25.05% of the total particles, respectively. We have also investigated the slowing down process of NBI fast ions.