An Adjoint Projection Formulation for Enforcing the divergence-free Constraint in Smoothed Particle Magnetohydrodynamics

Resumen Técnico: Una Formulación de Proyección Adjunta para Imponer la Restricción $\nabla \cdot \mathbf{B} = 0$ en la Magnetohidrodinámica de Partículas Suavizadas

Planteamiento del Problema
En la magnetohidrodinámica (MHD), la restricción de solenoidal $\nabla \cdot \mathbf{B} = 0$ es un requisito físico fundamental. Mientras que los códigos basados en mallas suelen emplear el transporte restringido (CT) para preservar esta restricción de forma discreta, los métodos sin malla como la Magnetohidrodinámica de Partículas Suavizadas (SPMHD) carecen de una topología natural de cara-arista para aplicar el CT directamente. En consecuencia, el enfoque principal en SPMHD ha sido la "limpieza de divergencia", que evoluciona un campo escalar auxiliar para propagar y amortiguar los errores de divergencia. Sin embargo, los métodos de limpieza introducen parámetros físicos (velocidad de propagación y escala de tiempo de amortiguación) que pueden afectar los resultados de la simulación si no se eligen de manera consistente, particularmente en sistemas con grandes rangos dinámicos. Los intentos de proyección previos en SPMHD, como los de Price y Monaghan, utilizaron sumatorias de tipo gravitatorio para las correcciones, pero sufrieron de inconsistencias entre el operador de divergencia utilizado para medir los errores y el operador de gradiente utilizado para computar las correcciones, lo que resultó en un desempeño deficiente ante el ruido de longitud de onda corta.

Metodología
El artículo propone un método de proyección que impone la restricción discreta libre de divergencia resolviendo un problema elíptico después de cada actualización de MHD. La innovación central es la construcción de un operador de corrección que es estrictamente el adjunto del operador de divergencia discreto utilizado para medir el error, definido dentro de una métrica ponderada por volumen.

1. Formulación: El método busca un campo magnético corregido $\mathbf{B}^{n+1}$ que minimice el cambio en la norma de la energía magnética discreta (definida por los volúmenes de las partículas $V_i = m_i/\rho_i$) sujeto a la restricción $D\mathbf{B}^{n+1} = 0$, donde $D$ es el operador de divergencia discreto de SPH.
2. Gradiente Adjunto: Mediante la introducción de un campo de multiplicador de Lagrange $\pi$, la corrección se deriva como $\mathbf{B}^{n+1} = \mathbf{B}^* - \mathcal{G}\pi$. El operador $\mathcal{G}$ se construye explícitamente como el adjunto de la métrica de volumen de $D$ ($\mathcal{G} = M_V^{-1} D^T$). Esto asegura que el espacio de corrección sea ortogonal al subespacio libre de divergencia con respecto al producto interno de volumen.
3. Discretización: El artículo deriva fórmulas explícitas de SPH tanto para el operador de divergencia $D$ como para su gradiente adjunto $\mathcal{G}$. Crucialmente, el gradiente adjunto tiene en cuenta la asimetría en los coeficientes de la derivada de SPH (derivada de las longitudes de suavizado variables y las densidades) utilizando la simetría específica de la lista de vecinos y la ponderación de volumen.
4. Solución Numérica: La ecuación de proyección resultante $L\pi = r$ (donde $L = D\mathcal{G}$) forma un sistema lineal simétrico definido positivo que permite resolver el sistema eficientemente utilizando el método de Gradiente de Conjugado Precondicionado (PCG) sin necesidad de ensamblar explícitamente la matriz global. La implementación distingue entre partículas activas (donde se resuelve $\pi$) y partículas de frontera fijas.

Contribuciones Clave

• Formulación Adjunta Consistente: El artículo proporciona la primera formulación de proyección en SPMHD donde el gradiente de corrección se construye rigurosamente como el adjunto del operador de divergencia discreto elegido bajo una métrica ponderada por volumen. Esto resuelve las inconsistencias a nivel de discretización encontradas en intentos de proyección previos.
• Propiedad de Minimización de Energía: Se demuestra que la métrica de volumen es la métrica natural para la energía magnética. Por consiguiente, la proyección produce una corrección que minimiza la energía y no aumenta la energía magnética discreta, una propiedad crítica para la fidelidad física.
• Solucionador Eficiente: La formulación conduce a un sistema simétrico definido positivo resoluble mediante PCG, evitando el costo de $O(N^2)$ de los métodos de sumatoria directa utilizados en enfoques anteriores.
• Criterios de Parada Prácticos: Los autores introducen criterios de parada prácticos basados en la acumulación del error de divergencia entre los pasos de proyección, en lugar de requerir la convergencia a la precisión de la máquina en cada paso, equilibrando la exactitud con el costo computacional.

Resultados
El método fue probado utilizando pruebas de divergencia de tipo Dedner en dos dimensiones y simulaciones de colapso gravitatorio magnetizado en tres dimensiones.

• Pruebas de Divergencia: En pruebas 2D con tanto redes perturbadas como distribuciones de partículas aleatorias, el método redujo exitosamente el residuo de divergencia discreta al nivel del error de redondeo de punto flotante ($\sim 10^{-15}$) con suficientes iteraciones de PCG. Notablemente, corrigió errores en longitudes de onda comparables a la longitud de suavizado ($h$), algo con lo que los métodos previos tuvieron dificultades.
• Simulaciones de Colapso: En ejecuciones de colapso 3D, el método de proyección superó significativamente a la limpieza de divergencia hiperbólica estándar.
  • Reducción de Error: Con un criterio de parada práctico (reduciendo el error acumulado a una fracción del error de actualización de MHD), el método de proyección suprimió los errores de divergencia normalizados muy por debajo de los obtenidos con la limpieza de divergencia.
  • Robustez: Variar el intervalo de proyección ($N_{interval} = 1, 10, 100$) produjo estructuras de densidad y de plasma-$\beta$ consistentes. En contraste, la ejecución de limpieza de divergencia mostró diferencias cuantitativas en la estructura de flujo de salida y en el plasma $\beta$.
  • Costo: El costo computacional de la proyección fue modesto, oscilando aproximadamente entre el 1% y el 10% del costo de la actualización de SPMHD para intervalos prácticos ($N_{interval} \approx 10$). Incluso en pruebas de estrés extremas, el esquema pudo llevar los errores a niveles de redondeo.

Significancia y Reivindicaciones
El artículo afirma que esta formulación de proyección adjunta es una alternativa robusta y atractiva para controlar los errores de $\nabla \cdot \mathbf{B}$ en SPMHD y esquemas de MHD relacionados con partículas o mallas móviles.

Las ventajas clave destacadas incluyen:

1. Independencia de Parámetros: A diferencia de la limpieza de divergencia, que depende de parámetros físicos (velocidad de onda, escala de amortiguación) que pueden alterar los resultados astrofísicos, los parámetros del método de proyección (criterio de parada, intervalo) son adimensionales y se demostró que tienen un impacto mínimo en la solución dentro de los rangos probados.
2. Consistencia Física: El método garantiza que la corrección no inyecta energía magnética artificial.
3. Escalabilidad: El enfoque no se limita a SPMHD; los autores sugieren que la idea de la proyección basada en operadores es aplicable a otros esquemas de malla móvil o sin malla (por ejemplo, GIZMO, AREPO) con las modificaciones apropiadas.

Los autores concluyen que, si bien el método requiere elegir parámetros libres para la terminación y la frecuencia, su capacidad para producir resultados consistentes y de alta fidelidad sin los problemas de sensibilidad de los métodos de limpieza lo convierte en una opción superior para el control de la divergencia numérica en la MHD basada en partículas.