Resumen Técnico: Esquema de volúmenes finitos para magnetohidrodinámica relativista viscorresistiva de primer orden

Planteamiento del Problema

La magnetohidrodinámica (MHD) relativista es esencial para modelar fluidos conductores en entornos de alta energía y astrofísicos, tales como plasmas de quarks y gluones y flujos de acreción. Un desafío persistente en este campo es la incorporación de coeficientes de transporte disipativos (viscosidad, conducción térmica y resistividad) de una manera que preserve la causalidad y la estabilidad. Las teorías tradicionales de primer orden producen ecuaciones parabólicas que violan la causalidad e inducen inestabilidades a velocidades altas. Aunque las teorías de Müller–Israel–Stewart (MIS) resuelven esto promoviendo las corrientes disipativas a variables dinámicas independientes, este enfoque aumenta significativamente la complejidad y el número de grados de libertad.

Desarrollos recientes en la teoría de Bemfica–Disconzi–Noronha–Kovtun (BDNK) ofrecen una alternativa de primer orden mínima que restaura la causalidad y la estabilidad mediante la inclusión de correcciones de derivadas superiores sin introducir nuevas variables. Sin embargo, extender la teoría BDNK a la MHD resistiva presenta una dificultad específica: mientras que la formulación de "una forma" (one-form) de la MHD (tratando el flujo magnético como una carga conservada) elimina la rigidez asociada con la ley de Ampère, vuelve la dinámica resistiva parabólica. Incorporar las correcciones BDNK en este marco de MHD de una forma acoplado, disipativo y de una forma dentro de un esquema numérico estable y eficiente, ha permanecido como una brecha abierta.

Metodología

Los autores presentan una implementación numérica de MHD relativista viscorresistiva de primer orden basada en el marco BDNK, utilizando el código BHAC. La metodología involucra:

1. Formulación Teórica:

   • Los autores adoptan el límite ultra-relativista (negligenciando la densidad de masa) y el enfoque de una forma, donde las líneas de campo magnético son tratadas como objetos similares a cuerdas conservados.
   • Las relaciones constitutivas para el tensor de energía-impulso ($T^{\mu\nu}$) y el tensor de intensidad de campo dual ($J^{\mu\nu}$) incluyen correcciones disipativas de primer orden para la conductividad térmica, la viscosidad y la resistividad.
   • Para asegurar la causalidad, las correcciones BDNK se aplican mediante una transformación on-shell de orden cero. Esto introduce tiempos de relajación ($\tau_\epsilon, \tau_u, \tau_b$) que convierten las ecuaciones de difusión parabólicas en ecuaciones de tipo telégrafo hiperbólicas.
   • Crucialmente, los autores identifican que tres términos BDNK son insuficientes para el sector acoplado de energía-impulso ultra-relativista. Se requiere un cuarto término, controlado por el coeficiente $\tau_X$, para asegurar que el sistema permanezca hiperbólico y causal, evitando el comportamiento elíptico en el sector del sonido.

2. Esquema Numérico:

   • Las ecuaciones se proyectan sobre hipersuperficies espaciotemporales y se descomponen en variables conservadas ($U_I$) y flujos ($F_I$), junto con variables primitivas ($P_I$) y términos fuente ($S_I$).
   • El esquema emplea una discretización de volúmenes finitos con reconstrucción espacial y pasos de tiempo TVD (Total Variation Diminishing) de segundo orden.
   • Los flujos disipativos se evalúan mediante diferencias centrales de segundo orden de los estados de interfaz reconstruidos, con flujos numéricos de tipo upwind computados mediante un resolvedor de Riemann aproximado.
   • Recuperación de Variables Primitivas: Una característica clave de la implementación es un método algebraico eficiente para recuperar las variables primitivas. Debido a que las correcciones BDNK hacen que las variables conservadas dependan linealmente de las derivadas temporales de las variables primitivas, la recuperación se reduce a una inversión de matriz local ($S_I = \sum (M^{-1})_{IJ} [U'_J - U_J(0)]$). Esto evita la necesidad de resolvedores no lineales iterativos típicamente requeridos en esquemas disipativos.
   • El código utiliza algoritmos de transporte restringido (FCT) para mantener la condición de divergencia nula del campo magnético ($\partial_i J^{it} = 0$) con precisión de máquina.

Contribuciones Clave

• Extensión de BDNK a la MHD Resistiva: El artículo construye un esquema numérico que combina la disipación viscosa y resistiva dentro del marco BDNK para la MHD relativista.
• Identificación de un Cuarto Término Necesario: Los autores demuestran que, si bien los términos mínimos de BDNK pueden crear ecuaciones de tipo telégrafo para variables individuales, un cuarto término ($\tau_X$) es estrictamente necesario para garantizar la hiperbolicidad del sistema acoplado de energía-impulso en el límite ultra-relativista.
• Recuperación Eficiente de Primitivas: El trabajo demuestra que las correcciones BDNK permiten un procedimiento de recuperación de variables primitivas algebraico y no iterativo, mejorando significativamente la eficiencia computacional en comparación con los esquemas disipativos tradicionales.
• Implementación Robusta: El esquema está implementado en el código BHAC y validado a través de una gama de regímenes disipativos, desde casi ideales hasta dominados por la difusión.

Resultados

La implementación fue validada a través de una serie de pruebas de referencia (benchmarks) en una y dos dimensiones:

1. Prueba Analítica: Se utilizó una solución de telégrafo con boost para la evolución del campo magnético para verificar la corrección; los resultados numéricos coincidieron con la solución analítica con convergencia de segundo orden.
2. Tubos de Choque 1D: Las simulaciones de tubos de choque con variaciones en la resistividad ($r_b$) y parámetros de difusión ($D_u, D_\epsilon$) mostraron que el esquema permanece estable incluso con datos iniciales discontinuos. Los resultados reprodujeron correctamente el comportamiento de MHD ideal para baja disipación y la dinámica dominada por difusión para alta disipación.
3. Inestabilidad de Kelvin-Helmholtz (KH): Las simulaciones 2D demostraron que aumentar la difusión de momento ($D_u$) suprime la inestabilidad, mientras que aumentar la resistividad ($r_b$) la potencia al relajar el acoplamiento entre la velocidad del fluido y las líneas de campo magnético. El esquema mantuvo velocidades de frente causales durante toda la evolución.
4. Vórtice de Orszag-Tang: Esta prueba resaltó el papel crítico del cuarto término de BDNK ($\tau_X$). Cuando $\tau_X = 0$, la simulación exhibió un comportamiento anti-difusivo y falló debido a velocidades de frente complejas. Con $\tau_X > 0$, el sistema permaneció estable con velocidades de frente reales y sublumínicas.
5. Hoja de Harris: Se simuló una hoja de Harris doble perturbada, que presenta capas de inversión magnética delgadas y grandes gradientes, utilizando refinamiento de malla adaptativa. El esquema resolvió con éxito la formación y deformación de la corriente, manteniendo una velocidad de frente máxima de $v_{max} \approx 0.999$, confirmando la evolución causal.

Significancia

El artículo sostiene que este trabajo proporciona un puente concreto entre los desarrollos teóricos recientes en hidrodinámica relativista (específicamente la teoría BDNK) y su aplicación en simulaciones realistas de flujos magnetizados. Al demostrar que la MHD relativista disipativa puede implementarse de manera numéricamente robusta, causal y eficiente sin introducir variables dinámicas adicionales, los autores establecen una base para simular fenómenos complejos en entornos como plasmas de quarks y gluones débilmente acoplados, el universo temprano y magnetosferas de magnetares. Se destaca la capacidad de controlar independientemente las escalas resistivas y viscosas (fijando el número de Prandtl magnético) como una capacidad crucial para estudiar fenómenos como los dinamos magnéticos y las llamaradas de magnetares. Los autores también señalan que la estructura algebraica de su recuperación de primitivas ofrece un camino prometedor para extender estas correcciones viscorresistivas causales a simulaciones de MHD de relatividad general.