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

Met behulp van schone numerieke simulaties (CNS) levert dit artikel numeriek bewijs dat de Navier-Stokes-vergelijkingen verschillende globale oplossingen kunnen hebben voor bijna identieke beginvoorwaarden, zelfs bij verschillen zo klein als $10^{-40}$, wat relevant is voor het Millennium Prize-probleem over uniciteit en bestaan.

De kern

De Navier-Stokes-riddle: Hoe een rimpel in een meer een tsunami kan worden

Stel je voor dat je twee identieke kopjes koffie hebt. Ze staan naast elkaar, ze hebben precies dezelfde temperatuur, dezelfde hoeveelheid suiker en ze zijn net zo stil. Je zou denken dat ze zich in de toekomst precies hetzelfde zouden gedragen, toch?

In de wereld van de wiskunde en de natuurkunde, waar stromende vloeistoffen (zoals water of lucht) worden beschreven door de beroemde Navier-Stokes-vergelijkingen, is het antwoord misschien wel verrassend: Nee.

Dit is wat de onderzoekers Liao en Qin in hun nieuwe paper ontdekken. Ze tonen aan dat je twee vloeistofstromen kunt hebben die op het begin bijna identiek zijn – zo bijna dat het verschil kleiner is dan een stofje in een universum – maar die toch totaal verschillende paden afleggen.

Hier is hoe ze dit hebben bewezen, vertaald in alledaagse taal:

1. Het probleem: De "Vlinder" en de computer

We weten al lang dat vloeistoffen chaotisch zijn. Een klein beetje beweging (een vlinder die met zijn vleugels slaat) kan later leiden tot een grote storm. Dit noemen we het vlindereffect.

Het probleem voor wiskundigen is dat computers niet perfect zijn. Computers maken altijd kleine rekenfoutjes (zoals als je een getal afrondt). Bij gewone computersimulaties van vloeistoffen worden deze kleine foutjes door het chaotische gedrag van de stroming enorm opgeblazen. Het is alsof je een foto probeert te kopiëren, maar elke keer dat je kopieert, wordt het beeld een beetje waziger. Na een tijdje is de kopie niet meer te onderscheiden van de echte foto.

Dit maakte het voor wetenschappers onmogelijk om te zeggen: "Kijk, dit is echt een ander pad dan dat andere, of is het gewoon een computerfout?"

2. De oplossing: De "Schoonmaak"-computer (CNS)

De onderzoekers gebruikten een speciale techniek die ze "Clean Numerical Simulation" (CNS) noemen. Je kunt dit zien als een superkrachtige rekenmachine die niet alleen heel nauwkeurig rekent, maar ook constant controleert of er geen "stof" (rekenfoutjes) op het lensje zit.

Ze gebruikten zo'n veel decimalen (cijfers achter de komma) dat de rekenfouten zo klein waren dat ze in de praktijk niet bestonden. Hierdoor konden ze een simulatie draaien die zo dicht bij de "ware" natuur zat, dat ze met zekerheid konden zeggen: "Dit is geen computerfout. Dit is echt wat er gebeurt."

3. Het experiment: Twee bijna-identieke startpunten

Ze startten twee simulaties:

• Simulatie A: Een perfecte, symmetrische start.
• Simulatie B: Dezelfde start, maar met een minuscule verstoring.

Deze verstoring was zo klein dat het verschil in grootte ongeveer 1 op de 10.000.000.000.000.000.000.000.000.000.000.000.000.000 was. Dat is zo klein dat het in de echte wereld onmeetbaar zou zijn. Het is alsof je aan een reus een stofje toevoegt.

Het resultaat?
In het begin leken de twee stromingen identiek. Maar na een tijdje (in hun simulatie na ongeveer 100 tijdseenheden) gebeurde er iets wonderlijks:

• Simulatie A bleef een mooie, symmetrische dans uitvoeren.
• Simulatie B begon totaal anders te bewegen. De symmetrie brak volledig. De stroming werd chaotisch en zag er compleet anders uit.

Het kleine stofje (de verstoring) was door de chaos van de stroming zo snel opgeblazen tot een enorme tsunami, dat het de hele dans van de vloeistof veranderde.

4. Waarom is dit belangrijk?

Dit is een enorme stap voor de wiskunde. Er is een beroemde prijs (de Millennium Prize) uitgeschreven voor iemand die kan bewijzen of de Navier-Stokes-vergelijkingen altijd maar één oplossing hebben voor een bepaalde start.

De onderzoekers zeggen met hun experiment: "Kijk, we hebben twee startpunten die bijna identiek zijn, maar ze leiden tot twee totaal verschillende oplossingen."

Dit suggereert dat de Navier-Stokes-vergelijkingen misschien niet uniek zijn. Dat wil zeggen: als je de natuur precies zo'n start geeft, is het niet altijd voorspelbaar welke uitkomst je krijgt. Het is alsof je twee identieke dobbelstenen gooit, maar ze landen altijd op een ander getal, puur door een onmeetbaar klein verschil in de luchtstroom.

Samenvatting

Dit papier toont aan dat in de wereld van turbulente vloeistoffen, een verschil dat zo klein is dat het onbestaanbaar lijkt, toch de hele toekomst van een stroming kan veranderen. Met hun "schoonmaak-computer" hebben ze bewezen dat de wiskunde misschien niet zo strak en voorspelbaar is als we hoopten. Het is een ontdekking die de deur opent naar een nieuw begrip van hoe chaos en orde in onze wereld werken.

Technische samenvatting

Probleemstelling

Het paper adresseert een van de zeven Millennium Prize Problems van het Clay Mathematics Institute: de uniekheid en bestaan van gladde oplossingen voor de Navier-Stokes-vergelijkingen (NS-vergelijkingen). Hoewel het algemeen wordt aanvaard dat turbulente stromingen door deze vergelijkingen worden beschreven, is het wiskundig nog niet bewezen of de oplossing uniek is voor gegeven beginvoorwaarden, zelfs als die beginvoorwaarden eindige energie hebben.

Bestaand onderzoek (zoals Coiculescu en Palasek) heeft aangetoond dat er verschillende globale oplossingen kunnen bestaan vanuit dezelfde beginvoorwaarden, maar deze constructies vereisen vaak data met oneindige energie of zijn puur theoretisch. De centrale vraag van dit paper is: Kunnen we numerieke voorbeelden geven van niet-uniekheid voor NS-vergelijkingen met eindige energie, waarbij de beginvoorwaarden willekeurig dicht bij elkaar liggen, maar de oplossingen fundamenteel verschillend zijn?

Traditionele numerieke methoden (zoals Direct Numerical Simulation, DNS) falen hierbij omdat numerieke ruis (afkappingsfouten en afrondingsfouten) door het "vlinder-effect" in chaotische systemen exponentieel groeit. Hierdoor verlaten DNS-resultaten snel de ware trajecten van de vergelijkingen, wat het onmogelijk maakt om echte niet-uniekheid te onderscheiden van numerieke artefacten.

Methodologie: Clean Numerical Simulation (CNS)

Om de beperkingen van traditionele methoden te overwinnen, maken de auteurs gebruik van Clean Numerical Simulation (CNS). Deze methode is ontworpen om numerieke ruis strikt verwaarloosbaar te houden over een eindig, maar lang genoeg tijdsinterval voor statistische analyse.

De kern van de methodologie omvat:

• Vergelijkingen: Een tweedimensionale, incompressibele Kolmogorov-stroming in een vierkante domein $[0, 2\pi] \times [0, 2\pi]$ met periodieke randvoorwaarden. De stroming wordt gestuurd door de dimensieloze NS-vergelijking in de vorm van een stroomfunctie $\psi$ met Reynoldsgetal $Re = 2000$ en Kolmogorov-schaal $n_K = 16$.
• Numerieke Techniek:
  • Ruimtelijke discretisatie: Fourier-pseudospectrale methode op een uniform rooster van $1024 \times 1024$.
  • Tijdsintegratie: Taylor-ontwikkeling van de $M$-de orde met een tijdstap $\Delta t = 10^{-3}$.
  • Meerdere precisie (Multiple-precision): In tegenstelling tot DNS, worden alle variabelen en parameters berekend met een groot aantal significante cijfers ($N_s$) om afrondingsfouten extreem te verkleinen.
• Validatie: Om de "kritische voorspelbare tijd" ($T_c$) te bepalen, worden simulaties uitgevoerd met nog kleinere numerieke ruis (hogere $M$ en $N_s$). Als de resultaten overeenkomen binnen het interval $t \in [0, 300]$, wordt de oplossing beschouwd als een convergente referentie die zeer dicht bij de exacte oplossing ligt.
• Computatie: De berekeningen zijn uitgevoerd op de "Tianhe New Generation Supercomputer" in China.

Experimenteel Opzet

De auteurs vergelijken vier simulaties met bijna identieke beginvoorwaarden:

1. Flow CNS: Beginvoorwaarde $\psi_1(x,y,0) = -\frac{1}{2}[\cos(x+y) + \cos(x-y)]$. Deze heeft zowel rotatie- als translatiesymmetrie.
2. Flow CNS'$_1$, CNS'$_2$, CNS'$_3$: Beginvoorwaarde $\psi_2(x,y,0) = \psi_1(x,y,0) + \delta \sin(x+y)$, waarbij $\delta$ respectievelijk $10^{-10}$, $10^{-20}$ en $10^{-40}$ is.
   • Hoewel $\delta$ extreem klein is (orde $10^{-40}$), breekt deze perturbatie de rotatiesymmetrie van de oorspronkelijke voorwaarde, hoewel de translatiesymmetrie behouden blijft.

Belangrijkste Resultaten

De simulaties tonen aan dat de Navier-Stokes-vergelijkingen distincte globale oplossingen toelaten vanuit beginvoorwaarden die willekeurig dicht bij elkaar liggen:

1. Afwijking in Trajecten: Hoewel de oplossingen aanvankelijk identiek lijken, beginnen ze na een bepaalde tijd (rond $t > 100$) fundamenteel van elkaar af te wijken door het chaotische karakter van de turbulentie. De "noise-expansion cascade" zorgt ervoor dat de microscopische verstoring $\delta$ exponentieel groeit tot een macroscopisch niveau.
2. Verschil in Statistieken:
   • De ruimtelijk gemiddelde kinetische energiedissipatie $\langle D \rangle_A$ van Flow CNS is na $t=100$ bijna twee keer zo groot als die van de andere drie stromingen.
   • De kansdichtheidsfuncties (PDF's) van de energiedissipatie en de ruimtetijd-gemiddelde dissipatie tonen duidelijke verschillen tussen de groep met de symmetrische beginvoorwaarde en de groep met de verstoorde beginvoorwaarde.
3. Verschil in Symmetrie:
   • Flow CNS behoudt de initiële rotatie- en translatiesymmetrie.
   • Flow CNS'$_{1,2,3}$ verliezen de rotatiesymmetrie en behouden alleen de translatiesymmetrie. Dit gebeurt ondanks dat de initiële verstoring $\delta$ verwaarloosbaar klein was ($10^{-40}$).
4. Eindige Energie: Alle gegenereerde oplossingen hebben een eindige energiedichtheid, wat hen onderscheidt van eerdere theoretische voorbeelden met oneindige energie.

Bijdrage en Betekenis

De belangrijkste bijdrage van dit paper is het leveren van expliciete numerieke bewijzen dat de Navier-Stokes-vergelijkingen niet-uniek kunnen zijn voor gladde oplossingen met eindige energie, zelfs bij beginvoorwaarden die willekeurig dicht bij elkaar liggen.

• Wiskundige Implicatie: De resultaten suggereren een nieuwe conjectuur: Er bestaan verschillende gladde globale oplossingen voor $t > 0$ vanuit beginvoorwaarden $U_1$ en $U_2$ waarbij $|U_1 - U_2| \leq \delta$ voor willekeurig klein $\delta$.
• Rol van CNS: Het paper demonstreert dat CNS een essentieel hulpmiddel is om de grenzen van de numerieke voorspelbaarheid van chaotische systemen te onderzoeken en om echte fysieke fenomenen te onderscheiden van numerieke ruis.
• Roadmap: Hoewel de resultaten numeriek zijn en beperkt tot een specifiek tijdsinterval en een 2D-systeem, bieden ze een "roadmap" voor wiskundigen om de uniekheid van de Navier-Stokes-vergelijkingen formeel te onderzoeken of te weerleggen.

Kortom, het paper stelt dat de gevoeligheid voor beginvoorwaarden in turbulente stromingen zo extreem is dat een infinitesimale verandering (zoals $10^{-40}$) kan leiden tot een volledig andere fysieke toestand met andere statistische eigenschappen en symmetrieën, wat de uniekheid van de oplossing in twijfel trekt.