Towards joint optimization of stellarator coils and support structures

Technische Samenvatting: Naar gezamenlijke optimalisatie van stellarator-spoelen en ondersteuningsstructuren

Probleemstelling
Het ontwerp van stellarator-fusiereactoren staat voor een cruciale uitdaging: het balanceren van de noodzaak voor een nauwkeurig magnetisch veld voor plasma-insluiting met de mechanische vereisten om brosse hoogtemperatuur-supergeleidende (HTS) geleiders te beschermen tegen enorme Lorentz-krachten en thermische spanningen. Momenteel is het ontwerpproces sequentieel en arbeidsintensief. Spoelgeometrieën worden in isolatie geoptimaliseerd voor magnetische prestaties, waarna ondersteuningsstructuren handmatig worden toegevoegd via iteratieve eindige-elementenanalyse (FEA) om aan technische spanningslimieten te voldoen. Deze "post-processing"-benadering houdt geen rekening met hoe de locatie van ondersteuning de spoelspanning beïnvloedt en hoe de deformatie van de spoel onder belasting magnetische veldfouten introduceert. Bestaande optimalisatietools vertrouwen ofwel op vereenvoudigde analytische proxies die ondersteuningsstructuren negeren, of maken gebruik van gesloten broncode commerciële FEA-tools die niet kunnen koppelen met gradiëntgebaseerde spoeloptimalisatielussen.

Methodologie
De auteurs presenteren coil-fem, een open-source, differentieerbare structurele mechanica-tool die geïntegreerd is in de stellarator-spoeloptimalisatielus. Gebouwd op de JAX-FEM bibliotheek, maakt coil-fem gelijktijdige optimalisatie van de spoelgeometrie en de locaties van de ondersteuningsstructuren mogelijk met behulp van gradiëntgebaseerde methoden.

Belangrijke methodologische componenten zijn:

• Differentieerbare FEA: De tool lost de lineaire elasticiteitsvergelijking $\nabla \cdot \sigma(u) + F_{body} = 0$ op met Robin (veersysteem) randvoorwaarden. Deze benadering vervangt vaste Dirichlet-randvoorwaarden, waardoor de locaties van de klemmen van de ondersteuning kunnen worden behandeld als optimaliseerbare variabelen zonder de ijle matrixstructuur van het FEA-systeem te wijzigen.
• Fysisch model: Het model houdt rekening met volumekrachten, waaronder zwaartekracht, onderlinge magnetische krachten (gemodelleerd als 1D-filamenten) en zelfkrachten (berekend via de Landreman-Hurwitz formule voor rechthoekige dwarsdoorsneden). Het incorporeert ook thermische rek die voortvloeit uit de temperatuurval van $\sim$200K van de spanningsvrije staat naar de operationele staat.
• Meshing en Geometrie: Spoelen worden gemodelleerd als eindige lichamen gegenereerd door een rechthoekige dwarsdoorsnede langs een middellijn te "sweepen" met een rotatie-minimaliserend frame. De mesh-topologie blijft vast tijdens de optimalisatie om differentieerbaarheid en computationele efficiëntie te waarborgen.
• Optimalisatielus: De auteurs formuleren een multi-objective penalty-functie die magnetische veldfout, mechanische belasting (volumintegraal van het kwadraat van de von Mises-spanning) en standaard geometrische beperkingen (spoellengte, tussenruimte, kromming) bevat. Het L-BFGS algoritme wordt gebruikt om deze doelstelling te minimaliseren met betrekking tot zowel de Fourier-coëfficiënten van de spoel als de parametrische locaties van de klemmen van de ondersteuning.

Belangrijkste Bijdragen

1. Eerste Gecombineerde Optimalisatie: Dit artikel presenteert het eerste geval in de literatuur van een gradiëntgestuurde, gecombineerde optimalisatie van de stellarator-spoelgeometrie en de locaties van de ondersteuningsstructuren.
2. coil-fem Tool: De release van een open-source tool die ondersteunende differentieerbare FEA direct integreert in de spoeloptimalisatielus, waarbij de locaties van de klemmen als optimaliseerbare variabelen worden behanden.
3. Demonstratie van het Belang van Ondersteuning: De studie levert numeriek bewijs dat de plaatsing van ondersteuningsstructuren een kritieke factor is bij het verminderen van mechanische spanning, een factor die vaak over het hoofd wordt gezien wanneer ondersteuningen vaststaan of handmatig worden geplaatst.

Resultaten
De auteurs hebben coil-fem getest tegen de DOLFINx FEA-suite en de commerciële-stijl surrogate-benadering, waarbij een overeenkomst van 1–5% werd aangetoond voor spanning- en verplaatsingsmetrieken. Ze hebben de tool vervolgens toegepast op een vereenvoudigde W7-X spoelset onder vier scenario's:

• Baseline: Originele W7-X spoelen met vaste bovenste/onderste klemmen.
• Geval A (Geoptimaliseerde Ondersteuning): Vaste spoelgeometrie, geoptimaliseerde klemlocaties.
• Geval B (Geoptimaliseerde Spoelen): Geoptimaliseerde spoelgeometrie, vaste klemlocaties.
• Geval C (Gezamenlijke Optimalisatie): Zowel de spoelgeometrie als de klemlocaties geoptimaliseerd.
• Geval D (Krachtminimalisatie): Geoptimaliseerde spoelgeometrie om de Lorentz-kracht te minimaliseren (standaard proxy) met vaste klemmen.

Belangrijkste Bevindingen:

• Spanningsreductie: Geval C bereikte een 2,4× reductie in de RMS von Mises-spanning vergeleken met de baseline.
• Rol van Ondersteuning: Alleen de spoelgeometrie optimaliseren met vaste ondersteuningen (Geval B) leverde slechts marginale verbeteringen in spanning op. Omgekeerd zorgde het optimaliseren van alleen de ondersteuningslocaties (Geval A) voor een significante vermindering van de spanning, wat bevestigt dat de plaatsing van ondersteuning een dominante factor is in de structurele mechanica.
• Beperkingen van Proxies: Geval D toonde aan dat het minimaliseren van Lorentz-krachten (een veelgebruikte proxy) niet noodzakelijkerwijs correleert met het minimaliseren van de von Mises-spanning, wat de noodzaak voor geïntegreerde structurele FEA rechtvaardigt.
• Thermische Effecten: Hoewel Lorentz-krachten de dominante stressfactor zijn, heeft thermische contractie een aanzienlijke impact op de maximale spanning (het verschuift de piek met $\sim$17%), wat het belang van het opnemen van thermische belastingen in de ontwerplus benadrukt.

Betekenis en Claims
Het artikel beweert dat coil-fem de iteratiesnelheid voor stellaratorontwerp aanzienlijk versnelt door het handmatige proces van trial-and-error bij het plaatsen van ondersteuningsstructuren te elimineren. Door de gelijktijdige vermindering van magnetische veldfouten en mechanische spanningen mogelijk te maken, biedt de tool een pad naar meer economische en robuuste fusiereactorontwerpen. De auteurs positioneren dit werk als een "proof of concept" die de haalbaarheid van differentieerbare structurele mechanica in fusie-spoelontwerp vaststelt. Ze stellen expliciet dat hoewel het huidige model gebruikmaakt van vereenvoudigende aannames (bijv. isotrope materialen, kooi-type ondersteuningen gemodelleerd als veren), het framework ruimte biedt voor toekomstige uitbreidingen naar hogere-getrouwheid thermo-elastische modellen en complexe ondersteuningsstructuren (bijv. balknetwerken of schalen). Het werk claimt niet alle stellaire mechanische uitdagingen te hebben opgelost, maar biedt het eerste computationele framework om spoelen en hun ondersteuningen gezamenlijk te optimaliseren.