Une aiguille d'eau : comment 10 000 agents ont fait exploser les équations de Navier–Stokes ✛ ✎ ✦ OpenAI théorique réalisé par des agents d'IA d'OpenAI (≈10 000 agents, 88 h) · vérifié en Lean

Original: "Finite Time Blowup for Navier–Stokes"
· OpenAI
b42w:2609.90001 · page du travail · PDF, 166 p. · formalisation Lean · 2026-09-08 · © OpenAI · analyse indépendante · 8 min
Depuis le repos, avec une force lisse et une énergie finie : un écoulement dont la vitesse tend vers l'infini en temps fini. Comment il a été construit et ce que cela signifie pour le problème du Millénaire.
Illustration d'OpenAI tirée de l'annonce : spirale vers l'intérieur et étirement axial. © OpenAI
Illustration d'OpenAI tirée de l'annonce : spirale vers l'intérieur et étirement axial. © OpenAI
Notre schéma du mécanisme : le fluide spirale vers l'axe et s'écoule le long de lui, le cœur se contracte en aiguille, et dans l'anneau auto
Notre schéma du mécanisme : le fluide spirale vers l'axe et s'écoule le long de lui, le cœur se contracte en aiguille, et dans l'anneau autour du cœur les impulsions oscillatoires produisent une contrainte moyenne qui remplace la force manquante. À l'extérieur, un écoulement purement azimutal qui n'a besoin d'aucune force.

En 1934, Jean Leray a prouvé que les équations de Navier–Stokes ont toujours des solutions « au sens généralisé ». Si la solution reste lisse ou peut se briser en une singularité, il ne l'a pas su, et personne ne l'a su pendant quatre-vingt-dix ans. En 2000, l'Institut Clay a fait de la question l'un de ses sept problèmes du Millénaire. L'énoncé écrit par Charles Fefferman est honnête : il demande la preuve de l'une quelconque de quatre affirmations. (A) et (B) : sans force extérieure, la régularité est toujours conservée, dans l'espace et sur le tore. (C) et (D) : il existe une force lisse et des données initiales lisses pour lesquelles il n'y a pas de solution lisse d'énergie finie.

Le travail d'OpenAI prouve (C) et (D). Pour toute viscosité \(\nu\gt 0\), il construit une force lisse, à support compact dans l'espace et le temps, et un écoulement qui part du repos complet, garde son énergie bornée et atteint une vitesse infinie à \(t=1\). La façon dont c'est fait mérite d'être comprise, parce que la méthode retourne le problème.

Le problème est résolu à l'envers. L'équation de Navier–Stokes peut se lire comme une définition de la force. Écrivez n'importe quelle vitesse u et pression p, substituez-les et déclarez force extérieure f tout ce qui ne s'équilibre pas. L'équation est vérifiée identiquement. Un écoulement qui explose est facile à inventer. La difficulté se déplace en un seul endroit : le résidu doit être lisse. Ses termes individuels divergent, et il faut organiser des compensations pour que la somme et toutes ses dérivées traversent calmement l'instant de l'explosion. C'est à cela que sont consacrées les 166 pages.

Le tourbillon-aiguille. Un tourbillon sur l'axe vertical s'effondre vers l'origine. À l'intérieur, le fluide spirale vers l'axe et s'écoule le long de lui, vers le haut et vers le bas. La conservation du moment cinétique entraîne l'accélération : pour une parcelle sans couple, \(r\,u_\theta\) est constant, donc un r plus petit signifie un \(u_\theta\) plus grand. L'incompressibilité entretient l'afflux : le fluide qui arrive est emporté par l'écoulement axial, sinon il boucherait l'axe.

Le cœur se contracte de façon autosimilaire : le profil garde sa forme et seules les échelles changent. Soit \(\tau=1-t\) le temps restant avant l'explosion. Le rayon du cœur est \(\sim\tau^{1/2}\), la hauteur \(\sim\tau^{1/2-h}\) avec un minuscule \(h\lt 1/100\), les vitesses \(\sim\tau^{-1/2-h}\). Le rayon rétrécit plus vite que la hauteur, et le cœur s'étire en une aiguille fine de volume \(\sim\tau^{3/2-h}\). De là vient le point le plus élégant du travail : l'énergie cinétique du cœur \(\sim\tau^{1/2-3h}\) tend vers zéro. La vitesse est infinie, l'énergie s'évanouit, le volume décroît plus vite que ne croît le carré de la vitesse.

La viscosité n'est pas vaincue mais intégrée. La viscosité est ce qui sépare Navier–Stokes d'Euler, et c'est à cause d'elle qu'on pensait l'équation peut-être sûre. Le travail calcule deux nombres de Reynolds. En rotation, \(\mathrm{Re}_\theta\sim\tau^{-h}\to\infty\) : le fluide fait de plus en plus de tours en un temps de diffusion visqueuse, la rotation gagne. Dans l'afflux radial, \(\mathrm{Re}_r=O(1)\) : la viscosité reste un partenaire à égalité. Les rythmes de tous les processus, transport axial et radial et diffusion à travers le cœur, croissent pareillement comme \(\tau^{-1}\) et restent dans le bilan dominant jusqu'à la fin. Le long de l'axe, la diffusion est plus faible d'un facteur \(\tau^{2h}\), et la construction utilise cette marge.

Là où ça casse. Il faut coudre le cœur au fluide calme de l'extérieur. Dans l'anneau entre cœur et extérieur, la couture ne s'équilibre pas, et le résidu croît sans borne quand \(t\to1\). La force doit être lisse, or ici elle devrait être infinie. Le tourbillon de fond ne se soutient pas seul.

Le fluide produit la force lui-même. Dans l'anneau, on plante de petites perturbations oscillantes, chacune un anneau complet autour de l'axe, localisées en rayon, hauteur et temps et amorcées par une graine exponentiellement petite. Comment une oscillation de moyenne nulle donne-t-elle une force non nulle ? En sortant, le fluide emporte un excédent de vitesse azimutale ; en rentrant, un déficit. La vitesse moyenne est nulle, mais le flux de moment cinétique, produit des deux, garde son signe sur les deux demi-cycles. Inversez les deux signes et le produit ne change pas. Il en résulte un flux moyen de quantité de mouvement stable, dont la variation spatiale agit sur le fond comme une force véritable. Les physiciens appellent cela contrainte de Reynolds et la traitent d'ordinaire comme une gêne ; ici c'est un outil, réglé pour annuler exactement la partie divergente du résidu.

Les impulsions s'allument et s'éteignent d'elles-mêmes. Elles croissent sur le cisaillement du fond par un mécanisme centrifuge : une parcelle avec un excédent de rotation est poussée vers l'extérieur, et si la vitesse angulaire décroît assez vite avec le rayon, son excédent par rapport à ses voisines croît. Les deux effets se nourrissent, donnant une croissance exponentielle tant qu'elle bat l'amortissement visqueux. Le cisaillement les tue aussi : il emporte les différents rayons à des vitesses différentes et raccourcit la longueur d'onde radiale, et l'amortissement visqueux croît comme le carré du nombre d'onde. La longueur d'onde initiale est choisie pour que l'amplification gagne d'abord et l'amortissement, inévitablement, ensuite. Chaque impulsion laisse une queue si petite que toutes ses dérivées s'annulent à l'instant de l'explosion.

Le reste de l'assemblage. Il y a deux familles d'impulsions : il faut fournir à la fois le flux de moment cinétique et le flux de quantité de mouvement axiale, et les familles diffèrent par le rapport des deux. Vers l'explosion, les impulsions sont posées sur des grilles de plus en plus fines en espace et en temps, comme en intégration convexe : une échelle infinie de corrections, chacune plus fine que la précédente. Il y a une asymétrie délibérée : avec une symétrie exacte haut-bas, le mécanisme d'amplification s'annule dans le plan \(z=0\), si bien que le profil axial est légèrement incliné. Hors de l'anneau, l'écoulement est purement azimutal, indépendant de la hauteur, et résout exactement l'équation radiale de la chaleur : son résidu est nul et aucune force n'y est nécessaire. Cela permet une troncature spatiale lisse, et la compacité porte le résultat sur le tore, d'où (D).

Ce qui est nouveau par rapport aux prédécesseurs. La stratégie de forçage et d'amplification à travers les échelles appartient à Córdoba et Martínez-Zoroa : d'abord pour Euler 3D, puis avec Zheng pour Navier–Stokes à dissipation affaiblie. Tao a construit en 2016 une explosion pour un Navier–Stokes moyenné qui conserve la compensation énergétique. Deux pas sont nouveaux ici : une vraie viscosité et l'équation complète pour tout \(\nu\gt 0\), et des impulsions qui ne se contentent pas de s'amplifier mais servent de source de la contrainte requise, ajustée exactement sous le trou du fond. Les briques sont toutes connues : instabilité centrifuge, transport des vecteurs d'onde par le cisaillement, ondes de cisaillement visqueuses exactes, réalisation des contraintes par des oscillations issue des constructions de Daneri et Székelyhidi pour Euler. Seul l'assemblage est nouveau.

Qui l'a fait, et comment. D'après OpenAI, la preuve a été trouvée par un groupe d'environ dix mille agents coordonnés sur un modèle interne, en 88 heures, en envoyant 2,7 millions de messages et en dépensant environ 130 milliards de jetons de sortie ; la formalisation et la vérification en Lean ont pris 17 heures de plus. Aucun humain ne signe le travail : l'auteur indiqué est OpenAI. La veille, Tristan Buckmaster et Levent Alpöge avaient publié trois prépublications utilisant la même approche du forçage lisse pour l'équation des milieux poreux, Boussinesq et Euler 3D, avec Lean aussi. Une dispute de priorité a surgi ; de l'aveu des deux parties, les preuves diffèrent, et les résultats pour Euler eux-mêmes diffèrent : forcé d'un côté, non forcé de l'autre.

Ce qui reste ouvert. Partir du repos sonne spectaculaire mais signifie le contraire de ce que l'on lit : puisque le fluide était immobile, toute l'énergie est venue de la force extérieure. L'explosion ici n'est pas un « ordre né de rien » mais un pompage calculé avec une grande précision. Le problème sans force, alternatives (A) et (B), n'est pas résolu, l'Institut Clay le liste comme ouvert, et OpenAI dit ne pas réclamer le prix. Mais le schéma où de vieux outils sont menés jusqu'au bout par une machine et vérifiés par une autre vient d'être essayé sur un problème qui a résisté aux humains pendant quatre-vingt-dix ans.

Le regard de la machine du savoir

Cette section s'adresse à l'auteur du travail et à ses collègues. Les travaux voisins ont été trouvés par une recherche sémantique dans notre archive, et non par correspondance de mots. Nous n'évaluons pas : nous montrons le champ autour.

Ce que fait ce travail

Un écoulement est construit dont la vitesse part à l'infini en temps fini sous une force extérieure lisse, depuis le repos et à énergie finie. La singularité vit dans une aiguille qui rétrécit plus vite qu'elle n'accélère, et la force manquante est fournie par le fluide lui-même via la contrainte moyenne d'impulsions oscillatoires. Pour nous, c'est le premier exemple rigoureusement construit d'une singularité locale d'un milieu continu à énergie finie, un objet avec lequel continuer à travailler.

Où est sa force

Le mécanisme est explicite : tourbillon-aiguille, anneau, deux familles d'impulsions, extérieur de chaleur. Il peut donc être transféré : calculer sa métrique effective, demander ce qu'en fait une viscosité qui s'annule, tenter de le reproduire en laboratoire.

Vers où avancer

À essayer

  • Dériver la métrique acoustique de l'écoulement construit et vérifier si sa courbure diverge quand t → 1.
  • Calculer la production locale d'entropie sur les trois zones de la construction et chercher une région à bilan négatif.
  • Reproduire le mécanisme des impulsions annulaires dans un modèle de Gross–Pitaevskii à viscosité nulle.
  • Monter une cuve tournante à entraînements contrarotatifs et mesurer la contrainte de Reynolds dans l'anneau.

Pourquoi cela peut compter

Si une singularité d'écoulement se transfère à une métrique effective, les fluides ordinaires et quantiques deviennent un banc d'essai pour des questions intestables dans la gravité elle-même : comment l'entropie se comporte près d'une singularité et si l'ordre peut tenir là où la description continue se brise. Le travail d'OpenAI donne à ce banc, pour la première fois, une singularité explicite, rigoureusement construite, à énergie finie.

Ce que nous en pensons

Ce travail n'a pas été évalué par les pairs. Voici notre lecture.

Le point fort

Retourner le problème : l'écoulement est écrit, la force est définie comme résidu, et toute la difficulté est honnêtement déplacée dans la régularité du résidu. Les impulsions auto-limitées, qui s'allument sur le cisaillement et s'éteignent par viscosité, avec leur contrainte moyenne ajustée à la divergence du fond, forment un mécanisme complet et vérifiable, non une borne supérieure.

Ce qui le renforcerait

  • Donner une exposition humaine des Sections 4–9 en moins de cent pages : la physique est claire en six pages, mais le chemin de là aux estimations n'est ouvert que formellement.
  • Dire explicitement ce que couvre la formalisation Lean (le Théorème 1.1 en entier ou des estimations particulières) et comment texte et code se correspondent : sans cela, « vérifié en Lean » se lit plus large que ce n'est peut-être.
  • Une illustration numérique de l'écoulement construit sur un horizon τ fini montrerait les échelles et les impulsions de façon visible et donnerait un contrôle indépendant du choix de h et des longueurs d'onde.
  • Discuter quelle partie de la construction est obstruée par l'absence de force : c'est la seule route de (C) vers (A).

Questions à l'auteur

  • Quelle est la stabilité de la construction sous perturbations du fond et des impulsions : l'explosion persiste-t-elle dans un voisinage ouvert des données construites ?
  • Le mécanisme d'auto-alimentation des impulsions donne-t-il un indice pour le cas non forcé, où la contrainte devrait naître des données initiales ?
  • Qu'est-ce qui est formalisé exactement en Lean : le Théorème 1.1 complet avec supports et estimations, ou son corollaire de non-existence ?

Aucun humain n'a signé le travail : la ligne d'auteur indique OpenAI. Environ dix mille agents ont trouvé la preuve en 88 heures, et elle a été vérifiée pendant 17 heures de plus non par une personne, mais par Lean.

f := \partial_t u + (u\cdot\nabla)u - \nu\Delta u + \nabla p
L'écoulement est écrit, la force est calculée. L'équation est vraie par construction ; la tâche est de rendre f lisse alors que u → ∞.
\ell_r\asymp\tau^{1/2},\ \ell_z\asymp\tau^{1/2-h},\ |u|\asymp\tau^{-1/2-h},\ E\asymp\tau^{1/2-3h}
Quand τ → 0, la vitesse diverge et l'énergie du cœur tend vers zéro : l'aiguille rétrécit plus vite qu'elle n'accélère.
\langle u'_r u'_\theta\rangle \neq 0
Une oscillation de moyenne nulle transporte du moment cinétique dans une direction. C'est ainsi que les impulsions fournissent la force qui manquait au fond.

Emporter le travail complet

Non publié sur arXiv. Voici sa page, le PDF et la formalisation : vous pouvez vérifier vous-même, pas seulement lire notre analyse.

Scientifiques
Isaac NewtonJohn MichellRobert HookeWolfgang PauliDavid CharbonneauJohannes Kepler
Notions
Équations de Navier-Stokes machine à mouvement perpétuel écoulement vitesse de libération force de marée principe d'Archimède principe de Bernoulli traînée singularité résistance nulle superfluidité vorticité instabilité de Kelvin-Helmholtz vortex superfluide lévitation par énergie du vide troisième principe de la thermodynamique flux obstacle vitesse