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

One of the scientific success stories of fusion research over the past decade is the development of the ExB shear stabilization model to explain the formation of transport barriers in magnetic confinement devices. This model was originally developed to explain the transport barrier formed at the plasma edge in tokamaks after the L (low) to H (high) transition. This concept has the universality needed to explain the edge transport barriers seen in limiter and divertor tokamaks, stellarators, and mirror machines. More recently, this model has been applied to explain the further confinement improvement from H (high)-mode to VH (very high)-mode seen in some tokamaks, where the edge transport barrier becomes wider. Most recently, this paradigm has been applied to the core transport barriers formed in plasmas with negative or low magnetic shear in the plasma core. These examples of confinement improvement are of considerable physical interest; it is not often that a system self-organizes to a higher energy state with reduced turbulence and transport when an additional source of free energy is applied to it. The transport decrease that is associated with ExB velocity shear effects also has significant practical consequences for fusion research. The fundamental physics involved in transport reduction is the effect of ExB shear on the growth, radial extent and phase correlation of turbulent eddies in the plasma. The same fundamental transport reduction process can be operational in various portions of the plasma because there are a number ways to change the radial electric field Er. An important theme in this area is the synergistic effect of ExB velocity shear and magnetic shear. Although the ExB velocity shear appears to have an effect on broader classes of microturbulence, magnetic shear can mitigate some potentially harmful effects of ExB velocity shear and facilitate turbulence stabilization.

/pdf/effects-of-e-b-velocity-shear-and-magnetic-shear-on-5f3wdz73f7.pdf

Effects of E×B velocity shear and magnetic shear on turbulence and transport in magnetic confinement devices

Progress in the area of MHD stability and disruptions, since the publication of the 1999 ITER Physics Basis document (1999 Nucl. Fusion 39 2137-2664), is reviewed. Recent theoretical and experimental research has made important advances in both understanding and control of MHD stability in tokamak plasmas. Sawteeth are anticipated in the ITER baseline ELMy H-mode scenario, but the tools exist to avoid or control them through localized current drive or fast ion generation. Active control of other MHD instabilities will most likely be also required in ITER. Extrapolation from existing experiments indicates that stabilization of neoclassical tearing modes by highly localized feedback-controlled current drive should be possible in ITER. Resistive wall modes are a key issue for advanced scenarios, but again, existing experiments indicate that these modes can be stabilized by a combination of plasma rotation and direct feedback control with non-axisymmetric coils. Reduction of error fields is a requirement for avoiding non-rotating magnetic island formation and for maintaining plasma rotation to help stabilize resistive wall modes. Recent experiments have shown the feasibility of reducing error fields to an acceptable level by means of non-axisymmetric coils, possibly controlled by feedback. The MHD stability limits associated with advanced scenarios are becoming well understood theoretically, and can be extended by tailoring of the pressure and current density profiles as well as by other techniques mentioned here. There have been significant advances also in the control of disruptions, most notably by injection of massive quantities of gas, leading to reduced halo current fractions and a larger fraction of the total thermal and magnetic energy dissipated by radiation. These advances in disruption control are supported by the development of means to predict impending disruption, most notably using neural networks. In addition to these advances in means to control or ameliorate the consequences of MHD instabilities, there has been significant progress in improving physics understanding and modelling. This progress has been in areas including the mechanisms governing NTM growth and seeding, in understanding the damping controlling RWM stability and in modelling RWM feedback schemes. For disruptions there has been continued progress on the instability mechanisms that underlie various classes of disruption, on the detailed modelling of halo currents and forces and in refining predictions of quench rates and disruption power loads. Overall the studies reviewed in this chapter demonstrate that MHD instabilities can be controlled, avoided or ameliorated to the extent that they should not compromise ITER operation, though they will necessarily impose a range of constraints.

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

Chapter 3: MHD stability, operational limits and disruptions

The ITER Physics Basis presents and evaluates the physics rules and methodologies for plasma performance projections, which provide the basis for the design of a tokamak burning plasma device whose goal is to demonstrate the scientific and technological feasibility of fusion energy for peaceful purposes. This Chapter summarizes the physics basis for burning plasma projections, which is developed in detail by the ITER Physics Expert Groups in subsequent chapters. To set context, the design guidelines and requirements established in the report of ITER Special Working Group 1 are presented, as are the specifics of the tokamak design developed in the Final Design Report of the ITER Engineering Design Activities, which exemplifies burning tokamak plasma experiments. The behaviour of a tokamak plasma is determined by the interaction of many diverse physics processes, all of which bear on projections for both a burning plasma experiment and an eventual tokamak reactor. Key processes summarized here are energy and particle confinement and the H-mode power threshold; MHD stability, including pressure and density limits, neoclassical islands, error fields, disruptions, sawteeth, and ELMs; power and particle exhaust, involving divertor power dispersal, helium exhaust, fuelling and density control, H-mode edge transition region, erosion of plasma facing components, tritium retention; energetic particle physics; auxiliary power physics; and the physics of plasma diagnostics. Summaries of projection methodologies, together with estimates of their attendant uncertainties, are presented in each of these areas. Since each physics element has its own scaling properties, an integrated experimental demonstration of the balance between the combined processes which obtains in a reactor plasma is inaccessible to contemporary experimental facilities: it requires a reactor scale device. It is argued, moreover, that a burning plasma experiment can be sufficiently flexible to permit operation in a steady state mode, with non-inductive plasma current drive, as well as in a pulsed mode where current is inductively driven. Overall, the ITER Physics Basis can support a range of candidate designs for a tokamak burning plasma facility. For each design, there will remain a significant uncertainty in the projected performance, but the projection methodologies outlined here do suffice to specify the major parameters of such a facility and form the basis for assuring that its phased operation will return sufficient information to design a prototype commercial fusion power reactor, thus fulfilling the goal of the ITER project.

https://iopscience.iop.org/article/10.1088/0029-5515/39/12/301/pdf

Chapter 1: Overview and summary

Nonlinear gyrokinetic equations play a fundamental role in our understanding of the long-time behavior of strongly magnetized plasmas. The foundations of modern nonlinear gyrokinetic the- ory are based on three important pillars: (1) a gyrokinetic Vlasov equation written in terms of a gyrocenter Hamiltonian with quadratic low-frequency ponderomotive-like terms; (2) a set of gyrokinetic Maxwell (Poisson-Ampere) equations written in terms of the gyrocenter Vlasov dis- tribution that contain low-frequency polarization (Poisson) and magnetization (Ampere) terms derived from the quadratic nonlinearities in the gyrocenter Hamiltonian; and (3) an exact energy conservationlaw for the gyrokineticVlasov-Maxwell equations that includes all the relevant linear and nonlinear coupling terms. The foundations of nonlinear gyrokinetic theory are reviewed with an emphasis on the rigorous applications of Lagrangian and Hamiltonian Lie-transform perturba- tion methods used in the variationalderivationof nonlineargyrokineticVlasov-Maxwell equations. The physical motivations and applications of the nonlinear gyrokinetic equations, which describe the turbulent evolution of low-frequency electromagnetic fluctuations in a nonuniform magnetized plasmas with arbitrary magnetic geometry, are also discussed.

/pdf/foundations-of-nonlinear-gyrokinetic-theory-3l78lnlpai.pdf

Foundations of nonlinear gyrokinetic theory

A physically comprehensive and theoretically based transport model tuned to three-dimensional (3-D) ballooning mode gyrokinetic instabilities and gyrofluid nonlinear turbulence simulations is formulated with global and local magnetic shear stabilization and E×B rotational shear stabilization. Taking no fit coefficients from experiment, the model is tested against a large transport profile database with good agreement. This model is capable of describing enhanced core confinement transport barriers in negative central shear discharges based on rotational shear stabilization. The model is used to make ignition projections from relative gyroradius scaling discharges.

A gyro-Landau-fluid transport model

A tokamak equilibrium model, local to a flux surface, is introduced which is completely described in terms of nine parameters including aspect ratio, elongation, triangularity, and safety factor. By allowing controlled variation of each of these nine parameters, the model is particularly suitable for localized stability studies such as those carried out using the ballooning mode representation of the gyrokinetic equations.

Noncircular, finite aspect ratio, local equilibrium model

The method of Hammett and Perkins [Phys. Rev. Lett. 64, 3019 (1990)] to model Landau damping has been recently applied to the moments of the gyrokinetic equation with curvature drift by Waltz, Dominguez, and Hammett [Phys. Fluids B 4, 3138 (1992)]. The higher moments are truncated in terms of the lower moments (density, parallel velocity, and parallel and perpendicular pressure) by modeling the deviation from a perturbed Maxwellian to fit the kinetic response function at all values of the kinetic parameters: k∥vth/ω, b=(k⊥ρ)2/2, and ωD/ω. Here the resulting gyro‐Landau fluid equations are applied to the simulation of ion temperature gradient (ITG) mode turbulence in toroidal geometry using a novel three‐dimensional (3‐D) nonlinear ballooning mode representation. The representation is a Fourier transform of a field line following basis (ky’,kx’,z’) with periodicity in toroidal and poloidal angles. Particular emphasis is given to the role of nonlinearly generated n=0 (ky’ = 0, kx’ ≠ 0) ‘‘radial modes’’ in s...

Toroidal gyro‐Landau fluid model turbulence simulations in a nonlinear ballooning mode representation with radial modes

Gyrokinetic simulation of tokamak transport has evolved sufficiently to allow direct comparison of numerical results with experimental data. It is to be emphasized that only with the simultaneous inclusion of many distinct and complex effects can this comparison realistically be made. Until now, numerical studies of tokamak microturbulence have been restricted to either (a) flux tubes or (b) electrostatic fluctuations. Using a newly developed global electromagnetic solver, we have been able to recover via direct simulation the Bohm-like scaling observed in DIII-D L-mode discharges. We also match, well within experimental uncertainty, the measured energy diffusivities.

R. E. Waltz

Papers

A gyro-Landau-fluid transport model

Noncircular, finite aspect ratio, local equilibrium model

Toroidal gyro‐Landau fluid model turbulence simulations in a nonlinear ballooning mode representation with radial modes

Anomalous Transport Scaling in the DIII-D Tokamak Matched by Supercomputer Simulation

An Eulerian gyrokinetic-Maxwell solver