The authors have written a code for the solution of the integro-differential equations describing coupling, propagation and absorption of high frequency waves in the ion cyclotron frequency range in tokamak plasmas, taking into account finite Larmor radius effects and parallel dispersion. The wave equations are discretized by a semi-spectral approach, using cubic Hermite finite elements in the radial direction and an expansion in Fourier modes in the poloidal direction. The latter permits analytic evaluation of the orbit integrals along magnetic field lines in the presence of a poloidal component of the static magnetic field. Thus, the code describes mode conversion to ion Bernstein waves, ion cyclotron damping at the fundamental and at the first harmonic, and electron transit time and Landau damping in full toroidal geometry. A few simulations of ion cyclotron heating of the ASDEX tokamak are presented. They show generally good agreement with the predictions of simpler models, namely plane layered models and ray tracing. Furthermore, the simulations can be used to identify a few interesting toroidal effects, particularly regarding the efficiency of mode conversion and the propagation of ion Bernstein waves.

In a medium with spatial dispersion the power Pabs dissipated per unit volume by an electromagnetic wave differs from ( J * E ) by the divergence of a vector which is interpreted as the correction askln to the Poynting flux due to the coherent particle motion in the field. We obtain expressions for PabE and askln in hot, inhomogeneous, weakly collisional plasmas by evaluating the time-averaged rate of change of the total kinetic energy density of charged particles, generalizing a method suggested by MCVEY et al. (1985, Phys. Rec. Lett. 55, 507).We prove that these expressions are consistent with the time averaged Poynting theorem for Maxwell equations. By specializing them to uniform and plane-stratified geometries, we derive some interesting properties of wave-plasma energy exchanges, and we check that the well known expressions for Pabs and

The code F ISIC solves the integro-differential wave equations for h.f. waves in the ion cyclotron frequency range in the fin ite Larmor radius approximation in tokamak geometry. The equations are discretised by a sernispectrai approach, using cubic Hermilc finite elements in the radial direction, and an expansion in Foufier modes in the poloidal direction. The latter allows the analytic evaluation of the orbit integrals along magnetic field lines in the presence of a poloidal component of the sLatic magnetic field. The code is thus able to describe mode conversion to ion Bernstein waves, ion cyclotron damping at the fundamental and at the first harmonic, and electron transit time and Landau damping, in full toroidal geometry. A detailed description of the model and the numerical method will be given elsewhere / 1/. First results from FISIC were presented in /2/. Since then a few improvements have been made, particularly in the description of collisional and electron damping.Here we present two simu lations of ICH in ASDEX (R = 167 cm, a = 40 cm,

The two dimensional electrostatic plasma particle in cell (PIC) code described an [l]