Abstract
Radial fluxes of parallel momentum due to and magnetic drifts are shown to be correlated in tokamak plasmas. This correlation comes from the onset of poloidal convective cells generated by turbulence. The entire process requires a symmetry breaking mechanism, e.g. a mean shear flow. An analytical calculation shows that anti-correlation between the poloidal and parallel components of the turbulent Reynolds stress results in anti-correlation of the fluxes of parallel momentum generated by and curvature drifts.
Export citation and abstract BibTeX RIS
Original content from this work may be used under the terms of the Creative Commons Attribution 3.0 licence. Any further distribution of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOI.
1. Introduction
Calculating particle, momentum and heat fluxes is of utmost importance in magnetized fusion plasmas. Indeed the confinement time is mainly determined by fluxes across the magnetic surfaces. These fluxes can be computed as velocity momenta of the distribution function for each species. At the microscopic level, the motion of a charged particle in a strong guide field can be separated into a cyclotron and gyrocenter motion. The velocity of particle gyrocenters perpendicular to the magnetic field can itself be split into two components. The first one is the E × B drift velocity, while the second one is the magnetic drift, due to field curvature and gradient of the magnetic field modulus. Hence gyrokinetic fluxes, defined here as fluxes of gyrocenters, contain contributions from the and magnetic drifts. Exact conservation laws can be derived in this framework [1]. Another approach, may be more intuitive, has been used in the past. It consists in separating the fluxes in collisional and turbulent parts. In toroidal fusion devices, collisional processes are amplified by trajectory effects. The resulting fluxes are called 'neoclassical'. Neoclassical and turbulent fluxes are often calculated separately and added up. Without turbulence, neoclassical fluxes match fluxes due to curvature drift, because the contribution of the E × B drift velocity is small [2, 3]. Cases exist however where the electric drift matters, like the neoclassical flux of toroidal momentum [4, 5], or impurity transport when poloidal asymmetries are large (e.g. due to the centrifugal force) [6]. Conversely, turbulent fluxes are usually computed by keeping only the fluctuating E × B velocity [7]. Hence a common methodology consists in using a code to calculate neoclassical fluxes (i.e. without turbulence), and a gyrokinetic code to compute turbulent fluxes, usually without direct contributions of the magnetic drift to fluxes. Both outputs are then added. This procedure can be justified under conditions of scale separation [8].
Recently several gyrokinetic codes have been upgraded to compute both collisional and turbulent transport, thanks to the implementation of accurate collisional operators. Results have been published for heat [9, 10], momentum [11–16] and particle [17] transport channels. These recent developments have raised the question of defining neoclassical and turbulent fluxes in a proper and practical way, and thus the issue of the interplay between collisional and turbulent processes. One way consists in defining the neoclassical flux as due to magnetic drift only, while turbulent fluxes are due to drift velocity fluctuations. This procedure finds some justification in the fact that it coincides with the conventional definition of a neoclassical flux in absence of turbulence, under the conditions given above, while the turbulent flux agrees with previous results that discarded the magnetic drift contribution to the flux. However this approach leads to inconsistencies, as shown in the references aforementioned. One reason is turbulent scattering in the velocity space, which plays a role similar to Coulomb collisions [15]. Another reason is the development of large scale poloidal asymmetries of the mean drift velocity due to turbulence. These asymmetries contribute to the neoclassical flux [9], whereas they are small in conventional neoclassical theory. These effects are sometimes called 'synergistic'. To avoid any misunderstanding and lengthy debates, fluxes are split in the present paper into contributions due to and curvature drifts, without using the words 'turbulent' and 'neoclassical'. We note that another possible definition would consist in separating fluxes due to large and small scales. Though the latter definition seems now to prevail, previous numerical works used the former one. Given these preliminaries, the question addressed in this paper is the following. Is the flux due to magnetic drift significant in a collisionless turbulent state? Since this is still a broad subject, it is limited here to toroidal momentum transport. This topic has triggered much attention due to the observation of intrinsic toroidal rotation in tokamaks, i.e. rotation without external torque (see [18] for an overview).
In several previous works, only the drift contribution to the parallel (or toroidal) momentum flux was calculated [14, 19–21]. However it was stressed in [22] that the contribution of the magnetic drift to the flux of momentum could be significant in non linear regime. Moreover two research teams recently computed both magnetic and drift contributions to the momentum flux. Both found that these two contributions are anti-correlated [13, 16], leading in some cases to an almost vanishing total flux. A mechanism is proposed here to explain the correlation between the contributions of magnetic and drifts to momentum transport, based on turbulent generation of poloidal convective cells. Poloidal convective cells are zonal structures with zero toroidal wavenumber, and finite poloidal wave numbers. We will make the difference with zonal flows by restricting the definition of zonal flows to structures with zero toroidal and poloidal wave numbers. The chain of causes starts with the well established correlation between radial fluxes of parallel and poloidal momentum due to mean shear flow or any other mechanism responsible for an up–down asymmetry of turbulence intensity. Turbulence is ballooned in a tokamak, i.e. is more intense on the low field side. The resulting asymmetry of the turbulent stress tensor generates poloidal convective cells, i.e. flows that are zonal in the toroidal rotation, but not in the poloidal direction. It appears that these flows are weakly damped at low frequency, and drive poloidal asymmetries of the distribution function, which contribute to a non zero flux of parallel momentum due to magnetic drift. Hence it appears that the magnetic drift component of the flux of parallel momentum is tied to the radial flux of poloidal momentum. Anti-correlation of the poloidal and parallel components of the turbulent Reynolds stress leads to anti-correlated radial fluxes of parallel momentum due to magnetic and drifts. This scheme is illustrated in figure 1. A schematic view of the circulation pattern in a tokamak poloidal plane is also shown in figure 2. The interplay between low wavenumber convective cells and turbulent transport has already been mentioned in the context of internal transport barrier formation [23].
Download figure:
Standard image High-resolution imageA simplified calculation is presented here, where the fluid parallel velocity and its gradient are small. In other words only the residual stress tensor is calculated. It appears that the sign of the correlation between the poloidal and parallel components of the Reynolds stress translates in a correlation of the curvature and drift contributions to the radial flux of parallel momentum. This general result comes from a simple relationship between the momentum flux due to curvature and the turbulent Reynolds stress (summarized in figure 2). Although poloidal convective cells do not appear explicitly in this relationship, they are instrumental to this mechanism. The overall process is illustrated by a quasi-linear calculation of stress tensors and fluxes based on linear ion temperature gradient (ITG) driven modes in the hydrodynamic limit. This calculation predicts similar amplitudes for the and magnetic drift components of the parallel momentum flux, but a sign of correlation that depends on plasma parameters. Typically for modes drifting in the ion diamagnetic direction, positive correlation is found for weak drive, and anti-correlation for strong drive. Since the hydrodynamic limit is valid only in the strong drive limit, this result can be considered as encouraging in view of previous numerical calculations, which found anti-correlation [13, 16].
The paper is organized as follows. General expressions of radial fluxes of parallel momentum are derived in section 2. The contribution from the magnetic drift is detailed in section 3. The and curvature driven fluxes of parallel momentum are compared and discussed in section 4. A conclusion follows. Most technical demonstrations can be found in the appendices, and may be skipped on first reading.
2. Fluxes of parallel momentum
Momentum flux is the same as the Reynolds stress tensor up to a mass density Nm, where N is the unperturbed density and m the ion mass. Both names will be used indistinctly throughout the paper. We calculate separately the contributions from the E × B drift and magnetic drift velocities to the flux of parallel momentum, i.e.
and
where F is the distribution function, is velocity that is aligned with the unperturbed magnetic field, is the E × B drift velocity and is the magnetic drift velocity due to curvature and gradient of the unperturbed magnetic field. The analysis is restricted to electrostatic turbulence, i.e.
where ϕ is the electric potential, and is the (unperturbed) magnetic field. Also we use a simplified expression of the magnetic drift velocity, valid at low values of the plasma β and normalized gyroradius ,
where m is the mass, e the charge, and is the modulus of the perpendicular velocity.
A simplified geometry of circular concentric magnetic surfaces is considered here. The spatial coordinates are , where r is the minor radius, θ and φ the poloidal and toroidal angles (see figure 2) and the magnetic field is
where the major radius, the safety factor, and B0 the magnetic field on the magnetic axis. The inverse aspect ratio is a small parameter throughout the paper.
2.1. Turbulence parametrization
2.1.1. Electric potential fluctuations
Since an evaluation of the turbulent Reynolds stress tensor is needed, it is necessary to postulate a spatial structure of the fluctuations of the electric potential. We follow a rationale that is close in spirit to the ballooning representation [24–29]. In a magnetized plasma with strong guide field, turbulent structures tend to align with the magnetic field. Given the structure of the magnetic field equation (5), it is convenient to use the variable instead of the toroidal angle. The periodicity in the toroidal direction allows using a a Fourier series
where each Fourier component can be written as a Fourier integral
valid locally near a reference magnetic surface. The periodicity condition in the poloidal direction imposes . In linear ballooning theory, the function is found to be separable in and θ for each n and pulsation ω, i.e. . Separability is not always rigorously demonstrated, but is a good proxy. The amplitude is called envelope, usually a localized around an angle called ballooning angle, with a narrow width of the order of , where Lp is a gradient length and the ion thermal gyroradius radius. The ballooning angle measures the departure from up–down symmetry on a given field line. Non zero values of result from symmetry breaking mechanisms, as will be seen later on. This leads to a slow radial dependence of in r. Hence 'fast' radial variations of come from the radial dependence of the safety factor in the exponential argument of equation (6), with a characteristic scale , whereas 'slow' radial variations come from the amplitude , with a mesoscale length . We insist here on the fact that a 'slow' radial scale is not a mean gradient length—the only constraint is that it should be larger than the 'fast' scale .
Calculating a mode envelope in turbulent regime is questionable since a global mode has presumably no enough time to form before turbulent decorrelation occurs, except very close to the threshold. Nevertheless simulations show that an instantaneous ballooning angle can still be identified [16]. We therefore adopt the following representation in non linear regime
The vector is a label for the couple so that the sum over designates a summation over the toroidal wave number n and also an integral over a distribution of angles . The value of the angle depends on time, as shown by simulations [16]. To some extent, this time dependence can be represented by such a distribution. Quantities labeled with a index will be called 'Fourier' harmonics, since can be seen as the Fourier counterpart of the radial variable . Each harmonic depends also slowly on the radius r. We will use in the following a 'slow' radial variable , where ε is a small parameter. All quantities depends on with a typical scale that is larger than a typical turbulent vortex size, but is smaller than a gradient length. The label is omitted to simplify the notations, unless specified otherwise. Equation (8) is in principle not acceptable since it is not periodic in θ. Nevertheless it is a reasonable proxy for a turbulence localized on the low field side of a tokamak (i.e. maximum near ) since is then small for . Hence each Fourier component contains the information on the poloidal localization of fluctuations, and therefore poloidal asymmetries. We postulate a Gaussian form for i.e. , exact in the linear framework
where and are complex numbers, and is a normalizing factor (, and depend on ).
2.1.2. Ballooning angle
The value of the ballooning angle is not easy to determine. In the ballooning formalism, it comes from a calculation of the mode envelope. However, as already mentioned, this calculation is questionable for a turbulent state. An estimate can be found by using a rapid distortion theory [30]. Indeed a structure with initially a zero ballooning angle will acquire an effective radial wave number due to a shear flow that is equal to after a time t. Here is the shear flow rate, i.e. the radial derivative of the mean velocity. Defining a correlation time , one gets the following estimate
Other expressions of have been given in the past [26–29, 31], with similar structures as equation (10) when envelopes are spatially localized. The exact expression of has no real importance here, since we are interested in the relative signs and amplitudes of the momentum fluxes. Sources of symmetry breaking different from shear flow will lead to different expressions of , but will cause correlated radial and parallel wavenumbers.
2.2. E × B drift stress tensor
We estimate now the radial fluxes of poloidal and parallel momentum due to drift, which are proportional to the corresponding components of the Reynolds stress
The distribution function is decomposed in the same way as the potential
where is the adiabatic invariant. The equations equations (11), (12) can then be recast as
where
is the Fourier component of the parallel fluid velocity and
is the Fourier component of the radial component of the drift velocity. The wavenumber is a reference poloidal wave number defined as . The poloidal structure function is defined as (see details in appendix
with the normalization
The amplitude does not depend on time for the structure equation (8). Note that at this stage, the fluxes depend on the poloidal angle and time (plus a slow radial variation in ). A rough estimate of is obtained by a rapid distortion argument [30], i.e. , where is a correlation time. The components of the Reynolds stress can be expanded as Fourier series in the poloidal direction
Using the expressions equations (23), (24), one finds
where is the cyclotron pulsation. The coefficients and measure the distortion of turbulent structures in the radial and parallel directions. They can be formally written as (see appendix
and
The wave number characterizes fluctuation poloidal symmetries. The useful values of are low, typically , and therefore different from the much larger turbulent poloidal wave numbers. The coefficients and can be seen as effective radial and parallel wavenumbers for a given set .
2.2.1. Wave numbers
Calculations in appendix
where and are short notations for and , i.e. components of the coefficients and . These coefficients are related to the parameters used for the turbulent field representation
where
is the magnetic shear calculated at , and . The symbol (resp. ) designates the real (resp. imaginary) part of a complex number.
In view of the field structure equation (9), the real part of must be positive, i.e. , since modes are spatially localized in the poloidal direction. As expected, the average radial and parallel wave numbers are correlated and proportional to , which measures the strength of the symmetry breaking mechanism. However a close inspection of equations (28), (29) indicates that the proportionality coefficients are not the same. For instance, in the limit (non propagative localized mode), it appears that is of the sign of , while is of the sign of . Hence the sign and amplitude of the correlation between the radial and parallel wave numbers is by no way trivial. This relationship has been discussed in [32], in the context of turbulent momentum transport in slab geometry. The components and are even functions of , so that and are even functions of θ. Here calculations are restricted to order one in .
Since we are interested in signs, a rapid distortion theory may not be accurate enough. A quasi-linear theory provides a more precise value of the time and poloidal average of . The calculation is done in appendix
3. Flux of parallel momentum due to magnetic drift
The expression equation (2) of the radial flux of parallel momentum due to magnetic drift shows that a finite flux can only be due to up–down symmetries of the distribution function. The rationale here is as follows: because turbulence is ballooned, the stress tensor is ballooned too, thus leading to poloidal asymmetries of the poloidal flow, and therefore of the electric potential. The resulting distorted distribution function is correlated with the Reynolds stress, and therefore with the radial flux of parallel momentum.
3.1. Generation of poloidal convective cells
Damping of time dependent and poloidally asymmetric flows is not primarily due to collisions, but rather wave particle resonant effects. An estimate can be obtained by solving the linear gyrokinetic equation (see [37, 38] for overviews on gyrokinetic theory)
coupled to the Poisson equation
where N is the unperturbed density, the gyroaverage operator and its inverse, the kinetic diamagnetic velocity, is the unshifted Maxwellian distribution function built with the density N and temperature T, and is the non adiabatic part of the distribution function. The distribution functions F and G are functions of the gyrocenter-center position, the adiabatic invariant μ and parallel velocity . The parallel non linear term has been neglected. The mirror force is also ignored since the quantity of interest is the ion parallel flux, to which passing particles mostly contribute. This nevertheless sets a lower bound in frequency, namely the trapped ion bounce frequency of the order of . It is reminded that the inverse aspect ratio is a small parameter in the present work. To simplify the calculation, we use the long wavelength limit of the gyroaverage operator , i.e. the polarization term is kept, but not the finite Larmor radius (FLR) effects. For non zonal modes , this gives an equation for the potential vorticity , i.e.
where is the ion parallel flux, and the diamagnetic velocity. The vorticity equation (33) has been derived by assuming that the electric potential wavelength is smaller than the gradient lengths of density and temperature. The axisymmetric electric potential n = 0 is Fourier expanded in poloidal angle and time
where since zonal flows are excluded. One can show the following variant of the Taylor identity [39] (see appendix
where
It is reminded that is the reference radius in the neighborhood of which the calculations are performed. Keeping FLR effects would replace one component of the drift velocity by in the stress tensor , i.e. becomes . This does not change the structure of the calculation as long as pressure fluctuations are proportional to potential fluctuations. Implications are discussed in section 4.2. A subsidiary small parameter is introduced at this point, which is the ratio of the thermal Larmor radius to the typical radial wavelength of the turbulent Reynolds stress tensor. Only terms of order are retained in the following.
The divergence of the parallel flux in the vorticity equation (33) can be calculated by computing the axisymmetric components of the perturbed distribution function, which is a linear solution of the gyrokinetic equation (31), i.e.
It appears that the component of the parallel flux divergence is
where is a frequency shift and is the flow damping rate (). The explicit expressions of and are
and
where is a thermal transit frequency. The last term in the vorticity equation (33) cancels out since the divergence of the diamagnetic current term is linear so that its toroidal average is zero. The other term is a particle flux which is equal to zero for an ITG turbulence. The vorticity equation (33) can then be solved by combining the Taylor identity equation (35) with the parallel flow divergence equation (38), thus providing the Fourier components of the electric potential
where the Reynolds stress tensor has been expanded in a Fourier series
and has to be understood as an operator that acts on , . It is reminded that only the leading term in is kept.
The expression of the electric potential Fourier components equation (41) is of central importance since it provides the structure of time-dependent poloidal convective cells which are driven by turbulent eddies. This result calls for several comments:
- 1.It appears that , as expected. However, even if the stress tensor is up–down symmetrical at lowest order in ballooning angle, i.e. (a consequence of the parity of and with ), this is not the case for the potential, which develops up–down asymmetries , due to damping.
- 2.Poloidal convective cells are Landau damped due to their finite poloidal wave number, as expected. However it appears that their damping rate is small at low frequencies , hence favoring their onset and sustainment. It appears that the scaling of the frequency spectrum plays an important role in the generation of parallel flow via the magnetic drift. Obviously the limit is an artifact. In fact two cut-off frequencies appear, that correspond to the approximations made in this calculation: the trapped ion bounce pulsation, and the ion collision frequency. The highest is the ion bounce pulsation, which remains nevertheless lower than the ion transit frequency in the limit of low aspect ratio. Collisions regularize the infrared singularity via resonance broadening.
- 3.The value of plays also an important role. First, it is assumed that is positive, i.e. we consider here the case of radial oscillations of the Reynolds stress. Secondly, the mechanism presented here requires values of that are not too small, i.e. corrugations of the Reynolds stress with a spatial scale of a few gyroradii [11, 13, 16]. This may push the present model to its validity limit, since it was assumed that the variation with r0 correspond to scales larger than the distance between resonant surfaces, itself of the order of the ion gyroradius. A more accurate calculation for is intricate as it requires a more precise procedure to solve the vorticity equation, as done for zonal flows ([40] and references therein, and also for convective cells driven by turbulence at low magnetic shear [23]).
- 4.The ordering of gyrokinetic theory is not broken since the amplitude of the normalized potential remains small, of the order of . This can be verified by using equation (41) , and noting that the turbulence stress tensor scales as , frequencies as so that . Macroscopic flows would lead to very small amplitudes of the potential turbulence on the electric potential with a negligible effect, as discussed in [12]. Numerical simulations suggest that is rather in the range of 0.1 [13].
3.2. Magnetic drift contribution to the flux of parallel momentum
The banana-plateau component of the neoclassical flux of parallel momentum is known to exhibit a Pfirsch-Schlüter scaling and thus vanishes in the collisionless regime [4, 5, 41] (see [42] for an overview). Therefore one is left with the cross-correlation between the magnetic drift and the perturbed distribution function due to the electric potential fluctuations driven by turbulent eddies. This cross-correlation is responsible for a finite flux of parallel momentum. The time dependent radial flux of parallel momentum due to magnetic drift reads
where is the radial component of the magnetic drift velocity. This flux can be calculated using the expression equation (41) of the amplitudes of poloidal convective cells (see appendix
For a ballooned turbulence, the Reynolds stress is ballooned too. If fluctuations are strongly localized on the low field side, the prefactor can be replaced by 1 in equation (44). This can also be demonstrated by using the structure of the stress tensor equation (21) and the expression of with (see equation (26)).
Given the simplicity of equation (44), one may actually wonder whether a more direct calculation of is possible, and indeed it is. Since the parallel velocity plays a subdominant role in the magnetic drift, essentially because the resonant velocity goes like the pulsation, the curvature driven flux of momentum equation (43) can be reformulated as
An integration by part allows a reformulation in terms of the parallel gradient of the parallel flux
Using the vorticity equation (33), and the Taylor identity equation (35), one finds equation (44). It is quite remarkable that this result does not depend on the details of the poloidal convective cells, which act as mediators. An example of turbulence self-organization via the generation of poloidal convective cells is shown in figure 3. This figure comes from simulations run with the GYSELA gyrokinetic code [43]. It turns out that poloidal convective cells play also a role in the interplay between turbulent and neoclassical impurity transport [17].
Download figure:
Standard image High-resolution image4. Comparison of and curvature driven momentum fluxes
4.1. General expressions of momentum fluxes fluxes
The expression equation (44) of the radial flux of parallel momentum due to the magnetic drift can be expressed in the Fourier space using equation (21)
This expression can be compared with the flux of momentum equation (77)
It appears that these two contributions to the total flux of parallel momentum are anti-correlated when and are of the same sign (for positive safety factor ). This is also the condition for anti-correlation of the poloidal and parallel components of the Reynolds stress. The final answer is clearly sensitive to processes that determine the radial and parallel wavenumbers, which are presumably determined by non linear effects. Nevertheless it is interesting to calculate the linear values of and , and the corresponding signs of fluxes.
4.2. Application to ITG driven turbulence
An instructive example is the case of an unstable linear toroidal ITG driven mode. General expressions of and are derived in the hydrodynamic limit in the appendix
The numbers and depend on the normalized frequencies and only, where and are the real and imaginary parts of the linear frequency in the plasma frame (frame where the local radial electric field is zero), and
is a thermal magnetic field curvature frequency. More precisely
and
The coefficient is always positive. These expressions are valid for any sign of and . We now restrict the discussion to the generic case and . These considerations lead to the following estimates of the momentum fluxes
and
For , the two numbers and are of the same order of magnitude. The two fluxes are therefore comparable when . This condition appears as reasonable since the intensity wave number spectrum typically peaks at , provided [44, 45].
Regarding the signs, it appears that for modes that drift in the ion diamagnetic direction , is negative for a weak drive , and positive for strong drive . Since is always positive, this means that and are positively correlated near threshold, and anti-correlated far from threshold. Consequently the two fluxes or parallel momentum are positively correlated for low drive and anti-correlated for strong drive . For modes drifting in the electron diamagnetic direction, , anti-correlation always occurs. For ITG modes, drift in the ion diamagnetic direction is expected, so that anti-correlation is expected only far enough from the instability threshold. It is not clear whether this finding agrees or not with previous numerical findings for momentum transport [13, 16]. Nevertheless it is stressed that the hydrodynamic limit that is being used here, is a rather poor approximation that becomes correct only well above the instability threshold, i.e. for . This is precisely the regime where an agreement with simulations is found, i.e. anticorrelation. In that regard some further analysis of the numerical simulations would be helpful. An encouraging observation though [16] is that the flux of parallel momentum is anti-correlated with the ballooning angle , which suggests positive in view of equation (56). Let us note that adding the FLR effects mentioned in section 3.1 may change the correlation sign. Indeed changing the stress tensor into is roughly equivalent to multiplying in equation (52) by a factor , where is the pressure diamagnetic frequency. Since the later tends to be larger than the mode pulsation for realistic parameters, this may change the sign of the coefficient , and therefore the relative signs of the E × B and magnetic drift components of the momentum flux. As a final note, it is possible (if not likely) that the two terms and are determined by non linear processes and not well captured by a linear analysis.
5. Conclusion
It is shown here that the Reynolds stress tensor generates poloidal asymmetries of the plasma flow due to turbulence ballooning. These poloidal convective cells are weakly damped at low frequency. Their radial scale is dictated by the turbulent Reynolds stress, and their poloidal wave numbers are small. These cells drive up–down asymmetries of the distribution function, which are responsible for a non-zero radial flux of parallel momentum due to the geodesic component of the particle magnetic drift. The entire process requires a symmetry breaking mechanism, for instance a mean shear flow. Since the turbulent Reynolds can be seen as a flux of momentum, it appears that the two components of the radial flux of parallel momentum due to curvature and drift are correlated. This general result comes from a simple relationship between the momentum flux due to magnetic drift and the turbulent Reynolds stress. Although poloidal convective cells do not appear explicitly in this relationship, they play an essential role in this mechanism. An analytic calculation shows that anti-correlation between the components of the turbulent Reynolds stress results in anti-correlation of the two contributions to the flux of parallel momentum that come from and magnetic drifts. A quasi-linear calculation of all quantities, based on ITG linear stability indicates that positive correlation is expected near threshold, and anti-correlation for strong drive. Hence no firm conclusion can be drawn as to the relevance of this mechanism to explain the numerical results. Nevertheless the hydrodynamic limit that has been used is a rather poor representation of the turbulence near threshold. In fact it is valid only well above the stability threshold, i.e. strong drive. This is encouraging since it is the case where anti-correlation is found, in agreement with numerical findings. It is likely that poloidal convective cells generated by turbulence have other consequences on turbulent transport and turbulence. Indeed they may participate in turbulence self-regulation via vortex shearing processes similar to the well documented effect of zonal flows on turbulence.
Appendix A.: Mode structure
The mode structure equation (9) lead to a turbulence intensity that reads
where is a form factor
and
The amplitude in equation (9) has been chosen such that
namely
The poloidal Fourier components of the radial and parallel wave numbers equation (25) and equation (26) can then be recast as
For a strongly ballooned turbulence , and small values of the ballooning angle , one gets at first order in
where . Hee and indicate the real and imaginary parts of a complex number z.
Appendix B.: Quasilinear expression of E × B momentum flux
We start from the expression of the E × B drift contribution to the momentum flux, and readily get its time average equation (12)
To calculate the distribution function versus the potential, we use a ballooning representation. The electric potential is written in the form
where is the ballooning mode. The single term p = 0 is kept for strongly ballooned fluctuations. A similar expansion is used for the non adiabatic part of the distribution function, i.e. . The gyrokinetic equation (31) reads
Here is an operator
where is a shifted poloidal angle. The kinetic magnetic drift frequency is defined as
where
is the thermal ion velocity, and is the poloidal wavenumber— is the safety factor at the reference radius . The definition of the kinetic diamagnetic frequency is the usual one
where
are the density and temperature diamagnetic frequencies and and are the density and pressure gradient lengths calculated at the reference radius . The gyroaverage is fairly well represented by a Bessel function with an argument , i.e. , where is the kinetic ion gyroradius and
The solution of the gyrokinetic equation (67) involves an integro-differential operator that relates to [46, 47]. A formal solution can be written in a Wentzel-Kramers-Brilloin (WKB) sense by dividing the rhs of equation (67) by the resonant term
We assume that the spectral turbulence intensity is of the form
which gives
In the hydrodynamic limit , one gets the result
This result is identical to equation (22) with . Strictly speaking there is also a contribution to the flux of parallel momentum that comes from the axisymmetric perturbations of the potential equation (41) and the distribution function response equation (37) . However this contribution is of second order in ballooning angle and will not be retained here. Hence the flux of momentum is in this peculiar case 'turbulent', i.e. produced by small scale fluctuations.
Appendix C.: Taylor identity
The purpose of this appendix is to demonstrate the identity
where is the vorticity, and the drift velocity. The demonstration is restricted to a geometry of concentric circular surfaces with large aspect ratio . We use the mode structure equation (8)
where , and the dependence on r of is made explicit for clarity throughout this section. The average over ϕ implies that the lhs of equation (78) appears as a sum over the index of operators acting on , namely
where
and is given by a similar expression. This suggests a change of variables for each , namely
The variable η plays the role of a coordinate along the field line, r0 represents the slow variation of the field in the radial electric field, and is a coordinate transverse to the field. We use the following ordering
which allows neglecting the slow variation of the field along the field line. We note that the unit vector along the magnetic field can be written,
where
For a strongly ballooned turbulence . Ignoring the derivatives with respect to η, one finds the following identity
Using , and neglecting again the derivatives along the field lines. Each Fourier harmonics of the vorticity can be written as
where it has been used that the variation of with is faster than the radial variation of . The largest term is the third one in the rhs of equation (89) , but it does not contribute to equation (88) for parity reasons. The largest contribution therefore comes from the second term of the rhs of equation (89) , i.e.
where
The expression within the brackets of the rhs of equation (90) is the component of the stress tensor, as can be verified from equation (14). Equation (90) demonstrates the Taylor identity equation (78).
Appendix D.: Explicit calculation of the flux of parallel momentum
The purpose of this section is to compute the radial flux of parallel momentum due to magnetic drift equation (43). The geodesic component of the magnetic drift velocity reads
where
We compute the momentum flux by using again the linear solution equation (37) of the gyrokinetic equation (31). The frequency Fourier components of the parallel momentum flux equation (43) read
Using equation (38) (or equivalently equations (40), (39)), the integral over the velocity can be expressed as a function of and , namely
We note that it is essential to deal with time dependent perturbations to get a finite flux. Static perturbations resonate at zero parallel velocity and therefore do not contribute to the flux of parallel momentum. This is also the basic reason why the banana-plateau neoclassical flux of parallel momentum is zero. In the present case, this cancellation is prevented by the increase of the potential amplitude at low frequency, as shown by the expression of equation (41). In other words, it is the time derivative of the electric potential that matters. Equation (95) combined with equation (41) leads to the following expression of the time Fourier transform of the momentum flux
where the up–down symmetry of the turbulent Reynolds stress has been used (consequence of the parity of and with ). The next step consists in taking the zero frequency limit of equation (96), which is equivalent to a time average. For frequencies lower than the transit frequency, , the frequency shift is close to unity, , while the damping rate vanishes. The expression equation (96) becomes equation (44).
Appendix E.: Linear calculation of radial and parallel wavenumber
Assuming an electron adiabatic response, and using the solution of the gyro-kinetic equation (74), the following electroneutrality equation is obtained for a non axisymmetric mode
where is the ratio of the ion to electron temperature at . This equation, which describes reasonably well toroidal ITG modes, can be easily extended to non adiabatic electrons [46] (see also [7] and references therein). In the hydrodynamic limit
the electro-neutrality condition equation (97) can be expanded at first order in , and , where is the thermal ion gyroradius. The equation that rules now reads [47]
where
is the η dependent local dispersion relation. The transit frequency is defined and is the pressure diamagnetic frequency
where is the pressure gradient length. The next step consists in expanding the magnetic drift frequency near up to and . After regrouping the various terms, one finds
where Λ is now a local dispersion relation independent of η and
where is the local magnetic drift frequency. Strictly speaking , but the analysis is restricted here to low values of . Calculations are run at first order in , so that the magnetic drift is the one defined in equation (51). The complex number is proportional to the ballooning angle
The smoothest solution of this differential equation is
with the conditions
The two equations (106), (107) provide the values of and of the complex mode frequency . Moreover, the solution equation (105) is acceptable only if , to guarantee mode spatial localization. The mode frequency is written as , where is the real part, and the growth rate. The later is positive above the instability threshold. We will focus on that case. A limit of interest is , and . The dispersion relation then becomes , i.e.
which yields
It is recovered that modes drift linearly in the ion diamagnetic direction for flat density profiles , while they rotate in the electron diamagnetic direction for strong density gradients [47]. We now turn to the quantities of interest, i.e. the radial and parallel wavenumber given by equations (28), (29). Since equation (107) depends on the detail of the local dispersion relation, we will use exclusively equation (106) and leave arbitrary, except at the very end, when making estimates for ITG modes. It appears that is given by the following relation
from which its imaginary part can be deduced
thus leading to
Using the relation , one gets the expression of
The calculation of is somewhat more delicate because of compensation effects. Using the rhs of equation (111), one finds
It appears readily that
where
The condition equation (106) can be used to provide a relation between and , i.e.
from which one gets the following result
An explicit expression of can be found by solving equation (106) , i.e.
It results from equation (119) that changes sign when
provided the rhs is positive. The positivity condition is fulfilled for modes drifting in the ion diamagnetic direction for low wave numbers and positive magnetic shear. This change of sign typically occurs when one moves from the situation of weak drive (near threshold) to a situation of strong drive .
To make these expressions more explicit, we concentrate on the limit of a strong magnetic shear and low wavenumbers . Using the relations equations (118) (120), explicit expressions of , can be found
from which and can be derived
where
The numbers and are always positive, and satisfy the useful relationships
The expression equation (125) of can be obtained directly from equation (106), thus providing a cross-check . Plugging equations (127), (128) in the expressions of and equations (113), (116), one obtains the following relations
where and are numbers that depend on the normalized frequencies and only, namely