A neoclassically optimized compact stellarator with simple coils has been designed. The magnetic field of the new stellarator is generated by only four planar coils including two interlocking coils of elliptical shape and two circular poloidal field coils. The interlocking coil topology is the same as that of the Columbia Non-neutral Torus (CNT). The new configuration was obtained by minimizing the effective helical ripple directly via the shape of the two interlocking coils. The optimized compact stellarator has very low effective ripple in the plasma core implying excellent neoclassical confinement. This is confirmed by the results of the drift-kinetic code SFINCS showing that the particle diffusion coefficient of the new configuration is one order of magnitude lower than CNTs.
Magnetic confinement devices for nuclear fusion can be large and expensive. Compact stellarators are promising candidates for costreduction, but introduce new difficulties: confinement in smaller volumes requires higher magnetic field, which calls for higher coil-currents and ultimately causes higher Laplace forces on the coils-if everything else remains the same. This motivates the inclusion of force reduction in stellarator coil optimization. In the present paper we consider a coil winding surface, we prove that there is a natural and rigorous way to define the Laplace force (despite the magnetic field discontinuity across the current-sheet), we provide examples of cost associated (peak force, surface-integral of the force squared) and discuss easy generalizations to parallel and normal force-components, as these will be subject to different engineering constraints. Such costs can then be easily added to the figure of merit in any multi-objective stellarator coil optimization code. We demonstrate this for a generalization of the REGCOIL code [1], which we rewrote in python, and provide numerical examples for the NCSX (now QUASAR) design. We present results for various definitions of the cost function, including peak force reductions by up to 40 %, and outline future work for further reduction.
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 distance from the magnetic axis. Our construction provides a parameterization of the space of such solutions, and the cylindrical reformulation and numerical method of Landreman and Sengupta (2018); Landreman et al. (2019), enables their efficient numerical construction.
We consider a problem relating to magnetic confinement devices known as stellarators. Plasma is confined by magnetic fields generated by current-carrying coils, and here we investigate how closely to the plasma they need to be positioned. Current-carrying coils are represented as singularities within the magnetic field and therefore this problem can be modelled mathematically as finding how far we can harmonically extend a vector field from the boundary of a domain. For this paper we consider two-dimensional domains with real analytic boundary, and prove that a harmonic extension exists if and only if the boundary data satisfies a combined compatibility and regularity condition. Our method of proof uses a generalisation of a result of Hadamard on the Cauchy problem for the Laplacian. We then provide a lower bound on how far we can harmonically extend the vector field from the boundary via the Cauchy--Kovalevskaya Theorem.
The nonlinear gyrokinetic code GS2 has been extended to treat non-axisymmetric stellarator geometry. Electromagnetic perturbations and multiple trapped particle regions are allowed. Here, linear, collisionless, electrostatic simulations of the quasi-axisymmetric, three-field period National Compact Stellarator Experiment (NCSX) design QAS3-C82 have been successfully benchmarked against the eigenvalue code FULL. Quantitatively, the linear stability calculations of GS2 and FULL agree to within ~10%.
One metric for comparing confinement properties of different magnetic fusion energy configurations is the linear critical gradient of drift wave modes. The critical gradient scale length determines the ratio of the core to pedestal temperature when a plasma is limited to marginal stability in the plasma core. The gyrokinetic turbulence code GS2 was used to calculate critical temperature gradients for the linear, collisionless ion temperature gradient (ITG) mode in the National Compact Stellarator Experiment (NCSX) and a prototypical shaped tokamak, based on the profiles of a JET H-mode shot and the stronger shaping of ARIES-AT. While a concern was that the narrow cross section of NCSX at some toroidal locations would result in steep gradients that drive instabilities more easily, it is found that other stabilizing effects of the stellarator configuration offset this so that the normalized critical gradients for NCSX are competitive with or even better than for the tokamak. For the adiabatic ITG mode, NCSX and the tokamak had similar critical gradients, though beyond marginal stability, NCSX had larger growth rates. However, for the kinetic ITG mode, NCSX had a higher critical gradient and lower growth rates until a/L_T is approximately 1.5 times a/L_{T,crit}, when it surpassed the tokamaks. A discussion of the results presented with respect to a/L_T vs. R/L_T is included.