Soluções Fracas, Fortes e Suaves para as Equações de Navier-Stokes

Aqui está a situação. O Prêmio do Milênio oferece um milhão de dólares para resolver se o Navier-Stokes 3D sempre tem soluções suaves que duram para sempre, ou se as coisas podem explodir. Suave significa que o campo de velocidade é perfeitamente bem comportado: sem saltos repentinos, sem velocidades infinitas, sem pontos onde a matemática falha. Mas o melhor resultado de existência que alguém já provou, em quase um século de tentativas, só garante algo mais fraco. Estas são chamadas soluções fracas.

Então, o que é uma solução fraca? Não é uma aproximação. Não é "quase certo". É uma solução exata para as equações, mas que segue regras relaxadas. Uma solução normal ("clássica") requer que a velocidade seja suave o suficiente para que você possa calcular sua taxa de mudança em cada ponto. Uma solução fraca pula esse requisito. Em vez de verificar as equações ponto a ponto, você as verifica "em média" através de regiões do espaço.

Aqui está uma analogia. Uma solução clássica é um estudante que resolve cada problema do exame mostrando todo o seu trabalho, passo a passo. Uma solução fraca é um estudante que não pode mostrar os passos intermediários, mas cujas respostas finais são comprovadamente corretas para cada possível pergunta que você poderia fazer. Você não pode vê-los trabalhar, mas as respostas sempre conferem.

Por que você aceitaria isso? Porque às vezes as equações são selvagens demais para soluções clássicas. O fluido pode desenvolver regiões onde a velocidade muda tão bruscamente que você simplesmente não pode calcular uma taxa de mudança ali. A matemática falha. Soluções fracas permitem que você continue onde as soluções clássicas desistem. Elas são a rede de segurança que mantém as equações vivas quando as coisas ficam difíceis.

O problema: soluções fracas podem não ser únicas. Você pode obter múltiplas soluções fracas começando do exato mesmo fluxo, e ninguém pode dizer qual é "a resposta real". Isso é um problema, porque a física diz que o fluido deve fazer uma coisa específica, não várias. E soluções fracas podem não ser suaves. Suavidade é o que o Prêmio do Milênio exige, e é o que ninguém pode provar.

Uma solução fraca para as equações de Navier-Stokes substitui a EDP pontual por uma formulação distribucional. Considere o sistema incompressível em $\mathbb{R}^3 \times (0,T)$:

$$\partial_t u + (u \cdot \nabla)u = -\nabla p + \nu \Delta u, \qquad \nabla \cdot u = 0.$$

Um campo vetorial livre de divergência $u \in L^2_{\mathrm{loc}}(\mathbb{R}^3 \times [0,T))$ é uma solução fraca se para toda função teste suave, com suporte compacto e livre de divergência $\varphi$:

$$\int_0^T \!\int_{\mathbb{R}^3} \bigl[-u \cdot \partial_t \varphi - (u \otimes u) : \nabla \varphi - \nu \, u \cdot \Delta \varphi \bigr] \, dx \, dt = \int_{\mathbb{R}^3} u_0 \cdot \varphi(x,0) \, dx.$$

A pressão desaparece desta formulação porque $\varphi$ é livre de divergência. Todas as derivadas foram movidas de $u$ para $\varphi$ via integração por partes. O ponto: $u$ não precisa ser classicamente diferenciável. Só precisa ser integrável o suficiente para que essas integrais convirjam.

A formulação fraca não é uma aproximação. Uma solução clássica satisfazendo a EDP pontualmente também satisfaz a formulação fraca para toda função teste (integre por partes na outra direção). A recíproca falha: uma solução fraca não precisa ser suave o suficiente para satisfazer a EDP pontualmente.

Isso importa porque: (1) soluções fracas existem globalmente para dados iniciais $L^2$ (Leray 1934), enquanto soluções clássicas são conhecidas por existir apenas localmente no tempo para 3D; (2) soluções fracas em geral não são conhecidas por serem únicas; (3) a questão de se soluções fracas são sempre suaves é equivalente ao Problema do Milênio de Clay para dados iniciais apropriados.

Em 1934, Jean Leray fez algo que ainda define o campo. Em um único artigo de 73 páginas, ele provou que soluções fracas para as equações de Navier-Stokes 3D existem para todo o tempo, começando de qualquer fluxo inicial razoável. Qualquer. Desde que a velocidade inicial não seja infinitamente energética ou fisicamente sem sentido, Leray garante que você obterá uma solução que dura para sempre. Esta foi a primeira vez que alguém provou um resultado de existência global para as equações 3D, e mais de noventa anos depois, ainda é o teorema de existência incondicional mais forte que temos.

Sua estratégia foi inteligente. As equações reais são difíceis demais para resolver diretamente por causa de como a velocidade do fluido se retroalimenta (essa é a não linearidade). Então Leray borrou as equações ligeiramente, como adicionar um pequeno filtro gaussiano a uma imagem. As equações borradas são mansas o suficiente para resolver. Então ele reduziu o borrão em direção a zero e mostrou que as soluções não se dispersam. Elas se estabelecem em algo que satisfaz as equações originais, não borradas, no sentido fraco.

Mas aqui está o que Leray NÃO provou. Unicidade. Seu método produz pelo menos uma solução fraca, mas pode haver outras começando do mesmo fluxo. Ele não pôde descartar isso. Ele também não provou suavidade. Suas soluções têm energia finita e satisfazem uma desigualdade de energia: o atrito pode drenar energia, mas a energia não pode aparecer espontaneamente do nada. Só isso. Nada mais.

O próprio Leray suspeitava que singularidades poderiam se formar. Ele esboçou como uma poderia parecer: o fluido colapsando em direção a um ponto, cada vez mais rápido, concentrando toda a sua energia em uma região cada vez menor, como um redemoinho encolhendo para um ponto em velocidade infinita. Em 1996, Nečas, Růžička e Šverák provaram que esse colapso auto-similar exato não pode acontecer. O palpite de Leray sobre a forma de uma potencial explosão estava errado. Se a explosão acontece de alguma forma? Ninguém sabe.

Em 1951, Eberhard Hopf estendeu a construção de Leray para fluidos em recipientes limitados (não apenas todo o espaço infinito), e a classe resultante ficou conhecida como soluções fracas de Leray-Hopf: soluções fracas que satisfazem a desigualdade de energia. Esta é a noção padrão. Quando pesquisadores dizem "soluções fracas" sem qualificação adicional, eles quase sempre significam isso.

Mais uma coisa. Mesmo dentro das soluções fracas de Leray-Hopf, há uma subclasse mais exigente chamada soluções fracas adequadas. Estas não apenas satisfazem a desigualdade de energia globalmente (a energia total não cresce). Elas a satisfazem localmente também: a energia não pode secretamente se acumular em um canto do fluido enquanto drena de outro. Caffarelli, Kohn e Nirenberg (CKN) provaram seu famoso resultado de regularidade parcial em 1982 especificamente para esta classe menor. Não confunda as duas: CKN se aplica a soluções fracas adequadas, não a todas as soluções de Leray-Hopf.

O artigo de Leray de 1934 estabeleceu o seguinte: para qualquer $u_0 \in L^2(\mathbb{R}^3)$ livre de divergência, existe pelo menos uma solução fraca $u$ para as equações de Navier-Stokes em $\mathbb{R}^3 \times (0,\infty)$ satisfazendo:

• $u \in L^\infty(0,\infty; L^2(\mathbb{R}^3)) \cap L^2(0,\infty; \dot{H}^1(\mathbb{R}^3))$
• A desigualdade de energia: $\|u(t)\|_{L^2}^2 + 2\nu \int_0^t \|\nabla u(s)\|_{L^2}^2 \, ds \leq \|u_0\|_{L^2}^2$ para q.t. $t > 0$

A construção procede por mollificação. Substitua a não linearidade $(u \cdot \nabla)u$ por $(u_\varepsilon \cdot \nabla)u$ onde $u_\varepsilon = J_\varepsilon * u$ é uma mollificação espacial. O sistema regularizado tem soluções suaves globais (a mollificação mata o pior das interações não lineares). Leray obteve limites de energia uniformes para as soluções regularizadas, então extraiu uma subsequência fracamente convergente. O limite satisfaz a formulação fraca e a desigualdade de energia.

Leray não estabeleceu unicidade. O argumento de compacidade dá existência de pelo menos um ponto de acumulação; diferentes subsequências podem convergir para diferentes limites. A unicidade de soluções fracas de Leray-Hopf em 3D permanece aberta até hoje.

Hopf (1951) adaptou a construção para domínios limitados $\Omega \subset \mathbb{R}^3$ com condições de contorno de Dirichlet, usando aproximação de Galerkin (projeção em subespaços de dimensão finita) em vez de mollificação. A classe resultante, soluções fracas satisfazendo a desigualdade de energia, carrega ambos os nomes: soluções fracas de Leray-Hopf.

O teorema de Caffarelli-Kohn-Nirenberg (1982) diz respeito a uma classe mais restritiva: soluções fracas adequadas, que adicionalmente satisfazem uma desigualdade de energia local da forma

$$\partial_t \left(\frac{|u|^2}{2}\right) + \nabla \cdot \left(\left(\frac{|u|^2}{2} + p\right)u\right) + \nu |\nabla u|^2 \leq \nu \Delta \left(\frac{|u|^2}{2}\right)$$

no sentido de distribuições. CKN provou que para qualquer solução fraca adequada, a medida de Hausdorff parabólica unidimensional do conjunto singular no espaço-tempo é zero. Isso significa que singularidades, se existirem, são extremamente esparsas (têm medida de Hausdorff parabólica unidimensional zero). Mas o teorema não diz nada sobre se singularidades realmente ocorrem, e se aplica apenas a soluções fracas adequadas, não a todas as soluções de Leray-Hopf.

O próprio Leray considerou a possibilidade de explosão auto-similar da forma $u(x,t) = (T-t)^{-1/2} U(x / (T-t)^{1/2})$. Nečas, Růžička e Šverák (1996) provaram que nenhuma explosão auto-similar desse tipo existe para soluções em $L^3(\mathbb{R}^3)$, e Tsai (1998) descartou certos cenários de explosão assintoticamente auto-similar sob hipóteses correspondentes. A forma de potenciais singularidades, se houver, permanece desconhecida.

Soluções fracas existem globalmente. Mas elas podem não ser únicas, e podem não ser suaves. Podemos fazer melhor?

Sim, mas apenas temporariamente. Soluções fortes são a melhoria: elas têm regularidade suficiente para que as equações valham quase em todo lugar, não apenas "em média". Soluções suaves, ou clássicas, são aquelas em que as equações valem ponto a ponto. Para dados iniciais suaves em 3D, soluções fortes existem por um curto período. Quão curto? Isso depende de quão selvagem é o fluxo inicial. Fluxos calmos e suaves obtêm garantias mais longas. Condições iniciais violentas e turbulentas? Microssegundos.

E ninguém pode provar que essas soluções fortes não eventualmente explodem.

Em 1962, James Serrin provou algo como uma regra de promoção. Funciona assim: se uma solução fraca acontece de permanecer bem comportada o suficiente (não muito grande, não concentrando sua energia em regiões cada vez menores), então ela era secretamente suave o tempo todo. Você pode promovê-la. E por um princípio chamado unicidade fraca-forte, ela também é a única solução fraca com essas condições iniciais. Uma solução, suave e única, caso encerrado. Mas se você não pode verificar que a solução permanece mansa? Nada. Você está preso.

Este é um resultado condicional. SE a solução não é muito selvagem, ENTÃO ela é perfeitamente bem comportada. Toda a dificuldade é provar o SE.

Em duas dimensões, as estimativas de energia são fortes o suficiente para que toda solução fraca automaticamente passe no teste de Serrin. Pronto. É por isso que 2D está resolvido. Em 3D, as estimativas ficam apenas um pouco aquém do que você precisaria, e fechar essa lacuna é todo o jogo.

Pesquisadores encontraram outros testes condicionais também, cada um um ângulo diferente de ataque: "Prove esta coisa específica sobre a solução, e eu lhe darei suavidade de graça". Provar qualquer um deles incondicionalmente resolveria o Problema do Milênio. Ninguém conseguiu. Para uma revisão das diferentes estratégias de prova que as pessoas tentaram, há uma página inteira sobre isso.

Uma solução forte para as equações de Navier-Stokes é aquela com regularidade suficiente para que a EDP valha pontualmente (q.t.p.) e o termo não linear $(u \cdot \nabla)u$ seja bem definido como uma função em vez de meramente uma distribuição. Tipicamente isso significa $u \in L^\infty(0,T; H^1(\mathbb{R}^3)) \cap L^2(0,T; H^2(\mathbb{R}^3))$. Uma estrutura relacionada mas distinta é Fujita-Kato (1964), que constrói soluções mild locais em espaços críticos.

A existência local de soluções fortes para dados de Sobolev suficientemente regulares está bem estabelecida. Em espaços críticos, a estrutura de Fujita-Kato (1964) constrói soluções mild locais via um argumento de ponto fixo:

$$u(t) = e^{\nu t \Delta} u_0 - \int_0^t e^{\nu(t-s)\Delta} \mathbb{P} \nabla \cdot (u \otimes u)(s) \, ds,$$

onde $\mathbb{P}$ é a projeção de Leray em campos livres de divergência. Esta equação integral tem uma solução local única pelo princípio da aplicação de contração em espaços de funções adequados. Para dados pequenos em espaços críticos ($L^3$, $\dot{H}^{1/2}$, $BMO^{-1}$), a solução é global.

A questão é se soluções fortes de dados grandes persistem para todo o tempo. O resultado condicional chave é de Serrin (1962): se uma solução fraca de Leray-Hopf satisfaz $u \in L^p(0,T; L^q(\mathbb{R}^3))$ com

$$\frac{2}{p} + \frac{3}{q} \leq 1, \qquad q > 3,$$

então $u$ é suave em $(0,T] \times \mathbb{R}^3$ e é a única solução fraca de Leray-Hopf com os dados iniciais dados. Estas são chamadas condições de Prodi-Serrin (Prodi 1959 estabeleceu um resultado relacionado).

O ponto final $q = 3$ ($p = \infty$) foi resolvido por Escauriaza, Seregin e Šverák (2003): se $u \in L^\infty(0,T; L^3(\mathbb{R}^3))$, então $u$ não explode no tempo $T$. Este é o ponto final da escala de Prodi-Serrin e um dos critérios de continuação mais precisos conhecidos.

Em duas dimensões, o limite de energia $u \in L^\infty_t L^2_x \cap L^2_t H^1_x$ combinado com a desigualdade de Ladyzhenskaya $\|f\|_{L^4}^2 \leq C\|f\|_{L^2}\|\nabla f\|_{L^2}$ (específica para 2D) dá $u \in L^4_t L^4_x$, que satisfaz a condição de Serrin $2/4 + 2/4 = 1$ (na versão 2D com $2/p + 2/q \leq 1$). Em 3D, o limite de energia dá $u \in L^{10/3}_t L^{10/3}_x$ por imersão de Sobolev, que satisfaz $2/(10/3) + 3/(10/3) = 3/2 > 1$. A condição de Serrin falha exatamente pela margem correspondente à lacuna de supercriticalidade. Veja Por Que Navier-Stokes É Difícil para mais sobre esta obstrução estrutural.

Soluções suaves são o padrão ouro. O campo de velocidade é perfeitamente bem comportado em todos os lugares, para todo o tempo. Sem saltos repentinos. Sem velocidades infinitas. Amplie o quanto quiser, e a solução continua sendo agradável.

O Problema do Prêmio do Milênio de Clay, formulado por Charles Fefferman em 2000, faz uma pergunta que cabe em uma ficha. Comece com qualquer campo de velocidade suave e fisicamente razoável preenchendo o espaço tridimensional. A equação de Navier-Stokes sempre produz uma solução suave que dura para sempre, ou você pode encontrar um fluxo inicial onde a solução eventualmente explode?

Qualquer resposta vale um milhão de dólares.

Aqui está onde estamos. Ninguém provou que soluções suaves sempre existem globalmente em 3D, e ninguém construiu uma explosão também. Estamos presos no meio desde o artigo de Leray de 1934, mais de noventa anos de uma das questões abertas mais difíceis em toda a matemática, e ainda não sabemos de que lado a resposta cai.

Curto prazo? Tudo bem. Para dados iniciais suaves, as equações produzem uma solução suave por algum período de tempo. O fluido começa a se mover, a matemática funciona, tudo está limpo. Mas o que acontece depois? A solução permanece suave para sempre, ou ela atinge um ponto onde a velocidade dispara para o infinito?

Se ela permanece suave, algo legal acontece. Essa solução suave automaticamente satisfaz as regras relaxadas para soluções fracas também, então é uma solução fraca. E pela unicidade fraca-forte, nenhuma outra solução fraca com essas condições iniciais pode existir. Então, se alguém provasse suavidade global, toda a hierarquia colapsaria: fraca, forte e suave acabariam sendo a mesma coisa, uma única solução única que é perfeitamente bem comportada para todo o tempo. É isso que torna este problema tão atraente e tão difícil. A lacuna entre o que podemos provar que existe (soluções fracas) e o que queremos (soluções suaves) é exatamente o conteúdo da questão de um milhão de dólares.

O Problema do Milênio de Clay (Fefferman 2000) pergunta: para qualquer $u_0 \in C^\infty(\mathbb{R}^3)$ livre de divergência satisfazendo $|\partial^\alpha u_0(x)| \leq C_{\alpha K}(1 + |x|)^{-K}$ para todos $\alpha, K$, existe $u \in C^\infty(\mathbb{R}^3 \times [0,\infty))$ e $p \in C^\infty(\mathbb{R}^3 \times [0,\infty))$ satisfazendo as equações de Navier-Stokes com $\int_{\mathbb{R}^3} |u(x,t)|^2 \, dx < C$ para todo $t \geq 0$?

A alternativa (também digna de prêmio): encontrar $u_0$ na classe acima tal que nenhuma solução suave desse tipo exista.

O que é conhecido:

• Existência local: Para $u_0 \in H^s(\mathbb{R}^3)$, $s \geq 1/2$, uma solução mild/forte local única existe em $[0,T^*)$ para algum $T^* > 0$ dependendo de $\|u_0\|_{H^s}$, e é suave para todo tempo positivo $t > 0$. Se $T^* < \infty$, então $\|u(t)\|_{H^s} \to \infty$ quando $t \to T^*$.
• Existência global para dados pequenos: Se $\|u_0\|_{\dot{H}^{1/2}}$ (ou $\|u_0\|_{L^3}$, ou $\|u_0\|_{BMO^{-1}}$) é menor que uma constante universal, a solução é global e suave. Referências chave: $\dot{H}^{1/2}$ (Fujita-Kato 1964), $BMO^{-1}$ (Koch-Tataru 2001), com resultados análogos em $L^3$ (Kato 1984).
• Unicidade fraca-forte: Se uma solução forte existe em $[0,T]$, então toda solução fraca de Leray-Hopf com os mesmos dados iniciais coincide com ela em $[0,T]$. Isso foi estabelecido por Serrin (1962) e refinado por trabalhos subsequentes. Significa que provar regularidade também resolve unicidade dentro da classe de Leray-Hopf.

A hierarquia colapsa para cima: suave $\Rightarrow$ forte $\Rightarrow$ fraca, e unicidade fraca-forte significa que uma solução suave, se existe, é a única solução fraca de Leray-Hopf. Então o Problema do Milênio é equivalente a perguntar: todas as soluções fracas de Leray-Hopf com dados iniciais suaves são elas mesmas suaves? A lacuna entre "solução fraca existe" (Leray 1934) e "solução suave existe" (aberto) é exatamente a questão do prêmio.

Se soluções fracas existem e descrevem o fluido, por que alguém deveria se importar com suavidade?

Três razões.

Primeiro, unicidade. A física exige uma resposta. Me dê o estado inicial de um fluido, e eu deveria poder dizer exatamente o que ele faz a seguir. Não "aqui estão várias possibilidades, escolha a que você quiser". Mas soluções fracas não garantem isso. Múltiplas soluções fracas podem emergir do mesmo fluxo inicial sem nenhuma maneira de dizer qual delas o fluido real segue. As equações se tornariam um menu em vez de uma receita. Isso não é física.

Segundo, confiabilidade numérica. Muitas simulações importantes de fluidos são baseadas em Navier-Stokes ou modelos intimamente relacionados: previsões meteorológicas, aerodinâmica, fluxo sanguíneo através de artérias, e mais. Torne a grade mais fina e a simulação deveria convergir para a resposta verdadeira. Sem uma garantia de suavidade e unicidade? Nenhum teorema diz que isso realmente acontece em todo cenário 3D. As simulações funcionam. Não podemos explicar completamente por quê.

Terceiro, física extrema. Se singularidades podem se formar, isso é a natureza nos enviando uma mensagem. As equações de Navier-Stokes, nosso melhor modelo de movimento de fluidos, teriam uma data de validade embutida: em alguma escala extrema o próprio modelo para de funcionar, e as equações estão nos dizendo, "Você precisa de nova física".

Isso não é uma tecnicalidade. É a linha de falha que atravessa tudo. Existência de um lado (Leray, 1934, feito). Suavidade do outro (aberto, um milhão de dólares). Por que cruzar é tão difícil? As estimativas de energia em 3D ficam apenas um pouco aquém do que é necessário, e noventa anos de esforço, centenas de artigos, carreiras inteiras gastas tentando, ninguém fechou essa lacuna.

Cada estratégia de prova sendo perseguida agora é uma tentativa de preencher essa divisão. Prove que soluções fracas são suaves. Ou prove que não são. Duas palavras ou três.

A hierarquia de soluções para as equações de Navier-Stokes incompressíveis 3D é:

$$\text{smooth} \subset \text{strong} \subset \text{Leray-Hopf weak} \subset \text{distributional weak}$$

O que é conhecido em cada nível:

O Problema do Milênio está na lacuna entre existência fraca global de Leray-Hopf e regularidade suave global. A dificuldade central: a desigualdade de energia fornece $u \in L^\infty_t L^2_x \cap L^2_t \dot{H}^1_x$, que está meia derivada abaixo do escalonamento crítico $\dot{H}^{1/2}$. Preencher essa lacuna de supercriticalidade é equivalente a provar regularidade global.

Unicidade fraca-forte implica que a hierarquia colapsa se soluções suaves existem globalmente. Mais precisamente: se para dado $u_0 \in C^\infty$ uma solução suave $u$ existe em $[0,T]$, então toda solução fraca de Leray-Hopf com os mesmos dados é igual a $u$ em $[0,T]$. Então suavidade global $\Rightarrow$ unicidade global dentro da classe de Leray-Hopf.

Inversamente, não unicidade de soluções fracas de Leray-Hopf para dados iniciais suaves admissíveis implicaria que soluções suaves globais não podem persistir para essa classe de dados (já que uma solução suave forçaria unicidade). Trabalho recente sobre integração convexa (baseando-se em De Lellis-Székelyhidi para Euler, estendido por Buckmaster-Vicol 2019 para construir soluções fracas não únicas de Navier-Stokes abaixo da regularidade de Leray-Hopf) mostra que soluções fracas distribucionais podem ser altamente não únicas. Se essa não unicidade se estende à classe de Leray-Hopf é uma questão aberta importante com implicações diretas para o Problema do Milênio.

Para o estado atual das estratégias de prova atacando essa lacuna, e para as razões estruturais pelas quais é tão resistente, veja Por Que Navier-Stokes É Difícil.