Por Que o Problema de Navier-Stokes É Difícil

Muitas das equações que as pessoas encontram primeiro em física são lineares: dobre a entrada e a resposta dobra. Navier-Stokes não é assim.

Navier-Stokes? Não linear. A velocidade do fluido afeta sua própria taxa de mudança, o que significa que o fluido empurra a si mesmo. Imagine tentar prever para onde uma multidão vai quando o movimento de cada pessoa depende do que todos ao redor estão fazendo, e o que essas pessoas estão fazendo depende de todos ao redor delas, em espiral para fora infinitamente. Essa é a situação que você está encarando.

O culpado é o termo de autointeração $(u \cdot \nabla)u$. Ele cria ciclos de retroalimentação onde pequenas perturbações se amplificam em grandes, e é por isso que a turbulência de fluidos é tão extremamente complexa (veja subproblemas para mais).

A não linearidade convectiva $(u \cdot \nabla)u$ é o obstáculo fundamental. Na formulação de vorticidade $\omega = \nabla \times u$, a equação torna-se

$$\partial_t \omega + (u \cdot \nabla)\omega = (\omega \cdot \nabla)u + \nu \Delta \omega$$

O termo de estiramento de vórtice $(\omega \cdot \nabla)u$ não tem sinal. Ele pode amplificar a vorticidade sem limite. Em 2D este termo desaparece (já que $\omega$ é um escalar perpendicular ao fluxo), razão pela qual a regularidade global em 2D é conhecida (Ladyzhenskaya, 1969). Em 3D, o estiramento de vórtice é o principal mecanismo candidato para explosão em tempo finito.

Aqui está o cerne: a não linearidade é quadrática em $u$. A estimativa de energia $H^1$ fornece $\|\nabla u\|_{L^2}$, mas controlar $(u \cdot \nabla)u$ em $L^2$ tipicamente requer informação mais forte, como $u \in L^\infty$ junto com $\nabla u \in L^2$, ou controle equivalente em escala crítica; a classe de energia sozinha não fornece isso.

As equações de Navier-Stokes têm uma simetria de escala. Amplie uma solução, torne tudo menor e mais rápido nas quantidades certas, e você obtém outra solução perfeitamente válida. Essa simetria é matematicamente natural, mas analiticamente perigosa.

Por quê? A única quantidade que podemos controlar de forma confiável é a energia total do fluido, e ela está em uma escala completamente errada, nos dizendo sobre o panorama geral mas não dizendo absolutamente nada sobre o que está acontecendo nas escalas microscópicas onde uma explosão realmente se formaria.

Pense em monitorar o uso total de eletricidade de uma cidade para detectar uma faísca. Útil? Claro. Detalhado o suficiente? Nem de longe. Essa lacuna é todo o problema.

Sob o escalonamento natural $u_\lambda(x,t) = \lambda u(\lambda x, \lambda^2 t)$, o espaço de Sobolev crítico é $\dot{H}^{1/2}(\mathbb{R}^3)$ (equivalentemente $L^3$). Uma norma é crítica em escala se permanece invariante sob essa transformação. Se uma norma controlada fica menor ao ampliar escalas pequenas, então a equação é supercrítica em relação a esse controle: o limite conhecido fica abaixo da escala necessária para enxergar concentração.

A desigualdade de energia fornece controle de $u \in L^\infty_t L^2_x \cap L^2_t \dot{H}^1_x$. No nível da energia, o problema 3D é supercrítico:

$$\|u_\lambda\|_{L^2} = \lambda^{-1/2} \|u\|_{L^2}, \quad \|u_\lambda\|_{L^2_t \dot{H}^1_x} = \lambda^{-1/2} \|u\|_{L^2_t \dot{H}^1_x}$$

Quando $\lambda \to \infty$ (ampliando escalas pequenas), esses limites enfraquecem em relação às quantidades invariantes $\|u\|_{L^3}$ e $\|u\|_{\dot{H}^{1/2}}$. Os termos não linear e viscoso escalam juntos; o problema é que a estimativa a priori conhecida fica abaixo da escala crítica. Passar da classe de energia para uma norma crítica é a dificuldade central.

Observe um rio. O movimento do fluido torna-se caótico. Turbulento. Grandes redemoinhos se fragmentam em menores, que se fragmentam em ainda menores, em cascata até escalas microscópicas onde a viscosidade finalmente suaviza as coisas.

Kolmogorov descreveu esta cascata de energia em 1941, e as equações de Navier-Stokes a capturam lindamente. Mas aqui está o que mantém as pessoas acordadas à noite: e se a energia se concentrar em regiões cada vez menores mais rápido do que a viscosidade pode dissipá-la? Isso é uma explosão.

Isso pode realmente acontecer? Ou a viscosidade sempre vence? Essa é a questão em aberto, ponto final. Para a ponte física do número de Reynolds para esta imagem de pequenas escalas, veja Número de Reynolds, Turbulência e Por Que Pequenas Escalas Importam.

A teoria K41 de Kolmogorov prevê um espectro de energia $E(k) \sim \varepsilon^{2/3} k^{-5/3}$ na faixa inercial $k_f \ll k \ll k_\eta$, onde $k_\eta \sim (\varepsilon/\nu^3)^{1/4}$ é o número de onda de dissipação de Kolmogorov. Fluxo constante de energia através das escalas. Essa é a imagem limpa.

A questão de regularidade pergunta algo mais sombrio: essa cascata pode degenerar? A escala de dissipação $k_\eta^{-1}$ pode encolher para zero em tempo finito, com $\|\nabla u\|_{L^2} \to \infty$ enquanto a energia total permanece finita?

A conjectura da anomalia de dissipação de Onsager (1949) diz que sim, em certo sentido: no limite de viscosidade desaparecendo $\nu \to 0$, a dissipação de energia persiste, e soluções fracas de Euler podem dissipar energia. Confirmado para expoentes de Hölder abaixo de $1/3$ (Isett, 2018; Buckmaster et al., 2018). Mas o que isso significa para a regularidade de Navier-Stokes permanece genuinamente incerto. Para intuição em nível de regime, veja Número de Reynolds, Turbulência e Por Que Pequenas Escalas Importam.

A pressão nas equações de Navier-Stokes é estranha. Ela não é uma variável independente; a velocidade a determina completamente através de uma única restrição: o fluido é incompressível, então não pode ser comprimido.

Isso torna a pressão não local. No modelo incompressível, mudar a velocidade em uma região afeta o campo de pressão globalmente, porque a pressão é determinada por todo o campo de velocidade naquele momento.

Para análise, isso é devastador. Você não pode estudar o que acontece em um único ponto sem considerar todo o fluido de uma vez. Raciocínio local? Esqueça.

A restrição de incompressibilidade $\nabla \cdot u = 0$ determina a pressão através da equação de Poisson

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

então $p = (-\Delta)^{-1} \partial_i \partial_j (u_i u_j)$, envolvendo transformadas de Riesz (operadores integrais singulares). A pressão é uma função não local da velocidade, e esta não localidade é o obstáculo central para estimativas pontuais ou locais no espaço.

Argumentos padrão de princípio do máximo falham aqui: embora o termo viscoso $\nu \Delta u$ seja dissipativo, o gradiente de pressão $-\nabla p$ pode concentrar energia de regiões distantes. A teoria de Caffarelli-Kohn-Nirenberg lida com isso via desigualdades de energia locais em cilindros parabólicos, mas extrair regularidade pontual destas permanece o passo difícil.

Para as equações de Navier-Stokes incompressíveis em 2D nas configurações padrão, soluções suaves globais são conhecidas; isso foi estabelecido em trabalho clássico incluindo o de Ladyzhenskaya (1969). Veja por que 2D é mais fácil.

Três dimensões? Tudo desmorona, e a razão se resume a um mecanismo: estiramento de vórtice. Em 2D, vórtices podem girar e se fundir, mas não podem se esticar. Em 3D, o fluido pode agarrar tubos de vórtice e puxá-los cada vez mais finos, potencialmente concentrando cada pedaço de energia em um filamento infinitamente fino.

Esse estiramento pode fugir ao infinito em tempo finito, ou a viscosidade sempre intervém? Essa é a questão de um milhão de dólares. Literalmente.

A dicotomia é nítida.

2D: Vorticidade $\omega$ é um escalar satisfazendo $\partial_t \omega + u \cdot \nabla \omega = \nu \Delta \omega$. O princípio do máximo fornece $\|\omega(t)\|_{L^\infty} \leq \|\omega_0\|_{L^\infty}$, BKM implica regularidade global, e o termo de estiramento de vórtice $(\omega \cdot \nabla)u$ é identicamente zero.

3D: Vorticidade $\omega \in \mathbb{R}^3$ satisfaz $\partial_t \omega + (u \cdot \nabla)\omega = (\omega \cdot \nabla)u + \nu \Delta \omega$, onde o termo de estiramento $(\omega \cdot \nabla)u$ pode amplificar $|\omega|$ superlinearmente (formalmente $\sim |\omega|^2$ via Biot-Savart). Nenhum princípio do máximo disponível.

A enstrofia $\|\omega\|_{L^2}^2$ satisfaz

$$\frac{d}{dt}\|\omega\|_{L^2}^2 \leq C\|\omega\|_{L^2}^2 \|\nabla u\|_{L^\infty} - 2\nu \|\nabla \omega\|_{L^2}^2$$

Controlar $\|\nabla u\|_{L^\infty}$ requer $\omega \in L^\infty$, que requer controlar $\|\nabla u\|_{L^\infty}$. Circular. Nenhuma técnica existente quebrou este ciclo.

Este artigo é parte de O Problema.

Esses obstáculos levaram matemáticos a decompor o problema em subproblemas e desenvolver abordagens especializadas para cada um.

Para contexto sobre por que o termo viscoso ajuda mas não é suficiente, veja Euler vs. Navier-Stokes.

Este artigo é parte de O Problema.

O problema de regularidade se divide em peças tratáveis; veja Subproblemas. As ferramentas analíticas construídas para atacar esses obstáculos são cobertas em Abordagens. Por que o termo viscoso ajuda mas ainda não fecha a lacuna? Euler vs. Navier-Stokes.