Evolution of passive scalar mixing layers in stratified and unstratified homogeneous turbulence

Simulações de grandes vórtices de alta resolução revelam que, embora a mistura de escalares passivos em turbulência homogênea estratificada se comporte de modo semelhante aos casos não estratificados na direção transversal, a estratificação inibe severamente a mistura vertical ao suprimir a agitação em grande escala, levando a limites de crescimento e intensidades de flutuação distintos que informam abordagens de modelagem específicas para o fluxo de escalares.

O essencial

Imagine que você está observando uma gota de tinta se espalhando através de um copo de água. Se você mexer a água suavemente, a tinta se distribui uniformemente. É assim que os cientistas geralmente pensam sobre a mistura em fluidos: a turbulência atua como uma colher gigante, misturando tudo até ficar uniforme.

Mas o que acontece se essa água estiver em camadas? Imagine que a água no fundo é pesada e salgada, enquanto a água no topo é leve e doce. Isso é chamado de estratificação. No mundo real, isso ocorre no oceano (onde a água profunda é mais densa) e na atmosfera (onde o ar fica mais fino e leve à medida que você sobe).

Este artigo é um experimento de computador de alta tecnologia que pergunta: Como essa camada altera a maneira como uma "mancha" (como poluição ou fumaça) se espalha através de um fluido turbulento?

Aqui está a história do que eles descobriram, dividida em conceitos simples.

O Cenário: Dois Tipos de "Manchas"

Os pesquisadores criaram um fluido virtual que inicialmente girava de forma caótica (turbulento). Em seguida, introduziram duas "manchas" diferentes (escalares passivos) para ver como elas se comportavam:

1. A Mancha Horizontal: Uma camada de tinta espalhada lateralmente, como uma folha de papel plana flutuando na água.
2. A Mancha Vertical: Uma camada de tinta espalhada de cima para baixo, como uma parede vertical de cor.

Eles executaram duas simulações: uma em água normal (não estratificada) e outra em água "em camadas" (estratificada).

A Grande Descoberta: A Diferença "Cima-Baixo" vs. "Lado a Lado"

1. A Mancha Lado a Lado (Camada Transversal)

O que aconteceu: Quando a mancha foi espalhada lateralmente, a estratificação não impediu sua dispersão. Na verdade, ela se espalhou ligeiramente mais rápido na água estratificada do que na água normal.
A Analogia: Imagine uma multidão de pessoas correndo em um corredor. Se o chão for perfeitamente plano (não estratificado), elas correm em todas as direções. Se o chão tiver uma leve, invisível inclinação (estratificada), elas ainda correm lateralmente muito bem, talvez até com um pouco mais de energia. A "tinta" se espalha amplamente em ambos os casos.
O Pulo do Gato: Embora a dispersão geral fosse semelhante, a "tinta" na água estratificada era mais "pontuda". Em vez de um gradiente suave, ela tinha bordas mais afiadas e irregulares. Era mais "intermitente", o que significa que havia bolsões de tinta pura e bolsões de água pura, com menos meio-termo suave.

2. A Mancha Cima-Baixo (Camada Vertical)

O que aconteceu: É aqui que a mágica (e a restrição) aconteceu. Na água normal, a mancha vertical se espalhava para cima e para baixo facilmente, assim como a lateral. Mas na água estratificada, a dispersão parou quase completamente.
A Analogia: Imagine tentar mexer um milk-shake espesso com uma colher. Se você tentar mover a colher para cima e para baixo, as camadas do milk-shake resistem a você. As coisas pesadas querem ficar no fundo, e as coisas leves querem ficar no topo. O movimento de "mexer" é esmagado.
O Resultado: A mancha vertical cresceu um pouquinho no início, mas depois atingiu um "teto". Ela não pôde ficar mais larga porque as camadas estáveis do fluido atuaram como uma tampa, impedindo que a turbulência misturasse as coisas verticalmente. O fluido ainda podia girar lateralmente, mas não podia misturar para cima e para baixo.

Por Que Isso Importa? (O "Porquê" por Trás do "O Quê")

Os pesquisadores descobriram que, na direção vertical, o fluido se comporta como uma mola. Assim que a turbulência tenta empurrar uma camada pesada para cima ou uma camada leve para baixo, a gravidade a puxa de volta. Isso interrompe o movimento de "mexer".

No entanto, o fluido ainda pode girar lateralmente. Assim, o "comprimento vertical" da turbulência fica travado em um tamanho específico (determinado pela força da gravidade e pelas camadas do fluido), e a mancha não pode crescer além desse tamanho.

A "Receita" para Previsão

O artigo também tentou criar uma "receita" matemática simples para prever como essas manchas se espalhariam sem precisar de um supercomputador.

• Se você conhece a forma da mancha: Você pode usar uma fórmula simples de um número para prever a velocidade com que ela se espalha lateralmente. Funciona muito bem.
• Se você não conhece a forma: Você tem que adivinhar a forma (assumindo que ela se parece com uma curva suave). Se fizer isso, precisa de uma fórmula de dois números. Isso funciona muito bem depois que a mancha teve tempo de se estabilizar em um ritmo com o fluido giratório.

O Resumo Final

• Mistura lateral: Camadas estáveis (como no oceano profundo ou na atmosfera superior) não impedem a mistura lateral; elas podem até torná-la um pouco mais intensa e irregular.
• Mistura vertical: Camadas estáveis atuam como um freio. Elas impedem que o fluido misture para cima e para baixo quase completamente.
• A distinção entre "Mexer" e "Misturar": O fluido ainda pode "mexer" (se mover) lateralmente, mas não pode "misturar" (combinar) verticalmente porque as camadas resistem à troca.

Os autores observam que seu experimento usou um tipo específico de propriedade do fluido (número de Prandtl de 0,7). Eles alertam que, se o fluido fosse "mais espesso" ou tivesse propriedades diferentes (número de Prandtl > 1), os resultados poderiam mudar devido a um efeito "inverso" onde a mistura cria sua própria flutuabilidade. Mas, para as condições que eles testaram, a regra "lateral é livre, cima-baixo é bloqueado" permanece verdadeira.

Resumo técnico

Resumo Técnico: Evolução de Camadas de Mistura de Escalares Passivos em Turbulência Homogênea Estratificada e Não Estratificada

Enunciação do Problema
A dispersão de partículas e escalares em ambientes estratificados estáveis, como a alta troposfera, a estratosfera e as camadas mistas oceânicas, constitui um desafio crítico nas ciências atmosféricas e oceânicas. Embora os efeitos da estratificação estável sobre o campo de densidade e o decaimento da turbulência sejam bem documentados, a dispersão de escalares passivos (por exemplo, aerossóis, poluentes, microplásticos) permanece menos compreendida. Especificamente, não está claro como uma camada de mistura de escalar passivo (PSML) evolui quando introduzida em turbulência homogênea isotrópica (HIT) que subsequentemente decai sob a influência de um gradiente de densidade estabilizante. A questão central é como a introdução de escalas de flutuabilidade e anisotropia altera a evolução auto-similar das estatísticas escalares observadas na HIT não estratificada.

Metodologia
Os autores empregaram simulações de grandes vórtices (LES) de alta resolução para investigar o decaimento tanto da turbulência homogênea não estratificada quanto da estratificada estável. As simulações utilizaram as equações de Navier-Stokes não divergentes com a aproximação de Boussinesq não hidrostática.

• Configuração do Escoamento: Duas simulações foram conduzidas. A primeira foi HIT não estratificada. A segunda foi idêntica nas condições iniciais, mas incluiu um gradiente de densidade estabilizante uniforme ($g \neq 0$), fazendo com que o escoamento transitasse para turbulência fortemente estratificada após aproximadamente um período de flutuabilidade.
• Escalares Passivos: Duas camadas de mistura de escalar passivo foram inicializadas: uma com um gradiente médio na direção vertical ($\xi_z$) e outra na direção transversal ($\xi_y$). Estas foram modeladas como frações de mistura ($0 \le \xi \le 1$).
• Detalhes Numéricos: Para minimizar efeitos difusivos que causam desvios da auto-similaridade, as simulações utilizaram hiperviscosidade e hiperdifusividade com coeficientes iguais. O número de Prandtl molecular foi definido como $Pr = 0.7$. Os números de Reynolds foram altos ($Re_h > 15.000$ e $Re_\lambda > 450$ para o caso não estratificado), garantindo a existência de intervalos inercial-convectivos e inercial-difusivos. A resolução foi suficiente para resolver a escala de comprimento de flutuabilidade ($L_b/\Delta \gg 1$).
• Análise: O estudo acompanhou a evolução das larguras dos perfis escalares, variância, momentos de ordem superior (assimetria, curtose), fluxos e eficiência de mistura. Também avaliou modelos de difusividade turbulenta para o fluxo escalar transversal sob duas condições: quando o perfil escalar médio é conhecido e quando deve ser assumido.

Principais Resultados

1. Mistura Transversal ($\xi_y$): Na direção perpendicular à gravidade, a evolução do escalar passivo no caso estratificado é amplamente semelhante ao caso não estratificado, embora se espalhe ligeiramente mais rápido.

   • As flutuações escalares são mais intensas no caso estratificado, com variância de pico aproximadamente 15% maior.
   • A interface turbulenta/não turbulenta é mais intermitente no caso estratificado.
   • Uma vez que o campo escalar atinge o quase-equilíbrio com o campo de velocidade, os perfis de fluxo normalizados para ambos os casos estratificado e não estratificado colapsam, indicando mecanismos de transporte subjacentes semelhantes apesar da estratificação.

2. Mistura Vertical ($\xi_z$): A estratificação inibe drasticamente a mistura vertical.

   • A largura da camada vertical ($\delta_z$) cresce inicialmente, mas rapidamente deixa de expandir após aproximadamente um período de flutuabilidade ($t^* \approx 3$).
   • Esta cessação do crescimento é atribuída à supressão da agitação vertical em grande escala pelas forças de flutuabilidade.
   • A largura final da camada é proporcional à escala de comprimento integral vertical da velocidade horizontal ($L_v$), que é restringida para manter um número de Froude vertical da ordem de um ($Fr_v \sim O(1)$).
   • Embora o perfil médio mostre crescimento monótono devido à difusão molecular, a análise de fluxo revela um fenômeno de "fluxo reverso" onde a agitação ocasionalmente varre os tentáculos escalares de volta para o seu lado original, consistente com um comportamento semelhante a ondas internas.

3. Eficiência de Mistura: A mistura cumulativa (dissipação da variância escalar) no caso estratificado eventualmente excede a do caso não estratificado em tempos tardios. Isso é atribuído ao fato de que a mistura segue a agitação, e a taxa de agitação horizontal no escoamento estratificado é proporcional à escala de comprimento horizontal, que cresce maior do que no caso não estratificado.

4. Modelagem de Fluxo:

   • Perfil Conhecido: Se o perfil escalar médio é conhecido, um modelo de difusividade turbulenta de constante única ($F_y = -c_1 \delta_y v' \partial \xi / \partial y$) prevê efetivamente o fluxo para ambos os casos estratificado e não estratificado, com $c_1 \approx 0.4$.
   • Perfil Assumido: Se a forma do perfil (assumida como uma função erro) e a escala de comprimento escalar devem ser estimadas a partir das estatísticas de velocidade, um modelo de duas constantes é necessário. Este modelo desempenha bem em tempos tardios quando o escalar está em quase-equilíbrio com o campo de velocidade, permitindo que a escala de comprimento escalar seja escalonada a partir da taxa de dissipação de energia cinética.

Significado e Alegações
O artigo afirma fornecer uma caracterização detalhada da mistura de escalares passivos em turbulência estratificada decaindo, preenchendo a lacuna entre previsões teóricas de difusão anisotrópica e observações numéricas.

• Anisotropia: O estudo confirma que a estratificação cria uma dicotomia nítida na mistura: a mistura horizontal permanece ativa e até mesmo aprimorada, enquanto a mistura vertical é efetivamente suprimida uma vez que o escoamento se torna fortemente estratificado.
• Utilidade de Modelagem: Os autores demonstram que abordagens padrão de difusividade turbulenta permanecem viáveis para escoamentos estratificados, desde que o escalar esteja em quase-equilíbrio com o campo de velocidade. Eles destacam que assumir um perfil conhecido simplifica significativamente os requisitos de modelagem em comparação com assumir uma forma de perfil.
• Limitações e Trabalho Futuro: Os autores modestamente observam que seus resultados são específicos para um número de Prandtl de 0,7. Eles afirmam explicitamente que a literatura recente sugere que números de Prandtl mais altos ($Pr > 1$) induzem fluxos de flutuabilidade reversos em pequenas escalas, o que poderia alterar as taxas de decaimento do escoamento e, consequentemente, as taxas de mistura de escalares passivos. Portanto, as descobertas atuais não se aplicam diretamente a cenários de alto número de Prandtl sem investigação adicional.

O estudo conclui que, embora a estratificação altere fundamentalmente a dinâmica vertical de uma camada de mistura, a mistura transversal retém muitas características da turbulência não estratificada, embora com maior intermitência e intensidade de flutuação.