Andrew Giuliani

None, since per your instructions on the web, Notes should not have abstracts.

The non-resonant divertor (NRD) offers a promising exhaust solution for stellarators, combining topological simplicity with resilience to magnetic field perturbations. To experimentally validate the robustness of NRDs in a quasi-axisymmetric (QA) configuration, we introduce STAR_Lite, a new stellarator experiment at Hampton University. This paper details the design and analysis of the first STAR_Lite coil configuration, STAR_Lite-A. The two-field-period configuration manifests an NRD through X-points with zero rotational transform, at the top and bottom of the device. The divertor legs extruding from the X-points are topologically similar to the poloidal divertors of tokamaks. To expand the experimental range, STAR_Lite-A is optimised for experimental flexibility, producing a wide range of distinct QA configurations by only varying the currents in the modular coils. The NRDs not only persist across these configurations, but numerical strike-line simulations confirm that heat exhaust remains resilient to changes in coil currents, with plasma following the divertor legs and creating a toroidal, discontinuous, strike pattern. We further examine the resilience of the NRD to magnetic perturbations caused by manufacturing errors in the modular coils. We find that quasi-symmetry and the existence of X-points are well preserved under these magnetic field changes, but the rotational transform may vary substantially and displacements of the divertor X-points may lead to one X-point having a dominant effect on edge transport. Overall, our analysis indicates that a compact, modular design can likely generate a resilient NRD structure while satisfying the practical constraints of a university-scale experiment.

We augment the ‘quasisisymmetric stellarator repository’ (QUASR) to include vacuum field stellarators with quasihelical symmetry using a globalized optimization workflow. The database now has over 300 000 quasisaxisymmetry and quasihelically symmetric devices along with coil sets, optimized for a variety of aspect ratios, rotational transforms and discrete rotational symmetries. This paper outlines a couple of ways to explore and characterize the data set. We plot devices on a near-axis quasisymmetry (QS) landscape, revealing close correspondence to this predicted landscape. We also use principal component analysis (PCA) to reduce the dimensionality of the data so that it can easily be visualized in two or three dimensions. The PCA also gives a mechanism to compare the new devices here with previously published ones in the literature. We are able to characterize the structure of the data, observe clusters and visualize the progression of devices in these clusters. The topology of the data are governed by the interplay of the design constraints and valleys of the QS objective. These techniques reveal that the data has structure, and that typically one, two or three principal components are sufficient to characterize it. The latest version of QUASR is archived at https://zenodo.org/doi/10.5281/zenodo.10050655 and can be explored online at quasr.flatironinstitute.org.

Single-stage optimization, also known as combined plasma-coil algorithms or direct coil optimization, has recently emerged as a possible method to expedite the design of stellarator devices by including, in a single step, confinement, stability, and engineering constraints. In this work, we show how such frameworks allow us to find new designs in a streamlined manner, yielding a broad range of new configurations. Examples are shown for stellarators with a small number of coils and quasisymmetric stellarators with only one to three coils per half field-period, with external trim coils, helical coils, and a single set of coils generating both a quasi-axisymmetric and a quasi-helical equilibrium.

We propose a practical finite volume method on cut cells using state redistribution. Our algorithm is provably monotone, total variation diminishing, and GKS (Gustafsson, Kreiss, Sundström) stable in many situations, and shuts off continuously as the cut cell size approaches a target value. Our analysis reveals why original state redistribution works so well: it results in a monotone scheme for most configurations, though at times subject to a slightly smaller CFL condition. Our analysis also explains why a premerging step is beneficial. We show computational experiments in two and three dimensions.

Many stellarator coil design problems are plagued by multiple minima, where the locally optimal coil sets can sometimes vary substantially in performance. As a result, solving a coil design problem a single time with a local optimization algorithm is usually insufficient and better optima likely do exist. To address this problem, we propose a global optimization algorithm for the design of stellarator coils and outline how to apply box constraints to the physical positions of the coils. The algorithm has a global exploration phase that searches for interesting regions of design space and is followed by three local optimization algorithms that search in these interesting regions (a ``global-to-local" approach). The first local algorithm (phase I), following the globalization phase, is based on near-axis expansions and finds stellarator coils that optimize for quasisymmetry in the neighborhood of a magnetic axis. The second local algorithm (phase II) takes these coil sets and optimizes them for nested flux surfaces and quasisymmetry on a toroidal volume. The final local algorithm (phase III) polishes these configurations for an accurate approximation of quasisymmetry. Using our global algorithm, we study the trade-off between coil length, aspect ratio, rotational transform, and quality of quasi-axisymmetry. The database of stellarators, which comprises approximately 200,000 coil sets, is available online and is called QUASR, for `QUAsi-symmetric Stellarator Repository'.

We propose a practical finite volume method on cut cells using state redistribution. Our algorithm is provably monotone, total variation diminishing, and GKS stable in many situations, and shuts off continuously as the cut cell size approaches a target value. Our analysis reveals why original state redistribution works so well: it results in a monotone scheme for most configurations, though at times subject to a slightly smaller CFL condition. Our analysis also explains why a pre-merging step is beneficial. We show computational experiments in two and three dimensions.

Many stellarator coil design problems are plagued by multiple minima, where the locally optimal coil sets can sometimes vary substantially in performance. As a result, solving a coil design problem a single time with a local optimization algorithm is usually insufficient and better optima likely do exist. To address this problem, we propose a global optimization algorithm for the design of stellarator coils and outline how to apply box constraints to the physical positions of the coils. The algorithm has a global exploration phase that searches for interesting regions of design space and is followed by three local optimization algorithms that search in these interesting regions (a ‘global-to-local’ approach). The first local algorithm (phase I), following the globalization phase, is based on near-axis expansions and finds stellarator coils that optimize for quasisymmetry in the neighbourhood of a magnetic axis. The second local algorithm (phase II) takes these coil sets and optimizes them for nested flux surfaces and quasisymmetry on a toroidal volume. The final local algorithm (phase III) polishes these configurations for an accurate approximation of quasisymmetry. Using our global algorithm, we study the trade-off between coil length, aspect ratio, rotational transform and quality of quasi-axisymmetry. The database of stellarators, which comprises approximately 200 000 coil sets, is available online and is called QUASR, for ‘quasi-symmetric stellarator repository’.

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.

It was recently shown in Wechsung et al (2022 Proc. Natl Acad. Sci. USA 119 e2202084119) 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 possible 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.