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

Sintesi Tecnica: Instabilità Centrifuga dei Flussi Compressibili e la Stabilità Idrodinamica dei Dischi di Accrezione

Definizione del Problema
L'instabilità centrifuga nei flussi rotanti è tradizionalmente governata dal criterio di Rayleigh, il quale stabilisce che l'instabilità sorge quando il momento angolare specifico diminuisce con il raggio. Questo criterio è stato derivato per fluidi incompressibili e inviscidi. Tuttavia, analisi recenti sulla rotazione relativistica compressibile supportata dalla pressione (Gourgouliatos & Komissarov 2018b) hanno suggerito una generalizzazione di questo criterio in cui l'instabilità dipende esplicitamente dal numero di Mach del flusso ($M$). Nel limite Newtoniano, questo "criterio KG" postula che un flusso sia instabile se $d \ln \Psi / d \ln R < M^2$, dove $\Psi = \rho (v_\phi R)^2$. Ciò implica che anche i flussi stabili secondo il classico criterio di Rayleigh (dove $d(v_\phi R)^2/dR > 0$) potrebbero diventare instabili se il numero di Mach fosse sufficientemente elevato.

Questa previsione teorica sfida la comprensione convenzionale dei dischi di accrezione astrofisici. I dischi Keplariani sottili, che ruotano a velocità ipersoniche, sono ampiamente considerati idrodinamicamente stabili perché soddisfano il criterio di Rayleigh. Gli autori mettono in discussione se gli elevati numeri di Mach in questi sistemi, combinati con la compressibilità, possano in realtà renderli suscettibili all'instabilità centrifuga, un fattore spesso trascurato nelle analisi di stabilità standard che assumono l'incompressibilità o ignorano la dipendenza esplicita dal numero di Mach.

Metodologia
Per validare la versione Newtoniana del criterio KG, gli autori hanno eseguito simulazioni numeriche assiali di flussi rotanti transonici non relativistici. Le simulazioni hanno risolto le equazioni della dinamica dei fluidi compressibili ideali (equazioni di continuità, quantità di moto ed energia) utilizzando il codice AMR-VAC con un'equazione di stato di gas ideale ($\gamma = 5/3$).

La configurazione della simulazione prevedeva:

• Equilibrio Iniziale: Una configurazione con densità uniforme ($\rho = \rho_0$) e un profilo di rotazione $v_\phi(R)$ che transita da una rotazione a corpo rigido in un nucleo centrale a un decadimento a legge di potenza ($v_\phi \propto R^{-n}$) all'esterno del nucleo.
• Bilancio di Pressione: La distribuzione della pressione è stata derivata per soddisfare l'equazione di bilancio della forza radiale $dp/dR = \rho v_\phi^2 / R$.
• Perturbazione: L'equilibrio è stato perturbato con una piccola perturbazione della velocità azimutale ($\delta v_\phi = \sin(kz)$) per innescare potenziali instabilità.
• Spazio dei Parametri: Gli autori hanno variato l'indice della legge di potenza $n$ (compreso tra 0,2 e 1,0) e la pressione centrale $p_0$ per creare modelli con diversi numeri di Mach ($M_c$) al confine del nucleo. Ciò ha permesso di testare configurazioni che sono stabili secondo l'originale criterio di Rayleigh ma instabili secondo il criterio KG.

Contributi Chiave e Risultati

1. Validazione del Criterio KG: I risultati numerici hanno dimostrato un perfetto accordo con il criterio teorico KG. I modelli predetti come instabili (dove $2(1-n) < M^2$) hanno mostrato una rapida crescita delle perturbazioni, mentre i modelli predetti come stabili lo sono rimasti. Ciò conferma che l'instabilità centrifuga nei flussi compressibili dipende fortemente dal numero di Mach, una caratteristica assente nel classico limite incompressibile.
2. Evoluzione Non Lineare: Le simulazioni hanno rivelato lo sviluppo non lineare dell'instabilità. Le perturbazioni sono cresciute più rapidamente al confine interno del dominio instabile. Man mano che l'instabilità evolveva, il flusso ha sviluppato correnti radiali (dita) che si allontanano dal nucleo. Queste correnti interagiscono con il fluido circostante per formare strutture a forma di fungo e "cocoon" che esibiscono rotoli caratteristici dell'instabilità di Kelvin-Helmholtz. Questo processo ha portato a uno strato turbolento con un forte mescolamento, che si estende significativamente oltre la regione inizialmente instabile.
3. Rivalutazione della Stabilità dei Dischi di Accrezione: Gli autori hanno esteso la loro analisi euristica di stabilità per includere sia le forze di pressione che la gravità di un oggetto compatto centrale. Hanno derivato un criterio di instabilità generale che mantiene la dipendenza esplicita dal numero di Mach.
   • Nel limite di gravità nulla, questo criterio si riduce al criterio KG per la rotazione supportata dalla pressione.
   • Nel limite di rotazione nulla, si riduce al criterio di Schwarzschild per l'instabilità convettiva.
   • Nel caso generale, il criterio derivato è matematicamente equivalente allo standard criterio di Solberg-Høiland, che non mostra esplicitamente la dipendenza dal numero di Mach.

Significato e Rivendicazioni
L'articolo sostiene che la stabilità idrodinamica dei dischi di accrezione viene spesso spiegata usando l'originale criterio di Rayleigh, anche quando la rotazione è altamente supersonica. Derivando un criterio che include esplicitamente il numero di Mach e mostrando la sua equivalenza con il criterio di Solberg-Høiland, gli autori chiariscono che la stabilità dei dischi Keplariani sottili non è semplicemente il risultato del soddisfacimento della condizione di Rayleigh incomprimibile, ma è una conseguenza del bilancio specifico tra forze centrifughe e convettive descritto da Solberg-Høiland.

Gli autori concludono che, sebbene il criterio KG fornisca una formulazione utile per i flussi supportati dalla pressione e per i getti (dove la gravità è trascurabile), la dipendenza esplicita dal numero di Mach viene effettivamente assorbita nei gradienti di entropia e momento angolare nella formulazione generale di Solberg-Høiland. Di conseguenza, per i flussi isentropici (dove l'entropia è costante), il criterio si riduce all'originale criterio di Rayleigh, in linea con l'analisi di Rayleigh (1917) per i fluidi barotropici. Lo studio conferma che elevati numeri di Mach possono destabilizzare flussi supportati dalla pressione che sarebbero altrimenti stabili, un risultato rilevante per la dinamica dei getti astrofisici e l'interazione di venti supersonici con il plasma circostante, sebbene gli autori notino che sono necessarie simulazioni 3D complete per caratterizzare appieno il mescolamento non lineare innescato da questa instabilità.