No evidence of vorticity production from initially irrotational turbulent gravitational collapse

Resumen Técnico: Ausencia de evidencia de producción de vorticidad a partir de un colapso gravitatorio turbulento inicialmente irrotacional

Planteamiento del problema
Se cree ampliamente que el colapso gravitatorio genera turbulencia, la cual podría potencialmente sembrar un dínamo de pequeña escala para amplificar campos magnéticos por debajo de la escala de Jeans. Sin embargo, los dínamos de pequeña escala requieren fundamentalmente turbulencia vortical (rotacional), mientras que el colapso gravitatorio produce predominantemente flujos irrotacionales (potenciales). En ausencia de campos magnéticos o fuerzas de excitación externas, y empleando una ecuación de estado barotrópica (donde los gradientes de presión y densidad son paralelos), el término baroclínico no puede generar vorticidad. Por consiguiente, el único mecanismo disponible para producir vorticidad a partir de un flujo inicialmente irrotacional es a través de la interacción con la viscosidad. Estudios previos han sugerido que la compresión inducida por el colapso podría impulsar la acción de un dínamo, pero no está claro si el flujo de colapso genera la vorticidad necesaria o si los efectos de dínamo observados son simplemente el decaimiento de la turbulencia inicial preexistente.

Metodología
Los autores realizan simulaciones numéricas directas (DNS) de colapso gravitatorio utilizando el código Pencil, empleando diferencias centradas de sexto orden y un esquema de paso de tiempo de tercer orden. El estudio utiliza un dominio cúbico periódico con una ecuación de estado isentrópica ($\gamma = 5/3$) y un perfil de densidad de una esfera de Bonnor-Ebert modificada de tipo isotérmico.

Las decisiones metodológicas clave incluyen:

• Condiciones Iniciales: El campo de velocidad se construye para ser estrictamente irrotacional en el espacio de Fourier, con fases aleatorias y una pendiente espectral específica. No se incluyen campos magnéticos para aislar los mecanismos de generación de vorticidad.
• Viscosidad: En lugar de utilizar viscosidad de choque (que puede oscurecer los resultados), los autores utilizan una viscosidad cinemática constante ($\nu$).
• Resolución y Parámetros: Las simulaciones se ejecutan en resoluciones que van desde $512^3$ hasta $2048^3$ puntos de malla. Los autores varían la amplitud de la velocidad inicial ($u_{ini}$) y la viscosidad ($\nu$) para determinar los límites numéricamente fiables.
• Análisis: El estudio rastrea la evolución de la entropía de giro ($\langle \omega^2 \rangle/2$) y descompone su derivada temporal en tres términos: un término tipo dínamo ($T_{dyn}$), un término de generación impulsado por la viscosidad y los gradientes de densidad ($T_{gen}$), y un término disipativo ($T_{dis}$). Los autores también analizan los espectros de velocidad, densidad y vorticidad, y examinan el escalamiento de la vorticidad con la amplitud del flujo inicial.

Resultos Clave

1. Dependencia de la Viscosidad: Se encuentra que la producción de vorticidad es un efecto secundario impulsado por la amplificación de las divergencias del flujo inicial mediante la viscosidad. El estudio halla que la relación entre la generación de vorticidad y la disipación ($R_{gen}$) no es proporcional a la viscosidad; más bien, la generación de vorticidad es más significativa durante la fase temprana, antes del colapso principal, y no como una consecuencia del flujo de colapso en sí mismo.
2. Escalamiento con las Condiciones Iniciales: La vorticidad cuadrática media ($\omega_{rms}$) escala con el cuadrado de la amplitud de la velocidad inicial ($u_{ini}^2$), mientras que la divergencia de la velocidad cuadrática media escala linealmente con $u_{ini}$. Este escalamiento persiste durante toda la simulación, lo que sugiere que el colapso gravitatorio no impulsa directamente la producción de vorticidad. En su lugar, el colapso amplifica las divergencias del flujo inicial, las cuales la viscosidad convierte posteriormente en vorticidad.
3. Fiabilidad Numérica: Los autores establecen que los resultados dejan de ser fiables cuando el número de Reynolds de malla ($Re_{mesh}$) supera valores entre 0.7 y 2.0, dependiendo de la ejecución específica. En el régimen fiable, el término tipo dínamo ($T_{dyn}$) permanece subdominante, contribuyendo con menos del 20% a la evolución de la entropía de giro.
4. Análisis Espectral: El espectro de vorticidad exhibe un rango subinercial de $k^4$, análogo al comportamiento de los campos magnéticos en MHD, pero los autores lo atribuyen a la inestabilidad de Jeans impulsando fluctuaciones de velocidad aleatorias, en lugar de a un proceso de dínamo. El espectro de velocidad muestra un rango inercial pronunciado sin un efecto de "cuello de botella" (bottleneck).
5. Vorticidad Crítica para el Dínamo: Revisando el número de Reynolds de vorticidad magnética crítica ($Pr_M Re_\omega$) para el inicio del dínamo, los autores confirman una bifurcación alrededor de un valor de 300. Esto implica que los flujos puramente irrotacionales, que no pueden sostener suficiente vorticidad sin excitación externa o condiciones iniciales específicas, no actuarían como dínamos.

Significancia y Afirmaciones
El artículo concluye que el colapso gravitatorio, bajo las condiciones de una ecuación de estado barotrópica y una turbulencia inicialmente irrotacional, no produce vorticidad de manera eficiente. La vorticidad observada es un efecto secundario resultante de la conversión viscosa de las divergencias del flujo inicial, no una consecuencia directa de la dinámica del colapso.

Los autores afirman que sus hallazgos desafían las afirmaciones previas sobre la acción de dínamos de pequeña escala inducidos por el colapso. Sugieren que las observaciones anteriores de la amplificación de campos magnéticos en flujos de colapso fueron probablemente alimentadas por el decaimiento de la turbulencia inicial en lugar del colapso mismo. El estudio enfatiza la necesidad de realizar simulaciones numéricas directas sin modelos de sub-malla para evitar la generación artificial de vorticidad. El trabajo implica que lograr un dínamo autosostenible en tales escenarios de colapso es difícil y que el flujo de colapso por sí solo es insuficiente para generar la turbulencia vortical requerida para la acción de un dínamo.