Luke Easy

Three dimensional simulations of plasma filaments in the scrape off layer: A comparison with models of reduced dimensionality

This paper presents simulations of isolated 3D filaments in a slab geometry obtained using a newly developed 3D reduced fluid code, written using the BOUT++ framework. First, systematic scans were performed to investigate how the dynamics of a filament are affected by its amplitude, perpendicular size, and parallel extent. The perpendicular size of the filament was found to have a strong influence on its motions, as it determined the relative importance of parallel currents to polarization and viscous currents, whilst drift-wave instabilities were observed if the initial amplitude of the blob was increased sufficiently. Next, the 3D simulations were compared to 2D simulations using different parallel closures; namely, the sheath dissipation closure, which neglects parallel gradients, and the vorticity advection closure, which neglects the influence of parallel currents. The vorticity advection closure was found to not replicate the 3D perpendicular dynamics and overestimated the initial radial acceleratio...

BOUT++: Recent and current developments

BOUT++ is a 3D nonlinear finite-difference plasma simulation code, capable of solving quite general systems of Partial Differential Equations (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 Magnetohydrodynamics (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 (the Portable, Extensible Toolkit for Scientific Computation) (Balay et al. 2014), 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.

Numerical investigation of isolated filament motion in a realistic tokamak geometry

This paper presents a numerical investigation of isolated filament dynamics in a simulation geometry representative of the scrape-off layer (SOL) of the Mega Amp Spherical Tokamak (MAST) previously studied in [N.R.Walkden et.al, Plasma Phys. Control. Fusion, 55 (2013) 105005]. This paper focuses on the evolution of filament cross-sections at the outboard midplane and investigates the scaling of the centre of mass velocity of the filament cross-section with filament width and electron temperature. By decoupling the vorticity equation into even and odd parity components about the centre of the filament in the bi-normal direction parallel density gradients are shown to drive large velocities in the bi-normal (approximately poloidal) direction which scale linearly with electron temperature. In this respect increasing the electron temperature causes a departure of the filament dynamics from 2D behaviours. Despite the strong impact of 3D effects the radial motion of the filament is shown to be relatively well predicted by 2D scalings. The radial velocity is found to scale positively with both electron temperature and cross-sectional width, suggesting an inertially limited nature. Comparison with the two-region model [J. R. Myra et.al, Phys.Plasmas, 13 (2006) 112502] achieves reasonable agreement when using a corrected parallel connection length due to the neglect of diamagnetic currents driven in the divertor region of the filament. Analysis of the transport of particles due to the motion of the filament shows that the background temperature has a weak overall effect on the radial particle flux whilst the filament width has a strong effect.

The effects of shape and amplitude on the velocity of scrape-off layer filaments

A complete model of the dynamics of scrape-off layer filaments will be rather complex, including temperature evolution, three dimensional geometry and finite Larmor radius effects. However, the basic mechanism of

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

\boldsymbol{E}\times\boldsymbol{B}

Verification of BOUT++ by the method of manufactured solutions

advection due to electrostatic potential driven by the diamagnetic current can be captured in a much simpler model; a complete understanding of the physics in the simpler model will then aid interpretation of more complex simulations, by allowing the new effects to be disentangled. Here we consider such a simple model, which assumes cold ions and isothermal electrons and is reduced to two dimensions. We derive the scaling with width and amplitude of the velocity of isolated scrape-off layer filaments, allowing for arbitrary elliptical cross-sections, where previously only circular cross-sections have been considered analytically. We also put the scaling with amplitude in a new and more satisfactory form. The analytical results are extensively validated with two dimensional simulations and also compared, with reasonable agreement, to three dimensional simulations having minimal variation parallel to the magnetic field.

Investigation of the effect of resistivity on scrape off layer filaments using three-dimensional simulations

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

Multi-code analysis of scrape-off layer filament dynamics in MAST

{{k}_{\parallel}}=0

Dynamics of 3D isolated thermal filaments

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.

On the interaction of scrape off layer filaments

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.