Centrifugal instability of compressible flows and the hydrodynamic stability of accretion disks

Technische Samenvatting: Centrifugale Instabiliteit van Samendrukbare Stromingen en de Hydrodynamische Stabiliteit van Accretieschijven

Probleemstelling
De centrifugale instabiliteit in roterende stromingen wordt traditioneel beheerst door het Rayleigh-criterium, dat stelt dat instabiliteit ontstaat wanneer het specifieke impulsmoment met de straal afneemt. Dit criterium werd afgeleid voor onviskeuze, incompressibele vloeistoffen. Echter, recente analyses van drukgestuurde samendrukbare relativistische rotatie (Gourgouliatos & Komissarov 2018b) suggereerden een generalisatie van dit criterium waarbij de instabiliteit expliciet afhankelijk is van het Mach-getal ($M$) van de stroming. In de Newtoniaanse limiet stelt dit "KG-criterium" dat een stroming instabiel is als $d \ln \Psi / d \ln R < M^2$, waarbij $\Psi = \rho (v_\phi R)^2$. Dit impliceert dat zelfs stromingen die stabiel zijn volgens het klassieke Rayleigh-criterium (waarbij $d(v_\phi R)^2/dR > 0$) instabiel kunnen worden als het Mach-getal voldoende hoog is.

Deze theoretische voorspelling daagt het conventionele begrip van astrofysische accretieschijven uit. Dunne Keplerse schijven, die met hypersonische snelheden roteren, worden breed beschouwd als hydrodynamisch stabiel omdat ze aan het Rayleigh-criterium voldoen. De auteurs stellen de vraag of de hoge Mach-getallen in deze systemen, in combinatie met samendrukbaarheid, hen in werkelijkheid vatbaar kunnen maken voor centrifugale instabiliteit, een factor die vaak over het hoofd wordt gezien in standaard stabiliteitsanalyses die incompressibiliteit veronderstellen of de expliciete Mach-afhankelijkheid negeren.

Methodologie
Om de Newtoniaanse versie van het KG-criterium te valideren, voerden de auteurs axiaal symmetrische numerieke simulaties uit van niet-relativistische, transsonische roterende stromingen. De simulaties losten de vergelijkingen van de ideale samendrukbare fluïdumdynamica op (continuïteit, impuls- en energievergelijkingen) met behulp van de AMR-VAC code met een ideale gaswet ($\gamma = 5/3$).

De simulatieopstelling omvatte:

• Initiële Evenwichtstoestand: Een configuratie met een uniforme dichtheid ($\rho = \rho_0$) en een rotatieprofiel $v_\phi(R)$ dat overgaat van een starre rotatie in een centrale kern naar een machtswetmatige afname ($v_\phi \propto R^{-n}$) buiten de kern.
• Drukbalans: De drukverdeling werd afgeleid om te voldoen aan de radiale krachtbalansvergelijking $dp/dR = \rho v_\phi^2 / R$.
• Perturbatie: Het evenwicht werd geperturbeerd met een kleine azimuthale snelheidsverstoring ($\delta v_\phi = \sin(kz)$) om potentiële instabiliteiten te triggeren.
• Parametruimte: De auteurs varieerden de machtswet-index $n$ (variërend van 0,2 tot 1,0) en de centrale druk $p_0$ om modellen te creëren met verschillende Mach-getallen ($M_c$) bij de grens van de kern. Dit maakte het mogelijk om configuraties te testen die stabiel zijn volgens het oorspronkelijke Rayleigh-criterium, maar instabiel volgens het KG-criterium.

Belangrijkste Bijdragen en Resultaten

1. Validatie van het KG-criterium: De numerieke resultaten vertoonden een perfecte overeenstemming met het theoretische KG-criterium. Modellen die voorspeld werden als instabiel (waar $2(1-n) < M^2$) vertoonden een snelle groei van perturbaties, terwijl modellen die voorspeld werden als stabiel dat ook bleven. Dit bevestigt dat de centrifugale instabiliteit in samendrukbare stromingen sterk afhankelijk is van het Mach-getal, een kenmerk dat afwezig is in de klassieke incompressibele limiet.
2. Nietlineaire Evolutie: De simulaties onthulden de nietlineaire ontwikkeling van de instabiliteit. Perturbaties groeiden het snelst bij de binnenste grens van het instabiele domein. Terwijl de instabiliteit evolueerde, ontwikkelde de stroming radiale stromen (vingers) die weg bewogen van de kern. Deze stromen interageerden met de omringende vloeistof om paddenstoelachtige structuren en kokons te vormen die rollen vertonen die karakteristiek zijn voor de Kelvin-Helmholtz-instabiliteit. Dit proces leidde tot een turbulente laag met sterke menging, die zich aanzienlijk buiten de initiële lineair instabiele regio uitstrekte.
3. Herwaardering van de Stabiliteit van Accretieschijven: De auteurs breidden hun heuristische stabiliteitsanalyse uit om zowel drukversterkende krachten als de zwaartekracht van een centraal compact object te bevatten. Ze leidden een algemeen instabiliteitscriterium af dat de expliciete afhankelijkheid van het Mach-getal behoudt.
   • In de limiet van verwaarloosbare zwaartekracht reduceert dit criterium tot het KG-criterium voor drukgestuurde rotatie.
   • In de limiet van verwaarloosbare rotatie reduceert het tot het Schwarzschild-criterium voor convectieve instabiliteit.
   • In het algemene geval is het afgeleide criterium wiskundig equivalent aan het standaard Solberg-Høiland-criterium, dat de expliciete Mach-getal afhankelijkheid niet vertoont.

Betekenis en Claims
Het artikel beweert dat de hydrodynamische stabiliteit van accretieschijven vaak wordt verklaard met het oorspronkelijke Rayleigh-criterium, zelfs wanneer de rotatie sterk supersonisch is. Door een criterium af te leiden dat expliciet het Mach-getal bevat en door de equivalentie met het Solberg-Høiland-criterium aan te tonen, verhelderen de auteurs dat de stabiliteit van dunne Keplerse schijven niet louter een gevolg is van het voldoen aan de incompressibele Rayleigh-conditie, maar een gevolg is van de specifieke balans tussen centrifugale en convectieve krachten zoals beschreven door Solberg-Høiland.

De auteurs concluderen dat hoewel het KG-criterium een nuttige formulering biedt voor drukgestuurde stromingen en jets (waar zwaartekracht verwaarloosbaar is), de expliciete Mach-getal afhankelijkheid effectief wordt geabsorbeerd in de entropie- en impulsmomentgradiënten in de algemene Solberg-Høiland-formulering. Voor isentropische stromingen (waar entropie constant is) reduceert het criterium derhalve tot het oorspronkelijke Rayleigh-criterium, consistent met de analyse van Rayleigh (1917) van barotrope vloeistoffen. De studie bevestigt dat hoge Mach-getallen drukgestuurde stromingen instabiel kunnen maken die anders stabiel zouden zijn, een bevinding die relevant is voor de dynamica van astrofysische jets en de interactie van supersonische winden met omliggend plasma, hoewel de auteurs opmerken dat volledige 3D-simulaties nodig zijn om de nietlineaire menging die door deze instabiliteit wordt getriggerd volledig te karakteriseren.