Towards joint optimization of stellarator coils and support structures

Resumen Técnico: Hacia la optimización conjunta de bobinas y estructuras de soporte de esteleradores

Planteamiento del problema
El diseño de reactores de fusión tipo estelerador enfrenta un desafío crítico: equilibrar la necesidad de un campo magnético preciso para el confinamiento del plasma con los requisitos mecánicos para proteger los conductores de superconductores de alta temperatura (HTS) —que son frágiles— de las enormes fuerzas de Lorentz y los esfuerzos térmicos. Actualmente, el proceso de diseño es secuencial y requiere mucha mano de obra. Las geometrías de las bobinas se optimizan para el rendimiento magnético de forma aislada, y posteriormente se añaden las estructuras de soporte de manera manual mediante análisis de elementos finitos (FEA) iterativos para cumplir con los límites de tensión de ingeniería. Este enfoque de "posprocesamiento" no tiene en cuenta cómo la ubicación de los soportes influye en el esfuerzo de la bobina y cómo la deformación de la bobina bajo carga introduce errores en el campo magnético. Las herramientas de optimización existentes dependen de aproximaciones analíticas simplificadas que ignoran las estructuras de soporte o utilizan herramientas de FEA comerciales de código cerrado que no pueden interactuar con los bucles de optimización de bobinas basados en gradientes.

Metodología
Los autores presentan coil-fem, una herramienta de mecánica estructural diferenciable y de código abierto integrada en el bucle de optimización de bobinas de esteleradores. Construida sobre la biblioteca JAX-FEM, coil-fem permite la optimización simultánea de la geometría de la bobina y la ubicación de las estructuras de soporte mediante métodos basados en gradientes.

Los componentes metodológicos clave incluyen:

• FEA Diferenciable: La herramienta resuelve la ecuación de elasticidad lineal $\nabla \cdot \sigma(u) + F_{body} = 0$ con condiciones de contorno de Robin (fundación de resorte). Esta aproximación reemplaza las condiciones de Dirichlet fijas, permitiendo que las ubicaciones de los cierres de soporte sean tratadas como variables optimizables sin alterar la estructura de la matriz dispersa del sistema FEA.
• Modelo de Física: El modelo contabiliza las fuerzas corporales, incluyendo la gravedad, las fuerzas magnéticas mutuas (modeladas como filamentos 1D) y las fuerzas propias (calculadas mediante la fórmula de Landreman-Hurwitz para secciones transversales rectangulares). También incorpora la deformación térmica resultante de la caída de temperatura de $\sim$200 K desde el estado libre de tensiones hasta la operación.
• Mallado y Geometría: Las bobinas se modelan como cuerpos de construcción finita generados mediante el barrido de una sección transversal rectangular a lo largo de una línea central utilizando un marco de rotación mínima. La topología de la malla se mantiene fija durante la optimización para garantizar la diferenciabilidad y la eficiencia computacional.
• Bucle de Optimización: Los autores formulan una función de penalización multiobjetivo que incluye el error del campo magnético, la carga mecánica (integral de volumen del cuadrado del esfuerzo de von Mises) y restricciones geométricas estándar (longitud de la bobina, espaciado, curvatura). Se utiliza el algoritmo L-BFGS para minimizar este objetivo con respecto tanto a los coeficientes de Fourier de la bobina como a las ubicaciones paramétricas de los cierres de soporte.

Contribuciones Clave

1. Primera Optimización Conjunta: Este artículo presenta el primer caso en la literatura de una optimación conjunta, impulsada por gradientes, de la geometría de la bobina del estelerador y la ubicación de la estructura de soporte.
2. Herramienta coil-fem: El lanzamiento de una herramienta de código abierto que integra la FEA de soporte diferenciable directamente en el bucle de optimización de la bobina, tratando las ubicaciones de los cierres de soporte como variables optimizables.
3. Demostración de la Importancia del Soporte: El estudio proporciona evidencia numérica de que la colocación de las estructuras de soporte es un factor crítico para reducir el esfuerzo mecánico, un factor que a menudo se pasa por alto cuando los soportes son fijos o se colocan manualmente.

Resultados
Los autores compararon coil-fem con la suite FEA DOLFINx y el enfoque de sustitución de estilo comercial, demostrando una concordancia dentro del 1–5% para las métricas de esfuerzo y desplazamiento. Luego aplicaron la herramienta a un conjunto de bobinas de W7-X simplificado bajo cuatro escenarios:

• Línea base (Baseline): Bobinas originales de W7-X con cierres superior e inferior fijos.
• Caso A (Soportes Optimizados): Geometría de bobina fija, ubicaciones de cierres optimizadas.
• Caso B (Bobinas Optimizadas): Geometría de bobina optimizada, ubicaciones de cierres fijas.
• Caso C (Optimización Conjunta): Tanto la geometría de la bobina como la ubicación de los cierres optimizadas.
• Caso D (Minimización de Fuerza): Geometría de bobina optimizada para minimizar la fuerza de Lorentz (proxy estándar) con cierres fijos.

Hallazgos Clave:

• Reducción de Esfuerzo: El Caso C logró una reducción de 2.4× en el esfuerzo de von Mises RMS en comparación con la línea base.
• Rol de los Soportes: Optimizar únicamente la geometría de la bobina con soportes fijos (Caso B) produjo solo mejoras marginales en el esfuerzo. Por el contrario, optimizar únicamente las ubicaciones de los soportes (Caso A) redujo significativamente el esfuerzo, confirmando que la colocación del soporte es un factor dominante en la mecánica estructural.
• Limitaciones de los Proxies: El Caso D demostró que minimizar las fuerzas de Lorentz (un proxy común) no necesariamente se correlaciona con la minimización del esfuerzo de von Mises, lo que justifica la necesidad de una FEA estructural integrada.
• Efectos Térmicos: Aunque las fuerzas de Lorentz son el principal estresor, la contracción térmica impacta significativamente el esfuerzo máximo (desplazando el pico en un ~17%), lo que resalta la necesidad de incluir cargas térmicas en el bucle de diseño.

Significancia y Reivindicaciones
El artículo afirma que coil-fem acelera significativamente la tasa de iteración para el diseño de esteleradores al eliminar el proceso manual de ensayo y error para la colocación de las estructuras de soporte. Al permitir la reducción simultánea de los errores del campo magnético y los esfuerzos mecánicos, la herramienta ofrece una vía hacia diseños de reactores de fusión más económicos y robustos. Los autores posicionan este trabajo como una "prueba de concepto" que establece la viabilidad de la mecánica estructural diferenciable en el diseño de bobinas de fusión. Expresan explícitamente que, si bien el modelo actual utiliza suposiciones simplificadoras (por ejemplo, materiales isotrópicos, soportes tipo jaula modelados como resortes), el marco permite futuras extensiones a modelos termoelásticos de mayor fidelidad y estructuras de soporte complejas (por ejemplo, redes de vigas o carcasas). El trabajo no pretende haber resuelto todos los desafíos mecánicos de los esteleradores, sino que proporciona el primer marco computacional para optimizar conjuntamente las bobinas y sus soportes.