The confinement of guiding center trajectories in a stellarator is determined by the variation of the magnetic field strength B in Boozer coordinates (r,θ,φ), but B(r,θ,φ) depends on the flux surface shape in a complicated way. Here we derive equations relating B(r,θ,φ) in Boozer coordinates and the…

Recently designed optimized stellarator experiments have suffered from very tight construction tolerances, but some level of deviation of the coil system is unavoidable during fabrication of the coils and assembly of the coil system. In this paper, we present a new approach that incorporates reduced sensitivity to construction tolerances of…

A method is given to rapidly compute quasisymmetric stellarator magnetic fields for plasma confinement, without the need to call a three-dimensional magnetohydrodynamic equilibrium code inside an optimization iteration. The method is based on direct solution of the equations of magnetohydrodynamic equilibrium and quasisymmetry using Garren and…

This paper describes a new and efficient method of defining an annular region of a curl-free magnetic field with specific physics and coil properties that can be used in stellarator design. Three statements define the importance: (1) Codes can follow an optimized curl-free initial state to a final full-pressure equilibrium. The large size of…

Quasisymmetric stellarators are appealing intellectually and as fusion reactor candidates since the guiding center particle trajectories and neoclassical transport are isomorphic to those in a tokamak, implying good confinement. Previously, quasisymmetric magnetic fields have been identified by applying black-box optimization algorithms to…

The condition of omnigenity is investigated, and applied to the near-axis expansion of Garren and Boozer (1991a). Due in part to the particular analyticity requirements of the near-axis expansion, we find that, excluding quasi-symmetric solutions, only one type of omnigenity, namely quasi-isodynamicity, can be satisfied at first order in the…

A method is presented to establish regions of phase space for 3D vector fields through which pass no co-oriented invariant 2D submanifolds transverse to a given oriented 1D foliation. Refinements are given for the cases of volume-preserving or Cartan–Arnol’d Hamiltonian flows and for boundaryless submanifolds.

We revisit the Hahm-Kulsrud-Taylor (HKT) problem, a classic prototype problem for studying resonant magnetic perturbations and 3D magnetohydrodynamical equilibria. We employ the boundary-layer techniques developed by Rosenbluth, Dagazian, and Rutherford (RDR) for the internal m=1 kink instability, while addressing the subtle…

It is shown that the resistive magnetohydrodynamic stability of a slab force-free current sheet can be calculated using the variational principle of multi-region relaxed magnetohydrodynamics and that the corresponding stability boundary is in exact agreement with linear tearing mode theory.

A direct construction of equilibrium magnetic fields with toroidal topology at arbitrary order in the distance from the magnetic axis is carried out, yielding an analytical framework able to explore the landscape of possible magnetic flux surfaces in the vicinity of the axis. This framework can provide meaningful analytical insight on the…

Most quasisymmetric stellarators to date have been designed by numerically optimizing the plasma boundary shape to minimize symmetry-breaking Fourier modes of the magnetic field strength B. At high aspect ratio, a faster approach is to directly construct the plasma shape from the equations of quasisymmetry near the magnetic axis derived…

We describe a method to construct smooth and compactly supported solutions of 3D incompressible Euler equations and related models. The method is based on localizable Grad-Shafranov equations and is inspired by the recent result.

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…

Analytic scaling relations are derived for a phenomenological model of the plasmoid instability in an evolving current sheet, including the effects of reconnection outflow. Two scenarios are considered, where the plasmoid instability can be triggered either by an injected initial perturbation or by the natural noise of the system (here…

In this work, we consider optimal control problems constrained by elliptic partial differential equations (PDEs) with lognormal random coefficients, which are represented by a countably infinite-dimensional random parameter with i.i.d. normal distribution. We approximate the optimal solution by a suitable truncation of its Hermite polynomial…

We present a boundary integral equation solver for computing Taylor relaxed states in non-axisymmetric solid and shell-like toroidal geometries. The computation of Taylor states in these geometries is a key element for the calculation of stepped pressure stellarator equilibria. The integral representation of the magnetic field in this work is…

Differential forms provide a coordinate-free way to express many quantities and relations in mathematical physics. In particular, they are useful in plasma physics. This tutorial gives a guide so that you can read the plasma physics literature that uses them and apply them yourself.

It is shown that the variational principle of multi-region relaxed magnetohydrodynamics (MRxMHD) can be used to predict the stability and nonlinear saturation of tearing modes in strong guide field configurations without resolving the dynamics and without explicit dependence on the plasma resistivity. While the magnetic helicity is not a good…