Non-uniqueness of smooth solutions of the Navier-Stokes equations from almost the same initial conditions

En utilisant la simulation numérique propre (CNS), cet article fournit des preuves numériques suggérant que les équations de Navier-Stokes admettent des solutions globales distinctes à partir de conditions initiales quasi identiques, offrant ainsi un éclairage nouveau sur le problème de l'unicité et de l'existence lié au prix du Millénaire.

L'essentiel

🌊 Le Mystère des Deux Rivières Qui Ne Se Ressemblent Pas

Imaginez que vous avez deux gouttes d'eau. Elles sont presque identiques. Si vous les laissez tomber dans un ruisseau calme, vous vous attendez à ce qu'elles suivent exactement le même chemin, non ? C'est ce que la physique classique nous a appris pendant des siècles : si vous connaissez le départ, vous connaissez l'arrivée.

Mais les mathématiciens et physiciens Liao et Qin, de l'Université Jiaotong de Shanghai, ont découvert quelque chose de fascinant (et de troublant) avec les équations de Navier-Stokes, qui décrivent comment les fluides (comme l'eau ou l'air) bougent.

Leur découverte ? Même si vous commencez avec deux gouttes d'eau presque indiscernables, elles peuvent finir par devenir deux rivières totalement différentes.

🦋 L'Effet Papillon : Un Souffle de Vent Qui Change le Monde

Pour comprendre cela, il faut parler de l'effet papillon. C'est l'idée qu'un tout petit battement d'ailes de papillon (une perturbation minuscule) peut, à la longue, provoquer une tempête à l'autre bout du monde.

Dans les fluides turbulents (comme l'eau qui dévale une cascade ou l'air autour d'une aile d'avion), les choses sont chaotiques.

• Le problème des ordinateurs classiques : Jusqu'à présent, les ordinateurs ne pouvaient pas prouver cela. Pourquoi ? Parce que les ordinateurs font des erreurs de calcul minuscules (comme un grain de poussière sur une balance). Dans un système chaotique, ce "grain de poussière" numérique grossit si vite qu'il fausse tout le résultat. C'est comme essayer de tracer la trajectoire d'un papillon avec un crayon qui tremble : le dessin devient illisible très vite.

🔍 La "Simulation Propre" (CNS) : Le Microscope Ultime

Pour résoudre ce problème, les auteurs ont utilisé une méthode spéciale qu'ils appellent la Simulation Numérique Propre (Clean Numerical Simulation ou CNS).

Imaginez que les ordinateurs classiques sont des lunettes de vue un peu rayées. Les auteurs ont créé des lunettes parfaitement nettes, avec une précision si incroyable qu'elles peuvent voir des différences de taille 10⁻⁴⁰.
(Pour se faire une idée : 10⁻⁴⁰, c'est comme comparer la taille d'un atome à la taille de tout l'univers observable, et trouver une différence encore plus petite !)

Grâce à cette précision extrême, ils ont pu simuler deux fluides qui commencent avec une différence infime, presque inexistante.

🎭 Le Résultat : Deux Mondes Différents

Voici ce qu'ils ont observé dans leur simulation :

1. Le Départ : Ils ont pris un fluide (appelons-le "Fluide A") et un autre fluide ("Fluide B") qui était identique au premier, sauf pour une toute petite modification, comme ajouter une pincée de sel invisible.
2. Le Chaos : Au début, les deux fluides bougent exactement de la même façon. Mais après un certain temps, la "pincée de sel" invisible s'amplifie grâce au chaos.
3. La Séparation : Soudainement, le Fluide A et le Fluide B prennent des chemins radicalement différents.
   • Le Fluide A garde une certaine symétrie (il tourne de manière régulière).
   • Le Fluide B perd cette symétrie et devient désordonné.
   • Leurs statistiques (leur "personnalité" globale) deviennent totalement différentes.

💡 Pourquoi est-ce important ?

Cette découverte est une bombe pour les mathématiques. Elle suggère que les équations de Navier-Stokes (qui sont l'un des sept problèmes du Prix du Millénaire, avec une récompense d'un million de dollars) pourraient avoir plusieurs solutions différentes pour un même départ.

C'est comme si vous lanciez deux balles de tennis avec la même force et le même angle, et que, par magie, l'une atterrissait dans le jardin du voisin et l'autre dans votre piscine, sans qu'aucun vent n'ait soufflé.

En résumé :
Les auteurs nous disent : "Attention ! Même si vous connaissez l'état initial d'un fluide avec une précision absolue, le futur n'est peut-être pas unique. Le chaos peut transformer une différence infime en deux mondes totalement opposés."

C'est une preuve numérique (très solide grâce à leur méthode "propre") qui pourrait aider les mathématiciens à résoudre l'un des plus grands mystères de la physique moderne.

Résumé technique

1. Problématique

L'article aborde l'un des problèmes du prix du millénaire de l'Institut Clay : l'existence et l'unicité des solutions lisses (régulières) des équations de Navier-Stokes (NS). Bien qu'il soit généralement admis que les écoulements turbulents sont régis par ces équations, la question de savoir si une condition initiale donnée conduit à une solution unique pour tout temps $t > 0$ reste ouverte.

Des travaux récents (Coiculescu et Palasek) ont démontré l'existence de deux solutions globales distinctes à partir de la même condition initiale dans l'espace critique $BMO^{-1}$, mais ces constructions nécessitent des données à énergie infinie. L'objectif de cet article est de fournir des exemples numériques explicites montrant que des conditions initiales presque identiques (avec une différence infinitésimale et une énergie finie) peuvent conduire à des solutions globales distinctes pour les équations de Navier-Stokes.

Le défi majeur réside dans la nature chaotique de la turbulence : les erreurs numériques traditionnelles (bruit numérique) s'amplifient exponentiellement (effet papillon), rendant impossible la distinction entre une divergence due à la physique du système et une divergence due aux artefacts numériques dans les simulations classiques (DNS - Direct Numerical Simulation).

2. Méthodologie

Pour surmonter les limitations des méthodes numériques traditionnelles, les auteurs utilisent la Simulation Numérique Propre (Clean Numerical Simulation - CNS).

• Principe de la CNS : Contrairement à la DNS, la CNS réduit rigoureusement les erreurs de troncature et d'arrondi à un niveau négligeable sur un intervalle de temps fini mais suffisamment long pour les statistiques. Cela est obtenu en utilisant :
  • Une précision multiple (nombre de chiffres significatifs $N_s$ très élevé).
  • Des développements de Taylor d'ordre élevé ($M$) pour l'approximation temporelle.
  • Une méthode pseudo-spectrale de Fourier pour l'approximation spatiale.
• Configuration du modèle :
  • Équations : Écoulement de Kolmogorov bidimensionnel incompressible dans un domaine carré périodique $[0, 2\pi] \times [0, 2\pi]$.
  • Paramètres : Nombre de Reynolds $Re = 2000$, échelle de forçage $n_K = 16$.
  • Conditions initiales :
    • $\psi_1$ : Une condition initiale possédant une symétrie de rotation et de translation.
    • $\psi_2$ : Une condition initiale perturbée par un terme $\delta \sin(x+y)$, où $\delta$ est une constante non nulle extrêmement petite ($10^{-10}$, $10^{-20}$, et $10^{-40}$).
  • Validation : Pour garantir la fiabilité, une solution de référence est calculée avec des paramètres de précision encore plus stricts ($M$ et $N_s$ accrus) pour confirmer que le bruit numérique de la solution principale est négligeable sur l'intervalle $t \in [0, 300]$.

3. Contributions Clés

• Preuve numérique de non-unicité : L'article fournit des preuves numériques montrant que des conditions initiales différant d'un ordre de grandeur aussi faible que $10^{-40}$ conduisent à des solutions globales distinctes.
• Utilisation de la CNS : Démonstration que la CNS permet d'isoler la sensibilité intrinsèque du système (chaos) du bruit numérique, permettant d'observer la divergence réelle des trajectoires.
• Analyse de la symétrie et des statistiques : Mise en évidence du fait que la perturbation infinitésimale détruit la symétrie spatiale initiale et modifie radicalement les propriétés statistiques de l'écoulement, bien que les conditions initiales soient quasi identiques.

4. Résultats

Les simulations ont été menées sur l'intervalle de temps $t \in [0, 300]$. Les résultats montrent :

• Divergence des trajectoires : Initialement, les écoulements issus de $\psi_1$ (Flow CNS) et de $\psi_2$ (Flow CNS'$_{1,2,3}$) sont indiscernables. Cependant, après un certain temps, ils divergent complètement.
• Taux de dissipation d'énergie cinétique : La moyenne spatiale du taux de dissipation d'énergie $\langle D \rangle_A$ pour le cas non perturbé (Flow CNS) devient environ deux fois plus grande que celle des cas perturbés après $t > 100$.
• Brisure de symétrie :
  • La solution non perturbée conserve la symétrie de rotation et de translation de sa condition initiale.
  • Les solutions perturbées, bien que la perturbation soit minuscule ($\delta \to 0$), perdent la symétrie de rotation et ne conservent que la symétrie de translation. Cette rupture est due à l'amplification exponentielle de la perturbation par le mécanisme de « cascade d'expansion du bruit » propre aux systèmes chaotiques.
• Différences statistiques : Les fonctions de densité de probabilité (PDF) du taux de dissipation et les moyennes spatio-temporelles diffèrent significativement entre le cas non perturbé et les cas perturbés, confirmant que les solutions appartiennent à des régimes statistiques distincts.
• Énergie finie : Toutes les solutions simulées possèdent une densité d'énergie finie, répondant ainsi à la contrainte de réalisme physique souvent absente dans les contre-exemples théoriques.

5. Signification et Conclusion

Ces résultats suggèrent fortement une conjecture de non-unicité pour les équations de Navier-Stokes :

Les équations de Navier-Stokes admettent des solutions globales lisses distinctes pour toutes $t > 0$, à partir de conditions initiales quasi identiques $U_1$ et $U_2$ telles que $|U_1 - U_2| \leq \delta$, pour un $\delta$ arbitrairement petit (ou tendant vers 0).

Implications :

• Ces exemples numériques servent de contre-exemples potentiels à l'unicité des solutions lisses à énergie finie.
• Ils offrent une « feuille de route » pour les mathématiciens cherchant à prouver rigoureusement la non-unicité dans le cadre du problème du prix du millénaire.
• Ils soulignent que la turbulence, en tant que système chaotique, est intrinsèquement sensible à des perturbations infinitésimales, au point que des conditions initiales indiscernables numériquement (mais mathématiquement distinctes) peuvent mener à des états physiques totalement différents.

L'article conclut que, grâce à la précision extrême de la CNS, ces divergences ne sont pas des artefacts numériques mais des propriétés fondamentales des équations de Navier-Stokes dans le régime turbulent.