BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

H. R. Wilson

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

ResearcherUnclaimed source record
physics.plasm-phphysics.comp-ph

Catalog footprint

What is connected

11works
2topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

11 published item(s)

preprint2020arXiv

The physics of a small-scale tearing mode in collisionless slab plasmas

Microtearing modes have been widely reported as a tearing parity electron temperature gradient driven plasma instability, which leads to fine scale tearing of the magnetic flux surfaces thereby resulting in reconnection of magnetic field lines and formation of magnetic islands. In slab geometry it has previously been shown that the drive mechanism requires a finite collision frequency. However, we find in linear gyrokinetic simulations that a collisionless fine-scale tearing parity instability exists even at low and zero collision frequency. Detailed studies reveal that these slab modes are also driven by electron temperature gradient but are sensitive to electron finite Larmor radius effects, and have a radial wavenumber much smaller than the binormal wavenumber, which is comparable to the ion Larmor radius. Furthermore, they exist even in the electrostatic limit and electromagnetic effects actually have a stabilising influence on this collisionless tearing mode. An analytic model shows that this collisionless small scale tearing mode is consistent with a tearing parity slab electron temperature gradient (ETG) mode, which can be more unstable than the twisting parity ETG mode that is often studied. This small-scale tearing parity mode can lead to magnetic islands, which, in turn, can influence turbulent transport in magnetised plasmas.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Using a local gyrokinetic code to study global ITG modes in tokamaks

In this paper the global mode structures of linear ion-temperature-gradient (ITG) modes in tokamak plasmas are obtained by combining results from the local gyrokinetic code GS2 with analytical theory. Local gyrokinetic calculations, using GS2, are performed for a range of radial flux surfaces, ${x}$, and ballooning phase angles, ${p}$, to map out the local complex mode frequency, ${Ω_{0}(x,p)=ω_{0}(x,p)+iγ_{0}(x,p)}$ for a single toroidal mode number, ${n}$. Taylor expanding ${Ω_{0}}$ about ${x=0}$, and employing the Fourier-ballooning representation leads to a second order ODE for the amplitude envelope, ${A\left(p\right)}$ , which describes how the local results are combined to form the global mode. We employ the so-called CYCLONE base case for circular Miller equilibrium model. Assuming radially varying profiles of ${a/L_{T}}$ and ${a/L_{n}}$, peaked at ${x=0}$, and with all other equilibrium profiles held constant, ${Ω_{0}(x,p)}$ is found to have a stationary point. The reconstructed global mode sits at the outboard mid-plane of the tokamak, with global growth rate, ${γ\sim}$Max${\left[γ_{0}\right]}$. Including the radial variation of other equilibrium profiles like safety factor and magnetic shear, leads to a mode that peaks away from the outboard mid-plane, with a reduced global growth rate. Finally, the influence of toroidal flow shear has also been investigated through the introduction of a Doppler shift, ${ω_{0} \rightarrow ω_{0} - nΩ_ϕ^{\prime} x}$, where ${Ω_ϕ}$ is the equilibrium toroidal flow, and a prime denotes the radial derivative. The equilibrium profile variations introduce an asymmetry into the global growth rate spectrum with respect to the sign of ${Ω_ϕ^{\prime}}$, such that the maximum growth rate is achieved with non-zero shearing, consistent with recent global gyrokinetic calculations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Explosive Instability and Erupting Flux Tubes in a Magnetised Plasma Atmosphere

The eruption of multiple flux tubes in a magnetised plasma atmosphere is proposed as a mechanism for explosive release of energy in plasmas. Linearly stable isolated flux tubes are shown to be metastable in a box model magnetised atmosphere in which ends of the field lines are embedded in conducting walls. The energy released by destabilising such field lines can be a significant fraction of the gravitational energy stored in the system. This energy can be released in a fast dynamical time.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Structure of Micro-instabilities in Tokamak Plasmas: Stiff Transport or Plasma Eruptions?

Solutions to a model 2D eigenmode equation describing micro-instabilities in tokamak plasmas are presented that demonstrate a sensitivity of the mode structure and stability to plasma profiles. In narrow regions of parameter space, with special plasma profiles, a maximally unstable mode is found that balloons on the outboard side of the tokamak. This corresponds to the conventional picture of a ballooning mode. However, for most profiles this mode cannot exist and instead a more stable mode is found that balloons closer to the top or bottom of the plasma. Good quantitative agreement with a 1D ballooning analysis is found provided the constraints associated with higher order profile effects, often neglected, are taken into account. A sudden transition from this general mode to the more unstable ballooning mode can occur for a critical flow shear, providing a candidate model for why some experiments observe small plasma eruptions (Edge Localised Modes, or ELMs) in place of large Type I ELMs.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Using the local gyrokinetic code, GS2, to investigate global ITG modes in tokamaks. (I) s-$α$ model with profile and flow shear effects

This paper combines results from a local gyrokinetic code with analytical theory to reconstruct the global eigenmode structure of the linearly unstable ion-temperature-gradient (ITG) mode with adiabatic electrons. The simulations presented here employ the s-$α$ tokamak equilibrium model. Local gyrokinetic calculations, using GS2 have been performed over a range of radial surfaces, x, and for ballooning phase angle, p, in the range -$π {\leq} p {\leqπ}$, to map out the complex local mode frequency, ${Ω_0(x, p) = ω_0(x, p) + iγ_0(x, p)}$. Assuming a quadratic radial profile for the drive, namely ${η_i = L_n/L_T}$, (holding constant all other equilibrium profiles such as safety factor, magnetic shear etc.), ${Ω_0(x, p)}$ has a stationary point. The reconstructed global mode then sits on the outboard mid plane of the tokamak plasma, and is known as a conventional or isolated mode, with global growth rate, $γ$ ~ Max[${γ_0(x, p)}$], where ${γ_0(x, p)}$ is the local growth rate. Taking the radial variation in other equilibrium profiles (e.g safety factor q(x)) into account, removes the stationary point in ${Ω_0(x, p)}$ and results in a mode that peaks slightly away from the outboard mid-plane with a reduced global growth rate. Finally, the influence of flow shear has also been investigated through a Doppler shift, ${ω_0 \rightarrow ω_0 + nΩ^{\prime}x}$, where n is the toroidal mode number and ${Ω^{\prime}}$ incorporates the effect of flow shear. The equilibrium profile variation introduces an asymmetry to the growth rate spectrum with respect to the sign of ${Ω^{\prime}}$, consistent with recent global gyrokinetic calculations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Microtearing modes at the top of the pedestal

Microtearing modes (MTMs) are unstable in the shallow gradient region just inside the top of the pedestal in the spherical tokamak experiment MAST, and may play an important role in the pedestal evolution. The linear properties of these instabilities are compared with MTMs deeper inside the core, and further detailed investigations expose the basic drive mechanism. MTMs near the MAST pedestal top are not well described by existing theories. In particular the growth rate of the dominant edge MTM does not peak at a finite collision frequency, as frequently reported for MTMs further into the core. Our study suggests that the edge MTM is driven by a collisionless trapped particle mechanism that is sensitive to magnetic drifts. This drive is enhanced in the outer region of MAST by a high magnetic shear and a high trapped particle fraction. Observations of similar modes on conventional aspect ratio devices suggests this drive mechanism may be somewhat ubiquitous towards the edge of current day and future hot tokamaks.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Collision frequency dependence of polarization current in neoclassical tearing modes

The neoclassical polarization current, generated when a magnetic island propagates through a tokamak plasma, is believed to influence the initial stage of the neoclassical tearing mode evolution. Understanding the strength of its contribution in the relevant plasma collision frequency regimes for future tokamaks such as ITER is crucial for the successful control and/or avoidance of the neoclassical tearing mode. A nonlinear drift kinetic theory is employed to determine the full collision frequency dependence of the neoclassical polarization current in the small island limit, comparable to the trapped ion orbit width. Focusing on the region away from the island separatrix (where a layer with a complex mix of physics processes exists), we evaluate for the first time the variation of the neoclassical ion polarization current in the transition regime between the analytically tractable collisionless and collisional limits. In addition, the island propagation frequency-dependence of the neoclassical polarization current and its contribution to the island evolution is revealed. For a range of propagation frequencies, we find that the neoclassical polarization current is a maximum in the intermediate collision frequency regime analyzed here - a new and unexpected result.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Kinetic instabilities that limit β in the edge of a tokamak plasma: a picture of an H-mode pedestal

Plasma equilibria reconstructed from the Mega-Amp Spherical Tokamak (MAST) have sufficient resolution to capture plasma evolution during the short period between edge-localized modes (ELMs). Immediately after the ELM steep gradients in pressure, P, and density, ne, form pedestals close to the separatrix, and they then expand into the core. Local gyrokinetic analysis over the ELM cycle reveals the dominant microinstabilities at perpendicular wavelengths of the order of the ion Larmor radius. These are kinetic ballooning modes (KBMs) in the pedestal and microtearing modes (MTMs) in the core close to the pedestal top. The evolving growth rate spectra, supported by gyrokinetic analysis using artificial local equilibrium scans, suggest a new physical picture for the formation and arrest of this pedestal.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Towards the construction of a model to describe the inter-ELM evolution of the pedestal on MAST

Pedestal profiles that span the ELM cycle have been obtained and used to test the idea that the pedestal pressure gradient in MAST is limited by the onset of Kinetic Ballooning Modes (KBMs). During the inter-ELM period of a regularly type I ELM-ing discharge on MAST, the pressure pedestal height and width increase together while the pressure gradient increases by only 15 % during the ELM cycle. Stability analyses show that the pedestal region over which infinite-n ballooning modes are unstable also broadens during the ELM cycle. To test the relationship between the width of the region that is unstable to n = \infty ideal magnetohydrodynamic ballooning modes and KBMs the gyrokinetic code, GS2, has been used for microstability analysis of the edge plasma region in MAST. The gyrokinetic simulations find that KBM modes with twisting parity are the dominant microinstabilities in the steep pedestal region, with a transition to more slowly growing tearing parity modes in the shallower pressure gradient core region immediately inside the pedestal top. The region over which KBMs are unstable increases during the ELM cycle, and a good correlation is found between the regions unstable to KBMs and infinite-n ideal ballooning modes.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Simulation of Edge Localised Modes using BOUT++

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

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

BOUT++: a framework for parallel plasma fluid simulations

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

0 citations0 reviewsNo code linkedOpen access