Pourquoi Navier-Stokes 2D est plus facile qu'en 3D

Les équations de Navier-Stokes décrivent le mouvement des fluides. Elles fonctionnent en 2D (plat, comme de l'eau qui s'étale sur une table) et en 3D (la vraie vie, comme les courants océaniques ou l'air autour d'une aile). Les équations sont presque identiques dans les deux cas.

Mais voici la surprise : en 2D, les mathématiciens peuvent prouver que les équations se comportent toujours bien. Les mathématiques ne se brisent jamais. En 3D ? Personne ne sait. Le fluide pourrait faire quelque chose de si violent et soudain que les mathématiques cesseraient complètement de fonctionner. Prouver si cela peut arriver est le problème du Millénaire de Clay, doté d'un million de dollars.

Ce n'est pas simplement « la 3D est plus difficile parce qu'il y a plus de choses ». Il existe un mécanisme spécifique en 3D qui n'existe pas en 2D, et il change tout.

Les équations de Navier-Stokes incompressibles sur $\mathbb{R}^n$ (ou un domaine périodique $\mathbb{T}^n$) avec $\nu > 0$ sont

$$\partial_t u + (u \cdot \nabla)u = -\nabla p + \nu \Delta u, \qquad \nabla \cdot u = 0.$$

Pour $n = 2$, l'existence globale et l'unicité des solutions classiques pour des données initiales lisses à divergence nulle est un théorème. Les références clés sont Ladyzhenskaya (1959), s'appuyant sur les travaux antérieurs de Leray (1934). Le résultat s'étend aux solutions lisses sur $\mathbb{R}^2$, $\mathbb{T}^2$ et les domaines bornés avec conditions aux limites standard.

Pour $n = 3$, l'existence globale de solutions classiques à partir de données lisses arbitraires est un problème ouvert. C'est le contenu du problème du Millénaire de Clay tel que formulé par Fefferman (2000). Leray (1934) a établi l'existence globale de solutions faibles dans $L^2(\mathbb{R}^3)$, mais l'unicité et la régularité de ces solutions restent non résolues.

L'écart entre $n = 2$ et $n = 3$ n'est pas une question de comptabilité. Il reflète une différence fondamentale dans la structure de l'équation de vorticité, les propriétés de changement d'échelle de la non-linéarité, et les estimations a priori disponibles. Chacun de ces points sera examiné dans les sections qui suivent.

La question à un million de dollars ne concerne que la 3D. Pourquoi ? Parce que la version 2D est déjà résolue. Les mathématiciens ont prouvé il y a des décennies que les solutions de Navier-Stokes 2D restent toujours lisses. Pas besoin de prix — c'est fait.

La vraie question n'est donc pas seulement « pourquoi la 3D est-elle difficile ? » C'est « pourquoi la 2D est-elle facile et la 3D difficile ? » Qu'est-ce qui se brise quand on ajoute cette troisième dimension ?

Le problème de Clay (Fefferman 2000) considère le problème de Cauchy pour Navier-Stokes incompressible sur $\mathbb{R}^3$ avec viscosité $\nu > 0$ et données initiales lisses à divergence nulle $u_0 \in C^\infty(\mathbb{R}^3)$ satisfaisant des conditions de décroissance appropriées. La question : existe-t-il une solution lisse $u \in C^\infty(\mathbb{R}^3 \times [0,\infty))$ à énergie bornée ?

L'énoncé analogue pour $\mathbb{R}^2$ est un théorème. Ladyzhenskaya a prouvé l'existence et l'unicité globales de solutions fortes pour Navier-Stokes 2D avec $u_0 \in H^1(\mathbb{R}^2)$, et le bootstrap donne la régularité $C^\infty$ pour des données lisses. La preuve repose sur des estimations a priori spécifiques à deux dimensions qui ne s'étendent pas à trois.

Le problème du Millénaire est donc purement tridimensionnel. Les résultats partiels qui existent en 3D — y compris les solutions faibles de Leray (1934), le théorème de régularité partielle de Caffarelli-Kohn-Nirenberg (1982), et divers critères de régularité conditionnelle — restent tous en deçà de la résolution complète de la question de régularité globale. Chaque résultat partiel met en lumière un écart spécifique dans notre contrôle des solutions 3D.

L'arme secrète en 2D est la vorticité — la mesure de la rotation du fluide en chaque point. Imaginez de petits tourbillons dispersés dans le fluide.

En 2D, la vorticité est simplement un nombre en chaque point : rotation horaire ou antihoraire, rapide ou lente. Et voici la magie — en 2D, ces petits tourbillons peuvent dériver et s'atténuer, mais ils ne peuvent jamais devenir plus forts. La rotation maximale du début est la rotation maximale pour toujours. C'est comme une règle que le fluide doit respecter : « pas de nouvelle rotation autorisée ».

C'est incroyablement puissant. Une fois qu'on sait que la rotation ne peut pas s'emballer, on peut prouver que tout le reste reste sous contrôle. La vitesse reste lisse, la pression reste lisse, et la solution fonctionne pour toujours. C'est comme une chaîne de dominos — renversez le premier (la vorticité est bornée) et tout le reste suit.

En 2D, la vorticité $\omega = \partial_1 u_2 - \partial_2 u_1$ est un scalaire satisfaisant

$$\partial_t \omega + (u \cdot \nabla)\omega = \nu \Delta \omega.$$

C'est une équation d'advection-diffusion pour la vorticité scalaire $\omega$. Bien que le système complet reste non linéaire parce que $u$ est reconstruit à partir de $\omega$ via la loi de Biot-Savart $u = \nabla^\perp (-\Delta)^{-1} \omega$, l'équation ne contient pas de terme source d'étirement tourbillonnaire : le membre de droite ne contient que le laplacien, et non le terme $(\omega \cdot \nabla)u$ qui apparaît en 3D.

Le principe du maximum scalaire s'applique directement : $\|\omega(t)\|_{L^\infty} \leq \|\omega_0\|_{L^\infty}$ pour tout $t \geq 0$. Simultanément, les normes $L^p$ de $\omega$ sont décroissantes pour tout $1 \leq p \leq \infty$.

À partir de la borne $L^\infty$ sur $\omega$, on récupère $u \in L^\infty(0,T; W^{1,p})$ pour tout $p < \infty$ via les estimations de Calderón-Zygmund sur le noyau de Biot-Savart. La régularité supérieure suit en différenciant l'équation de vorticité et en appliquant le bootstrap parabolique : chaque dérivée spatiale de $\omega$ satisfait une équation parabolique à coefficients contrôlables, de sorte que les bornes se propagent à tous les ordres.

Le principe du maximum de vorticité 2D a des antécédents plus anciens dans le cas non visqueux. Wolibner (1933) a prouvé l'existence globale pour l'Euler 2D dans les espaces de Hölder, et Yudovich (1963) a établi l'unicité pour les solutions d'Euler 2D à vorticité bornée. Avec viscosité ($\nu > 0$), le lissage parabolique ne fait que renforcer ces estimations. La preuve de Ladyzhenskaya de la régularité globale de Navier-Stokes 2D repose sur cette structure, combinée à l'inégalité d'interpolation de Ladyzhenskaya $\|f\|_{L^4}^2 \leq C \|f\|_{L^2} \|\nabla f\|_{L^2}$ (valide en 2D, avec une structure d'exposants différente de son homologue 3D).

En 3D, la vorticité n'est pas juste un nombre — c'est une direction et une intensité. Imaginez de petits tubes de tornades traversant le fluide, chacun pointant dans une direction et tournant à une certaine vitesse.

Et voici ce qui brise tout : en 3D, ces tubes de tornades peuvent être étirés. Imaginez qu'on tire un tube en rotation comme un caramel — il s'affine et tourne plus vite. C'est l'étirement tourbillonnaire, et c'est le méchant de l'histoire.

La belle règle de la 2D (« pas de nouvelle rotation autorisée ») est complètement détruite. En 3D, le fluide peut amplifier sa propre rotation. L'étirement injecte de l'énergie dans des échelles de plus en plus petites, créant potentiellement une rotation infiniment intense en un seul point. Ce serait une « explosion » — l'effondrement des mathématiques.

La viscosité (la friction interne du fluide) arrête-t-elle toujours l'étirement avant qu'il ne devienne infini ? Ou l'étirement peut-il parfois gagner ? Personne ne sait. C'est littéralement la question à un million de dollars.

Cette compétition entre étirement et friction est au cœur de pourquoi le problème est si difficile.

En 3D, la vorticité $\omega = \nabla \times u$ satisfait

$$\partial_t \omega + (u \cdot \nabla)\omega = (\omega \cdot \nabla)u + \nu \Delta \omega.$$

Le terme $(\omega \cdot \nabla)u$ est le terme d'étirement tourbillonnaire, absent en 2D. Il est quadratique au sens où $u$ est reconstruit à partir de $\omega$ via la loi de Biot-Savart 3D, de sorte que $(\omega \cdot \nabla)u$ est approximativement de l'ordre de $|\omega|^2$ dans son pire cas de changement d'échelle. Ce terme permet la croissance de $\|\omega\|_{L^\infty}$ et détruit le principe du maximum scalaire disponible en 2D.

Pour l'Euler 3D, le critère de Beale-Kato-Majda (1984) énonce qu'une solution lisse sur $[0, T)$ se prolonge au-delà du temps $T$ si et seulement si

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

Des critères de continuation analogues s'appliquent à Navier-Stokes 3D. L'explosion, si elle se produit, nécessite que $\|\omega\|_{L^\infty}$ devienne non intégrable en temps. Le terme d'étirement est le terme source dans l'équation de vorticité qui peut amplifier la vorticité et empêcher un argument de type principe du maximum. En 2D, $\|\omega\|_{L^\infty}$ est borné par les données initiales pour tout temps ; en 3D, contrôler cette norme est le problème ouvert central.

Progrès partiel : Caffarelli-Kohn-Nirenberg (1982) ont montré que pour toute solution faible appropriée de Navier-Stokes 3D, l'ensemble des points singuliers dans l'espace-temps a une mesure de Hausdorff parabolique unidimensionnelle nulle. Cela signifie que les singularités, si elles existent, sont extrêmement rares. Mais le théorème n'exclut pas leur existence.

Pour le cas non visqueux, Elgindi (2021) a prouvé la formation de singularités en temps fini pour l'Euler 3D avec des données initiales $C^{1,\alpha}$ ($\alpha$ petit), utilisant un mécanisme entraîné par l'étirement tourbillonnaire le long d'un axe de symétrie. Cela n'implique pas directement l'explosion pour Navier-Stokes (la viscosité pourrait encore régulariser), mais montre que le mécanisme d'étirement est assez puissant pour produire des singularités sans amortissement visqueux.

Il y a une raison plus profonde pour laquelle la 3D est difficile, au-delà de l'étirement tourbillonnaire. Cela a à voir avec ce qui se passe quand on « zoome » sur le fluide.

Les équations de Navier-Stokes ont une astuce de zoom : si on a une solution, on peut zoomer dans une région plus petite et accélérer le temps, et on obtient une autre solution valide. La question est : quand on zoome, l'énergie augmente-t-elle, diminue-t-elle, ou reste-t-elle la même ?

• En 2D, zoomer garde l'énergie identique. C'est ce qu'on appelle un changement d'échelle critique. Cela signifie que les estimations d'énergie fonctionnent à toutes les échelles — grandes ou petites. On ne perd jamais le contrôle.
• En 3D, zoomer fait croître l'énergie. C'est le changement d'échelle supercritique. Cela signifie qu'aux petites échelles, les effets non linéaires violents sont relativement plus forts que les effets visqueux calmants. Les outils mathématiques perdent leur prise exactement là où on en a le plus besoin.

Pensez-y ainsi : en 2D, votre lampe de poche est toujours assez brillante pour voir ce que fait le fluide. En 3D, plus on regarde petit, plus la lampe faiblit, tandis que le fluide devient plus chaotique. Finalement, on est dans le noir.

Ce n'est pas un simple inconvénient technique — c'est une barrière fondamentale. Les outils mathématiques standard ne sont tout simplement pas assez puissants pour contrôler Navier-Stokes 3D aux petites échelles. De nouvelles idées sont nécessaires.

Les équations de Navier-Stokes sont invariantes sous le changement d'échelle

$$u(x,t) \mapsto \lambda u(\lambda x, \lambda^2 t), \qquad p(x,t) \mapsto \lambda^2 p(\lambda x, \lambda^2 t),$$

pour tout $\lambda > 0$. Sous ce changement d'échelle, la norme $L^2$ se transforme comme $\|u_\lambda\|_{L^2(\mathbb{R}^n)} = \lambda^{1-n/2} \|u\|_{L^2}$.

• Pour $n = 2$ : $\|u_\lambda\|_{L^2} = \|u\|_{L^2}$. L'énergie est une quantité invariante d'échelle (critique). L'équation est énergie-critique.
• Pour $n = 3$ : $\|u_\lambda\|_{L^2} = \lambda^{-1/2} \|u\|_{L^2}$, qui croît quand $\lambda \to \infty$ (zoom avant). La norme d'énergie est supercritique : elle devient plus faible par rapport au changement d'échelle aux petites échelles.

L'espace de Sobolev critique pour Navier-Stokes 3D est $\dot{H}^{1/2}(\mathbb{R}^3)$, qui est invariant d'échelle sous le changement d'échelle naturel. Mais l'identité d'énergie ne contrôle $u$ que dans $L^\infty_t L^2_x \cap L^2_t \dot{H}^1_x$, soit une demi-dérivée en dessous de l'espace critique. C'est l'écart de supercriticalité.

En 2D, l'identité d'énergie $\frac{d}{dt}\|u\|_{L^2}^2 = -2\nu\|\nabla u\|_{L^2}^2$ fournit exactement le contrôle au niveau critique. Combinée au principe du maximum de vorticité, elle donne suffisamment de régularité pour bootstrapper jusqu'à $C^\infty$. En 3D, la même identité d'énergie fournit une borne plus faible que le niveau critique, et aucune estimation a priori supplémentaire n'est connue pour combler l'écart.

Tao a souligné la barrière de supercriticalité : des arguments fondés uniquement sur l'inégalité d'énergie et le changement d'échelle ne sont pas censés résoudre la régularité globale en 3D, de sorte que toute preuve réussie nécessitera probablement une structure supplémentaire. En particulier, les méthodes « aveugles à la supercriticalité » (traitant l'équation uniquement par son changement d'échelle et sa structure d'énergie) ne peuvent pas réussir. La structure algébrique spécifique de l'équation — telle que la condition de divergence nulle ou la structure antisymétrique de la non-linéarité — devrait jouer un rôle. Pour un traitement plus approfondi de ces barrières, voir Pourquoi Navier-Stokes est difficile.

Si la preuve 2D fonctionne grâce à une vorticité bornée et un changement d'échelle critique, et que la 3D n'a ni l'un ni l'autre… de quoi une preuve 3D aurait-elle besoin ?

Honnêtement ? Personne ne sait. Mais voici les principales idées sur lesquelles les gens travaillent :

• Trouver un nouveau « bouton de contrôle ». En 2D, la vorticité est le bouton de contrôle — elle reste bornée et tout le reste suit. En 3D, il faut autre chose qui reste sous contrôle et soit assez puissant pour maintenir la solution lisse. Personne ne l'a encore trouvé.
• Exploiter une structure cachée. Le fluide doit satisfaire une contrainte : il est incompressible (il ne peut pas être comprimé dans un volume plus petit). Cette contrainte limite ce que l'étirement tourbillonnaire peut faire. Peut-être y a-t-il un motif géométrique plus profond ici que personne n'a encore pleinement exploité.
• Prouver que ça casse vraiment. Il y a une autre possibilité : peut-être que les solutions 3D peuvent exploser. Le prouver serait tout aussi important. Il faudrait construire un exemple spécifique où l'étirement tourbillonnaire l'emporte sur la viscosité. Pour les équations d'Euler plus simples (Navier-Stokes sans friction), la formation de singularités a été démontrée dans des cadres apparentés, mais le cas visqueux reste ouvert.

Pour en savoir plus sur ce qui a été tenté, voir Sous-problèmes de Navier-Stokes.

Une preuve de régularité globale en 3D nécessiterait de combler l'écart de supercriticalité. Concrètement, il faut une estimation a priori de la forme $\|u(t)\|_X \leq C(\|u_0\|_Y, t)$ pour une norme $X$ au niveau critique ou au-dessus, où $C$ reste fini pour tout $t$. La borne d'énergie connue $u \in L^\infty_t L^2_x \cap L^2_t \dot{H}^1_x$ est une demi-dérivée en dessous du niveau critique et ne suffit pas.

Plusieurs programmes de recherche ciblent cet écart :

• Décomposition en profils et concentration-compacité. Adaptées du succès des équations dispersives critiques (Kenig-Merle 2006), ces méthodes cherchent à classifier les profils d'explosion. Pour Navier-Stokes, des résultats partiels existent (par ex. Gallagher-Koch-Planchon 2016), mais la nature supercritique de l'énergie rend le programme complet plus difficile à exécuter que dans les cas de l'onde énergie-critique ou de Schrödinger.
• Extensions de solutions mild. Le cadre de Fujita-Kato (1964) donne le caractère bien posé local dans $\dot{H}^{1/2}(\mathbb{R}^3)$ et le caractère bien posé global pour de petites données dans les espaces critiques ($L^3$, $\dot{H}^{1/2}$, $BMO^{-1}$). La question est de savoir si les solutions à grandes données peuvent être continuées globalement, ce qui nécessite de contrôler la norme critique.
• Critères de régularité. Au-delà de Beale-Kato-Majda ($\int_0^T \|\omega\|_{\infty} < \infty$), il y a les conditions de Prodi-Serrin ($u \in L^p_t L^q_x$ avec $2/p + 3/q = 1$, $q > 3$), Escauriaza-Seregin-Šverák ($u \in L^\infty_t L^3_x$, 2003), et d'autres critères aux points limites. Chacun réduit la régularité globale à une seule estimation a priori, mais prouver cette estimation reste ouvert.
• Construction d'explosion. Tao (2016) a construit une solution explosive pour un système de Navier-Stokes moyenné qui respecte l'identité d'énergie et le changement d'échelle mais pas la structure complète de divergence nulle. Cela indique que toute preuve de régularité doit utiliser la structure géométrique spécifique de la non-linéarité, pas seulement ses propriétés de changement d'échelle. La question de savoir si le vrai Navier-Stokes admet l'explosion est ouverte.

Pour le problème non visqueux, l'explosion $C^{1,\alpha}$ d'Elgindi pour l'Euler 3D (2021) montre que l'étirement tourbillonnaire peut produire des singularités en dessous de la régularité $C^\infty$. La question de l'explosion lisse ($C^\infty$) pour Euler reste ouverte, tout comme la question de savoir si la viscosité peut arrêter de tels mécanismes dans le cadre de Navier-Stokes.

Voici la différence fondamentale en un coup d'œil :

L'écart entre 2D et 3D n'est pas un petit détail technique. C'est un mur. La stratégie de preuve qui fonctionne parfaitement en 2D n'a pas juste « besoin d'un peu plus de travail » en 3D — elle ne peut fondamentalement pas fonctionner parce que la structure sur laquelle elle repose n'existe pas.

Pour les équations complètes, voir Que sont les équations de Navier-Stokes ?. Pour le problème ouvert précis, voir Existence et régularité de Navier-Stokes. Pour comprendre pourquoi c'est si difficile, voir Pourquoi Navier-Stokes est difficile.

Les contrastes suivants résument le fossé mathématique :

Le résultat 2D n'est pas un simple échauffement en dimension inférieure. C'est un théorème complet dont le mécanisme de preuve — le principe du maximum de vorticité combiné à la criticalité de l'énergie — n'a pas d'homologue connu en 3D. Toute résolution du problème 3D, qu'il s'agisse de régularité ou d'explosion, nécessitera des idées fondamentalement nouvelles. Pour l'état actuel des résultats partiels et des programmes de recherche, voir Sous-problèmes de Navier-Stokes et Pourquoi Navier-Stokes est difficile.