Source author record

M. J. Hole

M. J. Hole 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

22works
5topics
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

22 published item(s)

preprint2020arXiv

Sandpile modelling of dual location fuelling in fusion plasmas

We modify the Chapman sandpile model (Chapman \textit{et al} \textit{Physical Review Letters} 86, 2814 (2001)) to form comparisons with pellet pacing, which is used to reduce or eliminate ELMs in a fusion plasma. We employ a variation of that model in which a pedestal with feedback is introduced (Bowie and Hole \textit{Phys. Plasmas} 25, 012511 (2018)), which we further modify to provide for dual fuelling - sand is added both at the centre of the sandpile, and near the edge. We observe that when the additional sand is added at the top of the pedestal, MLEs are largely suppressed. While this suppression comes at a cost by way of reduction in total energy confinement, that reduction is lower than the reduction in MLE size. The trade-off between MLE suppression and reduction in energy confinement depends not only on the amount of extra sand, but also on its precise location relative to the top of the pedestal. We suggest that the approach of constant dual fuelling may be equally applicable to plasmas, and may suggest a strategy for ELM suppression in fusion plasmas. We observe that when the proposed amount of extra sand is added in 'pellets', using frequencies and amounts based on those proposed for ELM suppression for ITER, MLEs are similarly suppressed, although MLEs are not significantly suppressed when the pellet rate does not substantially exceed the MLE frequency. This suggests that pellet injection at the top of the pedestal at small pellet size and high frequency may represent a reasonable physical proxy for our proposed scheme. However, our results suggest that it is not the synchronisation of pellets to ELM frequencies which is the key factor for ELM suppression in this regime, but rather the introduction of additional fuelling at the top of the pedestal.

preprint2020arXiv

Stepped pressure equilibrium with relaxed flow and applications in reversed-field pinch plasmas

The Multi-region Relaxed MHD (MRxMHD) has been successful in the construction of equilibria in three-dimensional (3D) configurations. In MRxMHD, the plasma is sliced into sub-volumes separated by ideal interfaces, each undergoing relaxation, allowing the formation of islands and chaos. The resulting equilibrium has a stepped pressure profile across sub-volumes. The Stepped Pressure Equilibrium Code (SPEC) [S.R. Hudson et al., Phys. Plasmas 19, 112502 (2012)] was developed to calculate MRxMHD equilibria numerically. In this work, we have extended the SPEC code to compute MRxMHD equilibria with field-aligned flow and rotation, following the theoretical development to incorporate cross-helicity and angular momentum constraints. The code has been verified for convergence and compared to a Grad-Shafranov solver in 2D. We apply our new tool to study the flow profile change before and after the sawtooth crash of a reversed-field pinch discharge, in which data of the parallel flow is available. We find the promising result that under the constraints of cross-helicity and angular momentum, the parallel flow profile in post-crash SPEC equilibrium is flat in the plasma core and the amplitude of the flow matches experimental observations. Finally, we provide an example equilibrium with a 3D helical field structure as the favoured lower energy state. This will be the first 3D numerical equilibrium in which the flow effects are self-consistently calculated.

preprint2020arXiv

The impact of anisotropy on ITER scenarios

We report on the impact of anisotropy to tokamak plasma configuration and stability. Our focus is on analysis of the impact of anisotropy on ITER pre-fusion power operation 5~MA, $B=1.8$~T ICRH scenarios. To model ITER scenarios remapping tools are developed to distinguish the impact of pressure anisotropy from the change in magnetic geometry caused by an anisotropy-modified current profile. The remappings iterate the anisotropy-modified current profile to produce the same $q$ profile with matched thermal energy. The analysis is a step toward equilibria that are kinetically self-consistent for a prescribed scenario. We find characteristic detachment of flux surfaces from pressure surfaces, and an outboard (inboard) shift of peak density for $T_{\parallel}>T_\perp$ ( $T_{\parallel}<T_\perp$). Differences in the poloidal current profile are evident, albeit not as pronounced as for the spherical tokamak. We find that the incompressional continuum is largely unchanged in the presence of anisotropy, and the mode structure of gap modes is largely unchanged. The compressional branch however exhibits significant differences in the continuum. We report on the implication of these modifications.

preprint2020arXiv

Time-dependent relaxed magnetohydrodynamics -- inclusion of cross helicity constraint using phase-space action

A new formulation of time-dependent Relaxed Magnetohydrodynamics (RxMHD) is derived variationally from Hamilton's Action Principle using microscopic conservation of mass, and macroscopic conservation of total magnetic helicity, cross helicity and entropy, as the only constraints on variations of density, pressure, fluid velocity, and magnetic vector potential over a relaxation domain. A novel phase-space version of the MHD Lagrangian is derived, which gives Euler--Lagrange equations consistent with previous work on exact ideal and relaxed MHD equilibria with flow, but generalizes the relaxation concept from statics to dynamics. The application of the new dynamical formalism is illustrated for short-wavelength linear waves, and the interface connection conditions for Multiregion Relaxed MHD (MRxMHD) are derived. The issue of whether $\vec{E} + \vec{u}\times\vec{B} = 0$ should be a constraint is discussed.

preprint2019arXiv

Sandpile modelling of pellet pacing in fusion plasmas

Sandpile models have been used to provide simple phenomenological models without incorporating the detailed features of a fully featured model. The Chapman sandpile model Chapman et al Physical Review Letters 86, 2814 (2001) has been used as an analogue for the behaviour of a plasma edge, with mass loss events being used as analogues for ELMs. In this work we modify the Chapman sandpile model by providing for both increased and intermittent driving. We show that the behaviour of the sandpile, when continuously fuelled at very high driving, can be determined analytically by a simple algorithm. We observe that the size of the largest avalanches is better reduced by increasing constant driving than by the intermittent introduction of `pellets' of sand. Using the sandpile model as a reduced model of ELMing behaviour, we conject that ELM control in a fusion plasma may similarly prove more effective with increased total fuelling than with pellet addition.

preprint2015arXiv

Calculation of continuum damping of Alfvén eigenmodes in 2D and 3D cases

In ideal MHD, shear Alfvén eigenmodes may experience dissipationless damping due to resonant interaction with the shear Alfvén continuum. This continuum damping can make a significant contribution to the overall growth/decay rate of shear Alfvén eigenmodes, with consequent implications for fast ion transport. One method for calculating continuum damping is to solve the MHD eigenvalue problem over a suitable contour in the complex plane, thereby satisfying the causality condition. Such an approach can be implemented in three-dimensional ideal MHD codes which use the Galerkin method. Analytic functions can be fitted to numerical data for equilibrium quantities in order to determine the value of these quantities along the complex contour. This approach requires less resolution than the established technique of calculating damping as resistivity vanishes and is thus more computationally efficient. The complex contour method has been applied to the three-dimensional finite element ideal MHD code CKA . In this paper we discuss the application of the complex contour technique to calculate the continuum damping of global modes in tokamak as well as torsatron, W7X and H1-NF stellarator cases. To the authors' knowledge these stellarator calculations represent the first calculation of continuum damping for eigenmodes in fully three-dimensional equilibria. The continuum damping of global modes in W7X and H1-NF stellarator configurations investigated is found to depend sensitively on coupling to numerous poloidal and toroidal harmonics.

preprint2014arXiv

Magnetohydrodynamic normal mode analysis of plasma with equilibrium pressure anisotropy

In this work, we generalise linear magnetohydrodynamic (MHD) stability theory to include equilibrium pressure anisotropy in the fluid part of the analysis. A novel 'single-adiabatic' (SA) fluid closure is presented which is complementary to the usual 'double-adiabatic' (CGL) model and has the advantage of naturally reproducing exactly the MHD spectrum in the isotropic limit. As with MHD and CGL, the SA model neglects the anisotropic perturbed pressure and thus loses non-local fast-particle stabilisation present in the kinetic approach. Another interesting aspect of this new approach is that the stabilising terms appear naturally as separate viscous corrections leaving the isotropic SA closure unchanged. After verifying the self-consistency of the SA model, we re-derive the projected linear MHD set of equations required for stability analysis of tokamaks in the MISHKA code. The cylindrical wave equation is derived analytically as done previously in the spectral theory of MHD and clear predictions are made for the modification to fast-magnetosonic and slow ion sound speeds due to equilibrium anisotropy.

preprint2014arXiv

Multi-region relaxed magnetohydrodynamics with flow

We present an extension of the multi-region relaxed magnetohydrodynamics (MRxMHD) equilibrium model that includes plasma flow. This new model is a generalization of Woltjer's model of relaxed magnetohydrodynamics equilibria with flow. We prove that as the number of plasma regions becomes infinite our extension of MRxMHD reduces to ideal MHD with flow. We also prove that some solutions to MRxMHD with flow are not time-independent in the laboratory frame, and instead have 3D structure which rotates in the toroidal direction with fixed angular velocity. This capability gives MRxMHD potential application to describing rotating 3D MHD structures such as 'snakes' and long-lived modes.

preprint2014arXiv

Recent Developments in Bayesian Inference of Tokamak Plasma Equilibria and High-Dimensional Stochastic Quadratures

We present recent results and technical breakthroughs for the Bayesian inference of tokamak equilibria using force-balance as a prior constraint. Issues surrounding model parameter representation and posterior analysis are discussed and addressed. These points motivate the recent advancements embodied in the Bayesian Equilibrium Analysis and Simulation Tool (BEAST) software being presently utilised to study equilibria on the Mega-Ampere Spherical Tokamak (MAST) experiment in the UK (von Nessi \emph{et. al.} 2012 \emph{J. Phys. A} \textbf{46} 185501). State-of-the-art results of using BEAST to study MAST equilbria are reviewed, with recent code advancements being systematically presented though out the manuscript.

preprint2013arXiv

A minimally constrained model of self-organized helical states in reversed-field pinches

We show that the self-organized single-helical-axis (SHAx) and double-helical-axis (DAx) states in reversed field pinches can be reproduced in a minimally constrained equilibrium model using only five parameters. This is a significant reduction on previous representations of the SHAx which have required an infinite number of constraints. The DAx state, which has a non-trivial topology, has not been previously reproduced using an equilibrium model that preserves this topological structure. We show that both states are a consequence of transport barrier formation in the plasma core, in agreement with experimental results.

preprint2013arXiv

A Unified Method for Inference of Tokamak Equilibria and Validation of Force-Balance Models Based on Bayesian Analysis

A new method, based on Bayesian analysis, is presented which unifies the inference of plasma equilibria parameters in a Tokamak with the ability to quantify differences between inferred equilibria and Grad-Shafranov force-balance solutions. At the heart of this technique is the new method of observation splitting, which allows multiple forward models to be associated with a single diagnostic observation. This new idea subsequently provides a means by which the the space of GS solutions can be efficiently characterised via a prior distribution. Moreover, by folding force-balance directly into one set of forward models and utilising simple Biot-Savart responses in another, the Bayesian inference of the plasma parameters itself produces an evidence (a normalisation constant of the inferred posterior distribution) which is sensitive to the relative consistency between both sets of models. This evidence can then be used to help determine the relative accuracy of the tested force-balance model across several discharges/times. These ideas have been implemented in a code called BEAST (Bayesian Equilibrium Analysis and Simulation Tool), which uses a special implementation of Skilling's nested sampling algorithm [Skilling, Bayesian Analysis 1(4), 833--859 (2006)] to perform sampling and evidence calculations on high-dimensional, non-Gaussian posteriors. Initial BEAST equilibrium inference results are presented for two high-performance MAST discharges.

preprint2013arXiv

First evidence of Alfven wave activity in KSTAR plasmas

We report on first evidence of wave activity during neutral beam heating in KSTAR plasmas: 40 kHz magnetic fluctuations with a toroidal mode number of n=1. Our analysis suggests this a beta-induced Alfven eigenmode resonant with the q=1 surface. A kinetic analysis, when coupled with electron temperature measurements from electron cyclotron emission and ion/electron temperature ratios from crystallography, enables calculation of the frequency evolution, which is in agreement with observations. Complementary detailed MHD modelling of the magnetic configuration and wave modes supports the BAE mode conclusion, by locating an n=1 mode separated from the continuum in the core region. Finally, we have computed the threshold to marginal stability for a range of ion temperature profiles. These suggest the BAE can be driven unstable by energetic ions when the ion temperature radial gradient is sufficiently large. Our findings suggest that mode existence could be used as a form of inference for temperature profile consistency in the radial interval of the mode, thereby extending the tools of MHD spectroscopy.

preprint2013arXiv

Using Bayesian Analysis and Gaussian Processes to Infer Electron Temperature and Density Profiles on the MAST Experiment

A unified, Bayesian inference of midplane electron temperature and density profiles using both Thompson scattering (TS) and interferometric data is presented. Beyond the Bayesian nature of the analysis, novel features of the inference are the use of a Gaussian process prior to infer a mollification length-scale of inferred profiles and the use of Gauss-Laguerre quadratures to directly calculate the depolarisation term associated with the TS forward model. Results are presented from an application of the method to data from the high resolution TS system on the Mega-Ampere Spherical Tokamak, along with a comparison to profiles coming from the standard analysis carried out on that system.

preprint2012arXiv

Computation of multi-region relaxed magnetohydrodynamic equilibria

We describe the construction of stepped-pressure equilibria as extrema of a multi-region, relaxed magnetohydrodynamic (MHD) energy functional that combines elements of ideal MHD and Taylor relaxation, and which we call MRXMHD. The model is compatible with Hamiltonian chaos theory and allows the three-dimensional MHD equilibrium problem to be formulated in a well-posed manner suitable for computation. The energy-functional is discretized using a mixed finite-element, Fourier representation for the magnetic vector potential and the equilibrium geometry; and numerical solutions are constructed using the stepped-pressure equilibrium code, SPEC. Convergence studies with respect to radial and Fourier resolution are presented.

preprint2012arXiv

Evidence Cross-Validation and Bayesian Inference of MAST Plasma Equilibria

In this paper, current profiles for plasma discharges on the Mega-Ampere Spherical Tokamak (MAST) are directly calculated from pickup coil, flux loop and Motional-Stark Effect (MSE) observations via methods based in the statistical theory of Bayesian analysis. By representing toroidal plasma current as a series of axisymmetric current beams with rectangular cross-section and inferring the current for each one of these beams, flux-surface geometry and q-profiles are subsequently calculated by elementary application of Biot-Savart's law. The use of this plasma model in the context of Bayesian analysis was pioneered by Svensson and Werner on the Joint-European Tokamak (JET) [J. Svensson and A. Werner. Current tomography for axisymmetric plasmas. {\em Plasma Physics and Controlled Fusion}, 50(8):085002, 2008]. In this framework, linear forward models are used to generate diagnostic predictions, and the probability distribution for the currents in the collection of plasma beams was subsequently calculated directly via application of Bayes' formula. In this work, we introduce a new diagnostic technique to identify and remove outlier observations associated with diagnostics falling out of calibration or suffering from an unidentified malfunction. These modifications enable good agreement between Bayesian inference of the last closed flux-surface (LCFS) with other corroborating data, such as such as that from force balance considerations using EFIT++ [L. Appel et al., Proc. 33rd EPS Conf., Rome, Italy, 2006]. In addition, this analysis also yields errors on the plasma current profile and flux-surface geometry, as well as directly predicting the Shafranov shift of the plasma core.

preprint2012arXiv

Wave modelling in a cylindrical non-uniform helicon discharge

A radio frequency (RF) field solver based on Maxwell's equations and a cold plasma dielectric tensor is em- ployed to describe wave phenomena observed in a cylindrical non-uniform helicon discharge. The experiment is carried out on a recently built linear plasma-material interaction machine: the MAGnetized Plasma In- teraction Experiment (MAGPIE) [B. D. Blackwell, J. F. Caneses, C. Samuell, J. Wach, J. Howard, and C. S. Corr, submitted on 25 March 2012 to Plasma Sources Science and Technology], in which both plasma density and static magnetic field are functions of axial position. The field strength increases by a factor of 15 from source to target plate, and plasma density and electron temperature are radially non-uniform. With an enhancement factor of 9.5 to the electron-ion Coulomb collision frequency, 12% reduction in the antenna radius, and the same other conditions as employed in the experiment, the solver produces axial and radial profiles of wave amplitude and phase that are consistent with measurements. Ion-acoustic turbulence, which can happen if electron drift velocity exceeds the speed of sound in magnetized plasmas, may account for the factor of 9.5 used to match simulated results with experimental data. To overcome the single m vacuum solu- tion limitations of the RF solver, which can only compute the glass response to the same mode number of the antenna, we have adjusted the antenna radius to match the wave field strength in the plasma.(not finished because of the limited number of characters, please see the full paper)

preprint2011arXiv

A flowing plasma model to describe drift waves in a cylindrical helicon discharge

A two-fluid model developed originally to describe wave oscillations in the vacuum arc centrifuge, a cylindrical, rapidly rotating, low temperature and confined plasma column, is applied to interpret plasma oscillations in a RF generated linear magnetised plasma (WOMBAT), with similar density and field strength. Compared to typical centrifuge plasmas, WOMBAT plasmas have slower normalised rotation frequency, lower temperature and lower axial velocity. Despite these differences, the two-fluid model provides a consistent description of the WOMBAT plasma configuration and yields qualitative agreement between measured and predicted wave oscillation frequencies with axial field strength. In addition, the radial profile of the density perturbation predicted by this model is consistent with the data. Parameter scans show that the dispersion curve is sensitive to the axial field strength and the electron temperature, and the dependence of oscillation frequency with electron temperature matches the experiment. These results consolidate earlier claims that the density and floating potential oscillations are a resistive drift mode, driven by the density gradient. To our knowledge, this is the first detailed physics model of flowing plasmas in the diffusion region away from the RF source. Possible extensions to the model, including temperature non-uniformity and magnetic field oscillations, are also discussed.

preprint2011arXiv

Nonaxisymmetric, multi-region relaxed magnetohydrodynamic equilibrium solutions

We describe a magnetohydrodynamic (MHD) constrained energy functional for equilibrium calculations that combines the topological constraints of ideal MHD with elements of Taylor relaxation. Extremizing states allow for partially chaotic magnetic fields and non-trivial pressure profiles supported by a discrete set of ideal interfaces with irrational rotational transforms. Numerical solutions are computed using the Stepped Pressure Equilibrium Code, SPEC, and benchmarks and convergence calculations are presented.

preprint2010arXiv

Model Data Fusion: developing Bayesian inversion to constrain equilibrium and mode structure

Recently, a new probabilistic "data fusion" framework based on Bayesian principles has been developed on JET and W7-AS. The Bayesian analysis framework folds in uncertainties and inter-dependencies in the diagnostic data and signal forward-models, together with prior knowledge of the state of the plasma, to yield predictions of internal magnetic structure. A feature of the framework, known as MINERVA (J. Svensson, A. Werner, Plasma Physics and Controlled Fusion 50, 085022, 2008), is the inference of magnetic flux surfaces without the use of a force balance model. We discuss results from a new project to develop Bayesian inversion tools that aim to (1) distinguish between competing equilibrium theories, which capture different physics, using the MAST spherical tokamak; and (2) test the predictions of MHD theory, particularly mode structure, using the H-1 Heliac.

preprint2009arXiv

Magnetohydrodynamic Stability of Plasmas with Ideal and Relaxed Regions

A unified energy principle approach is presented for analysing the magnetohydrodynamic (MHD) stability of plasmas consisting of multiple ideal and relaxed regions. By choosing an appropriate gauge, we show that the plasma displacement satisfies the same Euler-Lagrange equation in ideal and relaxed regions, except in the neighbourhood of magnetic surfaces. The difference at singular surfaces is analysed in cylindrical geometry: in ideal MHD only Newcomb's [W. A. Newcomb (2006) Ann. Phys., 10, 232] small solutions are allowed, whereas in relaxed MHD only the odd-parity large solution and even-parity small solution are allowed. A procedure for constructing global multi-region solutions in cylindrical geometry is presented. Focussing on the limit where the two interfaces approach each other arbitrarily closely, it is shown that the singular-limit problem encountered previously [M.J. Hole et al. (2006) J. Plasma Phys., 77, 1167] in multi-region relaxed MHD is stabilised if the relaxed-MHD region between the coalescing interfaces is replaced by an ideal-MHD region. We then present a stable (k, pressure) phase space plot, which allows us to determine the form a stable pressure and field profile must take in the region between the interfaces. From this knowledge, we conclude that there exists a class of single interface plasmas that were found stable by Kaiser and Uecker [R. Kaiser et al (2004) Q. Jl Mech. Appl. Math., 57, 1], but are shown to be unstable when the interface is resolved.

preprint2009arXiv

Relaxed MHD states of a multiple region plasma

We calculate the stability of a multiple relaxation region MHD (MRXMHD) plasma, or stepped-Beltrami plasma, using both variational and tearing mode treatments. The configuration studied is a periodic cylinder. In the variational treatment, the problem reduces to an eigenvalue problem for the interface displacements. For the tearing mode treatment, analytic expressions for the tearing mode stability parameter $Δ'$, being the jump in the logarithm in the helical flux across the resonant surface, are found. The stability of these treatments is compared for $m=1$ displacements of an illustrative RFP-like configuration, comprising two distinct plasma regions. For pressure-less configurations, we find the marginal stability conclusions of each treatment to be identical, confirming analytic results in the literature. The tearing mode treatment also resolves ideal MHD unstable solutions for which $Δ' \to \infty$: these correspond to displacement of a resonant interface. Wall stabilisation scans resolve the internal and external ideal kink. Scans with increasing pressure are also performed: these indicate that both variational and tearing mode treatments have the same stability trends with $β$, and show pressure stabilisation in configurations with increasing edge pressure. Combined, our results suggest that MRXMHD configurations which are stable to ideal perturbations plus tearing modes are automatically in a stable state. Such configurations, and their stability properties, are of emerging importance in the quest to find mathematically rigorous solutions of ideal MHD force balance in 3D geometry.