Table of contents

A comprehensive review of zonal flow phenomena in plasmas is presented. While the emphasis is on zonal flows in laboratory plasmas, planetary zonal flows are discussed as well. The review presents the status of theory, numerical simulation and experiments relevant to zonal flows. The emphasis is on developing an integrated understanding of the dynamics of drift wave–zonal flow turbulence by combining detailed studies of the generation of zonal flows by drift waves, the back-interaction of zonal flows on the drift waves, and the various feedback loops by which the system regulates and organizes itself. The implications of zonal flow phenomena for confinement in, and the phenomena of fusion devices are discussed. Special attention is given to the comparison of experiment with theory and to identifying directions for progress in future research.

The behaviour of fast ions, generated by high energy neutral beams, is investigated in the MAST spherical tokamak using a spatially scanning neutral particle analyser (NPA), which has views that detect charge exchange neutrals from both passing and trapped ions. Almost featureless neutral particle spectra are produced, characteristic of thermalization of the neutral beam energy components, when viewing the charge exchange neutrals from passing fast ions in both L- and H-modes. Using NPA views with low tangency radius, R_NPA, both passing and trapped fast ions can be monitored. In L-mode, similar featureless spectra are observed in contrast to H-mode, where they are highly structured. All these aspects have been faithfully reproduced using particle transport codes showing, in particular, that the H-mode spectra found at low R_NPA arise from a small fraction of poorly confined trapped fast ions at the edge, an effect enhanced by the relatively high edge ion and neutral densities.

A fluid model of the magnetic presheath in a turbulent boundary plasma is presented. Turbulent transport corrections of the classical three-dimensional fluid transport equations, which can be used to study magnetic presheaths in various geometries, are derived by means of the ensemble averaging procedure from the statistical theory of plasma turbulence. Then, the magnetic presheath in front of an infinite plane surface is analysed in detail. The linearized planar magnetic presheath equations are applied to the plasma-presheath–magnetic-presheath boundary (i.e. the magnetic presheath edge), whereas the original non-linear planar magnetic presheath equations are used for the entire magnetic presheath, allowing for various sets of experimentally relevant free model parameters to be applied. Important new results of this study are, among others, new expressions for the fluid Bohm criterion at the Debye sheath edge and for the ion flux density perpendicular to the wall. These new results, which exhibit corrections due to the turbulent charged particle transport, can qualitatively explain the fact that whenever the angle between the magnetic field and the wall is very small (i.e. several degrees) or zero, electric currents, measured by Langmuir probes in the boundary regions of nuclear fusion devices and in various low-temperature plasmas, are anomalously enhanced in comparison with those expected or predicted by other theoretical models.

The plasma edge MHD stability is analysed for several JET discharges in the diagnostic optimized configuration. The stability analysis of Type I ELMy plasmas shows how after an edge localized mode (ELM) crash the plasma edge is deep in the stable region against low- to intermediate-n peeling–ballooning modes. As the pressure gradient steepens and the edge current builds up, the plasma reaches the low- to intermediate-n peeling–ballooning mode stability boundary just before the ELM crash. Increasing the plasma fuelling by gas puffing makes the second stability access against high-n ballooning modes narrower until it closes completely and the ELMs change from Type I to Type III. Reducing the plasma heating has a similar effect. Increasing the safety factor at the plasma edge improves the stability against low- to intermediate-n modes allowing steeper pressure gradients to develop before an ELM crash.