ELECTRIC TURBULENCE IN A PLASMA SUBJECT TO A STRONG MAGNETIC FIELD

G. Loeper * † A. Vasseur * †

November 27, 2006

Abstract
We consider in this paper a plasma subject to a strong deterministic magnetic field and we investigate the effect on this plasma of a stochastic electric field. We show that the limit behavior, which corresponds to the transfer of energy from the electric wave to the particles (Landau phenomena), is described by a Spherical Harmonics Expansion (SHE) model.

1 Introduction

This paper is concerned with the effect of a stochastic electric field on a plasma subject to a strong magnetic field. This is motivated by the study of the electric turbulence in a fusion machine as a Tokamak. Tokamaks are used to confine high energy plasmas in order to obtain the conditions needed for nuclear fusion reactions to take place. The plasma evolves in a toroidal reactor and is confined in the heart of the torus by the the mean of a strong magnetic field. A classical approximation is to suppose the ions to be at rest. Then only the electrons are moving. Another classical approximation argument in this type of study is the following: we are here interested only in interactions of particles over short distances of the order of the Larmor radius, moreover we suppose that at this scale the curvature of the magnetic field-lines can be neglected and that the plasma can be considered to be homogeneous along these field-lines. Thus we can restrict ourselves to a bidimensional problem. In this approach, the Vlasov equation describing the evolution of the repartition function f   of the electrons is:
m ( ∂ f ∂ t + v ⋅ ∇ x f ) + q ( B v ⊥ + ∇ V t u r b ( t , x ) ) ⋅ ∇ v f = 0 , (1)
m   stands for the electron's mass, q   its electric charge, f   the distribution function on ( t , x , v ) ∈ R + × R 2 × R 2   , with t   the time variable, x   the space variable and v   the velocity variable. B   is the (constant) norm of the transverse magnetic field, ∇ V t u r b   is the turbulent electric field, v ⊥   is the velocity vector after a rotation of π / 2   . We denote 1 / ε = q B m   the cyclotronic frequency, and we want to study the effect of ∇ V t u r b   in the limit ε   going to zero. In the deterministic case, the limits of related problems have been studied by several authors. In the case of the Vlasov-Poisson system (when the electric field is coupled with the density ∫ f d v   ), the limit has been studied, even in the 3D framework, by Frénod and Sonnendrucker [?] . In the 2D framework, using a slow time scale adapted to the problem, the convergence of the averaged motion, ∫ R 2 f ( t , x , v ) d v   to the 2D Euler system of equations has been performed simultaneously by Brenier [?] , Golse and Saint-Raymond [?]  and Frénod and Sonnendrucker [?] . A general result in 3D taking into account the two effects has been performed by Saint-Raymond in [?] .
In our case we neglect the Poisson non linear effect, concentrating on the stochastic behavior of the equation. Hamiltonian chaos method suggest that the modes of the turbulent electric field interacting with the electrons are those having a frequency of ω n = 2 π n ε   with n   an integer. This is roughly speaking the Landau resonance. Then the quasi-linear theory, (see Garbet [?] ) predicts a diffusive behavior with respect to the velocity variable. The diffusion coefficient obtained by this method being constant, it can not take into account the abnormal diffusion phenomena. In this paper we are interested in turbulent electric fields whose spectrum is spread around the Landau frequency and whose spatial fluctuations are of the scale of the Larmor radius (of order ε   ). We will show that the limit system is then governed by the following equation:
∂ t ρ − ∂ e ( a ( e ) ∂ e ρ ) = 0 (2)
where e = | v | 2 / 2   , ρ   is the average of f   over a sphere | v | 2 = 2 e   and the diffusion parameter a ( e )   is an explicit function of the correlation of V t u r b   , the turbulent electric potential, and of the energy, thus allowing abnormal diffusion. This diffusion parameter is undimensionally defined by ( 5 ). This equation is similar to the so-called Spherical Harmonics Expansion (SHE) model in high field limit modeling microelectronics semiconductor devices (see P.Degond [?]  or Ben Abdallah, Degond, Markowich and Schmeiser [?] ). It describes the Landau phenomena of transfer of energy from the electric wave to the particles. This work uses the techniques introduced by Poupaud and Vasseur [?]  to derive diffusive equation from transport in random media. This method works directly on the equation and, for this reason, is different from the method used in previous works (see Kesten and Papanicolaou [?] , [?]  and Fannjiang, Ryzhik and Papanicolaou [?] ). The paper is organized as follows: the precise result is stated in Section  2 . In Section  3 we show how we can compute explicitly the diffusion coefficients. Finally we give the proof of the theorem in Section  4 .

2 Results

In the remainder of the paper we fix n   and we denote
∇ V ε ( t , x ) = ε ∇ V t u r b ( 2 π n ε t , ε x ) (3)
for the stochastic potential. Equation ( 1 ) takes then the following undimensional form:
∂ f ε ∂ t + v ⋅ ∇ x f ε + ( v ⊥ ε + 1 ε ∇ V ε ( t 2 π n ε , x ε ) ) ⋅ ∇ v f ε = 0 . (4)
we denote E   the expectation value of any variable and make the following assumptions on the electrostatic potential:
( H 1 ) V ε ∈ L ∞ ( R + ; W 3 , ∞ ( R 2 ) ) and N ( ε ) : = E ( ( ∥ V ε ∥ L ∞ ( W 3 , ∞ ) ) 3 ) < ∞ ,
( H 2 ) E V ε ( t , x ) = 0 , for all t ∈ R + , x ∈ R 2 ,
( H 3 ) V ε ( t , x ) , V ε ( s , y ) are uncorrelated as soon as | t − s | ≥ 1 ,
( H 4 ) E ( V ε ( t , x ) V ε ( s , y ) ) = A ( t − s , x − y ) + g ε ( t , s , x , y ) ,
with:
∂ | α | ∂ x α A ∈ L ∞ ( R × R 2 ) , for α ∈ N 2 , | α | ≤ 3 ,
∥ | ∇ x , y 2 g ε | + | ∇ x , y , y 3 g ε | ∥ L ∞ ( ( R + ) 2 × R 4 ) ⟶ ε → 0 0 ,
where ∇ x , y 2 g ε   is the matrix ( ∂ x i ∂ y j g ε ) i , j   and ∇ x , y , y 3 g ε   is ( ∂ x i ∂ y j ∂ y k g ε ) i , j , k   .
Those assumptions are the same than in [?] . Hypothesis (H1) is an assumption on the regularity of V ε   for ε   fixed. Indeed the norm N ( ε )   can go to infinity when ε   goes to 0.
Hypothesis (H2) fixes the averaged potential at 0 which is not restrictive. In view of ( 3 ), Hypothesis (H3) determines 2 π n ε   as the decorrelated lapse of time for the turbulent electric field. Namely, it is the bigger lapse of time t − s   such that the electric fields at time t   and at time s   can be dependent on each other. Finally Hypothesis (H4), which is very classical, can be seen as an homogeneity property which takes place at the local scale ε   , since a quadratic quantity which depends on four variables ( t , x , s , y )   , at the limit, depends only on two variables ( t − s , x − y )   .
We denote R s ( v )   the rotation of angle s   with center 0 of v   . We consider the angular average of A   :
A ~ ( t , x ) = 1 2 π ∫ 0 2 π A ( t , R θ x ) d θ .   We then have the following result:
Theorem 2.1 Let V ε   be a stochastic potential satisfying assumptions ( H )   and independent of the initial data f ε 0 ∈ L 2 ( R 4 )   . Let a ( e )   be the function defined by:
a ( e ) = 1 2 π n 2 ∫ 0 + ∞ ( − ∂ t t 2 A ~ ) ( − s 2 π n , 2 2 e | sin s 2 | ) d s . (5)
This function is non negative. Assume that there is a constant C 0   such that ∥ f ε 0 ∥ L 2 ( R 4 ) ≤ C 0   and:
ε ( 1 + N ( ε ) 2 ) → 0 . (6)
Let ρ ε   be the gyro-average of f ε   defined by :
ρ ε ( t , x , e ) = 1 2 π ∫ 0 2 π f ε ( t , x , R θ v ) d θ , (7)
for every v   such that | v | 2 / 2 = e   . Then up to extraction of a subsequence, E f ε 0   converges weakly in L 2 ( R 4 )   to a function f 0 ∈ L 2 ( R 4 )   , E f ε   converges weakly in L 2   to a function f ∈ L ∞ ( R + , L 2 ( R 4 ) )   , E ρ ε   converges in C 0 ( [ 0 , T ] , L 2 ( R 2 ) − w )   for all T > 0   toward a function ρ ∈ L ∞ ( R + ; L 2 ( R 2 × R + ) ) ∩ C 0 ( R + ; L 2 ( R 2 × R + ) − W )   with ρ ( t = 0 , x , e ) = ∫ | v | 2 = 2 e f 0 d v   .
This function ρ   is solution to:
∂ t ρ − ∂ e ( a ( e ) ∂ e ρ ) = 0 t > 0 , x ∈ R 2 , e ∈ R + , (8)
in the distribution sense. Finally E f ε   converges weakly in L 2 ( R + × R 4 )   toward ρ ( t , x , | v | 2 / 2 )   .
Remark: depending on the regularity of the function a ( e )   the solution of the Cauchy problem for equation ( 8 ) may be unique. In this case, the whole sequence ρ ε   converges to ρ   the unique solution of ( 8 ).

3 Explicit computation of the diffusion parameter

In order to explicit the behavior of a ( e )   we must need the correlation function A   . We assume that it follows a “Richardson-like” law A ( t , x ) = f ( t ) | x | α .   Then we have
a ( e ) = − 1 2 π n 2 ∫ 0 + ∞ ∂ t t 2 A ~ ( − s 2 π n , 2 2 e | sin s 2 | ) d s
= − 2 3 2 α e α / 2 2 π n 2 ∫ 0 + ∞ f ′ ′ ( − s 2 π n ) | sin s 2 | α d s
= K e α / 2 .
A necessary condition for a function of the form ρ ( t , e ) = γ ( t ) ρ 0 ( e / t β )   to be an auto-similar solution of ( 2 ) is that β = 2 4 − α .   we then have an abnormal diffusion in e = t 2 4 − α .   For example if α = 4 / 3   we find a ( e ) = K e 2 / 3   and an abnormal diffusion in e = t 3 / 4   .

4 Proof of the result

We denote S ( R 4 )   the Schwartz space and S ′ ( R 4 )   its dual. We denote 〈 ⋅ ; ⋅ 〉   the duality brackets between those two spaces. We recall that L 2 ( R 4 ) ⊂ S ′ ( R 4 )   and by extension we will denote 〈 ⋅ ; ⋅ 〉   as well for the scalar product on L 2 ( R 4 )   . For every linear operator P   on S ( R 4 )   we will denote in the same way P   its extension on S ′ ( R 4 )   defined for every ψ ∈ S ( R 4 )   by: 〈 P ψ ; η 〉 = 〈 ψ ; P * η 〉 , η ∈ S ( R 4 ) .   Finally we will say that ψ n ∈ S ′ ( R 4 )   converges to ψ ∈ S ′ ( R 4 )   in S ′ ( R 4 )   if for every η ∈ S ( R 4 )   , 〈 ψ n ; η 〉   converges to 〈 ψ ; η 〉   . (This is the weak convergence for S ′ ( R 4 )   .) Let us rewrite equation ( 4 ) in the following way:
{ ∂ t f ε + C f ε + B f ε ε = − θ t ε f ε f ε | t = 0 = f ε 0 (9)
where C , B , θ t ε   are linear operators on S ( R 4 )   defined by:
C = v ⋅ ∇ x
B = v ⊥ ⋅ ∇ v
θ t ε = 1 ε ∇ V ε ( t 2 π n ε , x ε ) ⋅ ∇ v .
Notice that C   and B   are deterministic and non dependent on ε   nor on t   unlike θ t ε   .
We introduce the projection operator J   defined on S ( R 4 )   which averages the values of the function on the spheres | v | 2 / 2 = e   . Namely, for η ∈ S ( R 4 )   : J η ( x , v ) = 1 2 π ∫ 0 2 π η ( x , R θ v ) d θ .   We call J   the ”gyroaverage operator”. This operator is self adjoint for the L 2   scalar product. Applying the projection operator J   on ( 9 ) and taking its expectation value leads to: ∂ t E ( J f ε ) + E ( J C f ε ) + E ( J B f ε ) ε = − E ( J θ t ε f ε ) .   We first study some properties of the operators in order to pass to the limit in the left hand side of this equation. Then we investigate the limit of E ( J θ t ε f ε )   following the procedure of [?] . Finally we derive the SHE equation giving the explicit form of the diffusion coefficient a ( e )   .

4.1 Properties of the operators

We have the following properties on the operators C   , B   , θ t ε   and J   :
Lemma 4.1 Operators C   , B   and θ t ε   are skew adjoint for the L 2   scalar product. Operators C   , B   , and J   commute with the expectation operator E   . The operator J   is the restriction on S ( R 4 )   of the orthogonal projector of L 2 ( R 4 )   into K e r B   . In particular:
∥ J ∥ ℒ ( L 2 ( R 4 ) ) = 1 ,
J 2 = J ,
K e r B = I m J .
In addition: J B = J C J = 0 .  
Proof. –Operators C   , B   and θ t ε   can be rewritten as b . D   , where D   is a gradient operator and b   a regular function verifying D ⋅ b = 0   . For every functions η 1 , η 2 ∈ S ( R 4 )   we have:
〈 b ⋅ D η 1 ; η 2 〉 = − 〈 η 1 ; D ( b η 2 ) 〉
= − 〈 η 1 ; b ⋅ D η 2 〉 .
Hence they are skew adjoint operators for the L 2   scalar product.
–The operator J   is clearly the L 2   projection on L 2   functions which depends only on | v | 2 / 2   with respect to v   . In polar coordinates v = ( r sin θ , r cos θ )   , we have B = ∂ / ∂ θ   . So J   is the projection on K e r B   .
–Operator J   is the projector on K e r B   , hence B J = 0   . Since B   is skew adjoint and J   is self adjoint, we have ( B J ) * = − J B = 0   .
–Let us fix η ∈ S ( R 4 )   . We denote: J η ( x , v ) = ρ η ( x , | v | 2 2 )   . Hence J C J η ( x , v ) = ∇ x ⋅ ( ρ η ( x , | v | 2 2 ) 1 2 π ∫ 0 2 π R θ v d θ ) = 0 .   Finally J C J = 0   .
–Since C   , B   and J   are linear and deterministic, they commute with E   . □   .
From those properties we deduce the following proposition:
Proposition 4.2 For every ε   and every t ∈ R +   we have: ∥ f ε ( t ) ∥ L 2 ( R 4 ) = ∥ f ε 0 ∥ L 2 ( R 4 ) .   There exists a function f 0 ∈ L 2 ( R 4 )   and a function f ∈ L ∞ ( R + ; L 2 ( R 4 ) )   such that, up to a subsequence, E f ε 0   converges weakly in L 2 ( R 4 )   to f 0   , E f ε   converges weakly in L 2 ( [ 0 , T ] × R 4 )   to f   for every T > 0   . For every t > 0   , f   verifies J f ( t ) = f ( t )   . The function J E f ε   is solution to:
{ ∂ t J E f ε + E ( J θ t ε f ε ) = w ε J E f ε | t = 0 = J E f ε 0 , (10)
where w ε   converges to 0 in S ′   .
Proof. Since C   , B   and θ t ε   are skew adjoint operators, we have: ∂ t 〈 f ε ( t ) ; f ε ( t ) 〉 = 0 ,   which gives the first equality. By weak compactness there exists two functions f 0 ∈ L 2 ( R 4 )   and f ∈ L ∞ ( R + ; L 2 ( R 4 ) )   such that, up to a subsequence, E f ε 0   converges weakly in L 2 ( R 4 )   to f 0   , E f ε   converges weakly in L 2 ( [ 0 , T ] × R 4 )   to f   for every T > 0   . Notice that ε θ t ε   converges to 0 in S ′ ( R 4 )   . Multiplying Equation ( 9 ) by ε   , taking its expectation value, and letting ε   go to 0, we find: B f ( t ) = 0 o n R + ,   since B   and E   commute. Thanks to Lemma  4.1  f ( t ) ∈ I m J   , and since J 2 = J   , we have J f ( t ) = f ( t )   for almost every t > 0   . Since J B = 0   , applying the operator J   on equation ( 9 ) and taking its expectation value gives: ∂ t J E f ε + E ( J θ t ε f ε ) = w ε ,   with w ε = − E J C f ε   . This converges in S ′   to − J C f = − J C J f = 0   , thanks to Lemma  4.1 . □   Hence we are now concerned by the limit in S ′   of E ( J θ t ε f ε )   .

4.2 Computation of E ( J θ t ε f ε )  

Let us denote S t ε t ∈ R   the group on S ( R 4 )   generated by the operator C + B / ε   . Namely, for every h ∈ S   , S t ε h   is the unique solution on R   to:
{ ∂ t g + C g + B g ε = 0 g | t = 0 = h . (11)
The operator S t ε   can be explicitly given by: S t ε h ( x , v ) = h ( T ε ( t ) ( x , v ) ) ,   where T ε ( t ) ( x , v ) = ( x + ε v ⊥ − ε R − t / ε v ⊥ , R − t / ε v ) .   The function T ε ( t ) ( x , v )   gives the position at − t   of the particle being in x   with speed v   at time 0   and moving at constant speed | v |   on a circle of radius ε | v |   . In particular S t ε   is 2 π ε   periodic. Notice that the adjoint of S t ε   is S − t ε   .
Following the procedure of [?] , we use a 2 times iterated Duhamel formula. The first iteration gives:
f ε ( t ) = S 2 π n ε ε f ε ( t − 2 π n ε ) − ∫ 0 2 π n ε ( S σ ε θ t − σ ε S − σ ε ) S σ ε f ε ( t − σ ) d σ
and then we write the Duhamel formula for the f ε ( t − σ )   in the integral, and this yields:
f ε ( t ) = S 2 π n ε ε f ε ( t − 2 π n ε ) − ∫ 0 2 π n ε ( S σ ε θ t − σ ε S − σ ε ) S 4 π n ε ε f ε ( t − 4 π n ε ) d σ
+ ∫ 0 2 π n ε ∫ 0 4 π n ε − σ ( S σ ε θ t − σ ε S − σ ε ) ( S s + σ ε θ t − σ − s ε S − s − σ ε ) S s + σ ε f ε ( t − σ − s ) d s d σ .
We obtain:
E ( J θ t ε f ε ( t ) ) = J E ( θ t ε S 2 π n ε ε f ε ( t − 2 π n ε ) )
− ∫ 0 2 π n ε E ( J θ t ε ( S σ ε θ t − σ ε S − σ ε ) S 4 π n ε ε f ε ( t − 4 π n ε ) ) d σ + r t ε (12)
with
r t ε = ∫ 0 2 π n ε ∫ 0 4 π n ε − σ J E ( θ t ε ( S σ ε θ t − σ ε S − σ ε ) ( S s + σ ε θ t − σ − s ε S − s − σ ε ) S s + σ ε f ε ( t − σ − s ) ) d s d σ .
The function f ε 0   is independent of the operators θ t ε , t ∈ R   . In particular, in view of the assumption (H3), θ t ε   and f ε ( t − s )   are independent as soon as t ≥ s + 2 π n ε   and s ≥ 0   .
Combining this fact with (H2), Equation ( 12 ) becomes for t ≥ 4 π n ε  
E ( J θ t ε f ε ) = J E ( θ t ε ) E ( S 2 π n ε ε f ε ( t − 2 π n ε ) )
− ∫ 0 2 π n ε E ( J θ t ε ( S σ ε θ t − σ ε S − σ ε ) ) E ( S 4 π n ε ε f ε ( t − 4 π n ε ) ) d σ + r t ε ,
E ( J θ t ε f ε ) = − ∫ 0 2 π n ε E ( J θ t ε ( S σ ε θ t − σ ε S − σ ε ) ) E f ε ( t ) d σ + r t ε + e t ε , (13)
with e t ε = − ∫ 0 2 π n ε E ( J θ t ε ( S σ ε θ t − σ ε S − σ ε ) ) ( E S 4 π n ε ε f ε ( t − 4 π n ε ) − E f ε ( t ) ) d σ .
Since S t ε   is 2 π ε   periodic, S 4 π n ε ε f ε ( t − 4 π n ε ) = f ε ( t − 4 π n ε )   . We have: S s ε θ t − s ε S − s ε = 1 ε ( S s ε E ε ( t − s ) ) ⋅ D s ε   where we denote E ε ( t , x ) = ∇ V ε ( t 2 π n ε , x ε ) ,   and we define the differential operator D s ε   by D s ε = R − s / ε ∇ v + ε R − s / ε ∇ x ⊥ − ε ∇ x ⊥ .   Note that D s ε   is skew adjoint. Let us introduce the operator L t ε   on S ( R 4 )   (extended on S ′ ( R 4 ) )   defined for every η ∈ S ( R 4 )   by: L t ε η = − ∫ 0 2 π n ε E ( θ t ε ( S σ ε θ t − σ ε S − σ ε ) ) η d σ .   We can gather those results in the following way:
Lemma 4.3 We have the following equality: E ( J θ t ε f ε ) = J L t ε E f ε + r t ε + e t ε ,   where the operator L t ε   is defined for every η ∈ S ( R 4 )   by:
L t ε η ( x , v ) = − 1 ε ∫ 0 2 π n ε ∇ v ⋅ ( E ( S σ ε E ε ( t − σ ) ⊗ E ε ( t ) ) ⋅ D σ ε η ( x , v ) ) d σ , (14)
and the remainders are defined by:
e t ε = J L t ε ( E f ε ( t − 4 π n ε ) − E f ε ( t ) ) ,
r t ε = 1 ε ε ∫ 0 2 π n ε ∫ 0 4 π n ε − σ J E ( E ε ( t ) ⋅ ∇ v ( S σ ε E ε ( t − σ ) ⋅ D σ ε ( S s + σ ε E ε ( t − s − σ )
⋅ D s + σ ε ( S s + σ ε f ε ( t − σ − s ) ) ) ) ) d s d σ .
We can now show the following lemma:
Lemma 4.4 For every η ∈ S ( R 4 )   , the remainder r t ε   verifies: | 〈 r t ε ; η 〉 | ≤ C ( η ) ε N ( ε ) ,   and ( L t ε ) * η   converges in L 2 ( R 4 )   to: ∫ 0 2 π n R − s ∇ v ⋅ ( ∇ x x 2 A ( − s 2 π n , v ⊥ − R − s v ⊥ ) ∇ v η ) d s .  
Proof. We have:
( L t ε ) * η = − 1 ε ∫ 0 2 π n ε D σ ε ⋅ ( E ( S σ ε E ε ( t − σ ) ⊗ E ε ( t ) ) ⋅ ∇ v η ) d σ
= − ∫ 0 2 π n D ε σ ε ⋅ ( E ( S ε σ ε E ε ( t − ε σ ) ⊗ E ε ( t ) ) ⋅ ∇ v ) η d σ .
But thanks to the definition to T ε ( s )   , E ε   and Hypothesis (H3), the term E ( S ε σ ε E ε ( t − ε σ ) ⊗ E ε ( t ) ) = E ( ∇ V ε ( t − ε σ 2 π n ε , x / ε + v ⊥ − R − σ v ⊥ ) ⊗ ∇ V ε ( t 2 π n ε , x / ε ) )   converges strongly to ( − ∇ x x 2 A ) ( − σ / ( 2 π n ) , v ⊥ − R − σ v ⊥ )   in L ∞ ( ( R + ) 2 ; W 1 , ∞ ( R 2 ) )   . Hence thanks to the definition of D s ε   , ( L t ε ) * η   converges strongly to: ∫ 0 2 π n R − s ∇ v ⋅ ( ∇ x x 2 A ( − s 2 π n , v ⊥ − R − s v ⊥ ) ∇ v η ) d s   in L ∞ ( R + × R 4 )   . We recall that we have D ε s ε = R − s ∇ v + ε R − s ∇ x ⊥ − ε ∇ x ⊥   thus
∥ D ε s ′ ε Φ ∥ L 2 ( R + × R 4 ) ≤ C ∥ Φ ∥ W 1 , 2 ( R + × R 4 ) ,
∥ D ε s ′ ε D ε s Φ ∥ L 2 ( R + × R 4 ) ≤ C ∥ Φ ∥ W 2 , 2 ( R + × R 4 ) ,
hence
| 〈 r t ε ; η 〉 | ≤ C ( n ) ε ∥ Φ ∥ W 3 , 2 ∥ f ε ∥ L 2
sup s , s ′ { E ( | E ε ( T ( ε s ′ ) ) | ( | E ε ( T ( ε s ) ) | + | D ε s ′ ε [ E ε ( T ( ε s ) ) ] | )
( | E ε ( t , x ) | + | D ε s ε E ε ( t , x ) | + | D ε s ′ ε E ε ( t , x ) | + | D ε s ′ ε D ε s ε E ε ( t , x ) | ) ) } .
We then use the following bounds: (recall that E ε ( t , x ) = ∇ V ε ( t 2 π n ε , x ε ) .   )
| E ε | ≤ ∥ ∇ x V ε ∥ L ∞ ( R + × R 2 ) )
| D ε s ε E ε ( t , x ) | ≤ ∥ ∇ x x 2 V ε ∥ L ∞ ( R + × R 2 ) )
| D ε s ′ ε D ε s ε E ε ( t , x ) | ≤ ∥ ∇ x x x 3 V ε ∥ L ∞ ( R + × R 2 ) )
| D ε s ′ ε [ E ε ( T ( ε s ) ) ] | ≤ C ε | ∇ x E ε | ( T ( ε s ) ) ≤ C ∥ ∇ x x 2 V ε ∥ L ∞ ( R + × R 2 ) ) .
Then Hypothesis ( H 1 )   ensures that | 〈 r t ε ; η 〉 | ≤ C ε ∥ f 0 ε ∥ L 2 N ( ε ) ∥ Φ ∥ W 3 , 2 ,   which ends the proof of the lemma. □   We can now state the following proposition:
Proposition 4.5 Assume that ε ( N ( ε ) ) 2   converges to 0 when ε   goes to 0. Then the convergence (up to a subsequence) of J E f ε   to f   holds in C 0 ( R + ; L 2 ( R 4 ) − w )   , and f   is solution to:
∂ t f + J L t 0 J f = 0 , (15)
where the operator L t 0   is defined for every η ∈ S ( R 4 )   by: ( L t 0 ) * η = ∫ 0 2 π n R − s ∇ v ⋅ ( ∇ x x 2 A ( − s 2 π n , v ⊥ − R − s v ⊥ ) ∇ v η ) d s .  
Proof. Thanks to the previous lemma, for every test function η ∈ S ( R 4 )   : | 〈 r t ε ; η 〉 | ⟶ ε → 0 0 ,   and ( L t ε ) * η   converges strongly in L 2 ( R 4 )   to ( L t 0 ) * η   . But, thanks to Proposition  4.2 , f ε   converges weakly to f   in L 2 − w .   So 〈 J L t ε E f ε ; η 〉 = 〈 E f ε ; ( L t ε ) * J η 〉   converges to:
〈 f ; ( L t 0 ) * J η 〉 = 〈 J L t 0 f ; η 〉 .   The function f ε ( t ) − S 2 π n ε f ε ( t − 2 π n ε )   converges to 0 in L 2 − w   as well. So e t ε   converges to 0 in S ′ ( R 4 )   . Passing to the limit in equation ( 10 ) gives equation ( 15 ). This shows that ∂ t f ε   is uniformly bounded in time in a negative Sobolev space. Hence f ε   converges to f   in the space of continuous function in time with values in this Sobolev space. Finally since f ε   is bounded in L ∞ ( R + ; L 2 ( R 4 ) )   , the convergence holds in C 0 ( [ 0 , T ] ; L 2 ( R 4 ) − w )   for every T > 0   . □  

4.3 Convergence to the SHE model

Since J f = f   , we can introduce the gyroaverage function defined by: ρ ( t , x , e ) = f ( t , x , v ) ,   for every v   such that 2 e = | v | 2   . This subsection is devoted to the proof of the following lemma:
Lemma 4.6 The function ρ   lies in C 0 ( R + ; L 2 ( R 2 × R + ) − w ) ∩ L ∞ ( R + ; L 2 ( R 2 × R + ) )   .
It is solution to: { ∂ t ρ − ∂ e ( a ( e ) ∂ e ρ ) = 0 ρ | t = 0 = 1 2 π ∫ 0 2 π f 0 ( t , x , R θ v ) d θ   where the diffusion parameter is defined by: a ( e ) = ∫ 0 2 π n ∫ 0 2 π R θ v ⋅ ( − ∇ x x 2 A ) ( − s 2 π n , R θ v ⊥ − R − s + θ v ⊥ ) ⋅ R − s R θ v d s d θ ,   for every v   such that e = | v | 2 / 2   .
Proof. Let us first compute the operator J L t 0 J   . Let η 1 , η 1   be two test functions in S ( R 4 )   .
We have:
〈 η 1 ; J L t 0 J η 2 〉 = 〈 ( L t 0 ) * J η 1 ; J η 2 〉
= ∫ R 4 ∫ 0 2 π n ∇ v J η 1 ( − ∇ x x 2 A ) ( − s 2 π n , v ⊥ − R − s v ⊥ ) R − s ∇ v J η 2 d s d x d v .
Let us denote ρ η i   for i = 1 , 2   the functions defined by: ρ η i ( x , | v | 2 2 ) = J η i ( x , v ) .   Using polar coordinates and noticing that d v = d θ d e   we find:
〈 η 1 ; J L t 0 J η 2 〉 = ∫ R 2 ∫ 0 ∞ ∂ e ρ η 1 ( x , e ) ∂ e ρ η 2 ( x , e ) ∫ 0 2 π n ∫ 0 2 π R θ e ⃗ ⋅
( − ∇ x x 2 A ) ( − s 2 π n , − R − s + θ + π 2 e ⃗ + R θ + π 2 e ⃗ ) ⋅ R − s + θ e ⃗ d s d θ d e d x
= ∫ R 2 ∫ 0 ∞ ∂ e ρ η 1 ( x , e ) ∂ e ρ η 2 ( x , e ) a ( e ) d e d x ,
where e ⃗ = ( 2 e , 0 )   . Hence for every test function ρ η   , let us multiply it by Equation ( 15 ) and integrate with respect to x , v   . Since d e d θ = d v   we find: ∂ t ∫ R 2 ∫ 0 ∞ ρ ( t , x , e ) ρ η ( x , e ) d x d e = ∫ R 2 ∫ 0 ∞ ρ ( t , x , e ) ∂ e ( a ( e ) ∂ e ρ η ( x , e ) ) d e d x .   This, with Proposition  4.5 gives the desired result. □   Remark: we have a family of equations parametrized by x ∈ R 2   , and the solutions of two equations at two distinct x   do not interact.

4.4 Explicit computation of the diffusion coefficient

We derive in the following a suitable form to the diffusion coefficient a ( e )   . We will show, in particular, that a ( e )   is non negative. From ( H 4 )   the correlation function A ( t , x )   is even with respect to t   and x   . This with ( H 3 )   gives:
Lemma 4.7 The correlation function A   satisfies:
S u p p A ⊂ [ − 2 π n , 2 π n ] × R 2 ,
∇ x A ( 0 , 0 ) = 0 ,
∂ s A ( 0 , 0 ) = 0 .
This last subsection is devoted to the following proposition. Theorem  2.1 follows from this proposition, Proposition  4.2 and Proposition  4.5 .
Proposition 4.8 Let us denote A ~ ( t , x ) = 1 2 π ∫ 0 2 π A ( R θ x , t ) d θ .   Then a ( e )   is non negative and equal to: 1 2 π n 2 ∫ 0 ∞ ( − ∂ t t 2 A ~ ) ( − s 2 π n , 2 e 1 − cos s ) d s .  
Proof. Thanks to Lemma  4.6 and lemma  4.7 , we have a ( e ) = ∫ s = 0 ∞ ∫ 0 2 π R θ v ⋅ ( − ∇ x x 2 A ) ( − s 2 π n , R θ v ⊥ − R − s R θ v ⊥ ) ⋅ R − s R θ v d s d θ .   Since
− ∇ x x 2 A ( − s 2 π n , v ⊥ − R − s v ⊥ ) ⋅ R − s v
= 1 2 π n ∇ x ∂ s A ( − s 2 π n , v ⊥ − R − s v ⊥ ) + ∂ s ( ∇ x A ( − s 2 π n , v ⊥ − R − s v ⊥ ) ) ,
we find
a ( e ) = 1 2 π n ∫ s = 0 ∞ ∫ 0 2 π R θ v ⋅ ∇ x ∂ s A ( − s 2 π n , R θ v ⊥ − R θ − s v ⊥ ) d s d θ
− ∫ 0 2 π R θ v ⋅ ∇ x A ( 0 , 0 ) d θ
= 1 2 π n ∫ s = 0 ∞ ∫ 0 2 π R θ e ⃗ ⋅ ∇ x ∂ s A ( − s 2 π n , R θ + π / 2 e ⃗ − R − s + θ + π / 2 e ⃗ ) d s d θ
where e ⃗ = ( 2 e , 0 )   . Let us do the change of variables s ′ = θ − s   to get a ( e ) = 1 2 π n ∫ 0 2 π ∫ R 1 { s ≤ θ } R θ e ⃗ ⋅ ∇ x ∂ s A ( s − θ 2 π n , R θ + π / 2 e ⃗ − R s + π / 2 e ⃗ ) d θ d s .   Next we have
R θ e ⃗ ⋅ ∇ x ∂ s A ( s − θ 2 π n , R θ + π / 2 e ⃗ − R s + π / 2 e ⃗ ) = − 1 2 π n ∂ s s 2 A ( s − θ 2 π n , R θ + π / 2 e ⃗ − R s + π / 2 e ⃗ )
− ∂ θ { ∂ s A ( s − θ 2 π n , R θ + π / 2 e ⃗ − R s + π / 2 e ⃗ ) } .
Integrating by parts the second term of the RHS gives
1 2 π n ∫ 0 2 π ∫ R 1 { s ≤ θ } ∂ θ { ∂ s A ( s − θ 2 π n , R θ + π / 2 e ⃗ − R s + π / 2 e ⃗ ) } d s d θ
= 1 2 π n ∫ 0 2 π ∂ θ ∫ R 1 { s ≤ θ } ∂ s A ( s − θ 2 π n , R θ + π / 2 e ⃗ − R s + π / 2 e ⃗ ) d s d θ
− 1 2 π n ∫ 0 2 π ∫ R δ s = θ ∂ s A ( s − θ 2 π n , R θ + π / 2 e ⃗ − R s + π / 2 e ⃗ ) d s d θ
= 1 2 π n ∫ − ∞ 2 π ∂ s A ( s − 2 π 2 π n , R π / 2 e ⃗ − R s + π / 2 e ⃗ ) d s
− 1 2 π n ∫ − ∞ 0 ∂ s A ( s 2 π n , R π / 2 e ⃗ − R s + π / 2 e ⃗ ) d s
− 1 2 π n ∫ 0 2 π ∂ s A ( 0 , 0 ) d s
The first two lines cancel by doing the change of variables s ′ = s − 2 π   and the third line vanishes thanks to Lemma  4.7 , thus
a ( e ) = 1 ( 2 π n ) 2 ∫ 0 2 π ∫ R 1 { s ≤ θ } ( − ∂ s s 2 A ) ( s − θ 2 π n , R θ + π / 2 e ⃗ − R s + π / 2 e ⃗ ) d θ d s
Doing the change of variables s ′ = θ − s   gives
a ( e ) = 1 ( 2 π n ) 2 ∫ 0 2 π ∫ 0 ∞ ( − ∂ s s 2 A ) ( − s 2 π n , R θ + π / 2 ( ( I − R − s ) e ⃗ ) ) d θ d s
= 1 2 π n 2 ∫ 0 ∞ ( − ∂ s s 2 A ~ ) ( − s 2 π n , | ( I − R − s ) e ⃗ | ) d s .
Finally
| ( I − R s ) e ⃗ | = ( | 1 − cos s | 2 + sin 2 s ) 2 e
= 2 ( 1 − cos s ) 2 e
= 2 e 1 − cos s
= 2 2 e | sin ( s / 2 ) |
which ends the proof of the second assertion. □   Computation of the sign of the diffusion coefficient.
Here we check the non-negativity of the diffusion coefficient by expressing it in another form.
Thanks to lemma  4.6 and to hypothesis ( H 3 ) , ( H 4 )   we have
a ( e ) = 1 2 N ∫ s = 0 + ∞ ∫ − 2 π N 2 π N R θ v ⋅ ( − ∇ x x 2 A ) ( − s 2 π n , R θ v ⊥ − R θ − s v ⊥ ) ⋅ R θ − s v d s d θ
for all N   . Then doing the change of variable s : = θ − s   we find:
a ( e ) = 1 2 N ∫ s ∈ R ∫ − 2 π N 2 π N 1 { θ ≥ s } R θ v ⋅ ( − ∇ x x 2 A ) ( s − θ 2 π n , R θ v ⊥ − R s v ⊥ ) ⋅ R s v d s d θ
But we remind that thanks to hypothesis ( H 4 )  
− ∇ x x 2 A ( s − θ 2 π n , R θ v ⊥ − R s v ⊥ ) = lim ε → 0 E ( ∇ x V ε ( − s 2 π n , − R s v ⊥ ) ⊗ ∇ x V ε ( − θ 2 π n , − R θ v ⊥ ) )
Thus
a ( e ) = lim ε , N 1 2 N ⋅
∫ ∫ − 2 π N 2 π N 1 { θ ≥ s } E ( [ ∇ x V ε ( − s 2 π n , − R s v ⊥ ) ⋅ R s v ] [ ∇ x V ε ( − θ 2 π n , − R θ v ⊥ ) ⋅ R θ v ] ) d s d θ
Interverting s   and θ   we see that we can replace 1 { θ − s ≥ 0 }   by 1 { s − θ ≥ 0 }   and finally by adding both we obtain:
a ( e ) = lim ε , N 1 4 N E ( [ ∫ − 2 π N 2 π N ∇ x V ε ( − s 2 π n , − R s v ⊥ ) ⋅ R s v d s ] 2 )
which is a positive quantity. □   Acknowledgments: The idea of this paper was born as the authors attended a course on plasma's turbulence given by Xavier Garbet at the CEA center of of Cadarache. They are glad to thank him for having introduced them to the subject and for subsequent discussions.

1   Laboratoire J.A.Dieudonné, Université de Nice-Sophia-Antipolis, Parc Valrose, 06108 NICE Cedex 2.

Supported by

2   Work supported by the European Atomic Energy Community EURATOM in the Research Laboratory Agremented by CEA n o   01-24 UMR n o   6621 CNRS-Universite de Nice.