Decoupled iterative schemes for solving stationary MHD problems

Resumen Técnico: Esquemas Iterativos Desacoplados para Problemas de MHD Estacionarios

Planteamiento del Problema
El artículo aborda la solución numérica del sistema de magnetohidrodinámica (MHD) estacionario e incompresible. Este sistema modela la interacción entre fluidos conductores de electricidad y campos electromagnéticos, regido por el acoplamiento de las ecuaciones de Navier-Stokes y las ecuaciones de Maxwell. La formulación matemática consiste en hallar el campo de velocidad ($u$), el campo magnético ($B$), la presión modificada ($p$) y un multiplicador de Lagrange ($\lambda$) que satisfagan un conjunto de ecuaciones diferenciales parciales acopladas. El desafío principal reside en la fuerte no linealidad y el estrecho acoplamiento entre las variables de fluido y magnéticas, lo que complica la construcción de resolvedores numéricos eficientes y estables.

Metodología
Los autores proponen y analizan dos novedosos esquemas iterativos desacoplados inspirados en la estrategia de División de Viscosidad Incremental (IVS, por sus siglas en inglés), desarrollada originalmente para flujos de Navier-Stokes no estacionarios.

1. Esquema IVS Estándar (Algoritmo 5.1):

   • Estrategia de División: El método divide las ecuaciones de momento e inducción con respecto a sus términos difusivos.
   • Paso 1: Resuelve las variables intermedias ($\tilde{u}, \tilde{B}$) tratando los términos convectivos y de acoplamiento de forma explícita utilizando valores de la iteración anterior, mientras se retienen la presión y el multiplicador de Lagrange del paso previo. Esto resulta en un sistema lineal acoplado para los campos intermedios.
   • Pasos 2 y 3: Corrige los campos intermedios para obtener la velocidad final ($u$) y el campo magnético ($B$) mediante la resolución de dos sistemas de tipo Stokes separados y desacoplados. Estos sistemas imponen las restricciones de divergencia nula ($\nabla \cdot u = 0, \nabla \cdot B = 0$) y actualizan la presión y el multiplicador de Lagrange.
   • Ventaja en Álgebra Lineal: Una característica clave es que las matrices del complemento de Schur para las correcciones de Stokes son independientes del nivel de iteración ($k$). En consecuencia, deben ensamblarse y precondicionarse una sola vez, reduciendo la carga computacional en comparación con los métodos donde la matriz del sistema cambia en cada paso.

2. Esquema IVS de Elsässer (Algoritmo 6.1):

   • Reformulación: Para desacoplar aún más el sistema, los autores reformulan las ecuaciones de MHD utilizando las variables de Elsässer ($w = u + \sqrt{\kappa}B$ y $z = u - \sqrt{\kappa}B$).
   • Desacoplamiento: En esta formulación, el acoplamiento entre las nuevas variables ocurre a través de un término no lineal y uno lineal. El algoritmo propuesto resuelve $w$ y $z$ en pasos separados, reemplazando el gran sistema acoplado del IVS estándar por dos sistemas lineales más pequeños y totalmente desacoplados.
   • Aplicabilidad: Este enfoque es teóricamente ventajoso bajo condiciones de contorno de Dirichlet homogéneas, aunque los autores señalan que surgen desafíos con las condiciones de contorno mixtas.

Contribuciones Teóricas
El artículo proporciona un análisis matemático riguroso para ambos esquemas:

• Bien posedura: Los autores establecen resultados de existencia y unicidad para el problema de MHD continuo bajo un supuesto de "datos pequeños" (específicamente, un límite sobre los números de Reynolds y de Reynolds magnético en relación con la viscosidad).
• Pruebas de Convergencia:
  • Se analiza el esquema iterativo de Picard clásico y se demuestra su convergencia bajo el mismo supuesto de datos pequeños.
  • Ambos esquemas, IVS estándar e IVS de Elsässer, se demuestran uniformemente acotados y con convergencia geométrica hacia la solución única.
  • Las demostraciones utilizan estimaciones de energía, la condición inf-sup (LBB) y lemas específicos sobre la convergencia de secuencias no negativas.
• Análisis de Error: Se derivan cotas teóricas que muestran que los errores en la velocidad, el campo magnético, la presión y los multiplicadores de Lagrange decaen a cero a medida que aumenta el número de iteraciones.

Resultos Numéricos
Los autores validan sus hallazgos teóricos utilizando el software FreeFem++ con pares de elementos finitos $P_2-P_1$.

• Pruebas de Convergencia: Soluciones manufacturadas en 2D demuestran tasas de convergencia óptimas ($O(h^2)$ en $L^2$ y $O(h)$ en $H^1$) tanto para el campo de velocidad como para el magnético, confirmando la precisión de los esquemas.
• Problemas de Referencia (Benchmarks):
  • Flujo de Hartmann: Se resolvieron los problemas de flujo de Hartmann en 2D y 3D. Los resultados numéricos coincidieron con las soluciones analíticas con alta precisión para diversos números de Hartmann ($Ha$).
  • Cavidad Impulsada por Tapa (Lid-Driven Cavity): Los esquemas simularon con éxito patrones de flujo complejos para varios números de Reynolds ($Re$), de Reynolds magnético ($Rm$) y de acoplamiento ($\kappa$), incluyendo casos de alto acoplamiento ($\kappa=5000$).
  • Flujo sobre un Escalón: El método capturó con precisión la separación del flujo y las caídas de presión, demostrando robustez en el manejo de singularidades geométricas y variaciones en la fuerza de acoplamiento.

Significancia y Reivindicaciones
El artículo afirma que los esquemas iterativos desacoplados propuestos ofrecen una alternativa práctica y eficiente a los métodos iterativos clásicos monolíticos o totalmente acoplados (como los métodos de Newton o Picard estándar).

• Eficiencia: Al dividir los términos difusivos, el método reduce la complejidad de los sistemas lineales resueltos en cada paso. Específicamente, la formulación de Elsässer permite la resolución de dos sistemas más pequeños y desacoplados en lugar de un único sistema grande y acoplado.
• Reutilización: El enfoque IVS asegura que las matrices del complemento de Schur para las correcciones de Stokes sean estáticas, lo que permite un único paso de ensamblaje y precondicionamiento, siendo esto computacionalmente beneficioso para problemas en 3D.
• Robustez: Los esquemas mantienen la propiedad de divergencia nula de los campos de velocidad y magnético a nivel discreto y permanecen estables bajo las condiciones de datos pequeños establecidas.

Los autores concluyen que, si bien el presente trabajo se centra en problemas estacionarios, la metodología proporciona una base para futuras extensiones a problemas de MHD dependientes del tiempo (no estacionarios), con un enfoque en el análisis de estabilidad temporal y espacial.