Global stellarator coil optimization with quadratic constraints and objectives
Most present stellarator designs are produced by costly two-stage optimization: the first for an optimized equilibrium, and the second for a coil design reproducing its magnetic configuration. Few proxies for coil complexity and forces exist at the equilibrium stage. Rapid initial state finding for...
Saved in:
Main Authors: | , , , |
---|---|
Format: | Article |
Language: | English |
Published: |
IOP Publishing
2025-01-01
|
Series: | Nuclear Fusion |
Subjects: | |
Online Access: | https://doi.org/10.1088/1741-4326/ada810 |
Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
_version_ | 1832591476874280960 |
---|---|
author | Lanke Fu Elizabeth J. Paul Alan A. Kaptanoglu Amitava Bhattacharjee |
author_facet | Lanke Fu Elizabeth J. Paul Alan A. Kaptanoglu Amitava Bhattacharjee |
author_sort | Lanke Fu |
collection | DOAJ |
description | Most present stellarator designs are produced by costly two-stage optimization: the first for an optimized equilibrium, and the second for a coil design reproducing its magnetic configuration. Few proxies for coil complexity and forces exist at the equilibrium stage. Rapid initial state finding for both stages is a topic of active research. Most present convex coil optimization codes use the least square winding surface method by Merkel (NESCOIL), with recent improvements in conditioning, regularization, sparsity, and physics objectives. While elegant, the method is limited to modeling the norms of linear functions in coil current. We present QUADCOIL, a global coil optimization method that targets combinations of linear and quadratic functions of the current. It can directly constrain and/or minimize a wide range of physics objectives unavailable in NESCOIL and REGCOIL, including the Lorentz force, magnetic energy, curvature, field-current alignment, and the maximum density of a dipole array. QUADCOIL requires no initial guess and runs nearly $10^2\times$ faster than filament optimization. Integrating it in the equilibrium optimization stage can potentially exclude equilibria with difficult-to-design coils, without significantly increasing the computation time per iteration. QUADCOIL finds the exact, global minimum in a large parameter space when possible, and otherwise finds a well-performing approximate global minimum. It supports most regularization techniques developed for NESCOIL and REGCOIL. We demonstrate QUADCOIL’s effectiveness in coil topology control, minimizing non-convex penalties, and predicting filament coil complexity with three numerical examples. |
format | Article |
id | doaj-art-d31669e4a10f429ebcad2875a1556a44 |
institution | Kabale University |
issn | 0029-5515 |
language | English |
publishDate | 2025-01-01 |
publisher | IOP Publishing |
record_format | Article |
series | Nuclear Fusion |
spelling | doaj-art-d31669e4a10f429ebcad2875a1556a442025-01-22T10:56:46ZengIOP PublishingNuclear Fusion0029-55152025-01-0165202604510.1088/1741-4326/ada810Global stellarator coil optimization with quadratic constraints and objectivesLanke Fu0https://orcid.org/0000-0002-6845-387XElizabeth J. Paul1https://orcid.org/0000-0002-9355-5595Alan A. Kaptanoglu2Amitava Bhattacharjee3Princeton Plasma Physics Laboratory , Princeton, NJ, United States of America; Department of Astrophysical Sciences, Princeton University , Princeton, NJ, United States of AmericaDepartment of Applied Physics, Columbia University , New York, NY, United States of AmericaCourant Institute of Mathematical Sciences, New York University , New York, NY, United States of AmericaDepartment of Astrophysical Sciences, Princeton University , Princeton, NJ, United States of AmericaMost present stellarator designs are produced by costly two-stage optimization: the first for an optimized equilibrium, and the second for a coil design reproducing its magnetic configuration. Few proxies for coil complexity and forces exist at the equilibrium stage. Rapid initial state finding for both stages is a topic of active research. Most present convex coil optimization codes use the least square winding surface method by Merkel (NESCOIL), with recent improvements in conditioning, regularization, sparsity, and physics objectives. While elegant, the method is limited to modeling the norms of linear functions in coil current. We present QUADCOIL, a global coil optimization method that targets combinations of linear and quadratic functions of the current. It can directly constrain and/or minimize a wide range of physics objectives unavailable in NESCOIL and REGCOIL, including the Lorentz force, magnetic energy, curvature, field-current alignment, and the maximum density of a dipole array. QUADCOIL requires no initial guess and runs nearly $10^2\times$ faster than filament optimization. Integrating it in the equilibrium optimization stage can potentially exclude equilibria with difficult-to-design coils, without significantly increasing the computation time per iteration. QUADCOIL finds the exact, global minimum in a large parameter space when possible, and otherwise finds a well-performing approximate global minimum. It supports most regularization techniques developed for NESCOIL and REGCOIL. We demonstrate QUADCOIL’s effectiveness in coil topology control, minimizing non-convex penalties, and predicting filament coil complexity with three numerical examples.https://doi.org/10.1088/1741-4326/ada810stellaratorcoilswinding surfaceoptimizationconvex relaxation |
spellingShingle | Lanke Fu Elizabeth J. Paul Alan A. Kaptanoglu Amitava Bhattacharjee Global stellarator coil optimization with quadratic constraints and objectives Nuclear Fusion stellarator coils winding surface optimization convex relaxation |
title | Global stellarator coil optimization with quadratic constraints and objectives |
title_full | Global stellarator coil optimization with quadratic constraints and objectives |
title_fullStr | Global stellarator coil optimization with quadratic constraints and objectives |
title_full_unstemmed | Global stellarator coil optimization with quadratic constraints and objectives |
title_short | Global stellarator coil optimization with quadratic constraints and objectives |
title_sort | global stellarator coil optimization with quadratic constraints and objectives |
topic | stellarator coils winding surface optimization convex relaxation |
url | https://doi.org/10.1088/1741-4326/ada810 |
work_keys_str_mv | AT lankefu globalstellaratorcoiloptimizationwithquadraticconstraintsandobjectives AT elizabethjpaul globalstellaratorcoiloptimizationwithquadraticconstraintsandobjectives AT alanakaptanoglu globalstellaratorcoiloptimizationwithquadraticconstraintsandobjectives AT amitavabhattacharjee globalstellaratorcoiloptimizationwithquadraticconstraintsandobjectives |