Instabilidade de Kelvin-Helmholtz — Escoamento Cisalhado Bifásico

Duas camadas de fluido que deslizam uma contra a outra se desestabilizam na interface formando os característicos rolos KH. Ajuste a razão de densidades, o cisalhamento, a tensão superficial e ative a estabilização magnética (MHD).

Evolução da interface e campo de velocidade

O fluido inferior pesado (azul) desliza contra o fluido superior leve (laranja). Quando o cisalhamento supera a gravidade, a tensão e a rigidez magnética, a interface se enrola em rolos KH. As setas mostram o campo de velocidade perturbado.

Taxa de crescimento γ(k)

Parte imaginária de ω. Pico no número de onda mais instável k*; a interface amplifica preferencialmente esse comprimento de onda.

Dispersão ωᵣ(k)

Parte real da frequência. Onde γ=0 é uma onda estável gravidade-capilaridade-Alfvén; onde γ>0 a onda é sobrescrita por crescimento exponencial.

O mecanismo de Kelvin-Helmholtz

Quando duas camadas de fluido deslizam uma contra a outra, o salto de velocidade na interface é uma lâmina de vorticidade. Qualquer pequena corrugação da interface reduz a pressão na crista (Bernoulli) e eleva o vale, então a corrugação cresce — a clássica instabilidade KH. Para densidades iguais sem gravidade nem tensão, todos os comprimentos de onda são instáveis e a taxa de crescimento cresce com o número de onda k, então as ondas mais curtas vencem até que a tensão superficial as corte. A teoria linear dá a relação de dispersão (ρ₁+ρ₂)ω² − 2(ρ₁U₁+ρ₂U₂)kω + (ρ₁U₁²+ρ₂U₂² − B²/μ₀)k² − (ρ₁−ρ₂)gk − σk³ = 0, cujas raízes complexas ω = ωᵣ − iγ dão uma taxa de crescimento positiva γ onde o cisalhamento domina.

O que estabiliza a interface

Três efeitos combatem a instabilidade. (1) Estratificação de densidade estável — quando o fluido mais pesado está abaixo (ρ₁>ρ₂), a gravidade (ρ₁−ρ₂)gk atua como força restauradora em ondas longas; é por isso que a superfície do mar só se quebra em ondas acima de uma velocidade crítica do vento. (2) A tensão superficial σk³ mata os comprimentos de onda mais curtos, definindo um corte em alto k; é por isso que pequenas ondulações são suaves. (3) Um campo magnético horizontal adiciona tensão magnética B²k²/μ₀ (efeito Alfvén), que pode suprimir KH inteiramente — é por isso que a magnetopausa do vento solar pode permanecer nítida. Aumentar ΔU ou a razão de densidades empurra de volta para a instabilidade.

Onde aparecem os rolos KH

A instabilidade KH é onipresente. Vento sobre a água semeia ondas oceânicas e spray. A atmosfera em faixas de Júpiter é uma galeria de vórtices enrolados por KH, incluindo a fronteira de cisalhamento da Grande Mancha Vermelha. A magnetopausa terrestre — fronteira entre o vento solar e a magnetosfera — é uma camada de cisalhamento magnetizada cujas ondas KH deixam o plasma do vento solar entrar. KH também impulsiona a mistura em transbordamentos oceânicos, camadas de nuvens por cisalhamento (a clássica "nuvem em rolo"), escapes de motores turbofan e a interface combustível-ar na combustão scramjet. Em astrofísica molda bordas de remanescentes de supernova e turbulência de discos de acreção.

Como usar

Comece com Cisalhamento clássico: densidades iguais, cisalhamento puro — todos os comprimentos de onda são instáveis e os rolos se enrolam rápido. Agora adicione Gravidade (ρ₂<ρ₁) para imitar uma estratificação estável e observe as ondas longas serem suprimidas. Aumente a Tensão superficial para recortar as ondas mais curtas e veja a curva de crescimento ganhar um corte em alto k. Ative o Modo MHD e aumente a tensão magnética: a taxa de crescimento cai em todo k, e um campo forte o suficiente estabiliza a interface totalmente (caso magnetopausa). O painel de Taxa de crescimento destaca o modo mais instável k*; o de Dispersão mostra a frequência de onda real que sobrevive quando γ=0.