main concepts scientists sections
SimplePopulaireAvancé

Explosion en temps fini pour les équations de Navier–Stokes tridimensionnelles avec force lisse à support compact : anatomie de la construction ✛ ✎ ✦ 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 · 7 min
Analyse des EDP Dynamique des fluides
Pour toute viscosité : une force lisse à support compact et une solution depuis le repos, d'énergie finie, dont la vitesse devient infinie à t = 1, réglant les alternatives (C) et (D) de Fefferman. Schéma de la preuve, échelles, rôle des impulsions et limites du résultat.
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.

Contexte

Le problème de régularité demande si une solution des équations de Navier–Stokes incompressibles tridimensionnelles issue de données lisses reste lisse. Leray (1934) a construit des solutions faibles globales d'énergie finie ; Caffarelli, Kohn et Nirenberg (1982) ont montré que l'ensemble singulier d'une solution faible convenable a une mesure de Hausdorff parabolique unidimensionnelle nulle, ce qui autorise des singularités isolées ; Escauriaza, Seregin et Šverák (2003) ont donné la régularité sous une norme \(L^\infty_t L^3_x\) invariante d'échelle bornée. Tao (2016) a construit une explosion pour des équations moyennées conservant la compensation énergétique ; Buckmaster et Vicol (2019) ont obtenu la non-unicité des solutions faibles par intégration convexe ; Albritton, Brué et Colombo (2022) ont construit des solutions de Leray–Hopf distinctes avec données nulles et une même force, singulière à l'instant initial. Les prédécesseurs directs sont Córdoba et Martínez-Zoroa : explosion pour Euler 3D forcé (2023), pour les milieux poreux avec source lisse (2024) et, avec Zheng, pour Navier–Stokes hypodissipatif. Dans leurs constructions, la déformation à grande échelle amplifie la vorticité à petite échelle avec des auto-interactions supprimées. La veille du travail d'OpenAI, Buckmaster et Alpöge ont publié des explosions à forçage lisse pour les milieux poreux, Boussinesq 2D et Euler 3D avec formalisations Lean ; Tao a jugé le même jour leur approche « susceptible de s'étendre à Navier–Stokes ». L'énoncé de Fefferman accepte explicitement l'une quelconque de (A)–(D) comme solution ; la condition (5) sur la force est satisfaite pour toute \(f\in C_c^\infty\).

Méthodes

La méthode est inverse : le résidu \(R(u,p)=\partial_t u+(u\cdot\nabla)u-\Delta u+\nabla p\) (\( u\) = 1 après changement d'échelle) est déclaré force, et la tâche devient de construire u, p à vitesse non bornée et résidu lisse jusqu'à \(t=1\) avec toutes ses dérivées spatio-temporelles. La construction est en couches. (1) Fond axisymétrique autosimilaire \((u_B,p_B)\) : en coordonnées cylindriques, un cœur avec afflux radial, accélération azimutale et sortie axiale de part et d'autre d'une couche de séparation près de \(z=0\) ; échelles \(\ell_r\asymp\tau^{1/2},\ \ell_z\asymp\tau^{1/2-h},\ 0\lt h\lt 1/100\), vitesses \(|u_\theta|,\,|u_z|\asymp\tau^{-1/2-h},\ |u_r|=O(\tau^{-1/2})\) ; la pression décroît vers l'axe et fournit la force centripète ; la viscosité transporte le moment cinétique vers l'extérieur, et la croissance de la vitesse est le bilan de l'afflux contre la perte visqueuse. Les équations de profil imposent le bilan dominant dans le cœur ; le raccord à un extérieur lisse laisse dans l'anneau un résidu égal à moins la divergence d'une contrainte annulaire, à un terme près qui s'annule à tout ordre au point singulier (Proposition 5.5). Le profil axial est délibérément un peu asymétrique, avec un léger biais vers le haut et \(u_z(r,0,t)\neq0\) : avec une symétrie exacte, le transport axial de moment cinétique et le cisaillement radial de la vitesse axiale s'annuleraient tous deux en \(z=0\), et les impulsions ne pourraient s'y amplifier. (2) Impulsions oscillatoires : suites de perturbations spatialement oscillantes localisées en r, z, t et fermées en anneaux complets, deux familles aux rapports différents entre flux radial de moment cinétique et flux radial de quantité de mouvement axiale ; leurs produits quadratiques moyennés réalisent la contrainte dominante (Proposition 7.5). La graine est une force extérieure exponentiellement petite ; la croissance se fait sur le cisaillement du fond par le mécanisme centrifuge (critères de Leibovich–Stewartson, Billant–Gallaire) avec une contribution du cisaillement axial ; les orientations sont choisies pour que le cisaillement raccourcisse la longueur d'onde radiale, que l'amortissement \(\propto k^2\) l'emporte sur l'amplification, et que les queues avec toutes leurs dérivées s'annulent quand \(t\to1\). La construction s'appuie sur les solutions exactes de Craik–Criminale et les ondes de cisaillement visqueuses de Singh–Sridhar, et sur la technique de Daneri–Székelyhidi de réalisation des contraintes par des oscillations. (3) Des corrections successives améliorent la décroissance du résidu (Proposition 9.6) ; un tore auxiliaire et la séparation des supports oscillatoires (Section 6) ; des corrections moyennes à support compact (Section 8). (4) L'écoulement extérieur est purement azimutal, indépendant de z, et résout exactement l'équation radiale de la chaleur (« extérieur de chaleur ») : son résidu s'annule identiquement ; pour tout r > 0 fixé, les champs et toutes leurs dérivées ont des limites lisses quand \(t\uparrow1\), ce qui permet une troncature spatiale lisse (Section 10) ; la force de troncature reste lisse. Les estimations contrôlent toutes les dérivées cartésiennes du résidu. La formalisation en Lean 4 avec Mathlib (dépôt openai/NavierStokesAndEuler, Apache-2.0) couvre la non-existence d'une solution lisse globale à énergie bornée sur R³ et sur le tore.

Résultats

Pour tout \(\nu\gt 0\) : une solution lisse à support compact et énergie finie depuis le repos avec \(\|u(t)\|_{L^\infty}\to\infty\) quand \(t\uparrow1\), sous une force \(f\in C_c^\infty(\mathbb R^3\times(0,\infty))\). Échelles clés : volume du cœur \(\asymp\tau^{3/2-h}\), énergie cinétique du cœur \(\asymp\tau^{1/2-3h}\to0\) avec vitesse non bornée ; nombres de Reynolds \(\mathrm{Re}_\theta\asymp\tau^{-h}\to\infty\) et \(\mathrm{Re}_r=O(1)\) ; les rythmes \(|u_r|/\ell_r,\ |u_z|/\ell_z,\ \nu/\ell_r^2\) sont tous \(\asymp\tau^{-1}\), la diffusion axiale plus faible d'un facteur \(\tau^{2h}\). Corollaire 10.6 : transport sur T³. Les alternatives (C) et (D) de Fefferman sont revendiquées. OpenAI signale en outre une explosion pour les équations d'Euler non forcées (travail séparé, ~100 agents, ~50 heures). D'après OpenAI : ~10 000 agents coordonnés, 88 heures jusqu'à la solution, 17 heures de formalisation et vérification en Lean, 2,7 millions de messages, ~130 milliards de jetons de sortie sur ce problème.

Implications

Premier exemple d'explosion en temps fini pour les équations de Navier–Stokes tridimensionnelles complètes avec vraie viscosité et force lisse, depuis le repos et à énergie finie. Le mécanisme n'est pas un « dépassement » de la viscosité mais son inclusion dans un bilan autosimilaire à nombre de Reynolds de rotation divergent. Il est montré de façon constructive que le terme non linéaire peut lui-même fournir la contrainte manquante par les flux quadratiques moyens d'oscillations, autrement dit la contrainte de Reynolds agit comme un outil contrôlable et non comme une gêne. Pour le problème (A)/(B), le résultat ne donne pas de réponse directe : la force est essentielle et l'énergie est injectée de l'extérieur. Mais il déplace la frontière : la question sans force devient celle de savoir si une contrainte analogue peut naître des seules données initiales, sans pompage. Le fait méthodologique est à part entière significatif : une preuve de 166 pages assemblée à partir d'outils connus par un système d'agents et vérifiée formellement ; c'est le premier résultat de ce poids publié sans signature humaine.

Développements futurs

Directions naturelles : tenter de porter la construction au cas non forcé en remplaçant le pompage par une énergie stockée dans les données initiales ; étudier la stabilité de la construction sous perturbations et la taille de h ; une illustration numérique de l'écoulement construit ; une exposition humaine indépendante des Sections 4–9 ; une comparaison avec l'approche de Buckmaster–Alpöge pour voir quelles parties coïncident en substance. Pour les équations hypodissipatives et les modèles, le résultat simplifiera probablement les constructions existantes. Reste ouverte la question de savoir si un tel pompage est physiquement réalisable : la force est mathématique, compacte et lisse, mais sa contrepartie physique n'est pas discutée.

Impact

Dynamique des fluides mathématique (régularité, solutions faibles, intégration convexe), théorie des instabilités de tourbillons, vérification formelle des mathématiques, méthodologie de la recherche par IA. Indirectement, dynamique des fluides numérique : la construction fournit un cas test exact de concentration à vitesse divergente et énergie évanescente.

Prochaines étapes

Lire la Section 2 (image physique) et la Section 3 (schéma) de l'original ; confronter les conditions (4)–(7) à l'énoncé de Fefferman ; consulter le dépôt Lean et son README ; lire le billet de Tao du 7 septembre sur les travaux de Buckmaster–Alpöge et les prépublications de Córdoba–Martínez-Zoroa (arXiv:2309.08495, 2410.22920).

Problèmes ouverts clés

Le problème du Millénaire sur l'existence et la régularité de Navier–Stokes : alternatives (C) et (D) réglées, (A) et (B) ouvertes. Connexes : régularité d'Euler (OpenAI signale une explosion non forcée ; Buckmaster–Alpöge une explosion à forçage lisse), le critère de Beale–Kato–Majda, l'explosion conjecturée d'Euler en temps fini depuis des données lisses (indices numériques de Luo–Hou, preuve de Chen–Hou avec bord).

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 devient infinie 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. La viscosité n'est pas supprimée mais intégrée à un bilan autosimilaire. Pour nous, ce n'est pas seulement une réponse aux alternatives (C) et (D) de l'énoncé de Fefferman, mais le premier exemple rigoureusement construit d'une singularité locale d'un milieu continu à énergie finie : un objet avec lequel on peut continuer à travailler.

Où est sa force

Le mécanisme est entièrement explicite : tourbillon-aiguille, anneau, deux familles d'impulsions, extérieur de chaleur. Il peut donc non seulement être prouvé mais transféré : calculer sa métrique effective, demander ce qu'en fait une viscosité qui s'annule, tenter de le reproduire en laboratoire. Les travaux à construction explicite sont rares, et la valeur est précisément là.

Vers où avancer

À essayer

  • Dériver la métrique acoustique de l'écoulement construit et vérifier si sa courbure diverge quand t → 1 et comment cela s'accorde avec l'énergie évanescente du cœur.
  • Calculer la production locale d'entropie sur les trois zones de la construction, cœur, anneau, extérieur, et chercher une région à bilan négatif tenue par les impulsions.
  • Reproduire le mécanisme des impulsions annulaires dans un modèle de Gross–Pitaevskii à viscosité nulle : l'explosion survit-elle, et qu'est-ce qui remplace l'amortissement visqueux des impulsions ?
  • Monter une cuve tournante à entraînements contrarotatifs et mesurer la contrainte de Reynolds dans l'anneau en fonction de la vitesse : la façon la moins coûteuse de voir le mécanisme de ses propres mains.

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.

Travaux voisins dans notre archive

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 ?
Le graphe entier →

Le résidu du fond dans l'anneau de raccord diverge et ne peut être réparé par une force extérieure. C'est le fluide qui l'a réparé : des oscillations de vitesse moyenne nulle transportent de la quantité de mouvement dans une direction, et ce transport est ajusté exactement sous le trou. Ce que la théorie de la turbulence appelle une gêne est devenu un outil.

\partial_t u + (u\cdot\nabla)u - \nu\Delta u + \nabla p = f,\qquad \nabla\cdot u = 0
L'accélération de la parcelle (à gauche) égale la somme des forces : pression, frottement visqueux et force extérieure f. La seconde équation est l'incompressibilité. Dans le travail, f est lisse et à support compact ; u et p sont construits pour que l'équation soit vérifiée identiquement.
R(u,p) = \partial_t u + (u\cdot\nabla)u - \Delta u + \nabla p
Le geste central. Vitesse et pression sont écrites à la main et la force est déclarée égale au résidu R. L'équation est alors vraie par construction ; toute la tâche est de rendre R et toutes ses dérivées lisses quand t → 1 alors que u → ∞.
\ell_r \asymp \tau^{1/2},\qquad \ell_z \asymp \tau^{1/2-h},\qquad 0<h<\tfrac{1}{100}
Le rayon se contracte plus vite que la hauteur : ℓ_r/ℓ_z ≍ τ^h → 0. Le cœur s'étire en une aiguille fine de volume ~τ^{3/2−h}.
|u_\theta|,\,|u_z| \asymp \tau^{-1/2-h},\qquad E_{\text{core}} \asymp \tau^{1/2-3h} \to 0
La vitesse diverge tandis que l'énergie cinétique du cœur tend vers zéro : le volume décroît plus vite que ne croît le carré de la vitesse. La vitesse infinie vit dans une aiguille toujours plus fine.
\mathrm{Re}_\theta = \frac{|u_\theta|\,\ell_r}{\nu} \asymp \tau^{-h}\to\infty,\qquad \mathrm{Re}_r = \frac{|u_r|\,\ell_r}{\nu} = O(1)
La rotation dépasse la viscosité (Re_θ croît) tandis que l'afflux radial reste à son niveau (Re_r borné). La viscosité n'est pas vaincue : elle est intégrée au bilan dominant et y reste jusqu'à l'explosion.
r\,u_\theta = \text{const}\ \Rightarrow\ u_\theta \propto \frac{1}{r}
Une parcelle sans couple tourne plus vite en s'approchant de l'axe, comme une patineuse qui ramène ses bras. Dans le cœur réel, la viscosité emporte une partie du moment ; le bilan entre afflux et perte visqueuse fixe la croissance de la vitesse.
\langle u'_r\,u'_\theta\rangle \neq 0\quad\text{alors que}\quad \langle u'_r\rangle=\langle u'_\theta\rangle=0
Contrainte de Reynolds : si le fluide qui sort emporte un excédent de rotation et celui qui rentre un déficit, le produit garde son signe sur les deux demi-cycles. Ainsi les impulsions, sans vitesse moyenne, livrent une force moyenne qui bouche le trou du bilan de fond.

Chiffres clés

  • viscosité ν: quelconque, \(\nu\gt 0\)
  • instant d'explosion: \(t=1\), départ du repos u(·,0) = 0
  • rayon du cœur: \(\sim\tau^{1/2}\), \(\tau=1-t\)
  • hauteur du cœur: \(\sim\tau^{1/2-h}\), \(h\lt 1/100\)
  • vitesse dans le cœur: \(\sim\tau^{-1/2-h}\) → ∞
  • énergie du cœur: \(\sim\tau^{1/2-3h}\) → 0
  • nombre de Reynolds de rotation: \(\sim\tau^{-h}\to\infty\)
  • taille du travail: 166 pages + formalisation en Lean
  • agents dans le groupe: ≈ 10 000, 88 heures
  • vérification en Lean: 17 heures ; 2,7 M de messages, ≈130 milliards de jetons

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.

page du travailPDF, 166 p.formalisation Lean
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