F. J. Casson

Overview of toroidal momentum transport

Toroidal momentum transport mechanisms are reviewed and put in a broader perspective. The generation of a finite momentum flux is closely related to the breaking of symmetry (parity) along the field. The symmetry argument allows for the systematic identification of possible transport mechanisms. Those that appear to lowest order in the normalized Larmor radius (the diagonal part, Coriolis pinch, E x B shearing, particle flux, and up-down asymmetric equilibria) are reasonably well understood. At higher order, expected to be of importance in the plasma edge, the theory is still under development.

Influence of the centrifugal force and parallel dynamics on the toroidal momentum transport due to small scale turbulence in a tokamak

The paper derives the gyro-kinetic equation in the comoving frame of a toroidally rotating plasma, including both the Coriolis drift effect [A. G. Peeters et al., Phys. Rev. Lett. 98, 265003 (2007)] as well as the centrifugal force. The relation with the laboratory frame is discussed. A low field side gyro-fluid model is derived from the gyro-kinetic equation and applied to the description of parallel momentum transport. The model includes the effects of the Coriolis and centrifugal force as well as the parallel dynamics. The latter physics effect allows for a consistent description of both the Coriolis drift effect as well as the ExB shear effect [R. R. Dominguez and G. M. Staebler, Phys. Fluids B 5, 3876 (1993)] on the momentum transport. Strong plasma rotation as well as parallel dynamics reduce the Coriolis (inward) pinch of momentum and can lead to a sign reversal generating an outward pinch velocity. Also, the ExB shear effect is, in a similar manner, reduced by the parallel dynamics and stronger rotation.

Off-diagonal particle and toroidal momentum transport: a survey of experimental, theoretical and modelling aspects

In tokamaks, turbulent particle and toroidal momentum transport are both characterized by the presence of off-diagonal contributions which play an essential role in establishing the profile shapes of the density and the toroidal rotation under most conditions. In this paper similarities and differences between the two turbulent transport channels are pointed out and, thereby, interesting physical aspects which connect the two channels are identified. The main contributions to off-diagonal particle and toroidal momentum transport are reviewed by means of a rather simplified description, which aims at providing, when possible, a direct connection between theoretical, modelling and experimental research.

Tungsten transport in JET H-mode plasmas in hybrid scenario, experimental observations and modelling

The behaviour of tungsten in the core of hybrid scenario plasmas in JET with the ITER-like wall is analysed and modelled with a combination of neoclassical and gyrokinetic codes. In these discharges, good confinement conditions can be maintained only for the first 2?3?s of the high power phase. Later W accumulation is regularly observed, often accompanied by the onset of magneto-hydrodynamical activity, in particular neoclassical tearing modes (NTMs), both of which have detrimental effects on the global energy confinement. The dynamics of the accumulation process is examined, taking into consideration the concurrent evolution of the background plasma profiles, and the possible onset of NTMs. Two time slices of a representative discharge, before and during the accumulation process, are analysed with two independent methods, in order to reconstruct the W density distribution over the poloidal cross-section. The same time slices are modelled, computing both neoclassical and turbulent transport components and consistently including the impact of centrifugal effects, which can be significant in these plasmas, and strongly enhance W neoclassical transport. The modelling closely reproduces the observations and identifies inward neoclassical convection due to the density peaking of the bulk plasma in the central region as the main cause of the accumulation. The change in W neoclassical convection is directly produced by the transient behaviour of the main plasma density profile, which is hollow in the central region in the initial part of the high power phase of the discharge, but which develops a significant density peaking very close to the magnetic axis in the later phase. The analysis of a large set of discharges provides clear indications that this effect is generic in this scenario. The unfavourable impact of the onset of NTMs on the W behaviour, observed in several discharges, is suggested to be a consequence of a detrimental combination of the effects of neoclassical transport and of the appearance of an island.

Anomalous parallel momentum transport due to E×B flow shear in a tokamak plasma

Nondiffusive anomalous momentum transport in toroidal plasmas occurs through symmetry breaking mechanisms. In this paper the contribution of sheared E X B flows to parallel momentum transport [R. R. Dominguez and G. M. Staebler, Phys Fluids B 5, 3876 (1993)] is investigated with nonlinear gyrokinetic simulations in toroidal geometry. The background perpendicular shear is treated independently from the parallel velocity shear to isolate a nondiffusive, nonpinch contribution to the parallel momentum flux. It is found that the size of the term depends strongly on the magnetic shear, with the sign reversing for negative magnetic shear. Perpendicular shear flows are responsible for both symmetry breaking and suppression of turbulence, resulting in a shearing rate at which there is a maximum contribution to the momentum transport. The E X B momentum transport is shown to be quenched by increasing flow shear more strongly than the standard linear quench rule for turbulent heat diffusivity

Gyrokinetic simulations including the centrifugal force in a rotating tokamak plasma

Tokamak experiments operate with a rotating plasma, with toroidal velocity which can be driven externally but can also arise spontaneously. In the frame that corotates with the plasma, the effects of the centrifugal force are felt through a centrifugal drift and an enhanced mirror force [Peeters et al., Phys. Plasmas 16, 042310 (2009)]. These inertial terms become important in the case of strong rotation, as is common in spherical devices, and are also important for heavy impurity ions even at small toroidal velocities. In this work, the first gyrokinetic simulations including the centrifugal force in a strongly rotating plasma are presented. The enhanced mirror force redistributes density over a flux surface and modifies the trapping condition, destabilizing trapped electron modes. At intermediate scales this can result in promotion of the trapped electron mode over the ion temperature gradient (ITG) mode as the dominant instability, which under marginal conditions could result in an enhanced electron heat flux. The centrifugal drift acts to damp the residual zonal flow of the geoacoustic mode, while its frequency is increased. For nonlinear ITG dominated turbulence, increased trapped electron drive and reduced zonal flow lead to an increase in ion heat diffusivity if the increased rotation is not accompanied by rotational shear stabilization. An increased fraction of slow trapped electrons enhances the convective particle pinch, leading to an increase in the steady state density gradient with strong rotation. Linear ITG mode results show an increased pinch of heavy trace impurities due to their strong centrifugal trapping.

Gyrokinetic simulations of spherical tokamaks

This paper reviews transport and confinement in spherical tokamaks (STs) and our current physics understanding of this that is partly based on gyrokinetic simulations. Equilibrium flow shear plays an important role, and we show how this is consistently included in the gyrokinetic framework for flows that greatly exceed the diamagnetic velocity. The key geometry factors that influence the effectiveness of turbulence suppression by flow shear are discussed, and we show that toroidal equilibrium flow shear can sometimes entirely suppress ion scale turbulence in todays STs. Advanced nonlinear simulations of electron temperature gradient (ETG) driven turbulence, including kinetic ion physics, collisions and equilibrium flow shear, support the model that ETG turbulence can explain electron heat transport in many ST discharges.

The nonlinear coupling between gyroradius scale turbulence and mesoscale magnetic islands in fusion plasmas

The interaction between small scale turbulence (of the order of the ion Larmor radius) and mesoscale magnetic islands is investigated within the gyrokinetic framework. Turbulence, driven by background temperature and density gradients, over nonlinear mode coupling, pumps energy into long wavelength modes, and can result in an electrostatic vortex mode that coincides with the magnetic island. The strength of the vortex is strongly enhanced by the modified plasma flow response connected with the change in topology, and the transport it generates can compete with the parallel motion along the perturbed magnetic field. Despite the stabilizing effect of sheared plasma flows in and around the island, the net effect of the island is a degradation of the confinement. When density and temperature gradients inside the island are below the threshold for turbulence generation, turbulent fluctuations still persist through turbulence convection and spreading. The latter mechanisms then generate a finite transport flux ...

Intrinsic rotation driven by the electrostatic turbulence in up-down asymmetric toroidal plasmas

The transport of parallel momentum by small scale fluctuations is intrinsically linked to symmetry breaking in the direction of the magnetic field. In tokamaks, an up-down asymmetry in the equilibrium proves to be an efficient parallel symmetry breaking mechanism leading to the generation of a net radial flux of parallel momentum by the electrostatic turbulence [Y. Camenen et al., Phys. Rev. Lett. 102, 125001 (2009)]. This flux is neither proportional to the toroidal rotation nor to its gradient and arises from an incomplete cancellation of the local contributions to the parallel momentum flux under the flux surface average. The flux of parallel momentum then depends on the asymmetry of the curvature drift and on the extension of the fluctuations around the low field side midplane. In this paper, the mechanisms underlying the generation of the flux of parallel momentum are highlighted and the main dependences on plasma parameters investigated using linear gyrokinetic simulations.

Theoretical description of heavy impurity transport and its application to the modelling of tungsten in JET and ASDEX upgrade

The effects of poloidal asymmetries and heated minority species are shown to be necessary to accurately describe heavy impurity transport in present experiments in JET and ASDEX Upgrade. Plasma rotation, or any small background electrostatic field in the plasma, such as that generated by anisotropic external heating can generate strong poloidal density variation of heavy impurities. These asymmetries have recently been added to numerical tools describing both neoclassical and turbulent transport and can increase neoclassical tungsten transport by an order of magnitude. Modelling predictions of the steady-state two-dimensional tungsten impurity distribution are compared with tomography from soft x-ray diagnostics. The modelling identifies neoclassical transport enhanced by poloidal asymmetries as the dominant mechanism responsible for tungsten accumulation in the central core of the plasma. Depending on the bulk plasma profiles, turbulent diffusion and neoclassical temperature screening can prevent accumulation. Externally heated minority species can significantly enhance temperature screening in ICRH plasmas.