N. Walkden

Characterization of 3D filament dynamics in a MAST SOL flux tube geometry

Non-linear simulations of filament propagation in a realistic MAST SOL flux tube geometry using the BOUT++ fluid modelling framework show an isolation of the dynamics of the filament in the divertor region from the midplane region due to three features of the magnetic geometry; the variation of magnetic curvature along the field line, the expansion of the flux tube and strong magnetic shear. Of the three effects, the latter two lead to a midplane ballooning feature of the filament, whilst the former leads to a ballooning around the X-points. In simulations containing all three effects the filament is observed to balloon at the midplane, suggesting that the role of curvature variation is sub-dominant to the flux expansion and magnetic shear. The magnitudes of these effects are all strongest near the X-point which leads to the formation of parallel density gradients. The filaments simulated, which represent filaments in MAST, are identified as resistive ballooning, meaning that their motion is inertially limited, not sheath limited. Parallel density gradients can drive the filament towards a Boltzmann response when the collisionalityof the plasma is low. The results here show that the formation of parallel density gradients is a natural and inevitable consequence of a realistic magnetic geometry and therefore the transition to the Boltzmann response is a consequence of the use of realistic magnetic geometry and does not require initializing specifically varying background profiles as in slab simulations. The filaments studied here are stable to the linear resistive drift-wave instability but are subject to the non-linear effects associated with the Boltzmann response, particularly Boltzmann spinning. The Boltzmann response causes the filament to spin on an axis. In later stages of its evolution a non-linear turbulent state develops where the vorticity evolves into a turbulent eddy field on the same length scale as the parallel current. The transition from interchange motion to the Boltzmann response occurs with increasing temperature through a decrease in collisionality. This is confirmed by measuring the correlation between density and potential perturbations within the filament, which is low in the antisymmetric state associated with the interchange mechanism, but high in the Boltzmann regime. In the Boltzmann regime net radial transport is drastically reduced whilst a small net toroidal transport is observed. This suggests that only a subset of filaments, those driven by the interchange mechanism at the separatrix, can propagate into the far SOL. Filaments in the Boltzmann regime will be confined to the near separatrix region and quickly disperse. It is plausible that filaments in both regimes can contribute to the SOL transport observed in experiment; the former by propagating the filament into the far SOL and the latter by dispersion of the density within the filament.

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

Profile measurements in the plasma edge of mega amp spherical tokamak using a ball pen probe

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

A novel flexible field-aligned coordinate system for tokamak edge plasma simulation

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.

Influence of plasma background on 3D scrape-off layer filaments

The ball pen probe (BPP) technique is used successfully to make profile measurements of plasma potential, electron temperature, and radial electric field on the Mega Amp Spherical Tokamak. The potential profile measured by the BPP is shown to significantly differ from the floating potential both in polarity and profile shape. By combining the BPP potential and the floating potential, the electron temperature can be measured, which is compared with the Thomson scattering (TS) diagnostic. Excellent agreement between the two diagnostics is obtained when secondary electron emission is accounted for in the floating potential. From the BPP profile, an estimate of the radial electric field is extracted which is shown to be of the order ∼1 kV/m and increases with plasma current. Corrections to the BPP measurement, constrained by the TS comparison, introduce uncertainty into the ER measurements. The uncertainty is most significant in the electric field well inside the separatrix. The electric field is used to estimate toroidal and poloidal rotation velocities from E × B motion. This paper further demonstrates the ability of the ball pen probe to make valuable and important measurements in the boundary plasma of a tokamak.

Filamentary velocity scaling validation in the TCV tokamak

Abstract Tokamak plasmas are confined by a magnetic field that limits the particle and heat transport perpendicular to the field. Parallel to the field the ionised particles can move freely, so to obtain confinement the field lines are “closed” (i.e. form closed surfaces of constant poloidal flux) in the core of a tokamak. Towards, the edge, however, the field lines intersect physical surfaces, leading to interaction between neutral and ionised particles, and the potential melting of the material surface. Simulation of this interaction is important for predicting the performance and lifetime of future tokamak devices such as ITER. Field-aligned coordinates are commonly used in the simulation of tokamak plasmas due to the geometry and magnetic topology of the system. However, these coordinates are limited in the geometry they allow in the poloidal plane due to orthogonality requirements. A novel 3D coordinate system is proposed herein that relaxes this constraint so that any arbitrary, smoothly varying geometry can be matched in the poloidal plane while maintaining a field-aligned coordinate. This system is implemented in BOUT++ and tested for accuracy using the method of manufactured solutions. A MAST edge cross-section is simulated using a fluid plasma model and the results show expected behaviour for density, temperature, and velocity. Finally, simulations of an isolated divertor leg are conducted with and without neutrals to demonstrate the ion-neutral interaction near the divertor plate and the corresponding beneficial decrease in plasma temperature.

Modification of SOL profiles and fluctuations with line-average density and divertor flux expansion in TCV

This paper presents the effect of self-consistent plasma backgrounds including plasma-neutral interactions, on the dynamics of filament propagation. The principle focus is on the influence of the neutrals on the filament through both direct interactions and through their influence on the plasma background. Both direct and indirect interactions influence the motion of filaments. A monotonic increase of filament peak velocity with upstream electron temperature is observed, while a decrease with increasing electron density is observed. If ordered by the target temperature, the density dependence disappears and the filament velocity is only a function of the target temperature. Smaller filaments keep a density dependence, as a result of the density dependence of the plasma viscosity. The critical size

Influence of plasma background including neutrals on scrape-off layer filaments using 3D simulations

\delta^*