Dynamics of zonal flows in helical systems

How the birth and death of shear layers determine confinement evolution: from the L → H transition to the density limit

Electric field profile structure-especially its shear-is a natural order parameter for the edge plasma, and characterizes confinement regimes ranging from the H-mode (Wagner et al. 1982 Phys. Rev. Lett. 49, 1408-1412 (doi:10.1103/PhysRevLett.49.1408)) to the density limit (DL) (Greenwald et al. 1988 Nucl. Fusion 28, 2199-2207 (doi:10.1088/0029-5515/28/12/009)). The theoretical developments and lessons learned during 40 years of H-mode studies (Connor & Wilson 1999 Plasma Phys. Control. Fusion 42, R1-R74 (doi:10.1088/0741-3335/42/1/201); Wagner 2007 Plasma Phys. Control. Fusion 49, B1-B33 (doi:10.1088/0741-3335/49/12b/s01)) are applied to the shear layer collapse paradigm (Hong et al. 2017 Nucl. Fusion 58, 016041 (doi:10.1088/1741-4326/aa9626)) for the onset of DL phenomena. Results from recent experiments on edge shear layers and DL phenomenology are summarized and discussed in the light of L [Formula: see text] H transition physics. The theory of shear layer collapse is then developed. We demonstrate that shear layer physics captures both the well known current (Greenwald) scaling of the DL (Greenwald 2002 Plasma Phys. Control. Fusion 44, R27-R53 (doi:10.1088/0741-3335/44/8/201); Greenwald et al. 2014 Phys. Plasmas 21, 110501 (doi:10.1063/1.4901920)), as well as the emerging power scaling (Zanca, Sattin, JET Contributors 2019 Nucl. Fusion 59, 126011 (doi:10.1088/1741-4326/ab3b31)). The derivation of the power scaling theory exploits an existing model, originally developed for the L [Formula: see text] H transition (Diamond, Liang, Carreras, Terry 1994 Phys. Rev. Lett. 72, 2565-2568 (doi:10.1103/PhysRevLett.72.2565); Kim & Diamond 2003 Phys. Rev. Lett. 90, 185006 (doi:10.1103/PhysRevLett.90.185006)). We describe the enhanced particle transport events that occur following shear layer collapse. Open problems and future directions are discussed. This article is part of a discussion meeting issue 'H-mode transition and pedestal studies in fusion plasmas'.

Impact of Magnetic Field Configuration on Heat Transport in Stellarators and Heliotrons

We assess the magnetic field configuration in modern fusion devices by comparing experiments with the same heating power, between a stellarator and a heliotron. The key role of turbulence is evident in the optimized stellarator, while neoclassical processes largely determine the transport in the heliotron device. Gyrokinetic simulations elucidate the underlying mechanisms promoting stronger ion scale turbulence in the stellarator. Similar plasma performances in these experiments suggests that neoclassical and turbulent transport should both be optimized in next step reactor designs.

Vanishing neoclassical viscosity and physics of the shear layer in stellarators

The drift kinetic equation is solved for low density TJ-II plasmas employing slowly varying, time-dependent profiles. This allows us to simulate density ramp-up experiments and describe from first principles the formation and physics of the radial electric field shear layer. The main features of the transition are perfectly captured by the calculation, and good quantitative agreement is also found. The results presented here, that should be valid for other nonquasisymmetric stellarators, provide a fundamental explanation for a wealth of experimental observations connected to the shear layer emergence in TJ-II. The key quantity is the neoclassical viscosity, which is shown to go smoothly to zero when the critical density is approached from below. This makes it possible for turbulence-related phenomena, and particularly zonal flows, to arise in the neighborhood of the transition.

Energetic-particle-induced geodesic acoustic mode

A new energetic particle-induced geodesic acoustic mode (EGAM) is shown to exist. The mode frequency and mode structure are determined nonperturbatively by energetic particle kinetic effects. In particular the EGAM frequency is found to be substantially lower than the standard GAM frequency. The radial mode width is determined by the energetic particle drift orbit width and can be fairly large for high energetic particle pressure and large safety factor. These results are consistent with the recent experimental observation of the beam-driven n=0 mode in DIII-D.

Intrinsic ambipolarity and rotation in stellarators

It is shown that collisional plasma transport is intrinsically ambipolar only in quasiaxisymmetric or quasihelically symmetric magnetic configurations. Only in such fields can the plasma rotate freely, and then only in the direction of quasisymmetry. In a non-quasi-symmetric magnetic field, the average radial electric field is determined by parallel viscosity, which in turn is usually governed by collisional processes. Locally, the radial electric field may be affected by turbulent Reynolds stress producing zonal flows, but on a radial average taken over several ion gyroradii, it is determined by parallel viscosity, at least if the turbulence is electrostatic and obeys the conventional gyrokinetic orderings. This differs from the situation in a tokamak, where there is no flow damping by parallel viscosity in the symmetry direction and the turbulent Reynolds stress may affect the global radial electric field.

Neoclassical simulation of tokamak plasmas using the continuum gyrokinetic code TEMPEST

We present gyrokinetic neoclassical simulations of tokamak plasmas with a self-consistent electric field using a fully nonlinear (full- f ) continuum code TEMPEST in a circular geometry. A set of gyrokinetic equations are discretized on a five-dimensional computational grid in phase space. The present implementation is a method of lines approach where the phase-space derivatives are discretized with finite differences, and implicit backward differencing formulas are used to advance the system in time. The fully nonlinear Boltzmann model is used for electrons. The neoclassical electric field is obtained by solving the gyrokinetic Poisson equation with self-consistent poloidal variation. With a four-dimensional (psi,theta,micro) version of the TEMPEST code, we compute the radial particle and heat fluxes, the geodesic-acoustic mode, and the development of the neoclassical electric field, which we compare with neoclassical theory using a Lorentz collision model. The present work provides a numerical scheme for self-consistently studying important dynamical aspects of neoclassical transport and electric field in toroidal magnetic fusion devices.

Reduction of turbulent transport with zonal flows enhanced in helical systems

Gyrokinetic Vlasov simulations of the ion temperature gradient turbulence are performed in order to investigate effects of helical magnetic configurations on turbulent transport and zonal flows. The obtained results confirm the theoretical prediction that helical configurations optimized for reducing neoclassical ripple transport can simultaneously reduce the turbulent transport with enhancing zonal-flow generation. Stationary zonal-flow structures accompanied with transport reduction are clearly identified by the simulation for the neoclassically optimized helical geometry. The generation of the stationary zonal flow explains a physical mechanism for causing the confinement improvement observed in the inward-shifted plasma in the Large Helical Device [O. Motojima, Nucl. Fusion 43, 1674 (2003)].