Numerical thermalization in $n$-D particle-in-cell simulations

Résumé Technique : Thermalisation Numérique dans les Simulations de Particules dans une Cellule (PIC) en n-D

Énoncé du Problème
Les simulations de particules dans une cellule (PIC) sont largement utilisées pour modéliser le comportement cinétique des plasmas en représentant les particules physiques par des macroparticules dotées de fonctions de forme finies. Bien qu'elles soient souvent supposées résoudre l'équation de Vlasov sans collision grâce aux effets de lissage de l'interpolation sur la grille, les simulations PIC souffrent fréquemment d'une collisionalité artificiellement élevée. Cette « thermalisation numérique » survient lorsque le poids de la macroparticule (le rapport entre les particules physiques et les macroparticules) est important, amplifiant les forces inter-particulaires malgré le lissage de la grille. Dans les dimensions supérieures (2D et 3D), parvenir à un temps de thermalisation numérique ($\tau_{num}$) plus long que le temps de thermalisation physique ($\tau_{phys}$) est souvent intraitable, même lorsque la grille résout la longueur de Debye. Cette collisionalité artificielle distord les distributions de vitesse et affecte des quantités clés d'intérêt (par exemple, le transport, l'émission, la rupture) sur des échelles de temps de simulation $t \gtrsim \tau_{num}$, pouvant potentiellement conduire à une convergence fausse ou à des prédictions physiques incorrectes.

Méthodologie
Les auteurs valident une théorie cinétique des temps de collision numériques par rapport à des simulations PIC électrostatiques en 1D, 2D et 3D.

1. Simulations : Des simulations PIC électrostatiques à espèce unique ont été réalisées avec des conditions initiales uniformes. La fonction d'autocorrélation de la vitesse (VACF), $Z(t)$, a été calculée à partir des données de trajectoire des particules. Le temps de collision numérique a été défini comme le temps de décroissance exponentielle de $Z(t)$.
2. Paramètres : Les simulations ont fait varier le poids de la macroparticule ($W$), l'espacement de la grille ($\Delta$) et le rayon de la forme de la macroparticule ($R$). Des simulations non filtrées (où $R = \Delta/2$) et filtrées (où $R > \Delta/2$ via un filtrage dans l'espace $k$) ont été effectuées pour isoler les effets de forme du chauffage de la grille.
3. Théorie : L'étude emploie la théorie cinétique développée par Okuda, Birdsall et Langdon. Ce cadre intègre la fonction de forme et le poids de la macroparticule dans une équation cinétique de type Lenard-Balescu. La théorie dérive des opérateurs de collision pour les macroparticules façonnées, en tenant compte de la modification du noyau de force inter-particulaire par la fonction de forme $S(k)$ et le poids $W$.
4. Comparaison : Les prédictions théoriques pour le temps de collision ont été dérivées à l'aide d'approximations dynamiques, statiques et à faible vitesse de l'intégrale de traînée. Celles-ci ont été comparées directement aux taux de décroissance observés dans les simulations PIC.

Contributions Clés et Résultats

• Validation de la Théorie Cinétique : L'article démontre que la théorie cinétique d'Okuda, Birdsall et Langdon prédit avec précision le temps de collision numérique dans les simulations PIC 1D, 2D et 3D. La théorie capture correctement la dépendance de $\tau_{num}$ vis-à-vis du poids de la macroparticule, de la taille de la forme et de la dimensionnalité.
• Interprétation du Noyau de Force : L'étude clarifie que les collisions numériques sont pilotées par une force inter-particulaire effective où la fonction de forme finie affaiblit et lisse la force, tandis que le poids de la macroparticule l'amplifie. La théorie montre que le rapport entre les temps de collision numériques et physiques est indépendant des rapports de masse et de charge des espèces, dépendant principalement de $W$ et $R$.
• Effets de Dimensionnalité : Les résultats soulignent qu'atteindre le temps de thermalisation physique est nettement plus difficile en 3D qu'en 1D ou 2D pour des tailles de macroparticules qui résolvent la longueur de Debye. En 3D, le temps de collision évolue de manière logarithmique avec le paramètre de plasma $\Lambda$ dans la limite d'un petit rayon de particule ($R \ll \lambda_D$).
• Analyse des Compromis : L'article quantifie le compromis entre la réduction de la thermalisation numérique et le maintien de l'exactitude physique. Augmenter la taille de la macroparticule ($R$) pour réduire les collisions (en lissant les forces) affaiblit simultanément le champ de Vlasov moyen et altère les processus sans collision (par exemple, la dispersion et l'amortissement des ondes). Il existe une limite de « champ distordu » où l'augmentation supplémentaire de $R$ viole les contraintes physiques de résolution spatiale.
• Formes d'Ordre Supérieur : L'analyse s'étend aux fonctions de forme d'ordre supérieur (par exemple, les B-splines). Il est montré qu'augmenter l'ordre de la spline est mathématiquement équivalent à augmenter la largeur effective de la particule. Ainsi, l'utilisation de formes d'ordre supérieur ne modifie pas fondamentalement l'échelle de la thermalisation numérique, mais augmente effectivement la taille de la particule.
• Application à des Systèmes Réels : Les auteurs appliquent les contraintes dérivées à plusieurs études PIC existantes (par exemple, points chauds d'ICF, reconnexion magnétique, décharges de Penning, propulseurs à effet Hall, décharges de streamers et injecteurs de faisceaux neutres d'ITER). Dans de nombreux cas 3D, le temps de thermalisation numérique est trouvé plus court que le temps de collision physique, suggérant que la thermalisation numérique peut dominer ou influencer significativement les résultats, à moins que des stratégies de mitigation spécifiques (comme les collisions de Monte Carlo ou un raffinement extrême de la grille) ne soient employées.

Signification
La principale signification de ce travail est la fourniture d'un cadre validé et prédictif pour estimer l'échelle de temps des interactions coulombiennes auto-cohérentes dans les simulations PIC. En établissant que la théorie cinétique prédit avec précision $\tau_{num}$ en fonction de l'espacement de la grille et du poids de la macroparticule, l'article permet aux utilisateurs de PIC de :

1. Caractériser la Pertinence : Déterminer a priori si la thermalisation numérique impactera la physique d'intérêt spécifique (par exemple, si $\tau_{num}$ est plus court que le temps de simulation ou le temps de collision physique).
2. Définir des Contraintes : Définir les régimes de fonctionnement permis pour le poids de la macroparticule et la taille de la forme afin de garantir que les artefacts numériques n'obscurcissent pas les phénomènes physiques.
3. Éviter une Fausse Convergence : Reconnaître que la convergence par rapport au nombre de particules par cellule ne garantit pas l'élimination de la thermalisation numérique ; un nombre suffisant de macroparticules est requis pour résoudre le taux de collision physique, et non seulement les caractéristiques du bruit.

L'article conclut que bien que les méthodes PIC soient puissantes, l'hypothèse « sans collision » est souvent invalide dans les simulations 3D avec des paramètres standards. Les utilisateurs doivent explicitement prendre en compte la thermalisation numérique pour assurer la fidélité des processus de transport et de collision à long terme.