
The confinement of heat in the core of a magnetic fusion reactor is optimised
using a multidimensional optimisation algorithm. For the first time in such a
study, the loss of heat due to turbulence is modelled at every stage using
firstprinciples nonlinear simulations which accurately capture the turbulent
cascade and largescale zonal flows. The simulations utilise a novel approach,
with gyrofluid treatment of the smallscale drift waves and gyrokinetic
treatment of the largescale zonal flows. A simple nearcircular equilibrium
with standard parameters is chosen as the initial condition. The figure of
merit, fusion power per unit volume, is calculated, and then two control
parameters, the elongation and triangularity of the outer flux surface, are
varied, with the algorithm seeking to optimise the chosen figure of merit. A
twofold increase in the plasma power per unit volume is achieved by moving to
higher elongation and strongly negative triangularity.

This paper describes a new code for simulating astrophysical plasmas that
solves a hybrid model composed of gyrokinetic ions (GKI) and an isothermal
electron fluid (ITEF) [A. Schekochihin et al., Astrophys. J. Suppl.
\textbf{182}, 310 (2009)]. This model captures ion kinetic effects that are
important near the ion gyroradius scale while electron kinetic effects are
ordered out by an electronion mass ratio expansion. The code is developed by
incorporating the ITEF approximation into {\tt AstroGK}, an Eulerian $\delta f$
gyrokinetics code specialized to a slab geometry [R. Numata et al., J. Compute.
Pays. \textbf{229}, 9347 (2010)]. The new code treats the linear terms in the
ITEF equations implicitly while the nonlinear terms are treated explicitly. We
show linear and nonlinear benchmark tests to prove the validity and
applicability of the simulation code. Since the fast electron timescale is
eliminated by the mass ratio expansion, the CourantFriedrichsLewy condition
is much less restrictive than in full gyrokinetic codes; the present hybrid
code runs $\sim 2\sqrt{m_\mathrm{i}/m_\mathrm{e}} \sim 100$ times faster than
{\tt AstroGK}\ with a single ion species and kinetic electrons where
$m_\mathrm{i}/m_\mathrm{e}$ is the ionelectron mass ratio. The improvement of
the computational time makes it feasible to execute ion scale gyrokinetic
simulations with a high velocity space resolution and to run multiple
simulations to determine the dependence of turbulent dynamics on parameters
such as electronion temperature ratio and plasma beta.

Observational evidence in space and astrophysical plasmas with long
collisional mean free path suggests that more massive charged particles may be
preferentially heated. One possible mechanism for this is the turbulent cascade
of energy from injection to dissipation scales, where the energy is converted
to heat. Here we consider a simple system consisting of a magnetized plasma
slab of electrons and a single ion species with a crossfield density gradient.
We show that such a system is subject to an electron drift wave instability,
known as the universal instability, which is stabilized only when the electron
and ion thermal speeds are equal. For unequal thermal speeds, we find that the
instability gives rise to turbulent energy exchange between ions and electrons
that acts to equalize the thermal speeds. Consequently, this turbulent heating
tends to equalize the component temperatures of pair plasmas and to heat ions
to much higher temperatures than electrons for conventional massratio plasmas.

Tokamak plasmas rotate even without external injection of momentum. A Doppler
backscattering system installed at MAST has allowed this intrinsic rotation to
be studied in Ohmic Lmode and Hmode plasmas, including the first observation
of intrinsic rotation reversals in a spherical tokamak. Experimental results
are compared to a novel 1D model, which captures the collisionality dependence
of the radial transport of toroidal angular momentum due to the effect of
neoclassical flows on turbulent fluctuations. The model is able to accurately
reproduce the change in sign of core toroidal rotation, using experimental
density and temperature profiles from shots with rotation reversals as inputs
and no free parameters fit to experimental data.

The global ideal kink equation, for cylindrical geometry and zero beta, is
simplified in the high poloidal mode number limit and used to determine the
tearing stability parameter, $\Delta^\prime$. In the presence of a steep
monotonic current gradient, $\Delta^\prime$ becomes a function of a parameter,
$\sigma_0$, characterising the ratio of the maximum current gradient to
magnetic shear, and $x_s$, characterising the separation of the resonant
surface from the maximum of the current gradient. In equilibria containing a
current "spike", so that there is a nonmonotonic current profile,
$\Delta^\prime$ also depends on two parameters: $\kappa$, related to the ratio
of the curvature of the current density at its maximum to the magnetic shear,
and $x_s$, which now represents the separation of the resonance from the point
of maximum current density. The relation of our results to earlier studies of
tearing modes and to recent gyrokinetic calculations of current driven
instabilities, is discussed, together with potential implications for the
stability of the tokamak pedestal.

This paper presents a complete theoretical framework for plasma turbulence
and transport in tokamak plasmas. The fundamental scale separations present in
plasma turbulence are codified as an asymptotic expansion in the ratio of the
gyroradius to the equilibrium scale length. Proceeding orderbyorder in this
expansion, a framework for plasma turbulence is developed. It comprises an
instantaneous equilibrium, the fluctuations driven by gradients in the
equilibrium quantities, and the transporttimescale evolution of mean profiles
of these quantities driven by the fluctuations. The equilibrium distribution
functions are local Maxwellians with each flux surface rotating toroidally as a
rigid body. The magnetic equillibrium is obtained from the GradShafranov
equation for a rotating plasma and the slow (resistive) evolution of the
magnetic field is given by an evolution equation for the safety factor q.
Largescale deviations of the distribution function from a Maxwellian are given
by neoclassical theory. The fluctuations are determined by the highflow
gyrokinetic equation, from which we derive the governing principle for
gyrokinetic turbulence in tokamaks: the conservation and local cascade of free
energy. Transport equations for the evolution of the mean density, temperature
and flow velocity profiles are derived. These transport equations show how the
neoclassical corrections and the fluctuations act back upon the mean profiles
through fluxes and heating. The energy and entropy conservation laws for the
mean profiles are derived. Total energy is conserved and there is no net
turbulent heating. Entropy is produced by the action of fluxes flattening
gradients, Ohmic heating, and the equilibration of mean temperatures. Finally,
this framework is condensed, in the lowMachnumber limit, to a concise set of
equations suitable for numerical implementation.

The effect of small deviations from a Maxwellian equilibrium on turbulent
momentum transport in tokamak plasmas is considered. These nonMaxwellian
features, arising from diamagnetic effects, introduce a strong dependence of
the radial flux of cocurrent toroidal angular momentum on collisionality: As
the plasma goes from nearly collisionless to weakly collisional, the flux
reverses direction from radially inward to outward. This indicates a
collisionalitydependent transition from peaked to hollow rotation profiles,
consistent with experimental observations of intrinsic rotation.

In magnetized plasmas, a turbulent cascade occurs in phase space at scales
smaller than the thermal Larmor radius ("subLarmor scales") [Phys. Rev. Lett.
103, 015003 (2009)]. When the turbulence is restricted to two spatial
dimensions perpendicular to the background magnetic field, two independent
cascades may take place simultaneously because of the presence of two
collisionless invariants. In the present work, freely decaying turbulence of
twodimensional electrostatic gyrokinetics is investigated by means of
phenomenological theory and direct numerical simulations. A dual cascade
(forward and inverse cascades) is observed in velocity space as well as in
position space, which we diagnose by means of nonlinear transfer functions for
the collisionless invariants. We find that the turbulence tends to a
timeasymptotic state, dominated by a single scale that grows in time. A theory
of this asymptotic state is derived in the form of decay laws. Each case that
we study falls into one of three regimes (weakly collisional, marginal, and
strongly collisional), determined by a dimensionless number D*, a quantity
analogous to the Reynolds number. The marginal state is marked by a critical
number D* = D0 that is preserved in time. Turbulence initialized above this
value become increasingly inertial in time, evolving toward larger and larger
D*; turbulence initialized below D0 become more and more collisional, decaying
to progressively smaller D*.

Beam Emission Spectroscopy (BES) measurements of ionscale density
fluctuations in the MAST tokamak are used to show that the turbulence
correlation time, the drift time associated with ion temperature or density
gradients, the particle (ion) streaming time along the magnetic field and the
magnetic drift time are consistently comparable, suggesting a "critically
balanced" turbulence determined by the local equilibrium. The resulting
scalings of the poloidal and radial correlation lengths are derived and tested.
The nonlinear time inferred from the density fluctuations is longer than the
other times; its ratio to the correlation time scales as
$\nu_{*i}^{0.8\pm0.1}$, where $\nu_{*i}=$ ion collision rate/streaming rate.
This is consistent with turbulent decorrelation being controlled by a zonal
component, invisible to the BES, with an amplitude exceeding the drift waves'
by $\sim \nu_{*i}^{0.8}$.

Scaling laws for the transport and heating of trace heavy ions in
lowfrequency, magnetized plasma turbulence are derived and compared with
direct numerical simulations. The predicted dependences of turbulent fluxes and
heating on ion charge and mass number are found to agree with numerical results
for both stationary and differentially rotating plasmas. Heavy ion momentum
transport is found to increase with mass, and heavy ions are found to be
preferentially heated, implying a massdependent ion temperature for very
weakly collisional plasmas and for partiallyionized heavy ions in strongly
rotating plasmas.

Sheared toroidal flows can cause bifurcations to zeroturbulenttransport
states in tokamak plasmas. The maximum temperature gradients that can be
reached are limited by subcritical turbulence driven by the parallel velocity
gradient. Here it is shown that q/\epsilon (magnetic field pitch/inverse aspect
ratio) is a critical control parameter for sheared tokamak turbulence. By
reducing q/\epsilon, far higher temperature gradients can be achieved without
triggering turbulence, in some instances comparable to those found
experimentally in transport barriers. The zeroturbulence manifold is mapped
out, in the zeromagneticshear limit, over the parameter space (\gamma_E,
q/\epsilon, R/L_T), where \gamma_E is the perpendicular flow shear and R/L_T is
the normalised inverse temperature gradient scale. The extent to which it can
be constructed from linear theory is discussed.

Firstprinciples numerical simulations are used to describe a transport
bifurcation in a differentially rotating tokamak plasma. Such a bifurcation is
more probable in a region of zero magnetic shear than one of finite magnetic
shear because in the former case the component of the sheared toroidal flow
that is perpendicular to the magnetic field has the strongest suppressing
effect on the turbulence. In the zeromagneticshear regime, there are no
growing linear eigenmodes at any finite value of flow shear. However,
subcritical turbulence can be sustained, owing to the transient growth of modes
driven by the ion temperature gradient (ITG) and the parallel velocity gradient
(PVG). Nonetheless, in a parameter space containing a wide range of temperature
gradients and velocity shears, there is a sizeable window where all turbulence
is suppressed. Combined with the relatively low transport of momentum by
collisional (neoclassical) mechanisms, this produces the conditions for a
bifurcation from low to high temperature and velocity gradients. The path of
this bifurcation is mapped out using interpolation from a large number of
simulations. Numerical simulations are also used to construct a parametric
model which accurately describes the combined effect of the temperature
gradient and the flow gradient over a wide range of their values. Using this
parametric model, it is shown that in this reducedtransport state, heat is
transported almost neoclassically, while momentum transport is dominated by
subcritical PVG turbulence. It is further shown that for any given input of
torque, there is an optimum input of heat which maximises the temperature
gradient. The parametric model describes both the behaviour of the subcritical
turbulence and the complicated effect of the flow shear on the transport
stiffness. It may prove useful for transport modelling of tokamaks with sheared
flows.

Scaling laws for ion temperature gradient driven turbulence in magnetized
toroidal plasmas are derived and compared with direct numerical simulations.
Predicted dependences of turbulence fluctuation amplitudes, spatial scales, and
resulting heat fluxes on temperature gradient and magnetic field line pitch are
found to agree with numerical results in both the driving and inertial ranges.
Evidence is provided to support the critical balance conjecture that parallel
streaming and nonlinear perpendicular decorrelation times are comparable at all
spatial scales, leading to a scaling relationship between parallel and
perpendicular spatial scales. This indicates that even strongly magnetized
plasma turbulence is intrinsically threedimensional.

Injection and beam dumping is considered for a 16.5 TeV hadron accelerator in
the current LHC tunnel, with an injection energy in the range 1  1.3 TeV. The
present systems are described and the possible upgrade scenarios investigated
for higher beam rigidity. In addition to the required equipment performance,
the machine protection related aspects are explored. The expected constraints
on the machine layout are also given. The technological challenges for the
different equipment subsystems are detailed, and areas where R&D is necessary
are highlighted.

The effect of momentum injection on the temperature gradient in tokamak
plasmas is studied. A plausible scenario for transitions to reduced transport
regimes is proposed. The transition happens when there is sufficient momentum
input so that the velocity shear can suppress or reduce the turbulence.
However, it is possible to drive too much velocity shear and rekindle the
turbulent transport. The optimal level of momentum injection is determined. The
reduction in transport is maximized in the regions of low or zero magnetic
shear.

The effect of flow shear on turbulent transport in tokamaks is studied
numerically in the experimentally relevant limit of zero magnetic shear. It is
found that the plasma is linearly stable for all nonzero flow shear values,
but that subcritical turbulence can be sustained nonlinearly at a wide range of
temperature gradients. Flow shear increases the nonlinear temperature gradient
threshold for turbulence but also increases the sensitivity of the heat flux to
changes in the temperature gradient, except over a small range near the
threshold where the sensitivity is decreased. A bifurcation in the equilibrium
gradients is found: for a given input of heat, it is possible, by varying the
applied torque, to trigger a transition to significantly higher temperature and
flow gradients.

Twodimensional electrostatic turbulence in magnetized weaklycollisional
plasmas exhibits a cascade of entropy in phase space [Phys. Rev. Lett. 103,
015003 (2009)]. At scales smaller than the gyroradius, this cascade is
characterized by the dimensionless ratio D of the collision time to the eddy
turnover time measured at the scale of the thermal Larmor radius. When D >> 1,
a broad spectrum of fluctuations at subLarmor scales is found in both position
and velocity space. The distribution function develops structure as a function
of v_{perp}, the velocity coordinate perpendicular to the local magnetic field.
The cascade shows a localscale nonlinear interaction in both position and
velocity spaces, and Kolmogorov's scaling theory can be extended into phase
space.

Nonlinear gyrokinetic simulations have been conducted to investigate
turbulent transport in tokamak plasmas with rotational shear. At sufficiently
large flow shears, linear instabilities are suppressed, but transiently growing
modes drive subcritical turbulence whose amplitude increases with flow shear.
This leads to a local minimum in the heat flux, indicating an optimal E x B
shear value for plasma confinement. Local maxima in the momentum fluxes are
also observed, allowing for the possibility of bifurcations in the E x B shear.
The sensitive dependence of heat flux on temperature gradient is relaxed for
large flow shear values, with the critical temperature gradient increasing at
lower flow shear values. The turbulent Prandtl number is found to be largely
independent of temperature and flow gradients, with a value close to unity.

Direct coupling between a transport solver and local, nonlinear gyrokinetic
calculations using the multiscale gyrokinetic code TRINITY [M. Barnes, Ph.D.
thesis, arxiv:0901.2868] is described. The coupling of the microscopic and
macroscopic physics is done within the framework of multiscale gyrokinetic
theory, of which we present the assumptions and key results. An assumption of
scale separation in space and time allows for the simulation of turbulence in
small regions of the spacetime grid, which are embedded in a coarse grid on
which the transport equations are implicitly evolved. This leads to a reduction
in computational expense of several orders of magnitude, making
firstprinciples simulations of the full fusion device volume over the
confinement time feasible on current computing resources. Numerical results
from TRINITY simulations are presented and compared with experimental data from
JET and ASDEX Upgrade plasmas.

Electrostatic turbulence in weakly collisional, magnetized plasma can be
interpreted as a cascade of entropy in phase space, which is proposed as a
universal mechanism for dissipation of energy in magnetized plasma turbulence.
When the nonlinear decorrelation time at the scale of the thermal Larmor radius
is shorter than the collision time, a broad spectrum of fluctuations at
subLarmor scales is numerically found in velocity and position space, with
theoretically predicted scalings. The results are important because they
identify what is probably a universal Kolmogorovlike regime for kinetic
turbulence; and because any physical process that produces fluctuations of the
gyrophaseindependent part of the distribution function may, via the entropy
cascade, result in turbulent heating at a rate that increases with the
fluctuation amplitude, but is independent of the collision frequency.

A set of key properties for an ideal dissipation scheme in gyrokinetic
simulations is proposed, and implementation of a model collision operator
satisfying these properties is described. This operator is based on the exact
linearized testparticle collision operator, with approximations to the
fieldparticle terms that preserve conservation laws and an HTheorem. It
includes energy diffusion, pitchangle scattering, and finite Larmor radius
effects corresponding to classical (realspace) diffusion. The numerical
implementation in the continuum gyrokinetic code GS2 is fully implicit and
guarantees exact satisfaction of conservation properties. Numerical results are
presented showing that the correct physics is captured over the entire range of
collisionalities, from the collisionless to the strongly collisional regimes,
without recourse to artificial dissipation.