Explore the vast range of titles available.

The confinement of the guiding-centre trajectories in a stellarator is determined by the variation of the magnetic field strength $B$ in Boozer coordinates $(r,\\unicode[STIX]{x1D703},\\unicode[STIX]{x1D711})$ , but $B(r,\\unicode[STIX]{x1D703},\\unicode[STIX]{x1D711})$ depends on the flux surface shape in a complicated way. Here we derive equations relating $B(r,\\unicode[STIX]{x1D703},\\unicode[STIX]{x1D711})$ in Boozer coordinates and the rotational transform to the shape of flux surfaces in cylindrical coordinates, using an expansion in distance from the magnetic axis. A related expansion was done by Garren and Boozer (Phys. Fluids B, vol. 3, 1991a, 2805) based on the Frenet–Serret frame, which can be discontinuous anywhere the magnetic axis is straight, a situation that occurs in the interesting case of omnigenity with poloidally closed $B$ contours. Our calculation in contrast does not use the Frenet–Serret frame. The transformation between the Garren–Boozer approach and cylindrical coordinates is derived, and the two approaches are shown to be equivalent if the axis curvature does not vanish. The expressions derived here help enable optimized plasma shapes to be constructed that can be provided as input to VMEC and other stellarator codes, or to generate initial configurations for conventional stellarator optimization.

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.

Steady plasma flows have been studied almost exclusively in systems with continuous symmetry or in open domains. In the absence of continuous symmetry, the lack of a conserved quantity makes the study of flows intrinsically challenging. In a toroidal domain, the requirement of double periodicity for physical quantities adds to the complications. In particular, the magnetohydrodynamics (MHD) model of plasma steady state with the flow in a non-symmetric toroidal domain allows the development of singularities when the rotational transform of the magnetic field is rational, much like the equilibrium MHD model. In this work, we show that steady flows can still be maintained provided the rotational transform is close to rational and the magnetic shear is weak. We extend the techniques developed in carrying out perturbation methods to all orders for static MHD equilibrium by Weitzner (Phys. Plasmas, vol. 21, 2014, p. 022515) to MHD equilibrium with flows. We construct perturbative MHD equilibrium in a doubly periodic domain with nearly parallel flows by systematically eliminating magnetic resonances order by order. We then utilize an additional symmetry of the flow problem, first discussed by Hameiri (J. Math. Phys., vol. 22, 1981, pp. 2080–2088, § III), to obtain a generalized Grad–Shafranov equation for a class of non-symmetric three-dimensional MHD equilibrium with flows both parallel and perpendicular to the magnetic field. For this class of flows, we can obtain non-symmetric generalizations of integrals of motion, such as Bernoulli's function and angular momentum. Finally, we obtain the generalized Hamada conditions, which are necessary to suppress singular currents in such a system when the magnetic field lines are closed. We do not attempt to address the question of neoclassical damping of flows.

Quasisymmetric stellarators are appealing intellectually and as fusion reactor candidates since the guiding-centre 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 minimize symmetry-breaking Fourier modes of the field strength $B$ . Here, instead, we directly construct magnetic fields in cylindrical coordinates that are quasisymmetric to leading order in the distance from the magnetic axis, without using optimization. The method involves solution of a one-dimensional nonlinear ordinary differential equation, originally derived by Garren & Boozer (Phys. Fluids B, vol. 3, 1991, p. 2805). We demonstrate the usefulness and accuracy of this optimization-free approach by providing the results of this construction as input to the codes VMEC and BOOZ_XFORM, confirming the purity and scaling of the magnetic spectrum. The space of magnetic fields that are quasisymmetric to this order is parameterized by the magnetic axis shape along with three other real numbers, one of which reflects the on-axis toroidal current density, and another one of which is zero for stellarator symmetry. The method here could be used to generate good initial conditions for conventional optimization, and its speed enables exhaustive searches of parameter space.

Stellarators are generically small current and low plasma beta devices. Often the construction of vacuum magnetic fields with good magnetic surfaces is the starting point for an equilibrium calculation. Although in cases with some continuous spatial symmetry, flux functions can always be found for vacuum magnetic fields, an analogous function does not, in general, exist in three dimensions. This work examines several simple equilibria and vacuum magnetic field problems with the intent of demonstrating the possibilities and limitations in the construction of such states. Starting with a simple vacuum magnetic field with closed field lines in a topological torus (toroidal shell with a flat metric), we obtain a self-consistent formal perturbation series using the amplitude of the non-symmetric vacuum fields as a small parameter. We show that systems possessing stellarator symmetry allow the construction order by order. We further indicate the significance of stellarator symmetry in the amplitude expansion of the full ideal magnetohydrodynamics (MHD) problem as well. We then investigate the conditions that guarantee neighbouring flux surfaces given the data on one surface, by expanding in the distance from that surface. We show that it is much more difficult to find low shear vacuum fields with surfaces than force-free fields or ideal MHD equilibrium. Finally, we demonstrate the existence of a class of vacuum magnetic fields, analogous to ‘snakes’ observed in tokamaks, which can be expanded to all orders.

A first-order model is derived for quasisymmetric stellarators where the vacuum field due to coils is dominant, but plasma-current-induced terms are not negligible and can contribute to magnetic differential equations, with $\\beta$ of the order of the ratio induced to vacuum fields. Under these assumptions, it is proven that the aspect ratio must be large and a simple expression can be obtained for the lowest-order vacuum field. The first-order correction, which involves both vacuum and current-driven fields, is governed by a Grad–Shafranov equation and the requirement that flux surfaces exist. These two equations are not always consistent, and so this model is generally overconstrained, but special solutions exist that satisfy both equations simultaneously. One family of such solutions is the set of first-order near-axis solutions. Thus, the first-order near-axis model is a subset of the model presented here. Several other solutions outside the scope of the near-axis model are also found. A case study comparing one such solution to a VMEC-generated solution shows good agreement.

A systematic theory of the asymptotic expansion of the magnetohydrostatics (MHS) equilibrium in the distance from the magnetic axis is developed to include arbitrary smooth currents near the magnetic axis. Compared with the vacuum and the force-free system, an additional magnetic differential equation must be solved to obtain the pressure-driven currents. It is shown that there exist variables in which the rest of the MHS system closely mimics the vacuum system. Thus, a unified treatment of MHS fields is possible. The mathematical structure of the near-axis expansions to arbitrary order is examined carefully to show that the double-periodicity of physical quantities in a toroidal domain can be satisfied order by order. The essential role played by the leading-order Birkhoff–Gustavson normal form in solving the magnetic differential equations is highlighted. Several explicit examples of vacuum, force-free and MHS equilibrium in different geometries are presented.

We carry out expansions of non-symmetric toroidal ideal magnetohydrodynamic (MHD) equilibria with nested flux surfaces about a periodic cylinder, in which physical quantities are periodic of period $2\\unicode[STIX]{x03C0}$ in the cylindrical angle $\\unicode[STIX]{x1D703}$ and $z$ . The cross-section of a flux surface at a constant toroidal angle is assumed to be approximately circular, and data are given on the cylindrical flux surface $r=1$ . Furthermore, we assume that the magnetic field lines are closed on the lowest-order flux surface, and the magnetic shear is relatively small. We extend earlier work in a flat torus by Weitzner (Phys. Plasmas, vol. 23, 2016, 062512) and demonstrate that a power series expansion can be carried out to all orders using magnetic flux as an expansion parameter. The cylindrical metric introduces certain new features to the expansions compared to the flat torus. However, the basic methodology of dealing with resonance singularities remains the same. The results, even though lacking convergence proofs, once again support the possibility of smooth, low-shear non-symmetric toroidal MHD equilibria.

Quasisymmetry (QS), a hidden symmetry of the magnetic field strength, is known to support nested flux surfaces and provide superior particle confinement in stellarators. In this work, we study the ideal magnetohydrodynamic (MHD) equilibrium and stability of high-beta plasma in a large-aspect-ratio stellarator. In particular, we show that the lowest-order description of a near-axisymmetric equilibrium vastly simplifies the problem of three-dimensional quasisymmetric MHD equilibria, which can be reduced to a standard elliptic Grad–Shafranov equation for the flux function. We show that any large-aspect-ratio tokamak, deformed periodically in the vertical direction, is a stellarator with approximate volumetric QS. We discuss exact analytical solutions and numerical benchmarks. Finally, we discuss the ideal ballooning and interchange stability of some of our equilibrium configurations.

A systematic theory of the asymptotic expansion of the magnetohydrostatics (MHS) equilibrium in the distance from the magnetic axis is developed to include arbitrary smooth currents near the magnetic axis. Compared with the vacuum and the force-free system, an additional magnetic differential equation must be solved to obtain the pressure-driven currents. It is shown that there exist variables in which the rest of the MHS system closely mimics the vacuum system. Thus, a unified treatment of MHS fields is possible. The mathematical structure of the near-axis expansions to arbitrary order is examined carefully to show that the double-periodicity of physical quantities in a toroidal domain can be satisfied order by order. The essential role played by the leading-order Birkhoff–Gustavson normal form in solving the magnetic differential equations is highlighted. Several explicit examples of vacuum, force-free and MHS equilibrium in different geometries are presented.