Description
For the modeling of wave heating and current drive in magnetic fusion plasmas, stationary full wave analysis based on the finite element method (FEM) is a powerful and flexible scheme in various configurations. An issue remaining in the FEM analysis, however, is to describe kinetic effects due to thermal motion of particles without wave number. Usually, kinetic effects in dielectric tensor are expressed with wave number vector in a uniform plasma. Several approaches based on cold-plasma wave number, differential operator, or Fourier analysis, have been developed for inhomogeneous plasmas, but their applicability is limited. We have developed a systematic scheme using integral form of dielectric tensor for kinetic full wave analysis in inhomogeneous plasmas. The integral form is derived by following the particle orbit in inhomogeneous plasmas. In Maxwellian plasmas, two kinds of kernel functions have been derived: the plasma dispersion kernel function(PDKF) for parallel motion and the plasma gyro kernel function (PGKF) for perpendicular motion. One-dimension analyses were applied to various phenomena in inhomogeneous plasmas; laser-plasma resonant interaction, magnetic beach heating in the presence of linear magnetic inhomogeneity, localized cyclotron heating near the extremum of magnetic field, and mode conversion to the Bernstein waves. In tokamak plasmas, two-dimensional inhomogeneity has to be taken into account. The two-dimensional FEM full wave code TASK/WF2D was extended to use the integral form of dielectric tensor and applied to the analysis of ion cyclotron heating in the presence of energetic ions and coupling with Bernstein waves. It was also applied to the lower hybrid current drive and the O-X-B mode conversion of the electron cyclotron waves. The extension to include relativistic effects and arbitrary velocity distribution function will be discussed.