Source author record

S. R. Hudson

S. R. Hudson 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

12works
6topics
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

12 published item(s)

preprint2020arXiv

Bifurcations of the magnetic axis and the alternating-hyperbolic sawtooth

We present a sawtooth model that explains observations where the central safety factor, $q_0$, stays well below one, which is irreconcilable with current models that predict a reset to $q_0=1$ after the crash. We identify the structure of the field around the magnetic axis with elements of the Lie group $\mathrm{SL}(2,\mathbb{R})$ and find a transition to an alternating-hyperbolic geometry when $q_0=2/3$. This transition is driven by an ideal MHD instability and leads to a chaotic magnetic field near the axis.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.