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.

/pdf/zonal-flows-in-plasma-a-review-3byloag115.pdf

Zonal flows in plasma—a review

The impact of radially sheared poloidal flows on ambient edge turbulence in tokamaks is investigated analytically. In the regime where poloidal shearing exceeds turbulent radial scattering, a hybrid time scale weighted toward the former is found to govern the decorrelation process. The coupling between radial and poloidal decorrelation results in a suppression of the turbulence below its ambient value. The turbulence quench mechanism is found to be insensitive to the sign of either the radial electric field or its shear.

Influence of sheared poloidal rotation on edge turbulence

Collisionless electron-temperature-gradient-driven (ETG) turbulence in toroidal geometry is studied via nonlinear numerical simulations To this aim, two massively parallel, fully gyrokinetic Vlasov codes are used, both including electromagnetic effects Somewhat surprisingly, and unlike in the analogous case of ion-temperature-gradient-driven (ITG) turbulence, we find that the turbulent electron heat flux is significantly underpredicted by simple mixing length estimates in a certain parameter regime (ŝ∼1, low α) This observation is directly linked to the presence of radially highly elongated vortices (“streamers”) which lead to very effective cross-field transport The simulations therefore indicate that ETG turbulence is likely to be relevant to magnetic confinement fusion experiments

Electron temperature gradient driven turbulence

Progress, since the ITER Physics Basis publication (ITER Physics Basis Editors et al 1999 Nucl. Fusion 39 2137–2664), in understanding the processes that will determine the properties of the plasma edge and its interaction with material elements in ITER is described. Experimental areas where significant progress has taken place are energy transport in the scrape-off layer (SOL) in particular of the anomalous transport scaling, particle transport in the SOL that plays a major role in the interaction of diverted plasmas with the main-chamber material elements, edge localized mode (ELM) energy deposition on material elements and the transport mechanism for the ELM energy from the main plasma to the plasma facing components, the physics of plasma detachment and neutral dynamics including the edge density profile structure and the control of plasma particle content and He removal, the erosion of low- and high-Z materials in fusion devices, their transport to the core plasma and their migration at the plasma edge including the formation of mixed materials, the processes determining the size and location of the retention of tritium in fusion devices and methods to remove it and the processes determining the efficiency of the various fuelling methods as well as their development towards the ITER requirements. This experimental progress has been accompanied by the development of modelling tools for the physical processes at the edge plasma and plasma–materials interaction and the further validation of these models by comparing their predictions with the new experimental results. Progress in the modelling development and validation has been mostly concentrated in the following areas: refinement in the predictions for ITER with plasma edge modelling codes by inclusion of detailed geometrical features of the divertor and the introduction of physical effects, which can play a major role in determining the divertor parameters at the divertor for ITER conditions such as hydrogen radiation transport and neutral–neutral collisions, modelling of the ion orbits at the plasma edge, which can play a role in determining power deposition at the divertor target, models for plasma–materials and plasma dynamics interaction during ELMs and disruptions, models for the transport of impurities at the plasma edge to describe the core contamination by impurities and the migration of eroded materials at the edge plasma and its associated tritium retention and models for the turbulent processes that determine the anomalous transport of energy and particles across the SOL. The implications for the expected performance of the reference regimes in ITER, the operation of the ITER device and the lifetime of the plasma facing materials are discussed.

/pdf/chapter-4-power-and-particle-control-3l68lj0fxw.pdf

Chapter 4: Power and particle control

The understanding and predictive capability of transport physics and plasma confinement is reviewed from the perspective of achieving reactor-scale burning plasmas in the ITER tokamak, for both core and edge plasma regions. Very considerable progress has been made in understanding, controlling and predicting tokamak transport across a wide variety of plasma conditions and regimes since the publication of the ITER Physics Basis (IPB) document (1999 Nucl. Fusion 39 2137-2664). Major areas of progress considered here follow. (1) Substantial improvement in the physics content, capability and reliability of transport simulation and modelling codes, leading to much increased theory/experiment interaction as these codes are increasingly used to interpret and predict experiment. (2) Remarkable progress has been made in developing and understanding regimes of improved core confinement. Internal transport barriers and other forms of reduced core transport are now routinely obtained in all the leading tokamak devices worldwide. (3) The importance of controlling the H-mode edge pedestal is now generally recognized. Substantial progress has been made in extending high confinement H-mode operation to the Greenwald density, the demonstration of Type I ELM mitigation and control techniques and systematic explanation of Type I ELM stability. Theory-based predictive capability has also shown progress by integrating the plasma and neutral transport with MHD stability. (4) Transport projections to ITER are now made using three complementary approaches: empirical or global scaling, theory-based transport modelling and dimensionless parameter scaling (previously, empirical scaling was the dominant approach). For the ITER base case or the reference scenario of conventional ELMy H-mode operation, all three techniques predict that ITER will have sufficient confinement to meet its design target of Q = 10 operation, within similar uncertainties.

https://iopscience.iop.org/article/10.1088/0029-5515/47/6/S02/pdf

Chapter 2: Plasma confinement and transport

Tokamak plasmas are inherently comprised of multiple ion species. This is due to wall-bred impurities and, in future reactors, will result from fusion-born alpha particles. Relatively small densities nI, of highly charged non-hydrogenic impurities can strongly influence plasma transport properties whenever . The determination of the complete neoclassical Onsager matrix for a toroidally confined multispecies plasma, which provides the linear relation between the surface averaged radial fluxes and the thermodynamic forces (i.e. gradients of density and temperature, and the parallel electric field), is reviewed. A closed set of one-dimensional moment equations is presented for the time evolution of thermodynamic and magnetic field quantities which results from collisional transport of the plasma and two-dimensional motion of the magnetic flux surface geometry. The effects of neutral-beam injection on the equilibrium and transport properties of a toroidal plasma are consistently included.

https://iopscience.iop.org/article/10.1088/0029-5515/21/9/003/pdf

Neoclassical transport of impurities in tokamak plasmas

Three-dimensional free boundary calculations using a spectral Green's function method

A multi-species fluid model is described for the steady state parallel and radial force balance equations in axisymmetric tokamak plasmas. The bootstrap current, electrical resistivity, and particle and heat fluxes are evaluated in terms of the rotation velocities and friction and viscosity coefficients. A recent formulation of the neoclassical plasma viscosity for arbitrary shape and aspect ratio (including the unity aspect ratio limit), arbitrary collisionality, and orbit squeezing from strong radial electric fields is used to illustrate features of the model. The bootstrap current for the very low aspect ratio National Spherical Torus Experiment [J. Spitzer et al., Fusion Technol. 30, 1337 (1996)] is compared with other models; the largest differences occur near the plasma edge from treatment of the collisional contributions. The effects of orbit squeezing on bootstrap current, thermal and particle transport, and poloidal rotation are illustrated for an enhanced reverse shear plasma in the Tokamak Fusion Test Reactor [D. Meade and the TFTR Group, Plasma Physics and Controlled Nuclear Fusion Research, 1990 (International Atomic Energy Agency, Vienna, 1991), Vol. I, p. 9]. Multiple charge states of impurities are incorporated using the reduced ion charge state formalism for computational efficiency. Because the force balance equations allow for inclusion of external momentum and heat sources and sinks they can be used for general plasma rotation studies while retaining the multi-species neoclassical effects.

Bootstrap current and neoclassical transport in tokamaks of arbitrary collisionality and aspect ratio

A variational principle is developed for the linearized drift‐kinetic, Fokker–Planck equation, from which both upper and lower bounds for neoclassical transport coefficients can be calculated for plasmas in three‐dimensional toroidal confinement geometries. These bounds converge monotonically with the increasing phase‐space dimensionality of the assumed trial function. This property may be used to identify those portions of phase space that make dominant contributions to the transport process. A computer code based on this principle has been developed that uses Fourier–Legendre expansions for the poloidal, toroidal, and pitch‐angle dependences of the distribution function. Numerical calculations of transport coefficients for a plasma in the TJ‐II flexible heliac [Nucl. Fusion 28, 157 (1988)] are used to demonstrate the application of this procedure.

S.P. Hirshman

Papers

Neoclassical transport of impurities in tokamak plasmas

Three-dimensional free boundary calculations using a spectral Green's function method

Bootstrap current and neoclassical transport in tokamaks of arbitrary collisionality and aspect ratio

Variational bounds for transport coefficients in three-dimensional toroidal plasmas

Preconditioned descent algorithm for rapid calculations of magnetohydrodynamic equilibria