Entity -- Hardware-agnostic Particle-in-Cell Code for Plasma Astrophysics. III: Higher-order shape functions & generalized field stencils

Résumé technique : Entité – Code de particules dans une cellule indépendant du matériel pour l'astrophysique des plasmas. III : Fonctions de forme d'ordre supérieur et maillages de champ généralisés

Énoncé du problème
Les simulations de particules dans une cellule (PIC) sont essentielles pour modéliser les processus de plasma astrophysique, allant de la turbulence et de la reconnexion magnétique à l'accélération des rayons cosmiques. Cependant, ces simulations font face à des défis majeurs concernant le bruit numérique, le chauffage artificiel et l'instabilité de Tcherenkov numérique. Ces artefacts résultent souvent de schémas d'interpolation d'ordre faible et de maillages de différences finies standards, qui peuvent masquer des effets physiques ou conduire à une accélération parasite des particules. Bien que l'augmentation de la résolution numérique (cellules de grille ou particules par cellule) soit une stratégie d'atténuation courante, elle est contrainte par les ressources informatiques. Les auteurs identifient un besoin de méthodes d'ordre supérieur pour supprimer le bruit et améliorer la stabilité sans nécessairement engendrer des coûts informatiques prohibitifs, en particulier dans le contexte du code PIC Entity, conçu pour les clusters GPU à l'échelle de l'exa.

Méthodologie
Ce travail (Article III de la série Entity) met en œuvre et teste deux avancées principales au sein du cadre Entity :

1. Dépôt de courant d'ordre supérieur et interpolation de champ :
   Les auteurs mettent en œuvre le schéma de dépôt de courant conservant la charge proposé par Esirkepov (2001), en l'étendant pour prendre en charge des fonctions de forme d'ordre arbitraire jusqu'à l'ordre 11 ($O=11$).

   • Fonctions de forme : Le code utilise des fonctions de forme dérivées des B-splines, construites par convolution successive. Ces fonctions satisfont les conditions de symétrie, de normalisation, de positivité, d'étendue finie et de régularité.
   • Cohérence : Pour assurer la cohérence avec le dépôt de courant, l'interpolation du champ de la grille vers les positions des particules est également mise à jour pour utiliser les mêmes fonctions de forme d'ordre supérieur, remplaçant les anciens schémas d'interpolation indirecte.
   • Note d'implémentation : Les ordres $O=10$ et $O=11$ nécessitent une arithmétique en double précision en raison de petites valeurs numériques aux limites du maillage, limitant leur utilisation sur certaines architectures GPU.

2. Maillages de champ généralisés pour le solveur de Maxwell :
   Le solveur de Faraday (évoluant le champ magnétique via $\partial_t B = -c\nabla \times E$) est mis à jour pour utiliser des maillages de différences finies généralisés.

   • Optimisation : Au lieu du maillage Yee standard d'ordre deux, les auteurs implémentent des maillages utilisant des points de grille adjacents supplémentaires. Les coefficients de ces maillages sont optimisés pour minimiser l'écart entre la vitesse de phase numérique et la vitesse physique de la lumière dans l'espace de Fourier.
   • Contrôle de la dispersion : Cette approche vise à supprimer l'instabilité de Tcherenkov numérique, qui se produit lorsque des particules relativistes voyagent plus vite que la vitesse numérique de la lumière sur la grille.
   • Préservation de la divergence : Les auteurs s'assurent que des configurations de maillage spécifiques satisfont la condition $\beta_{ij} = \delta_j$ pour imposer automatiquement $\nabla \cdot B = 0$ tout au long de la simulation.

Contributions et résultats clés
Les auteurs effectuent des tests de convergence extensifs pour valider les nouveaux schémas par rapport aux implémentations standards (dépôt Zigzag et maillage Yee) et à diverses configurations de filtres.

• Conservation de la charge : Le schéma Esirkepov mis à jour démontre une excellente conservation de la charge sur tous les ordres testés ($O=1$ à $O=11$), la précision étant limitée principalement par l'erreur de troncature plutôt que par l'ordre de la fonction de forme.
• Pouvoir d'arrêt : Dans les tests mesurant la force électrostatique d'arrêt sur des particules relativistes, les fonctions de forme d'ordre supérieur réduisent considérablement la perte d'énergie numérique. Le passage de l'ordre 1 à l'ordre 9 réduit la puissance d'arrêt effective, équivalent à augmenter la résolution numérique d'un facteur 3 à 4.
• Chauffage numérique et stabilité :
  • Plasma froid : Pour les plasmas froids avec des longueurs de Debye sous-résolues, les schémas d'ordre supérieur améliorent considérablement la stabilité. Le passage de l'ordre 1 à l'ordre 2 réduit la résolution de grille requise d'un facteur $\sim4$ ; le passage à l'ordre 3 la réduit d'un facteur $\sim8$.
  • Modes obliques : Dans les tests impliquant des faisceaux de rayons cosmiques, les formes d'ordre supérieur suppriment la croissance du bruit d'aliasing numérique. Si l'ordre 6 semble suffisant pour la convergence dans ces tests spécifiques, des ordres plus élevés retardent davantage l'apparition de l'instabilité.
  • Turbulence en décroissance : Dans la turbulence en décroissance 2D, les fonctions de forme d'ordre supérieur améliorent la conservation de l'énergie, réduisant l'erreur d'énergie totale jusqu'à deux ordres de grandeur par rapport au dépôt Zigzag standard.
• Performance des maillages de champ :
  • Particule unique : Des tests avec un électron relativiste unique montrent que les maillages optimisés (par exemple Cowan, min3, min6) réduisent considérablement les artefacts de Tcherenkov numérique par rapport au maillage Yee standard.
  • Plasma en dérive : Dans les simulations de plasma chaud en dérive, les maillages min3 et min6 suppriment efficacement la croissance de l'instabilité de Tcherenkov et du chauffage numérique associé.
  • Chocs relativistes : Dans une simulation de choc de paires relativiste, la combinaison du maillage min3 et des fonctions de forme d'ordre 3 stabilise le plasma amont, empêchant le chauffage artificiel et la filamentation de courant observés avec le maillage Yee. Cela conduit à des structures de rampe de choc et des spectres d'accélération de particules plus précis.
• Compromis de performance :
  • 1D : L'impact sur les performances des schémas d'ordre supérieur est négligeable.
  • 2D : Le coût augmente modérément (par exemple $\sim2,5\times$ pour l'ordre 3, $\sim24\times$ pour l'ordre 11).
  • 3D : Le coût devient prohibitif pour les ordres $O > 3$ (par exemple $\sim5\times$ pour l'ordre 3).
  • Atténuation : Les auteurs démontrent que le coût informatique accru de l'interpolation d'ordre supérieur peut être compensé en réduisant le nombre de passages de filtre de courant ou le nombre de particules par cellule, car les schémas d'ordre supérieur offrent intrinsèquement une plus grande stabilité.

Signification et affirmations
L'article affirme que, bien qu'aucune méthode unique ne résolve tous les défis numériques des simulations PIC, l'intégration de fonctions de forme d'ordre supérieur et de maillages de champ généralisés dans Entity offre une boîte à outils robuste pour améliorer la fidélité des simulations. Les auteurs soulignent que :

1. Stabilité vs Résolution : Les méthodes d'ordre supérieur permettent des simulations stables de plasmas froids et de chocs relativistes à des résolutions inférieures à celles précédemment possibles, réduisant potentiellement les coûts informatiques totaux malgré la surcharge par étape.
2. Suppression du bruit : Ces méthodes suppriment efficacement les effets de Tcherenkov numériques et le bruit d'aliasing, ce qui est crucial pour l'étude des instabilités physiques comme l'instabilité de Bell ou le flux de rayons cosmiques.
3. Flexibilité : Les maillages généralisés permettent aux utilisateurs d'ajuster le solveur pour des relations de dispersion spécifiques ou des symétries directionnelles (par exemple, utiliser le maillage min3 pour une stabilité générale ou Cowan/min6 pour de grands nombres CFL).
4. Recommandations pratiques : Les auteurs recommandent l'utilisation de fonctions de forme d'ordre 2 ou 3 comme compromis pratique entre précision et performance, et l'utilisation des maillages min3 (2D) ou min3D2 (3D) pour des simulations à usage général.

Le travail conclut que ces améliorations sont désormais disponibles dans Entity v1.3.0, permettant une modélisation plus précise et stable de phénomènes complexes de plasma astrophysique sur du matériel à l'échelle de l'exa.