Error analysis of a divergence-preserving mixed finite element scheme for the incompressible Hall--magnetohydrodynamic equations

Resumen Técnico: Análisis de Error de un Esquema de Elementos Finitos Mixtos que Preserva la Divergencia para las Ecuaciones Hall–MHD Incompresibles

Planteamiento del Problema
El sistema de magnetohidrodinámica con efecto Hall (Hall–MHD) incompresible modela la interacción entre el movimiento de fluidos y los campos magnéticos en fluidos eléctricamente conductores, incorporando el efecto Hall, que se vuelve significativo a escalas comparables a la profundidad de piel iónica. Este efecto introduce un término rígido y altamente no lineal ($\eta \text{curl}(\text{curl} B \times B)$) en la ecuación de inducción, lo que conduce a un comportamiento dispersivo, gradientes pronunciados y oscilaciones de alta frecuencia (ondas silbantes). Estas características plantean desafíos analíticos y computacionales sustanciales. Específicamente, la existencia global de soluciones fuertes en tres dimensiones sigue siendo un problema abierto, y los métodos numéricos estándar a menudo luchan con la estabilidad y el requisito estricto de que el campo magnético permanezca solenoidal ($\text{div } B = 0$). Las violaciones de esta condición libre de divergencia pueden inducir monopolos magnéticos no físicos y desencadenar inestabilidades numéricas severas. Además, los enfoques numéricos existentes, como los basados en métodos de Variable Auxiliar Escalar (SAV) o técnicas espectrales, a menudo disipan funcionales de energía modificados en lugar de la energía física real, no logran imponer la condición libre de divergencia exactamente a nivel discreto, o carecen de un análisis de error riguroso y flexibilidad geométrica.

Metodología
Para abordar estos desafíos, los autores proponen un método de elementos finitos mixtos (FEM) totalmente discreto, lineal y que preserva la estructura para un sistema Hall–MHD regularizado mediante Voigt. La regularización de Voigt reemplaza la derivada temporal en la ecuación de momento con una versión filtrada espacialmente $(I - \alpha_1 \Delta)\partial_t u$ y retiene un término de inercia electrónica $\alpha_2 \partial_t J$ en la ley de Ohm generalizada. Esta modificación mejora la buena posición y la estabilidad numérica, preservando al mismo tiempo la estructura física esencial y las leyes de energía del sistema original.

El esquema numérico se construye utilizando los siguientes componentes:

1. Formulación Mixta: El problema se reformula utilizando el campo eléctrico $E$ y la densidad de corriente $J$ como variables primarias, lo que conduce a un sistema acoplado de ecuaciones de momento, inducción, ley de Ohm y ley de Ampère.
2. Discretización Espacial: El esquema emplea el Cálculo de Elementos Finitos Exterior (FEEC) para garantizar que la condición libre de divergencia se satisfaga exactamente.
   • La velocidad $u$ y la presión $p$ se aproximan utilizando el par de elementos finitos MINI (inf-sup estable de orden más bajo).
   • El campo magnético $B$ y la densidad de corriente $J$ se aproximan utilizando los espacios de Raviart–Thomas ($RT_0$) y Nédélec ($X_0$) de orden más bajo, respectivamente.
   • Estos espacios forman un complejo de de Rham discreto, asegurando que $\text{div } B_h = 0$ se preserve exactamente a nivel discreto.
3. Discretización Temporal: Se utiliza un esquema de paso de tiempo linealmente implícito y antisimétrico. Los términos convectivos y de Lorentz se tratan semi-implícitamente utilizando valores del paso de tiempo anterior para mantener la linealidad y la estabilidad energética incondicional.
4. Análisis de Error: Los autores establecen tasas de convergencia óptimas descomponiendo el error en errores de aproximación (utilizando proyecciones de Stokes y proyecciones ortogonales) y errores discretos. Utilizan desigualdades de Gaffney discretas y estimaciones de energía para manejar los términos de acoplamiento no lineales.

Contribuciones Clave

• Esquema que Preserva la Estructura: El artículo presenta un FEM totalmente discreto y lineal que impone la restricción solenoidal sobre el campo magnético exactamente, evitando la necesidad de pasos de proyección o métodos de penalización que pueden degradar la precisión o la estabilidad.
• Estabilidad Energética Incondicional: El esquema propuesto satisface una ley de disipación de energía discreta incondicionalmente, reflejando la ley de energía física continua, incluso en presencia del término Hall rígido.
• Estimaciones de Error Rigurosas: Los autores proporcionan un análisis de error completo para el problema regularizado por Voigt, estableciendo tasas de convergencia óptimas del orden $O(h^s + \tau)$ bajo suposiciones de regularidad adecuadas.
• Extensión al Sistema No Regularizado: Crucialmente, el análisis se extiende al sistema Hall–MHD no regularizado (donde $\alpha_1 = \alpha_2 = 0$). Los autores demuestran que la convergencia aún puede probarse para el sistema no regularizado, siempre que la resistividad $\sigma > 0$. En el régimen no resistivo ($\sigma = 0$), se muestra que la regularización de Voigt juega un papel crucial en las estimaciones de error.
• Implementación 2.5D y 3D: El marco se valida tanto en configuraciones 2.5D (invarianza traslacional) como totalmente 3D, demostrando la aplicabilidad del método a geometrías complejas y fenómenos físicos.

Resultados
Las simulaciones numéricas en 2.5D y 3D confirman los hallazgos teóricos:

• Convergencia: El esquema logra las tasas de convergencia espacial y temporal predichas para la velocidad, el campo magnético y la densidad de corriente. Si bien las tasas de convergencia para la velocidad y la densidad de corriente se reducen ligeramente en algunos casos (atribuido a la regularidad limitada de la solución debido al término Hall), el campo magnético mantiene una convergencia óptima.
• Preservación de la Divergencia: El error máximo de divergencia ($\max | \text{div } B_h |$) permanece al nivel de la tolerancia del solucionador durante toda la simulación, confirmando la preservación exacta de la restricción solenoidal.
• Estabilidad Energética: La energía total del sistema decae monótonamente con el tiempo, en consonancia con la ley teórica de disipación de energía.
• Fidelidad Física: Las simulaciones del flujo de Arnold–Beltrami–Childress (ABC), el vórtice de Orszag–Tang y la reconexión impulsada por Hall en una hoja de Harris capturan con éxito fenómenos físicos complejos, incluida la reconexión magnética, la formación de islas magnéticas, cuerdas de flujo y los patrones de corriente cuadrupolares característicos asociados con la dinámica Hall–MHD. El esquema resuelve estas características con precisión incluso en mallas relativamente gruesas.

Significado
El artículo afirma que este trabajo proporciona un marco matemáticamente riguroso y computacionalmente robusto para la simulación de sistemas Hall–MHD. Al combinar la regularización de Voigt con una formulación de elementos finitos mixtos que preserva la divergencia, los autores superan los desafíos duales de la no linealidad del término Hall y el requisito estricto de divergencia nula. El método ofrece una aproximación físicamente consistente que es estable, convergente y capaz de manejar geometrías complejas, llenando un vacío dejado por métodos anteriores que carecían de un análisis de error riguroso o una preservación exacta de la divergencia. Los resultados sugieren que este enfoque es una herramienta viable para estudiar fenómenos de plasma astrofísico y de laboratorio donde los efectos Hall son dominantes.