Search

Drift-reduced MHD models are widely used to study magnetised plasma phenomena, in particular for magnetically confined fusion applications, as well as in solar and astrophysical research. This letter discusses the choice of Ohm's law in these models, the resulting dispersion relations for the dynamics parallel to the magnetic field, and the implications for numerical simulations. We find that if electron pressure is included in Ohm's law, then both electromagnetic and finite electron mass effects must also be included in order to obtain physical dispersion relations. A simple modification to the plasma vorticity is also found which improves handling of low density regions, of particular relevance to the simulation of the boundary region of magnetised plasmas., Comment: 21 pages, 3 figures. To be submitted to PPCF

Due to their high cross-field mobility, neutrals can contribute to momentum transport even at the low relative densities found inside the separatrix and they can generate intrinsic rotation. We use a charge-exchange dominated solution to the neutral kinetic equation, coupled to neoclassical ions, to evaluate the momentum transport due to neutrals. Numerical solutions to the drift-kinetic equation allow us to cover the full range of collisionality, including the intermediate levels typical of the tokamak edge. In the edge there are several processes likely to contribute to momentum transport in addition to neutrals. Therefore, we present here an interpretive framework that can evaluate the momentum transport through neutrals based on radial plasma profiles. We demonstrate its application by analysing the neutral angular momentum flux for an L-mode discharge in the ASDEX Upgrade tokamak. The magnitudes of the angular momentum fluxes we find here due to neutrals of up to $1{-}2\;\mathrm{N\, m}$ are comparable to the net torque on the plasma from neutral beam injection, indicating the importance of neutrals for rotation in the edge., Comment: To be submitted to Nuclear Fusion

A theoretical framework is developed to clarify the relation between the profiles of density and temperature in the Scrape Off Layer (SOL) with the fluctuations (filaments) that generate them. The framework is based on the dynamics of independent filaments and on their statistical behaviour and can be used to rigorously understand the mechanisms that lead to the non-exponential nature of the radial SOL profiles as well as the increase of the relative fluctuation amplitude in the far SOL. Several models for the dynamics of the filaments, which can be applied to the framework, are derived and discussed for the purpose of identifying how different assumptions lead to the emergence of features in the profiles. It is found that multiple alternative models can explain the observations, thus motivating more stringent and focused experimental analysis. In particular, radially accelerating filaments, less efficient parallel exhaust and also a statistical distribution of the velocity of the filaments can all contribute to induce flatter profiles in the far SOL. A quite general result is the resiliency of the non-exponential nature of the profiles. At the same time, several of the models discussed can also capture the increase of the relative fluctuation amplitude observed in the far SOL. It is also shown that several scenarios are compatible with the broadening of the SOL, which could be caused by charge exchange interactions with neutral particles or by a significant radial acceleration of the filaments., Comment: 28 pages, 5 figures. This is an author-created, un-copyedited version of an article submitted for publication in Plasma Physics and Controlled Fusion. IoP Publishing Ltd is not responsible for any errors or omissions in this version of the manuscript or any version derived from it. arXiv admin note: text overlap with arXiv:1604.08790

Simulations have been carried out to establish how electron thermal physics, introduced in the form of a dynamic electron temperature, affects isolated filament motion and dynamics in 3D. It is found that thermal effects impact filament motion in two major ways when the filament has a significant temperature perturbation compared to its density perturbation: They lead to a strong increase in filament propagation in the bi-normal direction and a significant decrease in net radial propagation. Both effects arise from the temperature dependence of the sheath current which leads to a non-uniform floating potential, with the latter effect supplemented by faster pressure loss. The reduction in radial velocity can only occur when the filament cross-section loses angular symmetry. The behaviour is observed across different filament sizes and suggests that filaments with much larger temperature perturbations than density perturbations are more strongly confined to the near SOL region., Comment: 24 pages, 15 figures. This is an author-created, un-copyedited version of an article submitted for publication in Plasma Physics and Controlled Fusion. IoP Publishing Ltd is not responsible for any errors or omissions in this version of the manuscript or any version derived from it

A theoretical framework is developed to link the density profiles in the Scrape Off Layer (SOL) with the fluctuations (filaments) that generate them. The framework is based on the dynamics of independent filaments and their statistical behaviour and can be used to rigorously understand the mechanisms that lead to flattening and broadening of the SOL profiles as well as the radial increase of the relative fluctuation amplitude., Comment: 12 pages, 1 figure

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., Comment: Submitted to Nuclear Fusion

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 figures

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 $\boldsymbol{E}\times\boldsymbol{B}$ 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., Comment: This is an author-created, un-copyedited version of an article accepted for publication/published in Plasma Physics and Controlled Fusion. IOP Publishing Ltd is not responsible for any errors or omissions in this version of the manuscript or any version derived from it. The Version of Record is available online at http://dx.doi.org/10.1088/0741-3335/58/1/014030

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., Comment: 16 pages, 8 figures, submitted to Journal of Plasma Physics

Non-local closures allow kinetic effects on parallel transport to be included in fluid simulations. This is especially important in the scrape-off layer, but to be useful there the non-local model requires consistent kinetic boundary conditions at the sheath. A non-local closure scheme based on solution of a kinetic equation using a diagonalized moment expansion has been previously reported. We derive a method for imposing kinetic boundary conditions in this scheme and discuss their implementation in BOUT++. To make it feasible to implement the boundary conditions in the code, we are lead to transform the non-local model to a different moment basis, better adapted to describe parallel dynamics. The new basis has the additional benefit of enabling substantial optimization of the closure calculation, resulting in an O(10) speedup of the non-local code.