No evidence of vorticity production from initially irrotational turbulent gravitational collapse

Resumo Técnico: Ausência de Evidência de Produção de Vorticidade a partir de Colapso Gravitacional Turbulento Inicialmente Irrotacional

Problema
Acredita-se amplamente que o colapso gravitacional gere turbulência, o que poderia potencialmente semear um dínamo de pequena escala para amplificar campos magnéticos abaixo da escala Jeans. No entanto, os dínamos de pequena escala requerem fundamentalmente turbulência vortical (rotacional), ao passo que o colapso gravitacional produz predominantemente fluxos irrotacionais (potenciais). Na ausência de campos magnéticos ou forças externas, e ao empregar uma equação de estado barotrópica (onde os gradientes de pressão e densidade são paralelos), o termo baroclínico não pode gerar vorticidade. Consequentemente, o único mecanismo disponível para produzir vorticidade a partir de um fluxo inicialmente irrotacional é através da interação com a viscosidade. Estudos anteriores sugeriram que o colapso induzido pela compressão poderia impulsionar a ação de dínamo, mas permanece incerto se o próprio fluxo de colapso gera a vorticidade necessária ou se os efeitos de dínamo observados são meramente o decaimento de uma turbulência inicial pré-existente.

Metodologia
Os autores realizam simulações numéricas diretas (DNS) de colapso gravitacional utilizando o Pencil Code, empregando diferenças centradas de sexta ordem e um esquema de integração temporal de terceira ordem. O estudo utiliza um domínio cúbico periódico com uma equação de estado isentrópica ($\gamma = 5/3$) e um perfil de densidade inicial de uma esfera de Bonnor-Ebert modificada e isotérmica.

Principais escolhas metodológicas incluem:

• Condições Iniciais: O campo de velocidade é construído para ser estritamente irrotacional no espaço de Fourier, com fases aleatórias e uma inclinação espectral específica. Nenhum campo magnético é incluído para isolar os mecanismos de geração de vorticidade.
• Viscosidade: Em vez de viscosidade de choque (que pode obscurecer os resultados), os autores utilizam uma viscosidade cinemática constante ($\nu$).
• Resolução e Parâmetros: As simulações são executadas em resoluções variando de $512^3$ a $2048^3$ pontos de malha. Os autores variam a amplitude da velocidade inicial ($u_{ini}$) e a viscosidade ($\nu$) para determinar limites numericamente confiáveis.
• Análise: O estudo monitora a evolução da enstrofia ($\langle \omega^2 \rangle/2$) e decompõe sua derivada temporal em três termos: um termo do tipo dínamo ($T_{dyn}$), um termo de geração impulsionado pela viscosidade e gradientes de densidade ($T_{gen}$), e um termo dissipativo ($T_{dis}$). Os autores também analisam os espectros de velocidade, densidade e vorticidade e examinam o escalonamento da vorticidade com a amplitude do fluxo inicial.

Resultos Principais

1. Dependência da Viscosidade: A produção de vorticidade é identificada como um efeito secundário impulsionado pela amplificação das divergências do fluxo inicial via viscosidade. O estudo constata que a razão entre a geração de vorticidade e a dissipação ($R_{gen}$) não é proporcional à viscosidade; em vez disso, a geração de vorticidade é mais significativa durante a fase inicial, antes do colapso principal, e não como uma consequência do próprio fluxo de colapso.
2. Escalonamento com Condições Iniciais: A vorticidade raiz-média do quadrado ($\omega_{rms}$) escala com o quadrado da amplitude da velocidade inicial ($u_{ini}^2$), enquanto a divergência da velocidade raiz-média escala linearmente com $u_{ini}$. Este escalonamento persiste ao longo da simulação, sugerindo que o colapso gravitacional não impulsiona diretamente a produção de vorticidade. Em vez disso, o colapso amplifica as divergências do fluxo inicial, as quais a viscosidade então converte em vorticidade.
3. Confiabilidade Numérica: Os autores estabelecem que os resultados tornam-se não confiáveis quando o número de Reynolds de malha ($Re_{mesh}$) excede valores entre 0,7 e 2,0, dependendo da execução específica. No regime confiável, o termo do tipo dínamo ($T_{dyn}$) permanece subdominante, contribuindo com menos de 20% para a evolução da enstrofia.
4. Análise Espectral: O espectro de vorticidade exibe um intervalo subinercial de $k^4$, análogo ao comportamento do campo magnético em MHD, mas os autores atribuem isso à instabilidade de Jeans impulsionando flutuações de velocidade aleatórias, em vez de um processo de dínamo. O espectro de velocidade mostra um intervalo inercial íngreme sem um efeito de "gargalo" (bottleneck).
5. Vorticidade Crítica para Dínamo: Revisando o número de Reynolds de vorticidade magnética crítica ($Pr_M Re_\omega$) para o início do dínamo, os autores confirmam uma bifurcação em torno de um valor de 300. Isso implica que fluxos puramente irrotacionais, que não podem sustentar vorticidade suficiente sem condução externa ou condições iniciais específicas, não atuariam como dínamos.

Significância e Alegações
O artigo conclui que o colapso gravitacional, sob as condições de uma equação de estado barotrópica e turbulência inicialmente irrotacional, não produz vorticidade de forma eficiente. A vorticidade observada é um efeito secundário resultante da conversão viscosa das divergências do fluxo inicial, não uma consequência direta da dinâmica do colapso.

Os autores afirmam que suas descobertas desafiam alegações anteriores de ação de dínamo de pequena escala induzida pelo colapso. Eles sugerem que as observações anteriores de amplificação de campo magnético em fluxos de colapso foram provavelmente alimentadas pelo decaimento da turbulência inicial, e não pelo colapso em si. O estudo enfatiza a necessidade de simulações numéricas diretas sem modelagem de subgrade para evitar a geração artificial de vorticidade. O trabalho implica que alcançar um dínamo autossustentável em tais cenários de colapso é difícil e que o fluxo de colapso sozinho é insuficiente para gerar a turbulência vortical necessária para a ação de dínamo.