Improving divergence cleaning in cosmological SPMHD simulations

Resumen Técnico: Mejora de la Limpieza de la Divergencia en Simulaciones de SPMHD Cosmológicas

Planteamiento del Problema
En las simulaciones de magnetohidrodinámica por partículas suavizadas (SPMHD) cosmológicas, mantener la condición solenoidal ($\nabla \cdot \mathbf{B} = 0$) es un desafío numérico persistente. Los esquemas estándar en códigos como Gadget suelen depender del enfoque de "8 ondas" de Powell et al. (1999). Aunque este método advecta los errores de divergencia con el flujo del fluido, no los elimina activamente. En los códigos basados en SPH, los errores residuales de $\nabla \cdot \mathbf{B}$ pueden alcanzar órdenes de magnitud unitarios, lo que conduce a una amplificación no física del campo magnético y a la "inestabilidad de tensión", donde términos fuente espurios en la ecuación de cantidad de movimiento dominan cuando la presión magnética supera la presión del gas. Si bien el esquema de limpieza hiperbólico/parabólico propuesto por Dedner et al. (2002) ofrece una solución, su formulación original ha demostrado ser inestable en SPMHD, particularmente en entornos dinámicos de pasos de tiempo individuales con saltos de densidad fuertes, actuando a menudo como una fuente no física de energía magnética.

Metodología
Los autores implementan una versión restringida y de conservación de energía del algoritmo de limpieza tipo Dedner (Tricco & Price 2012; Tricco et al. 2016a) en el código de SPMHD cosmológico OpenGadget3. Las adaptaciones metodológicas clave incluyen:

• Conservación de la Energía: Las ecuaciones gobernantes se reformulan para garantizar la conservación de la energía total. Esto implica evolucionar la variable de campo de limpieza como $\psi/c_h$ (donde $\psi$ es el escalar auxiliar y $c_h$ es la velocidad de limpieza) en lugar de solo $\psi$, e incluir un término fuente específico $-\frac{1}{2}\psi(\nabla \cdot \mathbf{v})$ en la ecuación de evolución. Esto asegura que el amortiguamiento parabólico del campo de limpieza elimine la energía magnética sin introducir un crecimiento ilimitado.
• Adaptación Cosmológica: Las ecuaciones se adaptan para coordenadas comóviles y pasos de tiempo de partículas individuales. La velocidad de limpieza $c_h$ se establece como la velocidad magnetosónica rápida local ($\sqrt{v_A^2 + c_s^2}$), asegurando que la limpieza sea intrínsecamente estable en el sentido de Courant dentro de las restricciones de paso de tiempo estándar de SPH.
• Configuración de la Simulación: Los autores realizan una serie de simulaciones sobre un cúmulo de galaxias masivo ($M_{200c} \approx 10^{15} M_\odot$) sembrado con un campo magnético uniforme ($10^{-14}$ G comóvil). Se comparan cinco ejecuciones distintas:
  1. NIFTY-noclean: Ejecución de referencia que utiliza únicamente el esquema de advección de Powell.
  2. NIFTY-clean: Ejecución fiduciaria con el esquema de limpieza restringido completo (incluyendo la retroacción de la entropía).
  3. NIFTY-clean2: Igual que la anterior pero con un factor de sobre-limpieza $\zeta=2$ (limpieza a dos veces la velocidad de onda local).
  4. NIFTY-advect: Esquema de limpieza con el término fuente de divergencia de velocidad desactivado.
  5. NIFTY-heat: Esquema de limpieza donde la energía de limpieza disipada se descarta en lugar de retroalimentarse a la ecuación de entropía.

Resultos Clave

• Reducción del Error de Divergencia: El resultado principal es una reducción en los errores de divergencia relativa ($|\nabla \cdot \mathbf{B}|/B$) de 2 a 3 órdenes de magnitud en todo el volumen del cúmulo en comparación con el esquema de solo Powell. Esta reducción es uniforme tanto en el núcleo como en las periferias del cúmulo.
• Amplificación del Campo Magnético: En las periferias del cúmulo ($r \approx 1\text{--}3 \, h^{-1}$ Mpc), donde el campo magnético es más débil y la resolución es menor, el esquema de limpieza activo permite una amplificación del campo magnético de un factor de 5 a 10 en comparación con la ejecución de solo Powell. Los autores atribuyen esto a la eliminación de las fuerzas de Lorentz no físicas causadas por los errores de divergencia, que anteriormente suprimían el dínamo de pequeña escala en estas regiones.
• Estabilidad del Núcleo: La estructura de densidad global y la geometría del campo magnético en el núcleo del cúmulo permanecen mayormente inalteradas entre las ejecuciones con y sin limpieza, lo que indica que la limpieza no altera artificialmente el estado saturado del dínamo del núcleo.
• Consistencia Termodinámica: Se demuestra que la inclusión de la retroacción de la entropía (retroalimentar la energía de limpieza disipada como calor) es crítica para la consistencia termodinámica. Las ejecuciones sin esta retroacción (NIFTY-heat) resultan en una entropía central más baja y campos magnéticos centrales saturados más altos debido a la compresión adiabática de un núcleo más frío.
• Costo Computacional: La implementación no introduce ningún sobrecosto medible en el tiempo de ejecución; la mejora en la limpieza se logra esencialmente "gratis" dentro de las restricciones de paso de tiempo de Courant existentes.

Significado y Reivindicaciones
El artículo sostiene que el control preciso de los errores de divergencia es esencial para capturar la amplificación del campo magnético en las regiones de baja densidad y pobremente resueltas de los cúmulos de galaxias. Los autores argumentan que los errores de divergencia incontrolados en los esquemas estándar de Powell actúan como ruido numérico que suprime el dínamo en las regiones sub-saturadas. Al implementar el esquema de limpieza restringido y de conservación de energía, la simulación logra una realización más fiel de la física de MHD sin penalización computacional.

Los autores concluyen que el esquema de limpieza restringido:

1. Reduce drásticamente los errores numéricos de divergencia (2 a 3 órdenes de magnitud).
2. Permite una representación más precisa de la amplificación del campo en regiones sub-resueltas.
3. Produce estructuras de campo magnético más suaves y coherentes.
4. Estabiliza la evolución del campo magnético a través de diferentes resoluciones, mejorando la fidelidad física de los resultados incluso con resoluciones moderadas.

El trabajo recomienda que las aplicaciones de SPMHD cosmológica adopten este esquema de limpieza por defecto, utilizando la velocidad de limpieza natural ($\zeta=1$) y manteniendo activa la retroacción de energía parabólica para cerrar el balance de energía. La implementación ya está disponible públicamente en OpenGadget3.