Source author record

Ben Dudson

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

7works
2topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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.

Published work

7 published item(s)

preprint2022arXiv

Comparison between MAST-U conventional and Super-X configurations through SOLPS-ITER modelling

MAST-U has recently started operating with a Super-X divertor, designed to increase total flux expansion and neutral trapping, both predicted through simple analytic models and SOLPS calculations to reduce the plasma and impurity density detachment thresholds. In this study, utilising the SOLPS-ITER code, we are quantifying the possible gain allowed by the MAST-U Super-X and neutral baffling geometry, in terms of access to detachment. We show that a significant reduction of the upstream density detachment threshold (up to a factor 1.6) could be achieved in MAST-U, for the Super-X, as opposed to conventional divertor geometry, mainly through an increased total flux expansion, neutral trapping being found very similar between the different configurations. We also show that variations of the strike-point angle are complex to interpret in such a tightly baffled geometry, and that a case in which the target normal points more towards the separatrix does not necessarily imply a lower detachment threshold. As in previous calculations for TCV, we quantify the role of neutral effects through developing and applying a quantitative definition of neutral trapping.

preprint2016arXiv

Blob dynamics in TORPEX poloidal null configurations

Three dimensional blob dynamics are simulated in X-point magnetic configurations in the TORPEX device via a non-field-aligned coordinate system, using an isothermal model which evolves density, vorticity, parallel velocity and parallel current density. By modifying the parallel gradient operator to include perpendicular perturbations from poloidal field coils, numerical singularities associated with field aligned coordinates are avoided. A comparison with a previously developed analytical model is performed and an agreement is found with minimal modification. Experimental comparison determines that the null region can cause an acceleration of filaments due to increasing connection length, but this acceleration is small relative to other effects, which we quantify. Experimental measurements are reproduced, and the dominant acceleration mechanism is identified as that of a developing dipole in a moving background. Contributions from increasing connection length close to the null point are a small correction.

preprint2016arXiv

Towards nonaxisymmetry; initial results using the Flux Coordinate Independent method in BOUT++

Fluid simulation of stellarator edge transport is difficult due to the complexities of mesh generation; the stochastic edge and strong nonaxisymmetry inhibit the use of field aligned coordinate systems. The recent implementation of the Flux Coordinate Independent method for calculating parallel derivatives in BOUT++ has allowed for more complex geometries. Here we present initial results of nonaxisymmetric diffusion modelling as a step towards stellarator turbulence modelling. We then present initial (non-turbulent) transport modelling using the FCI method and compare the results with analytical calculations. The prospects for future stellarator transport and turbulence modelling are discussed.

preprint2016arXiv

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 coordinate systems and the shifted metric procedure which is used to handle highly sheared grids. The Flux Coordinate Independent (FCI) 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, 3-field reduced MHD in 3D toroidal coordinates, and 5-field reduced MHD in slab geometry.

preprint2015arXiv

On the validity of drift-reduced fluid models for tokamak plasma simulation

Drift-reduced plasma fluid models are commonly used in plasma physics for analytics and simulations; however, the validity of such models must be verified for the regions of parameter space in which tokamak plasmas exist. By looking at the linear behaviour of drift-reduced and full-velocity models one can determine that the physics lost through the simplification that the drift-reduction provides is important in the core region of the tokamak. It is more acceptable for the edge-region but one must determine specifically for a given simulation if such a model is appropriate.

preprint2013arXiv

Non-local approach to kinetic effects on parallel transport in fluid models of the scrape-off layer

By using a non-local model, fluid simulations can capture kinetic effects in the parallel electron heat-flux better than is possible using flux limiters in the usual diffusive models. Non-local and diffusive models are compared using a test case representative of an ELM crash in the JET SOL, simulated in one dimension. The non-local model shows substantially enhanced electron temperature gradients, which cannot be achieved using a flux limiter. The performance of the implementation, in the BOUT++ framework, is also analysed to demonstrate its suitability for application in three-dimensional simulations of turbulent transport in the SOL.

preprint2012arXiv

Improved Nonlinear Solvers in BOUT++

Challenging aspects of large-scale turbulent edge simulations in plasma physics include robust nonlinear solvers and efficient preconditioners. This paper presents recent advances in the scalable solution of nonlinear partial differential equations in BOUT++, with emphasis on simulations of edge localized modes in tokamaks. A six-field, nonlinear, reduced magnetohydrodynamics model containing the fast shear Alfven wave and electron and ion heat conduction along magnetic fields is solved by using Jacobian-free Newton-Krylov (JFNK) methods and nonlinear GMRES (NGMRES). Physics-based preconditioning based on analytic Schur factorization of a simplified Jacobian is found to result in an order of magnitude reduction in runtime over unpreconditioned JFNK, and NGMRES is shown to significantly reduce runtime while requiring only the nonlinear function evaluation. We describe in detail the preconditioning algorithm, and we discuss parallel performance of NGMRES and Newton-Krylov methods using the PETSc library.