The (Limited) Effect of Viscosity in Multiphase Turbulent Mixing

Enunciado del Problema
El gas multifásico es ubicuo en entornos astrofísicos como el medio interestelar (ISM), el medio circumgaláctico (CGM) y el medio intracúmulo (ICM). La dinámica de estos sistemas está gobernada por la interacción entre la turbulencia, el enfriamiento radiativo y los procesos de transporte. Si bien se ha estudiado el papel de los campos magnéticos en la supresión de la inestabilidad de Kelvin-Helmholtz (KHI) y en la alteración de la mezcla en capas de mezcla radiativa turbulenta (TRMLs), el impacto específico de la viscosidad sigue siendo objeto de debate. Se espera que los plasmas totalmente ionizados sean altamente viscosos (Spitzer 1962), sin embargo, la viscosidad efectiva neta en sistemas astrofísicos es incierta. Una pregunta crítica es si la viscosidad suprime la turbulencia y la mezcla necesarias para que se cumplan las relaciones de escala "universales" entre la velocidad turbulenta y la tasa de enfriamiento (Fielding et al. 2020; Tan et al. 2021), particularmente en regímenes donde el número de Karlovitz ($Ka < 1$) sugiere flujo laminar.

Metodología
Los autores emplearon simulaciones numéricas idealizadas utilizando el código Athena para aislar los efectos de la viscosidad. El estudio procedió en dos etapas principales:

1. Simulaciones 2D Adiabáticas: Para determinar la viscosidad crítica ($\eta_{Crit}$) requerida para suprimir la KHI, los autores simularon flujos de cizalla con sobredensidades ($\chi = \rho_{frío}/\rho_{caliente}$) y números de Mach ($M$) variables. Utilizaron una configuración de placa plana para probar la dependencia de la sobredensidad y una configuración de lámina plana para probar la dependencia del número de Mach, evitando artefactos de frontera. Derivaron una estimación analítica para $\eta_{Crit}$ basada en la comparación de la escala de tiempo viscosa y la escala de tiempo de crecimiento de la KHI, utilizando una media armónica para la viscosidad cinemática para tener en cuenta las diferencias de momento en fluidos de múltiples densidades.
2. Capas de Mezcla Radiativa Turbulenta 3D (TRMLs): Utilizando la viscosidad crítica derivada para una sobredensidad de $\chi=100$, los autores realizaron simulaciones 3D de TRMLs con viscosidad explícita. Variaron la fracción de viscosidad relativa a $\eta_{Crit}$ (desde 0% hasta 50%) y el número de Damköhler ($Da = t_{turb}/t_{enfriamiento}$) ajustando la longitud de onda de la perturbación. Esto les permitió explorar tanto regímenes de enfriamiento débil ($Da < 1$) como de enfriamiento fuerte ($Da > 1$). Las simulaciones incluyeron enfriamiento radiativo (algoritmo de Townsend) pero ninguna conducción térmica explícita, asegurando $Ka < 1$.

Contribuciones y Resultados Clave

• Escala de la Viscosidad Crítica: El estudio confirma que la viscosidad crítica requerida para suprimir la KHI escala linealmente con la sobredensidad ($\eta_{Crit} \propto \chi$) en el límite de alta sobredensidad. Esto se debe a que el medio frío y denso transporta la mayor parte del momento, requiriendo una fricción más fuerte para suprimir la inestabilidad. Los autores también verificaron que a altos números de Mach, los efectos de compresibilidad estabilizan la KHI incluso sin viscosidad, definiendo un número de Mach crítico ($M_{Crit}$) más allá del cual la inestabilidad es suprimida.
• Régimen de Enfriamiento Débil ($Da < 1$): En este régimen, las escalas de tiempo de enfriamiento son más largas que las de mezcla turbulenta. Los autores encontraron que la viscosidad suprime fuertemente la turbulencia, transicionando el flujo de turbulento a laminar. En consecuencia, la relación universal entre la velocidad turbulenta ($u'$) y el brillo superficial ($Q$) derivada para fluidos no viscosos se rompe. Sin embargo, contra intuitivamente, la luminosidad bolométrica total permanece en gran medida inalterada. Los autores atribuyen esto a un "efecto de cuello de botella" donde el largo tiempo de enfriamiento limita la pérdida total de energía independientemente del mecanismo de mezcla (difusión vs. convección).
• Régimen de Enfriamiento Fuerte ($Da > 1$): En este régimen, las escalas de tiempo de enfriamiento son más cortas que las escalas de tiempo viscosas. Los autores descubrieron que la viscosidad tiene un efecto limitado en la dinámica del sistema. A pesar de la presencia de viscosidad, el proceso de enfriamiento afina los gradientes de temperatura y velocidad en la interfaz, contrarrestando efectivamente el efecto de suavizado de la viscosidad. Esto permite que la KHI se desarrolle y mantenga niveles de turbulencia similares a los casos no viscosos. Como resultado, las relaciones de escala entre $u'$ y $Q$ permanecen intactas, y el sistema se comporta efectivamente como no viscoso.
• Competencia entre Viscosidad y Enfriamiento: Los autores derivaron una condición analítica ($Da > \zeta$, donde $\zeta$ es la fracción de viscosidad crítica) para cuando el enfriamiento domina sobre la viscosidad. Sus resultados numéricos se alinean con esta predicción. En el régimen de enfriamiento fuerte, la mayor cantidad de gas a temperatura intermedia en las corridas viscosas (debido a un frente ligeramente más ancho) conduce a un aumento modesto en el brillo superficial (hasta un factor de $\sim 2$) en comparación con las corridas no viscosas, pero las leyes de escala fundamentales se mantienen.

Significado e Implicaciones
El artículo afirma que el papel de la viscosidad en la mezcla turbulenta multifásica es más limitado de lo que se esperaba ingenuamente, particularmente en el régimen de enfriamiento fuerte que caracteriza muchos entornos astrofísicos (por ejemplo, flujos de salida galácticos, CGM, ICM).

• Diagnósticos Observacionales: Los hallazgos sugieren que las predicciones teóricas establecidas que relacionan los anchos de línea (turbulencia) con la emisividad (enfriamiento) aún pueden aplicarse a datos observacionales, incluso en sistemas donde la viscosidad podría ser significativa, siempre que el sistema se encuentre en el régimen de enfriamiento fuerte.
• Modelos de Subrejilla: Los resultados proporcionan una justificación física para los actuales modelos de subrejilla multifásicos utilizados en simulaciones cosmológicas a gran escala. Estos modelos a menudo se basan en simulaciones a pequeña escala no viscosas para informar las tasas de transferencia de masa y las prescripciones de enfriamiento. Los autores demuestran que estas prescripciones siguen siendo válidas incluso cuando está presente la viscosidad, ya que el régimen dominado por el enfriamiento niega efectivamente la supresión viscosa de la turbulencia y la mezcla.
• Contexto Astrofísico: El estudio destaca que, si bien la viscosidad puede influir en la morfología a gran escala (por ejemplo, en colas despojadas por presión de arrastre), no altera fundamentalmente las tasas de disipación de energía o transferencia de masa en el límite de enfriamiento rápido. Esto ofrece confianza en la interpretación de la evolución del gas multifásico en halos galácticos y cúmulos.

Los autores señalan advertencias respecto a su configuración idealizada, específicamente el uso de viscosidad constante en lugar de la viscosidad dependiente de la temperatura de Spitzer, la ausencia de campos magnéticos (que inducen viscosidad anisotrópica) y la falta de conducción térmica explícita. Sugieren estos aspectos como vías para trabajos futuros que refinen aún más la comprensión de los procesos de transporte en plasmas astrofísicos.