Quasi-axisymmetric magnetic fields: weakly non-axisymmetric case in a vacuum
An asymptotic expansion is performed to obtain quasi-axisymmetric magnetic configurations that are weakly non-axisymmetric. A large space of solutions is identified, which satisfy the condition of quasi-axisymmetry on a single magnetic flux surface, while (non-axisymmetric) globally quasi-axisymmetric solutions are shown to not exist, agreeing with the conclusions of previous theoretical work. The solutions found are shown to be geometrically constrained at low aspect ratio or high toroidal period number. Solutions satisfying the more general condition of omnigeneity (generalized quasi-axisymmetry) are also shown to exist, and it is found that quasi-axisymmetric deformations can be superposed with an omnigenous solution, while preserving the property of omnigeneity, effectively extending the space of "good" configurations. A numerical solution of the first order quasi-axisymmetry problem is demonstrated and compared with solutions found with a widely used MHD equilibrium solver, independently verifying that quasi-axisymmetry is satisfied at the appropriate order. It is thereby demonstrated that approximately quasi-axisymmetric solutions can be directly constructed, i.e. without using numerical search algorithms.
Constructing stellarators with quasisymmetry to high order
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 & Boozer’s expansion about the magnetic axis (Phys Fluids B, vol. 3, 1991, p. 2805), and it is several orders of magnitude faster than the conventional optimization approach. The work here extends the method of Landreman et al. (J. Plasma Phys., vol. 85, 2019, 905850103), which was limited to flux surfaces with elliptical cross-section, to higher order in the aspect-ratio expansion. As a result, configurations can be generated with strong shaping that achieve quasisymmetry to high accuracy. Using this construction, we give the first numerical demonstrations of Garren and Boozer’s ideal scaling of quasisymmetry breaking with the cube of the inverse aspect ratio. We also demonstrate a strongly non-axisymmetric configuration (vacuum rotational transform $\unicode[STIX]{x1D704}>0.4$ ) in which symmetry-breaking mode amplitudes throughout a finite volume are ${<}2\times 10^{-7}$ , the smallest ever reported. To generate boundary shapes of finite-minor-radius configurations, a careful analysis is given of the effect of substituting a finite minor radius into the near-axis expansion. The approach here can provide analytic insight into the space of possible quasisymmetric stellarator configurations, and it can be used to generate good initial conditions for conventional stellarator optimization.
Existence of three‐dimensional toroidal MHD equilibria with nonconstant pressure
We establish an existence result for the three‐dimensional MHD equations (∇ x B) x B = ∇p, ∇·B = 0, B·n|∂T = 0 with p ≠ const in tori T without symmetry. More precisely, our theorems insure the existence of sharp boundary solutions for tori whose departure from axisymmetry is sufficiently small; they allow for solutions to be constructed with an arbitrary number of pressure jumps.