Towards joint optimization of stellarator coils and support structures

Sintesi Tecnica: Verso l'Ottimizzazione Congiunta di Bobine Stellarator e Strutture di Supporto

Definizione del Problema
La progettazione di reattori a fusione stellarator affronta una sfida critica: bilanciare la necessità di un campo magnetico accurato per il confinamento del plasma con i requisiti meccanici per proteggere i fragili conduttori superconduttori ad alta temperatura (HTS) dalle enormi forze di Lorentz e dagli stress termici. Attualmente, il processo di progettazione è sequenziale e laborioso. Le geometrie delle bobine vengono ottimizzate per le prestazioni magnetiche in isolamento, e le strutture di supporto vengono successivamente aggiunte manualmente attraverso un'analisi iterativa agli elementi finiti (FEA) per soddisfare i limiti di stress ingegneristico. Questo approccio di "post-processing" non tiene conto di come la posizione dei supporti influenzi lo stress della bobina e di come la deformazione della bobina sotto carico introduca errori nel campo magnetico. Gli strumenti di ottimizzazione esistenti si affidano o a proxy analitici semplificati che ignorano le strutture di supporto, o utilizzano strumenti FEA commerciali a codice chiuso che non possono interfacciarsi con i cicli di ottimizzazione della bobina basati sul gradiente.

Metodologia
Gli autori presentano coil-fem, uno strumento di meccanica strutturale differenziabile open-source integrato nel ciclo di ottimizzazione della bobina stellarator. Costruito sulla libreria JAX-FEM, coil-fem consente l'ottimizzazione simultanea della geometria della bobina e della posizione delle strutture di supporto utilizzando metodi basati sul gradiente.

Componenti metodologici chiave:

• FEA Differenziabile: Lo strumento risolve l'equazione di elasticità lineare $\nabla \cdot \sigma(u) + F_{body} = 0$ con condizioni al contorno di Robin (fondazione a molla). Questa approssimazione sostituisce le condizioni Dirichlet fisse, permettendo alle posizioni dei morsetti di supporto di essere trattate come variabili ottimizzabili senza alterare la struttura della matrice sparsa del sistema FEA.
• Modello Fisico: Il modello tiene conto delle forze di corpo, inclusa la gravità, le forze magnetiche reciproche (modellate come filamenti 1D) e le forze auto-indotte (calcolate tramite la formula di Landreman-Hurwitz per sezioni rettangolari). Incorpora inoltre la deformazione termica derivante dal calo di circa 200K dallo stato privo di stress allo stato operativo.
• Meshing e Geometria: Le bobine sono modellate come corpi a costruzione finita generati mediante l'estrusione di una sezione trasversale rettangolare lungo una linea centrale utilizzando un frame minimizzante la rotazione (rotation-minimizing frame). La topologia della mesh rimane fissa durante l'ottimizzazione per garantire differenziabilità ed efficienza computazionale.
• Ciclo di Ottimizzazione: Gli autori formulano una funzione di penalità multi-obiettivo che include l'errore del campo magnetico, il carico meccanico (integrale volumetrico dello stress di von Mises al quadrato) e i vincoli geometrici standard (lunghezza della bobina, spaziatura, curvatura). L'algoritmo L-BFGS viene utilizzato per minimizzare questo obiettivo rispetto sia ai coefficienti di Fourier della bobina che alle posizioni parametriche dei morsetti di supporto.

Contributi Chiave

1. Prima Ottimizzazione Combinata: Questo articolo presenta il primo caso in letteratura di un'ottimizzazione guidata dal gradiente e combinata della geometria della bobina stellarator e della posizione delle strutture di supporto.
2. Strumento coil-fem: Il rilascio di uno strumento open-source che integra la FEA differenziabile del supporto direttamente nel ciclo di ottimizzazione della bobina, trattando le posizioni dei morsetti di supporto come variabili ottimizzabili.
3. Dimostrazione dell'Importanza del Supporto: Lo studio fornisce prove numeriche del fatto che il posizionamento delle strutture di supporto è un fattore critico per ridurre lo stress meccanico, un fattore spesso trascurato quando i supporti sono fissi o posizionati manualmente.

Risultati
Gli autori hanno confrontato coil-fem con la suite FEA DOLFINx e con l'approccio surrogate di tipo commerciale, dimostrando un accordo entro l'1–5% per le metriche di stress e spostamento. Hanno poi applicato lo strumento a un set di bobine W-XX semplificato sotto quattro scenari:

• Baseline: Bobine W-X originali con morsetti superiore/inferiore fissi.
• Caso A (Supporti Ottimizzati): Geometria della bobina fissa, posizioni dei morsetti ottimizzate.
• Caso B (Bobine Ottimizzate): Geometria della bobina ottimizzata, posizioni dei morsetti fisse.
• Caso C (Ottimizzazione Congiunta): Sia la geometria della bobina che la posizione dei morsetti sono ottimizzate.
• Caso D (Minimizzazione della Forza): Geometria della bobina ottimizzata per minimizzare la forza di Lorentz (proxy standard) con morsetti fissi.

Risultati Principali:

• Riduzione dello Stress: Il Caso C ha ottenuto una riduzione di 2,4 volte dello stress di von Mises RMS rispetto alla baseline.
• Ruolo dei Supporti: L'ottimizzazione della sola geometria della bobina con supporti fissi (Caso B) ha prodotto solo miglioramenti marginali nello stress. Al contrario, l'ottimizzazione della sola posizione dei supporti (Caso A) ha ridotto significativamente lo stress, confermando che il posizionamento dei supporti è un fattore dominante nella meccanica strutturale.
• Limiti dei Proxy: Il Caso D ha dimostrato che minimizzare le forze di Lorentz (un proxy comune) non correla necessariamente con la minimizzazione dello stress di von Mises, giustificando la necessità di una FEA strutturale integrata.
• Effetti Termici: Sebbene le forze di Lorentz siano lo stressore dominante, la contrazione termica impatta significativamente lo stress massimo (spostando il picco di circa il 17%), evidenziando la necessità di includere i carichi termici nel ciclo di progettazione.

Significato e Rivendicazioni
L'articolo afferma che coil-fem accelera significativamente il tasso di iterazione per la progettazione di stellarator eliminando il processo manuale di tentativi ed errori nel posizionamento delle strutture di supporto. Abilitando la riduzione simultanea degli errori di campo magnetico e degli stress meccanici, lo strumento offre una via verso progettazioni di reattori a fusione più economiche e robuste. Gli autori pongono questo lavoro come una "prova di concetto" che stabilisce la fattibilità della meccanica strutturale differenziabile nella progettazione di bobine per la fusione. Essi dichiarano esplicitamente che, sebbene l'attuale modello utilizzi assunzioni semplificative (ad esempio, materiali isotropi, supporti a gabbia modellati come molle), il framework consente estensioni future a modelli termo-elastici ad alta fedeltà e strutture di supporto più complesse (ad esempio, reti di travi o gusci). Il lavoro non pretende di aver risolto tutte le sfide meccaniche degli stellarator, ma piuttosto di fornire il primo framework computazionale per ottimizzare congiuntamente bobine e relativi supporti.