Quasilinear theory is developed by using canonical variables for a relativistic plasma. It is self-consistent, including momentum, pitch angle, and spatial diffusions. By assuming the wave field as a superposition of known toroidal and poloidal Fourier modes, the quasilinear diffusion coefficients are written in a form which can be directly evaluated using the output of a spectral full-wave solver of Maxwell equations in toroidal plasmas. The formalism is special for tokamas and, therefore, simple and suitable for simulations of cyclotron heating, current drive, and radio-frequency wave-induced radial transport in ITER.
1. Introduction
Interaction of radio-frequency (RF) wave with plasma in magnetic confinement devices has been a very important discipline of plasma physics. To approach more realistic description of wave-plasma interaction in a time scale longer than the kinetic time scales bounce-average is needed. The long-time evolution of the kinetic distribution can be treated with Fokker-Planck equation. The behavior of the plasma and the most interesting macroscopic effects are obtained by balancing the diffusion term with a collision term.
For the relativistic particles the action and angle variables initiated by Kaufman [1] are introduced. The technique of the area-conserved transformation proposed by Lichtenberg and Lieberman [2] is employed. A new invariant which actually is an implicit Hamiltonian is formed by using bounce average and from which the bounce frequency and precession frequency can be calculated. Using new action and angle variables the relativistic quasi-linear equation is derived, including spatial diffusion. It is different from Brizard and Chan [3] where magnetic field does not have the toroidal component. An elegant form of the quasi-linear kinetic equation including a term describing RF wave-induced radial diffusion has been derived by Eriksson and Helander [4]. This paper is an extension of their work with relativistic effects. In most of the above derivations for the quasi-linear operator it is assumed that the wave field in the vicinity of the resonances can be represented in an Eikonal form. Here, the field is supposed as a superposition of given toroidal and poloidal Fourier modes, instead. This is the representation used to solve the Maxwell equation in tokamak plasmas with a spectral full-wave code, for example, TORIC [5]. For the circulating particles, under the conditions of small Larmor radius and first harmonic resonance, the diffusion coefficient is compatible with the numerical code developed by Cardinali et al. [6].
The rest of the paper is organized as follows. In next section the exact guiding center variables are derived with Hamiltonian transformation. The center variables are derived with Hamiltonian transformation. The bounce-averaged quasi-linear equation is carried out in Section 3. A brief summary is presented in the last section.
2. Exact Guiding Center Variables
In tokamak configuration, the relativistic Hamiltonian of a charged particle can be expressed as
H=[(PR-ecAR)2+(PZ-ecAZ)2+𝔄/R2]c2+m02c4+eΦ,
where 𝔄 denotes (Pϕ-(e/c)RAϕ)2, AR, AZ, and Aϕ are the vector potential components of the magnetic field, and Φ is the electrical potential and assumed to be a function of Ψ. Here, Ψ is the poloidal flux of the magnetic field. m0 and e are the rest mass and charge. PR,Pϕ, and PZ are the canonical momenta conjugate to the cylindrical coordinates R, ϕ and Z, respectively.
The magnetic field can be expressed as
B=∇ϕ×∇Ψ+I∇ϕ,
where I is related to the poloidal current. Then, in tokamaks, we have
AR=0,AZ=-IlnRR0,Aϕ=-ΨR.
We introduce a generating function [7] for changing to the guiding center variables:
F1=-m0Ω0R022exp(Xm0Ω0R0)(lnRR0-Xm0Ω0R0)2tgα-ZX,
where
X=m0Ω0R0lnRCR0,
and Ω0 is the gyrofrequency in the toroidal field, ρ is the Larmor radius, α is the gyrophase, and subscripts 0 and c refer to the values at the magnetic axis and the guiding center position, respectively. X and α are the new coordinates conjugate to the momenta
PX=Z+ρsinα+ρ24RCsin2α,Pα=12m0ΩCρ2.
The other two canonical variables Pϕ and ϕ do not change in the new coordinates. The old coordinates are connected with the new ones through four equations:
PR=m0Ωcρe(ρ/Rc)cosαsinα,PZ=-X=-m0Ω0R0lnRCR0,R=RCexp(-ρcosαRC),Z=PX-ρsinα-ρ24RCsin2α,
where PX is actually the Z coordinate of the guiding center, Zc. Such fact that a momentum is turned to be a coordinate often occurs during area-conserved canonical transformation [2].
The Jacobian in the area-conserved transformation is unity [2], that is,
dτ=dPαdPXdPϕdαdXdϕ.
Substituting (7) and (8) into (1), the exact Hamiltonian for the relativistic particles is
H={2m0ΩcPα𝔅+1R2[Pϕ+eΨ]2}c2+m02c4+eΦ,
where 𝔅 denotes [(Rc/R)2sin2α+cos2α]. It is suitable for particle simulation. The equations of motion and Vlasov’s equation can be derived from the Hamiltonian.
3. QuasiLinear Equation
If the Larmor radius is smaller than the scale length of the system, a small parameter δ=ρ/b may be introduced where b is the scale length.
For the gyrokinetics the Hamiltonian in (12) could be averaged and reduced as
H=(2m0ΩcPα+m02uϕ2)c2+m02c4+eΦ
to the first order of δ, where uϕ=(Pϕ+eΨ)/m0R is the toroidal velocity.
We form a new invariant [2]:
Π=12π∮PxdX.
The action invariant is the toroidal flux enclosed by a particle orbit which actually is an implicit Hamiltonian from which the bounce frequency and precession frequency can be calculated.
For trapped particles in a large aspect ratio configuration, that is, ε=r/R0≪1, we get
Πt=8qR0m0(εΩ0Pα/m0)0.5π[E(k1)-(1-k12)K(k1)].
The bounce frequency and the precession frequency are obtained from (15) as in [8, 9]:
ωbt=π(εΩ0Pα/m0)0.52γqR0K(k1),ωζt=2Ω0PαγΩpm0R02[E(k1)K(k1)-12]+4Ω0PαŝγΩpm0R02[E(k1)K(k1)-(1-k12)],
where Ωp is the gyrofrequency in poloidal magnetic field, k12=uϕ02/(4εΩ0Pα), and ŝ is the magnetic shear. E(k1) and K(k1) are complete elliptic function of the first and second kinds.
For the circulating particles,
Πc=Ω0m0r22+2qR0m0uϕ0σπE(k),ωbc=πuϕ0σ2γqR0K(k),ωζc=qωbcσ+uϕ022γΩprR0[E(k)K(k)-(1-k22)]+uϕ02ŝγΩpR02E(k)K(k),
where σ represents direction of the circulating particle and k2=k1-2. The bounce-averaged gyrofrequency of the trapped particles is
Ω=Ω0γ*[1-2ε(E(k1)K(k1)-12)],
while for the circulating particles,
Ω=Ω0γ*{1-2εk2(1+ε)[E(k)K(k)-(1-k22)]},
where Ω is a nonlocal gyro-frequency (also seen in [4]).
Lamalle [10] seems to prove that the frequencies for the local and nonlocal scenarios in the resonances are compatible. If the momentum, Pi, is constants of motion, the frequency is ωi=dqi/dt=∂H/∂Pi which is bounce-averaged frequency [4]. The coordinates are qi=ωit+const. where qi could represent η, ξ, or α [2]. The new momenta Π,Pα, and Pζ are conjugate to the coordinates η,α, and ζ in which Pζ=-(e/c)Ψ0 is actually a position variable [8].
In the extended phase space the Hamiltonian is written as follows [2]:
H̅(p̅,q̅)=H(p,q,t)-H,
where p̅n+1=-H,q̅n+1=t, H is the particle energy.
According to Liouville’s theorem, the distribution function, f, satisfies Vlasov’s equation:
∂f∂t+η̇∂f∂η+ς̇∂f∂ς+α̇∂f∂α+Π̇∂f∂Π+Ṗς∂f∂Pς+Ṗα∂f∂Pα+Ḣ∂f∂H=0,
where f can be divided in two parts, the averaged part and oscillatory part:
f=f̅+f̃.
The energy perturbation H1 from (1) is
H1=eΦ1-eυ⃗⊥·A⃗⊥-eυϕAϕ,
where considering (9) and (10), Φ1 and A⇀1 are
f̃=-(ωH1(∂f̅/∂H)+nH1(∂f̅/∂Pς)+ℭ)(ω-nως-mωη-lΩ),
where ℭ denotes mH1(∂f̅/∂Π)+lH1(∂f̅/∂Pα), and f̅ in (22) is a function of H, Pα, and Pζ only, independent of η,ς,α, and Π and satisfies the equation:
∂f̅∂t-i(n∂f̃•∂Pς+l∂f̃•∂Pα+ω∂f̃•∂H)H1=0
from (21), where we have used Hamiltonian equations, Ṗi=-∂H̅/∂qi, Pi and qi are canonical variables, and f̃• is the conjugate to f̃. If H1 is a slow function of Pς,Pα, and H, we get the quasi-linear equation:
∂f̅∂τ+L̂DL̂f̅=C(f̅),
where τ=tν, L̂=∂/∂H+(n/ω)(∂/∂Pς)+(l/ω)(∂/∂Pα) is differential operator of actions. Here, Pς represents spatial diffusion while H and Pα represent velocity space diffusion, C is a collision operator normalized to ν which can take the form given by [11], ν is collision frequency, and D is diffusion coefficient:
D=πH12δ(ω-nως-mωη-lΩ},
which is similar to the one obtained in [4]. We have used the well-known Plemelj formula called as Landau’s rule in [4]:
1ω-mωη-nωξ-lΩ=P1ω-mωη-nωξ-lΩ-iπδ(ω-mωη-nωξ-lΩ).
For one harmonic from (24) and (29) we get
H1l=[(eΦk-ecυςAφk)Jl(k⊥ρ)-ecυ⊥AZkJl-1+Jl+12-ecυ⊥ARkJl-1-Jl+12i],ilexp(ik⃗·r⃗c-ilα0}exp(inς+imη+lΩ-iωt),
where A⇀ is the vector potential of incident electromagnetic wave, r⇀c is the guiding center position, and α0 is the angle between k⇀ and -R directions in the cylindrical coordinates. Equations (27), (28), and (30) are the main results in this paper.
For the relativistic circulating particles the resonant term in (28) is written as
δ[ωγ(γ-lγΩω-N//qR0γωbc)]=cqR0ωωbδ(Nres-N//),
where γ is the relativistic factor, and Nres=(2K(k)/π)(γ-γΩ/ω)/(βthP//) which is consistent with the code developed in [6] if k approaches zero. However, there are differences for trapped or barely circulating particles.
From (28) and (31) we get
D=2πK(k)P⊥2γP//D⊥δ(Nres-N//),
where P⊥ and P// are normalized with respect to the thermal momentum Pt=(m0T)0.5. The steady state is described by a relativistic Maxwellian with temperature Tt. From the definition of the δ-function, when ΔN//→0,
δ(Nres-N//)=π-0.51ΔN//exp[-(Nres-N//)2ΔN//2],D⊥=πωe2A⊥2βthνT2,L̂=γ∂P∂P+nTtω∂∂Ψ•+lΩ0ω(1+ε)∂P2μ∂μ,
where the chosen variables are related to the momentum P and μ=cosθ, with θ being the pitch angle.
4. Summary
The action and angle variables initiated by Kaufman are used [1]. The area-conserved transformation is employed [2]. The bounce-averaged quasi-linear Fokker-Planck equation for the relativistic particles is rigorously derived, including momentum, pitch angle, and spatial diffusion, that is, the three-dimensional diffusion. Equations (27), (28), and (30) are the main results of the work. The diffusion coefficient expressed in (28) is similar to the one obtained by Eriksson and Helander [4]. Equation (34) is only for circulating particles. Most of the past derivations of the quasi-linear operator assume that the wave field in the vicinity of the resonances can be represented in a Eikonal form. Here, the field is supposed, instead, to be given as a superposition of toroidal and poloidal Fourier modes. This is the representation used to solve Maxwell equation in tokamak plasmas with spectral full-wave code, for example, TORIC [5]. The formalism is special for tokamas; so it is simple and suitable for simulation of cyclotron heating, current drive, and radio-frequency wave-induced radial transport in ITER.
Acknowledgments
The authors would like to give special thanks to Dr. A. Cardinali for his report and his code left here during his visit to SWIP. Z. T. Wang would like to thank Dr. R. D. Hazeltine and Dr. P. Morrision for checking the canonical transformation procedure and many helpful discussions when he worked at the Institute for Fusion Studies University of Texas at Austin. This research is supported by the National Natural Science Foundations of China under Grant nos. 10475043, 10535020, 10375019, and 10135020.
KaufmanA. N.Quasilinear diffusion of an axisymmetric toroidal plasma19721561063106910.1063/1.16940312-s2.0-36849104388LichtenbergJ.LiebermanM. A.198338New York, NY, USASpringerBrizardA. J.ChanA. A.Relativistic quasilinear diffusion in axisymmetric magnetic geometry for arbitrary-frequency electromagnetic fluctuations2004119422042292-s2.0-003803031710.1063/1.1773554ErikssonL.-G.HelanderP.Monte Carlo operators for orbit-averaged Fokker-Planck equations1994123083142-s2.0-0001228979BrambilaM.1999411CardinaliA.MoriniL.ZoncaF.ConnerJ.SauterO.SindoniE.871Proceedings of the Joint Varenna-Lausanne International Workshop on Theory of Fusion Plasmas2006Varenna, ItalyAmerican Institute of Physics292WangZ.Plasma transport at magnetic axis in toroidal confinement systems199941A679A68610.1088/0741-3335/41/3A/061HazeltineR. D.MahajanS. M.HitchcockD. A.Quasi-linear diffusion and radial transport in tokamaks1981246116411792-s2.0-36749117333WangZ. T.LongY. X.DongJ. Q.WangL.ZoncaF.Fishbone instability excited by barely trapped electrons2006231158160LamalleP. U.A qualitative comparison of theoretical models of radiofrequency wave propagation and absorption in tokamak plasmas19984044654792-s2.0-36149033373BraamsB. J.KarneyC. F. F.Conductivity of a relativistic plasma198917135513682-s2.0-5844308118