Researcher profile

S. R. Hudson

S. R. Hudson contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
9works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

9 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.

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.

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.