Abordagens ao Problema de Navier-Stokes

A abordagem mais antiga começa com energia. Um fluido em movimento carrega energia cinética, e a viscosidade a consome, como o atrito que faz as coisas pararem. A energia total só pode diminuir ao longo do tempo, assumindo que nada está bombeando energia de fora.

Leray viu isso em 1934 e fez um movimento-chave: usar o limite de energia para provar que uma solução fraca global com energia cinética finita tem que existir. Construir soluções aproximadas, artificialmente suavizadas. Mostrar que todas obedecem ao limite de energia. Tomar um limite. Algo deve sobreviver nesse limite, e sobrevive.

Mas aqui está o problema. Limites de energia são instrumentos grosseiros. Eles garantem que o fluido tem energia total finita, claro, mas não podem dizer se a velocidade permanece finita em cada ponto do espaço e do tempo. Essa lacuna entre "energia finita" e "suave em toda parte" é exatamente o problema de regularidade, e está em aberto há noventa anos.

Links de artigos: Leray, Sur le mouvement d'un liquide visqueux emplissant l'espace (1934); Hopf, Über die Anfangswertaufgabe für die hydrodynamischen Grundgleichungen (1951).

A construção de Leray-Hopf procede em duas variantes principais: aproximação de Galerkin (projetando em subespaços de dimensão finita) ou mollificação de Friedrichs (suavizando a não-linearidade via convolução). Ambas compartilham o mesmo esqueleto de quatro etapas. A estratégia padrão tem quatro passos:

1. Aproximar: Resolver o sistema mollificado $\partial_t u_\varepsilon + (J_\varepsilon u_\varepsilon \cdot \nabla) u_\varepsilon = \nu \Delta u_\varepsilon - \nabla p_\varepsilon$ em subespaços de dimensão finita.
2. Limite de energia: A estimativa a priori $\|u_\varepsilon(t)\|_{L^2}^2 + 2\nu \int_0^t \|\nabla u_\varepsilon\|_{L^2}^2 \leq \|u_0\|_{L^2}^2$ vale uniformemente em $\varepsilon$.
3. Compacidade: Extrair uma subsequência fracamente convergente $u_\varepsilon \rightharpoonup u$ em $L^2_t \dot{H}^1_x$ usando o lema de Aubin-Lions.
4. Passar ao limite: O termo não-linear converge por convergência forte $L^2_{\text{loc}}$ de $u_\varepsilon$.

A solução fraca resultante satisfaz a desigualdade de energia, não a igualdade. A energia pode ser perdida em tempos irregulares. E esse é todo o problema: a lacuna entre a classe de energia $L^\infty_t L^2_x \cap L^2_t \dot{H}^1_x$ e a suavidade real é precisamente o que não conseguimos fechar.

Links de artigos: Leray, Sur le mouvement d'un liquide visqueux emplissant l'espace (1934); Hopf, Über die Anfangswertaufgabe für die hydrodynamischen Grundgleichungen (1951).

A abordagem de Caffarelli-Kohn-Nirenberg (1982) não tenta provar suavidade completa. Ela pergunta algo completamente diferente: quão ruins as singularidades podem realmente ser?

Quase nada ruins. Seu teorema de $\varepsilon$-regularidade diz que se certas quantidades locais invariantes por escala são pequenas o suficiente em uma pequena região espaço-temporal, a solução é automaticamente suave ali. E como a energia total é finita, simplesmente não há "orçamento" suficiente para muitos pontos singulares coexistirem.

Pense assim. Uma parede pode ter rachaduras. Mas o comprimento total de todas essas rachaduras combinadas é zero, significando que o conjunto singular é extremamente pequeno no sentido de medida parabólica teórica (medida de Hausdorff parabólica unidimensional zero).

Links de artigos: Caffarelli-Kohn-Nirenberg, Partial regularity of suitable weak solutions of the Navier-Stokes equations (1982); Albritton-Barker-Prange, Epsilon regularity for the Navier-Stokes equations via weak-strong uniqueness.

O teorema de $\varepsilon$-regularidade CKN: para soluções fracas adequadas, a pequenez de certas quantidades locais invariantes por escala (envolvendo tanto $u$ quanto $p$ em um cilindro parabólico $Q_r(z_0) = B_r(x_0) \times (t_0 - r^2, t_0)$) implica regularidade (continuidade Hölder) perto do ponto $z_0 = (x_0, t_0)$.

A prova combina a desigualdade de energia local com uma iteração tipo Campanato: se a energia invariante por escala é pequena, um argumento de bootstrap mostra que $u$ é limitada, depois Hölder, depois suave pela teoria clássica de Schauder.

A estimativa dimensional $\mathcal{P}^1(\Sigma) = 0$ segue por um argumento de cobertura de Vitali. Pontos singulares forçam cotas inferiores em quantidades localizadas invariantes por escala de dissipação/energia; a cobertura então controla a medida parabólica unidimensional do conjunto singular por uma medida local de energia, que se anula em pequenas escalas.

Links de artigos: Caffarelli-Kohn-Nirenberg, Partial regularity of suitable weak solutions of the Navier-Stokes equations (1982); Albritton-Barker-Prange, Epsilon regularity for the Navier-Stokes equations via weak-strong uniqueness.

Aqui está uma redução precisa de todo o problema. Beale, Kato e Majda provaram em 1984 que para as equações de Euler 3D, explosão só pode acontecer se o controle de vorticidade for perdido. Critérios análogos foram posteriormente estabelecidos para Navier-Stokes. É isso. Uma condição.

A vorticidade mede rotação local. O critério BKM diz: mantenha a rotação máxima limitada na norma certa, e a solução permanece suave. Todo o resto se alinha automaticamente.

Uma família de quantidades para controlar. Infelizmente, realmente controlá-las revelou-se exatamente tão difícil quanto o problema original. A redução é limpa. A execução permanece fora de alcance.

Links de artigos: Beale-Kato-Majda, Remarks on the breakdown of smooth solutions for the 3-D Euler equations (1984); Kozono-Taniuchi, Bilinear estimates in BMO and the Navier-Stokes equations (2000).

O teorema original de Beale-Kato-Majda (1984) é para as equações de Euler 3D. Para Navier-Stokes, critérios de continuação análogos implicam que uma solução suave $u$ em $[0, T^*)$ se estende além de $T^*$ sempre que

$$\int_0^{T^*} \|\omega(\cdot, t)\|_{L^\infty} \, dt < \infty,$$

onde $\omega = \nabla \times u$ é a vorticidade. Refinamentos incluem:

• Kozono-Taniuchi (2000): $\|\omega\|_{L^\infty}$ pode ser substituído por $\|\omega\|_{\mathrm{BMO}}$
• Variantes em espaços de Besov: controle crítico ou limítrofe em Besov também pode servir como critério de continuação
• Critérios com restrição de direção: condições tipo Serrin em componentes de $\nabla u$ também podem servir como critérios de continuação (ver por exemplo Beirão da Veiga, 1995)

Esses critérios conectam-se à imagem de estiramento de vórtice: qualquer singularidade em tempo finito deve forçar a vorticidade a acumular-se rápido demais para que a integral temporal acima permaneça finita.

Links de artigos: Beale-Kato-Majda, Remarks on the breakdown of smooth solutions for the 3-D Euler equations (1984); Kozono-Taniuchi, Bilinear estimates in BMO and the Navier-Stokes equations (2000); Chemin-Planchon, Self-improving bounds for the Navier-Stokes equations (2012).

Um ângulo mais moderno trabalha com espaços de funções (como $L^3$ ou $\dot{H}^{1/2}$) que ficam exatamente na fronteira do que a simetria de escala permite. Estes são espaços críticos, e é onde vivem os resultados de regularidade mais precisos.

A lógica é limpa: se você pode mostrar que uma solução permanece dentro de certos limites de espaço crítico, a suavidade segue automaticamente. Múltiplas equipes provaram isso, construindo um menu completo de critérios de regularidade (condições que garantem suavidade se você puder verificá-las).

O problema é a lacuna. Métodos de energia dão controle abaixo do nível crítico, em um regime supercrítico em energia. Precisamos de limites críticos. Essa lacuna é estreita, às vezes uma única derivada de regularidade, mas resistiu a todas as tentativas de fechá-la.

Links de artigos: Koch-Tataru, Well-posedness for the Navier-Stokes equations (2001); Kenig-Koch, An alternative approach to regularity for the Navier-Stokes equations in critical spaces; Gallagher-Koch-Planchon, A profile decomposition approach to the $L^\infty_t(L^3_x)$ Navier-Stokes regularity criterion. Para uma comparação detalhada de por que a criticalidade de energia funciona em 2D mas falha em 3D, veja Por Que 2D É Mais Fácil Que 3D.

Principais programas em regularidade de espaço crítico:

• Koch-Tataru (2001): Boa colocação global para dados pequenos em $\text{BMO}^{-1}$. A estimativa-chave é um limite de ponto fixo para a aplicação bilinear de Duhamel no espaço de soluções de fluxo de calor de Koch-Tataru construído sobre dados em $\text{BMO}^{-1}$, não uma estimativa estática nua de produto em $\text{BMO}^{-1}$. Este é um dos resultados mais precisos conhecidos para métodos perturbativos.
• Gallagher-Koch-Planchon (2013): Abordagem de decomposição de perfil ao critério de regularidade $L^\infty_t L^3_x$ de Navier-Stokes. Qualquer sequência de soluções com norma crítica limitada tem uma subsequência decompondo-se em perfis assintoticamente desacoplados.

A obstrução central: nenhum funcional coercivo conhecido é tanto controlado pela evolução quanto crítico em relação ao escalonamento de Navier-Stokes.

Links de artigos: Koch-Tataru, Well-posedness for the Navier-Stokes equations (2001); Kenig-Koch, An alternative approach to regularity for the Navier-Stokes equations in critical spaces; Gallagher-Koch-Planchon, A profile decomposition approach to the $L^\infty_t(L^3_x)$ Navier-Stokes regularity criterion.

A teoria moderna de EDP empresta muito da análise harmônica. A ideia central: quebrar uma função em ondas em diferentes frequências, da mesma forma que você dividiria um acorde musical em notas individuais. Exceto que aqui, as "notas" são oscilações espaciais da velocidade do fluido em escalas muito diferentes.

A decomposição de Littlewood-Paley faz exatamente isso. Cortar o campo de velocidade em componentes escala por escala. Rastrear como a energia flui entre elas. De repente, a intuição física informal de "cascata de energia" torna-se algo sobre o qual você pode realmente provar teoremas, e os teoremas são precisos. Esses métodos produziram muitos dos resultados mais precisos sobre critérios de regularidade e taxas de explosão.

Links de artigos: Cannone-Meyer, Littlewood-Paley decomposition and Navier-Stokes equations (1995); Gallagher-Koch-Planchon, A profile decomposition approach to the $L^\infty_t(L^3_x)$ Navier-Stokes regularity criterion.

A teoria de Littlewood-Paley decompõe $u = \sum_j \Delta_j u$ onde $\Delta_j$ localiza em frequências $|\xi| \sim 2^j$. Aplicada à não-linearidade de transporte $(u \cdot \nabla)u$:

A decomposição de paraproduto da não-linearidade $(u \cdot \nabla)u$ divide-se em interações de frequência baixa-alta, alta-baixa e alta-alta:

$$(u \cdot \nabla)u = T_u \nabla u + T_{\nabla u} u + R(u, \nabla u)$$

onde $T$ é o paraproduto e $R$ o resto. Cada peça tem propriedades de regularidade diferentes em espaços de Besov $\dot{B}^s_{p,q}$.

Resultados-chave usando essa maquinaria:

• Espaços de Chemin-Lerner: $\widetilde{L}^\rho_T \dot{B}^s_{p,q}$ fornecem a estrutura natural para boa colocação crítica: a forma bilinear de Navier-Stokes mapeia $\widetilde{L}^\infty_T \dot{B}^{-1+3/p}_{p,q} \times \widetilde{L}^1_T \dot{B}^{1+3/p}_{p,q} \to \widetilde{L}^1_T \dot{B}^{-1+3/p}_{p,q}$.
• Cannone-Meyer: Métodos de Littlewood-Paley dão uma formulação limpa em wavelet/Besov da teoria de dados pequenos.

Links de artigos: Cannone-Meyer, Littlewood-Paley decomposition and Navier-Stokes equations (1995); Gallagher-Koch-Planchon, A profile decomposition approach to the $L^\infty_t(L^3_x)$ Navier-Stokes regularity criterion.

Aqui está um instinto completamente diferente. Em vez de rastrear números (normas, energias), esses métodos estudam a forma da solução: como tubos de vórtice se dobram, como regiões de rotação intensa se organizam no espaço.

A percepção-chave é que explosão não é apenas sobre algo ficar grande. É sobre o fluido se organizar em uma configuração geométrica muito específica. Se você pode mostrar que essa configuração é impossível (porque contradiz a estrutura de dissipação de energia, ou incompressibilidade, ou ambas), você descartou a explosão sem nunca calcular uma norma.

Esse ponto de vista geométrico cresceu em uma perspectiva que inspirou vários critérios de regularidade rigorosos ao lado de métodos puramente analíticos. E parece diferente. Pergunta que forma o desastre toma? em vez de quão grande esse número pode ficar?

Links de artigos: Constantin-Fefferman, Geometric constraints on potentially singular solutions for the 3-D Euler equations (1993); Albritton-Barker-Prange, Localized smoothing and concentration for the Navier-Stokes equations in the half space.

Abordagens geométrico-topológicas exploram restrições estruturais que são invisíveis a métodos puramente analíticos:

• Geometria de linha de vórtice: Constantin e Fefferman (1993) mostraram que se o campo de direção de vorticidade $\hat{\omega} = \omega/|\omega|$ é Lipschitz em regiões de alta vorticidade, a solução é regular. Explosão requer que a direção de vorticidade desenvolva uma singularidade simultaneamente com a magnitude.
• Argumentos de incompatibilidade: se uma configuração de explosão é geometricamente restrita (por exemplo, via limites de empacotamento no número de regiões de concentração independentes que cabem dentro de orçamentos de energia e dissipação), pode-se derivar uma contradição sem estimar diretamente normas críticas.
• Partição de casos (especulativa/programática): uma estratégia proposta classificaria cada região espacial como pertencente a um de finitos cenários (por exemplo, localmente regular, tipo Tipo-I, tipo Tipo-II, densamente empacotado) e tentaria mostrar que cada cenário ou dá regularidade ou transfere o problema para um argumento de contagem limitada. Isso permanece um programa de pesquisa em vez de um resultado estabelecido.

A página de prova neste site explora argumentos desse tipo; nada aqui é apresentado como uma prova formal completa do problema do Milênio.

Links de artigos: Constantin-Fefferman, Geometric constraints on potentially singular solutions for the 3-D Euler equations (1993); Albritton-Barker-Prange, Localized smoothing and concentration for the Navier-Stokes equations in the half space.

Esta pegou as pessoas de surpresa. As soluções fracas do método de Leray (Seção 1) revelam-se não-únicas, pelo menos quando forçamento externo está presente.

A arma é a integração convexa, uma técnica originalmente construída para problemas de geometria e adaptada para equações de fluidos por De Lellis e Székelyhidi começando por volta de 2009. A ideia: construir soluções "selvagens" iterativamente empilhando correções de alta frequência que coletivamente satisfazem a equação mas se comportam erraticamente.

Para Navier-Stokes 3D, Buckmaster e Vicol (2019) provaram não-unicidade de soluções fracas de energia finita abaixo da classe de Leray-Hopf. A não-unicidade de Euler vem do programa de integração convexa de De Lellis-Székelyhidi, Isett e Buckmaster. Então em 2022, Albritton, Brué e Colombo provaram que mesmo soluções de Leray-Hopf de Navier-Stokes 3D são não-únicas quando força externa está presente. Se a não-unicidade persiste para as equações de Navier-Stokes sem forçamento permanece em aberto.

Por que isso importa? Porque "uma solução fraca existe" tem sido o resultado principal desde 1934. Agora sabemos que não determina uma única resposta. A questão se aguça: qual solução, se alguma, é a fisicamente correta?

Links de artigos: De Lellis-Székelyhidi, Dissipative continuous Euler flows (2013); Buckmaster-Vicol, Nonuniqueness of weak solutions to the Navier-Stokes equation (2019); Albritton-Brué-Colombo, Non-uniqueness of Leray solutions of the forced Navier-Stokes equations (2022).

A integração convexa para equações de fluidos origina-se no programa De Lellis-Székelyhidi (2009–2013), adaptando a técnica de imersão isométrica $C^1$ de Nash-Kuiper para construir soluções fracas das equações de Euler que dissipam energia. Estágios-chave:

• De Lellis-Székelyhidi (2013): existência de fluxos de Euler dissipativos contínuos ($C^0$) em $\mathbb{T}^3$ (regularidade Hölder abaixo de $1/5$ alcançada no trabalho subsequente de Buckmaster-De Lellis-Isett-Székelyhidi 2015; posteriormente melhorada para $< 1/3$ por Isett, 2018, resolvendo o lado flexível da conjectura de Onsager).
• Buckmaster-Vicol (2019): não-unicidade de soluções fracas de energia finita para Navier-Stokes 3D, com soluções em $C^0_t H^\beta_x$ para algum $\beta > 0$. A construção usa fluxos de Beltrami intermitentes como blocos de construção, adicionando correções oscilatórias em cada passo de iteração enquanto mantém controle do tensor de Reynolds. Isso está abaixo da classe de energia de Leray-Hopf, então não contradiz diretamente a unicidade de Leray.
• Albritton-Brué-Colombo (2022): não-unicidade de soluções de Leray-Hopf para as equações de Navier-Stokes 3D forçadas. A prova constrói uma solução auto-similar instável de fundo e usa um mecanismo de instabilidade para ramificar em soluções de Leray-Hopf distintas a partir dos mesmos dados iniciais. Isso mostra que a desigualdade de energia sozinha não seleciona uma solução única quando forçamento está presente.

A questão central em aberto é se a não-unicidade persiste para as equações de Navier-Stokes sem forçamento na classe de Leray-Hopf. O resultado forçado mostra que a desigualdade de energia não é um princípio de seleção suficiente, mas não resolve se a equação sem forçamento tem estrutura adicional que restaura a unicidade.

Links de artigos: De Lellis-Székelyhidi, Dissipative continuous Euler flows (2013); Buckmaster-Vicol, Nonuniqueness of weak solutions to the Navier-Stokes equation (2019); Albritton-Brué-Colombo, Non-uniqueness of Leray solutions of the forced Navier-Stokes equations (2022).

Podemos pelo menos descartar certas estratégias de prova? Terence Tao mostrou em 2016 que sim, podemos. E o resultado é sóbrio.

Tao construiu uma versão modificada das equações de Navier-Stokes, um sistema "médio", que mantém muitas características estruturais-chave das equações reais: a identidade de energia, a forma como a enstrofia (uma medida de intensidade de vorticidade) cresce, a simetria de escala. Mas neste sistema modificado, soluções explodem em tempo finito.

A implicação descarta amplas famílias de estratégias de prova. Qualquer prova de que suavidade global vale para as equações reais deve usar alguma propriedade estrutural específica da verdadeira não-linearidade que o sistema médio não tem. Você não pode provar regularidade usando apenas limites de energia, escala e crescimento de enstrofia. Essas ferramentas sozinhas são consistentes com explosão.

Isso não diz que as equações reais explodem. Diz que famílias inteiras de estratégias de prova são becos sem saída. A prova eventual (se regularidade vale) deve ser mais precisa que um argumento genérico de energia. Muito mais precisa.

Links de artigos: Tao, Finite time blowup for an averaged three-dimensional Navier-Stokes equation (2016).

Tao (2016) considera um sistema da forma

$$\partial_t u + \tilde{B}(u, u) = \nu \Delta u - \nabla p, \quad \nabla \cdot u = 0,$$

onde $\tilde{B}$ é um operador bilinear que concorda com a verdadeira não-linearidade de Navier-Stokes $(u \cdot \nabla)u$ nos seguintes sentidos:

• é um multiplicador de Fourier de ordem 1, preservando o escalonamento $u(x,t) \mapsto \lambda u(\lambda x, \lambda^2 t)$
• satisfaz a mesma identidade de energia: $\langle \tilde{B}(u,u), u \rangle = 0$
• reproduz a estrutura de crescimento de enstrofia

Para este sistema médio, Tao constrói dados iniciais suaves cuja solução explode em tempo finito. O mecanismo de explosão programa uma sequência de concentrações de enstrofia cada vez mais agudas em escalas cada vez menores, com cada estágio dobrando a enstrofia em uma cascata controlada.

A obstrução que isso cria: qualquer quantidade supercrítica que é (a) controlada pela evolução na classe de energia e (b) invariante sob o escalonamento de Navier-Stokes não pode por si só descartar explosão, porque também seria controlada no sistema médio, que explode. Uma prova de regularidade deve explorar a estrutura específica de cancelamento algébrico do verdadeiro termo de transporte $(u \cdot \nabla)u$ que o operador médio $\tilde{B}$ não compartilha.

Links de artigos: Tao, Finite time blowup for an averaged three-dimensional Navier-Stokes equation (2016).

Este artigo é parte de Progresso.

Da prova de existência de Leray de 1934 através da integração convexa e barreiras de prova de Tao, estas são as principais estratégias que as pessoas lançaram ao problema de Navier-Stokes 3D. Nenhuma resolveu o problema completo de regularidade 3D. Para contexto sobre como a viscosidade molda a matemática comparada às equações de Euler invíscidas, veja Euler vs. Navier-Stokes. Para o status atual, veja O Problema de Navier-Stokes Foi Resolvido? Para a declaração formal exata, retorne a O Problema do Milênio.

Este artigo é parte de Progresso.

As abordagens acima representam os principais fios rigorosos na literatura de regularidade e unicidade até 2022. Nenhuma combinação resolveu o problema 3D completo. O campo continua avançando. Métodos de prova assistidos por computador são uma área ativa coberta separadamente neste site.

Para a comparação entre sistemas viscosos e invíscidos (e por que a viscosidade ajuda mas não o suficiente), veja Euler vs. Navier-Stokes. Para os subproblemas que essas abordagens visam, veja Subproblemas. Para os obstáculos de escala, veja Por Que É Difícil.