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

Resumen Técnico: Inestabilidad Centrífuga de Flujos Compresibles y la Estabilidad Hidrodinámica de los Discos de Acreción

Planteamiento del Problema
La inestabilidad centrífuga en flujos rotatorios está gobernada tradicionalmente por el criterio de Rayleigh, el cual establece que la inestabilidad surge cuando el momento angular específico disminuye con el radio. Este criterio fue derivado para fluidos incompresibles e inviscosos. Sin embargo, análisis recientes de la rotación relativista compresible soportada por presión (Gourgouliatos & Komissarov 2018b) sugirieron una generalización de este criterio donde la inestabilidad depende explícitamente del número de Mach del flujo ($M$). En el límite Newtoniano, este "criterio KG" postula que un flujo es inestable si $d \ln \Psi / d \ln R < M^2$, donde $\Psi = \rho (v_\phi R)^2$. Esto implica que incluso los flujos estables bajo el criterio clásico de Rayleigh (donde $d(v_\phi R)^2/dR > 0$) podrían volverse inestables si el número de Mach es lo suficientemente alto.

Esta predicción teórica desafía la comprensión convencional de los discos de acreción astrofísicos. Los discos Keplerianos delgados, que rotan a velocidades hipersónicas, son ampliamente considerados hidrodinámicamente estables porque satisfacen el criterio de Rayleigh. Los autores cuestionan si los altos números de Mach en estos sistemas, combinados con la compresibilidad, podrían en realidad hacerlos susceptibles a la inestabilidad centrífuga, un factor que a menudo se pasa por alto en los análisis de estabilidad estándar que asumen incompresibilidad o ignoran la dependencia explícita del número de Mach.

Metodología
Para validar la versión Newtoniana del criterio KG, los autores realizaron simulaciones numéricas axisimétricas de flujos rotatorios transónicos no relativistas. Las simulaciones resolvieron las ecuaciones de la dinámica de fluidos compresibles ideales (ecuaciones de continuidad, momento y energía) utilizando el código AMR-VAC con una ecuación de estado de gas ideal ($\gamma = 5/3$).

La configuración de la simulación involucró:

• Equilibrio Inicial: Una configuración de densidad uniforme ($\rho = \rho_0$) y un perfil de rotación $v_\phi(R)$ que transiciona de una rotación de cuerpo rígido en un núcleo central a un decaimiento de ley de potencia ($v_\phi \propto R^{-n}$) fuera del núcleo.
• Balance de Presión: La distribución de presión se derivó para satisfacer la ecuación de balance de fuerza radial $dp/dR = \rho v_\phi^2 / R$.
• Perturbación: El equilibrio fue perturbado con una pequeña perturbación de velocidad azimutal ($\delta v_\phi = \sin(kz)$) para activar posibles inestabilidades.
• Espacio de Parámetros: Los autores variaron el índice de la ley de potencia $n$ (rango de 0.2 a 1.0) y la presión central $p_0$ para crear modelos con diferentes números de Mach ($M_c$) en el límite del núcleo. Esto permitió probar configuraciones que son estables según el criterio original de Rayleigh pero inestables según el criterio KG.

Contribuciones Clave y Resultados

1. Validación del Criterio KG: Los resultados numéricos demostraron una concordancia perfecta con el criterio teórico KG. Los modelos predichos como inestables (donde $2(1-n) < M^2$) exhibieron un crecimiento rápido de las perturbaciones, mientras que los modelos predichos como estables permanecieron así. Esto confirma que la inestabilidad centrífuga en flujos compresibles depende fuertemente del número de Mach, una característica ausente en el límite incompresible clásico.
2. Evolución No Lineal: Las simulaciones revelaron el desarrollo no lineal de la inestabilidad. Las perturbaciones crecieron más rápido en el límite interno del dominio inestable. A medida que la inestabilidad evolucionaba, el flujo desarrolló corrientes radiales (dedos) que se alejaban del núcleo. Estas corrientes interactuaron con el fluido circundante para formar estructuras similares a hongos y capullos que exhiben rollos característicos de la inestabilidad de Kelvin-Helmholtz. Este proceso condujo a una capa turbulenta con una fuerte mezcla, la cual se extendió significativamente más allá de la región inicialmente inestable.
3. Reevaluación de la Estabilidad de los Discos de Acreción: Los autores extendieron su análisis heurístico de estabilidad para incluir tanto las fuerzas de presión como la gravedad de un objeto compacto central. Derivaron un criterio de inestabilidad general que retiene la dependencia explícita en el número de Mach.
   • En el límite de gravedad evanescente, este criterio se reduce al criterio KG para la rotación soportada por presión.
   • En el límite de rotación evanescente, se reduce al criterio de Schwarzschild para la inestabilidad convectiva.
   • En el caso general, el criterio derivado es matemáticamente equivalente al estándar criterio de Solberg-Høiland, el cual no muestra explícitamente la dependencia del número de Mach.

Significancia y Reivindicaciones
El artículo afirma que la estabilidad hidrodinámica de los discos de acreción se explica a menudo utilizando el criterio original de Rayleigh, incluso cuando la rotación es altamente supersónica. Al derivar un criterio que incluye explícitamente el número de Mach y mostrar su equivalencia con el criterio de Solberg-Høiland, los autores aclaran que la estabilidad de los discos Keplerianos delgados no es simplemente un resultado de satisfacer la condición incompresible de Rayleigh, sino una consecuencia del balance específico entre las fuerzas centrífugas y convectivas descritas por Solberg-Høiland.

Los autores concluyen que, si bien el criterio KG proporciona una formulación útil para flujos soportados por presión y chorros (jets) (donde la gravedad es despreciable), la dependencia explícita del número de Mach es efectivamente absorbida por los gradientes de entropía y momento angular en la formulación general de Solberg-Høiland. Por consiguiente, para flujos isentrópicos (donde la entropía es constante), el criterio se reduce al criterio original de Rayleigh, consistente con el análisis de Rayleigh (1917) para fluidos barotrópicos. El estudio confirma que los altos números de Mach pueden desestabilizar flujos soportados por presión que de otro modo serían estables, un hallazgo relevante para la dinámica de los chorros astrofísicos y la interacción de vientos supersónicos con el plasma circundante, aunque los autores señalan que se requieren simulaciones completas en 3D para caracterizar plenamente la mezcla no lineal desencadenada por esta inestabilidad.