1,589 research outputs found
Optimisation of confinement in a fusion reactor using a nonlinear turbulence model
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
first-principles nonlinear simulations which accurately capture the turbulent
cascade and large-scale zonal flows. The simulations utilise a novel approach,
with gyrofluid treatment of the small-scale drift waves and gyrokinetic
treatment of the large-scale zonal flows. A simple near-circular 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
two-fold increase in the plasma power per unit volume is achieved by moving to
higher elongation and strongly negative triangularity.Comment: 32 pages, 8 figures, accepted to JP
Phase mixing vs. nonlinear advection in drift-kinetic plasma turbulence
A scaling theory of long-wavelength electrostatic turbulence in a magnetised,
weakly collisional plasma (e.g., ITG turbulence) is proposed, with account
taken both of the nonlinear advection of the perturbed particle distribution by
fluctuating ExB flows and of its phase mixing, which is caused by the streaming
of the particles along the mean magnetic field and, in a linear problem, would
lead to Landau damping. It is found that it is possible to construct a
consistent theory in which very little free energy leaks into high velocity
moments of the distribution function, rendering the turbulent cascade in the
energetically relevant part of the wave-number space essentially fluid-like.
The velocity-space spectra of free energy expressed in terms of Hermite-moment
orders are steep power laws and so the free-energy content of the phase space
does not diverge at infinitesimal collisionality (while it does for a linear
problem); collisional heating due to long-wavelength perturbations vanishes in
this limit (also in contrast with the linear problem, in which it occurs at the
finite rate equal to the Landau-damping rate). The ability of the free energy
to stay in the low velocity moments of the distribution function is facilitated
by the "anti-phase-mixing" effect, whose presence in the nonlinear system is
due to the stochastic version of the plasma echo (the advecting velocity
couples the phase-mixing and anti-phase-mixing perturbations). The partitioning
of the wave-number space between the (energetically dominant) region where this
is the case and the region where linear phase mixing wins its competition with
nonlinear advection is governed by the "critical balance" between linear and
nonlinear timescales (which for high Hermite moments splits into two
thresholds, one demarcating the wave-number region where phase mixing
predominates, the other where plasma echo does).Comment: 45 pages (single-column), 3 figures, replaced with version published
in JP
Kinetic Simulations of Magnetized Turbulence in Astrophysical Plasmas
This letter presents the first ab initio, fully electromagnetic, kinetic
simulations of magnetized turbulence in a homogeneous, weakly collisional
plasma at the scale of the ion Larmor radius (ion gyroscale). Magnetic and
electric-field energy spectra show a break at the ion gyroscale; the spectral
slopes are consistent with scaling predictions for critically balanced
turbulence of Alfven waves above the ion gyroscale (spectral index -5/3) and of
kinetic Alfven waves below the ion gyroscale (spectral indices of -7/3 for
magnetic and -1/3 for electric fluctuations). This behavior is also
qualitatively consistent with in situ measurements of turbulence in the solar
wind. Our findings support the hypothesis that the frequencies of turbulent
fluctuations in the solar wind remain well below the ion cyclotron frequency
both above and below the ion gyroscale.Comment: 4 pages, 3 figures, submitted to Physical Review Letter
Rotation and Neoclassical Ripple Transport in ITER
Neoclassical transport in the presence of non-axisymmetric magnetic fields
causes a toroidal torque known as neoclassical toroidal viscosity (NTV). The
toroidal symmetry of ITER will be broken by the finite number of toroidal field
coils and by test blanket modules (TBMs). The addition of ferritic inserts
(FIs) will decrease the magnitude of the toroidal field ripple. 3D magnetic
equilibria with toroidal field ripple and ferromagnetic structures are
calculated for an ITER steady-state scenario using the Variational Moments
Equilibrium Code (VMEC). Neoclassical transport quantities in the presence of
these error fields are calculated using the Stellarator Fokker-Planck Iterative
Neoclassical Conservative Solver (SFINCS). These calculations fully account for
, flux surface shaping, multiple species, magnitude of ripple, and
collisionality rather than applying approximate analytic NTV formulae. As NTV
is a complicated nonlinear function of , we study its behavior over a
plausible range of . We estimate the toroidal flow, and hence , using
a semi-analytic turbulent intrinsic rotation model and NUBEAM calculations of
neutral beam torque. The NTV from the ripple dominates
that from lower perturbations of the TBMs. With the inclusion of FIs, the
magnitude of NTV torque is reduced by about 75% near the edge. We present
comparisons of several models of tangential magnetic drifts, finding
appreciable differences only for superbanana-plateau transport at small .
We find the scaling of calculated NTV torque with ripple magnitude to indicate
that ripple-trapping may be a significant mechanism for NTV in ITER. The
computed NTV torque without ferritic components is comparable in magnitude to
the NBI and intrinsic turbulent torques and will likely damp rotation, but the
NTV torque is significantly reduced by the planned ferritic inserts
Gyrokinetic Simulations of Solar Wind Turbulence from Ion to Electron Scales
The first three-dimensional, nonlinear gyrokinetic simulation of plasma
turbulence resolving scales from the ion to electron gyroradius with a
realistic mass ratio is presented, where all damping is provided by resolved
physical mechanisms. The resulting energy spectra are quantitatively consistent
with a magnetic power spectrum scaling of as observed in \emph{in
situ} spacecraft measurements of the "dissipation range" of solar wind
turbulence. Despite the strongly nonlinear nature of the turbulence, the linear
kinetic \Alfven wave mode quantitatively describes the polarization of the
turbulent fluctuations. The collisional ion heating is measured at
sub-ion-Larmor radius scales, which provides the first evidence of the ion
entropy cascade in an electromagnetic turbulence simulation.Comment: 4 pages, 2 figures, submitted to Phys. Rev. Let
Transition to subcritical turbulence in a tokamak plasma
Tokamak turbulence, driven by the ion-temperature gradient and occurring in
the presence of flow shear, is investigated by means of local, ion-scale,
electrostatic gyrokinetic simulations (with both kinetic ions and electrons) of
the conditions in the outer core of the Mega-Ampere Spherical Tokamak (MAST). A
parameter scan in the local values of the ion-temperature gradient and flow
shear is performed. It is demonstrated that the experimentally observed state
is near the stability threshold and that this stability threshold is nonlinear:
sheared turbulence is subcritical, i.e. the system is formally stable to small
perturbations, but, given a large enough initial perturbation, it transitions
to a turbulent state. A scenario for such a transition is proposed and
supported by numerical results: close to threshold, the nonlinear saturated
state and the associated anomalous heat transport are dominated by long-lived
coherent structures, which drift across the domain, have finite amplitudes, but
are not volume filling; as the system is taken away from the threshold into the
more unstable regime, the number of these structures increases until they
overlap and a more conventional chaotic state emerges. Whereas this appears to
represent a new scenario for transition to turbulence in tokamak plasmas, it is
reminiscent of the behaviour of other subcritically turbulent systems, e.g.
pipe flows and Keplerian magnetorotational accretion flows.Comment: 16 pages, 5 figures, accepted to Journal of Plasma Physic
Multiscale Gyrokinetics for Rotating Tokamak Plasmas: Fluctuations, Transport and Energy Flows
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 order-by-order 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 transport-timescale 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 Grad-Shafranov
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.
Large-scale deviations of the distribution function from a Maxwellian are given
by neoclassical theory. The fluctuations are determined by the high-flow
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 low-Mach-number limit, to a concise set of
equations suitable for numerical implementation.Comment: 113 pages, 3 figure
Zero-Turbulence Manifold in a Toroidal Plasma
Sheared toroidal flows can cause bifurcations to zero-turbulent-transport
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 zero-turbulence manifold is mapped
out, in the zero-magnetic-shear 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.Comment: 5 Pages, 4 Figures, Submitted to PR
- …