D.1 Cold Plasma Model

In this section, we use the cold plasma model to calculate the wave properties. This approximation is valid only if the following conditions are satisfied:

In the first subsection, we discuss the modeling of RF k spectrum. In the second subsection, we use the cold plasma model to calculate wave properties. In the third and fourth subsections, we use further approximations to give analytical formulas for the properties of LHW and ECW. This leads to a simplified description of LHCD and ECCD and allows us to compare our models and results with other codes.

D.1.1 Wave Equation and Dispersion Tensor

We consider the two-fluids description of a non relativistic plasma in a constant magnetic field B0 = B0^z . In a continuous homogeneous linear medium, a Fourier component Eb(kb,ωb) verifies the wave equation [?]

Nb ×  (Nb ×  Eb)+ K ⋅Eb =  0
(D.2)

or equivalently

                2
(Nb ⋅Eb )Nb  - NbEb + K ⋅Eb =  0
(D.3)

where Nb = ckb∕ωb is the wave refractive index and K is the dielectric tensor. This equation can be expressed in tensorial form as

D ⋅Eb = 0
(D.4)

where the dispersion tensor is

             2
D =NbNb  - N bI+ K
(D.5)

By evaluating the susceptibilities in the cold plasma limit, the following dielectric tensor is obtained [?] in the frame (x,y,z):

    (                  )
       K ⊥   - iKxy  0
K = ( iKxy    K ⊥    0 )
        0      0    K ∥
(D.6)

whose elements in a plasma made of species s are given by

K = 1 - s    2
---ωps--
ω2 - ω2cs
K = 1 - sω2
--ps2
 ω
Kxy = - s   2
--ωpsωcs---
ω(ω2 - ω2cs) (D.7)

where

      ∘ ------
           2
ωps =   nsqs-
        ε0ms
(D.8)

is the plasma frequency and

ωcs = |qs|B0
       ms
(D.9)

the cyclotron frequency for the species s.

D.1.2 Dispersion Relation

In order to have a non-trivial solution to the wave equation, the dispersion relation must therefore be satisfied:

D  (kb,ωb) = |D| = 0
(D.10)

The axis ^x is chosen such that

Nb=Nb  ⊥x^+ Nb ∥z^
(D.11)

Taking the parallel component Nb as given, (D.10) lead to the following equation for Nb

AN 4b⊥ + BN 2b⊥ + C = 0
(D.12)

with

A = K (D.13)
B = (  2      )
 N b∥ - K⊥(        )
 K ⊥ + K ∥ + Kxy2 (D.14)
C = K[ (         )2      ]
   N 2b∥ - K⊥   - K2xy (D.15)

or equivalently

A = K (D.16)
B = (K   + K )
   ⊥     ∥Nb2 -(K  K   + K  K )
   R  L    ⊥   ∥ (D.17)
C = Kb(         )
 N 2b∥ - KR(         )
 N 2b∥ - KL (D.18)

with

KR = K + Kxy (D.19)
KL = K- Kxy (D.20)

We find the expression for Nb as a function of Nb and ω:

            √ ----------
N2  = --B-±---B2---4AC-
 b⊥          2A
(D.21)

We see immediately that the resonances occur for A = 0, that is,

K   = 0
  ⊥
(D.22)

and that cut-offs occur for C = 0, that is,

        2
KR  = Nb∥
(D.23)

K   = N 2
  L    b∥
(D.24)

or

K  =  0
  ∥
(D.25)

D.1.3 Polarization components

Once Nb has been evaluated, the components of the polarization of the electric field are simply given as the eigenvector of the wave equation (D.4)

D ⋅Eb = 0
(D.26)

D.1.4 Power flow

The relation between power flux and electric field was described using the vector Φb defined in (4.284)

Φb = ΦbP + ΦbT
(D.27)

with (4.287)

                        *
ΦbP = Re [Nb - (Nb ⋅eb)eb]
(D.28)

and (4.286)

        1 *  ∂KH
ΦbT = - 2eb ⋅∂N---⋅eb
                b
(D.29)

In the cold plasma model, the dielectric tensor K is independent of Nb and the contribution ΦbT , from the coherent motion of particles, vanishes.

ΦbT = 0
(D.30)

It is usually a very good approximation to neglect this kinetic power flux, even for quasi-electrostatic waves like Lower-Hybrid waves (See D.2.7).

D.1.5 Conclusion

The cold plasma model gives expressions (D.21), (D.26) and (D.28) for the calculation of N, eb and Φb respectively. These formulas can be generally used for the calculation of the diffusion coefficient in LHCD and ECCD.