Towards joint optimization of stellarator coils and support structures

Resumo Técnico: Rumo à Otimização Conjunta de Bobinas e Estruturas de Suporte de Stellarators

Definição do Problema
O projeto de reatores de fusão stellarator enfrenta um desafio crítico: equilibrar a necessidade de um campo magnético preciso para o confinamento do plasma com os requisitos mecânicos de proteção dos condutores de supercondutores de alta temperatura (HTS) frágeis contra enormes forças de Lorentz e tensões térmicas. Atualmente, o processo de projeto é sequencial e trabalhoso. As geometrias das bobinas são otimizadas para o desempenho magnético de forma isolada, e as estruturas de suporte são adicionadas posteriormente de forma manual através de análise de elementos finitos (FEA) iterativa para atender aos limites de tensão de engenharia. Essa abordagem de "pós-processamento" não leva em conta como a localização dos suportes influencia a tensão da bobina e como a deformação da bobina sob carga introduz erros no campo magnético. As ferramentas de otimização existentes ou dependem de proxies analíticos simplificados que ignoram as estruturas de suporte ou utilizam ferramentas de FEA comerciais de código fechado que não podem se integrar a loops de otimização de bobinas baseados em gradiente.

Metodologia
Os autores apresentam o coil-fem, uma ferramenta de mecânica estrutural diferenciável de código aberto integrada ao loop de otimização de bobinas de stellarator. Construído sobre a biblioteca JAX-FEM, o coil-fem permite a otimização simultânea da geometria da bobina e da localização das estruturas de suporte usando métodos baseados em gradiente.

Componentes metodológicos principais incluem:

• FEA Diferenciável: A ferramenta resolve a equação de elasticidade linear $\nabla \cdot \sigma(u) + F_{body} = 0$ com condições de contorno de Robin (fundação de mola). Esta aproximação substitui as condições de Dirichlet fixas, permitindo que as localizações de fixação dos suportes sejam tratadas como variáveis otimizáveis sem alterar a estrutura de matriz esparsa do sistema FEA.
• Modelo de Física: O modelo considera forças de corpo, incluindo gravidade, forças magnéticas mútuas (modeladas como filamentos 1D) e forças próprias (calculadas via fórmula de Landreman-Hurwitz para seções transversais retangulares). Também incorpora a deformação térmica resultante da queda de $\sim$200K do estado livre de tensão para a operação.
• Malha e Geometria: As bobinas são modeladas como corpos de construção finita gerados pelo varrimento de uma seção transversal retangular ao longo de uma linha central usando um frame de rotação mínima. A topologia da malha é fixa durante a otimização para garantir diferenciabilidade e eficiência computacional.
• Loop de Otimização: Os autores formulam uma função de penalidade multiobjetivo que inclui erro de campo magnético, carga mecânica (integral de volume do quadrado da tensão de von Mises) e restrições geométricas padrão (comprimento da bobina, espaçamento, curvatura). O algoritmo L-BFGS é usado para minimizar este objetivo em relação tanto aos coeficientes de Fourier da bobina quanto às localizações paramétricas de fixação dos suportes.

Principais Contribuições

1. Primeira Otimização Combinada: Este artigo apresenta o primeiro caso na literatura de uma otimização conjunta, impulsionada por gradiente, da geometria da bobina de stellarator e da localização da estrutura de suporte.
2. Ferramenta coil-fem: O lançamento de uma ferramenta de código aberto que integra o suporte de FEA diferenciável diretamente no loop de otimização da bobina, tratando as localizações de fixação dos suportes como variáveis otimizáveis.
3. Demonstração da Importância do Suporte: O estudo fornece evidência numérica de que o posicionamento das estruturas de suporte é um fator crítico para reduzir a tensão mecânica, um fator frequentemente negligenciado quando os suportes são fixos ou colocados manualmente.

Resultos
Os autores testaram o coil-fem contra o conjunto de FEA DOLFINx e a abordagem de surrogate de estilo comercial, demonstrando concordância de 1–5% para métricas de tensão e deslocamento. Eles então aplicaram a ferramenta a um conjunto de bobinas W-X simplificado sob quatro cenários:

• Baseline: Bobinas originais do W-X com grampos superior/inferior fixos.
• Caso A (Suportes Otimizados): Geometria da bobina fixa, localizações de grampo otimizadas.
• Caso B (Bobinas Otimizadas): Geometria da bobina otimizada, localizações de grampo fixas.
• Caso C (Otimização Conjunta): Tanto a geometria da bobina quanto as localizações de grampo otimizadas.
• Caso D (Minimização de Força): Geometria da bobina otimizada para minimizar a força de Lorentz (proxy padrão) com grampos fixos.

Principais Descobertas:

• Redução de Tensão: O Caso C alcançou uma redução de 2,4× na tensão RMS de von Mises em comparação com o baseline.
• Papel dos Suportes: Otimizar apenas a geometria da bobina com suportes fixos (Caso B) rendeu apenas melhorias marginais na tensão. Por outro lado, otimizar apenas as localizações dos suportes (Caso A) reduziu significativamente a tensão, confirmando que o posicionamento do suporte é um fator dominante na mecânica estrutural.
• Limitações de Proxies: O Caso D demonstrou que minimizar as forças de Lorentz (um proxy comum) não necessariamente correlaciona com a minimização da tensão de von Mises, justificando a necessidade de FEA estrutural integrada.
• Efeitos Térmicos: Embora as forças de Lorentz sejam o principal estressor, a contração térmica impacta significativamente a tensão máxima (deslocando o pico em ~17%), destacando a necessidade de incluir cargas térmicas no loop de projeto.

Significância e Alegações
O artigo afirma que o coil-fem acelera significativamente a taxa de iteração para o projeto de stellarators ao eliminar o processo manual de tentativa e erro de posicionamento de estruturas de suporte. Ao permitir a redução simultânea de erros de campo magnético e tensões mecânicas, a ferramenta oferece um caminho para projetos de reatores de fusão mais econômicos e robustos. Os autores posicionam este trabalho como uma "prova de conceito" que estabelece a viabilidade da mecânica estrutural diferenciável no projeto de bobinas de fusão. Eles declaram explicitamente que, embora o modelo atual utilize suposições simplificadoras (ex: materiais isotrópicos, suportes do tipo gaiola modelados como molas), o framework permite extensões futuras para modelos termoelásticos de maior fidelidade e estruturas de suporte complexas (ex: redes de vigas ou cascas). O trabalho não alega ter resolvido todos os desafios mecânicos de stellarator, mas sim fornecer o primeiro framework computacional para otimizar conjuntamente bobinas e seus suportes.