Showing 1–50 of 89 results for author: Landreman, M

Flexible Stellarator Physics Facility

Authors: F. I. Parra, S. -G. Baek, M. Churchill, D. R. Demers, B. Dudson, N. M. Ferraro, B. Geiger, S. Gerhardt, K. C. Hammond, S. Hudson, R. Jorge, E. Kolemen, D. M. Kriete, S. T. A. Kumar, M. Landreman, C. Lowe, D. A. Maurer, F. Nespoli, N. Pablant, M. J. Pueschel, A. Punjabi, J. A. Schwartz, C. P. S. Swanson, A. M. Wright

Abstract: We propose to build a Flexible Stellarator Physics Facility to explore promising regions of the vast parameter space of disruption-free stellarator solutions for Fusion Pilot Plants (FPPs). We propose to build a Flexible Stellarator Physics Facility to explore promising regions of the vast parameter space of disruption-free stellarator solutions for Fusion Pilot Plants (FPPs). △ Less

Submitted 4 July, 2024; originally announced July 2024.

Comments: White paper submitted to FESAC subcommittee on Facilities, 8 pages

A family of quasi-axisymmetric stellarators with varied rotational transform

Authors: Stefan Buller, Matt Landreman, John Kappel, Rahul Gaur

Abstract: We apply a continuation method to recently optimized stellarator equilibria with excellent quasi-axisymmetry (QA) to generate new equilibria with a wide range of rotational transform profiles. Using these equilibria, we investigate how the rotational transform affects fast-particle confinement, the maximum coil-plasma distance, the maximum growth rate in linear gyrokinetic ion-temperature gradient… ▽ More We apply a continuation method to recently optimized stellarator equilibria with excellent quasi-axisymmetry (QA) to generate new equilibria with a wide range of rotational transform profiles. Using these equilibria, we investigate how the rotational transform affects fast-particle confinement, the maximum coil-plasma distance, the maximum growth rate in linear gyrokinetic ion-temperature gradient (ITG) simulations, and the ion heat flux in corresponding nonlinear simulations. We find values of two-term quasisymmetry error comparable to or lower than the similar Landreman-Paul (Phys. Rev. Lett. 128, 035001) configuration for values of the mean rotational transform $\barι$ between $0.12$ and $0.75$. The fast-particle confinement improves with $\barι$ until $\barι = 0.73$, at which point the degradation in quasisymmetry outweighs the benefits of further increasing $\barι$. The required coil-plasma distance only varies by about $\pm 10\%$ for the configurations under consideration, and is between $2.8\,\mathrm{m}$ to $3.3\,\mathrm{m}$ when the configuration is scaled up to reactor size. The maximum growth rate from linear gyrokinetic simulations increases with $\barι$, but also shifts towards higher $k_y$ values. The maximum linear growth rate is sensitive to the choice of flux tube at rational $ι$, but this can be compensated for by taking the maximum over several flux tubes. The corresponding ion heat fluxes from nonlinear simulations display a non-monotonic relation to $ι$. Sufficiently large positive shear is destabilizing. This is reflected both in linear growth rates and nonlinear heat fluxes. △ Less

Submitted 17 January, 2024; originally announced January 2024.

Comments: 20 pages, 15 figures

MONKES: a fast neoclassical code for the evaluation of monoenergetic transport coefficients

Authors: F. J. Escoto, J. L. Velasco, I. Calvo, M. Landreman, F. I. Parra

Abstract: MONKES is a new neoclassical code for the evaluation of monoenergetic transport coefficients in stellarators. By means of a convergence study and benchmarks with other codes, it is shown that MONKES is accurate and efficient. The combination of spectral discretization in spatial and velocity coordinates with block sparsity allows MONKES to compute monoenergetic coefficients at low collisionality,… ▽ More MONKES is a new neoclassical code for the evaluation of monoenergetic transport coefficients in stellarators. By means of a convergence study and benchmarks with other codes, it is shown that MONKES is accurate and efficient. The combination of spectral discretization in spatial and velocity coordinates with block sparsity allows MONKES to compute monoenergetic coefficients at low collisionality, in a single core, in approximately one minute. MONKES is sufficiently fast to be integrated into stellarator optimization codes for direct optimization of the bootstrap current and to be included in predictive transport suites. The code and data from this paper are available at https://github.com/JavierEscoto/MONKES/ △ Less

Submitted 20 December, 2023; v1 submitted 19 December, 2023; originally announced December 2023.

Coil Optimization for Quasi-helically Symmetric Stellarator Configurations

Authors: Alexander Vyacheslav Wiedman, Stefan Buller, Matt Landreman

Abstract: Filament-based coil optimizations are performed for several quasihelical stellarator configurations, notably the one from [M. Landreman, E. Paul, PRL 128, 035001, 2022], demonstrating that precise quasihelical symmetry can be achieved with realistic coils. Several constraints are placed on the shape and spacing of the coils, such as low curvature and sufficient plasma-coil distance for neutron shi… ▽ More Filament-based coil optimizations are performed for several quasihelical stellarator configurations, notably the one from [M. Landreman, E. Paul, PRL 128, 035001, 2022], demonstrating that precise quasihelical symmetry can be achieved with realistic coils. Several constraints are placed on the shape and spacing of the coils, such as low curvature and sufficient plasma-coil distance for neutron shielding. The coils resulting from this optimization have a maximum curvature 0.8 times that of the coils of the Helically Symmetric eXperiment (HSX) and a mean squared curvature 0.4 times that of the HSX coils when scaled to the same plasma minor radius. When scaled up to reactor size and magnetic field strength, no fast particle losses were found in the free-boundary configuration when simulating 5000 alpha particles launched at 3.5 MeV on the flux surface with normalized toroidal flux of s=0.5. An analysis of the tolerance of the coils to manufacturing errors is performed using a Gaussian process model, and the coils are found to maintain low particle losses for smooth, large-scale errors up to amplitudes of about 0.15 m. Another coil optimization is performed for the Landreman-Paul configuration with the additional constraint that the coils are purely planar. Visual inspection of the Poincaré plot of the resulting magnetic field-lines reveal that the planar modular coils alone do a poor job of reproducing the target equilibrium. Additional non-planar coil optimizations are performed for the quasihelical configuration with 5% volume-averaged plasma beta from [M. Landreman, S. Buller, M. Drevlak, PoP 29, 082501, 2022], and a similar configuration also optimized to satisfy the Mercier criterion. The finite beta configurations had larger fast-particle losses, with the free-boundary Mercier-optimized configuration performing the worst, losing about 5.5% of alpha particles launched at s=0.5. △ Less

Submitted 27 November, 2023; originally announced November 2023.

Comments: 20 pages, 17 figures

Grad-Shafranov equilibria via data-free physics informed neural networks

Authors: Byoungchan Jang, Alan A. Kaptanoglu, Rahul Gaur, Shaowu Pan, Matt Landreman, William Dorland

Abstract: A large number of magnetohydrodynamic (MHD) equilibrium calculations are often required for uncertainty quantification, optimization, and real-time diagnostic information, making MHD equilibrium codes vital to the field of plasma physics. In this paper, we explore a method for solving the Grad-Shafranov equation by using Physics-Informed Neural Networks (PINNs). For PINNs, we optimize neural netwo… ▽ More A large number of magnetohydrodynamic (MHD) equilibrium calculations are often required for uncertainty quantification, optimization, and real-time diagnostic information, making MHD equilibrium codes vital to the field of plasma physics. In this paper, we explore a method for solving the Grad-Shafranov equation by using Physics-Informed Neural Networks (PINNs). For PINNs, we optimize neural networks by directly minimizing the residual of the PDE as a loss function. We show that PINNs can accurately and effectively solve the Grad-Shafranov equation with several different boundary conditions. We also explore the parameter space by varying the size of the model, the learning rate, and boundary conditions to map various trade-offs such as between reconstruction error and computational speed. Additionally, we introduce a parameterized PINN framework, expanding the input space to include variables such as pressure, aspect ratio, elongation, and triangularity in order to handle a broader range of plasma scenarios within a single network. Parametrized PINNs could be used in future work to solve inverse problems such as shape optimization. △ Less

Submitted 22 November, 2023; originally announced November 2023.

Optimization of Nonlinear Turbulence in Stellarators

Authors: Patrick Kim, Stefan Buller, Rory Conlin, William Dorland, Daniel W. Dudt, Rahul Gaur, Rogerio Jorge, Egemen Kolemen, Matt Landreman, Noah R. Mandell, Dario Panici

Abstract: We present new stellarator equilibria that have been optimized for reduced turbulent transport using nonlinear gyrokinetic simulations within the optimization loop. The optimization routine involves coupling the pseudo-spectral GPU-native gyrokinetic code GX with the stellarator equilibrium and optimization code DESC. Since using GX allows for fast nonlinear simulations, we directly optimize for r… ▽ More We present new stellarator equilibria that have been optimized for reduced turbulent transport using nonlinear gyrokinetic simulations within the optimization loop. The optimization routine involves coupling the pseudo-spectral GPU-native gyrokinetic code GX with the stellarator equilibrium and optimization code DESC. Since using GX allows for fast nonlinear simulations, we directly optimize for reduced nonlinear heat fluxes. To handle the noisy heat flux traces returned by these simulations, we employ the simultaneous perturbation stochastic approximation (SPSA) method that only uses two objective function evaluations for a simple estimate of the gradient. We show several examples that optimize for both reduced heat fluxes and good quasisymmetry as a proxy for low neoclassical transport. Finally, we run full transport simulations using the T3D stellarator transport code to evaluate the changes in the macroscopic profiles. △ Less

Submitted 22 February, 2024; v1 submitted 28 October, 2023; originally announced October 2023.

Efficient calculation of self magnetic field, self-force, and self-inductance for electromagnetic coils. II. Rectangular cross-section

Authors: Matt Landreman, Siena Hurwitz, Thomas M Antonsen Jr

Abstract: For designing high-field electromagnets, the Lorentz force on coils must be computed to ensure a support structure is feasible, and the inductance should be computed to evaluate the stored energy. Also, the magnetic field and its variation inside the conductor is of interest for computing stress and strain, and due to superconducting quench limits. For these force, inductance, energy, and internal… ▽ More For designing high-field electromagnets, the Lorentz force on coils must be computed to ensure a support structure is feasible, and the inductance should be computed to evaluate the stored energy. Also, the magnetic field and its variation inside the conductor is of interest for computing stress and strain, and due to superconducting quench limits. For these force, inductance, energy, and internal field calculations, the coils cannot be naively approximated as infinitesimally thin filaments due to divergences when the source and evaluation points coincide, so more computationally demanding calculations are usually required, resolving the finite cross-section of the conductors. Here, we present a new alternative method that enables the internal magnetic field vector, self-force, and self-inductance to be computed rapidly and accurately within a 1D filament model. The method is applicable to coils for which the curve center-line can have general noncircular shape, as long as the conductor width is small compared to the radius of curvature. This paper extends a previous calculation for circular-cross-section conductors [Hurwitz et al, arXiv:2310.09313 (2023)] to consider the case of rectangular cross-section. The reduced model is derived by rigorous analysis of the singularity, regularizing the filament integrals such that they match the true high-dimensional integrals at high coil aspect ratio. The new filament model exactly recovers analytic results for a circular coil, and is shown to accurately reproduce full finite-cross-section calculations for a non-planar coil of a stellarator magnetic fusion device. Due to the efficiency of the model here, it is well suited for use inside design optimization. △ Less

Submitted 18 October, 2023; originally announced October 2023.

Efficient calculation of the self magnetic field, self-force, and self-inductance for electromagnetic coils

Authors: Siena Hurwitz, Matt Landreman, Thomas M. Antonsen

Abstract: The design of electromagnetic coils may require evaluation of several quantities that are challenging to compute numerically. These quantities include Lorentz forces, which may be a limiting factor due to stresses; the internal magnetic field, which is relevant for determining stress as well as a superconducting coil's proximity to its quench limit; and the inductance, which determines stored magn… ▽ More The design of electromagnetic coils may require evaluation of several quantities that are challenging to compute numerically. These quantities include Lorentz forces, which may be a limiting factor due to stresses; the internal magnetic field, which is relevant for determining stress as well as a superconducting coil's proximity to its quench limit; and the inductance, which determines stored magnetic energy and dynamics. When computing the effect on one coil due to the current in another, these quantities can often be approximated quickly by treating the coils as infinitesimally thin. When computing the effect on a coil due to its own current (e.g., self-force or self-inductance), evaluation is difficult due to the presence of a singularity; coils cannot be treated as infinitesimally thin as each quantity diverges at zero conductor width. Here, we present novel and well-behaved methods for evaluating these quantities using non-singular integral formulae of reduced dimensions. These formulae are determined rigorously by dividing the domain of integration of the magnetic vector potential into two regions, exploiting appropriate approximations in each region, and expanding in high aspect ratio. Our formulae show good agreement to full finite-thickness calculations even at low aspect ratio, both analytically for a torus and numerically for a non-planar coil of a stellarator fusion device, the Helically Symmetric eXperiment (HSX). Because the integrands of these formulae develop fine structure as the minor radius becomes infinitely thin, we also develop a method of evaluating the self-force and self-inductance with even greater efficiency by integrating this sharp feature analytically. We demonstrate with this method that the self-force can be accurately computed for the HSX coil with as few as 12 grid points. △ Less

Submitted 12 October, 2023; originally announced October 2023.

The Magnetic Gradient Scale Length Explains Why Certain Plasmas Require Close External Magnetic Coils

Authors: John Kappel, Matt Landreman, Dhairya Malhotra

Abstract: The separation between the last closed flux surface of a plasma and the external coils that magnetically confine it is a limiting factor in the construction of fusion-capable plasma devices. This plasma-coil separation must be large enough so that components such as a breeding blanket and neutron shielding can fit between the plasma and the coils. Plasma-coil separation affects reactor size, engin… ▽ More The separation between the last closed flux surface of a plasma and the external coils that magnetically confine it is a limiting factor in the construction of fusion-capable plasma devices. This plasma-coil separation must be large enough so that components such as a breeding blanket and neutron shielding can fit between the plasma and the coils. Plasma-coil separation affects reactor size, engineering complexity, and particle loss due to field ripple. For some plasmas it can be difficult to produce the desired flux surface shaping with distant coils, and for other plasmas it is infeasible altogether. Here, we seek to understand the underlying physics that limits plasma-coil separation and explain why some configurations require close external coils. In this paper, we explore the hypothesis that the limiting plasma-coil separation is set by the shortest scale length of the magnetic field as expressed by the $\nabla$B tensor. We tested this hypothesis on a database of > 40 stellarator and tokamak configurations. Within this database, the coil-to-plasma distance compared to the minor radius varies by over an order of magnitude. The magnetic scale length is well correlated to the coil-to-plasma distance of actual coil designs generated using the REGCOIL method [Landreman, Nucl. Fusion 57, 046003 (2017)]. Additionally, this correlation reveals a general trend that larger plasma-coil separation is possible with a small number of field periods. △ Less

Submitted 20 September, 2023; originally announced September 2023.

Comments: 16 pages, 11 figures; fixed typos, added supplemental material Zenodo DOI

Topology optimization for inverse magnetostatics as sparse regression: application to electromagnetic coils for stellarators

Authors: Alan A. Kaptanoglu, Gabriel P. Langlois, Matt Landreman

Abstract: Topology optimization, a technique to determine where material should be placed within a predefined volume in order to minimize a physical objective, is used across a wide range of scientific fields and applications. A general application for topology optimization is inverse magnetostatics; a desired magnetic field is prescribed, and a distribution of steady currents is computed to produce that ta… ▽ More Topology optimization, a technique to determine where material should be placed within a predefined volume in order to minimize a physical objective, is used across a wide range of scientific fields and applications. A general application for topology optimization is inverse magnetostatics; a desired magnetic field is prescribed, and a distribution of steady currents is computed to produce that target field. In the present work, electromagnetic coils are designed by magnetostatic topology optimization, using volume elements (voxels) of electric current, constrained so the current is divergence-free. Compared to standard electromagnet shape optimization, our method has the advantage that the nonlinearity in the Biot-Savart law with respect to position is avoided, enabling convex cost functions and a useful reformulation of topology optimization as sparse regression. To demonstrate, we consider the application of designing electromagnetic coils for a class of plasma experiments known as stellarators. We produce topologically-exotic coils for several new stellarator designs and show that these solutions can be interpolated into a filamentary representation and then further optimized. △ Less

Submitted 21 June, 2023; originally announced June 2023.

Understanding Trade-offs in Stellarator Design with Multi-objective Optimization

Authors: David Bindel, Matt Landreman, Misha Padidar

Abstract: In designing stellarators, any design decision ultimately comes with a trade-off. Improvements in particle confinement, for instance, may increase the burden on engineers to build more complex coils, and the tightening of financial constraints may simplify the design and worsen some aspects of transport. Understanding trade-offs in stellarator designs is critical in designing high performance devi… ▽ More In designing stellarators, any design decision ultimately comes with a trade-off. Improvements in particle confinement, for instance, may increase the burden on engineers to build more complex coils, and the tightening of financial constraints may simplify the design and worsen some aspects of transport. Understanding trade-offs in stellarator designs is critical in designing high performance devices that satisfy the multitude of physical, engineering, and financial criteria. In this study we show how multi-objective optimization (MOO) can be used to investigate trade-offs and develop insight into the role of design parameters. We discuss the basics of MOO, as well as practical solution methods for solving MOO problems. We apply these methods to bring insight into the selection of two common design parameters: the aspect ratio of an ideal magnetohydrodynamic equilibrium, and the total length of the electromagnetic coils. △ Less

Submitted 17 April, 2023; originally announced April 2023.

Direct Optimization of Fast-Ion Confinement in Stellarators

Authors: David Bindel, Matt Landreman, Misha Padidar

Abstract: Confining energetic ions such as alpha particles is a prime concern in the design of stellarators. However, directly measuring alpha confinement through numerical simulation of guiding-center trajectories has been considered to be too computationally expensive and noisy to include in the design loop, and instead has been most often used only as a tool to assess stellarator designs post hoc. In its… ▽ More Confining energetic ions such as alpha particles is a prime concern in the design of stellarators. However, directly measuring alpha confinement through numerical simulation of guiding-center trajectories has been considered to be too computationally expensive and noisy to include in the design loop, and instead has been most often used only as a tool to assess stellarator designs post hoc. In its place, proxy metrics, simplified measures of confinement, have often been used to design configurations because they are computationally more tractable and have been shown to be effective. Despite the success of proxies, it is unclear what is being sacrificed by using them to design the device rather than relying on direct trajectory calculations. In this study, we optimize stellarator designs for improved alpha particle confinement without the use of proxy metrics. In particular, we numerically optimize an objective function that measures alpha particle losses by simulating alpha particle trajectories. While this method is computationally expensive, we find that it can be used successfully to generate configurations with low losses. △ Less

Submitted 22 February, 2023; originally announced February 2023.

Single-Stage Stellarator Optimization: Combining Coils with Fixed Boundary Equilibria

Authors: R. Jorge, A. Goodman, M. Landreman, J. Rodrigues, F. Wechsung

Abstract: We introduce a novel approach for the simultaneous optimization of plasma physics and coil engineering objectives using fixed-boundary equilibria that is computationally efficient and applicable to a broad range of vacuum and finite plasma pressure scenarios. Our approach treats the plasma boundary and coil shapes as independently optimized variables, penalizing the mismatch between the two using… ▽ More We introduce a novel approach for the simultaneous optimization of plasma physics and coil engineering objectives using fixed-boundary equilibria that is computationally efficient and applicable to a broad range of vacuum and finite plasma pressure scenarios. Our approach treats the plasma boundary and coil shapes as independently optimized variables, penalizing the mismatch between the two using a quadratic flux term in the objective function. Four use cases are presented to demonstrate the effectiveness of the approach, including simple and complex stellarator geometries. As shown here, this method outperforms previous 2-stage approaches, achieving smaller plasma objective function values when coils are taken into account. △ Less

Submitted 25 April, 2023; v1 submitted 21 February, 2023; originally announced February 2023.

Comments: 23 pages, 7 figures

An adjoint-based method for optimizing MHD equilibria against the infinite-n, ideal ballooning mode

Authors: Rahul Gaur, Stefan Buller, Maximilian E. Ruth, Matt Landreman, Ian G. Abel, William D. Dorland

Abstract: We demonstrate a fast adjoint-based method to optimize tokamak and stellarator equilibria against a pressure-driven instability known as the infinite-$n$ ideal ballooning mode. We present three finite-$β$ (the ratio of thermal to magnetic pressure) equilibria: one tokamak equilibrium and two stellarator equilibria that are unstable against the ballooning mode. Using the self-adjoint property of id… ▽ More We demonstrate a fast adjoint-based method to optimize tokamak and stellarator equilibria against a pressure-driven instability known as the infinite-$n$ ideal ballooning mode. We present three finite-$β$ (the ratio of thermal to magnetic pressure) equilibria: one tokamak equilibrium and two stellarator equilibria that are unstable against the ballooning mode. Using the self-adjoint property of ideal MHD, we construct a technique to rapidly calculate the change in the growth rate, a measure of ideal ballooning instability. Using the~\texttt{SIMSOPT} framework, we then implement our fast adjoint gradient-based optimizer to minimize the growth rate and find stable equilibria for each of the three initially unstable equilibria. △ Less

Submitted 15 February, 2023; originally announced February 2023.

Comments: 24 pages, 8 tables, 9 figures

Direct Microstability Optimization of Stellarator Devices

Authors: R. Jorge, W. Dorland, P. Kim, M. Landreman, N. R. Mandell, G. Merlo, T. Qian

Abstract: Turbulent transport is regarded as one of the key issues in magnetic confinement nuclear fusion, both for tokamaks in stellarators. In this letter, we show that a significant decrease in the turbulent heat flux can be obtained in an efficient manner by coupling stellarator optimization with linear gyrokinetic simulations. This is accomplished by computing the quasi-linear heat flux at each step of… ▽ More Turbulent transport is regarded as one of the key issues in magnetic confinement nuclear fusion, both for tokamaks in stellarators. In this letter, we show that a significant decrease in the turbulent heat flux can be obtained in an efficient manner by coupling stellarator optimization with linear gyrokinetic simulations. This is accomplished by computing the quasi-linear heat flux at each step of the optimization process, as well as the deviation from quasisymmetry, and minimizing their sum, leading to a balance between neoclassical and turbulent transport. △ Less

Submitted 23 January, 2023; originally announced January 2023.

Comments: 6 pages, 3 figures

Neoclassical transport in strong gradient regions of large aspect ratio tokamaks

Authors: Silvia Trinczek, Felix I. Parra, Peter J. Catto, Iván Calvo, Matt Landreman

Abstract: We present a new neoclassical transport model for large aspect ratio tokamaks where the gradient scale lengths are of the size of the poloidal gyroradius. Previous work on neoclassical transport across transport barriers assumed large density and potential gradients but a small temperature gradient, or neglected the gradient of the mean parallel flow. Using large aspect ratio and low collisionalit… ▽ More We present a new neoclassical transport model for large aspect ratio tokamaks where the gradient scale lengths are of the size of the poloidal gyroradius. Previous work on neoclassical transport across transport barriers assumed large density and potential gradients but a small temperature gradient, or neglected the gradient of the mean parallel flow. Using large aspect ratio and low collisionality expansions, we relax these restrictive assumptions. We define a new set of variables based on conserved quantities, which simplifies the drift kinetic equation whilst keeping strong gradients, and derive equations describing the transport of particles, parallel momentum and energy by ions in the banana regime. The poloidally varying parts of density and electric potential are included. Studying contributions from both passing and trapped particles, we show that the resulting transport is dominated by trapped particles. We find that a non-zero neoclassical particle flux requires parallel momentum input which could be provided through interaction with turbulence or impurities. We derive upper and lower bounds for the energy flux across a transport barrier in both temperature and density and present example profiles and fluxes. △ Less

Submitted 19 April, 2023; v1 submitted 17 January, 2023; originally announced January 2023.

Comments: 52 pages, 13 figures, submitted to Journal of Plasma Physics

Constructing precisely quasi-isodynamic magnetic fields

Authors: Alan Goodman, Katia Camacho Mata, Sophia A Henneberg, Rogerio Jorge, Matt Landreman, Gabriel Plunk, Hakan Smith, Ralf Mackenbach, Per Helander

Abstract: We present a novel method for numerically finding quasi-isodynamic stellarator magnetic fields with excellent fast-particle confinement and extremely small neoclassical transport. The method works particularly well in configurations with only one field period. We examine the properties of these newfound quasi-isodynamic configurations, including their bootstrap currents, particle confinement, and… ▽ More We present a novel method for numerically finding quasi-isodynamic stellarator magnetic fields with excellent fast-particle confinement and extremely small neoclassical transport. The method works particularly well in configurations with only one field period. We examine the properties of these newfound quasi-isodynamic configurations, including their bootstrap currents, particle confinement, and available energy for trapped-electron driven turbulence, as well as the degree to which they change when a finite pressure profile is added. We finally discuss the differences between the magnetic axes of the optimized solutions and their respective initial conditions, and conclude with the prospects for future quasi-isodynamic optimization. △ Less

Submitted 17 November, 2022; originally announced November 2022.

Comments: 25 pages, 10 figures

Direct stellarator coil optimization for nested magnetic surfaces with precise quasi-symmetry

Authors: Andrew Giuliani, Florian Wechsung, Antoine Cerfon, Matt Landreman, Georg Stadler

Abstract: We present a robust optimization algorithm for the design of electromagnetic coils that generate vacuum magnetic fields with nested flux surfaces and precise quasi-symmetry. The method is based on a bilevel optimization problem, where the outer coil optimization is constrained by a set of inner least-squares optimization problems whose solutions describe magnetic surfaces. The outer optimization o… ▽ More We present a robust optimization algorithm for the design of electromagnetic coils that generate vacuum magnetic fields with nested flux surfaces and precise quasi-symmetry. The method is based on a bilevel optimization problem, where the outer coil optimization is constrained by a set of inner least-squares optimization problems whose solutions describe magnetic surfaces. The outer optimization objective targets coils that generate a field with nested magnetic surfaces and good quasi-symmetry. The inner optimization problems identify magnetic surfaces when they exist, and approximate surfaces in the presence of magnetic islands or chaos. We show that this formulation can be used to heal islands and chaos, thus producing coils that result in magnetic fields with precise quasi-symmetry. We show that the method can be initialized with coils from the traditional two stage coil design process, as well as coils from a near axis expansion optimization. We present a numerical example where island chains are healed and quasi-symmetry is optimized up to surfaces with aspect ratio 6. Another numerical example illustrates that the aspect ratio of nested flux surfaces with optimized quasi-symmetry can be decreased from 6 to approximately 4. The last example shows that our approach is robust and a cold-start using coils from a near-axis expansion optimization. △ Less

Submitted 13 March, 2023; v1 submitted 6 October, 2022; originally announced October 2022.

Mapping the space of quasisymmetric stellarators using optimized near-axis expansion

Authors: Matt Landreman

Abstract: A method is demonstrated to rapidly calculate the shapes and properties of quasi-axisymmetric and quasi-helically symmetric stellarators. In this approach, optimization is applied to the equations of magnetohydrodynamic equilibrium and quasisymmetry, expanded in the small distance from the magnetic axis, as formulated by Garren and Boozer [Phys. Fluids B, 3, 2805 (1991)]. Due to the reduction of t… ▽ More A method is demonstrated to rapidly calculate the shapes and properties of quasi-axisymmetric and quasi-helically symmetric stellarators. In this approach, optimization is applied to the equations of magnetohydrodynamic equilibrium and quasisymmetry, expanded in the small distance from the magnetic axis, as formulated by Garren and Boozer [Phys. Fluids B, 3, 2805 (1991)]. Due to the reduction of the equations by the expansion, the computational cost is significantly reduced, to times on the order of 1 cpu-second, enabling wide and high-resolution scans over parameter space. In contrast to traditional stellarator optimization, here the cost function serves to maximize the volume in which the expansion is accurate. A key term in the cost function is $|| \nabla\mathbf B ||$, to maximize scale lengths in the field. Using this method, a database of $5\times 10^5$ optimized configurations is calculated and presented. Quasisymmetric configurations are observed to exist in continuous bands, varying in the ratio of the magnetic axis length to average major radius. Several qualitatively new types of configuration are found, including quasi-helically symmetric fields in which the number of field periods is two or more than six. △ Less

Submitted 8 December, 2022; v1 submitted 23 September, 2022; originally announced September 2022.

Comments: Version accepted in the Journal of Plasma Physics. Minor revisions compared to v1

Greedy permanent magnet optimization

Authors: Alan A. Kaptanoglu, Rory Conlin, Matt Landreman

Abstract: A number of scientific fields rely on placing permanent magnets in order to produce a desired magnetic field. We have shown in recent work that the placement process can be formulated as sparse regression. However, binary, grid-aligned solutions are desired for realistic engineering designs. We now show that the binary permanent magnet problem can be formulated as a quadratic program with quadrati… ▽ More A number of scientific fields rely on placing permanent magnets in order to produce a desired magnetic field. We have shown in recent work that the placement process can be formulated as sparse regression. However, binary, grid-aligned solutions are desired for realistic engineering designs. We now show that the binary permanent magnet problem can be formulated as a quadratic program with quadratic equality constraints (QPQC), the binary, grid-aligned problem is equivalent to the quadratic knapsack problem with multiple knapsack constraints (MdQKP), and the single-orientation-only problem is equivalent to the unconstrained quadratic binary problem (BQP). We then provide a set of simple greedy algorithms for solving variants of permanent magnet optimization, and demonstrate their capabilities by designing magnets for stellarator plasmas. The algorithms can a-priori produce sparse, grid-aligned, binary solutions. Despite its simple design and greedy nature, we provide an algorithm that outperforms the state-of-the-art algorithms while being substantially faster, more flexible, and easier-to-use. △ Less

Submitted 22 August, 2022; originally announced August 2022.

Energetic particle loss mechanisms in reactor-scale equilibria close to quasisymmetry

Authors: E. J. Paul, A. Bhattacharjee, M. Landreman, D. Alex, J. L. Velasco, R. Nies

Abstract: Collisionless physics primarily determines the transport of fusion-born alpha particles in 3D equilibria. Several transport mechanisms have been implicated in stellarator configurations, including stochastic diffusion due to class transitions, ripple trapping, and banana drift-convective orbits. Given the guiding center dynamics in a set of six quasihelical and quasiaxisymmetric equilibria, we per… ▽ More Collisionless physics primarily determines the transport of fusion-born alpha particles in 3D equilibria. Several transport mechanisms have been implicated in stellarator configurations, including stochastic diffusion due to class transitions, ripple trapping, and banana drift-convective orbits. Given the guiding center dynamics in a set of six quasihelical and quasiaxisymmetric equilibria, we perform a classification of trapping states and transport mechanisms. In addition to banana drift convection and ripple transport, we observe substantial non-conservation of the parallel adiabatic invariant which can cause losses through diffusive banana tip motion. Furthermore, many lost trajectories undergo transitions between trapping classes on longer time scales, either with periodic or irregular behavior. We discuss possible optimization strategies for each of the relevant transport mechanisms. We perform a comparison between fast ion losses and metrics for the prevalence of mechanisms such as banana-drift convection [1], transitioning orbits, and wide orbit widths. Quasihelical configurations are found to have natural protection against ripple-trapping and diffusive banana tip motion leading to a reduction in prompt losses. △ Less

Submitted 10 October, 2022; v1 submitted 3 August, 2022; originally announced August 2022.

Stellarator coil optimization supporting multiple magnetic configurations

Authors: Brandon F Lee, Elizabeth J Paul, Georg Stadler, Matt Landreman

Abstract: We present a technique that can be used to design stellarators with a high degree of experimental flexibility. For our purposes, flexibility is defined by the range of values the rotational transform can take on the magnetic axis of the vacuum field while maintaining satisfactory quasisymmetry. We show that accounting for configuration flexibility during the modular coil design improves flexibilit… ▽ More We present a technique that can be used to design stellarators with a high degree of experimental flexibility. For our purposes, flexibility is defined by the range of values the rotational transform can take on the magnetic axis of the vacuum field while maintaining satisfactory quasisymmetry. We show that accounting for configuration flexibility during the modular coil design improves flexibility beyond that attained by previous methods. Careful placement of planar control coils and the incorporation of an integrability objective enhance the quasisymmetry and nested flux surface volume of each configuration. We show that it is possible to achieve flexibility, quasisymmetry, and nested flux surface volume to reasonable degrees with a relatively simple coil set through an NCSX-like example. This example coil design is optimized to achieve three rotational transform targets and nested flux surface volumes in each magnetic configuration larger than the NCSX design plasma volume. Our work suggests that there is a tradeoff between flexibility, quasisymmetry, and volume of nested flux surfaces. △ Less

Submitted 4 November, 2022; v1 submitted 1 August, 2022; originally announced August 2022.

Comments: Minor phrasing change for clarity

Permanent magnet optimization for stellarators as sparse regression

Authors: Alan A. Kaptanoglu, Tony Qian, Florian Wechsung, Matt Landreman

Abstract: A common scientific inverse problem is the placement of magnets that produce a desired magnetic field inside a prescribed volume. This is a key component of stellarator design, and recently permanent magnets have been proposed as a potentially useful tool for magnetic field shaping. Here, we take a closer look at possible objective functions for permanent magnet optimization, reformulate the probl… ▽ More A common scientific inverse problem is the placement of magnets that produce a desired magnetic field inside a prescribed volume. This is a key component of stellarator design, and recently permanent magnets have been proposed as a potentially useful tool for magnetic field shaping. Here, we take a closer look at possible objective functions for permanent magnet optimization, reformulate the problem as sparse regression, and propose an algorithm that can efficiently solve many convex and nonconvex variants. The algorithm generates sparse solutions that are independent of the initial guess, explicitly enforces maximum strengths for the permanent magnets, and accurately produces the desired magnetic field. The algorithm is flexible, and our implementation is open-source and computationally fast. We conclude with two new permanent magnet configurations for the NCSX and MUSE stellarators. Our methodology can be additionally applied for effectively solving permanent magnet optimizations in other scientific fields, as well as for solving quite general high-dimensional, constrained, sparse regression problems, even if a binary solution is required. △ Less

Submitted 18 July, 2022; originally announced July 2022.

A single-field-period quasi-isodynamic stellarator

Authors: R. Jorge, G. G. Plunk, M. Drevlak, M. Landreman, J. -F. Lobsien, K. Camacho Mata, P. Helander

Abstract: A single-field-period quasi-isodynamic stellarator configuration is presented. This configuration, which resembles a twisted strip, is obtained by the method of direct construction, that is, it is found via an expansion in the distance from the magnetic axis. Its discovery, however, relied on an additional step involving numerical optimization, performed within the space of near-axis configuration… ▽ More A single-field-period quasi-isodynamic stellarator configuration is presented. This configuration, which resembles a twisted strip, is obtained by the method of direct construction, that is, it is found via an expansion in the distance from the magnetic axis. Its discovery, however, relied on an additional step involving numerical optimization, performed within the space of near-axis configurations defined by a set of adjustable magnetic-field parameters. This optimization, completed in 30 seconds on a single cpu core using the SIMSOPT code, yields a solution with excellent confinement, as measured by the conventional figure of merit for neoclassical transport, effective ripple, at a modest aspect ratio of eight. The optimization parameters that led to this configuration are described, its confinement properties are assessed, and a set of magnetic-field coils is found. The resulting transport at low collisionality is much smaller than that of W7-X, and the device needs significantly fewer coils thanks to the reduced number of field periods. △ Less

Submitted 1 September, 2022; v1 submitted 11 May, 2022; originally announced May 2022.

Comments: 13 pages, 9 figures

Optimization of quasisymmetric stellarators with self-consistent bootstrap current and energetic particle confinement

Authors: Matt Landreman, Stefan Buller, Michael Drevlak

Abstract: Quasisymmetry can greatly improve the confinement of energetic particles and thermal plasma in a stellarator. The magnetic field of a quasisymmetric stellarator at high plasma pressure is significantly affected by the bootstrap current, but the computational cost of accurate stellarator bootstrap calculations has precluded use inside optimization. Here, a new efficient method is demonstrated for o… ▽ More Quasisymmetry can greatly improve the confinement of energetic particles and thermal plasma in a stellarator. The magnetic field of a quasisymmetric stellarator at high plasma pressure is significantly affected by the bootstrap current, but the computational cost of accurate stellarator bootstrap calculations has precluded use inside optimization. Here, a new efficient method is demonstrated for optimization of quasisymmetric stellarator configurations such that the bootstrap current profile is consistent with the geometry. The approach is based on the fact that all neoclassical phenomena in quasisymmetry are isomorphic to those in axisymmetry. Therefore accurate formulae for the bootstrap current in tokamaks, which can be evaluated rapidly, can be applied also in stellarators. The deviation between this predicted parallel current and the actual parallel current in the magnetohydrodynamic equilibrium is penalized in the objective function, and the current profile of the equilibrium is included in the parameter space. Quasisymmetric configurations with significant pressure are thereby obtained with self-consistent bootstrap current and excellent confinement. In a comparison of fusion-produced alpha particle confinement across many stellarators, the new configurations have significantly lower alpha energy losses than many previous designs. △ Less

Submitted 26 July, 2022; v1 submitted 5 May, 2022; originally announced May 2022.

Comments: Version to appear in Physics of Plasmas

Stochastic and a posteriori optimization to mitigate coil manufacturing errors in stellarator design

Authors: Florian Wechsung, Andrew Giuliani, Matt Landreman, Antoine Cerfon, Georg Stadler

Abstract: It was recently shown in [Wechsung et. al., Proc. Natl. Acad. Sci. USA, 2022] that there exist electromagnetic coils that generate magnetic fields which are excellent approximations to quasi-symmetric fields and have very good particle confinement properties. Using a Gaussian process based model for coil perturbations, we investigate the impact of manufacturing errors on the performance of these c… ▽ More It was recently shown in [Wechsung et. al., Proc. Natl. Acad. Sci. USA, 2022] that there exist electromagnetic coils that generate magnetic fields which are excellent approximations to quasi-symmetric fields and have very good particle confinement properties. Using a Gaussian process based model for coil perturbations, we investigate the impact of manufacturing errors on the performance of these coils. We show that even fairly small errors result in noticeable performance degradation. While stochastic optimization yields minor improvements, it is not able to mitigate these errors significantly. As an alternative to stochastic optimization, we then formulate a new optimization problem for computing optimal adjustments of the coil positions and currents without changing the shapes of the coil. These a-posteriori adjustments are able to reduce the impact of coil errors by an order of magnitude, providing a new perspective for dealing with manufacturing tolerances in stellarator design. △ Less

Submitted 30 July, 2022; v1 submitted 18 March, 2022; originally announced March 2022.

Direct computation of magnetic surfaces in Boozer coordinates and coil optimization for quasi-symmetry

Authors: Andrew Giuliani, Florian Wechsung, Matt Landreman, Georg Stadler, Antoine Cerfon

Abstract: We propose a new method to compute magnetic surfaces that are parametrized in Boozer coordinates for vacuum magnetic fields. We also propose a measure for quasi-symmetry on the computed surfaces and use it to design coils that generate a magnetic field that is quasi-symmetric on those surfaces. The rotational transform of the field and complexity measures for the coils are also controlled in the d… ▽ More We propose a new method to compute magnetic surfaces that are parametrized in Boozer coordinates for vacuum magnetic fields. We also propose a measure for quasi-symmetry on the computed surfaces and use it to design coils that generate a magnetic field that is quasi-symmetric on those surfaces. The rotational transform of the field and complexity measures for the coils are also controlled in the design problem. Using an adjoint approach, we are able to obtain analytic derivatives for this optimization problem, yielding an efficient gradient-based algorithm. Starting from an initial coil set that presents nested magnetic surfaces for a large fraction of the volume, our method converges rapidly to coil systems generating fields with excellent quasi-symmetry and low particle losses. In particular for low complexity coils, we are able to significantly improve the performance compared to coils obtained from the standard two-stage approach, e.g.~reduce losses of fusion-produced alpha particles born at half-radius from $17.7\%$ to $6.6\%$. We also demonstrate 16-coil configurations with alpha loss < $1\%$ and neoclassical transport magnitude $ε_{\mathrm{eff}}^{3/2}$ less than approximately $5\times 10^{-9}.$ △ Less

Submitted 29 April, 2022; v1 submitted 7 March, 2022; originally announced March 2022.

Stellarator optimization for nested magnetic surfaces at finite $β$ and toroidal current

Authors: A. Baillod, J. Loizu, J. P. Graves, M. Landreman

Abstract: Good magnetic surfaces, as opposed to magnetic islands and chaotic field lines, are generally desirable for stellarators. In previous work, M. Landreman et al. [Phys. of Plasmas 28, 092505 (2021)] showed that equilibria computed by the Stepped-Pressure Equilibrium Code (SPEC) [S. P. Hudson et al., Phys. Plasmas 19, 112502 (2012)] could be optimized for good magnetic surfaces in vacuum. In this pap… ▽ More Good magnetic surfaces, as opposed to magnetic islands and chaotic field lines, are generally desirable for stellarators. In previous work, M. Landreman et al. [Phys. of Plasmas 28, 092505 (2021)] showed that equilibria computed by the Stepped-Pressure Equilibrium Code (SPEC) [S. P. Hudson et al., Phys. Plasmas 19, 112502 (2012)] could be optimized for good magnetic surfaces in vacuum. In this paper, we build upon their work to show the first finite-$β$, fixed- and free-boundary optimization of SPEC equilibria for good magnetic surfaces. The objective function is constructed with the Greene's residue of selected rational surfaces and the optimization is driven by the SIMSOPT framework [M. Landreman et al., J. Open Source Software 6, 3525 (2021)]. We show that the size of magnetic islands and the consequent regions occupied by chaotic field lines can be minimized in a classical stellarator geometry by optimizing either the injected toroidal current profile, the shape of a perfectly conducting wall surrounding the plasma (fixed-boundary case), or the coils (free-boundary case), in a reasonable amount of computational time. This work shows that SPEC can be used as an equilibrium code both in a two-step or single-step stellarator optimization loop. △ Less

Submitted 30 November, 2021; originally announced November 2021.

Comments: 6 pages, 4 figures, to be submitted to Physics of Plasmas

Journal ref: Physics of Plasmas 29, 042505 (2022)

Magnetic fields with precise quasisymmetry for plasma confinement

Authors: Matt Landreman, Elizabeth Paul

Abstract: Quasisymmetry is an unusual symmetry that can be present in toroidal magnetic fields, enabling confinement of charged particles and plasma. Here it is shown that both quasi-axisymmetry and quasi-helical symmetry can be achieved to a much higher precision than previously thought over a significant volume, resulting in exceptional confinement. For a 1 Tesla mean field far from axisymmetry (vacuum ro… ▽ More Quasisymmetry is an unusual symmetry that can be present in toroidal magnetic fields, enabling confinement of charged particles and plasma. Here it is shown that both quasi-axisymmetry and quasi-helical symmetry can be achieved to a much higher precision than previously thought over a significant volume, resulting in exceptional confinement. For a 1 Tesla mean field far from axisymmetry (vacuum rotational transform $>$ 0.4), symmetry-breaking mode amplitudes throughout a volume of aspect ratio 6 can be made as small as the typical $\sim 50$ $μ$T geomagnetic field. △ Less

Submitted 23 November, 2021; v1 submitted 8 August, 2021; originally announced August 2021.

Comments: Version to appear in Physical Review Letters, including new Supplemental Material

Stellarator optimization for good magnetic surfaces at the same time as quasisymmetry

Authors: Matt Landreman, Bharat Medasani, Caoxiang Zhu

Abstract: A method is demonstrated to optimize a stellarator's geometry to eliminate magnetic islands and achieve other desired physics properties at the same time. For many physics quantities that have been used in stellarator optimization, including quasisymmetry, neoclassical transport, and magnetohydrodynamic stability, it is convenient to use a magnetic equilibrium representation that assures the exist… ▽ More A method is demonstrated to optimize a stellarator's geometry to eliminate magnetic islands and achieve other desired physics properties at the same time. For many physics quantities that have been used in stellarator optimization, including quasisymmetry, neoclassical transport, and magnetohydrodynamic stability, it is convenient to use a magnetic equilibrium representation that assures the existence of magnetic surfaces. However, this representation hides the possible presence of magnetic islands, which are typically undesirable. To include both surface-based objectives and island widths in a single optimization, two fixed-boundary equilibrium calculations are run at each iteration of the optimization: one that enforces the existence of magnetic surfaces (VMEC [S. P. Hirshman and J. C. Whitson, Phys. Fluids 26, 3553 (1983)]), and one that does not (SPEC [S. R. Hudson, et al, Phys. Plasmas 19, 112502 (2012)]). By penalizing the island residues in the objective function, the two magnetic field representations are brought into agreement during the optimization. An example is presented in which, particularly on the surface where quasisymmetry was targeted, quasisymmetry is achieved more accurately than in previously published examples. △ Less

Submitted 28 June, 2021; originally announced June 2021.

Single-stage gradient-based stellarator coil design: stochastic optimization

Authors: Florian Wechsung, Andrew Giuliani, Matt Landreman, Antoine Cerfon, Georg Stadler

Abstract: We extend the single-stage stellarator coil design approach for quasi-symmetry on axis from [Giuliani et al, 2020] to additionally take into account coil manufacturing errors. By modeling coil errors independently from the coil discretization, we have the flexibility to consider realistic forms of coil errors. The corresponding stochastic optimization problems are formulated using a risk-neutral a… ▽ More We extend the single-stage stellarator coil design approach for quasi-symmetry on axis from [Giuliani et al, 2020] to additionally take into account coil manufacturing errors. By modeling coil errors independently from the coil discretization, we have the flexibility to consider realistic forms of coil errors. The corresponding stochastic optimization problems are formulated using a risk-neutral approach and risk-averse approaches. We present an efficient, gradient-based descent algorithm which relies on analytical derivatives to solve these problems. In a comprehensive numerical study, we compare the coil designs resulting from deterministic and risk-neutral stochastic optimization and find that the risk-neutral formulation results in more robust configurations and reduces the number of local minima of the optimization problem. We also compare deterministic and risk-neutral approaches in terms of quasi-symmetry on and away from the magnetic axis, and in terms of the confinement of particles released close to the axis. Finally, we show that for the optimization problems we consider, a risk-averse objective using the Conditional Value-at-Risk leads to results which are similar to the risk-neutral objective. △ Less

Submitted 22 June, 2021; originally announced June 2021.

A neoclassically optimized compact stellarator with four planar coils

Authors: Guodong Yu, Zhichen Feng, Peiyou Jiang, Neil Pomphrey, Matt Landreman, GuoYong Fu

Abstract: 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 th… ▽ More 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 CNT's. △ Less

Submitted 10 June, 2021; originally announced June 2021.

Energetic particle transport in optimized stellarators

Authors: A Bader, D T Anderson, M Drevlak, B J Faber, C C Hegna, S Henneberg, M Landreman, J C Schmitt, Y Suzuki, A Ware

Abstract: Nine stellarator configurations, three quasiaxisymmetric, three quasihelically symmetric and three non-quasisymmetric are scaled to ARIES-CS size and analyzed for energetic particle content. The best performing configurations with regard to energetic particle confinement also perform the best on the neoclassical Γc metric, which attempts to align contours of the second adiabatic invariant with flu… ▽ More Nine stellarator configurations, three quasiaxisymmetric, three quasihelically symmetric and three non-quasisymmetric are scaled to ARIES-CS size and analyzed for energetic particle content. The best performing configurations with regard to energetic particle confinement also perform the best on the neoclassical Γc metric, which attempts to align contours of the second adiabatic invariant with flux surfaces. Quasisymmetric configurations that simultaneously perform well on Γc and quasisymmetry have the best overall confinement, with collisional losses under 3%, approaching the performance of ITER with ferritic inserts. △ Less

Submitted 1 June, 2021; originally announced June 2021.

Comments: 18 pages, 10 figures, 1 Table

Ion-temperature-gradient stability near the magnetic axis of quasisymmetric stellarators

Authors: R. Jorge, M. Landreman

Abstract: The stability of the ion-temperature gradient mode in quasisymmetric stellarators is assessed. This is performed using a set of analytical estimates together with linear gyrokinetic simulations. The peak growth rates, their corresponding real frequencies and wave-vectors are identified. A comparison is made between a first-order near-axis expansion model and eleven realistic designs obtained using… ▽ More The stability of the ion-temperature gradient mode in quasisymmetric stellarators is assessed. This is performed using a set of analytical estimates together with linear gyrokinetic simulations. The peak growth rates, their corresponding real frequencies and wave-vectors are identified. A comparison is made between a first-order near-axis expansion model and eleven realistic designs obtained using numerical optimization methods. It is found that while the near-axis expansion is able to replicate the growth rates, real frequencies and perpendicular wave-vector at the inner core (both using simplified dispersion relations and first-principle gyrokinetic simulations), it leads to an overestimation of the growth rate at larger radii. An approximate analytic solution of the ITG dispersion relation for the non-resonant limit suggests growth rates could be systematically higher in quasi-axisymmetric (QA) configurations compared to quasi-helically (QH) symmetric ones. However except for very close to the axis, linear gyrokinetic simulations do not show systematic differences between QA and QH configurations. △ Less

Submitted 24 February, 2021; originally announced February 2021.

An adjoint method for determining the sensitivity of island size to magnetic field variations

Authors: Alessandro Geraldini, Matt Landreman, Elizabeth Paul

Abstract: An adjoint method to calculate the gradient of island width in stellarators is presented and applied to a set of magnetic field configurations. The underlying method of calculation of the island width is that of Cary & Hanson (1991) (with a minor modification), and requires that the residue of the island centre be small. Therefore, the gradient of the residue is calculated in addition. Both the is… ▽ More An adjoint method to calculate the gradient of island width in stellarators is presented and applied to a set of magnetic field configurations. The underlying method of calculation of the island width is that of Cary & Hanson (1991) (with a minor modification), and requires that the residue of the island centre be small. Therefore, the gradient of the residue is calculated in addition. Both the island width and the gradient calculations are verified using an analytical magnetic field configuration introduced in Reiman & Greenside (1986). The method is also applied to the calculation of the shape gradient of the width of a magnetic island in an NCSX vacuum configuration with respect to positions on a coil. A gradient-based optimization is applied to a magnetic field configuration studied in Hanson & Cary (1984) to minimize stochasticity by adding perturbations to a pair of helical coils. Although only vacuum magnetic fields and an analytical magnetic field model are considered in this work, the adjoint calculation of the island width gradient could also be applied to a magnetohydrodynamic (MHD) equilibrium if the derivative of the magnetic field with respect to the equilibrium parameters was known. Using the island width gradient calculation presented here, more general gradient-based optimization methods can be applied to design stellarators with small magnetic islands. Moreoever, the sensitivity of the island size may itself be optimized to ensure that coil tolerances with respect to island size are kept as high as possible. △ Less

Submitted 8 February, 2021; originally announced February 2021.

Comments: 41 pages, 8 figures

Gradient-based optimization of 3D MHD equilibria

Authors: Elizabeth Paul, Matt Landreman, Thomas Antonsen Jr

Abstract: Using recently developed adjoint methods for computing the shape derivatives of functions that depend on MHD equilibria (Antonsen et al. 2019; Paul et al. 2020), we present the first example of analytic gradient-based optimization of fixed-boundary stellarator equilibria. We take advantage of gradient information to optimize figures of merit of relevance for stellarator design, including the rotat… ▽ More Using recently developed adjoint methods for computing the shape derivatives of functions that depend on MHD equilibria (Antonsen et al. 2019; Paul et al. 2020), we present the first example of analytic gradient-based optimization of fixed-boundary stellarator equilibria. We take advantage of gradient information to optimize figures of merit of relevance for stellarator design, including the rotational transform, magnetic well, and quasisymmetry near the axis. With the application of the adjoint method, we reduce the number of equilibrium evaluations by the dimension of the optimization space ($\sim 50-500$) in comparison with a finite-difference gradient-based method. We discuss regularization objectives of relevance for fixed-boundary optimization, including a novel method that prevents self-intersection of the plasma boundary. We present several optimized equilibria, including a vacuum field with very low magnetic shear throughout the volume. △ Less

Submitted 1 February, 2021; v1 submitted 17 December, 2020; originally announced December 2020.

Figures of merit for stellarators near the magnetic axis

Authors: Matt Landreman

Abstract: A new paradigm for rapid stellarator configuration design has been recently demonstrated, in which the shapes of quasisymmetric or omnigenous flux surfaces are computed directly using an expansion in small distance from the magnetic axis. To further develop this approach, here we derive several other quantities of interest that can be rapidly computed from this near-axis expansion. First, the… ▽ More A new paradigm for rapid stellarator configuration design has been recently demonstrated, in which the shapes of quasisymmetric or omnigenous flux surfaces are computed directly using an expansion in small distance from the magnetic axis. To further develop this approach, here we derive several other quantities of interest that can be rapidly computed from this near-axis expansion. First, the $\nabla\vec{B}$ and $\nabla\nabla\vec{B}$ tensors are computed, which can be used for direct derivative-based optimization of electromagnetic coil shapes to achieve the desired magnetic configuration. Moreover, if the norm of these tensors is large compared to the field strength for a given magnetic field, the field must have a short length scale, suggesting it may be hard to produce with coils that are suitably far away. Second, we evaluate the minor radius at which the flux surface shapes would become singular, providing a lower bound on the achievable aspect ratio. This bound is also shown to be related to an equilibrium beta limit. Finally, for configurations that are constructed to achieve a desired magnetic field strength to first order in the expansion, we compute the error field that arises due to second order terms. △ Less

Submitted 5 March, 2021; v1 submitted 1 December, 2020; originally announced December 2020.

Journal ref: J. Plasma Phys. 87 (2021) 905870112

Single-stage gradient-based stellarator coil design: Optimization for near-axis quasi-symmetry

Authors: Andrew Giuliani, Florian Wechsung, Antoine Cerfon, Georg Stadler, Matt Landreman

Abstract: We present a new coil design paradigm for magnetic confinement in stellarators. Our approach directly optimizes coil shapes and coil currents to produce a vacuum quasi-symmetric magnetic field with a target rotational transform on the magnetic axis. This approach differs from the traditional two-stage approach in which first a magnetic configuration with desirable physics properties is found, and… ▽ More We present a new coil design paradigm for magnetic confinement in stellarators. Our approach directly optimizes coil shapes and coil currents to produce a vacuum quasi-symmetric magnetic field with a target rotational transform on the magnetic axis. This approach differs from the traditional two-stage approach in which first a magnetic configuration with desirable physics properties is found, and then coils to approximately realize this magnetic configuration are designed. The proposed single-stage approach allows us to find a compromise between confinement and engineering requirements, i.e., find easy-to-build coils with good confinement properties. Using forward and adjoint sensitivities, we derive derivatives of the physical quantities in the objective, which is constrained by a nonlinear periodic differential equation. In two numerical examples, we compare different gradient-based descent algorithms and find that incorporating approximate second-order derivative information through a quasi-Newton method is crucial for convergence. We also explore the optimization landscape in the neighborhood of a minimizer and find many directions in which the objective is mostly flat, indicating ample freedom to find simple and thus easy-to-build coils. △ Less

Submitted 15 March, 2022; v1 submitted 1 October, 2020; originally announced October 2020.

Comments: 28 pages

Calculation of permanent magnet arrangements for stellarators: A linear least-squares method

Authors: Matt Landreman, Caoxiang Zhu

Abstract: A problem arising in several engineering areas is to design magnets outside a volume that produce a desired magnetic field inside it. One instance of this problem is stellarator design, where it has recently been shown that permanent magnets can provide the required shaping of the magnetic field. Here we demonstrate a robust and efficient algorithm REGCOIL_PM to calculate the spatial distribution… ▽ More A problem arising in several engineering areas is to design magnets outside a volume that produce a desired magnetic field inside it. One instance of this problem is stellarator design, where it has recently been shown that permanent magnets can provide the required shaping of the magnetic field. Here we demonstrate a robust and efficient algorithm REGCOIL_PM to calculate the spatial distribution of these permanent magnets. The procedure involves a small number of fixed-point iterations, with a linear least-squares problem solved at each step. The method exploits the Biot-Savart Law's exact linearity in magnetization density and approximate linearity in magnet size, for magnets far from the target region. No constraint is placed on the direction of magnetization, so Halbach solutions are found naturally, and the magnitude of the magnetization can be made uniformly equal to a target value. △ Less

Submitted 14 September, 2020; originally announced September 2020.

The Use of Near-Axis Magnetic Fields for Stellarator Turbulence Simulations

Authors: R. Jorge, M. Landreman

Abstract: The design of turbulence optimized stellarators has so far relied on three-dimensional equilibrium codes such as VMEC in order to find the minimum of a given objective function. In this work, we propose a complimentary approach based on the near-axis expansion to compute the geometry parameters of neoclassicaly optimized stellarators used in turbulence studies. As shown here, the near-axis expansi… ▽ More The design of turbulence optimized stellarators has so far relied on three-dimensional equilibrium codes such as VMEC in order to find the minimum of a given objective function. In this work, we propose a complimentary approach based on the near-axis expansion to compute the geometry parameters of neoclassicaly optimized stellarators used in turbulence studies. As shown here, the near-axis expansion can be a reasonable approximation of the geometric parameters relevant for turbulence and stability simulations of the core of existing optimized stellarator designs. In particular, we examine the geometry coefficients that appear in the gyrokinetic equation, the drift-reduced fluid equations and the ideal ballooning equation. This approach may allow for the development of new stellarator optimization techniques significantly faster than conventional methods. △ Less

Submitted 20 August, 2020; originally announced August 2020.

Comments: 29 pages, 14 figures

Magnetic well and Mercier stability of stellarators near the magnetic axis

Authors: Matt Landreman, Rogerio Jorge

Abstract: We have recently demonstrated that by expanding in small distance from the magnetic axis compared to the major radius, stellarator shapes with low neoclassical transport can be generated efficiently. To extend the utility of this new design approach, here we evaluate measures of magnetohydrodynamic interchange stability within the same expansion. In particular, we evaluate magnetic well, Mercier's… ▽ More We have recently demonstrated that by expanding in small distance from the magnetic axis compared to the major radius, stellarator shapes with low neoclassical transport can be generated efficiently. To extend the utility of this new design approach, here we evaluate measures of magnetohydrodynamic interchange stability within the same expansion. In particular, we evaluate magnetic well, Mercier's criterion, and resistive interchange stability near a magnetic axis of arbitrary shape. In contrast to previous work on interchange stability near the magnetic axis, which used an expansion of the flux coordinates, here we use the inverse expansion in which the flux coordinates are the independent variables. Reduced expressions are presented for the magnetic well and stability criterion in the case of quasisymmetry. The analytic results are shown to agree with calculations from the VMEC equilibrium code. Finally, we show that near the axis, Glasser, Greene, & Johnson's stability criterion for resistive modes approximately coincides with Mercier's ideal condition. △ Less

Submitted 2 September, 2020; v1 submitted 26 June, 2020; originally announced June 2020.

A New Optimized Quasihelically SymmetricStellarator

Authors: A. Bader, B. J. Faber, J. C. Schmitt, D. T. Anderson, M. Drevlak, J. M. Duff, H. Frerichs, C. C. Hegna, T. G. Kruger, M. Landreman, I. J. McKinney, L. Singh, J. M. Schroeder, P. W. Terry, A. S. Ware

Abstract: A new optimized quasihelically symmetric configuration is described that has the desir-able properties of improved energetic particle confinement, reduced turbulent transportby 3D shaping, and non-resonant divertor capabilities. The configuration presented in thispaper is explicitly optimized for quasihelical symmetry, energetic particle confinement,neoclassical confinement, and stability near the… ▽ More A new optimized quasihelically symmetric configuration is described that has the desir-able properties of improved energetic particle confinement, reduced turbulent transportby 3D shaping, and non-resonant divertor capabilities. The configuration presented in thispaper is explicitly optimized for quasihelical symmetry, energetic particle confinement,neoclassical confinement, and stability near the axis. Post optimization, the configurationwas evaluated for its performance with regard to energetic particle transport, idealmagnetohydrodynamic (MHD) stability at various values of plasma pressure, and iontemperature gradient instability induced turbulent transport. The effect of discrete coilson various confinement figures of merit, including energetic particle confinement, aredetermined by generating single-filament coils for the configuration. Preliminary divertoranalysis shows that coils can be created that do not interfere with expansion of thevessel volume near the regions of outgoing heat flux, thus demonstrating the possibilityof operating a non-resonant divertor. △ Less

Submitted 23 April, 2020; originally announced April 2020.

Comments: 13 figures, 3 tables

Construction of Quasisymmetric Stellarators Using a Direct Coordinate Approach

Authors: R. Jorge, W. Sengupta, M. Landreman

Abstract: Optimized stellarator configurations and their analytical properties are obtained using a near-axis expansion approach. Such configurations are associated with good confinement as the guiding center particle trajectories and neoclassical transport are isomorphic to those in a tokamak. This makes them appealing as fusion reactor candidates. Using a direct coordinate approach, where the magnetic fie… ▽ More Optimized stellarator configurations and their analytical properties are obtained using a near-axis expansion approach. Such configurations are associated with good confinement as the guiding center particle trajectories and neoclassical transport are isomorphic to those in a tokamak. This makes them appealing as fusion reactor candidates. Using a direct coordinate approach, where the magnetic field and flux surface functions are found explicitly in terms of the position vector at successive orders in the distance to the axis, the set of ordinary differential equations for first and second order quasisymmetry is derived. Examples of quasi-axisymmetric shapes are constructed using a pseudospectral numerical method. Finally, the direct coordinate approach is used to independently verify two hypotheses commonly associated with quasisymmetric magnetic fields, namely that the number of equations exceeds the number of parameters at third order in the expansion and that the near-axis expansion does not prohibit exact quasisymmetry from being achieved on a single flux surface. △ Less

Submitted 13 March, 2020; originally announced March 2020.

Comments: 28 pages, 5 figures

Impurity temperature screening in stellarators close to quasisymmetry

Authors: Mike F Martin, Matt Landreman

Abstract: Impurity temperature screening is a favorable neoclassical phenomenon involving an outward radial flux of impurity ions from the core of fusion devices. Quasisymmetric magnetic fields are intrinsically ambipolar and give rise to temperature screening for low enough $η^{-1}\equiv d\ln n/d\ln T$. In contrast, neoclassical fluxes in generic stellarators will depend on the radial electric field, which… ▽ More Impurity temperature screening is a favorable neoclassical phenomenon involving an outward radial flux of impurity ions from the core of fusion devices. Quasisymmetric magnetic fields are intrinsically ambipolar and give rise to temperature screening for low enough $η^{-1}\equiv d\ln n/d\ln T$. In contrast, neoclassical fluxes in generic stellarators will depend on the radial electric field, which is predicted to be inward for ion-root plasmas, potentially leading to impurity accumulation. Here we examine the impurity particle flux in a number of approximately quasisymmetric stellarator configurations and parameter regimes while varying the amount of symmetry-breaking in the magnetic field. The amount of symmetry-breaking (from perfect quasisymmetry) required to lose temperature screening is dependent on the particular configuration and flux surface. Neoclassical fluxes have been obtained using the SFINCS drift-kinetic equation solver with the electrostatic potential $Φ=Φ(ψ)$, where $ψ$ is a flux surface label. Results indicate that achieving temperature screening is possible, but unlikely, at reactor-relevant conditions in the core. Thus, the small departures from symmetry in nominally quasisymmetric stellarators are large enough to significantly alter the neoclassical impurity transport. A further look at these fluxes when compared to a gyro-Bohm turbulence estimate suggests that neoclassical fluxes are small in optimized configurations compared to turbulent fluxes. Therefore, although neoclassical impurity accumulation is expected in most situations, the strength of the turbulence may render this irrelevant. △ Less

Submitted 22 February, 2020; originally announced February 2020.

Near-Axis Expansion of Stellarator Equilibrium at Arbitrary Order in the Distance to the Axis

Authors: R. Jorge, W. Sengupta, M. Landreman

Abstract: 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 character of high-aspect-ratio stellarator shapes, such as… ▽ More 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 character of high-aspect-ratio stellarator shapes, such as the dependence of the rotational transform and the plasma beta-limit on geometrical properties of the resulting flux surfaces. The approach developed here is based on an asymptotic expansion on the inverse aspect-ratio of the ideal MHD equation. The analysis is simplified by using an orthogonal coordinate system relative to the Frenet-Serret frame at the magnetic axis. The magnetic field vector, the toroidal magnetic flux, the current density, the field line label and the rotational transform are derived at arbitrary order in the expansion parameter. Moreover, a comparison with a near-axis expansion formalism employing an inverse coordinate method based on Boozer coordinates (the so-called Garren-Boozer construction) is made, where both methods are shown to agree at lowest order. Finally, as a practical example, a numerical solution using a W7-X equilibrium is presented, and a comparison between the lowest order solution and the W7-X magnetic field is performed. △ Less

Submitted 28 January, 2020; v1 submitted 6 November, 2019; originally announced November 2019.

Comments: 32 pages, 6 figures, 1 table

Journal ref: J. Plasma Phys. 86 (2020) 905860106

Adjoint approach to calculating shape gradients for three-dimensional magnetic confinement equilibria. Part II: Applications

Authors: Elizabeth J. Paul, Thomas Antonsen, Jr., Matt Landreman, W. Anthony Cooper

Abstract: The shape gradient is a local sensitivity function that provides the change in a figure of merit associated with a perturbation to the shape of the object. The shape gradient can be used for gradient-based optimization, sensitivity analysis, and tolerance calculations. However, it is generally expensive to compute from finite-difference derivatives for shapes which are described by many parameters… ▽ More The shape gradient is a local sensitivity function that provides the change in a figure of merit associated with a perturbation to the shape of the object. The shape gradient can be used for gradient-based optimization, sensitivity analysis, and tolerance calculations. However, it is generally expensive to compute from finite-difference derivatives for shapes which are described by many parameters, as is the case for stellarator geometry. In an accompanying work (Antonsen et al. 2019), generalized self-adjointness relations are obtained for MHD equilibria. These describe the relation between perturbed equilibria due to changes in the rotational transform or toroidal current profiles, displacements of the plasma boundary, modifications of currents in the vacuum region, or the addition of bulk forces. These are applied to efficiently compute the shape gradient of functions of magnetohydrodynamic (MHD) equilibria with an adjoint approach. In this way, the shape derivative with respect to any perturbation applied to the plasma boundary or coil shapes can be computed with only one additional MHD equilibrium solution. We demonstrate that this approach is applicable for several figures of merit of interest for stellarator configuration optimization: the magnetic well, the magnetic ripple on axis, the departure from quasisymmetry, the effective ripple in the low-collisionality $1/ν$ regime ($ε_{\text{eff}}^{3/2})$ (Nemov et al. 1999), and several finite collisionality neoclassical quantities. Numerical verification is demonstrated for the magnetic well figure of merit with the VMEC code (Hirshman & Whitson 1983) and for the magnetic ripple with modification of the ANIMEC code (Cooper et al. 1992). Comparisons with the direct approach demonstrate that to obtain agreement within several percent, the adjoint approach provides a factor of $\mathcal{O}(10^3)$ in computational savings. △ Less

Submitted 13 December, 2019; v1 submitted 30 October, 2019; originally announced October 2019.

Direct construction of optimized stellarator shapes. III. Omnigenity near the magnetic axis

Authors: G. G. Plunk, M. Landreman, P. Helander

Abstract: 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 provid… ▽ More 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. △ Less

Submitted 19 September, 2019; originally announced September 2019.

Constructing stellarators with quasisymmetry to high order

Authors: Matt Landreman, Wrick Sengupta

Abstract: 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 Boozer's expansion about the magnetic axis (Phys Fluid… ▽ More 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 Boozer's expansion about the magnetic axis (Phys Fluids B 3, 2805 (1991)), and it is several orders of magnitude faster than the conventional optimization approach. The work here extends the method of Landreman, Sengupta and Plunk (J Plasma Phys 85, 905850103 (2019)), 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 inverse aspect ratio. We also demonstrate a strongly nonaxisymmetric configuration (vacuum $ι> 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. △ Less

Submitted 14 October, 2019; v1 submitted 27 August, 2019; originally announced August 2019.

Journal ref: J. Plasma Phys. 85 (2019) 815850601

An adjoint method for neoclassical stellarator optimization

Authors: Elizabeth Paul, Ian Abel, Matt Landreman, William Dorland

Abstract: Stellarators are a promising route to steady-state fusion power. However, to achieve the required confinement, the magnetic geometry must be highly optimized. This optimization requires navigating high-dimensional spaces, often necessitating the use of gradient-based methods. The gradient of the neoclassical fluxes is expensive to compute with classical methods, requiring $O(N)$ flux computations,… ▽ More Stellarators are a promising route to steady-state fusion power. However, to achieve the required confinement, the magnetic geometry must be highly optimized. This optimization requires navigating high-dimensional spaces, often necessitating the use of gradient-based methods. The gradient of the neoclassical fluxes is expensive to compute with classical methods, requiring $O(N)$ flux computations, where $N$ is the number of parameters. To reduce the cost of the gradient computation, we present an adjoint method for computing the derivatives of moments of the neoclassical distribution function for stellarator optimization. The linear adjoint method allows derivatives of quantities which depend on solutions of a linear system, such as moments of the distribution function, to be computed with respect to many parameters from the solution of only two linear systems. This reduces the cost of computing the gradient to the point that the finite-collisionality neoclassical fluxes can be used within an optimization loop. With the neoclassical adjoint method, we compute solutions of the drift kinetic equation and an adjoint drift kinetic equation to obtain derivatives of neoclassical quantities with respect to geometric parameters. When the number of parameters in the derivative is large ($\mathcal{O}(10^2)$), this adjoint method provides up to a factor of 200 reduction in cost. We demonstrate adjoint-based optimization of the field strength to obtain minimal bootstrap current on a surface. With adjoint-based derivatives, we also compute the local sensitivity to magnetic perturbations on a flux surface and identify regions where tight tolerances on error fields are required for control of the bootstrap current or radial transport. Furthermore, the solve for the ambipolar electric field is accelerated using a Newton method with derivatives obtained from the adjoint method. △ Less

Submitted 1 June, 2019; v1 submitted 12 April, 2019; originally announced April 2019.

Optimized quasisymmetric stellarators are consistent with the Garren-Boozer construction

Authors: Matt Landreman

Abstract: 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 by Garren and Boozer [Phys Fluids B 3, 2805 (1991)]. Here w… ▽ More 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 by Garren and Boozer [Phys Fluids B 3, 2805 (1991)]. Here we show that the core shape and rotational transform of many optimization-based configurations can be accurately described by this direct-construction approach. This consistency supports use of the near-axis construction as an accurate analytical model for modern stellarator configurations. △ Less

Submitted 11 May, 2019; v1 submitted 5 February, 2019; originally announced February 2019.