(TH/P846) Global Gyrokinetic Simulations of ρ^{*} and ν^{*} Scalings of Turbulent Transport
Y. Sarazin^{1)},
G. DifPradalier^{1)},
V. Grandgirard^{1)},
P. Angelino^{1)},
X. Garbet^{1)},
Ph. Ghendrih^{1)},
R. Belaouar^{2)},
N. Crouseilles^{2)},
G. Latu^{2)},
E. Sonnendrücker^{2)},
S. Jolliet^{3)},
B.F. McMillan^{3)},
T.M. Tran^{3)},
L. Villard^{3)}
^{1)} CEA, IRFM, F13108 SaintPaullezDurance, France
^{2)} IRMAUniversité Louis Pasteur and INRIA Lorraine, 7 rue René Descartes, 67084 Strasbourg Cedex, France
^{3)} CRPP, Association EuratomConfédération Suisse, EPFL, 1015 Lausanne, Switzerland
Abstract. Predicting the confinement time in next step controlled fusion devices such as Iter requires to know how the turbulent transport depends on key scaling parameters, more specifically on ρ^{*}, the ratio of the ion Larmor radius to the machine size a, and on the collisionality nustar. Such a scan is performed with the global gyrokinetic code GYSELA, which solves the standard gyrokinetic equation for the full distribution function of the ions with adiabatic electrons. The code achieves more than 80% efficiency on 4 096 processors.
Gyrokinetic simulations have predicted the transition from Bohm to
gyroBohm scaling when decreasing ρ^{*}, the open issues remaining
the precise value of this transition and the underlying physics.
Such a scan is performed with GYSELA. While the correlation time is
roughly independent of ρ^{*}, the correlation length λ_{c} looks
consistent with the gyroBohm scaling at small ρ^{*} and well above
the linear threshold. Conversely, at large values of ρ^{*} and
close to the threshold, λ_{c} exhibits a dependence on the system size,
consistent with the Bohm scaling. When examining the latter regime,
one finds that the gyroBohm scaling is more relevant at small
ρ^{*}, suggesting that the transition value between Bohm and
gyroBohm scaling is close to
ρ^{*} = 0.01  0.02 and might depend on
the distance to the threshold.
Interplay between neoclassical and turbulent transport is addressed with ν^{*} scans. A FokkerPlanck operator has been added to GYSELA, acting in the parallel velocity space only. It is shown analytically that the neoclassical transport coefficients are recovered with such an ionion collision operator in the three collisional regimes, banana, plateau and PfirschSchlütter. Below the critical Ion Temperature Gradient for the onset of turbulence, the poloidal equilibrium velocity is found to reverse sign when increasing ν^{*} from the banana to the PfirschSchlütter regime, in good agreement with theoretical predictions. Collisions modify the response of the zonal flows, as well as the turbulent heat flux. The dependence of the turbulent heat flux on ν^{*} is presently being analyzed for these simulations. It is found to combine the reduction of the linear growth rate, the damping effect on zonal flows as well as the spreading into the stable boundary regions and the profile relaxation.
Full paper available (PDF)
