Self-consistent full-wave and Fokker-Planck calculations for ion cyclotron heating in non-Maxwellian plasmas

High-performance burning plasma devices such as ITER will contain significant concentrations of non-thermal plasma particles arising from fusion reactions, neutral beam injection, and wave-driven diffusion in velocity space. Initial studies in 1-D [1] and experimental results [2] show that non-thermal energetic ions can significantly alter wave propagation and absorption in the ion cyclotron range of frequencies. In addition, these ions can absorb power at high harmonics of the cyclotron frequency where conventional 2-D global-wave models are not valid. In this work, the all-orders, full-wave solver AORSA [3] is generalized to treat non-Maxwellian velocity distributions. Quasi-linear diffusion coefficients are derived directly from the global wave fields and used to calculate the energetic ion velocity distribution with the CQL3D Fokker-Planck code [4]. Alternately, the quasi-linear coefficients can be calculated numerically by integrating the Lorentz force equations along particle orbits. Self-consistency between the wave electric field and resonant ion distribution function is achieved by iterating between the full-wave and Fokker-Planck solutions.\newline [1] R. J. Dumont, C. K. Phillips and D. N. Smithe, Phys. Plasmas \textbf{12}, 042508 (2005).\newline [2] A. L. Rosenberg, J. E. Menard, J. R. Wilson, \textit{et al}., Phys. Plasmas \textbf{11}, 2441(2004).\newline [3] E. F. Jaeger, L. A. Berry, J. R. Myra, \textit{et al}., Phys. Rev. Lett. \textbf{90}, 195001-1 (2003).\newline [4] R. W. Harvey and M. G. McCoy, in \textit{Proceedings of the IAEA Technical Committee Meeting on Advances in Simulation and Modeling of Thermonuclear Plasmas} (IAEA, Montreal, 1992).