Modeling of ion losses in an NSTX/HHFW scenario with the RF Monte-Carlo code MCGO coupled to TORIC

The 5D and 6D Monte-Carlo code MCGO [1] has been recently upgraded with an RF kick operator based on full wave fields from AORSA, and was tested against CQL3D in a scenario of minority ICRH in C-Mod [2]. The operator utilizes local quasilinear diffusion coefficients that include all components of RF electric field and all-orders terms for finite Larmor radius effects. The RF operator acts on both perpendicular and parallel components of particle momentum. In the present study, the simulations are performed for the case of High Harmonic Fast wave (HHFW) heating in NSTX. The RF field data needed for the RF kick operator is obtained from TORIC runs. The MCGO runs are made in both zero-orbit-width and finite-orbit-width configurations. Of particular interest are the fast ion loss spectra to the 2D wall. A procedure for accumulation of particle loss distributions to the wall of test ions and neutrals is included as a function of poloidal and toroidal distance along the wall, particle energy, and the incident angle. With the applied RF power, the count of lost ions is increased only moderately, however, the heat load to the wall is increased by 60%.