Showing 1–37 of 37 results for author: Dudson, B

Coupling Fluid Plasma and Kinetic Neutral Models using Correlated Monte Carlo Methods

Authors: Gregory J. Parker, Maxim V. Umansky, Benjamin D. Dudson

Abstract: While boundary plasmas in present day tokamaks generally fall in a fluid regime, neutral species near the boundary often require kinetic models due to long mean-free-paths compared to characteristic spatial scales in the region. Monte-Carlo (MC) methods provide a complete, high-fidelity approach to solving kinetic models, and must be coupled to fluid plasma models to simulate the full plasma-neutr… ▽ More While boundary plasmas in present day tokamaks generally fall in a fluid regime, neutral species near the boundary often require kinetic models due to long mean-free-paths compared to characteristic spatial scales in the region. Monte-Carlo (MC) methods provide a complete, high-fidelity approach to solving kinetic models, and must be coupled to fluid plasma models to simulate the full plasma-neutrals system. The statistical nature of MC methods, however, prevents convergence of coupled fluid-kinetic simulations to an exact self-consistent steady-state. Moreover, this forces the use of explicit methods that can suffer from numerical errors and require huge computational resources. Correlated Monte-Carlo (CMC) methods are expected to alleviate these issues, but have historically enjoyed only mixed success. Here, a fully implicit method for coupled plasma-neutral systems is demonstrated in 1D using the UEDGE plasma code and a homemade CMC code. In particular, it is shown that ensuring the CMC method is a differentiable function of the background plasma is sufficient to employ a Jacobian-Free Newton-Krylov solver for implicit time steps. The convergence of the implicit coupling method is explored and compared with explicit coupling and uncorrelated methods. It is shown that ensuring differentiability by controlling random seeds in the MC is sufficient to achieve convergence, and that the use of implicit time-stepping methods has the potential for improved stability and runtimes over explicit coupling methods. △ Less

Submitted 15 July, 2024; originally announced July 2024.

Comments: 6 pages, 6 figures. Comments welcome!

Flexible Stellarator Physics Facility

Authors: F. I. Parra, S. -G. Baek, M. Churchill, D. R. Demers, B. Dudson, N. M. Ferraro, B. Geiger, S. Gerhardt, K. C. Hammond, S. Hudson, R. Jorge, E. Kolemen, D. M. Kriete, S. T. A. Kumar, M. Landreman, C. Lowe, D. A. Maurer, F. Nespoli, N. Pablant, M. J. Pueschel, A. Punjabi, J. A. Schwartz, C. P. S. Swanson, A. M. Wright

Abstract: We propose to build a Flexible Stellarator Physics Facility to explore promising regions of the vast parameter space of disruption-free stellarator solutions for Fusion Pilot Plants (FPPs). We propose to build a Flexible Stellarator Physics Facility to explore promising regions of the vast parameter space of disruption-free stellarator solutions for Fusion Pilot Plants (FPPs). △ Less

Submitted 4 July, 2024; originally announced July 2024.

Comments: White paper submitted to FESAC subcommittee on Facilities, 8 pages

Detachment scalings derived from 1D scrape-off-layer simulations

Authors: Thomas Body, Thomas Eich, Adam Q Kuang, Thomas Looby, Mike Kryjak, Benjamin D Dudson, Matthew Reinke

Abstract: Fusion power plants will require active detachment control to mitigate sputtering. Although detachment control has been demonstrated on current fusion experiments, achieving this on fusion power plants will be more challenging due to their reduced diagnostic set. Real-time scrape-off-layer (SOL) modeling will help to fill this "diagnostic gap". To determine the minimal set of physics required for… ▽ More Fusion power plants will require active detachment control to mitigate sputtering. Although detachment control has been demonstrated on current fusion experiments, achieving this on fusion power plants will be more challenging due to their reduced diagnostic set. Real-time scrape-off-layer (SOL) modeling will help to fill this "diagnostic gap". To determine the minimal set of physics required for accurate SOL modeling, we use the configurable Hermes-3 edge modeling framework to perform simple, fixed-fraction-impurity 1D Braginskii simulations. We show that these simulations reproduce distinctive attached and detached conditions, as well as the characteristic ion-flux rollover. We also perform scans of the input heat flux and impurity concentration, and show that the results predict similar trends to the lower-fidelity Lengyel model. This allows us to also indirectly compare to SOLPS simulations, and we find that the Lengyel model has similar "over-estimation" factors to Hermes-1D and SOLPS. This indirect comparison highlights a series of extensions to bring our Hermes-1D results closer to the higher-fidelity 2D modeling, with the goal of developing time-dependent detachment control suitable for next-generation fusion devices. △ Less

Submitted 24 June, 2024; originally announced June 2024.

Comments: Conference paper for PSI-26 (26th Conference on Plasma Surface Interactions, Marseille)

Global fluid turbulence simulations in the SOL of a stellarator island divertor

Authors: Brendan Shanahan, David Bold, Ben Dudson

Abstract: Isothermal fluid turbulence simulations have been performed in the edge and scrape-off-layer (SOL) of an analytic stellarator configuration with an island divertor, thereby providing numerical insight into edge turbulence in regions around islands in a stellarator. The steady-state transport follows the 1/R curvature drive toward the outboard side, but large fluctuations are present throughout the… ▽ More Isothermal fluid turbulence simulations have been performed in the edge and scrape-off-layer (SOL) of an analytic stellarator configuration with an island divertor, thereby providing numerical insight into edge turbulence in regions around islands in a stellarator. The steady-state transport follows the 1/R curvature drive toward the outboard side, but large fluctuations are present throughout the island divertor region, with the average wavelength of similar size to the island width. The system exhibits a prominent m=2, n=5 mode, although other modes are present. The amplitude and radial extent of the density fluctuations are similar throughout the edge and SOL, but can decrease near island O-points. The fluctuations exhibit a predominantly positive skewness on the outboard midplane, indicating blob-like perturbations for the transport into the outer SOL. It is determined that a point on the separatrix is generally more correlated with regions outside of the SOL than a nearby reference point which does not lie on the separatrix. △ Less

Submitted 26 March, 2024; originally announced March 2024.

Comments: 19 pages, 11 Figures

Two-stage Crash Process in Resistive Drift Ballooning Mode Driven ELM Crash

Authors: Haruki Seto, Xueqiao Xu, Benjamin D. Dudson, Masatoshi Yagi

Abstract: We report a two-stage crash process in edge localized mode (ELM) driven by resistive drift-ballooning modes (RDBMs) numerically simulated in a full annular torus domain. In the early nonlinear phase, the first crash is triggered by linearly unstable RDBMs and m/n = 2/1 magnetic islands are nonlinearly excited via nonlinear couplings of RDBMs. Simultaneously, middle-n RDBM turbulence develops but i… ▽ More We report a two-stage crash process in edge localized mode (ELM) driven by resistive drift-ballooning modes (RDBMs) numerically simulated in a full annular torus domain. In the early nonlinear phase, the first crash is triggered by linearly unstable RDBMs and m/n = 2/1 magnetic islands are nonlinearly excited via nonlinear couplings of RDBMs. Simultaneously, middle-n RDBM turbulence develops but is poloidally localized around X-points of the magnetic islands, leading to the small energy loss. Here m is the poloidal mode number, n is the toroidal mode number, the q = 2 rational surface exists at the pressure gradient peak, and q is the safety factor, respectively. The second crash occurs in the late nonlinear phase. Low-n magnetic islands are also excited around the q = 2 surface via nonlinear couplings among the middle-n turbulence. Since the turbulence develops from the X-points of higher harmonics of m/n = 2/1 magnetic islands, it expands out poloidally. The second crash is triggered when the turbulence covers the whole poloidal region. A scan of toroidal wedge number N, where full torus is divided into N segments in the toroidal direction, also reveals that the first crash process becomes more prominent with the higher toroidal wedge number where the RDBMs play a dominant role. These results indicate that nonlinear interactions of all channels in the full torus domain can significantly affect the trigger dynamics of ELMs driven by the RDBMs. △ Less

Submitted 7 October, 2023; originally announced October 2023.

Comments: 27 pages, 15 figures. The article has been submitted to Physics of Plasmas

Hermes-3: Multi-component plasma simulations with BOUT++

Authors: Ben Dudson, Mike Kryjak, Hasan Muhammed, Peter Hill, John Omotani

Abstract: A new open source tool for fluid simulation of multi-component plasmas is presented, based on a flexible software design that is applicable to scientific simulations in a wide range of fields. This design enables the same code to be configured at run-time to solve systems of partial differential equations in 1D, 2D or 3D, either for transport (steady-state) or turbulent (time-evolving) problems, w… ▽ More A new open source tool for fluid simulation of multi-component plasmas is presented, based on a flexible software design that is applicable to scientific simulations in a wide range of fields. This design enables the same code to be configured at run-time to solve systems of partial differential equations in 1D, 2D or 3D, either for transport (steady-state) or turbulent (time-evolving) problems, with an arbitrary number of ion and neutral species. To demonstrate the capabilities of this tool, applications relevant to the boundary of tokamak plasmas are presented: 1D simulations of diveror plasmas evolving equations for all charge states of neon and deuterium; 2D transport simulations of tokamak equilibria in single-null X-point geometry with plasma ion and neutral atom species; and simulations of the time-dependent propagation of plasma filaments (blobs). Hermes-3 is publicly available on Github under the GPL-3 open source license. The repository includes documentation and a suite of unit, integrated and convergence tests. △ Less

Submitted 16 June, 2023; v1 submitted 21 March, 2023; originally announced March 2023.

Comments: Submitted to Computer Physics Communications

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

Authors: Alexandre Fil, Bruce Lipschultz, David Moulton, Andrew Thornton, Ben Dudson, Omkar Myatra, Kevin Verhaegh, EUROFusion MST1 team

Abstract: 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 geome… ▽ More 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. △ Less

Submitted 13 July, 2022; v1 submitted 25 February, 2022; originally announced February 2022.

On Ohm's law in reduced plasma fluid models

Authors: B. D. Dudson, S. L. Newton, J. T. Omotani, J. Birch

Abstract: 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 e… ▽ More 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. △ Less

Submitted 11 May, 2021; originally announced May 2021.

Comments: 21 pages, 3 figures. To be submitted to PPCF

Comparing Two- and Three-Dimensional Models of Scrape-Off-Layer Turbulent Transport

Authors: T. E. G. Nicholas, J. Omotani, F. Riva, F. Militello, B. Dudson

Abstract: There exists a large body of previous work using reduced two-dimensional models of the SOL, which model fluctuations in the drift-plane but approximate parallel transport with effective loss terms. Full size three-dimensional simulations of SOL turbulence in experimental geometries are now possible, but are far more computationally expensive than 2D models. We therefore use a flux-tube geometry mo… ▽ More There exists a large body of previous work using reduced two-dimensional models of the SOL, which model fluctuations in the drift-plane but approximate parallel transport with effective loss terms. Full size three-dimensional simulations of SOL turbulence in experimental geometries are now possible, but are far more computationally expensive than 2D models. We therefore use a flux-tube geometry model of the scrape-off layer to compare the results of 2D simulations to 3D simulations with a similar setup, looking for systematic differences. Overall there is good agreement in the basic radial profiles, probability distribution functions, and power spectra of fluctuations. However, the average temperature is over-predicted in 2D relative to 3D, and we explain the difference in terms of the effect of geometrical simplifications of devices at low power. Varying geometric parameters, we find that supersonic flow in the divertor leg, which occurs because our simulations do not include neutrals and so represent low-recycling conditions, means that the divertor leg length only has a weak effect on the output. Finally, we examine the effect of altering the magnitude of source and sink terms in 2D, concluding that they cannot easily be used to recreate both the density and temperature profiles observed in 3D simultaneously. △ Less

Submitted 19 March, 2021; v1 submitted 17 March, 2021; originally announced March 2021.

Comments: 14 pages, 24 figures, submitted to Plasma Physics and Controlled Fusion

Dynamics of scrape-off layer filaments in detached conditions

Authors: David Schwörer, Nick R. Walkden, Benjamin D. Dudson, Fulvio Militello, Huw Leggate, Miles M. Turner

Abstract: The here presented work studies the dynamics of filaments using 3D fluid simulations in the presence of detached background profiles. It was found that evolving the neutrals on the time-scale of the filament did not have a significant impact on the dynamics of the filament. In general a decreasing filament velocity with increasing plasma background density has been observed, with the exception of… ▽ More The here presented work studies the dynamics of filaments using 3D fluid simulations in the presence of detached background profiles. It was found that evolving the neutrals on the time-scale of the filament did not have a significant impact on the dynamics of the filament. In general a decreasing filament velocity with increasing plasma background density has been observed, with the exception of detachment onset, where a temporarily increase in radial velocity occurs. The decreasing trend with temporary increase was found for filaments around the critical size and larger, while smaller filaments where less affected by detachment. With detachment the critical filament size increased, as larger filaments were faster in detached conditions. This breaks the trend of attached conditions, where the critical size decreases with increasing density. △ Less

Submitted 16 March, 2020; v1 submitted 26 November, 2019; originally announced November 2019.

The role of particle, energy and momentum losses in 1D simulations of divertor detachment

Authors: B D Dudson, J Allen, T Body, B Chapman, C Lau, L Townley, D Moulton, J Harrison, B Lipschultz

Abstract: A new 1D divertor plasma code, SD1D, has been used to examine the role of recombination, radiation, and momentum exchange in detachment. Neither momentum or power losses by themselves are found to be sufficient to produce a reduction in target ion flux in detachment (flux rollover); radiative power losses are required to a) limit and reduce the ionization source and b) access low-target temperatur… ▽ More A new 1D divertor plasma code, SD1D, has been used to examine the role of recombination, radiation, and momentum exchange in detachment. Neither momentum or power losses by themselves are found to be sufficient to produce a reduction in target ion flux in detachment (flux rollover); radiative power losses are required to a) limit and reduce the ionization source and b) access low-target temperature, T_target, conditions for volumetric momentum losses. Recombination is found to play little role at flux rollover, but as T_target drops to temperatures around 1eV, it becomes a strong ion sink. In the case where radiative losses are dominated by hydrogen, the detachment threshold is identified as a minimum gradient of the energy cost per ionisation with respect to T_target. This is also linked to thresholds in T_target and in the ratio of upstream pressure to power flux. A system of determining the detached condition is developed such that the divertor solution at a given T_target (or lack of one) is determined by the simultaneous solution of two equations for target ion current - one dependent on power losses and the other on momentum. Depending on the detailed momentum and power loss dependence on temperature there are regions of T_target where there is no solution and the plasma 'jumps' from high to low T_target states. The novel analysis methods developed here provide an intuitive way to understand complex detachment phenomena, and can potentially be used to predict how changes in the seeding impurity used or recycling aspects of the divertor can be utilised to modify the development of detachment. △ Less

Submitted 21 December, 2018; originally announced December 2018.

Comments: 19 pages, 11 figures. Submitted to Plasma Phys. Control. Fusion

The effects of non-uniform drive on plasma filaments

Authors: Brendan Shanahan, Ben Dudson, Peter Hill

Abstract: Wendelstein 7-X core fueling is primarily achieved through pellet injection. The trajectory of plasmoids from an ablating pellet is an ongoing research question, which is complicated by the complex magnetic geometry of W7-X; curvature drive varies significantly toroidally, including a change in the drift drive direction. Here we use the Hermes model in BOUT++ to simulate cold plasma filaments in s… ▽ More Wendelstein 7-X core fueling is primarily achieved through pellet injection. The trajectory of plasmoids from an ablating pellet is an ongoing research question, which is complicated by the complex magnetic geometry of W7-X; curvature drive varies significantly toroidally, including a change in the drift drive direction. Here we use the Hermes model in BOUT++ to simulate cold plasma filaments in slab geometries where the magnetic drift drive is non-uniform along the field line. It is shown that if the field-line-averaged curvature drive is non-zero, a filament will propagate coherently in the direction of average drive. It is also shown that a non-uniform drive will provide a non- uniform propagation; an effect which is reduced at higher temperatures due to an increased sound speed along the field line. Finally, simulations with curvature similar to that found in Wendelstein 7-X are performed which indicate a slow radial propagation of plasma filaments despite a highly non-uniform drive profile. △ Less

Submitted 10 October, 2018; originally announced October 2018.

Fluid simulations of plasma filaments in stellarator geometries with BSTING

Authors: B Shanahan, B Dudson, P Hill

Abstract: Here we present first results simulating plasma filaments in non-axisymmetric geometries, using a fluid turbulence extension of the \boutxx~framework. This is made possible by the implementation of the Flux Coordinate Independent scheme for parallel derivatives, an extension of the metric tensor components which allows them to vary in three dimensions, and development of grid generation. Tests hav… ▽ More Here we present first results simulating plasma filaments in non-axisymmetric geometries, using a fluid turbulence extension of the \boutxx~framework. This is made possible by the implementation of the Flux Coordinate Independent scheme for parallel derivatives, an extension of the metric tensor components which allows them to vary in three dimensions, and development of grid generation. Tests have been performed to confirm that the extension to three dimensional metric tensors does not compromise the accuracy and stability of the associated numerical operators. Recent changes to the FCI grid generator in \boutxx, including a curvilinear grid system which allows for potentially more efficient computation, are also presented. Initial simulations of seeded plasma filaments in a non-axisymmetric geometry are reported. We characterize filaments propagating in the closed-field-line region of a low-field-period, rotating ellipse equilibrium as inertially-limited by examining the velocity scaling and currents associated with the filament propagation. Finally, it is shown that filaments in a non-axisymmetric rotating ellipse equilibrium propagate in a toroidally nonuniform fashion, and it is determined that the long connection lengths in the scrape-off-layer enable parallel gradients to establish, which has consequences for interpretation of experimental data. △ Less

Submitted 27 August, 2018; originally announced August 2018.

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

Authors: David Schwörer, Nicholas R Walkden, Huw Leggate, Ben D Dudson, Fulvio Militello, Turlough Downes, Miles M Turner

Abstract: 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 increas… ▽ More 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 $δ^*$, where filaments are fastest, is shifted to larger sizes for higher densities, due to the plasma viscosity. If the density dependence of the plasma viscosity is removed, $δ^*$ has no temperature dependence, but rather a density dependence. △ Less

Submitted 19 June, 2018; originally announced June 2018.

Generation of synthetic magnetized plasma turbulence by representation as a stochastic process

Authors: Jarrod Leddy, Chris Bowman, Kieran Gibson, Ben Dudson

Abstract: Plasma turbulence simulations are often computationally expensive with delicate numerical stability. Yet, long simulations are needed to generate uncorrelated turbulence data for studies such as microwave scattering through density perturbations. For this reason, alternative methods of producing accurate synthetic turbulence profiles via statistical methods is of interest. Such a method is propose… ▽ More Plasma turbulence simulations are often computationally expensive with delicate numerical stability. Yet, long simulations are needed to generate uncorrelated turbulence data for studies such as microwave scattering through density perturbations. For this reason, alternative methods of producing accurate synthetic turbulence profiles via statistical methods is of interest. Such a method is proposed where the two-point covariance function of the desired turbulence is used to construct a multi-variate normal distribution. Sampling from this distribution produces random fields which are both qualitatively and quantitatively similar to the input turbulence data set. The resulting `synthetic' turbulent profiles are uncorrelated in `time' so it is useful only for scenarios that do not require such physical evolution. △ Less

Submitted 16 November, 2017; originally announced November 2017.

Comments: 5 pages, 5 figures, submitting to CPC

Intrinsic suppression of turbulence in linear plasma devices

Authors: Jarrod Leddy, Ben Dudson

Abstract: Plasma turbulence is the dominant transport mechanism for heat and particles in magnetized plasmas in linear devices and tokamaks, so the study of turbulence is important in limiting and controlling this transport. Linear devices provide an axial magnetic field that serves to confine a plasma in cylindrical geometry as it travels along the magnetic field from the source to the strike point. Due to… ▽ More Plasma turbulence is the dominant transport mechanism for heat and particles in magnetized plasmas in linear devices and tokamaks, so the study of turbulence is important in limiting and controlling this transport. Linear devices provide an axial magnetic field that serves to confine a plasma in cylindrical geometry as it travels along the magnetic field from the source to the strike point. Due to perpendicular transport, the plasma density and temperature have a roughly Gaussian radial profile with gradients that drive instabilities, such as resistive drift-waves and Kelvin-Helmholtz. If unstable, these instabilities cause perturbations to grow resulting in saturated turbulence, increasing the cross-field transport of heat and particles. When the plasma emerges from the source, there is a time, $τ_{\parallel}$, that describes the lifetime of the plasma based on parallel velocity and length of the device. As the plasma moves down the device, it also moves azimuthally according to $E\times B$ and diamagnetic velocities. There is a balance point in these parallel and perpendicular times that sets the stabilisation threshold. We simulate plasmas with a variety of parallel lengths and magnetic fields to vary the parallel and perpendicular lifetimes, respectively, and find that there is a clear correlation between the saturated RMS density perturbation level and the balance between these lifetimes. The threshold of marginal stability is seen to exist where $τ_{\parallel}\approx11τ_{\perp}$. This is also associated with the product $τ_{\parallel}γ_*$, where $γ_*$ is the drift-wave linear growth rate, indicating that the instability must exist for roughly 100 times the growth time for the instability to enter the non-linear growth phase. We explore the root of this correlation and the implications for linear device design. △ Less

Submitted 14 June, 2017; originally announced June 2017.

Comments: 11 pages, 5 figures

Journal ref: J Leddy and B Dudson 2017 Plasma Phys. Control. Fusion 59 125011

Testing nonlocal models of electron thermal conduction for magnetic and inertial confinement fusion applications

Authors: Jonathan Peter Brodrick, Robert J. Kingham, Michael M. Marinak, Mehul V. Patel, Alex V. Chankin, John Omotani, Maxim Umansky, Dario Del Sorbo, Ben Dudson, Joseph Thomas Parker, Gary D. Kerbel, Mark Sherlock, Christopher P Ridgers

Abstract: Three models for nonlocal electron thermal transport are here compared against Vlasov-Fokker-Planck (VFP) codes to assess their accuracy in situations relevant to both inertial fusion hohlraums and tokamak scrape-off layers. The models tested are (i) a moment-based approach using an eigenvector integral closure (EIC) originally developed by Ji, Held and Sovinec; (ii) the non-Fourier Landau-fluid (… ▽ More Three models for nonlocal electron thermal transport are here compared against Vlasov-Fokker-Planck (VFP) codes to assess their accuracy in situations relevant to both inertial fusion hohlraums and tokamak scrape-off layers. The models tested are (i) a moment-based approach using an eigenvector integral closure (EIC) originally developed by Ji, Held and Sovinec; (ii) the non-Fourier Landau-fluid (NFLF) model of Dimits, Joseph and Umansky; and (iii) Schurtz, Nicolaï and Busquet's multigroup diffusion model (SNB). We find that while the EIC and NFLF models accurately predict the damping rate of a small-amplitude temperature perturbation (within 10% at moderate collisionalities), they overestimate the peak heat flow by as much as 35% and do not predict preheat in the more relevant case where there is a large temperature difference. The SNB model, however, agrees better with VFP results for the latter problem if care is taken with the definition of the mean free path. Additionally, we present for the first time a comparison of the SNB model against a VFP code for a hohlraum-relevant problem with inhomogeneous ionisation and show that the model overestimates the heat flow in the helium gas-fill by a factor of ~2 despite predicting the peak heat flux to within 16%. △ Less

Submitted 6 September, 2017; v1 submitted 28 April, 2017; originally announced April 2017.

Comments: 15 pages, 8 figures. This article may be downloaded for personal use only. Any other use requires prior permission of the author and AIP Publishing. The following article appeared in Brodrick et al., Physics of Plasmas 24, 092309 (2017) and may be found at http://dx.doi.org/10.1063/1.5001079

Journal ref: Brodrick et al., Physics of Plasmas 24, 092309 (2017)

Overview of recent physics results from MAST

Authors: A Kirk, J Adamek, RJ Akers, S Allan, L Appel, F Arese Lucini, M Barnes, T Barrett, N Ben Ayed, W Boeglin, J Bradley, P K Browning, J Brunner, P Cahyna, M Carr, F Casson, M Cecconello, C Challis, IT Chapman, S Chapman, S Conroy, N Conway, WA Cooper, M Cox, N Crocker , et al. (138 additional authors not shown)

Abstract: New results from MAST are presented that focus on validating models in order to extrapolate to future devices. Measurements during start-up experiments have shown how the bulk ion temperature rise scales with the square of the reconnecting field. During the current ramp up models are not able to correctly predict the current diffusion. Experiments have been performed looking at edge and core turbu… ▽ More New results from MAST are presented that focus on validating models in order to extrapolate to future devices. Measurements during start-up experiments have shown how the bulk ion temperature rise scales with the square of the reconnecting field. During the current ramp up models are not able to correctly predict the current diffusion. Experiments have been performed looking at edge and core turbulence. At the edge detailed studies have revealed how filament characteristic are responsible for determining the near and far SOL density profiles. In the core the intrinsic rotation and electron scale turbulence have been measured. The role that the fast ion gradient has on redistributing fast ions through fishbone modes has led to a redesign of the neutral beam injector on MAST Upgrade. In H-mode the turbulence at the pedestal top has been shown to be consistent with being due to electron temperature gradient modes. A reconnection process appears to occur during ELMs and the number of filaments released determines the power profile at the divertor. Resonant magnetic perturbations can mitigate ELMs provided the edge peeling response is maximised and the core kink response minimised. The mitigation of intrinsic error fields with toroidal mode number n>1 has been shown to be important for plasma performance. △ Less

Submitted 18 November, 2016; originally announced November 2016.

Comments: 34 pages, 10 figures. This is an author-created, un-copyedited version of an article submitted for publication in Nuclear Fusion. IoP Publishing Ltd is not responsible for any errors or omissions in this version of the manuscript or any version derived from it

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

Authors: Brendan Shanahan, Peter Hill, Ben Dudson

Abstract: 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 nonaxi… ▽ More 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. △ Less

Submitted 21 September, 2016; originally announced September 2016.

Hermes: Global plasma edge fluid turbulence simulations

Authors: Ben Dudson, Jarrod Leddy

Abstract: The transport of heat and particles in the relatively collisional edge regions of magnetically confined plasmas is a scientifically challenging and technologically important problem. Understanding and predicting this transport requires the self-consistent evolution of plasma fluctuations, global profiles and flows, but the numerical tools capable of doing this in realistic (diverted) geometry are… ▽ More The transport of heat and particles in the relatively collisional edge regions of magnetically confined plasmas is a scientifically challenging and technologically important problem. Understanding and predicting this transport requires the self-consistent evolution of plasma fluctuations, global profiles and flows, but the numerical tools capable of doing this in realistic (diverted) geometry are only now being developed. Here a 5-field reduced 2-fluid plasma model for the study of instabilities and turbulence in magnetised plasmas is presented, built on the BOUT++ framework. This cold ion model allows the evolution of global profiles, electric fields and flows on transport timescales, with flux-driven cross-field transport determined self-consistently by electromagnetic turbulence. Developments in the model formulation and numerical implementation are described, and simulations are performed in poloidally limited and diverted tokamak configurations. △ Less

Submitted 12 September, 2016; originally announced September 2016.

Comments: Joint Varenna-Lausanne international workshop on Theory of Fusion Plasmas

Dirichlet boundary conditions for arbitrary-shaped boundaries in stellarator-like magnetic fields for the Flux-Coordinate Independent method

Authors: Peter Hill, Brendan Shanahan, Ben Dudson

Abstract: We present a technique for handling Dirichlet boundary conditions with the Flux Coordinate Independent (FCI) parallel derivative operator with arbitrary-shaped material geometry in general 3D magnetic fields. The FCI method constructs a finite difference scheme for $\nabla_\Vert$ by following field lines between planes and interpolating within planes. Doing so removes the need for field-aligned co… ▽ More We present a technique for handling Dirichlet boundary conditions with the Flux Coordinate Independent (FCI) parallel derivative operator with arbitrary-shaped material geometry in general 3D magnetic fields. The FCI method constructs a finite difference scheme for $\nabla_\Vert$ by following field lines between planes and interpolating within planes. Doing so removes the need for field-aligned coordinate systems that suffer from singularities in the metric tensor at null points in the magnetic field (or equivalently, when $q \to \infty$). One cost of this method is that as the field lines are not on the mesh, they may leave the domain at any point between neighbouring planes, complicating the application of boundary conditions. The Leg Value Fill (LVF) boundary condition scheme presented here involves an extrapolation/interpolation of the boundary value onto the field line end point. The usual finite difference scheme can then be used unmodified. We implement the LVF scheme in BOUT++ and use the Method of Manufactured Solutions to verify the implementation in a rectangular domain, and show that the error scaling of the finite difference scheme isn't modified. The use of LVF for arbitrary wall geometry is outlined. We also demonstrate the feasibility of using the FCI approach in non-axisymmetric configurations for a simple diffusion model in a "straight stellarator" magnetic field. A Gaussian blob diffuses along the field lines, tracing out flux surfaces. Dirichlet boundary conditions impose a last closed flux surface (LCFS) that confines the density. Including a poloidal limiter moves the LCFS to a smaller radius. The expected scaling of the numerical perpendicular diffusion, which is a consequence of the FCI method, in stellarator-like geometry is recovered. A novel technique for increasing the parallel resolution during post-processing is also described. △ Less

Submitted 8 August, 2016; originally announced August 2016.

Blob dynamics in TORPEX poloidal null configurations

Authors: Brendan Shanahan, Ben Dudson

Abstract: 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 wi… ▽ More 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. △ Less

Submitted 2 August, 2016; v1 submitted 3 May, 2016; originally announced May 2016.

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

Authors: Jarrod Leddy, Ben Dudson, Michele Romanelli, Brendan Shanahan, Nick Walkden

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" (ie. form closed surfaces of constant poloidal flux) in the core of a tokamak. Towards, the edge, however, the field lines begin to intersect physical surfaces, lea… ▽ More 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" (ie. form closed surfaces of constant poloidal flux) in the core of a tokamak. Towards, the edge, however, the field lines begin to 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. △ Less

Submitted 20 April, 2016; originally announced April 2016.

Comments: 11 pages, 12 figures, submitted to Computer Physics Communications

Journal ref: Computer Physics Communications, Volume 212, 2017, Pages 59-68

Verification of BOUT++ by the Method of Manufactured Solutions

Authors: Ben Dudson, Jens Madsen, John Omotani, Peter Hill, Luke Easy, Michael Løiten

Abstract: 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… ▽ More 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. △ Less

Submitted 22 February, 2016; originally announced February 2016.

Investigation of the Effect of Resistivity on Scrape Off Layer Filaments using Three Dimensional Simulations

Authors: Luke Easy, Fulvio Militello, John Omotani, Nick Walkden, Benjamin Dudson

Abstract: The propagation of filaments in the Scrape Off Layer (SOL) of tokamaks largely determine the plasma profiles in the region. In a conduction limited SOL, parallel temperature gradients are expected, such that the resistance to parallel currents is greater at the target than further upstream. Since the perpendicular motion of an isolated filament is largely determined by balance of currents that flo… ▽ More The propagation of filaments in the Scrape Off Layer (SOL) of tokamaks largely determine the plasma profiles in the region. In a conduction limited SOL, parallel temperature gradients are expected, such that the resistance to parallel currents is greater at the target than further upstream. Since the perpendicular motion of an isolated filament is largely determined by balance of currents that flow through it, this may be expected to affect filament transport. 3D simulations have thus been used to study the influence of enhanced parallel resistivity on the dynamics of filaments. Filaments with the smallest perpendicular length scales, which were inertially limited at low resistivity (meaning that polarization rather than parallel currents determine their radial velocities), were unaffected by resistivity. For larger filaments, faster velocities were produced at higher resistivities, due to two mechanisms. Firstly parallel currents were reduced and polarization currents were enhanced, meaning that the inertial regime extended to larger filaments, and secondly a potential difference formed along the parallel direction so that higher potentials were produced in the region of the filament for the same amount of current to flow into the sheath. These results indicate that broader SOL profiles could be produced at higher resistivities. △ Less

Submitted 17 August, 2015; originally announced August 2015.

Comments: 14 Pages, 15 Figures

Numerical investigation of isolated filament motion in a realistic tokamak geometry

Authors: N. R. Walkden, B. D. Dudson, L. Easy, G. Fishpool, J. T. Omotani

Abstract: 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 th… ▽ More 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. △ Less

Submitted 17 August, 2015; originally announced August 2015.

Comments: Submitted to Nuclear Fusion

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

Authors: Jarrod Leddy, Ben Dudson, Michele Romanelli

Abstract: 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… ▽ More 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. △ Less

Submitted 11 June, 2015; v1 submitted 10 June, 2015; originally announced June 2015.

Comments: submitted to PPCF

Journal ref: 2015 Plasma Phys. Control. Fusion 57 125016

Profile measurements in the plasma edge of MAST using a ball pen probe

Authors: N. R. Walkden, J. Adamek, S. Allan, B. D. Dudson, S. Elmore, G. Fishpool, J. Harrison, A. Kirk, M. Komm

Abstract: 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 (MAST). 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 elec… ▽ More 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 (MAST). 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 ~1kV/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 ExB 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. △ Less

Submitted 26 November, 2014; originally announced November 2014.

Comments: 9 pages, 10 figures, submitted to Review of Scientific Instruments

3D Simulations of Plasma Filaments in the Scrape Off Layer: A Comparison with Models of Reduced Dimensionality

Authors: Luke Easy, Fulvio Militello, John Omotani, Benjamin Dudson, Eva Havlíčková, Patrick Tamain, Volker Naulin, Anders H. Nielsen

Abstract: This paper presents simulations of isolated 3D filaments in a slab geometry obtained using a 3D reduced fluid code. 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 import… ▽ More This paper presents simulations of isolated 3D filaments in a slab geometry obtained using a 3D reduced fluid code. 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 polarisation 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 acceleration of all the filaments studied. In contrast, a more satisfactory comparison with the sheath dissipation closure was obtained, even in the presence of significant parallel gradients, where the closure is no longer valid. Specifically it captured the contrasting dynamics of filaments with different perpen- dicular sizes that were observed in the 3D simulations which the vorticity advection closure failed to replicate. However, neither closure successfully replicated the Boltzmann spinning effects and associated poloidal drift of the blob that was observed in the 3D simulations. Although the sheath dissipation closure was concluded to be more successful in replicating the 3D dynamics, it is emphasised that the vorticity closure may still be relevant for situations where the parallel current is inhibited from closing through the sheath due to effects such as increased resistivity near the targets. △ Less

Submitted 8 October, 2014; originally announced October 2014.

Comments: 16 pages, 21 Figures

BOUT++: Recent and current developments

Authors: B. D. Dudson, A. Allen, G. Breyiannis, E. Brugger, J. Buchanan, L. Easy, S. Farley, I. Joseph, M. Kim, A. D. McGann, J. T. Omotani, M. V. Umansky, N. R. Walkden, T. Xia, X. Q. Xu

Abstract: 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 str… ▽ More 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. △ Less

Submitted 30 May, 2014; originally announced May 2014.

Comments: 16 pages, 8 figures, submitted to Journal of Plasma Physics

Non-local parallel transport in BOUT++

Authors: J. T. Omotani, B. D. Dudson, E. Havlickova, M. Umansky

Abstract: 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 de… ▽ More 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. △ Less

Submitted 23 May, 2014; originally announced May 2014.

Characterization of 3D filament dynamics in a MAST SOL fluxtube geometry

Authors: Nicholas Walkden, Benjamin Dudson, Geoff Fishpool

Abstract: 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 th… ▽ More 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. The Boltzmann response causes the filament to self-organise and spin on an axis. 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 anti-symmetric state associated with the interchange mechanism, but high in the Boltzmann regime. △ Less

Submitted 19 July, 2013; originally announced July 2013.

Comments: 28 pages, 9 figures

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

Authors: John Omotani, Ben Dudson

Abstract: 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… ▽ More 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. △ Less

Submitted 15 July, 2013; v1 submitted 4 February, 2013; originally announced February 2013.

Comments: Accepted version. Corrected normalization in (6). This is an author-created, un-copyedited version of an article accepted 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. The Version of Record is available online at doi:10.1088/0741-3335/55/5/055009

Journal ref: Plasma Physics and Controlled Fusion 55 055009, 2013

Improved Nonlinear Solvers in BOUT++

Authors: Ben Dudson, Sean Farley, Lois Curfman McInnes

Abstract: 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… ▽ More 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. △ Less

Submitted 10 September, 2012; originally announced September 2012.

Comments: 15 pages, 4 figures. Submitted to Comp. Phys. Comm

Energy dynamics in a simulation of LAPD turbulence

Authors: B. Friedman, T. A. Carter, M. V. Umansky, D. Schaffner, B. Dudson

Abstract: Energy dynamics calculations in a 3D fluid simulation of drift wave turbulence in the linear Large Plasma Device (LAPD) [W. Gekelman et al., Rev. Sci. Inst. 62, 2875 (1991)] illuminate processes that drive and dissipate the turbulence. These calculations reveal that a nonlinear instability dominates the injection of energy into the turbulence by overtaking the linear drift wave instability that do… ▽ More Energy dynamics calculations in a 3D fluid simulation of drift wave turbulence in the linear Large Plasma Device (LAPD) [W. Gekelman et al., Rev. Sci. Inst. 62, 2875 (1991)] illuminate processes that drive and dissipate the turbulence. These calculations reveal that a nonlinear instability dominates the injection of energy into the turbulence by overtaking the linear drift wave instability that dominates when fluctuations about the equilibrium are small. The nonlinear instability drives flute-like ($k_\parallel = 0$) density fluctuations using free energy from the background density gradient. Through nonlinear axial wavenumber transfer to $k_\parallel \ne 0$ fluctuations, the nonlinear instability accesses the adiabatic response, which provides the requisite energy transfer channel from density to potential fluctuations as well as the phase shift that causes instability. The turbulence characteristics in the simulations agree remarkably well with experiment. When the nonlinear instability is artificially removed from the system through suppressing $k_\parallel=0$ modes, the turbulence develops a coherent frequency spectrum which is inconsistent with experimental data. △ Less

Submitted 3 January, 2013; v1 submitted 10 May, 2012; originally announced May 2012.

Journal ref: B. Friedman, T. A. Carter, M. V. Umansky, D. Schaffner, and B. Dudson, Phys. Plasmas 19, 102307 (2012)

Simulation of Edge Localised Modes using BOUT++

Authors: B. D. Dudson, X. Q. Xu, M. V. Umansky, H. R. Wilson, P. B. Snyder

Abstract: 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 visc… ▽ More 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. △ Less

Submitted 26 August, 2010; originally announced August 2010.

Comments: 14 pages, 14 figures. Theory of Fusion Plasmas workshop, Varenna 2010. Submitted to Plas. Phys. Control. Fusion

BOUT++: a framework for parallel plasma fluid simulations

Authors: B. D. Dudson, M. V. Umansky, X. Q. Xu, P. B. Snyder, H. R. Wilson

Abstract: 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… ▽ More 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. △ Less

Submitted 5 November, 2008; v1 submitted 31 October, 2008; originally announced October 2008.

Comments: Submitted to Computer Physics Communications. Revised to reduce page count. 18 pages, 16 figures

Journal ref: Computer Physics Communications 180 (2009), pp. 1467-1480