-
An impulsive geomagnetic effect from an early-impulsive flare
Authors:
Hugh S. Hudson,
Edward. W. Cliver,
Lyndsay Fletcher,
Declan A. Diver,
Peter T. Gallagher,
Ying Li,
Christopher M. J. Osborne,
Craig Stark,
Yang Su
Abstract:
The geomagnetic "solar flare effect" (SFE) results from excess ionization in the Earth's ionosphere, famously first detected at the time of the Carrington flare in 1859. This indirect detection of a flare constituted one of the first cases of "multimessenger astronomy," whereby solar ionizing radiation stimulates ionospheric currents. Well-observed SFEs have few-minute time scales and perturbation…
▽ More
The geomagnetic "solar flare effect" (SFE) results from excess ionization in the Earth's ionosphere, famously first detected at the time of the Carrington flare in 1859. This indirect detection of a flare constituted one of the first cases of "multimessenger astronomy," whereby solar ionizing radiation stimulates ionospheric currents. Well-observed SFEs have few-minute time scales and perturbations of >10 nT, with the greatest events reaching above 100 nT. In previously reported cases the SFE time profiles tend to resemble those of solar soft X-ray emission, which ionizes the D-region; there is also a less-well-studied contribution from Lyman-alpha. We report here a specific case, from flare SOL2024-03-10 (M7.4), in which an impulsive SFE deviated from this pattern. This flare contained an "early impulsive" component of exceptionally hard radiation, extending up to gamma-ray energies above 1 MeV, distinctly before the bulk of the flare soft X-ray emission. We can characterize the spectral distribution of this early-impulsive component in detail, thanks to the modern extensive wavelength coverage. A more typical gradual SFE occurred during the flare's main phase. We suggest that events of this type warrant exploration of the solar physics in the "impulse response" limit of very short time scales.
△ Less
Submitted 12 July, 2024;
originally announced July 2024.
-
Hard X-rays from the deep solar atmosphere: An unusual UV burst with flare properties
Authors:
L. P. Chitta,
I. G. Hannah,
L. Fletcher,
H. S. Hudson,
P. R. Young,
S. Krucker,
H. Peter
Abstract:
Explosive transient events occur throughout the solar atmosphere. The differing manifestations range from coronal mass ejections to Ellermann bombs. The former may have negligible signatures in the lower atmosphere, and the latter may have negligible nonthermal emissions such as hard X-radiation. A solar flare generally involves a broad range of emission signatures. Using a suite of four space-bor…
▽ More
Explosive transient events occur throughout the solar atmosphere. The differing manifestations range from coronal mass ejections to Ellermann bombs. The former may have negligible signatures in the lower atmosphere, and the latter may have negligible nonthermal emissions such as hard X-radiation. A solar flare generally involves a broad range of emission signatures. Using a suite of four space-borne telescopes, we report a solar event that combines aspects of simple UV bursts and hard X-ray emitting flares at the same time. The event is a compact C-class flare in active region AR11861, SOL2013-10-12T00:30. By fitting a combined isothermal and nonthermal model to the hard X-ray spectrum, we inferred plasma temperatures in excess of 15\,MK and a nonthermal power of about $3\times10^{27}$\,erg\,s$^{-1}$ in this event. Despite these high temperatures and evidence for nonthermal particles, the flare was mostly confined to the chromosphere. However, the event lacked clear signatures of UV spectral lines, such as the Fe\,{\sc xii} 1349\,Å and Fe\,{\sc xxi} 1354\,Å emission lines, which are characteristic of emission from hotter plasma with a temperature over 1\,MK. Moreover, the event exhibited very limited signatures in the extreme-UV wavelengths. Our study indicates that a UV burst -- hard X-ray flare hybrid phenomenon exists in the low solar atmosphere. Plasma that heats to high temperatures coupled with particle acceleration by magnetic energy that is released directly in the lower atmosphere sheds light on the nature of active region core heating and on inferences of flare signatures.
△ Less
Submitted 8 July, 2024;
originally announced July 2024.
-
Flexible Stellarator Physics Facility
Authors:
F. I. Parra,
S. -G. Baek,
M. Churchill,
D. R. Demers,
B. Dudson,
N. M. Ferraro,
B. Geiger,
S. Gerhardt,
K. C. Hammond,
S. Hudson,
R. Jorge,
E. Kolemen,
D. M. Kriete,
S. T. A. Kumar,
M. Landreman,
C. Lowe,
D. A. Maurer,
F. Nespoli,
N. Pablant,
M. J. Pueschel,
A. Punjabi,
J. A. Schwartz,
C. P. S. Swanson,
A. M. Wright
Abstract:
We propose to build a Flexible Stellarator Physics Facility to explore promising regions of the vast parameter space of disruption-free stellarator solutions for Fusion Pilot Plants (FPPs).
We propose to build a Flexible Stellarator Physics Facility to explore promising regions of the vast parameter space of disruption-free stellarator solutions for Fusion Pilot Plants (FPPs).
△ Less
Submitted 4 July, 2024;
originally announced July 2024.
-
Constructing nested coordinates inside strongly shaped toroids using an action principle
Authors:
Zeno Tecchiolli,
Stuart Hudson,
Joaquim Loizu,
Robert Köberl,
Florian Hindenlang,
Brenno De Lucca
Abstract:
A new approach for constructing polar-like boundary-conforming coordinates inside a toroid with strongly shaped cross-sections is presented. A coordinate mapping is obtained through a variational approach, which involves identifying extremal points of a proposed action in the mapping space from [0, 2π] x [0, 2π] x [0, 1] to a toroidal domain in R3. This approach employs an action built on the squa…
▽ More
A new approach for constructing polar-like boundary-conforming coordinates inside a toroid with strongly shaped cross-sections is presented. A coordinate mapping is obtained through a variational approach, which involves identifying extremal points of a proposed action in the mapping space from [0, 2π] x [0, 2π] x [0, 1] to a toroidal domain in R3. This approach employs an action built on the squared Jacobian and radial length. Extensive testing is conducted on general toroidal boundaries using a global Fourier-Zernike basis via action minimization. The results demonstrate successful coordinate construction capable of accurately describing strongly shaped toroidal domains. The coordinate construction is successfully applied to the computation of 3D MHD equilibria in the GVEC code where the use of traditional coordinate construction by interpolation from the boundary failed.
△ Less
Submitted 13 May, 2024;
originally announced May 2024.
-
Application of Lagrangian techniques for calculating the on-axis rotational transform
Authors:
S. Guinchard,
W. Sengupta,
S. R. Hudson
Abstract:
The Floquet exponents of periodic field lines are studied through the variations of the magnetic action on the magnetic axis, which is assumed to be elliptical. The near-axis formalism developed by Mercier, Solov'ev and Shafranov is combined with a Lagrangian approach. The on-axis Floquet exponent is shown to coincide with the on-axis rotational transform, and this is a coordinate-independent resu…
▽ More
The Floquet exponents of periodic field lines are studied through the variations of the magnetic action on the magnetic axis, which is assumed to be elliptical. The near-axis formalism developed by Mercier, Solov'ev and Shafranov is combined with a Lagrangian approach. The on-axis Floquet exponent is shown to coincide with the on-axis rotational transform, and this is a coordinate-independent result. A discrete solution suitable for numerical implementation is introduced, which gives the Floquet exponents as solutions to an eigenvalue problem. This discrete formalism expresses the exponents as the eigenvalues of a $6\times 6$ matrix.
△ Less
Submitted 26 April, 2024;
originally announced April 2024.
-
Exploration of the parameter space of quasisymmetric stellarator vacuum fields through adjoint optimisation
Authors:
Richard Nies,
Elizabeth J. Paul,
Dario Panici,
Stuart R. Hudson,
Amitava Bhattacharjee
Abstract:
Optimising stellarators for quasisymmetry leads to strongly reduced collisional transport and energetic particle losses compared to unoptimised configurations. Though stellarators with precise quasisymmetry have been obtained in the past, it remains unclear how broad the parameter space is where good quasisymmetry may be achieved. We study the range of aspect ratio and rotational transform values…
▽ More
Optimising stellarators for quasisymmetry leads to strongly reduced collisional transport and energetic particle losses compared to unoptimised configurations. Though stellarators with precise quasisymmetry have been obtained in the past, it remains unclear how broad the parameter space is where good quasisymmetry may be achieved. We study the range of aspect ratio and rotational transform values for which stellarators with excellent quasisymmetry on the boundary can be obtained. A large number of Fourier harmonics is included in the boundary representation, which is made computationally tractable by the use of adjoint methods to enable fast gradient-based optimisation, and by the direct optimisation of vacuum magnetic fields, which converge more robustly compared to solutions from magnetohydrostatics. Several novel configurations are presented, including stellarators with record levels of quasisymmetry on a surface; three field period quasiaxisymmetric stellarators with substantial magnetic shear, and compact quasisymmetric stellarators at low aspect ratios similar to tokamaks.
△ Less
Submitted 23 May, 2024; v1 submitted 2 April, 2024;
originally announced April 2024.
-
DAS-N2N: Machine learning Distributed Acoustic Sensing (DAS) signal denoising without clean data
Authors:
Sacha Lapins,
Antony Butcher,
J. -Michael Kendall,
Thomas S. Hudson,
Anna L. Stork,
Maximilian J. Werner,
Jemma Gunning,
Alex M. Brisbourne
Abstract:
This article presents a weakly supervised machine learning method, which we call DAS-N2N, for suppressing strong random noise in distributed acoustic sensing (DAS) recordings. DAS-N2N requires no manually produced labels (i.e., pre-determined examples of clean event signals or sections of noise) for training and aims to map random noise processes to a chosen summary statistic, such as the distribu…
▽ More
This article presents a weakly supervised machine learning method, which we call DAS-N2N, for suppressing strong random noise in distributed acoustic sensing (DAS) recordings. DAS-N2N requires no manually produced labels (i.e., pre-determined examples of clean event signals or sections of noise) for training and aims to map random noise processes to a chosen summary statistic, such as the distribution mean, median or mode, whilst retaining the true underlying signal. This is achieved by splicing (joining together) two fibres hosted within a single optical cable, recording two noisy copies of the same underlying signal corrupted by different independent realizations of random observational noise. A deep learning model can then be trained using only these two noisy copies of the data to produce a near fully-denoised copy. Once the model is trained, only noisy data from a single fibre is required. Using a dataset from a DAS array deployed on the surface of the Rutford Ice Stream in Antarctica, we demonstrate that DAS-N2N greatly suppresses incoherent noise and enhances the signal-to-noise ratios (SNR) of natural microseismic icequake events. We further show that this approach is inherently more efficient and effective than standard stop/pass band and white noise (e.g., Wiener) filtering routines, as well as a comparable self-supervised learning method based on masking individual DAS channels. Our preferred model for this task is lightweight, processing 30 seconds of data recorded at a sampling frequency of 1000 Hz over 985 channels (approx. 1 km of fiber) in $<$1 s. Due to the high noise levels in DAS recordings, efficient data-driven denoising methods, such as DAS-N2N, will prove essential to time-critical DAS earthquake detection, particularly in the case of microseismic monitoring.
△ Less
Submitted 24 November, 2023; v1 submitted 17 April, 2023;
originally announced April 2023.
-
Structure of pressure-gradient-driven current singularity in ideal magnetohydrodynamic equilibrium
Authors:
Yi-Min Huang,
Yao Zhou,
Joaquim Loizu,
Stuart Hudson,
Amitava Bhattacharjee
Abstract:
Singular currents typically appear on rational surfaces in non-axisymmetric ideal magnetohydrodynamic equilibria with a continuum of nested flux surfaces and a continuous rotational transform. These currents have two components: a surface current (Dirac $δ$-function in flux surface labeling) that prevents the formation of magnetic islands and an algebraically divergent Pfirsch--Schlüter current de…
▽ More
Singular currents typically appear on rational surfaces in non-axisymmetric ideal magnetohydrodynamic equilibria with a continuum of nested flux surfaces and a continuous rotational transform. These currents have two components: a surface current (Dirac $δ$-function in flux surface labeling) that prevents the formation of magnetic islands and an algebraically divergent Pfirsch--Schlüter current density when a pressure gradient is present across the rational surface. At flux surfaces adjacent to the rational surface, the traditional treatment gives the Pfirsch--Schlüter current density scaling as $J\sim1/Δι$, where $Δι$ is the difference of the rotational transform relative to the rational surface. If the distance $s$ between flux surfaces is proportional to $Δι$, the scaling relation $J\sim1/Δι\sim1/s$ will lead to a paradox that the Pfirsch--Schlüter current is not integrable. In this work, we investigate this issue by considering the pressure-gradient-driven singular current in the Hahm\textendash Kulsrud\textendash Taylor problem, which is a prototype for singular currents arising from resonant magnetic perturbations. We show that not only the Pfirsch--Schlüter current density but also the diamagnetic current density are divergent as $\sim1/Δι$. However, due to the formation of a Dirac $δ$-function current sheet at the rational surface, the neighboring flux surfaces are strongly packed with $s\sim(Δι)^{2}$. Consequently, the singular current density $J\sim1/\sqrt{s}$, making the total current finite, thus resolving the paradox.
△ Less
Submitted 9 February, 2024; v1 submitted 3 March, 2023;
originally announced March 2023.
-
Periodic Korteweg-de Vries soliton potentials generate magnetic field strength with excellent quasisymmetry
Authors:
W. Sengupta,
N. Nikulsin,
S. Buller,
R. Madan,
E. J. Paul,
R. Nies,
A. A. Kaptanoglu,
S. R. Hudson,
A. Bhattacharjee
Abstract:
Quasisymmetry (QS) is a hidden symmetry of the magnetic field strength, $B$, that confines charged particles effectively in a three-dimensional toroidal plasma equilibrium. Here, we show that QS has a deep connection to the underlying symmetry that makes solitons possible. Our approach uncovers a hidden lower dimensionality of $B$ on a magnetic flux surface, which could make stellarator optimizati…
▽ More
Quasisymmetry (QS) is a hidden symmetry of the magnetic field strength, $B$, that confines charged particles effectively in a three-dimensional toroidal plasma equilibrium. Here, we show that QS has a deep connection to the underlying symmetry that makes solitons possible. Our approach uncovers a hidden lower dimensionality of $B$ on a magnetic flux surface, which could make stellarator optimization schemes significantly more efficient. Recent numerical breakthroughs have yielded configurations with excellent volumetric QS and surprisingly low magnetic shear. Our approach elucidates why the magnetic shear is low in these configurations. Furthermore, we deduce an upper bound on the maximum toroidal volume that can be quasisymmetric and verify it for the Landreman-Paul precise quasiaxisymmetric (QA) stellarator configuration. In the neighborhood of the outermost surface, we show that the B approaches the form of the 1-soliton reflectionless potential. We present three independent approaches to demonstrate that quasisymmetric $B$ is described by well-known integrable systems such as the Korteweg-de Vries (KdV) equation. The first approach is weakly nonlinear multiscale perturbation theory, which highlights the crucial role that magnetic shear plays in QS. We show that the overdetermined problem of finding quasisymmetric vacuum fields admits solutions for which the rotational transform is not free but highly constrained. We obtain the KdV equation (and, more specifically, Gardner's equation for certain choices of parameters). Our second approach is non-perturbative and based on ensuring single-valuedness of $B$, which directly leads to its Painlevé property shared by the KdV equation. Our third approach uses machine learning, trained on a large dataset of numerically optimized quasisymmetric stellarators. We robustly recover the KdV (and Gardner's) equation from the data.
△ Less
Submitted 24 May, 2024; v1 submitted 27 February, 2023;
originally announced February 2023.
-
Prospects of Detecting Non-thermal Protons in Solar Flares via Lyman Line Spectroscopy: Revisiting the Orrall-Zirker Effect
Authors:
Graham S. Kerr,
Joel C. Allred,
Adam F. Kowalski,
Ryan O. Milligan,
Hugh S. Hudson,
Natalia Zambrana Prado,
Therese A. Kucera,
Jeffrey W. Brosius
Abstract:
Solar flares are efficient particle accelerators, with a substantial fraction of the energy released manifesting as non-thermal particles. While the role that non-thermal electrons play in transporting flare energy is well studied, the properties and importance of non-thermal protons is rather less well understood. This is in large part due to the paucity of diagnostics, particularly at the lower-…
▽ More
Solar flares are efficient particle accelerators, with a substantial fraction of the energy released manifesting as non-thermal particles. While the role that non-thermal electrons play in transporting flare energy is well studied, the properties and importance of non-thermal protons is rather less well understood. This is in large part due to the paucity of diagnostics, particularly at the lower-energy (deka-keV) range of non-thermal proton distributions in flares. One means to identify the presence of deka-keV protons is by an effect originally described by \cite{1976ApJ...208..618O}. In the Orrall-Zirker effect, non-thermal protons interact with ambient neutral hydrogen, and via charge exchange produce a population of energetic neutral atoms (ENAs) in the chromosphere. These ENAs subsequently produce an extremely redshifted photon in the red wings of hydrogen spectral lines. We revisit predictions of the strength of this effect using modern interaction cross-sections, and numerical models capable of self-consistently simulating the flaring non-equilibrium ionization stratification, and the non-thermal proton distribution (and, crucially, their feedback on each other). We synthesize both the thermal and non-thermal emission from \lya\ and \lyb, the most promising lines that may exhibit a detectable signal. These new predictions are are weaker and more transient than prior estimates, but the effects should be detectable in fortuitous circumstances. We degrade the \lyb\ emission to the resolution of the Spectral Imaging of the Coronal Environment (SPICE) instrument on board Solar Orbiter, demonstrating that though likely difficult, it should be possible to detect the presence of non-thermal protons in flares observed by SPICE.
△ Less
Submitted 8 February, 2023; v1 submitted 3 February, 2023;
originally announced February 2023.
-
Bayesian optimization of laser-plasma accelerators assisted by reduced physical models
Authors:
A. Ferran Pousa,
S. Jalas,
M. Kirchen,
A. Martinez de la Ossa,
M. Thévenet,
S. Hudson,
J. Larson,
A. Huebl,
J. -L. Vay,
R. Lehe
Abstract:
Particle-in-cell simulations are among the most essential tools for the modeling and optimization of laser-plasma accelerators, since they reproduce the physics from first principles. However, the high computational cost associated with them can severely limit the scope of parameter and design optimization studies. Here, we show that a multitask Bayesian optimization algorithm can be used to mitig…
▽ More
Particle-in-cell simulations are among the most essential tools for the modeling and optimization of laser-plasma accelerators, since they reproduce the physics from first principles. However, the high computational cost associated with them can severely limit the scope of parameter and design optimization studies. Here, we show that a multitask Bayesian optimization algorithm can be used to mitigate the need for such high-fidelity simulations by incorporating information from inexpensive evaluations of reduced physical models. In a proof-of-principle study, where a high-fidelity optimization with FBPIC is assisted by reduced-model simulations with Wake-T, the algorithm demonstrates an order-of-magnitude speedup. This opens a path for the cost-effective optimization of laser-plasma accelerators in large parameter spaces, an important step towards fulfilling the high beam quality requirements of future applications.
△ Less
Submitted 23 December, 2022;
originally announced December 2022.
-
On the relationship between the multi-region relaxed variational principle and resistive inner layer theory
Authors:
A. Kumar,
J. Loizu,
M. J. Hole,
Z. Qu,
S. R. Hudson,
R. L Dewar
Abstract:
We show that the variational energy principle of multi-region relaxed magnetohydrodynamic (MRxMHD) model can be used to predict finite-pressure linear tearing instabilities. In this model, the plasma volume is sliced into sub-volumes separated by "ideal interfaces", and in each volume the magnetic field relaxes to a Taylor state where the pressure gradient $\nabla p = 0$. The MRxMHD model is imple…
▽ More
We show that the variational energy principle of multi-region relaxed magnetohydrodynamic (MRxMHD) model can be used to predict finite-pressure linear tearing instabilities. In this model, the plasma volume is sliced into sub-volumes separated by "ideal interfaces", and in each volume the magnetic field relaxes to a Taylor state where the pressure gradient $\nabla p = 0$. The MRxMHD model is implemented in the SPEC code so that the equilibrium solution in each region is computed while the preserving force balance across the interfaces. As SPEC computes the Hessian matrix (a discretized stability matrix), the stability of an MRxMHD equilibrium can also be computed with SPEC. In this article, using SPEC, we investigate the effect of local pressure gradients and the $\nabla p = 0$ in the vicinity of the resonant surface of a tearing mode. For low beta plasma, we have been able to illustrate a relationship between the resistive singular layer theory [Coppi et al. (1966) Nucl. Fusion 6 101, Glasser et al. The Physics of Fluids 18, 875-888 (1975)], and the MRxMHD model. Within the singular layer, the volume-averaged magnetic helicity and the flux-averaged toroidal flux are shown to be the invariants for the linear tearing modes in SPEC simulations. Our technique to compute MRxMHD stability is first tested numerically in cylindrical tokamak and its application in toroidal geometry is demonstrated. We demonstrate an agreement between the stability boundary obtained with SPEC simulation and the resistive inner layer theories.
△ Less
Submitted 11 September, 2022;
originally announced September 2022.
-
Comparison of Multiobjective Optimization Methods for the LCLS-II Photoinjector
Authors:
Nicole Neveu,
Tyler H. Chang,
Paris Franz,
Stephen Hudson,
Jeffrey Larson
Abstract:
Particle accelerators are among some of the largest science experiments in the world and can consist of thousands of components with a wide variety of input ranges. These systems can easily become unwieldy optimization problems during design and operations studies. Starting in the early 2000s, searching for better beam dynamics configurations became synonymous with heuristic optimization methods i…
▽ More
Particle accelerators are among some of the largest science experiments in the world and can consist of thousands of components with a wide variety of input ranges. These systems can easily become unwieldy optimization problems during design and operations studies. Starting in the early 2000s, searching for better beam dynamics configurations became synonymous with heuristic optimization methods in the accelerator physics community. Genetic algorithms and particle swarm optimization are currently the most widely used. These algorithms can take thousands of simulation evaluations to find optimal solutions for one machine prototype. For large facilities such as the Linac Coherent Light Source (LCLS) and others, this equates to a limited exploration of many possible design configurations. In this paper, the LCLS-II photoinjector is optimized with three optimization algorithms. All optimizations were started from both a uniform random and Latin hypercube sample. In all cases, the optimizations started from Latin hypercube samples outperformed optimizations started from uniform samples. All three algorithms were able to optimize the photoinjector, with the model-based methods approximating the Pareto front in fewer simulation evaluations. This work, in combination with previous optimization observations, indicates objective penalties have a strong impact on the efficiency of such methods. In general, we recommend heuristic methods for initial optimizations and model-based methods when information about the objective space is available.
△ Less
Submitted 31 October, 2022; v1 submitted 8 June, 2022;
originally announced June 2022.
-
Adjoint methods for quasisymmetry of vacuum fields on a surface
Authors:
Richard Nies,
Elizabeth J. Paul,
Stuart R. Hudson,
Amitava Bhattacharjee
Abstract:
Adjoint methods can speed up stellarator optimisation by providing gradient information more efficiently compared to finite-difference evaluations. Adjoint methods are herein applied to vacuum magnetic fields, with objective functions targeting quasisymmetry and a rotational transform value on a surface. To measure quasisymmetry, a novel way of evaluating approximate flux coordinates on a single f…
▽ More
Adjoint methods can speed up stellarator optimisation by providing gradient information more efficiently compared to finite-difference evaluations. Adjoint methods are herein applied to vacuum magnetic fields, with objective functions targeting quasisymmetry and a rotational transform value on a surface. To measure quasisymmetry, a novel way of evaluating approximate flux coordinates on a single flux surface without the assumption of a neighbourhood of flux surfaces is proposed. The shape gradients obtained from the adjoint formalism are evaluated numerically and verified against finite-difference evaluations.
△ Less
Submitted 25 August, 2021;
originally announced August 2021.
-
Numerical study of $δ$-function current sheets arising from resonant magnetic perturbations
Authors:
Yi-Min Huang,
Stuart R. Hudson,
Joaquim Loizu,
Yao Zhou,
Amitava Bhattacharjee
Abstract:
General three-dimensional toroidal ideal magnetohydrodynamic equilibria with a continuum of nested flux surfaces are susceptible to forming singular current sheets when resonant perturbations are applied. The presence of singular current sheets indicates that, in the presence of non-zero resistivity, magnetic reconnection will ensue, leading to the formation of magnetic islands and potentially reg…
▽ More
General three-dimensional toroidal ideal magnetohydrodynamic equilibria with a continuum of nested flux surfaces are susceptible to forming singular current sheets when resonant perturbations are applied. The presence of singular current sheets indicates that, in the presence of non-zero resistivity, magnetic reconnection will ensue, leading to the formation of magnetic islands and potentially regions of stochastic field lines when islands overlap. Numerically resolving singular current sheets in the ideal MHD limit has been a significant challenge. This work presents numerical solutions of the Hahm-Kulsrud-Taylor (HKT) problem, which is a prototype for resonant singular current sheet formation. The HKT problem is solved by two codes: a Grad-Shafranov (GS) solver and the SPEC code. The GS solver has built-in nested flux surfaces with prescribed magnetic fluxes. The SPEC code implements multi-region relaxed magnetohydrodynamics (MRxMHD), where the solution relaxes to a Taylor state in each region while maintaining force balance across the interfaces between regions. As the number of regions increases, the MRxMHD solution approaches the ideal MHD solution assuming a continuum of nested flux surfaces. We demonstrate excellent agreement between the numerical solutions obtained from the two codes through a thorough convergence study.
△ Less
Submitted 1 September, 2023; v1 submitted 20 August, 2021;
originally announced August 2021.
-
Heat conduction in an irregular magnetic field: Part II. Heat transport as a measure of the effective non-integrable volume
Authors:
Elizabeth J. Paul,
Stuart R. Hudson,
Per Helander
Abstract:
Given the large anisotropy of transport processes in magnetized plasmas, the magnetic field structure can strongly impact heat diffusion: magnetic surfaces and cantori form barriers to transport while chaotic layers and island structures can degrade confinement. When a small but non-zero amount of perpendicular diffusion is included, the structure of the magnetic field becomes less important, allo…
▽ More
Given the large anisotropy of transport processes in magnetized plasmas, the magnetic field structure can strongly impact heat diffusion: magnetic surfaces and cantori form barriers to transport while chaotic layers and island structures can degrade confinement. When a small but non-zero amount of perpendicular diffusion is included, the structure of the magnetic field becomes less important, allowing pressure gradients to be supported across chaotic regions and island chains. We introduce a metric for the effective volume over which the local parallel diffusion dominates based on the solution to the anisotropic heat diffusion equation. To validate this metric, we consider model fields with a single island chain and a strongly chaotic layer for which analytic predictions of the relative parallel and perpendicular transport can be made. We also analyze critically chaotic fields produced from different sets of perturbations, highlighting the impact of the mode number spectrum on the heat transport. Our results indicate that this metric coincides with the effective volume of non-integrability in the limit $κ_{\perp} \rightarrow 0$. We propose that this metric be used to assess the impact of non-integrability on the heat transport in stellarator equilibria.
△ Less
Submitted 6 February, 2022; v1 submitted 13 August, 2021;
originally announced August 2021.
-
On the non-existence of stepped-pressure equilibria far from symmetry
Authors:
Z. S. Qu,
S. R. Hudson,
R. L. Dewar,
J. Loizu,
M. J. Hole
Abstract:
The Stepped Pressure Equilibrium Code (SPEC) [Hudson et al., Phys. Plasmas 19, 112502 (2012)] has been successful in the construction of equilibria in 3D configurations that contain a mixture of flux surfaces, islands and chaotic magnetic field lines. In this model, the plasma is sliced into sub-volumes separated by ideal interfaces, and in each volume the magnetic field is a Beltrami field. In th…
▽ More
The Stepped Pressure Equilibrium Code (SPEC) [Hudson et al., Phys. Plasmas 19, 112502 (2012)] has been successful in the construction of equilibria in 3D configurations that contain a mixture of flux surfaces, islands and chaotic magnetic field lines. In this model, the plasma is sliced into sub-volumes separated by ideal interfaces, and in each volume the magnetic field is a Beltrami field. In the cases where the system is far from possessing a continuous symmetry, such as in stellarators, the existence of solutions to a stepped-pressure equilibrium with given constraints, such as a multi-region relaxed MHD minimum energy state, is not guaranteed but is often taken for granted. Using SPEC, we have studied two different scenarios in which a solution fails to exist in a slab with analytic boundary perturbations. We found that with a large boundary perturbation, a certain interface becomes fractal, corresponding to the break up of a Kolmogorov-Arnold-Moser (KAM) surface. Moreover, an interface can only support a maximum pressure jump while a solution of the magnetic field consistent with the force balance condition can be found. An interface closer to break-up can support a smaller pressure jump. We discovered that the pressure jump can push the interface closer to being non-smooth through force balance, thus significantly decreasing the maximum pressure it can support. Our work shows that a convergence study must be performed on a SPEC equilibrium with interfaces close to break-up. These results may also provide insights into the choice of interfaces and have applications in finding out the maximum pressure a machine can support.
△ Less
Submitted 14 September, 2021; v1 submitted 4 August, 2021;
originally announced August 2021.
-
Computation of the Biot-Savart line integral with higher-order convergence using straight segments
Authors:
Nick McGreivy,
Caoxiang Zhu,
Lee Gunderson,
Stuart R. Hudson
Abstract:
One common approach to computing the magnetic field produced by a filamentary current-carrying coil is to approximate the coil as a series of straight segments. The Biot-Savart field from each straight segment is known analytically. However, if the endpoints of the straight segments are chosen to lie on the coil, then the accuracy of the Biot-Savart computation is generally only second-order in th…
▽ More
One common approach to computing the magnetic field produced by a filamentary current-carrying coil is to approximate the coil as a series of straight segments. The Biot-Savart field from each straight segment is known analytically. However, if the endpoints of the straight segments are chosen to lie on the coil, then the accuracy of the Biot-Savart computation is generally only second-order in the number of endpoints. We propose a simple modification: shift each endpoint off the coil in the outwards normal direction by an amount proportional to the local curvature. With this modification, the Biot-Savart accuracy increases to fourth order and the numerical error is dramatically reduced for a given number of discretization points.
△ Less
Submitted 27 July, 2021; v1 submitted 26 May, 2021;
originally announced May 2021.
-
Combined plasma-coil optimization algorithms
Authors:
Sophia A. Henneberg,
Stuart R. Hudson,
David Pfefferlé,
Per Helander
Abstract:
Combined plasma-coil optimization approaches for designing stellarators are discussed and a new method for calculating free-boundary equilibria is proposed. Four distinct categories of stellarator optimization, two of which are novel approaches, are the fixed-boundary optimization, the generalized fixed-boundary optimization, the quasi free-boundary optimization, and the free-boundary (coil) optim…
▽ More
Combined plasma-coil optimization approaches for designing stellarators are discussed and a new method for calculating free-boundary equilibria is proposed. Four distinct categories of stellarator optimization, two of which are novel approaches, are the fixed-boundary optimization, the generalized fixed-boundary optimization, the quasi free-boundary optimization, and the free-boundary (coil) optimization. These are described using the multi-region relaxed magnetohydrodynmics (MRxMHD) energy functional, the Biot-Savart integral, the coil-penalty functional and the virtual casing integral, and their derivatives. The proposed free-boundary equilibrium calculation differs from existing methods in how the boundary-value problem is posed, and for the new approach it seems that there is not an associated energy minimization principle because a non-symmetric functional arises. We propose to solve the weak formulation of this problem using a spectral-Galerkin method, and this will reduce the free-boundary equilibrium calculation to something comparable to a fixed-boundary calculation. In our discussion of combined plasma-coil optimization algorithms, we emphasize the importance of the stability matrix.
△ Less
Submitted 16 December, 2020;
originally announced December 2020.
-
Solar Flare Energy Partitioning and Transport -- the Impulsive Phase (a Heliophysics 2050 White Paper)
Authors:
Graham S. Kerr,
Meriem Alaoui,
Joel C. Allred,
Nicholas H. Bian,
Brian R. Dennis,
A. Gordon Emslie,
Lyndsay Fletcher,
Silvina Guidoni,
Laura A. Hayes,
Gordon D. Holman,
Hugh S. Hudson,
Judith T. Karpen,
Adam F. Kowalski,
Ryan O. Milligan,
Vanessa Polito,
Jiong Qiu,
Daniel F. Ryan
Abstract:
Solar flares are a fundamental component of solar eruptive events (SEEs; along with solar energetic particles, SEPs, and coronal mass ejections, CMEs). Flares are the first component of the SEE to impact our atmosphere, which can set the stage for the arrival of the associated SEPs and CME. Magnetic reconnection drives SEEs by restructuring the solar coronal magnetic field, liberating a tremendous…
▽ More
Solar flares are a fundamental component of solar eruptive events (SEEs; along with solar energetic particles, SEPs, and coronal mass ejections, CMEs). Flares are the first component of the SEE to impact our atmosphere, which can set the stage for the arrival of the associated SEPs and CME. Magnetic reconnection drives SEEs by restructuring the solar coronal magnetic field, liberating a tremendous amount of energy which is partitioned into various physical manifestations: particle acceleration, mass and magnetic-field eruption, atmospheric heating, and the subsequent emission of radiation as solar flares. To explain and ultimately predict these geoeffective events, the heliophysics community requires a comprehensive understanding of the processes that transform and distribute stored magnetic energy into other forms, including the broadband radiative enhancement that characterises flares. This white paper, submitted to the Heliophysics 2050 Workshop, discusses the flare impulsive phase part of SEEs, setting out the questions that need addressing via a combination of theoretical, modelling, and observational research. In short, by 2050 we must determine the mechanisms of particle acceleration and propagation, and must push beyond the paradigm of energy transport via nonthermal electron beams, to also account for accelerated protons & ions and downward directed Alfven waves.
△ Less
Submitted 17 September, 2020;
originally announced September 2020.
-
Optimized finite-build stellarator coils using automatic differentiation
Authors:
Nick McGreivy,
Stuart R. Hudson,
Caoxiang Zhu
Abstract:
A new stellarator coil design code is introduced that optimizes the position and winding pack orientation of finite-build coils. The new code, called FOCUSADD, performs gradient-based optimization in a high-dimensional, non-convex space. The derivatives with respect to parameters of finite-build coils are easily and efficiently computed using automatic differentiation. FOCUSADD parametrizes coil p…
▽ More
A new stellarator coil design code is introduced that optimizes the position and winding pack orientation of finite-build coils. The new code, called FOCUSADD, performs gradient-based optimization in a high-dimensional, non-convex space. The derivatives with respect to parameters of finite-build coils are easily and efficiently computed using automatic differentiation. FOCUSADD parametrizes coil positions in free space using a Fourier series and uses a multi-filament approximation to the coil winding pack. The orientation of the winding pack is parametrized with a Fourier series and can be optimized as well. Optimized finite-build coils for a W7-X-like stellarator are found, and compared with filamentary coil results. The final positions of optimized finite-build W7-X-like coils are shifted, on average, by approximately 2.5mm relative to optimized filamentary coils. These results suggest that finite-build effects should be accounted for in the optimization of stellarators with low coil tolerances.
△ Less
Submitted 30 November, 2020; v1 submitted 31 August, 2020;
originally announced September 2020.
-
Bifurcations of the magnetic axis and the alternating-hyperbolic sawtooth
Authors:
C. B. Smiet,
G. J. Kramer,
S. R. Hudson
Abstract:
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…
▽ More
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.
△ Less
Submitted 8 July, 2020;
originally announced July 2020.
-
Stepped pressure equilibrium with relaxed flow and applications in reversed-field pinch plasmas
Authors:
Z. S. Qu,
R. L. Dewar,
F. Ebrahimi,
J. K. Anderson,
S. R. Hudson,
M. J. Hole
Abstract:
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 Cod…
▽ More
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.
△ Less
Submitted 20 January, 2020;
originally announced January 2020.
-
Identification of important error fields in stellarators using Hessian matrix method
Authors:
Caoxiang Zhu,
David A. Gates,
Stuart R. Hudson,
Haifeng Liu,
Yuhong Xu,
Akihiro Shimizu,
Shoichi Okamura
Abstract:
Error fields are predominantly attributed to inevitable coil imperfections. Controlling error fields during coil fabrication and assembly is crucial for stellarators. Excessively tight coil tolerance increases time and cost, and, in part, led to the cancellation of NCSX and delay of W7-X. In this paper, we improve the recently proposed Hessian matrix method to rapidly identify important coil devia…
▽ More
Error fields are predominantly attributed to inevitable coil imperfections. Controlling error fields during coil fabrication and assembly is crucial for stellarators. Excessively tight coil tolerance increases time and cost, and, in part, led to the cancellation of NCSX and delay of W7-X. In this paper, we improve the recently proposed Hessian matrix method to rapidly identify important coil deviations. Two of the most common figures of merit, magnetic island size and quasi-symmetry, are analytically differentiated over coil parameters. By extracting the eigenvectors of the Hessian matrix, we can directly identify sensitive coil deviations in the order of the eigenvalues. The new method is applied to the upcoming CFQS configuration. Important perturbations that enlarge n/m=4/11 islands and deteriorate quasi-axisymmetry of the magnetic field are successfully determined. The results suggest each modular coil should have separate tolerance and some certain perturbation combinations will produce significant error fields. By relaxing unnecessary coil tolerance, this method will hopefully lead to a substantial reduction in time and cost.
△ Less
Submitted 4 September, 2019; v1 submitted 8 April, 2019;
originally announced April 2019.
-
Flow Resistance and Structures in Viscoelastic Channel Flows at Low Re
Authors:
Boyang Qin,
Paul F. Salipante,
Steven D. Hudson,
Paulo E. Arratia
Abstract:
The flow of viscoelastic fluids in channels and pipes remain poorly understood, particularly at low Reynolds numbers. Here, we investigate the flow of polymeric solutions in straight channels using pressure measurements and particle tracking. The law of flow resistance is established by measuring the flow friction factor $f_η$ versus flow rate. Two regimes are found: a transitional regime marked b…
▽ More
The flow of viscoelastic fluids in channels and pipes remain poorly understood, particularly at low Reynolds numbers. Here, we investigate the flow of polymeric solutions in straight channels using pressure measurements and particle tracking. The law of flow resistance is established by measuring the flow friction factor $f_η$ versus flow rate. Two regimes are found: a transitional regime marked by rapid increase in drag, and a turbulent-like regime characterized by a sudden decrease in drag and a weak dependence on flow rate. Lagrangian trajectories show finite transverse modulations not seen in Newtonian fluids. These curvature perturbations far downstream can generate sufficient hoop stresses to sustain the flow instabilities in the parallel shear flow.
△ Less
Submitted 2 July, 2018;
originally announced July 2018.
-
Non-planar elasticae as optimal curves for the magnetic axis of stellarators
Authors:
David Pfefferlé,
Lee Gunderson,
Stuart R. Hudson,
Lyle Noakes
Abstract:
The problem of finding an optimal curve for the target magnetic axis of a stellarator is addressed. Euler-Lagrange equations are derived for finite length three-dimensional curves that extremise their bending energy while yielding fixed integrated torsion. The obvious translational and rotational symmetry is exploited to express solutions in a preferred cylindrical coordinate system in terms of el…
▽ More
The problem of finding an optimal curve for the target magnetic axis of a stellarator is addressed. Euler-Lagrange equations are derived for finite length three-dimensional curves that extremise their bending energy while yielding fixed integrated torsion. The obvious translational and rotational symmetry is exploited to express solutions in a preferred cylindrical coordinate system in terms of elliptic Jacobi functions. These solution curves, which, up to similarity transformations, depend on three dimensionless parameters, do not necessarily close. Two closure conditions are obtained for the vertical and toroidal displacement (the radial coordinate being trivially periodic) to yield a countably infinite set of one-parameter families of closed non-planar curves. The behaviour of the integrated torsion (Twist of the Frenet frame), the Linking of the Frenet frame and the Writhe of the solution curves is studied in light of the \Calugareanu theorem. A refreshed interpretation of Mercier's formula for the on-axis rotational transform of stellarator magnetic field-lines is proposed.
△ Less
Submitted 6 September, 2018; v1 submitted 6 May, 2018;
originally announced May 2018.
-
First NuSTAR Limits on Quiet Sun Hard X-Ray Transient Events
Authors:
Andrew J. Marsh,
David M. Smith,
Lindsay Glesener,
Iain G. Hannah,
Brian W. Grefenstette,
Amir Caspi,
Säm Krucker,
Hugh S. Hudson,
Kristin K. Madsen,
Stephen M. White,
Matej Kuhar,
Paul J. Wright,
Steven E. Boggs,
Finn E. Christensen,
William W. Craig,
Charles J. Hailey,
Fiona A. Harrison,
Daniel Stern,
William W. Zhang
Abstract:
We present the first results of a search for transient hard X-ray (HXR) emission in the quiet solar corona with the \textit{Nuclear Spectroscopic Telescope Array} (\textit{NuSTAR}) satellite. While \textit{NuSTAR} was designed as an astrophysics mission, it can observe the Sun above 2~keV with unprecedented sensitivity due to its pioneering use of focusing optics. \textit{NuSTAR} first observed qu…
▽ More
We present the first results of a search for transient hard X-ray (HXR) emission in the quiet solar corona with the \textit{Nuclear Spectroscopic Telescope Array} (\textit{NuSTAR}) satellite. While \textit{NuSTAR} was designed as an astrophysics mission, it can observe the Sun above 2~keV with unprecedented sensitivity due to its pioneering use of focusing optics. \textit{NuSTAR} first observed quiet Sun regions on 2014 November 1, although out-of-view active regions contributed a notable amount of background in the form of single-bounce (unfocused) X-rays. We conducted a search for quiet Sun transient brightenings on time scales of 100 s and set upper limits on emission in two energy bands. We set 2.5--4~keV limits on brightenings with time scales of 100 s, expressed as the temperature T and emission measure EM of a thermal plasma. We also set 10--20~keV limits on brightenings with time scales of 30, 60, and 100 s, expressed as model-independent photon fluxes. The limits in both bands are well below previous HXR microflare detections, though not low enough to detect events of equivalent T and EM as quiet Sun brightenings seen in soft X-ray observations. We expect future observations during solar minimum to increase the \textit{NuSTAR} sensitivity by over two orders of magnitude due to higher instrument livetime and reduced solar background.
△ Less
Submitted 14 November, 2017;
originally announced November 2017.
-
Theory and discretization of ideal magnetohydrodynamic equilibria with fractal pressure profiles
Authors:
B. F. Kraus,
S. R. Hudson
Abstract:
In three-dimensional ideal magnetohydrodynamics, closed flux surfaces cannot maintain both rational rotational-transform and pressure gradients, as these features together produce unphysical, infinite currents. A proposed set of equilibria nullifies these currents by flattening the pressure on sufficiently wide intervals around each rational surface. Such rational surfaces exist at every scale, wh…
▽ More
In three-dimensional ideal magnetohydrodynamics, closed flux surfaces cannot maintain both rational rotational-transform and pressure gradients, as these features together produce unphysical, infinite currents. A proposed set of equilibria nullifies these currents by flattening the pressure on sufficiently wide intervals around each rational surface. Such rational surfaces exist at every scale, which characterizes the pressure profile as self-similar and thus fractal. The pressure profile is approximated numerically by considering a finite number of rational regions and analyzed mathematically by classifying the irrational numbers that support gradients into subsets. Applying these results to a given rotational-transform profile in cylindrical geometry, we find magnetic field and current density profiles compatible with the fractal pressure.
△ Less
Submitted 29 September, 2017;
originally announced September 2017.
-
New method to design stellarator coils without the winding surface
Authors:
Caoxiang Zhu,
Stuart R. Hudson,
Yuntao Song,
Yuanxi Wan
Abstract:
Finding an easy-to-build coils set has been a critical issue for stellarator design for decades. Conventional approaches assume a toroidal "winding" surface. We'll investigate if the existence of winding surface unnecessarily constrains the optimization, and a new method to design coils for stellarators is presented. Each discrete coil is represented as an arbitrary, closed, one-dimensional curve…
▽ More
Finding an easy-to-build coils set has been a critical issue for stellarator design for decades. Conventional approaches assume a toroidal "winding" surface. We'll investigate if the existence of winding surface unnecessarily constrains the optimization, and a new method to design coils for stellarators is presented. Each discrete coil is represented as an arbitrary, closed, one-dimensional curve embedded in three-dimensional space. A target function to be minimized that covers both physical requirements and engineering constraints is constructed. The derivatives of the target function are calculated analytically. A numerical code, named FOCUS, has been developed. Applications to a simple configuration, the W7-X, and LHD plasmas are presented.
△ Less
Submitted 8 August, 2017; v1 submitted 5 May, 2017;
originally announced May 2017.
-
Multi-region relaxed magnetohydrodynamics in plasmas with slowly changing boundaries --- resonant response of a plasma slab
Authors:
Robert L. Dewar,
Stuart R. Hudson,
Amitava Bhattacharjee,
Zensho Yoshida
Abstract:
The adiabatic limit of a recently proposed dynamical extension of Taylor relaxation, \emph{multi-region relaxed magnetohydrodynamics} (MRxMHD) is summarized, with special attention to the appropriate definition of relative magnetic helicity. The formalism is illustrated using a simple two-region, sheared-magnetic-field model similar to the Hahm--Kulsrud--Taylor (HKT) rippled-boundary slab model. I…
▽ More
The adiabatic limit of a recently proposed dynamical extension of Taylor relaxation, \emph{multi-region relaxed magnetohydrodynamics} (MRxMHD) is summarized, with special attention to the appropriate definition of relative magnetic helicity. The formalism is illustrated using a simple two-region, sheared-magnetic-field model similar to the Hahm--Kulsrud--Taylor (HKT) rippled-boundary slab model. In MRxMHD a linear Grad--Shafranov equation applies, even at finite ripple amplitude. The adiabatic switching on of boundary ripple excites a shielding current sheet opposing reconnection at a resonant surface. The perturbed magnetic field as a function of ripple amplitude is calculated by invoking conservation of magnetic helicity in the two regions separated by the current sheet. At low ripple amplitude "half islands" appear on each side of the current sheet, locking the rotational transform at the resonant value. Beyond a critical amplitude these islands disappear and the rotational transform develops a discontinuity across the current sheet.
△ Less
Submitted 17 March, 2017; v1 submitted 4 September, 2016;
originally announced September 2016.
-
Multi-region relaxed Hall magnetohydrodynamics with flow
Authors:
Manasvi Lingam,
Hamdi M. Abdelhamid,
Stuart R. Hudson
Abstract:
The recent formulations of multi-region relaxed magnetohydrodynamics (MRxMHD) have generalized the famous Woltjer-Taylor states by incorporating a collection of "ideal barriers" that prevent global relaxation, and flow. In this paper, we generalize MRxMHD with flow to include Hall effects (MRxHMHD), and thereby obtain the partially relaxed counterparts of the famous double Beltrami states as a spe…
▽ More
The recent formulations of multi-region relaxed magnetohydrodynamics (MRxMHD) have generalized the famous Woltjer-Taylor states by incorporating a collection of "ideal barriers" that prevent global relaxation, and flow. In this paper, we generalize MRxMHD with flow to include Hall effects (MRxHMHD), and thereby obtain the partially relaxed counterparts of the famous double Beltrami states as a special subset. The physical and mathematical consequences arising from the introduction of the Hall term are also presented. We demonstrate that our results (in the ideal MHD limit) constitute an important subset of ideal MHD equilibria, and we compare our approach against other variational principles proposed for deriving the partially relaxed states.
△ Less
Submitted 3 August, 2016; v1 submitted 18 July, 2016;
originally announced July 2016.
-
A rheological signature of frictional interactions in shear thickening suspensions
Authors:
John R. Royer,
Daniel L. Blair,
Steven D. Hudson
Abstract:
Colloidal shear thickening presents a significant challenge because the macroscopic rheology becomes increasingly controlled by the microscopic details of short ranged particle interactions in the shear thickening regime. Our measurements here of the first normal stress difference over a wide range of particle volume fraction elucidate the relative contributions from hydrodynamic lubrication and f…
▽ More
Colloidal shear thickening presents a significant challenge because the macroscopic rheology becomes increasingly controlled by the microscopic details of short ranged particle interactions in the shear thickening regime. Our measurements here of the first normal stress difference over a wide range of particle volume fraction elucidate the relative contributions from hydrodynamic lubrication and frictional contact forces, which have been debated. At moderate volume fractions we find $N_1<0$, consistent with hydrodynamic models, however at higher volume fractions and shear stresses these models break down and we instead observe dilation ($N_1>0$), indicating frictional contact networks. Remarkably, there is no signature of this transition in the viscosity, instead this change in the sign of $N_1$ occurs while the shear thickening remains continuous. These results suggest a scenario where shear thickening is driven primarily by the formation of frictional contacts, with hydrodynamic forces playing a supporting role at lower concentrations. Motivated by this picture, we introduce a simple model which combines these frictional and hydrodynamic contributions and accurately fits the measured viscosity over a wide range of particle volume fraction and shear stress.
△ Less
Submitted 1 March, 2016;
originally announced March 2016.
-
Variational formulation of relaxed and multi-region relaxed magnetohydrodynamics
Authors:
Robert L. Dewar,
Zensho Yoshida,
Amitava Bhattacharjee,
Stuart R. Hudson
Abstract:
Ideal magnetohydrodynamics (IMHD) is strongly constrained by an infinite number of microscopic constraints expressing mass, entropy and magnetic flux conservation in each infinitesimal fluid element, the latter preventing magnetic reconnection. By contrast, in the Taylor relaxation model for formation of macroscopically self-organized plasma equilibrium states, all these constraints are relaxed sa…
▽ More
Ideal magnetohydrodynamics (IMHD) is strongly constrained by an infinite number of microscopic constraints expressing mass, entropy and magnetic flux conservation in each infinitesimal fluid element, the latter preventing magnetic reconnection. By contrast, in the Taylor relaxation model for formation of macroscopically self-organized plasma equilibrium states, all these constraints are relaxed save for global magnetic fluxes and helicity. A Lagrangian variational principle is presented that leads to a new, fully dynamical, \emph{relaxed magnetohydrodynamics} (RxMHD), such that all static solutions are Taylor states but also allows flow. By postulating that some long-lived macroscopic current sheets can act as barriers to relaxation, separating the plasma into multiple relaxation regions, a further generalization, \emph{multi-region relaxed magnetohydrodynamics} (MRxMHD) is developed.
△ Less
Submitted 10 November, 2015; v1 submitted 1 September, 2015;
originally announced September 2015.
-
Multi-region relaxed magnetohydrodynamics with anisotropy and flow
Authors:
Graham R. Dennis,
Stuart R. Hudson,
Robert L. Dewar,
Matthew J. Hole
Abstract:
We present an extension of the multi-region relaxed magnetohydrodynamics (MRxMHD) equilibrium model that includes pressure anisotropy and general plasma flows. This anisotropic extension to our previous isotropic model is motivated by Sun and Finn's model of relaxed anisotropic magnetohydrodynamic equilibria. We prove that as the number of plasma regions becomes infinite, our anisotropic extension…
▽ More
We present an extension of the multi-region relaxed magnetohydrodynamics (MRxMHD) equilibrium model that includes pressure anisotropy and general plasma flows. This anisotropic extension to our previous isotropic model is motivated by Sun and Finn's model of relaxed anisotropic magnetohydrodynamic equilibria. We prove that as the number of plasma regions becomes infinite, our anisotropic extension of MRxMHD reduces to anisotropic ideal MHD with flow. The continuously nested flux surface limit of our MRxMHD model is the first variational principle for anisotropic plasma equilibria with general flow fields.
△ Less
Submitted 30 April, 2014;
originally announced May 2014.
-
Multi-region relaxed magnetohydrodynamics with flow
Authors:
G. R. Dennis,
S. R. Hudson,
R. L. Dewar,
M. J. Hole
Abstract:
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 wi…
▽ More
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.
△ Less
Submitted 3 March, 2014; v1 submitted 14 January, 2014;
originally announced January 2014.
-
A minimally constrained model of self-organized helical states in reversed-field pinches
Authors:
G. R. Dennis,
S. R. Hudson,
D. Terranova,
P. Franz,
R. L. Dewar,
M. J. Hole
Abstract:
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 bee…
▽ More
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.
△ Less
Submitted 18 June, 2013; v1 submitted 21 February, 2013;
originally announced February 2013.
-
The Infinite Interface Limit of Multiple-Region Relaxed MHD
Authors:
G. R. Dennis,
S. R. Hudson,
R. L. Dewar,
M. J. Hole
Abstract:
We show the stepped-pressure equilibria that are obtained from a generalization of Taylor relaxation known as multi-region, relaxed MHD (MRXMHD) are also generalizations of ideal MHD. We show this by proving that as the number of plasma regions becomes infinite, MRXMHD reduces to ideal MHD. Numerical convergence studies demonstrating this limit are presented.
We show the stepped-pressure equilibria that are obtained from a generalization of Taylor relaxation known as multi-region, relaxed MHD (MRXMHD) are also generalizations of ideal MHD. We show this by proving that as the number of plasma regions becomes infinite, MRXMHD reduces to ideal MHD. Numerical convergence studies demonstrating this limit are presented.
△ Less
Submitted 19 December, 2012;
originally announced December 2012.
-
Computation of multi-region relaxed magnetohydrodynamic equilibria
Authors:
S. R. Hudson,
R. L. Dewar,
G. Dennis,
M. J. Hole,
M. McGann,
G. von Nessi,
S. Lazerson
Abstract:
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 computa…
▽ More
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.
△ Less
Submitted 13 November, 2012;
originally announced November 2012.
-
Generalised action-angle coordinates defined on island chains
Authors:
Robert L. Dewar,
Stuart R. Hudson,
Ashley M. Gibson
Abstract:
Straight-field-line coordinates are very useful for representing magnetic fields in toroidally confined plasmas, but fundamental problems arise regarding their definition in 3-D geometries because of the formation of islands and chaotic field regions, ie non-integrability. In Hamiltonian dynamical systems terms these coordinates are a form of action-angle variables, which are normally defined only…
▽ More
Straight-field-line coordinates are very useful for representing magnetic fields in toroidally confined plasmas, but fundamental problems arise regarding their definition in 3-D geometries because of the formation of islands and chaotic field regions, ie non-integrability. In Hamiltonian dynamical systems terms these coordinates are a form of action-angle variables, which are normally defined only for integrable systems. In order to describe 3-D magnetic field systems, a generalisation of this concept was proposed recently by the present authors that unified the concepts of ghost surfaces and quadratic-flux-minimising (QFMin) surfaces. This was based on a simple canonical transformation generated by a change of variable $θ= θ(Θ,ζ)$, where $θ$ and $ζ$ are poloidal and toroidal angles, respectively, with $Θ$ a new poloidal angle chosen to give pseudo-orbits that are a) straight when plotted in the $ζ,Θ$ plane and b) QFMin pseudo-orbits in the transformed coordinate. These two requirements ensure that the pseudo-orbits are also c) ghost pseudo-orbits. In the present paper, it is demonstrated that these requirements do not \emph{uniquely} specify the transformation owing to a relabelling symmetry. A variational method of solution that removes this lack of uniqueness is proposed.
△ Less
Submitted 30 November, 2012; v1 submitted 1 April, 2012;
originally announced April 2012.
-
Nonaxisymmetric, multi-region relaxed magnetohydrodynamic equilibrium solutions
Authors:
S. R. Hudson,
R. L. Dewar,
M. J. Hole,
M. McGann
Abstract:
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 ar…
▽ More
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.
△ Less
Submitted 5 October, 2011; v1 submitted 26 July, 2011;
originally announced July 2011.
-
Action-gradient-minimizing pseudo-orbits and almost-invariant tori
Authors:
R. L. Dewar,
S. R. Hudson,
A. M. Gibson
Abstract:
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 surfa…
▽ More
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.
△ Less
Submitted 29 April, 2011;
originally announced April 2011.
-
Hamilton--Jacobi theory for continuation of magnetic field across a toroidal surface supporting a plasma pressure discontinuity
Authors:
M. McGann,
S. R. Hudson,
R. L. Dewar,
G. von Nessi
Abstract:
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…
▽ More
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.
△ Less
Submitted 18 February, 2010;
originally announced February 2010.
-
Unified Theory of Ghost and Quadratic-Flux-Minimizing Surfaces
Authors:
R. L. Dewar,
S. R. Hudson,
A. M. Gibson
Abstract:
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…
▽ More
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.
△ Less
Submitted 2 August, 2011; v1 submitted 4 January, 2010;
originally announced January 2010.
-
Are ghost surfaces quadratic-flux-minimizing?
Authors:
S. R. Hudson,
R. L. Dewar
Abstract:
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-orbit…
▽ More
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.
△ Less
Submitted 2 January, 2010; v1 submitted 11 September, 2009;
originally announced September 2009.
-
Relaxed MHD states of a multiple region plasma
Authors:
M. J. Hole,
R. Mills,
S. R. Hudson,
R. L. Dewar
Abstract:
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 stabil…
▽ More
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.
△ Less
Submitted 19 February, 2009;
originally announced February 2009.
-
Relaxed plasma equilibria and entropy-related plasma self-organization principles
Authors:
R. L. Dewar,
M. J. Hole,
M. McGann,
R. Mills,
S. R. Hudson
Abstract:
The concept of plasma relaxation as a constrained energy minimization is reviewed. Recent work by the authors on generalizing this approach to partially relaxed three-dimensional plasma systems in a way consistent with chaos theory is discussed, with a view to clarifying the thermodynamic aspects of the variational approach used. Other entropy-related approaches to finding long-time steady state…
▽ More
The concept of plasma relaxation as a constrained energy minimization is reviewed. Recent work by the authors on generalizing this approach to partially relaxed three-dimensional plasma systems in a way consistent with chaos theory is discussed, with a view to clarifying the thermodynamic aspects of the variational approach used. Other entropy-related approaches to finding long-time steady states of turbulent or chaotic plasma systems are also briefly reviewed.
△ Less
Submitted 5 October, 2008;
originally announced October 2008.