Résumé technique : Schéma de volume fini pour la magnétohydrodynamique relativiste viscorésistive de premier ordre

Énoncé du problème

La magnétohydrodynamique (MHD) relativiste est essentielle pour modéliser les fluides conducteurs dans les environnements de haute énergie et astrophysiques, tels que les plasmas quarks-gluons et les flux d'accrétion. Un défi persistant dans ce domaine est l'incorporation des coefficients de transport dissipatifs (viscosité, conduction thermique et résistivité) d'une manière qui préserve la causalité et la stabilité. Les théories traditionnelles de premier ordre produisent des équations paraboliques qui violent la causalité et induisent des instabilités à des vitesses élevées. Bien que les théories de Müller–Israel–Stein (MIS) résolvent ce problème en promouvant les courants dissipatifs au rang de variables dynamiques indépendantes, cette approche augmente considérablement la complexité et le nombre de degrés de liberté.

Les développements récents de la théorie de Bemfica–Disconzi–Noronha–Kovtun (BDNK) offrent une alternative de premier ordre minimale qui restaure la causalité et la stabilité en incluant des corrections de dérivées supérieures aux courants conservés sans introduire de nouvelles variables. Cependant, l'extension de la théorie BDNK à la MHD résistive présente une difficulté spécifique : bien que la formulation en « une-forme » de la MHD (traitant le flux magnétique comme une charge conservée) élimine la rigidité associée à la loi d'Ampère, elle rend la dynamique résistive parabolique. L'incorporation des corrections BDNK dans ce cadre de MHD en une-forme couplé et dissipatif, au sein d'un schéma numérique stable et efficace, est restée une lacune ouverte.

Méthodologie

Les auteurs présentent une implémentation numérique de la MHD relativiste viscorésistive de premier ordre basée sur le cadre BDNK, utilisant le code BHAC. La méthodologie implique :

1. Formulation théorique :

   • Les auteurs adoptent la limite ultra-relativiste (négligeant la densité de masse) et l'approche en une-forme, où les lignes de champ magnétique sont traitées comme des objets filamenteux conservés.
   • Les relations de comportement pour le tenseur énergie-impulsion ($T^{\mu\nu}$) et le tenseur de densité de flux dual ($J^{\mu\nu}$) incluent des corrections dissipatives de premier ordre pour la conductivité thermique, la viscosité et la résistivité.
   • Pour garantir la causalité, les corrections BDNK sont appliquées via une transformation de l'ordre zéro sur couche (on-shell). Cela introduit des temps de relaxation ($\tau_\epsilon, \tau_u, \tau_b$) qui convertissent les équations de diffusion paraboliques en équations de type télégraphe hyperboliques.
   • Crucialement, les auteurs identifient que trois termes BDNK sont insuffisants pour le secteur couplé de l'énergie-impulsion ultra-relativiste. Un quatrième terme, contrôlé par le coefficient $\tau_X$, est nécessaire pour garantir que le système reste hyperbolique et causal, empêchant un comportement elliptique dans le secteur acoustique.

2. Schéma numérique :

   • Les équations sont projetées sur des hypersurfaces espace-temps et décomposées en variables conservées ($U_I$) et flux ($F_I$), ainsi qu'en variables primitives ($P_I$) et termes sources ($S_I$).
   • Le schéma utilise une discrétisation par volumes finis avec un pas de temps TVD (Total Variation Diminishing) de second ordre et une reconstruction spatiale.
   • Les flux dissipatifs sont évalués par un calcul de différences centrales de second ordre des états d'interface reconstruits, avec des flux numériques de type upwind calculés via un solveur de Riemann approché.
   • Récupération des variables primitives : Une caractéristique clé de l'implémentation est une méthode algébrique efficace pour la récupération des variables primitives. Comme les corrections BDNK rendent les variables conservées dépendantes linéairement des dérivées temporelles des variables primitives, la récupération se réduit à une inversion de matrice locale ($S_I = \sum (M^{-1})_{IJ} [U'_J - U_J(0)]$). Cela évite le recours à des solveurs non linéaires itératifs typiquement requis dans les schémas dissipatifs.
   • Le code utilise des algorithmes de transport contraint (FCT) pour maintenir la condition de divergence nulle du champ magnétique ($\partial_i J^{it} = 0$) à la précision de la machine.

Contributions clés

• Extension de la BDNK à la MHD résistive : L'article construit un schéma numérique qui combine la dissipation visqueuse et résistive au sein du cadre BDNK pour la MHD relativiste.
• Identification d'un quatrième terme nécessaire : Les auteurs démontrent que si les termes BDNK minimaux peuvent créer des équations de type télégraphe pour des variables individuelles, un quatrième terme ($\tau_X$) est strictement nécessaire pour garantir l'hyperbolicité du système couplé énergie-impulsion dans la limite ultra-relativiste.
• Récupération efficace des variables primitives : Le travail démontre que les corrections BDNK permettent une procédure de récupération des variables primitives algébrique et non itérative, améliorant considérablement l'efficacité de calcul par rapport aux schémas dissipatifs traditionnels.
• Implémentation robuste : Le schéma est implémenté dans le code BHAC et validé à travers une gamme de régimes dissipatifs, allant du régime quasi-idéal au régime dominé par la diffusion.

Résultats

L'implémentation a été validée par une série de tests de référence en une et deux dimensions :

1. Test analytique : Une solution de télégraphe boostée pour l'évolution du champ magnétique a été utilisée pour vérifier l'exactitude. Les résultats numériques concordent avec la solution analytique avec une convergence de second ordre.
2. Tubes de choc 1D : Les simulations de tubes de choc avec des résistivités ($r_b$) et des paramètres de diffusion ($D_u, D_\epsilon$) variables ont montré que le schéma reste stable même avec des données initiales discontinues. Les résultats reproduisent correctement le comportement MHD idéal pour une faible dissipation et une dynamique dominée par la diffusion pour une forte dissipation.
3. Instabilité de Kelvin-Helmholtz (KH) : Des simulations 2D ont démontré que l'augmentation de la diffusion de quantité de mouvement ($D_u$) supprime l'instabilité, tandis que l'augmentation de la résistivité ($r_b$) l'accentue en relaxant le couplage entre la vitesse du fluide et les lignes de champ magnétique. Le schéma a maintenu des vitesses de front causales tout au long de l'évolution.
4. Vortex d'Orszag-Tang : Ce test a mis en évidence le rôle critique du quatrième terme BDNK ($\tau_X$). Lorsque $\tau_X = 0$, la simulation présentait un comportement anti-diffusif et plantait en raison de vitesses de front complexes. Avec $\tau_X > 0$, le système est resté stable avec des vitesses de front réelles et subluminales.
5. Feuille de Harris : Un test de feuille de Harris double perturbée, présentant des couches de renversement magnétique fines et de grands gradients, a été simulé en utilisant le raffinement de maillage adaptatif. Le schéma a réussi à résoudre la formation et la déformation de la couche de courant, maintenant une vitesse de front maximale de $v_{max} \approx 0,999$, confirmant une évolution causale.

Signification

L'article affirme que ce travail fournit un pont concret entre les développements théoriques récents en hydrodynamique relativiste (spécifiquement la théorie BDNK) et leur application aux simulations réalistes de flux magnétisés relativistes. En démontisant que la MHD relativiste dissipative peut être implémentée de manière numériquement robuste, causale et efficace sans introduire de variables dynamiques supplémentaires, les auteurs établissent une base pour simuler des phénomènes complexes dans des environnements tels que les plasmas quarks-gluons faiblement couplés, l'univers primordial et les magnétosphères de magnétars. La capacité de contrôler indépendamment les échelles résistives et visqueuses (en fixant le nombre de Prandtl magnétique) est soulignée comme une capacité cruciale pour étudier des phénomènes tels que les dynamos magnétiques et les éruptions de magnétars. Les auteurs notent également que la structure algébrique de leur récupération de variables primitives offre une voie prometteuse pour étendre ces corrections viscorésistives causales aux simulations de MHD en relativité générale.