F. Jenko

Electron Temperature Gradient Driven 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.

Chapter 2: Plasma confinement and transport

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.

Gyrokinetic simulation of collisionless trapped-electron mode turbulence

Some basic properties of collisionless, trapped-electron mode turbulence in tokamaks are investigated by means of massively parallel gyrokinetic Vlasov simulations. In particular, the spatial structure and wave number spectra of various fluctuating plasma quantities are presented and discussed. An analysis of several cross phase relations supports the view that the transport-dominating scales may be interpreted in terms of remnants of linear modes. In a few test cases, zonal flows are artificially suppressed, demonstrating that their influence on the transport level is small. Finally, the dependence of the latter on several plasma parameters is studied.

Critical Gradient Formula for Toroidal Electron Temperature Gradient Modes

Under certain conditions, the electron heat transport induced by electron temperature gradient (ETG) streamers is sufficiently large and sensitive with respect to the normalized electron temperature gradient to represent a possible cause for electron temperature profile consistency (“stiffness”). Here, linear gyrokinetic simulations of toroidal ETG modes in tokamak core and edge plasmas are presented. An algebraic formula for the threshold of the linear instability is derived from the numerical solutions of the linear gyrokinetic equations which recovers previous analytical results in the appropriate limits.

The global version of the gyrokinetic turbulence code GENE

The understanding and prediction of transport due to plasma microturbulence is a key open problem in modern plasma physics, and a grand challenge for fusion energy research. Ab initio simulations of such small-scale, low-frequency turbulence are to be based on the gyrokinetic equations, a set of nonlinear integro-differential equations in reduced (five-dimensional) phase space. In the present paper, the extension of the well-established and widely used gyrokinetic code GENE [F. Jenko, W. Dorland, M. Kotschenreuther, B.N. Rogers, Electron temperature gradient driven turbulence, Phys. Plasmas 7 (2000) 1904-1910] from a radially local to a radially global (nonlocal) version is described. The necessary modifications of both the basic equations and the employed numerical methods are detailed, including, e.g., the change from spectral methods to finite difference and interpolation techniques in the radial direction and the implementation of sources and sinks. In addition, code verification studies and benchmarks are presented

Characterizing electron temperature gradient turbulence via numerical simulation

Numerical simulations of electron temperature gradient (ETG) turbulence are presented that characterize the ETG fluctuation spectrum, establish limits to the validity of the adiabatic ion model often employed in studying ETG turbulence, and support the tentative conclusion that plasma-operating regimes exist in which ETG turbulence produces sufficient electron heat transport to be experimentally relevant. We resolve prior controversies regarding simulation techniques and convergence by benchmarking simulations of ETG turbulence from four microturbulence codes, demonstrating agreement on the electron heat flux, correlation functions, fluctuation intensity, and rms flow shear at fixed simulation cross section and resolution in the plane perpendicular to the magnetic field. Excellent convergence of both continuum and particle-in-cell codes with time step and velocity-space resolution is demonstrated, while numerical issues relating to perpendicular (to the magnetic field) simulation dimensions and resolution are discussed. A parameter scan in the magnetic shear, s, demonstrates that the adiabatic ion model is valid at small values of s (s < 0.4 for the parameters used in this scan) but breaks down at higher magnetic shear. A proper treatment employing gyrokinetic ions reveals a steady increase in the electron heat transport with increasing magnetic shear, reaching electron heat transport rates consistent with analyses of experimental tokamak discharges. (c) 2006 American Institute of Physics.

Gyrokinetic turbulence simulations at high plasma beta

Electromagnetic gyrokinetic turbulence simulations employing Cyclone Base Case parameters are presented for β values up to and beyond the kinetic ballooning threshold. The β scaling of the turbulent transport is found to be linked to a complex interplay of linear and nonlinear effects. Linear investigation of the kinetic ballooning mode is performed in detail, while nonlinearly, it is found to dominate the turbulence only in a fairly narrow range of β values just below the respective ideal limit. The magnetic transport scales like β2 and is well described by a Rechester–Rosenbluth-type ansatz.Electromagnetic gyrokinetic turbulence simulations employing Cyclone Base Case parameters are presented for β values up to and beyond the kinetic ballooning threshold. The β scaling of the turbulent transport is found to be linked to a complex interplay of linear and nonlinear effects. Linear investigation of the kinetic ballooning mode is performed in detail, while nonlinearly, it is found to dominate the turbulence only in a fairly narrow range of β values just below the respective ideal limit. The magnetic transport scales like β2 and is well described by a Rechester–Rosenbluth-type ansatz.

Clarifications to the limitations of the s-α equilibrium model for gyrokinetic computations of turbulence

In the context of gyrokinetic flux-tube simulations of microturbulence in magnetized toroidal plasmas, different treatments of the magnetic equilibrium are examined. Considering the Cyclone DIII-D base case parameter set [Dimits et al., Phys. Plasmas 7, 969 (2000)], significant differences in the linear growth rates, the linear and nonlinear critical temperature gradients, and the nonlinear ion heat diffusivities are observed between results obtained using either an

Nonlinear stabilization of tokamak microturbulence by fast ions.

s-\alpha

Interaction of energetic particles with large and small scale instabilities

or an MHD equilibrium. Similar disagreements have been reported previously [Redd et al., Phys. Plasmas 6, 1162 (1999)]. In this paper it is shown that these differences result primarily from the approximation made in the standard implementation of the