Navier-Stokes Incompressível vs. Compressível

"Incompressível vs compressível" resume-se à densidade. Ela permanece constante ou muda?

Experimente. Encha uma seringa com água e empurre o êmbolo. A água se move, mas em condições cotidianas ela não comprime visivelmente. A água resiste à compressão tão fortemente que tratá-la como incompressível é uma excelente aproximação. Incompressível. Agora encha essa seringa com ar e sele a extremidade. Empurre o êmbolo e você sentirá o ar ceder, a mesma massa de ar agora compactada em menos volume à medida que se comprime sob seu polegar. Isso é escoamento compressível.

No escoamento incompressível, a densidade $\rho$ é constante em todo o fluido, e cada pequena parcela mantém seu volume à medida que se move pelo espaço. O escoamento compressível é diferente. A densidade torna-se uma variável, livre para mudar de lugar para lugar e de momento a momento. Ar ao redor de um motor a jato, gás em uma explosão, a atmosfera em grandes escalas: todos compressíveis, todos impulsionados por variações de densidade.

Por que se importar? Porque essa distinção reformula como as equações de Navier-Stokes se parecem, o que elas preveem e quão difíceis são de analisar e resolver.

A hipótese de incompressibilidade define o campo de densidade $\rho$ como uma constante positiva em todo o domínio do escoamento. Fisicamente, isso significa que cada elemento de volume material preserva seu volume sob o mapa de escoamento. A configuração compressível promove $\rho(x,t)$ a uma incógnita completa, governada por sua própria equação de evolução.

Estes não são duas notações para o mesmo sistema. Eles diferem no número de incógnitas, na estrutura das equações de restrição e no caráter da pressão. O sistema incompressível tem $d+1$ incógnitas escalares ($u$ e $p$ em $\mathbb{R}^d$); o sistema compressível tipicamente tem $d+2$ campos primários (velocidade, densidade e energia interna ou temperatura), com a pressão determinada através de uma equação de estado.

A distinção também altera o tipo de EDP da pressão: elíptica no modelo incompressível, onde a pressão é determinada por uma resolução espacial instantânea, versus uma estrutura mista hiperbólico-parabólica no caso compressível viscoso. Modos acústicos na parte invíscida ou hiperbólica se propagam à velocidade finita do som, enquanto viscosidade e condução de calor introduzem difusão parabólica matematicamente.

As equações de Navier-Stokes incompressíveis descrevem fluidos cuja densidade é constante. Elas são a versão que aparece no Problema do Milênio Clay e a versão em que este site se concentra.

O sistema tem duas partes. A equação de momento:

$$\partial_t u + (u \cdot \nabla)u = -\nabla p + \nu \Delta u + f$$

e a restrição de incompressibilidade:

$$\nabla \cdot u = 0$$

A restrição $\nabla \cdot u = 0$ diz que o campo de velocidade é livre de divergência: o fluido não se acumula nem se dilui em lugar algum. O que quer que flua para dentro de uma pequena região deve fluir para fora na mesma taxa. Esta única condição substitui toda a equação de densidade. A densidade não muda, então você não precisa de uma equação para rastreá-la.

A pressão desempenha um papel especial aqui. Ela não é determinada por uma lei termodinâmica (como a lei dos gases ideais). Em vez disso, ela se ajusta instantaneamente em todos os lugares para manter o escoamento livre de divergência. Matematicamente, $p$ resolve uma equação de Poisson derivada da restrição. Mudanças de pressão se propagam infinitamente rápido. Não há "velocidade do som" no escoamento incompressível.

O sistema de Navier-Stokes incompressível tem dois campos desconhecidos: velocidade $u$ e pressão $p$. Essa simplicidade é enganosa. O termo não linear $(u \cdot \nabla)u$ ainda torna o sistema extremamente difícil em três dimensões.

O sistema de Navier-Stokes incompressível em $\mathbb{R}^3$ com viscosidade cinemática $\nu > 0$:

$$\partial_t u + (u \cdot \nabla)u = -\nabla p + \nu \Delta u + f, \qquad x \in \mathbb{R}^3,\; t > 0,$$

$$\nabla \cdot u = 0, \qquad u(x,0) = u_0(x).$$

A condição livre de divergência $\nabla \cdot u = 0$ é uma restrição pontual, não uma equação de evolução. Ela codifica a preservação local de volume: o mapa de escoamento $\Phi_t$ satisfaz $\det(D\Phi_t) = 1$ para todo $t$.

Aplicando o operador divergência à equação de momento e usando a incompressibilidade obtém-se a equação de Poisson para a pressão:

$$-\Delta p = \partial_i \partial_j (u_i u_j) - \nabla \cdot f.$$

Esta é uma equação elíptica para $p$ em cada tempo fixo. A pressão não é uma variável termodinâmica independente; na interpretação padrão de EDP, ela atua como um multiplicador de Lagrange impondo a restrição livre de divergência, determinada global e instantaneamente pelo campo de velocidade. A informação se propaga a velocidade infinita através da pressão, uma diferença estrutural do sistema compressível que não pode ser disfarçada.

As incógnitas são $u : \mathbb{R}^3 \times [0,T) \to \mathbb{R}^3$ e $p : \mathbb{R}^3 \times [0,T) \to \mathbb{R}$. Para a formulação Clay (Fefferman, 2000), $u_0 \in C^\infty(\mathbb{R}^3)$ é livre de divergência e a questão é se $u$ permanece em $C^\infty(\mathbb{R}^3 \times [0,\infty))$ com energia limitada.

As equações de Navier-Stokes compressíveis governam escoamentos onde a densidade varia. Sistema maior. Mais incógnitas. Mais equações.

Você ainda tem uma equação de momento, mas agora a densidade $\rho$ aparece explicitamente:

$$\partial_t (\rho u) + \nabla \cdot (\rho u \otimes u) = -\nabla p + \nabla \cdot \tau + \rho f$$

A restrição $\nabla \cdot u = 0$ desaparece. Em seu lugar, você obtém uma equação de continuidade que rastreia como a densidade evolui:

$$\partial_t \rho + \nabla \cdot (\rho u) = 0$$

Isso diz que a massa é conservada: a densidade muda porque o escoamento comprime ou expande parcelas de fluido.

O sistema também precisa de uma equação de energia e uma equação de estado, uma relação termodinâmica como $p = \rho R T$ (a lei dos gases ideais) que liga a pressão à densidade e temperatura. A pressão não é mais um impositor passivo de uma restrição. Sua parte acústica viaja à velocidade do som, enquanto o sistema viscoso completo também contém efeitos difusivos.

O sistema compressível é essencial para aerodinâmica em altas velocidades, dinâmica de gases astrofísicos, combustão e qualquer escoamento onde mudanças de densidade importam. Mas é um objeto matemático genuinamente diferente das equações incompressíveis. Mais incógnitas, mais equações, estrutura de EDP completamente diferente.

O sistema de Navier-Stokes compressível acopla a velocidade $u(x,t)$, densidade $\rho(x,t)$, pressão $p(x,t)$ e energia interna específica $e(x,t)$ (ou temperatura $\theta$). Na forma de conservação:

Continuidade: $$\partial_t \rho + \nabla \cdot (\rho u) = 0$$

Momento: $$\partial_t (\rho u) + \nabla \cdot (\rho u \otimes u) + \nabla p = \nabla \cdot \tau + \rho f$$

Energia: $$\partial_t (\rho E) + \nabla \cdot ((\rho E + p)u) = \nabla \cdot (\tau \cdot u) + \nabla \cdot (\kappa \nabla \theta) + \rho f \cdot u$$

onde $E = e + \tfrac{1}{2}|u|^2$ é a energia específica total, $\tau$ é o tensor de tensão viscosa (para um fluido newtoniano, $\tau = \mu(\nabla u + \nabla u^T) + \lambda (\nabla \cdot u)I$ com viscosidade volumétrica $\lambda$), e $\kappa$ é a condutividade térmica.

O fechamento requer uma equação de estado, por exemplo, $p = (\gamma - 1)\rho e$ para um gás ideal com índice adiabático $\gamma$.

Na parte invíscida ou hiperbólica do sistema compressível, perturbações acústicas se propagam a velocidade finita, com velocidade do som característica $c = \sqrt{\partial p / \partial \rho |_s}$. O sistema viscoso completo de Navier-Stokes-Fourier é misto hiperbólico-parabólico, então difusão viscosa e térmica introduzem propagação infinita matematicamente. As equações compressíveis suportam ondas de choque, rarefações e descontinuidades de contato que não têm análogo no escoamento incompressível.

Quando a compressibilidade importa? Um número decide: o número de Mach.

$$\text{Ma} = \frac{|u|}{c}$$

$|u|$ é a velocidade do escoamento. $c$ é a velocidade do som. Sua razão diz quão rápido o escoamento se move em comparação com a velocidade na qual perturbações de pressão podem se propagar através do meio, e essa comparação determina se você pode ignorar com segurança mudanças de densidade ou se elas dominarão a física.

Quando $\text{Ma} < 0.3$, a densidade muda menos de cerca de 5%. As equações incompressíveis funcionam. Ar em uma sala, água em um cano, vento ao redor de um edifício: todos escoamentos de baixo Mach onde perturbações de pressão viajam muito mais rápido que o próprio escoamento, de modo que a densidade mal se move.

Acima de $\text{Ma} \approx 0.3$, a compressibilidade começa a morder, e em torno de $\text{Ma} \approx 1$ você atinge o regime transônico onde bolsões supersônicos locais aparecem e ondas de choque se formam. Jatos de combate. Bocais de foguetes. Naves espaciais em reentrada.

Não é um interruptor binário. A maioria dos escoamentos de fluidos cotidianos, e o Problema do Milênio Clay, situam-se firmemente no regime de baixo Mach onde as equações incompressíveis se aplicam.

O número de Mach $\text{Ma} = |u|/c$ parametriza a importância da compressibilidade, onde $c = \sqrt{\partial p / \partial \rho |_s}$ é a velocidade do som isentrópica. Formalmente, as equações incompressíveis surgem como o limite de baixo Mach do sistema compressível.

A expansão assintótica em potências de $\text{Ma}^2$ (ver Klainerman & Majda, 1981, 1982; Schochet, 1986) mostra que à medida que $\text{Ma} \to 0$ com dados iniciais adequados, as soluções compressíveis convergem para a solução incompressível. Em adimensionalizações padrão de baixo Mach com dados bem preparados, frequentemente se escreve a pressão como um fundo termodinâmico quase espacialmente uniforme mais uma correção dinâmica menor que impõe a restrição de incompressibilidade no limite.

Este é um limite singular: a velocidade do som $c \to \infty$ e o caráter hiperbólico do sistema compressível degenera para a equação de pressão elíptica do escoamento incompressível. Os modos acústicos tornam-se infinitamente rápidos e se desacoplam da dinâmica vorticosa.

O regime $\text{Ma} < 0.3$ é uma heurística de engenharia. Reflete a observação empírica de que a variação relativa de densidade $\delta\rho / \rho \sim \text{Ma}^2 / 2$ permanece abaixo de $\sim$5% nesta faixa. A justificativa matemática é o teorema de convergência para o limite de baixo Mach, que requer dados iniciais bem preparados e condições de contorno compatíveis.

O Problema do Milênio Clay faz uma pergunta precisa: dada uma velocidade inicial suave e livre de divergência em $\mathbb{R}^3$, o sistema de Navier-Stokes incompressível sempre produz uma solução suave que existe para todo tempo?

Por que incompressível especificamente? Três razões.

Primeiro, já é difícil o suficiente. As equações incompressíveis 3D resistiram à prova de regularidade global desde o trabalho fundamental de Leray em 1934. Adicionar densidade variável, termodinâmica e ondas de choque tornaria o problema vastamente mais difícil, não mais tratável.

Segundo, a dificuldade é mecânica de fluidos pura. O sistema incompressível isola o desafio matemático central, a competição entre advecção não linear $(u \cdot \nabla)u$ e dissipação viscosa $\nu \Delta u$, sem complicações termodinâmicas ou acústicas. É a arena mais limpa para fazer a questão de regularidade.

Terceiro, a física é limpa. As equações incompressíveis modelam os escoamentos cotidianos mais comuns. Se elas podem produzir singularidades a partir de dados suaves é uma questão fundamental sobre a consistência matemática da mecânica de fluidos clássica.

O sistema compressível tem seus próprios problemas profundos em aberto (existência de soluções globais com dados grandes, formação e interação de choques), mas esses são problemas diferentes com estruturas diferentes. O prêmio Clay visa o caso incompressível porque essa é a questão específica de regularidade que Fefferman formulou para Navier-Stokes 3D.

A formulação oficial Clay (Fefferman, 2000) especifica o sistema incompressível em $\mathbb{R}^3$:

$$\partial_t u + (u \cdot \nabla)u = -\nabla p + \nu \Delta u, \qquad \nabla \cdot u = 0, \qquad u|_{t=0} = u_0,$$

com $u_0 \in C^\infty(\mathbb{R}^3)$ livre de divergência, e a questão é se $u \in C^\infty(\mathbb{R}^3 \times [0,\infty))$ com $\int_{\mathbb{R}^3} |u(x,t)|^2\,dx$ limitada para todo $t \geq 0$.

A escolha do sistema incompressível é matematicamente motivada. O principal problema em aberto, a lacuna entre soluções fracas de Leray-Hopf (que existem globalmente mas podem não ser únicas ou suaves) e soluções clássicas suaves (que existem localmente mas podem explodir), é específico das equações incompressíveis 3D. Em 2D, a regularidade global para soluções suaves de Navier-Stokes incompressíveis é conhecida; o caso 3D permanece em aberto.

O sistema compressível introduz dificuldades qualitativamente diferentes: formação de choques (que ocorre mesmo para equações de Euler com dados suaves), estados de vácuo ($\rho \to 0$) e o acoplamento entre vorticidade e modos acústicos. Estes são problemas importantes em aberto, mas são estruturalmente distintos da questão de regularidade incompressível.

O problema incompressível isola a competição entre a não linearidade supercrítica em energia e a dissipação viscosa. A estimativa de energia natural dá $u \in L^\infty_t L^2_x \cap L^2_t \dot{H}^1_x$, que fica aquém do espaço crítico de escala $L^\infty_t \dot{H}^{1/2}_x$ por meia derivada em 3D. Fechar essa lacuna, ou provar que não pode ser fechada, é o coração do Problema do Milênio.

Comece aqui. Quer cada termo no sistema incompressível desmembrado, com o significado físico e o papel matemático de cada peça explicados do zero? O Que São as Equações de Navier-Stokes?

De onde vem este sistema? Derivação das Equações de Navier-Stokes.

Elimine a viscosidade e você obtém as equações de Euler, que são um século mais antigas, parecem mais simples no papel e, de certa forma, são ainda mais difíceis de entender matematicamente porque você perde o efeito de suavização do termo de difusão. Euler vs. Navier-Stokes.

O prêmio. O Problema de Existência e Suavidade de Navier-Stokes.

Caminhos sugeridos a partir daqui:

• O Que São as Equações de Navier-Stokes?, o sistema incompressível completo com análise termo a termo
• Derivação das Equações de Navier-Stokes, dos axiomas da mecânica do contínuo ao sistema de EDP
• Euler vs. Navier-Stokes, o limite invíscido $\nu \to 0$ e o papel da viscosidade na regularidade
• O Problema de Existência e Suavidade de Navier-Stokes, a formulação precisa Clay, teoria de Leray-Hopf e a lacuna entre soluções fracas e fortes