Turbulence in hot magnetized plasmas is shown to generate permeable localized transport barriers that globally organize into the so-called "ExB staircase" [G. Dif-Pradalier et al., Phys. Rev. E, 82, 025401(R) (2010)]. Its domain of existence and dependence with key plasma parameters is discussed theoretically. Based on these predictions, staircases are observed experimentally in the Tore Supra tokamak by means of high-resolution fast-sweeping X-mode reflectometry. This observation strongly emphasizes the critical role of mesoscale self-organization in plasma turbulence and may have far-reaching consequences for turbulent transport models and their validation. A puzzling result in recent years in plasma turbulence has arguably been the discovery of the quasiregular pattern of E × B flows and interacting avalanches that we have come to call the "E × B staircase," or the "plasma staircase" in short [1]. This structure may be defined as a spontaneously formed, self-organizing pattern of quasiregular, long-lived, localized shear flow and stress layers coinciding with similarly long-lived pressure corrugations and interspersed between regions of turbulent avalanching. The plasma staircase exemplifies how a systematic organization of turbulent fluctuations may lead to the onset of strongly correlated flows on magnetic flux surfaces.Flow patterning is a prominent topic in many fluidrelated systems and hot magnetized plasmas are no exception to that. In fact the "staircase" name is borrowed from the vast literature in planetary flows motivated by the desire to explain the banded structure of observed atmospheres in our Solar System-including Earth [2] or Jupiter [3]-and of terrestrial oceans [4]. Just as in the geophysical or astrophysical systems where the planetary staircase strongly influences the general circulation, the plasma staircase plays an important role in organizing the heat transport [1]: avalanches and the staircase interplay, statistically interrupting at mesoscales the long-range radial avalanching that could otherwise expand over the whole system. The nonlocal heat transport thus remains contained at the mesoscale staircase step spacing, resulting in a beneficial scaling of confinement with machine size. This flow patterning is primarily a spontaneous mean zonal shear patterning. "Zonal" denotes the axisymmetric n ¼ m ¼ 0 component of the E × B flows [5], n and m respectively being the toroidal and poloidal mode numbers while "mean" refers to the ensemble-averaged part of the zonal flows. Remarkably, the plasma spontaneously generates robust shear patterns that endure despite the strong background turbulence and retain their coherence over long (several milliseconds) to very long (hundreds of milliseconds) periods of time. The results presented throughout this Letter are based on state-of-the-art flux-driven gyrokinetic [6] computations using the GYSELA code [7] with realistic tokamak plasma parameters. Systematic features of the plasma staircase can be inferred from extensive computational scans, see ...

International audienceIn a dedicated collisionality scan in Tore Supra, the geodesic acoustic mode (GAM) is detected and identified with the Doppler backscattering technique. Observations are compared to the results of a simulation with the gyrokinetic code GYSELA. We found that the GAM frequency in experiments is lower than predicted by simulation and theory. Moreover, the disagreement is higher in the low collisionality scenario. Bursts of non harmonic GAM oscillations have been characterized with filtering techniques, such as the Hilbert-Huang transform. When comparing this dynamical behaviour between experiments and simulation, the probability density function of GAM amplitude and the burst autocorrelation time are found to be remarkably similar. In the simulation, where the radial profile of GAM frequency is continuous, we observed a phenomenon of radial phase mixing of the GAM oscillations, which could influence the burst autocorrelation time

Abstract. For a comprehensive comparison with theoretical models and advanced numerical turbulence simulations, a large spectrum of fluctuation parameters were measured on the devices ASDEX Upgrade, TCV, and Tore-Supra. Radial profiles of scale-resolved turbulence levels in H-mode discharges are measured and compared with GENE simulations in the transition range from ion-temperature-gradient to trappedelectron-mode turbulence. Correlation reflectometry is used to study the microscopic structure of turbulence and GAMs in discharges where poloidal flow damping was varied by means of variations of the shape of the poloidal plasma cross-section and collisionality. The aspects of turbulence spreading and non-local transport in response of local changes in density and temperature gradients are studied during phases with ECRH power modulation by means of radial correlation measurements and an ultrafast reflectometer. Full-wave codes and synthetic diagnostics are applied for the interpretation of the data.

Simultaneous measurements of binormal velocity of density fluctuations using two separate Doppler backscattering systems at the low field side and at the top of the plasma show significant poloidal asymmetry. The measurements are performed in the core region between the radii 0.7 < ρ  <  0.95, over a limited number of L-mode discharges covering a wide range of plasma conditions in the Tore Supra tokamak. A possible generation mechanism by the ballooned structure of the underlying turbulence, in the form of convective cells, is proposed for explaining the observation of these poloidally asymmetric mean flows.

Various indications on the weakly nonlocal character of turbulent plasma transport both from experimental fluctuation measurements from Tore Supra and observations from the full-f, flux-driven gyrokinetic code GYSELA are reported. A simple Fisher equation model of this weakly nonlocal dynamics can be formulated in terms of an evolution equation for the turbulent entropy density, which contains the basic phenomenon of radial turbulence spreading in addition to avalanche-like dynamics via coupling to profile modulations. A derivation of this model, which contains the so-called beach effect, a diffusive and convective flux components for the flux of turbulence intensity, in addition to linear group propagation is given, starting from the drift-kinetic equation. The proposed model has the form of a transport equation for turbulence intensity, and may be considered as an addition to transport modelling. The kinetic fluxes given, can be computed using model closures, or local gyrokinetics. The model is also used in a particular setup that represents the near edge region as a relatively stable zone between the core and edge region where the energy injection is locally more substantial. It is observed that with constant, physical coefficients, the model gives a convincing qualitative profile of fluctuation intensity when the turbulence is coming from the core region with either a group velocity or a convective flux.

Gyrokinetic simulations of fusion plasmas give extensive information in 5D on turbulence and transport. This paper highlights a few of these challenging physics in global, flux driven simulations using experimental inputs from Tore Supra shot TS45511. The electrostatic gyrokinetic code GYSELA is used for these simulations. The 3D structure of avalanches indicate that these structures propagate radially at localised toroidal angles and then expand along the field line at sound speed to form the filaments. Analysing the poloidal mode structure of the potential fluctuations (at a given toroidal location), one finds that the low modes m = 0 and m = 1 exhibit a global structure; the magnitude of the m = 0 mode is much larger than that of the m = 1 mode. The shear layers of the corrugation structures are thus found to be dominated by the m = 0 contribution, that are comparable to that of the zonal flows. This global mode seems to localise the m = 2 mode but has little effect on the localisation of the higher mode numbers. However when analysing the pulsation of the latter modes one finds that all modes exhibit a similar phase velocity, comparable to the local zonal flow velocity. The consequent dispersion like relation between the modes pulsation and the mode numbers provides a means to measure the zonal flow. Temperature fluctuations and the turbulent heat flux are localised between the corrugation structures. Temperature fluctuations are found to exhibit two scales, small fluctuations that are localised by the corrugation shear layers, and appear to bounce back and forth radially, and large fluctuations, also readily observed on the flux, which are associated to the disruption of the corrugations. The radial ballistic velocity of both avalanche events if of the order of 0.5ρ * c0 where ρ * = ρ0/a, a being the tokamak minor radius and ρ0 being the characteristic Larmor radius, ρ0 = c0/Ω0. c0 is the reference ion thermal velocity and Ω0 = qiB0/mi the reference ion Larmor frequency for the characteristic amplitude of the magnetic field B0, qi and mi being respectively the ion charge and mass. The electric drift velocity is also found to exhibit a poloidal pattern, with maximum amplitude of the fluctuations either in the top or in the bottom regions of the machine depending on the sign of the zonal flow shear. This effect is found to be correlated to the stopping capability of the corrugation structures. The neoclassical properties stemming from the trapped particle drifts lead to large distortion of the distribution function. As expected, these prevail at the outer part of the simulation region despite the large collisionality. The distribution function fluctuations appear to be aligned along the v = constant lines at constant poloidal angle. A specific symmetry is observed regarding the interplay of turbulence with the trapped-passing region.

Starting from the Braginskii equations, relevant for the tokamak edge region, a complete set of nonlinear equations for the geodesic acoustic modes (GAM) has been derived which includes collisionality, plasma beta and external sources of particle, momentum and heat. Local linear analysis shows that the GAM frequency increases with collisionality at low radial wave number k r and decreases at high k r . GAM frequency also decreases with plasma beta. Radial profiles of GAM frequency for two Tore Supra shots, which were part of a collisionality scan, are compared with these calculations. Discrepency between experiment and theory is observed, which seems to be explained by a finite k r for the GAM when flux surface averaged density n and temperature T are assumed to vanish. It is shown that this agreement is incidental and self-consistent inclusion of n and T responses enhances the disagreement more with k r at high k r . So the discrepancy between the linear GAM calculation, (which persist also for more "complete" linear models such as gyrokinetics) can probably not be resolved by simply adding a finite k

In parallel to the direct contribution to the procurement phase of ITER and Broader Approach, CEA has initiated research & development programmes, accompanied by experiments together with a significant modelling effort, aimed at ensuring robust operation, plasma performance, as well as mitigating the risks of the procurement phase. This overview reports the latest progress in both fusion science and technology including many areas, namely the mitigation of superconducting magnet quenches, disruption-generated runaway electrons, edge-localized modes (ELMs), the development of imaging surveillance, and heating and current drive systems for steady-state operation. The WEST (W Environment for Steady-state Tokamaks) project, turning Tore Supra into an actively cooled W-divertor platform open to the ITER partners and industries, is presented.