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

R. L. Dewar

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

ResearcherUnclaimed source record
physics.plasm-phnlin.CDmath-phmath.MPastro-ph.SRphysics.ao-phphysics.class-phphysics.data-anphysics.flu-dynphysics.space-ph

Catalog footprint

What is connected

19works
10topics
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

19 published item(s)

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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Plasmoid solutions of the Hahm--Kulsrud--Taylor equilibrium model

The Hahm--Kulsrud (HK) [T. S. Hahm and R. M. Kulsrud, Phys. Fluids {\bf 28}, 2412 (1985)] solutions for a magnetically sheared plasma slab driven by a resonant periodic boundary perturbation illustrate fully shielded (current sheet) and fully reconnected (magnetic island) responses. On the global scale, reconnection involves solving a magnetohydrodynamic (MHD) equilibrium problem. In systems with a continuous symmetry such MHD equilibria are typically found by solving the Grad--Shafranov equation, and in slab geometry the elliptic operator in this equation is the 2-D Laplacian. Thus, assuming appropriate pressure and poloidal current profiles, a conformal mapping method can be used to transform one solution into another with different boundary conditions, giving a continuous sequence of solutions in the form of partially reconnected magnetic islands (plasmoids) separated by Syrovatsky current sheets. The two HK solutions appear as special cases.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Action-gradient-minimizing pseudo-orbits and almost-invariant tori

Transport in near-integrable, but partially chaotic, $1 1/2$ degree-of-freedom Hamiltonian systems is blocked by invariant tori and is reduced at \emph{almost}-invariant tori, both associated with the invariant tori of a neighboring integrable system. "Almost invariant" tori with rational rotation number can be defined using continuous families of periodic \emph{pseudo-orbits} to foliate the surfaces, while irrational-rotation-number tori can be defined by nesting with sequences of such rational tori. Three definitions of "pseudo-orbit," \emph{action-gradient--minimizing} (AGMin), \emph{quadratic-flux-minimizing} (QFMin) and \emph{ghost} orbits, based on variants of Hamilton's Principle, use different strategies to extremize the action as closely as possible. Equivalent Lagrangian (configuration-space action) and Hamiltonian (phase-space action) formulations, and a new approach to visualizing action-minimizing and minimax orbits based on AGMin pseudo-orbits, are presented.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Dressed test particles, oscillation centres and pseudo-orbits

A general semi-analytical method for accurate and efficient numerical calculation of the dielectrically screened ("dressed") potential around a non-relativistic test particle moving in an isotropic, collisionless, unmagnetised plasma is presented. The method requires no approximations and is illustrated using results calculated for two cases taken from the MSc thesis of the first author: test particles with velocities above and below the ion sound speed in plasmas with Maxwellian ions and warm electrons. The idea that the fluctuation spectrum of a plasma can be described as a superposition of the fields around \emph{non-interacting} dressed test particles is an expression of the quasiparticle concept, which has also been expressed in the development of the oscillation-centre and pseudo-orbit formalisms.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Unified Theory of Ghost and Quadratic-Flux-Minimizing Surfaces

A generalized Hamiltonian definition of ghost surfaces (surfaces defined by an action-gradient flow) is given and specialized to the usual Lagrangian definition. Numerical calculations show uncorrected quadratic-flux-minimizing (QFMin) and Lagrangian ghost surfaces give very similar results for a chaotic magnetic field weakly perturbed from an integrable case in action-angle coordinates, described by $L = L_0 + εL_1$, where $L_0(\dotθ)$ (with $\dotθ$ denoting $dθ/dζ$) is an integrable field-line Lagrangian and $ε$ is a perturbation parameter. This is explained using a perturbative construction of the auxiliary poloidal angle $Θ$ that corrects QFMin surfaces so they are also ghost surfaces. The difference between the corrected and uncorrected surfaces is $O(ε^2)$, explaining the observed smallness of this difference. An alternative definition of ghost surfaces is also introduced, based on an action-gradient flow in $Θ$, which appears to have superior properties when unified with QFMin surfaces.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Are ghost surfaces quadratic-flux-minimizing?

Two candidates for "almost-invariant" toroidal surfaces passing through magnetic islands, namely quadratic-flux-minimizing (QFMin) surfaces and ghost surfaces, use families of periodic pseudo-orbits (i.e. paths for which the action is not exactly extremal). QFMin pseudo-orbits, which are coordinate-dependent, are field lines obtained from a modified magnetic field, and ghost-surface pseudo-orbits are obtained by displacing closed field lines in the direction of steepest descent of magnetic action, $\oint \vec{A}\cdot\mathbf{dl}$. A generalized Hamiltonian definition of ghost surfaces is given and specialized to the usual Lagrangian definition. A modified Hamilton's Principle is introduced that allows the use of Lagrangian integration for calculation of the QFMin pseudo-orbits. Numerical calculations show QFMin and Lagrangian ghost surfaces give very similar results for a chaotic magnetic field perturbed from an integrable case, and this is explained using a perturbative construction of an auxiliary poloidal angle for which QFMin and Lagrangian ghost surfaces are the same up to second order. While presented in the context of 3-dimensional magnetic field line systems, the concepts are applicable to defining almost-invariant tori in other $1{1/2}$ degree-of-freedom nonintegrable Lagrangian/Hamiltonian systems.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Hamilton--Jacobi theory for continuation of magnetic field across a toroidal surface supporting a plasma pressure discontinuity

The vanishing of the divergence of the total stress tensor (magnetic plus kinetic) in a neighborhood of an equilibrium plasma containing a toroidal surface of discontinuity gives boundary and jump conditions that strongly constrain allowable continuations of the magnetic field across the surface. The boundary conditions allow the magnetic fields on either side of the discontinuity surface to be described by surface magnetic potentials, reducing the continuation problem to that of solving a Hamilton--Jacobi equation. The characteristics of this equation obey Hamiltonian equations of motion, and a necessary condition for the existence of a continued field across a general toroidal surface is that there exist invariant tori in the phase space of this Hamiltonian system. It is argued from the Birkhoff theorem that existence of such an invariant torus is also, in general, sufficient for continuation to be possible. An important corollary is that the rotational transform of the continued field on a surface of discontinuity must, generically, be irrational.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Nonlinear Simulation of Drift Wave Turbulence

In a two-dimensional version of the modified Hasegawa-Wakatani (HW) model, which describes electrostatic resistive drift wave turbulence, the resistive coupling between vorticity and density does not act on the zonal components ($k_{y}=0$). It is therefore necessary to modify the HW model to treat the zonal components properly. The modified equations are solved numerically, and visualization and analysis of the solutions show generation of stable zonal flows, through conversion of turbulent kinetic energy, and the consequent turbulence and transport suppression. It is demonstrated by comparison that the modification is essential for generation of zonal flows.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Zonal flow generation by modulational instability

This paper gives a pedagogic review of the envelope formalism for excitation of zonal flows by nonlinear interactions of plasma drift waves or Rossby waves, described equivalently by the Hasegawa-Mima (HM) equation or the quasigeostrophic barotropic potential vorticity equation, respectively. In the plasma case a modified form of the HM equation, which takes into account suppression of the magnetic-surface-averaged electron density response by a small amount of rotational transform, is also analyzed. Excitation of zonal mean flow by a modulated wave train is particularly strong in the modified HM case. A local dispersion relation for a coherent wave train is calculated by linearizing about a background mean flow and used to find the nonlinear frequency shift by inserting the nonlinearly excited mean flow. Using the generic nonlinear Schroedinger equation about a uniform carrier wave, the criterion for instability of small modulations of the wave train is found, as is the maximum growth rate and phase velocity of the modulations and zonal flows, in both the modified and unmodified cases.

0 citations0 reviewsNo code linkedOpen access