Source author record

B. D. Dudson

B. D. Dudson appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

physics.plasm-ph physics.comp-ph

Catalog footprint

What is connected

5works

2topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2015arXiv

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.

preprint2014arXiv

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.

preprint2014arXiv

Non-local parallel transport in BOUT++

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.

preprint2010arXiv

Simulation of Edge Localised Modes using BOUT++

The BOUT++ code is used to simulate ELMs in a shifted circle equilibrium. Reduced ideal MHD simulations are first benchmarked against the linear ideal MHD code ELITE, showing good agreement. Diamagnetic drift effects are included finding the expected suppression of high toroidal mode number modes. Nonlinear simulations are performed, making the assumption that the anomalous kinematic electron viscosity is comparable to the anomalous electron thermal diffusivity. This allows simulations with realistically high Lundquist numbers S = 1e8, finding ELM sizes of 5-10% of the pedestal stored thermal energy. Scans show a strong dependence of the ELM size resistivity at low Lundquist numbers, with higher resistivity leading to more violent eruptions. At high Lundquist numbers relevant to high-performance discharges, ELM size is independent of resistivity as hyper-resistivity becomes the dominant dissipative effect.

preprint2008arXiv

BOUT++: a framework for parallel plasma fluid simulations

A new modular code called BOUT++ is presented, which simulates 3D fluid equations in curvilinear coordinates. Although aimed at simulating Edge Localised Modes (ELMs) in tokamak X-point geometry, the code is able to simulate a wide range of fluid models (magnetised and unmagnetised) involving an arbitrary number of scalar and vector fields, in a wide range of geometries. Time evolution is fully implicit, and 3rd-order WENO schemes are implemented. Benchmarks are presented for linear and non-linear problems (the Orszag-Tang vortex) showing good agreement. Performance of the code is tested by scaling with problem size and processor number, showing efficient scaling to thousands of processors. Linear initial-value simulations of ELMs using reduced ideal MHD are presented, and the results compared to the ELITE linear MHD eigenvalue code. The resulting mode-structures and growth-rate are found to be in good agreement (BOUT++ = 0.245, ELITE = 0.239). To our knowledge, this is the first time dissipationless, initial-value simulations of ELMs have been successfully demonstrated.