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

Technische Samenvatting: Simulatie-gebaseerde parameterruimte voor schokken in transonische, sub-Kepleriaanse accretiestromen

Probleemstelling
Niet-dissipatieve, transonische, sub-Kepleriaanse accretiestromen naar zwarte gaten worden beheerst door twee geconserveerde parameters: specifieke energie ($\epsilon$) en specifiek impulsmoment ($l$). Hoewel analytische modellen gebaseerd op Rankine-Hugoniot (RH)-analyse en verticaal evenwicht een beperkte parameterruimte hebben geïdentificeerd waar stationaire schokken kunnen ontstaan, hebben eerdere meerdimensionale numerieke simulaties gesuggereerd dat schokken ook buiten deze theoretische limieten kunnen voorkomen. Deze simulaties hebben ook aangegeven dat schokken dynamisch (oscillerend) kunnen zijn in plaats van stationair. Echter, een uitgebreide, simulatie-gebaseerde mapping van de schok-parameterruimte ontbrak. Dit werk adresseert de kloof tussen analytische voorspellingen en de numerieke realiteit door de $(l, \epsilon)$ parameterruimte systematisch te scannen om te identificeren waar schokvorming plaatsvindt, het karakter van deze schokken te karakteriseren (stationair versus oscillerend), en de resulterende stromingsdynamica te analyseren.

Methodologie
De auteurs voerden meer dan 280 twee-dimensionale, ideale hydrodynamica-simulaties uit in een cilindrische $R-Z$ domein, waarbij zij uitgingen van axiale symmetrie. De simulaties maakten gebruik van een pseudo-Newtoniaans potentiaal ($\Phi(r) = -1/2(r-1)$) om een niet-roterend zwart gat te modelleren.

• Numeriek schema: De code maakte gebruik van een eindige volumemethode met tweede-orde nauwkeurigheid, gebruikmakend van een van Leer-limiter voor ruimtelijke reconstructie, een HLL Riemann-solver voor interface-fluxen, en een twee-traps strong stability preserving Runge-Kutta schema voor tijdintegratie.
• Randvoorwaarden: Inflow werd opgelegd aan de buitenste radiale grens ($R_{out}$) met constante primitieve variabelen afgeleid van analytische oplossingen voor specifieke $(l, \epsilon)$ paren, wat supersonische injectie garandeerde om akoestische feedback te voorkomen. Outflow-condities werden toegepast op de bovenste verticale grens ($Z_{out}$). Reflectiecondities werden gebruikt op de rotatieas ($R=0$) en het evenichtsvlak ($Z=0$). Materie binnen een straal $r_{BH} = 1.8$ werd geabsorbeerd.
• Parameterscan: De studie scande een rechthoekig domein van $l \in [1.4, 1.9]$ en $\epsilon \in [0.0001, 0.0061]$ met behulp van 286 punten. De simulatiedomeinen werden voor elk punt aangepast op basis van de locatie van het buitenste X-type sonische punt. Simulaties werden gedraaid tot $t=50000$ om een stationaire of periodieke toestand te waarborgen.

Belangrijkste bijdragen en resultaten
De primaire bijdrage is de identificatie van een schokvormings-parameterruimte die aanzienlijk groter is dan de analytisch voorspelde regio. De auteurs categoriseerden de parameterruimte in vijf onderscheidende regio's op basis van oplossingskenmerken:

1. Regio A (Stationaire Schokken): Deze regio bevat parameters waar stationaire, staande schokken ontstaan. De schok creëert een centrifugale barrière ondersteunde grenslaag (CENBOL). De auteurs merken op dat de numeriek geïdentificeerde grens voor deze regio zich uitstrekt naar lagere $l$-waarden dan de analytische voorspellingen. Deze uitbreiding wordt toegeschreven aan verticale bewegingen in de simulatie die materie nabij de evenicht accumuleren, waardoor de druk toeneemt en schokvorming vergemakkelijkt, zelfs wanneer de aannames van verticaal evenwicht in analytische modellen dit zouden verbieden.
2. Regio E (Geen Schok): Gelegen bij lagere $l$-waarden, vertonen deze oplossingen een centrifugale barrière die een dichtheidsophoping en een supersonische-naar-supersonische sprong veroorzaakt, maar er vormt zich geen echte schok (supersonische-naar-subsonische transitie). De stroming blijft supersonisch gedurende de gehele accretie naar het zwarte gat.
3. Regio B en C (Oscillerende Schokken): Deze regio's bevatten parameters waar schokken ontstaan maar oscilleren.
   • Regio B: Wordt gekenmerkt door hoogfrequente, kleine-amplitude oscillaties die voornamelijk beperkt zijn tot het onderste deel van het schokoppervlak (nabij de evenicht). Het bovenste deel van de CENBOL blijft statisch.
   • Regio C: Wordt gekenmerkt door grootschalige, laagfrequente oscillaties die het gehele schokoppervlak omvatten. Deze oplossingen vertonen significante turbulente vortex-structuren in de post-schok regio.
4. Regio D (Outflow-gedomineerd): Bij hoge $l$-waarden is de centrifugale barrière zo sterk dat accretie wordt verhinderd. De stroming ontwikkelt een naar buiten toe propagerende schok die de domeingrens bereikt, waarbij een subsonische stroming achterblijft. Meer dan 86% van de geïnjecteerde materie wordt uitgestoten als een outflow.

Parameterafhankelijkheid

• Schoklocatie: Voor een vaste $\epsilon$ beweegt de schoklocatie over het algemeen verder van het zwarte gat af naarmate $l$ toeneemt. De afhankelijkheid van $\epsilon$ voor een vaste $l$ is niet-monotoon; in sommige gevallen neemt de schoklocatie toe met $\epsilon$, terwijl deze in andere gevallen (met name bij lagere $l$) constant blijft of licht afneemt.
• Outflow: Outflow wordt gegenereerd in bijna alle simulaties door de gecombineerde effecten van de centrifugale barrière en thermische drukgradiënten. De ratio van outflow ten opzichte van inflow neemt significant toe naarmate $l$ stijgt, waarbij een transitie plaatsvindt van absorptie-gedomineerd (Regio's A en E) naar outflow-gedomineerd (Regio D).

Betekenis en Claims
Het artikel claimt dat meerdimensionale numerieke simulaties een bredere parameterruimte voor schokvorming onthullen dan voorheen begrepen via analytische methoden. Specifiek staat de inclusie van verticale dynamica toe dat schokken kunnen vormen bij lagere impulsmomenten dan voorspeld door modellen van verticaal evenwicht.

De auteurs benadrukken de astrofysische relevantie van deze bevindingen:

• Corona-vorming: De door schokken gecomprimeerde, geometrisch dikke post-schok stroming (CENBOL) fungeert als een hete, optisch dunne regio die in staat is tot inverse-Compton-verstrooiing van laagenergetische fotonen van een Kepleriaanse schijf (of zelf-Comptonisering indien de schijf afwezig is). Dit biedt een natuurlijk mechanisme voor de vorming van de corona die wordt waargenomen in röntgen-dubbelsterren en actieve sterrenstelsels in harde spectrale toestanden.
• Quasi-Periodieke Oscillaties (QPO's): De identificatie van oscillerende schokregio's (B en C) ondersteunt de theoretische link tussen schokoscillaties en waargenomen QPO's in hoogenergetische straling.
• Outflows: De simulaties produceren consistent outflows, die frequent worden waargenomen in harde spectrale toestanden, wat suggereert dat het schokmechanisme een verenigde drijfveer is voor zowel de corona als de outflow.

De auteurs blijven bescheiden over hun reikwijdte en merken op dat de resultaten specifiek zijn voor niet-roterende zwarte gaten met gebruik van een pseudo-Newtoniaanse potentiaal en ideale hydrodynamica. Zij erkennen dat het opnemen van viscositeit, radiatieve verliezen, magnetische velden of algemene relativistische effecten voor roterende zwarte gaten de grenzen van de parameterruimte zou kunnen veranderen. Zij stellen echter dat het huidige werk een fundamentele, simulatie-gebaseerde kaart van schokvorming vaststelt die de huidige analytische limieten overstijgt.