Euler contre Navier-Stokes : quelle est la différence ?

Les équations d'Euler décrivent le mouvement d'un fluide sans friction interne — sans viscosité. Les équations de Navier-Stokes décrivent le même fluide avec viscosité.

Mathématiquement, la différence entre les équations d'Euler et de Navier-Stokes se résume à un seul terme : $\nu \Delta u$, le terme de diffusion visqueuse. Supprimez-le et Navier-Stokes devient Euler. Conservez-le et l'équation acquiert un mécanisme de lissage qui modifie fondamentalement à la fois la physique et l'analyse.

Ce seul terme explique pourquoi la fumée se dissipe, pourquoi des couches limites se forment le long des surfaces et pourquoi le problème du Millénaire de Navier-Stokes a un caractère différent de la question correspondante pour Euler.

Les équations d'Euler incompressibles sur $\mathbb{R}^3$ s'écrivent

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

Les équations de Navier-Stokes incompressibles s'écrivent

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

avec $\nu > 0$ la viscosité cinématique. Formellement, poser $\nu = 0$ dans Navier-Stokes redonne Euler. Mais cette substitution formelle dissimule le fait que la limite $\nu \to 0$ est singulière : le terme visqueux $\nu \Delta u$ porte la dérivée spatiale d'ordre le plus élevé du système, et le supprimer change le type de l'EDP et les espaces fonctionnels dans lesquels vivent les solutions.

Voici les deux équations sous leur forme incompressible standard, écrites pour rendre la comparaison évidente :

Euler :

$$\partial_t u + (u \cdot \nabla)u = -\nabla p$$

Navier-Stokes :

$$\partial_t u + (u \cdot \nabla)u = -\nabla p + \nu \Delta u$$

Les deux équations partagent le même membre de gauche : le taux de variation de la vitesse plus le terme non linéaire d'auto-transport $(u \cdot \nabla)u$. Toutes deux imposent l'incompressibilité par $\nabla \cdot u = 0$. La seule différence structurelle est le terme visqueux $\nu \Delta u$ dans le membre de droite de Navier-Stokes.

Le paramètre $\nu$ est la viscosité cinématique — une constante physique du fluide. Le miel a un $\nu$ élevé. L'air a un $\nu$ faible. Les équations d'Euler correspondent à $\nu = 0$ : une idéalisation parfaitement sans friction qui peut approcher certains écoulements à grand nombre de Reynolds loin des parois, mais qui n'est jamais exactement réalisée dans les fluides ordinaires.

Les deux systèmes partagent la forme bilinéaire $(u \cdot \nabla)u$ et la pression déterminée implicitement par la contrainte de divergence nulle. En prenant la divergence de l'équation de quantité de mouvement et en utilisant $\nabla \cdot u = 0$, on obtient l'équation de Poisson pour la pression

$$-\Delta p = \nabla \cdot ((u \cdot \nabla)u) = \partial_i \partial_j (u_i u_j),$$

qui est identique pour les deux systèmes. La pression est une fonctionnelle non locale de $u$ dans les deux cas.

Le terme visqueux $\nu \Delta u$ est un opérateur elliptique linéaire du second ordre agissant sur chaque composante de la vitesse. Il s'agit d'une régularisation parabolique : sa présence fait de Navier-Stokes un système parabolique semi-linéaire, tandis qu'Euler incompressible est une équation de transport non linéaire du premier ordre avec un couplage non local par la pression. Cette différence de type d'EDP est à l'origine de presque toutes les différences ultérieures en théorie de la régularité.

La viscosité est la friction interne entre couches voisines de fluide. Lorsqu'une couche rapide se trouve à côté d'une couche lente, la viscosité transfère la quantité de mouvement entre elles, lissant la différence de vitesse.

Cela entraîne trois conséquences physiques majeures qui séparent Navier-Stokes d'Euler :

• Dissipation. L'énergie cinétique se convertit en chaleur. Remuez une tasse de café et arrêtez — le café finit par s'immobiliser. Les équations d'Euler ne peuvent pas prédire cela. Dans le modèle idéalisé d'Euler, il n'existe aucun mécanisme visqueux pour dissiper le mouvement en chaleur.
• Couches limites. Les fluides réels adhèrent aux surfaces (la condition de non-glissement). Cela crée de fines couches de variation rapide de la vitesse près des parois — les couches limites — qui engendrent la traînée sur les avions, le frottement dans les tuyaux et la turbulence à grande vitesse. Les écoulements d'Euler satisfont une condition de glissement aux parois plutôt que de non-glissement, et manquent donc la traînée visqueuse et l'essentiel de la physique des couches limites observée dans les fluides réels.
• Lissage à petite échelle. La viscosité supprime les gradients de vitesse les plus abrupts. Sans elle, rien n'empêche l'écoulement de développer une structure infiniment fine. Ce lissage est précisément ce qui rend la question de régularité pour Navier-Stokes différente du cas d'Euler.

L'identité énergétique pour Navier-Stokes sur $\mathbb{R}^3$ (ou un domaine périodique) s'écrit

$$\frac{d}{dt}\frac{1}{2}\|u\|_{L^2}^2 = -\nu \|\nabla u\|_{L^2}^2,$$

de sorte que l'énergie cinétique est dissipée de façon monotone. Pour Euler ($\nu = 0$), la norme $L^2$ de $u$ est formellement conservée.

Au niveau des conditions aux limites, Navier-Stokes avec la condition de non-glissement ($u|_{\partial \Omega} = 0$) est un problème aux limites parabolique bien posé, tandis qu'Euler ne requiert que la condition d'imperméabilité ($u \cdot n = 0$). L'incompatibilité entre ces conditions aux limites lorsque $\nu \to 0$ engendre la couche limite de Prandtl, un phénomène de perturbation singulière étudié depuis Prandtl (1904).

Physiquement, la viscosité agit comme un filtre passe-bas : elle amortit les modes de Fourier au taux $\nu |k|^2$, supprimant préférentiellement les petites échelles. Cet amortissement spectral est le mécanisme sous-jacent à l'échelle de dissipation de Kolmogorov $\eta \sim (\nu^3 / \varepsilon)^{1/4}$ en turbulence. Voir Nombre de Reynolds et turbulence pour le tableau complet des lois d'échelle.

Formellement, oui. Posez $\nu = 0$ dans les équations de Navier-Stokes et vous obtenez les équations d'Euler. Mais cela est trompeur si l'on s'arrête là.

La limite $\nu \to 0$ est singulière, pas lisse. La viscosité est responsable des dérivées d'ordre le plus élevé de l'équation. La supprimer ne se réduit pas à un petit changement — cela modifie le caractère fondamental de l'EDP. Les couches limites ne s'amincissent pas progressivement ; elles peuvent devenir turbulentes. Des solutions qui étaient lisses sous Navier-Stokes peuvent présenter un comportement totalement différent sous Euler.

Une analogie : supprimer le frottement d'un système mécanique ne fait pas simplement glisser les choses plus facilement. Cela peut rendre le système qualitativement différent — des orbites stables deviennent chaotiques, des contraintes imposées par le frottement sont entièrement relâchées. Il en va de même ici.

Ainsi, bien qu'Euler et Navier-Stokes partagent le même ADN mathématique, la limite à viscosité nulle est l'un des problèmes les plus profonds et les plus subtils de la dynamique des fluides, et non une simple substitution.

La limite non visqueuse $\nu \to 0$ est une perturbation singulière : $\nu \Delta u$ porte les dérivées spatiales d'ordre le plus élevé du système, de sorte que poser $\nu = 0$ abaisse l'ordre de l'EDP. Sur les domaines avec bord, la limite est liée à la validité du développement en couche limite de Prandtl, qui peut échouer de manière spectaculaire (Grenier 2000, Gérard-Varet & Dormy 2010).

Sur $\mathbb{R}^3$ ou $\mathbb{T}^3$ (sans bord), la situation est plus claire mais reste non triviale. Si la solution d'Euler $u^E$ demeure lisse sur $[0,T]$, alors les solutions de Navier-Stokes $u^\nu$ convergent vers $u^E$ dans $L^2$ lorsque $\nu \to 0$ (Kato 1972). Mais savoir si les solutions d'Euler restent lisses globalement — si une explosion en temps fini se produit pour des données lisses — est lui-même un problème ouvert en 3D.

La limite interagit également avec la théorie de la turbulence. Le modèle de Kolmogorov nécessite $\nu > 0$ pour définir une échelle de dissipation, pourtant la dissipation anomale — la persistance de la dissipation d'énergie lorsque $\nu \to 0$ — est une conjecture centrale. La conjecture d'Onsager (désormais un théorème : Isett 2018, avec optimalité par Buckmaster–De Lellis–Székelyhidi–Vicol 2019) caractérise les conditions sous lesquelles les solutions d'Euler peuvent dissiper de l'énergie en l'absence de viscosité.

Les équations d'Euler sont utilisées lorsque la viscosité est négligeable par rapport aux autres forces de l'écoulement. Cela se produit plus souvent qu'on ne le pense :

• Aérodynamique à grande vitesse loin des surfaces. Loin de la surface d'une aile, l'écoulement d'air se comporte de manière presque non visqueuse. Les ingénieurs utilisent des solveurs d'Euler pour l'écoulement en masse et ajoutent des modèles de couche limite près des parois.
• Écoulements astrophysiques. Le gaz interstellaire, les intérieurs stellaires et les disques d'accrétion opèrent à des échelles si vastes que la viscosité moléculaire est négligeable (bien que la viscosité turbulente effective puisse ne pas l'être).
• Dynamique des gaz compressibles. Les ondes de choc, les détonations et les écoulements supersoniques sont souvent modélisés par les équations d'Euler compressibles, où la physique dominante est la pression et l'inertie, et non le frottement.
• Analyse théorique. Euler est parfois étudié comme une étape vers Navier-Stokes, ou comme un système d'EDP intéressant en soi — avec des connexions profondes à la géométrie, à la dynamique des tourbillons et à la structure de la turbulence.

Mais pour tout ce qui implique le frottement, la traînée ou le comportement aux parois — écoulement dans les tuyaux, aérodynamique des véhicules près des surfaces, circulation sanguine, météorologie à l'échelle humaine — Navier-Stokes est le bon modèle.

Les équations d'Euler régissent le comportement à l'ordre dominant à grand nombre de Reynolds $\mathrm{Re} = UL/\nu \gg 1$, où les effets visqueux sont confinés à de fines couches limites et couches de cisaillement internes. Dans les écoulements en masse à grand nombre de Reynolds, loin des parois, Euler donne souvent l'approximation à l'ordre dominant, tandis que les corrections visqueuses se concentrent dans de fines couches limites ou de cisaillement.

Les équations d'Euler compressibles — un système hyperbolique à vitesse de propagation finie — constituent le modèle standard de la dynamique des gaz, y compris la formation de chocs et les problèmes de Riemann. Elles diffèrent des équations d'Euler incompressibles discutées ci-dessus : Euler compressible est véritablement hyperbolique, tandis qu'Euler incompressible possède un couplage non local par la pression et une vitesse de propagation infinie.

En analyse mathématique, Euler sert à la fois d'objet limite pour le problème de viscosité évanescente et de système d'EDP riche possédant sa propre théorie de régularité, ses quantités conservées (hélicité, Casimirs via le cadre d'Euler-Arnold sur le groupe des difféomorphismes) et ses connexions à la mécanique géométrique.

La question de la régularité est l'endroit où l'écart entre Euler et Navier-Stokes compte le plus.

Le problème du Millénaire de Navier-Stokes demande : si l'on part d'un écoulement lisse et bien posé en trois dimensions, la solution de Navier-Stokes reste-t-elle lisse pour tout temps ? Personne ne le sait.

La même question pour Euler est également ouverte en 3D. Mais les deux problèmes ont des caractères fondamentalement différents :

• Navier-Stokes bénéficie de la viscosité — qui lisse toujours, dissipe toujours l'énergie, supprime toujours les gradients les plus abrupts. La question est de savoir si ce lissage est suffisamment fort pour empêcher le terme non linéaire de créer une singularité.
• Euler n'a aucun lissage. Le terme non linéaire peut amplifier les gradients sans rien pour l'arrêter. Savoir si Euler en 3D développe une singularité en temps fini à partir de données initiales lisses est un problème ouvert.

En deux dimensions, les deux équations sont globalement bien posées pour des données initiales lisses. L'histoire en 2D est réglée. C'est en 3D que réside le mystère, pour les deux équations, de manières différentes.

Le paysage de la régularité :

2D : L'existence globale et l'unicité de solutions lisses sont connues pour les deux systèmes. Pour Euler 2D avec des données lisses, Wolibner (1933) a prouvé l'existence globale dans les espaces de Hölder ; Yudovich (1963) a établi l'unicité pour des données à vorticité bornée. Pour Navier-Stokes 2D, la régularité globale découle de l'inégalité de Ladyzhenskaya et du principe du maximum pour la vorticité.

Navier-Stokes 3D : Leray (1934) a prouvé l'existence globale de solutions faibles dans $L^2$, mais l'unicité et la régularité restent ouvertes. Le théorème de Caffarelli–Kohn–Nirenberg (1982) montre que l'ensemble singulier a une mesure de Hausdorff parabolique unidimensionnelle nulle — toute explosion, si elle se produit, est extrêmement clairsemée. Le terme visqueux fournit l'estimation a priori cruciale $\int_0^T \|\nabla u\|_{L^2}^2 \, dt \leq C(u_0)$, mais ce contrôle $H^1$ est sous-critique pour le changement d'échelle 3D et insuffisant pour fermer un argument de bootstrap. Voir Pourquoi Navier-Stokes est difficile pour l'écart de surcriticalité.

Euler 3D : Il n'existe aucune théorie globale pour des données lisses. Le caractère bien posé local dans les espaces de Sobolev $H^s$, $s > 5/2$, est classique (Kato 1972, Kato–Ponce 1988). Le critère de Beale–Kato–Majda (1984) réduit la détection d'explosion à $\int_0^T \|\omega\|_{L^\infty} \, dt$ : la solution reste lisse sur $[0,T]$ si et seulement si $\int_0^T \|\omega\|_{L^\infty} \, dt < \infty$. L'explosion requiert que le supremum de la vorticité croisse suffisamment vite pour être non intégrable en temps. Elgindi (2021, Annals of Mathematics) a prouvé la formation de singularité en temps fini pour des données $C^{1,\alpha}$ — une percée, mais en deçà du seuil lisse ($C^\infty$). Savoir si les solutions lisses d'Euler explosent en 3D reste ouvert.

La turbulence est l'endroit où la comparaison Euler-Navier-Stokes devient la plus physiquement parlante.

Dans un écoulement turbulent, l'énergie entre à grande échelle (la taille du tuyau, de l'aile, de la tempête) et cascade vers des tourbillons de plus en plus petits. C'est la cascade d'énergie. Au bas de la cascade, la viscosité finit par convertir l'énergie cinétique en chaleur, arrêtant le processus.

Les équations d'Euler peuvent modéliser le haut et le milieu de la cascade — la dynamique à grande échelle et dans la zone inertielle. Mais elles ne peuvent pas modéliser le bas, là où la viscosité agit. Sans viscosité, il n'y a pas d'échelle la plus petite naturelle, et l'énergie n'a nulle part où aller.

C'est pourquoi la modélisation de la turbulence fait presque toujours appel à Navier-Stokes, et non à Euler. Le nombre de Reynolds $\mathrm{Re} = UL/\nu$ contrôle l'étendue de la cascade : un $\mathrm{Re}$ élevé signifie de nombreuses décades d'échelles entre l'injection d'énergie et la dissipation visqueuse. La turbulence réelle vit dans l'interaction entre la cascade non visqueuse et la coupure visqueuse.

La théorie de Kolmogorov de 1941 prédit un spectre d'énergie $E(k) \sim \varepsilon^{2/3} k^{-5/3}$ dans la zone inertielle, où ni le forçage à grande échelle ni la dissipation visqueuse ne dominent. L'échelle de dissipation $\eta = (\nu^3/\varepsilon)^{1/4}$ fixe le bas de cette zone.

Les équations d'Euler décrivent la limite $\nu = 0$, où la zone inertielle s'étend à toutes les échelles. La question de savoir si la dissipation d'énergie persiste dans cette limite — la dissipation anomale — est la conjecture d'Onsager. Le côté rigide (pas de dissipation pour $u \in C^{0,\alpha}$, $\alpha > 1/3$) a été prouvé par Constantin–E–Titi (1994). Le côté flexible (existence de solutions dissipatives d'Euler en dessous de $C^{1/3}$) a été complété par Isett (2018), s'appuyant sur le programme d'intégration convexe de De Lellis–Székelyhidi.

Pour Navier-Stokes, le tableau de la cascade est inscrit dans le bilan d'énergie : $\varepsilon = \nu \langle |\nabla u|^2 \rangle$. La question ouverte clé est de savoir si les solutions de Navier-Stokes restent lisses suffisamment longtemps pour que cette théorie statistique soit mathématiquement fondée. Le problème d'existence et de régularité demande, en partie, si le mécanisme dissipatif de la viscosité est assez fort pour régulariser la cascade à toutes les échelles.

La différence entre les équations d'Euler et de Navier-Stokes se résume à un terme : $\nu \Delta u$. Mais ce terme change tout.

Euler n'est pas un Navier-Stokes simplifié. C'est un système fondamentalement différent qui partage l'essentiel de sa structure. La comparaison a également une portée pratique : choisir le mauvais modèle — Euler là où la viscosité compte, ou Navier-Stokes là où elle ne compte pas — peut faire la différence entre une simulation utile et une simulation inutile. Pour les équations complètes, voir Que sont les équations de Navier-Stokes ?. Pour les obstacles, voir Pourquoi c'est difficile. Pour le prix, voir Le problème du Millénaire. Pour incompressible contre compressible, voir Incompressible contre compressible.

Le terme visqueux $\nu \Delta u$ est une régularisation parabolique qui convertit un système non linéaire du premier ordre (Euler) en un système parabolique semi-linéaire (Navier-Stokes). Il fournit la dissipation d'énergie, des estimations a priori d'ordre supérieur et la structure de semi-groupe analytique qui sous-tend une grande partie de la théorie de Navier-Stokes (Leray 1934, Fujita–Kato 1964).

Pourtant, cette régularisation n'est pas suffisante en 3D pour fermer l'argument de régularité globale. Le problème du Millénaire de Clay demande précisément si le lissage parabolique dans Navier-Stokes peut contrôler le transfert non linéaire d'énergie à travers toutes les échelles et pour tout temps. La question parallèle pour Euler — si l'absence de tout lissage conduit à une explosion en temps fini à partir de données lisses — reste également ouverte et également fondamentale.

Les deux problèmes se situent au cœur de la théorie mathématique des fluides. La comparaison entre eux clarifie ce que la viscosité apporte et ce qu'elle n'apporte pas, et pourquoi la question de la régularité 3D, pour l'un ou l'autre système, figure parmi les problèmes ouverts les plus difficiles de l'analyse.