Augmented Lagrangian methods produce cutting-edge magnetic coils for stellarator fusion reactors

Technische Samenvatting: Augmented Lagrangian-methoden voor de optimalisatie van stellarator-spoelen

Probleemstelling
Het ontwerp van magnetische spoelen voor stellarator-fusiereactoren vormt een kritieke flessenhals bij het realiseren van gecontroleerde kernfusie. Hoewel stellarators een stationaire werking bieden zonder de noodzaak van geïnduceerde plasmastromen (in tegen tegenstelling tot tokamaks), vereisen hun complexe driedimensionale magnetische velden ingewikquete externe spoelsystemen. De optimalisatie van deze spoelen is een slecht gedefinieerd (ill-posed), sterk niet-convex invers probleem. Traditionele ontwerpprocessen, die de optimalisatie van de plasmabegrenzing (Fase I) scheiden van de optimalisatie van de spoelvorm (Fase II), worstelen met deze complexiteit.

Huidige methodologieën vertrouwen doorgaans op het transformeren van het beperkte optimalisatieprobleem naar een onbeperkt multi-objectief probleem door handmatige gewichten toe te kennen aan diverse doelstellingen en beperkingen (bijv. veldaccuratesse, spoelkromming, afstand tussen spoelen). Deze aanpak lijdt onder verschillende beperkingen:

1. Handmatige afstemming: De gewichten missen fysieke betekenis en vereisen uitgebreide, computationeel dure handmatige afstemming om haalbare oplossingen te vinden.
2. Lokale minima: Het niet-convexe landschap vangt optimalisatie-algoritmen vaak in suboptimale configuraties.
3. Technische onhaalbaarheid: Veel gegenereerde spoelsets zijn onpraktisch vanwege excessieve kromming, verstrengeling of krachten die de materiaaltoleranties overschrijden, wat leidt tot vertragingen in de bouw en kostenoverschrijdingen in projecten zoals Wendelstein 7-X en de annulering van NCSX.

Methodologie
De auteurs stellen een Augmented Lagrangian (AL) methode voor, geïmplementeerd binnen de SIMSOPT-optimalisatiesuite om deze uitdagingen aan te pakken. In tegen tegenstelling tot traditionele gewogen-som-benaderingen, behandelt de AL-methode beperkingen dynamisch.

• Formalisme: De optimalisatie minimaliseert een Lagrange-functie $L_A(x, \lambda, \mu) = f(x) - \lambda^\top c(x) + \frac{1}{2} \|\sqrt{\mu} \circ c(x)\|_2^2$, waarbij $f(x)$ de doelstelling is, $c(x)$ de beperkingen vertegenwoordigt, $\lambda$ de Lagrange-multiplicatoren zijn en $\mu$ de strafparameters (penalty parameters) zijn.
• Dynamische updates: In plaats van vaste, door de gebruiker gedefinieerde gewichten, worden de Lagrange-multiplicatoren ($\lambda$) en strafparameters ($\mu$) automatisch bijgewerkt tijdens het optimalisatieproces op basis van schendingen van de beperkingen. Dit elimineert de noodzaak voor handmatige gewichtstuning.
• Behandeling van beperkingen: De methode behandelt de gekwadrateerde fluxfout (een maat voor de afwijking van het magnetisch veld) als een beperking in plaats van de primaire doelstelling. Door een drempelwaarde in te stellen (bijv. $f_{SF} \lesssim 10^{-6}$), geeft het algoritme prioriteit aan het voldoen aan technische beperkingen zodra de veldaccuratesse voldoende is, in plaats van na te streven naar een theoretisch perfecte maar praktisch onmogelijke oplossing met nul-fout.
• Technische beperkingen: Het raamwerk bevat standaardmetrieken voor spoelcomplexiteit, waaronder minimale afstand tussen spoelen en tussen spoel en oppervlak, totale spoellengte, krommingslimieten, gemiddelde kwadratische kromming, Gauss-verknopingsgetallen (om verstrengeling te voorkomen) en elektromagnetische krachten.
• Optimalisatiealgoritme: De innerlijke minimalisatiestap maakt gebruik van de L-BFGS-B quasi-Newton methode. Gradiënten worden berekend met een combinatie van analytische afgeleiden en automatische differentiatie via vector-Jacobiaan producten.

Belangrijkste bijdragen en resultaten
De paper demonstreert de effectiviteit van de AL-methode door spoelsets te genereren voor vijf verschillende stellarator-configuraties, variërend van theoretische concepten tot volledig gebouwde apparaten. In alle gevallen produceerde de methode Pareto-optimale oplossingen die voorheen gepubliceerde sets evenaren of overtreffen zonder handmatige gewichtstuning.

1. Quasi-axiale symmetrische (QA) stellarator:

   • De methode bracht de Pareto-front in kaart tussen spoelkrachten en veldaccuratesse voor een Landreman&Paul QA-configuratie.
   • Het genereerde een 3-spoelen per half-veldperiode oplossing die een referentie 4-spoelen oplossing (Wechsung et al.) overtrof in veldaccuratesse en afstand tussen spoelen, ondanks de verwachting dat minder spoelen de ripple zouden vergroten.
   • Het traceerde succesvol nieuwe Pareto-grenzen door de beperkingen op totale lengte te versoepelen, waarbij veldaccuratessen tot een orde van grootte beter werden bereikt dan eerdere scans.

2. Quasi-helicale symmetrische (QH) stellarator:

   • Toegepast op een Landreman&Paul QH-configuratie, produceerde de methode een 4-spoelen oplossing die de fysica-prestaties van een referentie 5-spoelen oplossing (Wiedman et al.) evenaarde, terwijl de technische metrieken verbeterden (bijv. gereduceerde spoellengte, grotere afstand).
   • Een derde configuratie, die puur werd geoptimaliseerd voor technische eigenschappen, toonde een toename van 38% in de afstand tussen spoel en oppervlak en een reductie van 10% in HTS-materiaal kosten, zij het met een hogere quasisymmetrische fout, wat de mogelijkheid aantoont om fysica-getrouwheid in te ruilen voor maakbaarheid.

3. Stellaris (SQuID) configuratie:

   • Voor het Proxima Fusion Stellaris-concept vond de methode alternatieven met 5 en 6 spoelen per half-veldperiode.
   • Deze oplossingen verbeterden de magnetische veldfout met 25% en vergrootten de afstand tussen de spoelen met 13 cm vergeleken met het oorspronkelijke ontwerp, terwijl de piekkrachten op HTS-wikkelingen met 29% werden verminderd.
   • Een configuratie die geoptimaliseerd was voor technische eigenschappen behaalde een reductie van 6% in de geschatte HTS-kosten en aanzienlijk eenvoudigere geometrieën.

4. Wendelstein 7-X (W7-X):

   • De methode genereerde alternatieve spoelsets voor de bestaande W7-X reactor.
   • Configuratie #1: Verminderde het aantal spoelen van 50 naar 40 (5 naar 4 per half-veldperiode) terwijl de veldaccuratesse verbeterde en de totale spoellengte afnam.
   • Configuratie #2: Bereikte een reductie van 20% in elektromagnetische krachten met slechts een 5% toename in totale spoellengte.
   • Configuratie #3: Verbeterde de veldfout en kromming drastisch met verwaarloosbare veranderingen in lengte of krachten.
   • Alle alternatieven vertoonden lagere kromming en gladdere geometrieën dan het oorspronkelijke ontwerp.

5. Helically Symmetric eXperiment (HSX):

   • Voor het bestaande HSX-apparaat produceerde de methode een 4-spoelen oplossing die de totale spoellengte met 25% verminderde en de kromming met 45% vergeleken met het originele 6-spoelen ontwerp, terwijl een redelijke veldaccuratesse behouden bleef. Dit benadrukt het potentieel om de fabricage van onderzoeksschaal-apparaten te vereenvoudigen.

Betekenis en claims
De paper beweert dat de Augmented Lagrangian-benadering de stand van de techniek in stellarator-spoelontwerp fundamenteel vooruithelpt door:

• Automatisering van gewichtstuning: Het verwijderen van de afhankelijkheid van heuristische, handmatige gewichtsaanpassingen, die een belangrijke bron van inefficiëntie en suboptimaliteit zijn in huidige workflows.
• Robuustheid: Het succesvol navigeren door het niet-convexe optimalisatielandschap om haalbare, Pareto-optimale oplossingen te vinden voor diverse magnetische configuraties (QA, QH, QI).
• Veelzijdigheid: Het aantonen van toepasbaarheid op zowel conceptuele ontwerpen als bestaande, volledig gefabriceerde reactoren (W7-X, HSX), wat bewijst dat de methode gevestigde technische oplossingen kan verbeteren.
• Efficiëntie: Het behalen van deze resultaten met aanzienlijk minder computationele overhead dan de enorme parameterscans (bijv. 8.500 runs) die vereist zijn door eerdere stochastische of gewogen-som-methoden.

De auteurs concluderen dat deze methodologie een computationeel efficiënte en robuuste route biedt voor het genereren van spoelsets die de eisen van de plasmafysica in evenwicht brengen met de strikte technische beperkingen die nodig zijn voor integratie op reactorniveau.