Systematic comparison of VMEC and HINT equilibrium calculations for finite-beta LHD plasmas

Este artículo compara sistemáticamente los cálculos de equilibrio de VMEC y HINT para plasmas del Dispositivo Helicoidal Grande, revelando que, si bien ambos códigos coinciden a valores bajos de beta, divergen a valores de beta más altos, ya que HINT captura la estocasticidad del borde y la ruptura de las superficies de flujo que el supuesto de superficies de flujo anidadas de VMEC no puede representar.

Lo esencial

Imagina intentar hornear un pastel perfectamente redondo dentro de un horno muy extraño y retorcido. En el mundo de la energía de fusión, los científicos utilizan máquinas llamadas estelaradores (como el Dispositivo Helicoidal Grande, o LHD) para atrapar plasma supercaliente. Para mantener este plasma estable, necesitan calcular exactamente cómo deben verse las "paredes" magnéticas que lo contienen.

Este artículo compara dos "reposteros" diferentes (programas informáticos) que intentan determinar la forma de estas paredes magnéticas cuando el plasma se vuelve muy caliente y presurizado.

Los dos reposteros: VMEC y HINT

1. VMEC (El Arquitecto Estricto): Este programa es como un arquitecto que insiste en que cada capa del pastel debe ser una cebolla perfecta, suave y anidada. Supone que las paredes magnéticas nunca se rompen ni se tocan entre sí. Es excelente para situaciones simples y de baja presión, pero tiene un punto ciego: se niega a creer que las paredes puedan volverse desordenadas o romperse.
2. HINT (El Observador Realista): Este programa es como un científico que observa cómo se hornea realmente el pastel. No asume que las capas sean perfectas. En su lugar, deja que la física suceda de forma natural. Si el calor aumenta demasiado, permite que las paredes magnéticas se vuelvan inestables, se rompan o se conviertan en un caos.

El experimento: Subiendo la temperatura

Los investigadores probaron estos dos programas en la máquina LHD con tres formas diferentes de "horno" magnético (algunas desplazadas hacia adentro, otras hacia afuera). Aumentaron lentamente la presión del plasma (el "calor" del pastel) del 0% al 5%.

¿Qué pasó con la presión baja?
Cuando el plasma estaba frío y tranquilo, ambos reposteros estaban de acuerdo. Las paredes magnéticas se mantenían suaves y anidadas, tal como predijo el Arquitecto Estricto (VMEC). Todo estaba bien.

¿Qué pasó cuando subió el calor?
Una vez que la presión cruzó cierto "punto crítico", los dos reposteros empezaron a discrepar.

• VMEC seguía dibujando capas de cebolla perfectas, suaves y en expansión. Pensaba que el plasma simplemente se estaba haciendo más grande y redondo.
• HINT vio algo diferente. Notó que las paredes magnéticas empezaban a volverse "estocásticas".

El desastre "estocástico": Una analogía creativa

Piensa en las líneas del campo magnético como un manojo de espaguetis.

• En un estado perfecto (presión baja), las hebras de espagueti están ordenadamente agrupadas y corren paralelas entre sí.
• A medida que la presión aumenta, la corriente Pfirsch-Schlüter (un tipo de corriente eléctrica que se forma naturalmente en el plasma) actúa como una mano caótica que mezcla los espaguetis.
• Eventualmente, las hebras comienzan a superponerse y enredarse. Esto se llama islas magnéticas y estocasticidad. Las ordenadas "capas de cebolla" se rompen.

Debido a que HINT permite este enredo, ve cómo la "jaula magnética" se encoge. La mezcla caótica en el borde del plasma hace que el volumen efectivo sea menor. VMEC, sin embargo, sigue dibujando la cebolla perfecta y en expansión, por lo que piensa que el volumen está aumentando.

Los hallazgos clave

1. El "punto de inflexión": Existe un nivel de presión específico donde las ordenadas capas de cebolla se rompen. Una vez que pasas este punto, VMEC ya no es preciso porque no puede ver las paredes rotas.
2. La forma importa: El "punto de inflexión" ocurre antes (a una presión más baja) si la máquina está desplazada hacia afuera.
   • Analogía: Imagina que la máquina desplazada hacia afuera es como una mesa tambaleante. Es más fácil volcarla (crear el caos) que una mesa desplazada hacia adentro que es robusta. La forma hacia afuera crea más "ondulaciones" en el campo magnético, haciendo que el espagueti se enrede más rápido.
3. Pérdida de volumen: En las configuraciones desplazadas hacia afuera y estándar, a medida que la presión se vuelve muy alta, el volumen real del plasma (según el modelo realista de HINT) comienza a encogerse debido a que las paredes magnéticas se rompen. VMEC pasa esto por alto por completo, pensando que el volumen sigue creciendo.

La conclusión fundamental

Este artículo demuestra que para plasmas de fusión de alta presión, no podemos confiar únicamente en el modelo de la "cebolla perfecta" (VMEC). Necesitamos al "observador realista" (HINT) para ver cuándo las paredes magnéticas se están rompiendo y volviendo caóticas. Esto es especialmente cierto para las máquinas desplazadas hacia afuera, donde el campo magnético es más sensible a estos efectos tridimensionales desordenados. El estudio confirma que, a medida que presionamos para obtener más energía, la suposición de capas magnéticas perfectas y suaves se vuelve cada vez menos válida.

Resumen técnico

Resumen Técnico: Comparación Sistemática de Cálculos de Equilibrio VMEC y HINT para Plasmas LHD de $\beta$ Finito

Planteamiento del Problema
En la física de estelaradores, el cálculo del equilibrio magnetohidrodinámico (MHD) tridimensional se realiza frecuentemente utilizando el solver VMEC. VMEC se basa en el principio variacional y asume estrictamente la existencia de superficies de flujo perfectamente anidadas. Sin embargo, en plasmas de $\beta$ finito, el aumento de la presión del plasma puede alterar la topología del campo magnético, generando islas magnéticas y regiones estocásticas. Estos fenómenos, particularmente en dispositivos de corriente neta cero donde son impulsados por corrientes de Pfirsch-Schlüter (PS), pueden comprometer la validez del supuesto de superficies de flujo anidadas. Cuando las islas magnéticas se solapan, la última superficie de flujo cerrado (LCFS) se rompe, reduciendo el volumen de plasma encerrado —una limitación que se observa más severa que el desplazamiento de Shafranov. Mientras que el código HINT resuelve ecuaciones dinámicas para campos magnéticos y presión sin asumir superficies anidadas, ha sido necesaria una comparación sistemática entre VMEC y HINT a través de diversas configuraciones y valores de beta en el Dispositivo Helicoidal Grande (LHD) para aclarar la influencia de estos supuestos.

Metodología
El estudio llevó a cabo una comparación sistemática de los cálculos de equilibrio entre VMEC y HINT para plasmas LHD. La investigación se centró en tres configuraciones específicas de eje magnético de vacío: $R_{ax}^V = 3.53$ m (desplazado hacia adentro), $3.60$ m (estándar) y $3.85$ m (desplazado hacia afuera). El análisis cubrió valores de beta en el eje ($\beta_0$) que oscilaron desde el $0.0\%$ al $5.0\%$ en incrementos de $0.5\%$.

Ambos códigos utilizaron un perfil de presión inicial definido como $p = p_0(1-s)$, donde $s$ es el flujo toroidal normalizado. El estudio comparó tres métricas primarias:

1. La posición del eje magnético ($\langle R_{ax} \rangle$).
2. El transform de rotación en el eje ($\iota/2\pi$).
3. El volumen encerrado por la LCFS clara.

El código HINT, que emplea un método de relajación para resolver islas magnéticas y regiones estocásticas, se contrastó con VMEC, el cual impone la restricción de superficie de flujo anidada.

Resultos Clave

• Régimen de Bajo $\beta$: Para valores bajos de $\beta_0$, VMEC y HINT produjeron soluciones de equilibrio consistentes. Esto indica que el supuesto de superficies de flujo anidadas sigue siendo ampliamente válido en este régimen, y las respuestas MHD 3D son comparables entre los dos códigos.
• $\beta_0$ Crítico y Divergencia: Se identificó un $\beta_0$ crítico dependiente de la configuración. Más allá de este umbral, las soluciones de los dos códigos comenzaron a divergir. En los cálculos de HINT, el eje magnético exhibió un desplazamiento de Shafrander mayor y un transform de rotación más alto en comparación con VMEC.
• Ruptura de Superficie de Flujo y Volumen: En HINT, el volumen encerrado por la LCFS clara disminuyó en valores altos de beta para las tres configuraciones. Esta reducción se atribuye a la evolución de los campos magnéticos estocásticos cerca del borde del plasma, lo que rompe las superficies de flujo. Por el contrario, VMEC, restringido por su supuesto de superficie anidada, mostró un aumento monotónico en el volumen encerrado a medida que $\beta_0$ aumentaba, fallando en representar la ruptura de la superficie de flujo.
• Dependencia de la Configuración: El valor de $\beta_0$ crítico disminuyó desde las configuraciones desplazadas hacia adentro hasta las desplazadas hacia afuera. Se encontró que la configuración desplazada hacia afuera ($R_{ax}^V = 3.85$ m) era la más sensible a las respuestas de equilibrio 3D. Esta sensibilidad surge porque el rizado helicoidal aumenta con $R_{ax}^V$, impulsando corrientes PS más fuertes, mientras que el cizallamiento magnético disminuye. La combinación de corrientes PS más fuertes y un menor cizallamiento conduce a islas magnéticas más grandes y una estocasticidad más temprana a valores de beta más bajos.

Significancia y Reivindicaciones
El artículo establece que las respuestas de equilibrio 3D en los plasmas LHD se vuelven cada vez más significativas a medida que la configuración se desplaza de hacia adentro hacia hacia afuera. El estudio demuestra que el supuesto de superficie de flujo anidada, central en VMEC, se ve comprometido a beta finito debido a las perturbaciones impulsadas por la corriente PS y la resultante estocasticidad en el borde.

Los autores afirman que, si bien VMEC y HINT son comparables para configuraciones de bajo $\beta_0$, HINT proporciona una descripción físicamente más completa en beta alto al resolver la ruptura de las superficies de flujo y la reducción asociada del volumen de plasma. La divergencia entre los códigos resalta que la limitación del supuesto de superficie de flujo anidada no es uniforme, sino que depende fuertemente de la configuración magnética de vacío, específicamente de la interacción entre la toroidicidad (que impulsa las corrientes PS) y el cizallamiento magnético. Los resultados subrayan la importancia de utilizar códigos capaces de resolver la estocasticidad, como HINT, al analizar equilibrios de alto beta en estelaradores, particularmente en configuraciones desplazadas hacia afuera donde el umbral de beta crítico es menor.