Prof Michael Barnes

Abstract:

A novel approach to global gyrokinetic simulation is implemented in the flux-tube code stella. This is done by using a subsidiary expansion of the gyrokinetic equation in the perpendicular scale length of the turbulence, originally derived by Parra and Barnes [Plasma Phys. Controlled Fusion, 57 054003, 2015], which allows the use of Fourier basis functions while enabling the effect of radial profile variation to be included in a perturbative way. Radial variation of the magnetic geometry is included by utilizing a global extension of the Grad-Shafranov equation and the Miller equilibrium equations which is obtained through Taylor expansion. Radial boundary conditions that employ multiple flux-tube simulations are also developed, serving as a more physically motivated replacement to the conventional Dirichlet radial boundary conditions that are used in global simulation. It is shown that these new boundary conditions eliminate much of the numerical artefacts generated near the radial boundary when expressing a non-periodic function using a spectral basis. We then benchmark the new approach both linearly and non-linearly using a number of standard test cases.

Abstract:

Nonlinear multiscale gyrokinetic simulations of a Joint European Torus edge pedestal are used to show that electron-temperature-gradient (ETG) turbulence has a rich three-dimensional structure, varying strongly according to the local magnetic-field configuration. In the plane normal to the magnetic field, the steep pedestal electron temperature gradient gives rise to anisotropic turbulence with a radial (normal) wavelength much shorter than in the binormal direction. In the parallel direction, the location and parallel extent of the turbulence are determined by the variation in the magnetic drifts and finite-Larmor-radius (FLR) effects. The magnetic drift and FLR topographies have a perpendicular-wavelength dependence, which permits turbulence intensity maxima near the flux-surface top and bottom at longer binormal scales, but constrains turbulence to the outboard midplane at shorter electron-gyroradius binormal scales. Our simulations show that long-wavelength ETG turbulence does not transport heat efficiently, and significantly decreases overall ETG transport—in our case by ∼40%—through multiscale interactions.

Abstract:

<jats:title>Abstract</jats:title> <jats:p>Nonlinear multiscale gyrokinetic simulations of a Joint European Torus edge pedestal are used to show that electron-temperature-gradient (ETG) turbulence has a rich three-dimensional structure, varying strongly according to the local magnetic-field configuration. In the plane normal to the magnetic field, the steep pedestal electron temperature gradient gives rise to anisotropic turbulence with a radial (normal) wavelength much shorter than in the binormal direction. In the parallel direction, the location and parallel extent of the turbulence are determined by the variation in the magnetic drifts and finite-Larmor-radius (FLR) effects. The magnetic drift and FLR topographies have a perpendicular-wavelength dependence, which permits turbulence intensity maxima near the flux-surface top and bottom at longer binormal scales, but constrains turbulence to the outboard midplane at shorter electron-gyroradius binormal scales. Our simulations show that long-wavelength ETG turbulence does not transport heat efficiently, and significantly decreases overall ETG transport—in our case by ∼40%—through multiscale interactions.</jats:p>

Abstract:

The first experimental campaigns have proven that, due to the optimization of the magnetic configuration with respect to neoclassical transport, the contribution of turbulence is essential to understand and predict the total particle and energy transport in Wendelstein 7-X (W7-X). This has spurred much work on gyrokinetic modelling for the interpretation of the available experimental results and for the preparation of the next campaigns. At the same time, new stellarator gyrokinetic codes have just been or are being developed. It is therefore desirable to have a sufficiently complete, documented and verified set of gyrokinetic simulations in W7-X geometry against which new codes or upgrades of existing codes can be tested and benchmarked. This paper attempts to provide such a set of simulations in the form of a comprehensive benchmark between the recently developed code stella and the well-established code GENE. The benchmark consists of electrostatic gyrokinetic simulations in the W7-X magnetic geometry and includes different flux tubes, linear ion-temperature-gradient (ITG) and trapped-electron-mode stability analyses, computation of linear zonal-flow responses and calculation of ITG-driven heat fluxes.

Abstract:

The theory of magnetohydrodynamic (MHD) turbulence predicts that Alfvénic and slow-mode-like compressive fluctuations are energetically decoupled at small scales in the inertial range. The partition of energy between these fluctuations determines the nature of dissipation, which, in many astrophysical systems, happens on scales where plasma is collisionless. However, when the magnetorotational instability (MRI) drives the turbulence, it is difficult to resolve numerically the scale at which both types of fluctuations start to be decoupled because the MRI energy injection occurs in a broad range of wavenumbers, and both types of fluctuations are usually expected to be coupled even at relatively small scales. In this study, we focus on collisional MRI turbulence threaded by a near-azimuthal mean magnetic field, which is naturally produced by the differential rotation of a disc. We show that, in such a case, the decoupling scales are reachable using a reduced MHD model that includes differential-rotation effects. In our reduced MHD model, the Alfvénic and compressive fluctuations are coupled only through the linear terms that are proportional to the angular velocity of the accretion disc. We numerically solve for the turbulence in this model and show that the Alfvénic and compressive fluctuations are decoupled at the small scales of our simulations as the nonlinear energy transfer dominates the linear coupling below the MRI-injection scale. In the decoupling scales, the energy flux of compressive fluctuations contained in the small scales is almost double that of Alfvénic fluctuations. Finally, we discuss the application of this result to prescriptions of ion-to-electron heating ratio in hot accretion flows.