Simulation based parameter space for shock in transonic, sub-Keplerian accretion flow onto non-rotating black holes

Résumé Technique : Espace de Paramètres par Simulation pour les Chocs dans un Écoulement d'Accrétion Transsonique et Sub-Keplerien

Énoncé du Problème
Les écoulements d'accrétion transsoniques, sub-kepleriens et non dissipatifs sur les trous noirs sont régis par deux paramètres conservés : l'énergie spécifique ($\epsilon$) et le moment angulaire spécifique ($l$). Bien que les modèles analytiques basés sur l'analyse de Rankine-Hugoniot (RH) et l'équilibre vertical aient identifié un espace de paramètres limité où des chocs stationnaires peuvent se former, des simulations numériques multidimensionnelles antérieures ont suggéré que les chocs peuvent se produire en dehors de ces limites théoriques. Ces simulations ont également indiqué que les chocs peuvent être dynamiques (oscillatoires) plutôt que stationnaires. Cependant, une cartographie systématique de l'espace de paramètres des chocs basée sur des simulations faisait défaut. Ce travail comble l'écart entre les prédictions analytiques et la réalité numérique en scannant systématiquement l'espace de paramètres $(l, \epsilon)$ pour identifier où la formation de chocs se produit, caractériser la nature de ces chocs (stationnaires vs oscillatoires) et analyser la dynamique de l'écoulement qui en résulte.

Méthodologie
Les auteurs ont réalisé plus de 280 simulations d'hydrodynamique idéale en deux dimensions dans un domaine cylindrique $R-Z$, en supposant l'axi-symétrie. Les simulations ont utilisé un potentiel pseudo-newtonien ($\Phi(r) = -1/2(r-1)$) pour modéliser un trou noir non rotatif.

• Schéma Numérique : Le code a employé une méthode de volume fini de second ordre, utilisant un limiteur de van Leer pour la reconstruction spatiale, un solveur de Riemann HLL pour les flux interfaciaux, et un schéma Runge-Kutta à deux étapes de type "strong stability preserving" pour l'intégration temporelle.
• Conditions aux Limites : L'injection a été imposée à la limite radiale extérieure ($R_{out}$) avec des variables primitives constantes dérivées de solutions analytiques pour des paires $(l, \epsilon)$ spécifiques, assurant une injection supersonique pour empêcher le feedback acoustique. Des conditions de sortie ont été appliquées à la limite verticale supérieure ($Z_{out}$). Des conditions de réflexion ont été utilisées sur l'axe de rotation ($R=0$) et sur le plan équatorial ($Z=0$). La matière à l'intérieur d'un rayon $r_{BH} = 1,8$ a été absorbée.
• Balayage de Paramètres : L'étude a balayé un domaine rectangulaire de $l \in [1,4, 1,9]$ et $\epsilon \in [0,0001, 0,0061]$ en utilisant 286 points. Les domaines de simulation ont été ajustés pour chaque point en fonction de l'emplacement du point sonique de type X le plus externe. Les simulations ont été exécutées jusqu'à $t=50000$ pour garantir un comportement stationnaire ou périodique.

Contributions Clés et Résultats
La contribution principale est l'identification d'un espace de paramètres de formation de chocs qui est significativement plus large que la région prédite analytiquement. Les auteurs ont classé l'espace de paramètres en cinq régions distinctes selon les caractéristiques des solutions :

1. Région A (Chocs Stationnaires) : Cette région contient des paramètres où des chocs stationnaires se forment. Le choc crée une couche limite supportée par une barrière centrifuge (CENBOL). Les auteurs notent que la limite identifiée numériquement pour cette région s'étend vers des valeurs de $l$ plus faibles que les prédictions analytiques. Cette extension est attribuée aux mouvements verticaux dans la simulation qui accumulent la matière près de l'équateur, augmentant la pression et facilitant la formation du choc, même lorsque les hypothèses d'équilibre vertical des modèles analytiques l'interdiraient.
2. Région E (Pas de Choc) : Située à des valeurs de $l$ plus faibles, ces solutions présentent une barrière centrifuge causant un empilement de densité et un saut supersonique-vers-supersonique, mais aucun véritable choc (transition supersonique-vers-subsonique) ne se forme. L'écoulement reste supersonique tout au long du processus et accrète sur le trou noir.
3. Régions B et C (Chocs Oscillatoires) : Ces régions contiennent des paramètres où des chocs se forment mais oscillent.
   • Région B : Caractérisée par des oscillations de haute fréquence et de faible amplitude, confinées principalement à la partie inférieure de la surface du choc (près de l'équateur). La partie supérieure du CENBOL reste statique.
   • Région C : Caractérisée par des oscillations de grande amplitude et de basse fréquence impliquant toute la surface du choc. Ces solutions présentent des vortex turbulents importants dans la région post-choc.
4. Région D (Dominée par les Écoulements Sortants) : À des valeurs de $l$ élevées, la barrière centrifuge est si forte qu'elle empêche l'accrétion. L'écoulement développe un choc se propageant vers l'extérieur qui atteint la limite du domaine, laissant derrière lui un écoulement subsonique. Plus de 86 % de la matière injectée est éjectée sous forme d'écoulement sortant (outflow).

Dépendance des Paramètres

• Localisation du Choc : Pour un $\epsilon$ fixé, la position du choc s'éloigne généralement du trou noir à mesure que $l$ augmente. La dépendance vis-à-vis de $\epsilon$ pour un $l$ fixé est non monotone ; dans certains cas, la position du choc augmente avec $\epsilon$, tandis que dans d'autres (particulièrement à bas $l$), elle reste constante ou diminue légèrement.
• Écoulement Sortant (Outflow) : Un écoulement sortant est généré dans presque toutes les simulations en raison des effets combinés de la barrière centrifuge et des gradients de pression thermique. Le rapport entre l'écoulement sortant et l'écoulement entrant augmente considérablement à mesure que $l$ augmente, passant d'une dominance par absorption (Régions A et E) à une dominance par écoulement sortant (Région D).

Signification et Revendications
L'article affirme que les simulations numériques multidimensionnelles révèlent un espace de paramètres pour la formation de chocs plus large que ce qui était compris précédemment via les méthodes analytiques. Spécifiquement, l'inclusion de la dynamique verticale permet aux chocs de se former à des moments angulaires plus faibles que ceux prédits par les modèles d'équilibre vertical.

Les auteurs soulignent la pertinence astrophysique de ces découvertes :

• Formation de la Couronne : L'écoulement post-choc, compressé par le choc et géométriquement épais (CENBOL), agit comme une région chaude et optiquement mince capable de Comptoniser inversement les photons de basse énergie provenant d'un disque Keplerien (ou auto-Comptoniser si le disque est absent). Cela fournit un mécanisme naturel pour la formation de la couronne observée dans les binaires X et les noyaux actifs de galaxies dans les états spectraux durs.
• Oscillations Quasi-Périodiques (QPO) : L'identification des régions de chocs oscillatoires (B et C) soutient le lien théorique entre les oscillations de choc et les QPO observées dans le rayonnement de haute énergie.
• Écoulements Sortants : Les simulations produisent de manière auto-cohérente des écoulements sortants, qui sont fréquemment observés dans les états spectraux durs, suggérant que le mécanisme de choc est un moteur unifié pour la fois la couronne et l'écoulement sortant.

Les auteurs font preuve de modestie quant à leur portée, notant que les résultats sont spécifiques aux trous noirs non rotatifs utilisant un potentiel pseudo-newtonien et l'hydrodynamique idéale. Ils reconnaissent que l'inclusion de la viscosité, des pertes radiatives, des champs magnétiques ou des effets de relativité générale pour les trous noirs rotatifs pourrait modifier les limites de l'espace de paramètres. Cependant, ils affirment que le présent travail établit une carte fondamentale, basée sur des simulations, de la formation de chocs qui s'étend au-delà des limites analytiques actuelles.