Équations de Navier–Stokes

1. Pourquoi ce problème est important (problème du millénaire)

L’eau qui contourne un bateau et l’air qui s’écoule autour d’une aile obéissent au même bilan local : la masse se conserve et les forces modifient la quantité de mouvement. Les équations de Navier–Stokes traduisent cette idée dans un modèle de milieu continu. Écrire ces équations est bien plus facile que démontrer ce que toute solution peut faire.

État de la recherche vérifié le 13 octobre 2026. Le 8 septembre, OpenAI a annoncé une construction d’explosion en temps fini pour les équations de Navier–Stokes incompressibles en trois dimensions avec une force extérieure lisse, accompagnée d’un article et de code Lean. Elle vise les alternatives C et D du problème du millénaire de Clay. Il s’agit d’une avancée majeure revendiquée, mais elle n’établit pas l’explosion pour les équations de Navier–Stokes sans force extérieure. Le résultat distinct sur les équations d’Euler sans force extérieure concerne une viscosité nulle : il ne faut pas les confondre. Annonce d’OpenAI, mise à jour le 10 septembre.

La distinction entre un résultat annoncé, un énoncé formel vérifié, une démonstration lisible et l’acceptation par la communauté est ici essentielle. La section 14 présente l’état daté de l’évaluation, y compris les examens critiques d’octobre. Les notions physiques ci-dessous expliquent comment une énergie bornée peut coexister avec un mouvement local de plus en plus intense.


2. Champs continus et lois de conservation

Notons la masse volumique, la vitesse, la pression et le tenseur des contraintes de Cauchy. La force volumique par unité de masse est , par exemple la pesanteur.

Un champ attribue une valeur à chaque position et à chaque instant ; il représente une moyenne sur un grand nombre de molécules. Le modèle suppose que les échelles de longueur pertinentes sont suffisamment grandes pour que cette moyenne et une loi locale des contraintes aient un sens. Une singularité mathématique concerne les équations du milieu continu ; elle ne prédit pas qu’une molécule réelle atteigne une vitesse infinie.

2.1 Conservation de la masse

Le premier terme mesure l’accumulation locale de masse ; la divergence mesure le flux net de masse sortant. À masse volumique constante, la conservation de la masse se réduit à une divergence nulle de la vitesse. Une parcelle de fluide peut s’étirer et changer de forme tout en conservant son volume.

2.2 Conservation de la quantité de mouvement (Cauchy)

Décomposons les contraintes en une pression isotrope et une partie visqueuse :

La dérivée matérielle suit une parcelle en mouvement. Même un champ de vitesse stationnaire peut accélérer une parcelle : l’eau qui entre dans un passage plus étroit atteint une région où la vitesse est différente. C’est ce que décrit le terme de transport non linéaire .


3. Loi de comportement (fluide newtonien)

Pour un fluide newtonien isotrope,

où est la viscosité dynamique et le second coefficient de viscosité ; la viscosité volumique vaut .

Entre une plaque inférieure immobile et une plaque supérieure mobile, la vitesse d’un écoulement laminaire stationnaire croît linéairement avec la hauteur. Le schéma illustre le cisaillement visqueux et la condition d’adhérence.

Pour des coefficients de viscosité constants dans l’espace, leur substitution dans le bilan de quantité de mouvement donne les équations de Navier–Stokes compressibles :

On définit la viscosité cinématique lorsque est constante.

La viscosité répond à une déformation, pas à une translation uniforme. Ses dimensions sont et . Dans le schéma d’écoulement entre plaques, un écart et une vitesse de la plaque supérieure donnent et une contrainte de cisaillement . Cette solution exacte simple est lisse ; une annonce de singularité porte sur l’existence de données particulières, pas sur tous les écoulements.


4. Forme incompressible

Pour une masse volumique constante et

on obtient

L’équation indique que l’accélération résulte de la pression, de la diffusion visqueuse et de l’accélération imposée :

TermeSignification
Variation mesurée en un point fixe
Variation due au déplacement d’une parcelle dans le champ de vitesse
Les forces de pression redistribuent la quantité de mouvement
La viscosité atténue les variations spatiales
Force extérieure appliquée par unité de masse

L’incompressibilité n’impose ni une vitesse nulle ni une accélération nulle. La pression joue le rôle d’une contrainte qui maintient la divergence nulle de la vitesse ; dans ce modèle à masse volumique constante, elle n’est pas une loi d’état indépendante.

Prendre la divergence de l’équation de quantité de mouvement donne une équation de Poisson pour la pression, avec des conditions aux limites appropriées :


5. Vorticité et hélicité

Définissons la vorticité .
Pour un écoulement incompressible () :

En composantes :

Le terme représente l’étirement des tourbillons (en 3D seulement).
L’hélicité renseigne sur la topologie des lignes de vorticité.

Un tube tourbillonnaire étiré le long de son axe peut s’amincir et devenir plus intense. Dans l’équation de la vorticité, l’étirement peut amplifier la rotation tandis que la viscosité la diffuse. Aucun de ces effets, pris isolément, ne détermine si la dynamique couplée complète reste lisse.

Pour un écoulement véritablement plan , la vorticité est perpendiculaire au plan et les champs ne varient pas dans cette direction. Ainsi, . Cela supprime un mécanisme majeur d’amplification en deux dimensions ; cela ne signifie pas que tout tourbillon tridimensionnel doit exploser.

Une circulation dans le sens inverse des aiguilles d’une montre, dans le plan, possède une vorticité dirigée vers le lecteur, représentée par un point cerclé.


6. Bilan d’énergie (cas incompressible)

Pour une solution lisse, multiplions l’équation de quantité de mouvement par et intégrons sur un domaine dont la frontière est immobile et impose l’adhérence. La même identité vaut dans une boîte périodique ou dans tout l’espace avec une décroissance suffisante :

Pour , l’énergie cinétique décroît de façon monotone ; la dissipation visqueuse vaut .

Cette énergie est normalisée par la masse volumique constante : l’énergie cinétique physique vaut fois l’intégrale du carré de la vitesse. Le transport et la pression ne contribuent que par des flux aux frontières qui s’annulent sous ces hypothèses ; une intégration par parties donne le terme visqueux non négatif. Une paroi mobile ou une frontière ouverte ajoute des termes de travail ou de flux.

6.1 Pourquoi une énergie finie ne garantit pas la régularité

Une borne intégrale contrôle le mouvement total, pas sa plus grande valeur locale. Pour illustrer le changement d’échelle, choisissons un champ fixe , lisse, à support compact et à divergence nulle, puis posons

Alors croît lorsque , tandis que décroît. Ce sont des instantanés redimensionnés, pas une solution ni une preuve d’explosion. Ils montrent pourquoi des vitesses toujours plus grandes dans un volume toujours plus petit peuvent respecter une borne d’énergie. Le contrôle des dérivées et de la concentration constitue la difficulté supplémentaire du problème de régularité.


7. Adimensionnement et paramètres de contrôle

Choisissons des échelles caractéristiques (longueur) et (vitesse), puis définissons , , , .

Les équations incompressibles deviennent

avec le nombre de Reynolds

Ici, . Le temps de diffusion vaut et le temps d’advection ; leur rapport est . Un grand nombre de Reynolds indique que l’advection est rapide par rapport au lissage visqueux à l’échelle choisie. Ce n’est ni un seuil universel de turbulence ni un critère de singularité : la géométrie, les perturbations et l’échelle de mesure comptent aussi. Par exemple, , et donnent .

Autres nombres sans dimension, selon les besoins :

  • Mach (compressibilité),
  • Froude (gravité),
  • Prandtl (diffusion de quantité de mouvement comparée à la diffusion thermique).

8. Navier–Stokes–Fourier compressible (aperçu)

On ajoute l’énergie interne ou la température , avec le flux thermique de Fourier (la conductivité thermique est distincte de la diffusivité thermique ci-dessus) :

Une loi d’état, par exemple celle d’un gaz parfait, ferme le système : .

Le travail de compression et la conversion de l’énergie cinétique en énergie interne interviennent alors explicitement. Le problème du millénaire incompressible ne résout ni la formation des chocs, ni le transport de chaleur, ni la théorie de la régularité de ce système différent.


9. Théorie mathématique (existence et unicité)

L’existence locale signifie qu’une solution existe pendant un certain intervalle après l’instant initial ; l’existence globale signifie qu’elle existe pour tous les temps futurs. Une solution classique possède suffisamment de dérivées pour satisfaire les équations en chaque point ; lisse signifie ici différentiable à tout ordre. Une solution faible satisfait une formulation intégrée dans laquelle les dérivées sont transférées sur des fonctions test lisses. Ces notions d’existence répondent à des exigences différentes.

9.1 Solutions faibles (Leray–Hopf)

Pour des données initiales à divergence nulle et d’énergie finie dans le cadre usuel de l’espace entier ou périodique, avec un forçage admissible, la théorie de Leray–Hopf fournit des solutions faibles globales satisfaisant une inégalité d’énergie. Une solution lisse est unique au sein de la classe d’énergie correspondante tant qu’elle reste lisse : c’est l’unicité faible–forte. Une éventuelle perte de régularité ne signifie donc pas que toutes les solutions généralisées cessent d’exister. Présentation officielle du problème par Fefferman.

L’unicité dépend de la classe de solutions, des données initiales et de la force. Albritton, Brué et Colombo établissent une non-unicité pour certaines forces. Hou, Wang et Yang annoncent une non-unicité assistée par ordinateur sans force extérieure pour des données initiales singulières. Ces résultats ne doivent pas être présentés comme une explosion à partir de données initiales lisses sans force extérieure.

9.2 Critères de régularité (3D)

  • Critère de prolongement par la vorticité, de type Beale–Kato–Majda : pour une solution classique dans le cadre usuel suffisamment régulier de l’espace entier ou périodique, avec un forçage suffisamment régulier jusqu’à l’instant inclus, une perte de régularité à l’instant exige

La présentation de Luo des critères de prolongement pour Navier–Stokes énonce le critère sans force extérieure ; la condition sur le forçage ci-dessus exclut une perte de régularité imposée par une entrée extérieure singulière.

  • Prodi–Serrin : une solution faible dans la classe d’énergie, avec un forçage approprié, est régulière si avec

la solution est alors lisse sur .

  • Caffarelli–Kohn–Nirenberg : les solutions faibles dites « convenables » (suitable), sous les hypothèses du théorème sur le forçage, ont un ensemble singulier de mesure de Hausdorff parabolique unidimensionnelle nulle. C’est un résultat de régularité partielle ; il ne prouve pas que l’ensemble singulier est vide. Contexte et références dans la présentation officielle du problème.

Ce sont des énoncés conditionnels. Ils précisent quel contrôle empêcherait une perte de régularité, sans fournir ce contrôle pour chaque écoulement tridimensionnel. Les exposants ci-dessus reflètent le changement d’échelle naturel ; en trois dimensions, la norme d’énergie seule n’a pas l’échelle critique.

9.3 Régularité globale en 2D

Pour des données lisses compatibles dans le cadre usuel de l’espace entier ou périodique, avec un forçage suffisamment régulier, l’écoulement incompressible bidimensionnel reste lisse à tout instant fini. Sans force extérieure, l’enstrophie est décroissante au sens large ; le forçage ajoute un terme source. L’absence d’étirement des tourbillons permet un contrôle bien plus fort qu’en 3D.

9.4 Écoulements compressibles

Des solutions faibles globales d’énergie finie existent pour les gaz à loi en , sous certaines conditions, mais la régularité et l’unicité restent largement ouvertes.

9.5 Ce que demandent réellement les quatre alternatives de Clay

La formulation officielle accepte une démonstration de l’un quelconque de quatre énoncés, avec une viscosité strictement positive :

AlternativeDomaineForce et question
AForce nulle ; solutions lisses globales d’énergie finie pour toutes les données initiales lisses autorisées
BDomaine périodique Force nulle ; solutions lisses globales pour toutes les données initiales périodiques lisses autorisées
CCertaines données lisses autorisées et une force lisse à décroissance rapide n’admettent aucune solution lisse globale d’énergie uniformément bornée
DDomaine périodique Certaines données périodiques lisses autorisées et une force périodique lisse n’admettent aucune solution lisse globale

Dans l’espace entier, les données initiales et le forçage satisfont les conditions de décroissance spécifiées ; la force périodique obéit aussi à des conditions de décroissance temporelle. C et D autorisent une force non nulle ; démontrer l’une de ces alternatives ne réfuterait donc pas logiquement les énoncés sans force A ou B. « Résoudre une alternative de Clay » et « régler toutes les versions du problème de régularité » sont deux affirmations différentes.


10. Turbulence et spectres

Dans une turbulence homogène et isotrope à grand , la théorie K41 donne le spectre d’énergie de la zone inertielle

où est le taux moyen de dissipation et la constante de Kolmogorov.

Micro-échelle de Kolmogorov :

L’estimation K41 donne . Résoudre cette échelle en trois dimensions suggère un nombre de points de grille de l’ordre de , avant même le coût de l’intégration temporelle. Cette estimation heuristique motive différentes approches de résolution et de modélisation : la simulation numérique directe (DNS), la simulation des grandes échelles (LES) et les équations de Navier–Stokes moyennées au sens de Reynolds (RANS). LES et RANS nécessitent des hypothèses de fermeture pour les corrélations non résolues.

Le spectre en décrit un modèle statistique dans une zone inertielle, pas tous les écoulements ni toutes les échelles. L’intermittence et les frontières compliquent cette image. Une solution singulière spécialement construite ne fournirait pas un modèle universel de turbulence et ne supprimerait pas le besoin d’estimer les erreurs numériques.


11. Conditions aux limites et géométrie

  • Paroi avec adhérence :
  • Glissement libre : ,
  • Boîte périodique (tore) pour la théorie et la DNS
  • Entrée et sortie avec des prescriptions appropriées de contrainte ou de pression

Ce choix influe sur le caractère bien posé du problème, le bilan d’énergie et la stabilité numérique.

Une boîte périodique identifie les faces opposées ; elle n’a pas de parois physiques. Dans l’espace entier, la décroissance à l’infini remplace une condition de paroi. Un résultat dans ces géométries ne se transpose pas automatiquement à la surface d’un avion ou à une conduite : les conditions aux limites font partie du problème mathématique.


12. Méthodes numériques (guide condensé)

  • Méthodes de projection : calculer , résoudre l’équation de Poisson pour , puis projeter pour obtenir .
  • Schémas stables en énergie : une forme antisymétrique du terme convectif préserve l’énergie discrète.
  • Contrainte CFL pour les schémas explicites d’advection et de diffusion : .
  • Contrôle de la divergence : bases spectrales à divergence nulle, ou éléments finis mixtes avec couples vitesse–pression satisfaisant une condition de stabilité inf–sup. Une divergence exactement nulle nécessite des espaces adaptés ; la stabilité inf–sup seule ne la garantit pas.

Une simulation peut révéler un pic croissant, mais une grille finie ne suffit pas à établir une limite infinie. Des études de raffinement, un contrôle des erreurs de discrétisation et un argument analytique sont nécessaires pour distinguer une concentration réelle d’une instabilité numérique. Les preuves assistées par ordinateur fournissent des bornes rigoureuses ; les assistants de preuve formelle vérifient des déductions à partir de définitions et d’hypothèses encodées. Ni l’une ni l’autre démarche ne se réduit à une animation à haute résolution.


13. Liens et prolongements

  • Équations d’Euler : poser supprime la diffusion visqueuse. Le passage à la limite peut lui-même être délicat, notamment près des parois. Le résultat sans force extérieure annoncé séparément par OpenAI pour Euler concerne cette équation non visqueuse, pas le problème sans force à viscosité strictement positive. Annonce et article distinct sur Euler.
  • Hydrodynamique relativiste : une transposition naïve de Navier–Stokes est acausale ; les théories du second ordre d’Israel–Stewart peuvent rétablir la causalité sous des conditions appropriées sur les coefficients de transport et les états du fluide (un exemple rigoureux avec viscosité volumique).
  • Fluides quantiques : analogues superfluides de Navier–Stokes, avec modèles à deux fluides et tourbillons quantifiés.

14. État des connaissances

Au 13 octobre 2026 : la construction tridimensionnelle avec force extérieure lisse a été annoncée et fait l’objet d’un examen ; la régularité des équations de Navier–Stokes en 3D sans force extérieure reste ouverte. La régularité en deux dimensions et l’utilité quotidienne des solutions d’écoulements lisses ne sont pas remises en cause.

14.1 L’énoncé précis et son mécanisme

Le théorème 1.1 de l’article d’OpenAI énonce que, pour tout , il existe une force lisse à support compact dans l’espace et dans le temps et un écoulement incompressible lisse sur , initialement au repos, tels que

La vitesse et la pression ont un support spatial compact fixe avant l’instant 1 ; la force se prolonge de manière lisse au-delà de cet instant. L’article revendique l’alternative C et, grâce à sa construction périodique du corollaire 10.6, l’alternative D. Sa convention de pression absorbe la masse volumique constante.

Le mécanisme proposé concentre un tourbillon qui s’enroule vers l’intérieur tout en s’étirant suivant son axe. Des corrections oscillantes et leur transport de quantité de mouvement annulent des forces résiduelles qui seraient autrement singulières. Choisir un champ de vitesse qui explose et définir une force à partir de l’équation ne suffit pas : la force et toutes ses dérivées doivent rester lisses au passage de l’instant singulier. Il s’agit ici d’une esquisse conceptuelle de la construction revendiquée, pas d’une démonstration indépendante.

« Pour toute viscosité, il existe une force » autorise la force à dépendre de la viscosité. Cela ne signifie pas qu’une seule force fixe l’emporte sur toute viscosité, qu’un forçage arbitrairement faible suffise, ou qu’un fluide initialement au repos explose sans apport d’énergie.

14.2 Lean, examen mathématique et prix

Le dépôt Lean publié contient les énoncés cibles dans l’espace entier et dans le cas périodique, ainsi que des instructions de vérification indépendante avec Comparator. Un vérificateur formel valide les déductions pour l’énoncé exact et les hypothèses encodées. Savoir si ces définitions expriment le problème visé, et comment le raisonnement éclaire la preuve écrite, constitue une autre question. Ces notes ne prétendent pas avoir recompilé ni vérifié indépendamment le projet Lean.

Le 6 octobre, Bastounis, Circelli et Hansen ont publié une critique de la correspondance entre l’exposé en langage naturel et Lean : leurs exemples comprennent des exigences différentes sur les dérivées dans le lemme 8.6 et un argument différent pour l’estimation de pression du lemme 10.5. Cela remet en question l’idée que la formalisation certifierait automatiquement chaque étape de la preuve écrite. Cela ne réfute pas, à lui seul, le théorème final d’explosion avec forçage et ne montre pas que l’énoncé formel final est erroné.

Le 11 septembre, la réponse de Clay a salué cette avancée apparente tout en soulignant le rythme délibérément posé de son processus d’évaluation. Le règlement du prix exige une publication admissible, un délai d’au moins deux ans après cette publication et une acceptation générale par la communauté mathématique avant tout examen en vue du prix. La date de l’annonce ne suffit pas à établir le début de ce délai. Aucune attribution du prix n’est affirmée ici ; OpenAI indique ne pas avoir l’intention de le réclamer.

14.3 Ce que les travaux ultérieurs établissent, ou non

Les travaux de septembre qui ont suivi doivent eux aussi être interprétés selon leur portée exacte. Petrillo et Glimm distinguent explicitement le résultat avec forçage du problème sans force extérieure encore ouvert. L’exposé de Lei et Ren étudie la construction du profil, en reportant l’annulation oscillante ultérieure ; ce n’est pas une vérification indépendante complète. Constantin, Ignatova et Vicol obtiennent un résultat de régularité qui exclut une singularité sous des hypothèses précises de géométrie, de croissance et de forçage. Ce résultat concerne les tentatives de supprimer la force de ce mécanisme ; ce n’est pas un théorème de régularité globale pour tous les écoulements soumis à un forçage analytique.

Les questions fécondes sont désormais précises : les hypothèses et les conclusions revendiquées sont-elles vérifiées, le mécanisme peut-il être expliqué de manière transparente, et que subsiste-t-il lorsqu’on supprime le forçage ? La prévision de la turbulence, l’unicité des solutions faibles dans toutes les classes et les écoulements compressibles restent des problèmes distincts.

15. Vérifiez votre compréhension

  1. Translation uniforme : pour une vitesse constante, une pression constante et une force nulle, quels termes de l’équation incompressible s’annulent ? Toutes les dérivées spatiales et temporelles s’annulent ; ce mouvement est donc compatible avec toute viscosité strictement positive.
  2. Énergie et vitesse maximale : dans la section 6.1, effectuons le changement de variables . Pourquoi l’énergie varie-t-elle comme ? Le carré de l’amplitude apporte un facteur et le volume tridimensionnel un facteur .
  3. Changement d’échelle : à et fixés, que devient le nombre de Reynolds si est divisé par deux ? Il est divisé par deux ; l’importance relative de la viscosité augmente à la plus petite échelle.
  4. Lire un théorème : pourquoi l’alternative C peut-elle être vraie alors que A reste indécise ? C autorise le choix d’une force lisse non nulle, tandis que A concerne uniquement une force nulle.
  5. Vérificateur de preuve : que faut-il encore vérifier après que Lean a accepté une preuve ? Que les définitions encodées, les hypothèses et la conclusion finale correspondent à la question mathématique visée ; la correspondance avec un exposé distinct doit elle aussi être examinée.