Implementation and Benchmark of TORIC-FP4D for Ion Cyclotron Resonance Heating in Tokamak Devices

Ion Cyclotron Resonance Heating (ICRH) in a tokamak is evaluated by Fokker-Planck code, such as CQL3D [1], coupled with wave code like TORIC [2]. Typically, Fokker-Planck codes assume that changes in energy and magnetic momentum due to RF waves with bounce motion are negligible. However, this assumption may not hold in the presence of strong poloidal inhomogeneous sources. To address this, a new Fokker-Planck code, FP4D [3], has been developed, which does not use bounce averaging and thus considers the poloidal angle in real space. The RF heating term in TORIC-FP4D is coupled to the Kennel-Engelmann quasi-linear diffusion form [4], and the coefficients are calculated using TORIC. In this study, we performed the simulation for ITER minority heating and compared the results of TORIC-FP4D for m=0 mode with TORIC-CQL3D by mode decomposition for poloidal angle. We also compare the magnitudes of m=0 and m/=0 modes to discuss effects not seen in the bounce-average code. Additionally, we present preliminary results on neoclassical impurity transport by ICRH.