Characterisation of the L-mode Scrape Off Layer in MAST: decay lengths

This work presents a detailed characterisation of the MAST Scrape Off Layer in L-mode. Scans in line averaged density, plasma current and toroidal magnetic field were performed. A comprehensive and integrated study of the SOL was allowed by the use of a wide range of diagnostics. In agreement with previous results, an increase of the line averaged density induced a broadening of the midplane density profile.Comment: 30 pages, 11 figure

Verification of BOUT++ by the method of manufactured solutions

BOUT++ is a software package designed for solving plasma fluid models. It has been used to simulate a wide range of plasma phenomena ranging from linear stability analysis to 3D plasma turbulence and is capable of simulating a wide range of drift-reduced plasma fluid and gyro-fluid models. A verification exercise has been performed as part of a EUROfusion Enabling Research project, to rigorously test the correctness of the algorithms implemented in BOUT++, by testing order-of-accuracy convergence rates using the Method of Manufactured Solutions (MMS). We present tests of individual components including time-integration and advection schemes, non-orthogonal toroidal field-aligned coordinate systems and the shifted metric procedure which is used to handle highly sheared grids. The flux coordinate independent approach to differencing along magnetic field-lines has been implemented in BOUT++ and is here verified using the MMS in a sheared slab configuration. Finally, we show tests of three complete models: 2-field Hasegawa-Wakatani in 2D slab, 3-field reduced magnetohydrodynamics (MHD) in 3D field-aligned toroidal coordinates, and 5-field reduced MHD in slab geometry

Kinetic and fluid simulations of parallel electron transport during equilibria and transients in the scrape-off layer

We present the first parallel electron transport results obtained using the newly developed 1D transport code SOL-KiT. With the capability to switch between consistent kinetic and fluid models for the electrons, we explore and report the differences in both equilibrium and transient simulations. Significant kinetic effects are found during transients, especially in the behaviour of the electron sheath heat transmission coefficient, which shows up to an eightfold increase. Equilibria are obtained for an input power scan with parameters relevant to medium size tokamaks. Detached equilibria are found to persist to higher input powers when electrons are treated kinetically. Furthermore, non-monotonic behaviour of the electron sheath heat transmission coefficient is observed in the power scan, with values being up to 40% above the classical value. We discuss the implications of the p

BOUT++ : Recent and current developments

BOUT++ is a 3D nonlinear finite-difference plasma simulation code, capable of solving quite general systems of PDEs, but targeted particularly on studies of the edge region of tokamak plasmas. BOUT++ is publicly available, and has been adopted by a growing number of researchers worldwide. Here we present improvements which have been made to the code since its original release, both in terms of structure and its capabilities. Some recent applications of these methods are reviewed, and areas of active development are discussed. We also present algorithms and tools which have been developed to enable creation of inputs from analytic expressions and experimental data, and for processing and visualisation of output results. This includes a new tool Hypnotoad for the creation of meshes from experimental equilibria. Algorithms have been implemented in BOUT++ to solve a range of linear algebraic problems encountered in the simulation of reduced MHD and gyro-fluid models: A preconditioning scheme is presented which enables the plasma potential to be calculated efficiently using iterative methods supplied by the PETSc library, without invoking the Boussinesq approximation. Scaling studies are also performed of a linear solver used as part of physics-based preconditioning to accelerate the convergence of implicit time-integration schemes

Hermes : global plasma edge fluid turbulence simulations

The transport of heat and particles in the relatively collisional edge regions of magnetically confined plasmas is a scientifically challenging and technologically important problem. Understanding and predicting this transport requires the self-consistent evolution of plasma fluctuations, global profiles and flows, but the numerical tools capable of doing this in realistic (diverted) geometry are only now being developed. Here a 5-field reduced 2-fluid plasma model for the study of instabilities and turbulence in magnetised plasmas is presented, built on the BOUT++ framework. This cold ion model allows the evolution of global profiles, electric fields and flows on transport timescales, with flux-driven cross-field transport determined self-consistently by electromagnetic turbulence. Developments in the model formulation and numerical implementation are described, and simulations are performed in poloidally limited and diverted tokamak configurations

Blob dynamics in the TORPEX experiment:a multi-code validation

Three-dimensional and two-dimensional seeded blob simulations are performed with five different fluid models, all based on the drift-reduced Braginskii equations, and the numerical results are compared among themselves and validated against experimental measurements provided by the TORPEX device (Fasoli et al 2006 Phys. Plasmas 13 055902). The five models are implemented in four simulation codes, typically used to simulate the plasma dynamics in the tokamak scrape-off layer, namely BOUT++ (Dudson et al 2009 Comput. Phys. Commun. 180 1467), GBS (Ricci et al 2012 Plasma Phys. Control. Fusion 54 124047), HESEL (Nielsen et al 2015 Phys. Lett. A 379 3097), and TOKAM3X (Tamain et al 2014 Contrib. Plasma Phys. 54 555). Three blobs with different velocities and different stability properties are simulated. The differences observed among the simulation results and the different levels of agreement with experimental measurements are investigated, increasing our confidence in our simulation tools and shedding light on the blob dynamics. The comparisons demonstrate that the radial blob dynamics observed in the three-dimensional simulations is in good agreement with experimental measurements and that, in the present experimental scenario, the two-dimensional model derived under the assumption of ${{k}_{\parallel}}=0$ is able to recover the blob dynamics observed in the three-dimensional simulations. Moreover, it is found that an accurate measurement of the blob temperature is important to perform reliable seeded blob simulations