Parte 1As equações de Navier–Stokes

Toda vez que você vê fumaça subindo de um incenso, água girando no ralo ou um furacão na imagem de satélite, o movimento obedece às mesmas duas equações, escritas por Claude-Louis Navier (1822) e George Gabriel Stokes (1845).

Primeiro: o que é um fluido para a matemática?

Em vez de acompanhar trilhões de moléculas, tratamos o fluido como um meio contínuo. Em cada ponto do espaço x e em cada instante t existe uma velocidade u(x, t) (uma setinha: direção e rapidez), uma pressão p(x, t) e uma densidade ρ. Esse conjunto de setas é o campo de velocidade — ligue “Setas do campo de velocidade” no simulador para vê-lo.

Derivada material: acompanhando uma gotinha

Juntando os dois primeiros termos obtemos a derivada material, a aceleração sentida por quem “viaja” junto com o fluido:

Mesmo num escoamento que não muda com o tempo (∂u/∂t = 0), uma gota que gira no redemoinho está sempre acelerando, porque sua direção muda o tempo todo.

Parte 2Como o computador resolve as equações

Não existe uma fórmula que dê a resposta pronta para um escoamento qualquer. Então o computador divide a tela numa grade de células e avança o tempo em passos pequenos (cerca de 60 por segundo), usando o método Stable Fluids de Jos Stam (1999), muito usado em efeitos visuais e jogos.

Cada passo aplica os termos da equação um de cada vez — a técnica se chama divisão de operadores (operator splitting). Os cálculos rodam na placa de vídeo (GPU), que processa todas as células em paralelo.

O truque mais importante: a projeção

Depois dos passos 1–3, o campo de velocidade u* em geral não tem divergente zero (ele “cria” ou “some” com fluido). O teorema de Helmholtz–Hodge garante que qualquer campo pode ser separado em uma parte sem divergente mais o gradiente de algo. Esse “algo” é justamente a pressão:

Ou seja: a pressão é a “força” que o fluido precisa fazer para continuar incompressível. É uma ideia profunda — a pressão não tem equação própria; ela se ajusta instantaneamente para impor ∇·u = 0.

Um detalhe sobre o ralo

Numa simulação 2D (vista de cima), a água que desce pelo ralo vai para uma terceira dimensão que não existe na tela. Para representar isso, a projeção impõe ∇·u < 0 no ralo (fluido sumindo) e um valor positivo minúsculo espalhado pela superfície inteira, que representa o nível da água baixando por igual. O total se conserva.

Parte 3Anatomia de um redemoinho

Vorticidade: a medida do giro

A vorticidade é o rotacional da velocidade:

Imagine um cata-ventinho minúsculo solto no fluido: a vorticidade diz quão rápido ele gira (é o dobro da velocidade angular local). Na visualização “Vorticidade ω”, vermelho é giro anti-horário e azul é horário. Repare que quase toda a vorticidade fica concentrada no núcleo: fora dele o fluido gira em volta do centro, mas o cata-vento não gira em torno de si mesmo!

A circulação Γ é a soma da vorticidade dentro de uma região — o “tamanho total” do redemoinho. É o controle principal do simulador.

O modelo usado: vórtice de Lamb–Oseen

O motor do redemoinho empurra o fluido para esta velocidade tangencial, que é uma solução exata das equações de Navier–Stokes:

Perto do centro (r ≪ r_c) isso vira v ≈ Ω r: o núcleo gira como um disco rígido. Longe do centro vira v ≈ Γ/(2πr): um vórtice livre, que gira mais devagar quanto mais longe. O gráfico abaixo do simulador compara essa teoria com a velocidade realmente medida na simulação — um ótimo teste de que o programa obedece à física.

O padrão de corante “Raios” mostra isso muito bem: no núcleo os raios giram retos (disco rígido); fora dele se transformam em espirais, porque o fluido de dentro dá voltas mais rápido que o de fora (rotação diferencial).

Por que a superfície afunda no meio?

Para uma gota andar em círculo, algo precisa puxá-la para o centro. Esse algo é a diferença de pressão (balanço ciclostrófico):

A pressão cresce para fora, então no centro ela é mínima. Num líquido com superfície livre, pressão baixa significa coluna de água mais baixa: a superfície afunda formando o famoso “funil”. No modo Água, o sombreamento da superfície é calculado a partir desse campo de pressão. No olho de um furacão acontece o mesmo: é a região de menor pressão atmosférica.

A bailarina e o ralo: conservação do momento angular

Quando uma patinadora fecha os braços, ela gira mais rápido. No ralo acontece a mesma coisa: a sucção traz fluido de longe para perto do centro e, conservando o momento angular (r·v ≈ constante), um giro fraquinho vira um redemoinho forte. Teste: cenário “Ralo da pia”, zere Γ, dê um pequeno giro com o mouse e aumente a sucção.

2D × 3D: o esticamento de vórtices

Tomando o rotacional de Navier–Stokes, obtém-se a equação da vorticidade. Em 3D ela tem um termo a mais:

O termo destacado é o esticamento de vórtices: quando um tubo de vórtice é esticado (como um tornado sendo puxado para cima, ou o funil do ralo), ele fica mais fino e gira mais rápido. Em 2D esse termo é sempre zero, e é por isso que a matemática 2D está resolvida. A grande pergunta é: em 3D, esse esticamento pode ficar infinito em tempo finito?

Parte 4Modo 3D: o redemoinho em volume e os objetos

No modo 3D o computador resolve Navier–Stokes numa grade de 64 × 64 × 64 = 262.144 células. Aqui aparece o ingrediente que falta no plano: o esticamento de vórtices.

Como cabe um volume numa placa de vídeo?

A GPU é ótima para processar imagens. Então o volume é “fatiado” em 64 camadas horizontais, e cada camada vira um ladrilho de uma única imagem de 512 × 512 pixels (um atlas). Cada etapa do método (forças, advecção, viscosidade, projeção) processa todas as 262 mil células de uma vez. Em 3D, cada célula conversa com 6 vizinhas (esquerda, direita, cima, baixo, frente, trás) e a vorticidade deixa de ser um número e vira um vetor.

O esticamento em ação: vórtice de Burgers

No tanque, um jato puxa a água para baixo pelo eixo (controle “🚀 Acelerador da descida”); perto do fundo ela se abre para os lados, sobe pelas paredes e volta — a água recircula, como num tanque com hélice. Isso estica o tubo de vórtice: ele fica mais fino e gira mais rápido (a bailarina fechando os braços, agora na vertical). Existe uma solução exata de Navier–Stokes em 3D que descreve exatamente esse equilíbrio entre esticamento e viscosidade, o vórtice de Burgers:

O raio do núcleo é √(4ν/a): quanto mais forte o esticamento a, mais fino e rápido fica o vórtice. A viscosidade é o que impede que ele fique infinitamente fino. O Problema do Milênio pergunta, no fundo, se em algum escoamento o esticamento pode vencer a viscosidade e produzir uma velocidade infinita.

A superfície em funil

A grade 3D tem uma “tampa” fixa no topo. A altura da superfície é calculada pela pressão logo abaixo dela: η = p/(ρ g) (aproximação de tampa rígida). Pressão baixa no eixo = superfície afundada. O controle “Exagero da superfície” amplia esse efeito, que na realidade é pequeno.

Por que cada objeto faz uma coisa diferente?

Os objetos obedecem a uma versão simplificada da equação de Maxey–Riley, usada por cientistas para estudar partículas em escoamentos:

• (u − v)/τ — arrasto: o fluido tenta arrastar o objeto para a sua própria velocidade. τ pequeno (folha, bolha) = obedece rápido; τ grande (aço) = tem “teimosia” (inércia).
• β Du/Dt — a mesma diferença de pressão que acelera o fluido também empurra o objeto. No redemoinho, ela aponta para o centro.
• (1 − β) g — peso menos empuxo, já incluindo a massa adicional (o objeto precisa empurrar água para se mover).

O paradoxo das folhas de chá

Mexa uma xícara de chá: as folhas no fundo vão para o centro, e não para a borda como numa centrífuga. Einstein explicou isso em 1926: perto do fundo, o atrito freia a água, a força centrífuga fica fraca ali e a pressão empurra a água do fundo para o centro (escoamento secundário). O tanque do simulador tem esse atrito no fundo, mas a bolinha de aço é pesada demais (7,8 vezes a densidade da água) e a “centrífuga” vence: ela vai para a borda. Uma folha de chá de verdade é quase tão densa quanto a água, por isso o escoamento secundário consegue levá-la ao centro.

E no tornado?

Na câmara de tornado, um exaustor no teto puxa o ar para cima enquanto o ar entra girando pela base (como na câmara de Ward, usada em laboratórios). No chão, o vento forte produz sustentação (efeito Bernoulli: ar rápido por cima = pressão menor em cima; é assim que telhados são arrancados). No ar, o objeto sobe se a corrente ascendente for mais rápida que a sua velocidade terminal g·τ. A folha, o isopor e o patinho (τ pequeno) sobem; o cubo de madeira só pula e o aço nunca sai do chão. A poeira (partículas) mostra o caminho do ar: entra rente ao chão e sobe em espiral pelo funil.

Parte 5Fumaça × Líquido: mesmas equações, fluidos diferentes

Ar e água obedecem exatamente às mesmas equações. O que muda são os números (densidade e viscosidade) e algumas forças extras.

Surpresa: o ar é “mais grudento” que a água em termos de ν, a viscosidade que aparece na equação. A água tem mais atrito (μ), mas é tão pesada que o atrito quase não consegue freá-la.

O que é específico da fumaça

A fumaça é ar quente carregando partículas. Ar quente é menos denso e sobe (empuxo). Na aproximação de Boussinesq, isso entra como uma força externa proporcional à temperatura: f = β (T − T₀) g. No cenário “Fumaça subindo”, a pluma quente sobe, fica instável e se enrola em redemoinhos em forma de cogumelo (instabilidade de Kelvin–Helmholtz).

O que é específico dos líquidos

Líquidos têm superfície livre, que afunda onde a pressão é baixa (o funil do ralo), e tintas se misturam muito devagar porque a difusão molecular é lenta — por isso a tinta desenha filamentos tão finos.

Número de Reynolds: o número mais importante da mecânica dos fluidos

Re compara o termo azul com o termo verde da equação. Re pequeno: escoamento calmo, organizado (laminar). Re grande: caótico (turbulento). Para um redemoinho usa-se Re_Γ = Γ/ν; para o escoamento ao redor de um obstáculo, Re_D = U·D/ν. O simulador mostra o valor no canto superior direito.

Parte 6O Problema do Milênio (e o anúncio da OpenAI)

Em 2000, o Clay Mathematics Institute escolheu sete grandes problemas em aberto e ofereceu US$ 1 milhão por cada solução. Um deles é sobre as equações desta página.

A pergunta

Em três dimensões, começando com um escoamento perfeitamente suave e com energia finita, a solução das equações de Navier–Stokes continua suave para sempre? Ou pode acontecer uma explosão (blow-up): em um tempo finito, a velocidade ou a vorticidade ficarem infinitas num ponto?

Fisicamente parece absurdo — água de verdade não fica com velocidade infinita. Mas ninguém conseguiu provar que as equações proíbem isso. Se a explosão existir, significa que Navier–Stokes deixa de descrever o fluido naquele ponto.

O enunciado oficial (escrito por Charles Fefferman) aceita quatro tipos de resposta: provar que as soluções são sempre suaves (A e B, sem força externa, em todo o espaço ou num domínio periódico) ou mostrar um exemplo de explosão (C e D), e nestes últimos é permitido usar uma força externa suave f.

Linha do tempo

O que a OpenAI anunciou

Segundo a própria empresa e reportagens da imprensa:

• Um sistema com pelo menos 10 mil agentes de IA trabalhou por 88 horas (a partir de 1º de setembro de 2026) e produziu um artigo de 166 páginas, acompanhado de uma verificação formal no assistente de provas Lean.
• O resultado é uma construção de explosão com força externa suave — ou seja, uma resposta do tipo (C) do enunciado do Clay.
• A OpenAI declarou que não vai reivindicar o prêmio de US$ 1 milhão.

Como é o escoamento que “explode” (e como o simulador o mostra)

Importante: a OpenAI não publicou uma simulação de computador — publicou uma prova matemática. Mas a prova descreve com precisão a forma do escoamento que explode, e pesquisadores já o analisaram (por exemplo, o artigo arXiv 2609.17642). O cenário 3D “Colapso da OpenAI” mostra esse perfil com as mesmas leis de escala:

• Um vórtice que se contrai sozinho. Tudo é escrito em termos de τ = T − t, o tempo que falta para a explosão. As variáveis de similaridade da prova são X = r²/(2τ) e η = z/τ^1/2−h, com um expoente pequeno 0 < h < 1/100.
• Ele fica cada vez mais fino: o raio encolhe como τ^1/2 e o comprimento como τ^1/2−h — a largura diminui mais depressa que o comprimento.
• Espiral para dentro, saída pelo eixo: o fluido entra girando e sai ao longo do eixo, dividido perto de z = 0 (parte sobe, parte desce), com o giro valendo zero no próprio eixo.
• A força age só num anel em volta do núcleo, em pulsos oscilantes (fase proporcional a log X). Ela fornece exatamente o que falta para o perfil obedecer às equações.
• O resultado: a velocidade cresce como |u| ~ τ^−1/2−h → ∞, enquanto a energia do núcleo diminui como τ^1/2−3h → 0 e o número de Reynolds cresce devagar, Re ~ τ^−h. É uma explosão da velocidade com energia finita.

No simulador, o perfil exato da prova (que é calculado numericamente) é trocado por fórmulas fechadas com o mesmo comportamento — giro do tipo Burgers com raio √τ e um escoamento de entrada e saída exatamente incompressível. As leis de escala, a geometria e o papel da força são os da construção. O tempo avança em escala logarítmica, como o próprio colapso autossimilar. Na análise do arXiv, os autores observam que um fluido real cavitaria (água) ou formaria ondas de choque (ar) antes de chegar à explosão.

Por que o resultado é polêmico

• Força “artificial”: para muitos matemáticos, a pergunta realmente importante é se o fluido pode explodir sozinho, sem força externa. A força usada é matematicamente válida pelo enunciado, mas fisicamente artificial. Como resumiu um especialista à Scientific American, o problema do Clay pode estar resolvido, mas o problema principal das equações não. Depois do anúncio, outros matemáticos mostraram que o método não funciona sem essa força.
• Disputa de autoria: cerca de 12 horas antes do anúncio, o matemático Tristan Buckmaster afirmou que ideias recentes dele e de Levent Alpöge (sobre as equações de Euler) teriam vazado para a OpenAI e inspirado os comandos dados aos agentes.
• Verificação: até as últimas notícias, a comunidade ainda não tinha concluído uma verificação independente completa da construção.

Parte 7Roteiro para a feira de ciências

Apresentação em ~6 minutos

1. Gancho (20 s): abra o cenário Ralo da pia. “Isto é água descendo pelo ralo — calculada, não filmada. O computador está resolvendo as equações de Navier–Stokes 60 vezes por segundo.”
2. A equação (60 s): mostre a barra da equação. Explique que é F = m·a para cada gotinha. Desligue a pressão: o fluido fica “irreal”. Religue.
3. Como funciona o giro (60 s): mude para Pressão (centro com pressão baixa) e depois para Vorticidade (giro concentrado no núcleo). Mostre o gráfico: teoria × medido.
4. Fluidos diferentes (40 s): Fumaça subindo (empuxo e redemoinhos em cogumelo) e depois a água trocada por Mel (viscosidade alta mata o giro).
5. 3D: o colapso da OpenAI (60 s): clique em 3D · volume + objetos. O primeiro cenário mostra o vórtice que encolhe e gira cada vez mais rápido até “explodir”. Aponte o gráfico: a velocidade vai ao infinito enquanto a energia vai a zero, e os anéis laranja são a força externa da polêmica.
6. 3D com objetos (40 s): no “Redemoinho no tanque”, clique em “🎲 Um de cada” e aumente o 🚀 acelerador. Pergunte ao público: “quem desce primeiro?” A bolha vai para o eixo e o aço é jogado para fora. Termine com o tornado levantando a folha, mas não a bolinha de aço.
7. Beleza da física (30 s): Esteira de von Kármán — os mesmos vórtices que aparecem atrás de ilhas em fotos de satélite e fazem fios cantarem com o vento.
8. Fechamento (40 s): o Problema do Milênio, US$ 1 milhão, e o anúncio da OpenAI em setembro de 2026 — com a polêmica da força externa e a lembrança de que a prova ainda está sendo verificada.

Perguntas que os jurados podem fazer

Referências

• Jos Stam, “Stable Fluids”, SIGGRAPH 1999.
• Mark Harris, “Fast Fluid Dynamics Simulation on the GPU”, GPU Gems, cap. 38 (2004).
• Fedkiw, Stam & Jensen, “Visual Simulation of Smoke”, SIGGRAPH 2001 (confinamento de vorticidade).
• Charles Fefferman, “Existence and Smoothness of the Navier–Stokes Equation” — enunciado oficial do Clay Mathematics Institute (2000).
• OpenAI, “On the Navier–Stokes Millennium Prize Problem” (set. 2026) — openai.com/index/navier-stokes-solution
• CNBC, “OpenAI claims to have solved the 90-year-old Navier-Stokes math problem in 88 hours” (9 set. 2026) — cnbc.com
• Axios, “OpenAI's historic math solution overshadowed by credit controversy” (8 set. 2026) — axios.com
• Scientific American, “Did OpenAI solve the wrong Navier-Stokes problem?” — scientificamerican.com
• “Self-similar swirl between contracting porous walls: the GD1998 exact Navier–Stokes solution revisited in the similarity variables of the OpenAI 2026 forced blow-up construction” — arXiv 2609.17642 (variáveis de similaridade e leis de escala usadas no cenário 3D).
• K. Adler, visualização do colapso da OpenAI em 3D — github.com/keithadler/superfx-navier-stokes.
• Wikipédia, “Navier–Stokes existence and smoothness”.