-
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 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.