In this self-contained document, we aim to present the basic theoretical building blocks to understand modeling of stellarator magnetic fields, some of the challenges associated with modeling, and optimization for designing stellarators. As often as possible, the ideas will be presented using equations and pictures, and references to other relevant introductory material will be included. This document is accessible to those who may not have a physics background but are interested in applications of mathematical and computational tools to stellarator research.

Small-amplitude, symmetry-breaking magnetic field perturbations, including resonant magnetic perturbations (RMPs) and error fields, can profoundly impact plasma properties in both tokamaks and stellarators. In this work, we perform the first comparison between the Stepped Pressure Equilibrium Code (SPEC) (a comparatively fast and efficient equilibrium code based on energy-minimisation principles) and M3D-C $^{1}$ (a high-fidelity albeit computationally expensive initial-value extended-magnetohydro- dynamic (MHD) code) to assess the conditions under which SPEC can be used to model the nonlinear, non-ideal plasma response to an externally applied $(m=2,n=1)$ RMP field in an experimentally relevant geometry. We find that SPEC is able to capture the plasma response in the weakly nonlinear regime – meaning perturbation amplitudes below the threshold for break up of the separatrix and onset of secondary magnetic island formation – when around half of the total toroidal flux is enclosed in the volume containing the $q=2$ resonant surface. The observed dependence of SPEC solutions on input parameters, including toroidal flux and the number of volumes into which the plasma is partitioned, indicates that additional exploration of the underlying Multi-Region Relaxed MHD physics model is needed to constrain the choice of parameters. Nonetheless, this work suggests promising applications of SPEC to optimisation and fusion plasma design.

We present recent highlights from the most recent operation phases of Wendelstein 7-X, the most advanced stellarator in the world. Stable detachment with good particle exhaust, low impurity content, and energy confinement times exceeding 100 ms, have been maintained for tens of seconds. Pellet fueling allows for plasma phases with reduced ion-temperature-gradient turbulence, and during such phases, the overall confinement is so good (energy confinement times often exceeding 200 ms) that the attained density and temperature profiles would not have been possible in less optimized devices, since they would have had neoclassical transport losses exceeding the heating applied in W7-X. This provides proof that the reduction of neoclassical transport through magnetic field optimization is successful. W7-X plasmas generally show good impurity screening and high plasma purity, but there is evidence of longer impurity confinement times during turbulence-suppressed phases.

Over the last decade, a variational principle based on a generalisation of Taylor’s relaxation, referred to as multi-region relaxed magnetohydrodynamics (MRxMHDs) has been developed. The numerical solutions of the MRxMHD equilibria have been constructed using the Stepped Pressure Equilibrium Code (SPEC) (Hudson et al 2012 Phys. Plasmas 19 112502). In principle, SPEC could also be established to describe the MRxMHD stability of an equilibrium. In this work, a theoretical framework is developed to relate the second variation of the energy functional to the so-called Hessian matrix, enabling the prediction of MHD linear instabilities of cylindrical plasmas, and is implemented in SPEC. The negative and positive eigenvalues of the Hessian matrix predict the stability of an equilibrium. Verification studies of SPEC stability results with the M3D-C 1 code and the tearing mode Δ′ criterion have been conducted for ideal and resistive MHD instabilities, respectively, in a pressureless cylindrical tokamak, and have shown good agreement. Our stability analysis is a critical step towards understanding the MHD stability of three-dimensional MHDs where nested flux surfaces, magnetic islands and stochastic regions co-exist.

Equilibria with extended regions of weak magnetic shear, including some tokamak scenarios and stellarators, can be susceptible to pressure-driven internal MHD instabilities even though there is no mode rational surface in the plasma. Nonresonant modes, in particular, can have properties that are unattractive for confinement, including displacing substantial volumes of the plasma and leading to more efficient pressure gradient flattening in the nonlinear regime. The onset and linear properties of the low shear stability regime are examined using the initial-value, extended-MHD code M3D-C1. For monotonic q-profiles, we demonstrate a clear correlation between the convergents associated with the continued fraction representation of q0 and the spectrum of unstable modes. Nonresonant modes are observed to be destabilized preferentially to any other resonant instability with the same toroidal mode number when n > 1. Using the observed connection between the spectrum associated with q0 and the overall equilibrium stability properties, we suggest a technique for reducing the uncertainty on both q0 and magnetic shear in the core region, obtained either via measurement or through the analysis and reconstruction of experimental results.

Nonresonant internal modes can be difficult to anticipate as there is no resonant surface in the plasma. However, equilibria that are unstable to multiple nonresonant magnetohydrodynamic (MHD) modes may be more prone to global loss of confinement since these instabilities generate spatially extended linear displacements, potentially enhancing magnetic field line chaos via nonlinear interactions. Here, we successfully predict the unstable nonresonant pressure-driven modes for equilibria with zero shear in the plasma core, irrational q on axis, and a central pressure gradient, which is consistent with pre-crash profiles in sawtoothing tokamak plasmas in the large-aspect-ratio limit. A criterion for identifying nonresonant modes most likely to be unstable is developed from the convergents of the continued fraction representation of q0. A higher-order analysis of the standard Energy Principle reveals the conditions under which these modes are expected to dominate. Linear growth rate spectra, as a function of toroidal mode number (up to n = 30), calculated using the initial-value extended-MHD code, M3D-C1, recover the characteristic dependence observed for ideal infernal modes. Nonresonant modes have also been invoked in some ideal sawtooth crash models. This work provides a mechanism to predict the mode numbers of infernal modes and, potentially, the width of some post-sawtooth-crash profiles.

A new detector array, Solenogam, has been developed at the Australian National University Heavy Ion Accelerator Facility. Coupled initially to the SOLITAIRE 6.5 T, gas-filled, solenoidal separator, and later to an 8 T solenoid, the system enables the study of long-lived nuclear states through γ-ray and conversion-electron spectroscopy in a low-background environment. The detector system is described and results from the commissioning experiments are presented.

In this self-contained document, we aim to present the basic theoretical building blocks to understand modeling of stellarator magnetic fields, some of the challenges associated with modeling, and optimization for designing stellarators. As often as possible, the ideas will be presented using equations and pictures, and references to other relevant introductory material will be included. This document is accessible to those who may not have a physics background but are interested in applications of mathematical and computational tools to stellarator research.

As a step toward understanding 3D magnetohydrodynamic (MHD) equilibria, for which smooth solutions may not exist, we develop a simple cylindrical model to investigate the resistive stability of MHD equilibria with alternating regions of constant and nonuniform pressure, producing states with continuous total pressure (i.e., no singular current sheets) but discontinuities in the parallel current density. We examine how the resistive stability characteristics of the model change as we increase the localization of pressure gradients at fixed radii, which approaches a discontinuous pressure profile in the zero-width limit. Equilibria with continuous pressure are found to be unstable to moderate/high-m modes and apparently tend toward ideal instability in some cases. We propose that additional geometric degrees of freedom or symmetry breaking via island formation may increase the parameter space on which equilibria of our model are physically realizable, while preserving the radial localization of pressure gradients. This is consistent with the possibility of realizing, in practice, 3D MHD equilibria which support both continuously nested flux surfaces (where ∇p ≠ 0) and chaotic field regions (where ∇p = 0).

A steady-state single fluid MHD model which describes the equilibrium of plasma parameters in a collisional, rotating plasma column with temperature gradients and a non-uniform externally applied magnetic field is developed. Two novel methods of simplifying the governing equations are introduced. Specifically, a 'radial transport constraint' and an ordering argument are applied. The reduced system is subsequently solved to yield the equilibrium of macroscopic plasma parameters in the bulk region of the plasma. The model is benchmarked by comparing these solutions to experimental measurements of axial velocity and density for a hydrogen plasma in the converging-field experiment MAGPIE and overall a good agreement is observed. The plasma equilibrium is determined by the interaction of a density gradient, due to a temperature gradient, with an electric field. The magnetic field and temperature gradient are identified as key parameters in determining the flow profile, which may be important considerations in other applications.

Using the Er-168(B-10,5n) reaction at a beam energy of 68 MeV, new data have been obtained for the population and decay of a T-1/2 = 148 ns, K-pi = 21/2(-) three-quasiparticle isomer at 1717 keV in Ta-173. Revised decay energies and intensities have been determined, together with newly observed members of a rotational band associated with the isomer. By comparison with other isomers in the A approximate to 180 deformed region, the 173Ta isomer properties help to specify the key degrees of freedom that determine K- forbidden transition rates. In particular, when all three quasiparticles are of the same nucleon type, there is a strong dependence of the E2 reduced hindrance factor on the isomer excitation energy.

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