-
Exploration of the parameter space of quasisymmetric stellarator vacuum fields through adjoint optimisation
Authors:
Richard Nies,
Elizabeth J. Paul,
Dario Panici,
Stuart R. Hudson,
Amitava Bhattacharjee
Abstract:
Optimising stellarators for quasisymmetry leads to strongly reduced collisional transport and energetic particle losses compared to unoptimised configurations. Though stellarators with precise quasisymmetry have been obtained in the past, it remains unclear how broad the parameter space is where good quasisymmetry may be achieved. We study the range of aspect ratio and rotational transform values…
▽ More
Optimising stellarators for quasisymmetry leads to strongly reduced collisional transport and energetic particle losses compared to unoptimised configurations. Though stellarators with precise quasisymmetry have been obtained in the past, it remains unclear how broad the parameter space is where good quasisymmetry may be achieved. We study the range of aspect ratio and rotational transform values for which stellarators with excellent quasisymmetry on the boundary can be obtained. A large number of Fourier harmonics is included in the boundary representation, which is made computationally tractable by the use of adjoint methods to enable fast gradient-based optimisation, and by the direct optimisation of vacuum magnetic fields, which converge more robustly compared to solutions from magnetohydrostatics. Several novel configurations are presented, including stellarators with record levels of quasisymmetry on a surface; three field period quasiaxisymmetric stellarators with substantial magnetic shear, and compact quasisymmetric stellarators at low aspect ratios similar to tokamaks.
△ Less
Submitted 23 May, 2024; v1 submitted 2 April, 2024;
originally announced April 2024.
-
Stellarator Optimization with Constraints
Authors:
Rory Conlin,
Patrick Kim,
Daniel W. Dudt,
Dario Panici,
Egemen Kolemen
Abstract:
In this work we consider the problem of optimizing a stellarator subject to hard constraints on the design variables and physics properties of the equilibrium. We survey current numerical methods for handling these constraints, and summarize a number of methods from the wider optimization community that have not been used extensively for stellarator optimization thus far. We demonstrate the utilit…
▽ More
In this work we consider the problem of optimizing a stellarator subject to hard constraints on the design variables and physics properties of the equilibrium. We survey current numerical methods for handling these constraints, and summarize a number of methods from the wider optimization community that have not been used extensively for stellarator optimization thus far. We demonstrate the utility of new methods of constrained optimization by optimizing a QA stellarator for favorable physics properties while preventing strong shaping of the plasma boundary which can be difficult to create with external current sources.
△ Less
Submitted 16 March, 2024;
originally announced March 2024.
-
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.
-
Magnetic Fields with General Omnigenity
Authors:
Daniel W. Dudt,
Alan G. Goodman,
Rory Conlin,
Dario Panici,
Egemen Kolemen
Abstract:
Omnigenity is a desirable property of toroidal magnetic fields that ensures confinement of trapped particles. Confining charged particles is a basic requirement for any fusion power plant design, but it can be difficult to satisfy with the non-axisymmetric magnetic fields used by the stellarator approach. Every ideal magnetohydrodynamic equilibrium previously found to approximate omnigenity has be…
▽ More
Omnigenity is a desirable property of toroidal magnetic fields that ensures confinement of trapped particles. Confining charged particles is a basic requirement for any fusion power plant design, but it can be difficult to satisfy with the non-axisymmetric magnetic fields used by the stellarator approach. Every ideal magnetohydrodynamic equilibrium previously found to approximate omnigenity has been either axisymmetric, quasi-symmetric or has poloidally closed contours of magnetic field strength $B$. However, general omnigenous equilibria are a much larger design space than these subsets. A new model is presented and employed in the DESC stellarator optimization suite to represent and discover the full parameter space of omnigenous equilibria. Although exact omnigenity aside from quasi-symmetry is impossible, these results reveal that excellent particle confinement can be achieved in practice. Examples far from quasi-symmetry with poloidally, helically and toroidally closed $B$ contours are attained with DESC and shown to have low neoclassical collisional transport and fast particle losses.
△ Less
Submitted 3 April, 2024; v1 submitted 13 May, 2023;
originally announced May 2023.
-
The DESC Stellarator Code Suite Part III: Quasi-symmetry optimization
Authors:
Daniel Dudt,
Rory Conlin,
Dario Panici,
Egemen Kolemen
Abstract:
The DESC stellarator optimization code takes advantage of advanced numerical methods to search the full parameter space much faster than conventional tools. Only a single equilibrium solution is needed at each optimization step thanks to automatic differentiation, which efficiently provides exact derivative information. A Gauss-Newton trust-region optimization method uses second-order derivative i…
▽ More
The DESC stellarator optimization code takes advantage of advanced numerical methods to search the full parameter space much faster than conventional tools. Only a single equilibrium solution is needed at each optimization step thanks to automatic differentiation, which efficiently provides exact derivative information. A Gauss-Newton trust-region optimization method uses second-order derivative information to take large steps in parameter space and converges rapidly. With just-in-time compilation and GPU portability, high-dimensional stellarator optimization runs take orders of magnitude less computation time with DESC compared to other approaches. This paper presents the theory of the DESC fixed-boundary local optimization algorithm along with demonstrations of how to easily implement it in the code. Example quasi-symmetry optimizations are shown and compared to results from conventional tools. Three different forms of quasi-symmetry objectives are available in DESC, and their relative advantages are discussed in detail. In the examples presented, the triple product formulation yields the best optimization results in terms of minimized computation time and particle transport. This paper concludes with an explanation of how the modular code suite can be extended to accommodate other types of optimization problems.
△ Less
Submitted 31 October, 2022; v1 submitted 31 March, 2022;
originally announced April 2022.
-
The DESC Stellarator Code Suite Part I: Quick and accurate equilibria computations
Authors:
Dario Panici,
Rory Conlin,
Daniel W. Dudt,
Kaya Unalmis,
Egemen Kolemen
Abstract:
3D equilibrium codes are vital for stellarator design and operation, and high-accuracy equilibria are also necessary for stability studies. This paper details comparisons of two three-dimensional equilibrium codes, VMEC, which uses a steepest-descent algorithm to reach a minimum-energy plasma state, and DESC, which minimizes the magnetohydrodynamic (MHD) force error in real space directly. Accurac…
▽ More
3D equilibrium codes are vital for stellarator design and operation, and high-accuracy equilibria are also necessary for stability studies. This paper details comparisons of two three-dimensional equilibrium codes, VMEC, which uses a steepest-descent algorithm to reach a minimum-energy plasma state, and DESC, which minimizes the magnetohydrodynamic (MHD) force error in real space directly. Accuracy as measured by satisfaction of MHD force balance is presented for each code, along with the computation time. It is shown that DESC is able to achieve more accurate solutions, especially near-axis. The importance of higher accuracy equilibria is shown in DESC's better agreement of stability metrics with asymptotic formulae. DESC's global Fourier-Zernike basis also yields the solution with analytic derivatives explicitly everywhere in the plasma volume, provides improved accuracy in the radial direction versus conventional finite differences, and allows for exponential convergence. Further, DESC can compute the same accuracy solution as VMEC in an order of magnitude less time.
△ Less
Submitted 31 March, 2023; v1 submitted 31 March, 2022;
originally announced March 2022.
-
The DESC Stellarator Code Suite Part II: Perturbation and continuation methods
Authors:
Rory Conlin,
Daniel W. Dudt,
Dario Panici,
Egemen Kolemen
Abstract:
A new perturbation and continuation method is presented for computing and analyzing stellarator equilibria. The method is formally derived from a series expansion about the equilibrium condition $F \equiv J \times B - \nabla p = 0$, and an efficient algorithm for computing solutions to 2nd and 3rd order perturbations is developed. The method has been implemented in the DESC stellarator equilibrium…
▽ More
A new perturbation and continuation method is presented for computing and analyzing stellarator equilibria. The method is formally derived from a series expansion about the equilibrium condition $F \equiv J \times B - \nabla p = 0$, and an efficient algorithm for computing solutions to 2nd and 3rd order perturbations is developed. The method has been implemented in the DESC stellarator equilibrium code, using automatic differentiation to compute the required derivatives. Examples are shown demonstrating its use for computing complicated equilibria, perturbing a tokamak into a stellarator, and performing parameter scans in pressure, rotational transform and boundary shape in a fraction of the time required for a full solution.
△ Less
Submitted 4 April, 2023; v1 submitted 29 March, 2022;
originally announced March 2022.