Two pathways to diapycnal mixing in strongly stratified flows with no initial vertical shear

Este estudo combina a teoria linear e simulações numéricas diretas para revelar que, em fluxos fortemente estratificados sem cisalhamento vertical inicial, instabilidades de cisalhamento horizontal inevitavelmente impulsionam a mistura diapicnal através de dois caminhos distintos — seja via emergência direta de cisalhamento vertical ou através de uma evolução não linear para vórtices colunares — ambos os quais, em última análise, desencadeiam instabilidades de Kelvin-Helmholtz de pequena escala, mas produzem diferentes eficiências de mistura devido à excitação de distintas escalas verticais.

O essencial

Imagine o oceano ou a atmosfera como um bolo gigante e em camadas. As camadas são feitas de fluidos com diferentes densidades (como diferentes sabores de bolo) e não gostam de se misturar facilmente. Normalmente, os cientistas estudam o que acontece quando você empurra essas camadas lateralmente e para cima ou para baixo ao mesmo tempo (cisalhamento vertical). Mas este artigo faz uma pergunta diferente: o que acontece se você empurrar as camadas apenas lateralmente, sem movimento inicial para cima e para baixo, em um ambiente fortemente estratificado?

Os pesquisadores descobriram que, mesmo que você comece com um fluxo perfeitamente plano e horizontal, a natureza possui duas "receitas" ou caminhos distintos para criar o caos (turbulência) e misturar as camadas. A receita que a natureza escolhe depende inteiramente de como você "semeia" o experimento — essencialmente, qual pequeno empurrãozinho você dá ao fluido no início.

Aqui está a divisão dos dois caminhos usando analogias simples:

A Configuração: O Rio Calmo

Imagine um rio largo e calmo fluindo horizontalmente. A água é estratificada em camadas como uma pilha de panquecas (estratificação forte). No início, o fluxo é suave e bidimensional (ele só se move para a esquerda e para a direita, não para cima e para baixo).

Caminho 1: O Efeito "Sala Lotada" (A Rota Direta)

Como começa: Você dá ao rio um pequeno empurrão aleatório em todos os lugares ao mesmo tempo (como jogar um punhado de confetes no ar).
O que acontece:

1. A Ondulação: Devido às camadas, o fluido não apenas ondula para a esquerda e para a direita; ele imediatamente começa a ondular para cima e para baixo em muitos tamanhos diferentes ao mesmo tempo. Pense nisso como uma multidão de pessoas em uma sala, todas tentando se mover ao mesmo tempo, criando um sobressalto caótico e multidirecional.
2. O Cisalhamento: Essas ondulações criam fortes correntes verticais (cisalhamento) muito rapidamente.
3. A Ruptura: Essas correntes verticais tornam-se tão fortes que se rompem, criando pequenos redemoinhos violentos (como minúsculos tornados). Esta é a instabilidade de "Kelvin-Helmholtz", que se parece com as ondas quebrando que você vê quando o vento sopra sobre a água.
   O Resultado: A mistura ocorre de forma eficiente. Como a energia é espalhada por muitas diferentes escalas de ondulação, o "atrito" (dissipação viscosa) é menor, tornando o processo de mistura relativamente eficiente.

Caminho 2: O Efeito "Dança Sincronizada" (A Rota Indireta)

Como começa: Você dá ao rio um empurrão muito específico e organizado (como um maestro agitando uma batuta para fazer todos se moverem em um padrão específico).
O que acontece:

1. O Vórtice: Em vez de ondulações caóticas, o fluido se organiza em longas colunas verticais de água giratória (como gigantes tornados em pé no rio). Por muito tempo, o fluxo permanece perfeitamente bidimensional, consistindo apenas nessas grandes colunas giratórias.
2. O Bamboleio: Eventualmente, essas gigantescas colunas tornam-se instáveis. Elas começam a bambolear de uma forma muito específica e de alta frequência. Os pesquisadores chamam isso de "instabilidade hiperbólica". Imagine um pião que começa a balançar violentamente logo antes de cair.
3. A Ruptura: Esse bamboleio violento cria camadas verticais muito finas e nítidas de cisalhamento. Essas camadas finas então se rompem em pequenos redemoinhos violentos, exatamente como no Caminho 1.
   O Resultado: A mistura acontece, mas é menos eficiente. Por quê? Porque este caminho cria camadas verticais extremamente finas e nítidas. É necessário muita energia (atrito) para criar e quebrar essas camadas minúsculas e afiadas. É como tentar cortar um bloco grosso de queijo com uma faca cega (Caminho 1) versus uma lâmina de barbear (Caminho 2); a lâmina de barbear cria um corte muito mais nítido e intensivo em termos de energia.

A Grande Conclusão

O artigo prova que o cisalhamento vertical (movimento para cima e para baixo) não precisa estar presente desde o início. Ele é um subproduto inevitável do cisalhamento horizontal em fluidos fortemente estratificados, desde que o fluido seja espesso o suficiente (número de Reynolds alto).

• Se você começar com ruído aleatório: Você obtém o Caminho 1 (Mistura direta e eficiente).
• Se você começar com um padrão específico: Você obtém o Caminho 2 (Mistura indireta e menos eficiente).

Os pesquisadores usaram simulações computacionais poderosas para mostrar que esses dois caminhos são reais, distintos, e que a "receita" que você escolhe no início determina quanta energia é desperdiçada como calor versus quanto é usada para realmente misturar as camadas.

Em resumo: Mesmo em um fluido estratificado perfeitamente calmo, um empurrão horizontal eventualmente criará o caos vertical. Mas, dependendo de como você inicia o empurrão, esse caos será diferente e misturará as camadas com diferentes níveis de eficiência.

Resumo técnico

Enunciado do Problema
A estabilidade de fluxos horizontais de grande escala e fortemente estratificados para perturbações tridimensionais (3D) é um mecanismo crítico para a transferência descendente de energia e o transporte diapicnal em sistemas geofísicos e astrofísicos. Embora fluxos estratificados com cisalhamento vertical sejam bem estudados, seus equivalentes com cisalhamento horizontal (sem cisalhamento vertical inicial) receberam menos atenção. Uma questão central permanece: como fluxos de cisalhamento horizontal inicialmente bidimensionais (2D) e fortemente estratificados ($Fr \ll 1$) transitam para a turbulência e geram cisalhamento vertical? A literatura anterior sugere dois caminhos distintos: um envolvendo o crescimento direto de modos primários 3D que levam a instabilidades de Kelvin-Helmholtz (KH), e outro envolvendo a evolução não linear em vórtices colunares seguidos por instabilidades secundárias. No entanto, permanecia incerto se esses caminhos são mutuamente exclusivos devido a diferenças nas condições iniciais (ex: tangente hiperbólico vs. perfis sinusoidais) ou na formulação do problema (Problema de Valor Inicial vs. forçado por corpo), ou se representam rotas alternativas disponíveis dentro do mesmo sistema físico. Além disso, o impacto da seleção do caminho escolhido na eficiência da mistura diapicnal não foi sistematicamente quantificado.

Metodologia
Os autores combinam a teoria de estabilidade linear com Simulações Numéricas Diretas (DNS) para investigar essas transições em um domínio triperiodicamente periódico. O estudo foca em um fluxo de Kolmogorov plano, horizontal e verticalmente invariante ($\bar{u}(y) = \sin(y)\mathbf{e}_x$) em um fluido linearmente estratificado com frequência de flutuabilidade constante $N$. As simulações são conduzidas em altos números de Reynolds ($Re = 10.000$) e Péclet ($Pe = 10.000$), com baixos números de Froude ($Fr = 0,15$e$Fr = 0,1$), resultando em números de Reynolds de flutuabilidade ($Re_b = Fr^2 Re$) de 225 e 100, respectivamente.

Para isolar o efeito das condições iniciais na seleção do caminho, dois conjuntos de DNS são realizados:

1. DNS1 (Caminho P1): Inicializado com ruído branco em escala de grade, permitindo que todos os modos instáveis da instabilidade de cisalhamento horizontal primária cresçam.
2. DNS2 (Caminho P2): Inicializado com ruído de escala de grade mais uma "semente" específica para o modo eigen do modo 2D ($k_z=0$) de crescimento mais rápido da instabilidade primária.

Adicionalmente, os autores realizam uma análise de estabilidade linear de um modelo analítico "congelado" do fluxo colunar que emerge em DNS2. Esse modelo aproxima o fluxo 2D evoluído não linearmente como um arranjo de vórtices de Taylor-Green sobreposto a fluxos sinusoidais, permitindo o cálculo de eigenmodos e taxas de crescimento de perturbações 3D secundárias.

Principais Contribuições e Resultados

1. Confirmação de Dois Caminhos em uma Única Configuração:
   O estudo demonstra que ambos os caminhos para a turbulência podem ocorrer dentro da exata mesma configuração de modelo, distinguidos unicamente pelo espectro de perturbação inicial:

• Caminho P1 (Crescimento 3D Direto): Quando inicializado com ruído branco, os eigenmodos de $k_z \neq 0$ da instabilidade de cisalhamento horizontal primária crescem diretamente. Esses modos geram cisalhamento vertical em uma ampla gama de escalas, que subsequentemente torna-se instável a instabilidades KH secundárias de pequena escala quando o $Re_b$ é suficientemente grande. Este caminho alinha-se com observações em fluxos de Kolmogorov forçados por corpo.
• Caminho P2 (Evolução Colunar): Quando o modo $k_z=0$ é semeado, ele domina o fluxo, conduzindo o fundo para um fluxo colunar (vortical) 2D de longa duração e dependente do tempo. Este estado 2D é subsequentemente instável a modos 3D secundários. Os autores identificam esses modos secundários como instabilidades hiperbólicas decorrentes dos pontos de estagnação hiperbólicos no arranjo de vórtices. Essas instabilidades possuem altos números de onda verticais ($k_z \sim Fr^{-1}$) e impulsionam a emergência de cisalhamento vertical, que eventualmente desencadeia instabilidades KH terciárias.

2. Caracterização da Instabilidade Secundária no Caminho P2:
   Através da análise de estabilidade linear do fluxo colunar congelado, os autores mostram que a instabilidade secundária é uma instabilidade hiperbólica de alto $k_z$. A taxa de crescimento desta instabilidade escala com a taxa de deformação local dos pontos hiperbólicos ($\Delta_h$) e é aproximadamente $\lambda \approx 0,46 \Delta_h$, sendo amplamente independente da força da estratificação para grandes amplitudes. A instabilidade manifesta-se como "ondulações" no campo de velocidade vertical centradas em torno de pontos hiperbólicos, distintas da tradicional instabilidade "zigzag" observada em pares de vórtices isolados.

3. Impacto na Eficiência de Mistura:
   O estudo quantifica a mistura diapicnal e a eficiência de mistura ($\eta$) para ambos os caminhos.

• Dissipação de Flutuabilidade ($\chi$): A taxa de pico de dissipação da variância da flutuabilidade é comparável em ambos os caminhos, pois é controlada por movimentos de pequena escala ($k_z \gtrsim Fr^{-1}$) em ambos os casos.
• Dissipação Viscosa ($\epsilon$): O Caminho P2 exibe uma dissipação viscosa significativamente maior do que o Caminho P1. Isso ocorre porque a instabilidade hiperbólica em P2 excita um espectro estreito de alto $k_z$, levando a uma dissipação intensa em camadas de cisalhamento vertical finas. Em contraste, o P1 excita um espectro amplo de escalas verticais, distribuindo a dissipação de forma mais ampla.
• Eficiência de Mistura: Consequentemente, a eficiência de mistura de pico é significativamente maior no Caminho P1 ($\eta_{max} \approx 0,25$) do que no Caminho P2 ($\eta_{max} \approx 0,165$). A diferença surge porque o P2 dissipa mais energia cinética em relação à quantidade de flutuabilidade misturada.

Significância e Alegações
O artigo alega que a emergência de cisalhamento vertical impulsionando instabilidades KH de pequena escala é um subproduto inevitável das instabilidades de cisalhamento horizontal em fluxos fortemente estratificados para valores suficientemente grandes de $Re_b$, desde que o fluxo possua um ponto de inflexão. A seleção entre os dois caminhos é determinada pelas condições iniciais: especificamente, a razão da amplitude inicial do modo $k_z=0$ em relação aos modos $k_z \neq 0$.

Os autores sugerem que a prevalência do Caminho P2 na literatura anterior (especificamente em estudos de camadas de cisalhamento de tangente hiperbólica) pode ser um artefato de práticas de inicialização onde uma semente $k_z=0$ é frequentemente adicionada para acelerar os cálculos, thereby suprimindo o crescimento direto de modos 3D (Caminho P1). Por outro vez, a ausência de P1 em alguns contextos pode ser devida a limites de estratificação mais fraca, onde o intervalo de modos $k_z \neq 0$ instáveis é limitado.

Em última análise, o trabalho estabelece que as instabilidades de cisalhamento horizontal são atores fundamentais na geração de turbulência estratificada e mistura diapicnal, mesmo na ausência de cisalhamento vertical pré-existente. Ele destaca que a rota específica para a turbulência (crescimento 3D direto vs. transição 2D-para-3D) altera fundamentalmente o espectro de dissipação de energia e, consequentemente, a eficiência da mistura, um fator que deve ser considerado ao parametrizar a turbulência em modelos climáticos e geofísicos. Os autores observam que, embora seus resultados sejam robustos para fluxos idealizados, a interação com a rotação e campos de ondas ambientes em cenários geofísicos realistas permanece uma questão aberta.