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

Resumo Técnico: Instabilidade Centrífuga de Fluxos Compressíveis e a Estabilidade Hidrodinâmica de Discos de Acreção

Enunciado do Problema
A instabilidade centrífuga em fluxos rotativos é tradicionalmente governada pelo critério de Rayleigh, que estabelece que a instabilidade surge quando o momento angular específico diminui com o raio. Este critério foi derivado para fluidos incompressíveis e não viscosos. No entanto, análises recentes de rotação relativística compressível suportada por pressão (Gourgouliatos & Komissarov 2018b) sugeriram uma generalização deste critério onde a instabilidade depende explicitamente do número de Mach do fluxo ($M$). No limite Newtoniano, este "critério KG" postula que um fluxo é instável se $d \ln \Psi / d \ln R < M^2$, onde $\Psi = \rho (v_\phi R)^2$. Isso implica que mesmo fluxos estáveis sob o critério clássico de Rayleigh (onde $d(v_\phi R)^2/dR > 0$) poderiam tornar-se instáveis se o número de Mach fosse suficientemente alto.

Esta previsão teórica desafia a compreensão convencional dos discos de acreção astrofísicos. Discos Keplerianos finos, que rotacionam em velocidades hipersônicas, são amplamente considerados hidrodinamicamente estáveis porque satisfazem o critério de Rayleigh. Os autores questionam se os altos números de Mach nestes sistemas, combinados com a compressibilidade, poderiam de fato torná-los suscetíveis à instabilidade centrífuga, um fator frequentemente negligenciado em análises de estabilidade padrão que assumem incompressibilidade ou ignoram a dependência explícita do número de Mach.

Metodologia
Para validar a versão Newtoniana do critério KG, os autores realizaram simulações numéricas axisimétricas de fluxos rotativos transônicos não relativísticos. As simulações resolveram as equações da dinâmica de fluidos compressíveis ideais (equações de continuidade, momento e energia) usando o código AMR-VAC com uma equação de estado de gás ideal ($\gamma = 5/3$).

A configuração da simulação envolveu:

• Equilíbrio Inicial: Uma configuração com densidade uniforme ($\rho = \rho_0$) e um perfil de rotação $v_\phi(R)$ que transiciona de rotação de corpo rígido em um núcleo central para um decaimento de lei de potência ($v_\phi \propto R^{-n}$) fora do núcleo.
• Equilíbrio de Pressão: A distribuição de pressão foi derivada para satisfazer a equação de balanço de força radial $dp/dR = \rho v_\phi^2 / R$.
• Perturbação: O equilíbrio foi perturbado com uma pequena perturbação de velocidade azimutal ($\delta v_\phi = \sin(kz)$) para desencadear potenciais instabilidades.
• Espaço de Parâmetros: Os autores variaram o índice da lei de potência $n$ (variando de 0,2 a 1,0) e a pressão central $p_0$ para criar modelos com diferentes números de Mach ($M_c$) na fronteira do núcleo. Isso permitiu testar configurações que são estáveis de acordo com o critério de Rayleigh original, mas instáveis de acordo com o critério KG.

Principais Contribuições e Resultados

1. Validação do Critério KG: Os resultados numéricos demonstraram concordância perfeita com o critério teórico KG. Modelos previstos como instáveis (onde $2(1-n) < M^2$) exibiram crescimento rápido de perturbações, enquanto modelos previstos como estáveis permaneceram assim. Isso confirma que a instabilidade centrífuga em fluxos compressíveis é fortemente dependente do número de Mach, uma característica ausente no limite incompressível clássico.
2. Evolução Não Linear: As simulações revelaram o desenvolvimento não linear da instabilidade. As perturbações cresceram mais rapidamente na fronteira interna do domínio instável. À medida que a instabilidade evoluía, o fluxo desenvolveu correntes radiais (dedos) movendo-se para longe do núcleo. Estas correntes interagiram com o fluido circundante para formar estruturas semelhantes a cogumelos e casulos que exibem rollos característicos da instabilidade de Kelvin-Helmholtz. Este processo levou a uma camada turbulenta com mistura forte, que se estendeu significativamente além da região inicialmente instável.
3. Reavaliação da Estabilidade de Discos de Acreção: Os autores estenderam sua análise heurística de estabilidade para incluir tanto forças de pressão quanto a gravidade de um objeto compacto central. Eles derivaram um critério de instabilidade geral que retém a dependência explícita do número de Mach.
   • No limite de gravidade nula, este critério reduz-se ao critério KG para rotação suportada por pressão.
   • No limite de rotação nula, reduz-se ao critério de Schwarzschild para instabilidade convectiva.
   • No caso geral, o critério derivado é matematicamente equivalente ao critério padrão de Solberg-Høiland, que não exibe explicitamente a dependência do número de Mach.

Significância e Alegações
O artigo afirma que a estabilidade hidrodinâmica de discos de acreção é frequentemente explicada usando o critério de Rayleigh original, mesmo quando a rotação é altamente supersônica. Ao derivar um critério que inclui explicitamente o número de Mach e mostrar sua equivalência ao critério de Solberg-Høiland, os autores esclarecem que a estabilidade dos discos Keplerianos finos não é meramente um resultado de satisfazer a condição incompressível de Rayleigh, mas sim uma consequência do equilíbrio específico entre forças centrífugas e convectivas descrito por Solberg-Høiland.

Os autores concluem que, embora o critério KG forneça uma formulação útil para fluxos suportados por pressão e jatos (onde a gravidade é desprezível), a dependência explícita do número de Mach é efetivamente absorvida nos gradientes de entropia e momento angular na formulação geral de Solberg-Høiland. Consequentemente, para fluxos isentrópicos (onde a entropia é constante), o critério reduz-se ao critério de Rayleigh original, consistente com a análise de Rayleigh (1917) para fluidos barotrópicos. O estudo confirma que altos números de Mach podem desestabilizar fluxos suportados por pressão que seriam de outra forma estáveis, um achado relevante para a dinâmica de jatos astrofísicos e a interação de ventos supersônicos com o plasma circundante, embora os autores notem que simulações 3D completas são necessárias para caracterizar totalmente a mistura não linear desencadeada por esta instabilidade.