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

Résumé Technique : Une formulation de projection par adjoint pour l'application de la contrainte $\nabla \cdot \mathbf{B} = 0$ en Magnétohydrodynamique à Particules lissées (SPMHD)

Énoncé du Problème
En magnétohydrodynamique (MHD), la contrainte solénoidale $\nabla \cdot \mathbf{B} = 0$ est une exigence physique fondamentale. Bien que les codes basés sur des mailles utilisent souvent le transport contraint (CT) pour préserver cette contrainte de manière discrète, les méthodes sans maillage comme la SPMHD ne possèdent pas de topologie naturelle face-arête pour appliquer directement le CT. Par conséquent, l'approche dominante en SPMHD a été le « nettoyage de la divergence » (divergence cleaning), qui fait évoluer un champ scalaire auxiliaire pour propager et amortir les erreurs de divergence. Cependant, les méthodes de nettoyage introduisent des paramètres physiques (vitesse de propagation et échelle de temps d'amortissement) qui peuvent affecter les résultats de la simulation si ceux-ci ne sont pas choisis de manière cohérente, particulièrement dans les systèmes présentant de larges plages dynamiques. Les tentatives de projection antérieures en SPMHD, telles que celles de Price et Monaghan, utilisaient des sommes de type gravitationnel pour les corrections, mais souffraient d'incohérences entre l'opérateur de divergence utilisé pour mesurer les erreurs et l'opérateur de gradient utilisé pour calculer les corrections, entraînant de mauvaises performances sur le bruit à courtes longueurs d'onde.

Méthodologie
Le présent article propose une méthode de projection qui impose la contrainte discrète de divergence nulle en résolvant un problème elliptique après chaque mise à jour MHD. L'innovation centrale réside dans la construction d'un opérateur de correction qui est strictement l'adjoint de l'opérateur de divergence discrète utilisé pour mesurer l'erreur, défini au sein d'une métrique pondérée par le volume.

1. Formulation : La méthode cherche un champ magnétique corrigé $\mathbf{B}^{n+1}$ qui minimise le changement de la norme d'énergie magnétique discrète (définie par les volumes de particules $V_i = m_i/\rho_i$) sous la contrainte $D\mathbf{B}^{n+1} = 0$, où $D$ est l'opérateur de divergence SPH discret.
2. Gradient Adjoint : En introduisant un champ de multiplicateurs de Lagrange $\pi$, la correction est dérivée de $\mathbf{B}^{n+1} = \mathbf{B}^* - \mathcal{G}\pi$. L'opérateur $\mathcal{G}$ est explicitement construit comme l'adjoint de $D$ sous la métrique volumique ($\mathcal{G} = M_V^{-1} D^T$). Cela garantit que l'espace de correction est orthogonal au sous-espace sans divergence par rapport au produit scalaire volumique.
3. Discrétisation : L'article dérive des formules SPH explicites pour l'opérateur de divergence $D$ et son gradient adjoint $\mathcal{G}$. Crucialement, le gradient adjoint tient compte de l'asymétrie des coefficients de dérivée SPH (résultant de la variation des longueurs de lissage et des densités) en utilisant la symétrie spécifique de la liste de voisins et la pondération volumique.
4. Solution Numérique : L'équation de projection résultante $L\pi = r$ (où $L = D\mathcal{G}$) forme un système linéaire symétrique semi-défini positif. Cela permet de résoudre le système efficacement à l'aide de la méthode de Gradient Conjugué Préconditionné (PCG) sans avoir à assembler explicitement la matrice globale. L'implémentation distingue les particules actives (où $\pi$ est résolu) et les particules de bord fixes.

Contributions Clés

• Formulation Adjointe Cohérente : L'article fournit la première formulation de projection SPMHD où la correction de gradient est rigoureusement construite comme l'adjoint de l'opérateur de divergence discret choisi sous une métrique pondérée par le volume. Cela résout les incohérences de discrétisation trouvées dans les tentatives de projection antérieures.
• Propriété de Minimisation de l'Énergie : Il est démontré que la métrique volumique est la métrique naturelle pour l'énergie magnétique. Par conséquent, la projection produit une correction minimisant l'énergie qui ne augmente pas l'énergie magnétique discrète, une propriété critique pour la fidélité physique.
• Solveur Efficace : La formulation conduit à un système symétrique semi-défini positif résoluble par PCG, évitant le coût en $O(N^2)$ des méthodes de sommation directe utilisées dans les approches précédentes.
• Critères d'Arrêt Pratiques : Les auteurs introduisent des critères d'arrêt pratiques basés sur l'accumulation de l'erreur de divergence entre les étapes de projection, plutôt que d'exiger une convergence à la précision machine à chaque étape, équilibrant ainsi précision et coût de calcul.

Résultats
La méthode a été testée à l'aide de tests de divergence de type Dedner en 2D et de simulations de l'effondrement gravitationnel magnétisé en 3D.

• Tests de Divergence : Dans des tests 2D avec des réseaux perturbés et des distributions de particules aléatoires, la méthode a réussi à réduire le résidu de divergence discrète au niveau de la précision de virgule flottante ($\sim 10^{-15}$) avec un nombre suffisant d'itérations PCG. Notamment, elle a corrigé les erreurs à des longueurs d'onde comparables à la longueur de lissage ($h$), là où les méthodes précédentes peinaient.
• Simulations d'Effondrement : Dans les simulations d'effondrement 3D, la méthode de projection a nettement surpassé le nettoyage de divergence hyperbolique standard.
  • Réduction de l'Erreur : Avec un critère d'arrêt pratique (réduisant l'erreur accumulée à une fraction de l'erreur de mise à jour MHD), la méthode de projection a supprimé les erreurs de divergence normalisées bien en dessous de celles obtenues avec le nettoyage de divergence.
  • Robustesse : La variation de l'intervalle de projection ($N_{interval} = 1, 10, 100$) a produit des structures de densité et de plasma-$\beta$ cohérentes. En revanche, la simulation avec nettoyage de divergence a montré des différences quantitatives dans la structure de l'écoulement sortant et le plasma $\beta$.
  • Coût : Le coût de calcul de la projection était modeste, allant d'environ 1 % à 10 % du coût de la mise à jour SPMHD pour des intervalles pratiques ($N_{interval} \approx 10$). Même lors de tests de stress extrêmes, le schéma pouvait conduire les erreurs au niveau de la précision machine.

Signification et Revendications
L'article affirme que cette formulation de projection par adjoint est une alternative robuste et attractive au nettoyage de divergence pour contrôler les erreurs de $\nabla \cdot \mathbf{B}$ en SPMHD et dans les schémas MHD à particules ou à mailles mobiles apparentés.

Les principaux avantages mis en avant incluent :

1. Indépendance des Paramètres : Contrairement au nettoyage de divergence, qui dépend de paramètres physiques (vitesse d'onde, échelle d'amortissement) pouvant altérer les résultats astrophysiques, les paramètres de la méthode de projection (critère d'arrêt, intervalle) sont adimensionnels et ont montré un impact minimal sur la solution dans les plages testées.
2. Cohérence Physique : La méthode garantit que la correction n'injecte pas d'énergie magnétique artificielle.
3. Scalabilité : L'approche n'est pas limitée à la SPMHD ; les auteurs suggèrent que l'idée de projection basée sur les opérateurs est applicable à d'autres schémas sans maillage ou à mailles mobiles (ex: GIZMO, AREPO) avec des modifications appropriées.

Les auteurs concluent que, bien que la méthode nécessite de choisir des paramètres libres pour la terminaison et la fréquence, sa capacité à produire des résultats cohérents et de haute fidélité sans les problèmes de sensibilité des méthodes de nettoyage en fait un choix supérieur pour le contrôle de la divergence numérique en MHD par particules.