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

Utilizzando la simulazione numerica pulita (CNS), l'articolo fornisce evidenze numeriche che le equazioni di Navier-Stokes ammettono soluzioni globali distinte a partire da condizioni iniziali quasi identiche, offrendo nuovi spunti sul problema del millennio relativo all'esistenza e all'unicità di tali soluzioni.

L'essenziale

Il Titolo: "Due mondi diversi partendo dallo stesso punto"

Immagina di avere due palline da biliardo quasi identiche. Sono così simili che, se le guardi con un microscopio, sembrano la stessa identica cosa. Se le lanci sul tavolo con la stessa forza e nella stessa direzione, ci aspetteremmo che facciano esattamente lo stesso percorso.

Questo è il modo in cui la fisica classica ci ha insegnato a pensare: stesso inizio = stesso risultato.

Tuttavia, questo articolo di Shijun Liao e Shijie Qin suggerisce che le equazioni che governano i fluidi turbolenti (come l'acqua che scorre veloce o l'aria che si muove in una tempesta) potrebbero non funzionare così. Potrebbero permettere che due palline quasi identiche finiscano in due posti completamente diversi, prendendo strade opposte.

Il Problema: Il "Farfalla" che distrugge i computer

Per capire perché nessuno aveva scoperto questo prima, dobbiamo parlare dei computer.
I fluidi turbolenti sono come un caos totale. Sono sensibili al minimo dettaglio. C'è un famoso detto: "Il battito d'ali di una farfalla in Brasile può scatenare un tornado in Texas". Questo è il "effetto farfalla".

Il problema è che i computer tradizionali (chiamati nel testo DNS) sono come bambini che cercano di fare calcoli complessi ma hanno sempre un po' di "rumore" di fondo. Immagina di dover copiare una frase perfetta, ma ogni volta che scrivi una lettera, fai un minuscolo errore di battitura (anche solo un punto di troppo).
In un sistema normale, questo errore rimane piccolo. Ma in un sistema caotico come la turbolenza, questo piccolo errore si ingrandisce esplosivamente, come una valanga che parte da un singolo sassolino e finisce per distruggere un villaggio.

Per decenni, i ricercatori hanno pensato: "Se il computer ci dà due risultati diversi partendo da quasi lo stesso punto, è colpa nostra! È solo un errore di calcolo, non è la realtà."

La Soluzione: La "Simulazione Pulita" (CNS)

Gli autori di questo articolo hanno usato una nuova tecnica chiamata CNS (Clean Numerical Simulation), ovvero "Simulazione Numerica Pulita".

Immagina la CNS non come un bambino che scrive, ma come un robot con una precisione sovrumana.

• I computer normali usano numeri con circa 15 cifre decimali.
• La CNS usa numeri con 260 cifre decimali.

È come se invece di scrivere "3,14", il computer scrivesse "3,1415926535..." fino a 260 posizioni dopo la virgola. In questo modo, il "rumore" o l'errore è così piccolo (ordine di $10^{-40}$, cioè un numero con 40 zeri dopo la virgola prima di un numero significativo) che è praticamente zero. È come se avessimo un microscopio così potente da vedere gli atomi, ma senza che la lente sia sporca.

L'Esperimento: Due fluidi, un destino diverso

Gli scienziati hanno creato un esperimento virtuale in un quadrato (come una vasca da bagno quadrata) con un fluido che scorre in modo turbolento.

1. Scenario A: Hanno lanciato il fluido con una certa configurazione iniziale (chiamiamola "Pallina A").
2. Scenario B: Hanno lanciato il fluido con una configurazione quasi identica alla A, ma con una differenza minuscola, quasi invisibile (chiamiamola "Pallina B"). La differenza era così piccola che era come aggiungere una goccia d'acqua a un oceano.

Cosa è successo?

• All'inizio, le due palline si muovevano insieme.
• Dopo un po' di tempo, la "Pallina B" ha iniziato a deviare.
• Dopo un po' di più, le due palline erano su percorsi completamente diversi.

Non solo si sono spostate in posti diversi, ma hanno assunto forme diverse:

• La "Pallina A" ha mantenuto una simmetria rotazionale (come una ruota che gira perfettamente).
• La "Pallina B" ha rotto questa simmetria e ha assunto una forma diversa.

E la cosa più incredibile? Anche le statistiche (la media di quanto il fluido si muoveva) erano diverse. Non era solo un piccolo errore di percorso; erano due stati fisici completamente differenti.

Perché è importante? (Il Premio da un Milione di Dollari)

Questo risultato è una bomba perché tocca uno dei problemi più famosi della matematica moderna: il Problema del Millennio del Clay Institute.

Uno dei sette problemi irrisolti della matematica chiede: "Le equazioni di Navier-Stokes hanno sempre una soluzione unica e liscia?"
In parole povere: "Se conosco lo stato di un fluido oggi, posso prevedere esattamente come sarà domani?"

Se la risposta è "Sì, c'è sempre una sola soluzione", allora il mondo è prevedibile (almeno in teoria).
Se la risposta è "No, possono esserci più soluzioni partendo dallo stesso punto", allora la prevedibilità ha un limite fondamentale, anche con computer perfetti.

Questo articolo non dice "abbiamo risolto il problema", ma dice: "Ehi, guardate qui! Con i nostri computer perfetti, abbiamo visto che due fluidi quasi identici fanno cose diverse. Forse la risposta al problema del Millennio è che la soluzione unica non esiste sempre."

In Sintesi

Immagina di avere due fogli di carta quasi identici. Se li lanci in una stanza piena di correnti d'aria caotiche, ci aspetteremmo che volino insieme. Questo studio dice che, se misuriamo con una precisione assoluta, potremmo scoprire che un foglio finisce sul soffitto e l'altro sul pavimento, non perché c'era un vento diverso, ma perché le leggi della fisica stessa permettono che due inizi quasi uguali portino a destini completamente diversi.

È una scoperta che potrebbe cambiare il modo in cui pensiamo alla prevedibilità del nostro universo.

Sommario tecnico

Titolo: Non-univocità delle soluzioni lisce delle equazioni di Navier-Stokes da condizioni iniziali quasi identiche

Autori: Shijun Liao e Shijie Qin
Affiliazione: Shanghai Jiaotong University, Cina

---

1. Il Problema

Il lavoro affronta uno dei problemi del Millennio dell'Istituto Clay: l'esistenza e l'univocità delle soluzioni lisce (regolari) per le equazioni di Navier-Stokes (NS) in tre dimensioni (e in questo caso specifico, in un contesto bidimensionale turbolento).
Sebbene sia ampiamente accettato che i flussi turbolenti siano governati dalle equazioni di NS, la questione se queste ammettano soluzioni uniche per dati iniziali assegnati rimane aperta.

• Contesto teorico: Recenti lavori (es. Coiculescu e Palasek) hanno dimostrato l'esistenza di dati iniziali nello spazio critico $BMO^{-1}$ che generano due soluzioni globali distinte, ma tali costruzioni richiedono dati con energia infinita.
• La sfida numerica: L'obiettivo degli autori è fornire esempi numerici espliciti di non-univocità per soluzioni con energia finita, partendo da condizioni iniziali quasi identiche (differenza dell'ordine di $10^{-40}$).
• Ostacolo tradizionale: I metodi numerici classici, come la Simulazione Numerica Diretta (DNS), falliscono in questo compito a causa della "sensibilità alle condizioni iniziali" (effetto farfalla) nei sistemi caotici. Il rumore numerico (errori di troncamento e di arrotondamento) cresce esponenzialmente, inquinando rapidamente la traiettoria esatta e rendendo impossibile distinguere tra divergenza fisica e errore numerico.

2. Metodologia: Simulazione Numerica Pulita (CNS)

Per superare le limitazioni dei metodi tradizionali, gli autori utilizzano la Clean Numerical Simulation (CNS), una metodologia sviluppata per garantire la precisione delle traiettorie spaziotemporali in intervalli di tempo finiti ma sufficientemente lunghi per l'analisi statistica.

• Equazioni Goveranti: Flusso di Kolmogorov bidimensionale incomprimibile in un dominio quadrato $[0, 2\pi] \times [0, 2\pi]$ con condizioni al contorno periodiche.
  • Numero di Reynolds: $Re = 2000$.
  • Scala di forzante di Kolmogorov: $n_K = 16$.
• Condizioni Iniziali:
  • Caso A (Flusso CNS): $\psi_1(x, y, 0) = -\frac{1}{2}[\cos(x+y) + \cos(x-y)]$. Questa condizione possiede simmetria di rotazione e traslazione.
  • Caso B (Flussi CNS'): $\psi_2(x, y, 0) = \psi_1(x, y, 0) + \delta \sin(x+y)$.
  • Vengono testati tre valori estremamente piccoli per $\delta$: $10^{-10}$, $10^{-20}$ e $10^{-40}$.
• Implementazione CNS:
  • Metodo spaziale: Pseudo-spettrale di Fourier su una griglia uniforme $1024 \times 1024$.
  • Metodo temporale: Espansione di Taylor di ordine $M$ con passo temporale $\Delta t = 10^{-3}$.
  • Precisione: Utilizzo di aritmetica a precisione multipla (con $N_s$ cifre significative) per ridurre gli errori di arrotondamento a livelli trascurabili.
  • Validazione: Per garantire l'affidabilità, vengono eseguite simulazioni di controllo con parametri ancora più rigorosi ($M$ e $N_s$ maggiori) per definire il "tempo prevedibile critico" ($T_c$). Nel lavoro, $T_c$ è stato stabilito fino a $t=300$.
  • Risorse computazionali: Calcoli eseguiti sul supercomputer "Tianhe".

3. Risultati Chiave

Le simulazioni CNS hanno rivelato divergenze significative tra le soluzioni derivanti da condizioni iniziali quasi indistinguibili:

1. Divergenza Temporale:

   • Inizialmente, i flussi derivanti da $\psi_1$ e $\psi_2$ (con $\delta$ piccolissimo) appaiono identici.
   • Dopo un certo periodo (intorno a $t > 100$), le soluzioni si separano drasticamente.
   • Il tasso di dissipazione dell'energia cinetica media spazialmente $\langle D \rangle_A$ del flusso CNS (condizione iniziale simmetrica) diventa circa due volte superiore rispetto ai flussi CNS' (condizioni con perturbazione $\delta$).

2. Rottura della Simmetria Spaziale:

   • La soluzione CNS mantiene la simmetria di rotazione e traslazione della sua condizione iniziale.
   • Le soluzioni CNS' perdono la simmetria di rotazione, conservando solo quella di traslazione.
   • Questo dimostra che la piccola perturbazione $\delta \sin(x+y)$, a causa del caos intrinseco della turbolenza, cresce fino a un livello macroscopico ("noise-expansion cascade"), distruggendo completamente la simmetria originale.

3. Differenze Statistiche:

   • Le distribuzioni di probabilità (PDF) del tasso di dissipazione dell'energia e le medie spazio-temporali sono statisticamente diverse tra il caso CNS e i casi CNS', nonostante le condizioni iniziali differiscano solo di $10^{-40}$.

4. Contributi Principali

• Evidenza Numerica di Non-Univocità: Forniscono i primi esempi numerici espliciti in cui le equazioni di Navier-Stokes ammettono soluzioni globali lisce e distinte partendo da condizioni iniziali con differenza arbitrariamente piccola ($\delta \to 0$) e energia finita.
• Validazione della CNS: Dimostrano l'efficacia della CNS nel distinguere tra rumore numerico e comportamento fisico reale in sistemi caotici su scale temporali lunghe, superando i limiti della DNS tradizionale.
• Meccanismo di Rottura: Illustrano come il "cascata di espansione del rumore" nei sistemi caotici possa trasformare micro-perturbazioni in macro-differenze strutturali (come la rottura di simmetria).

5. Significato e Implicazioni

Questo lavoro ha profonde implicazioni per la matematica pura e la fisica dei fluidi:

• Contro-esempio alla Congettura di Unicità: Sebbene i risultati siano numerici e limitati a un intervallo di tempo finito e a un caso 2D specifico, fungono da potenti contro-esempi che suggeriscono che la congettura di unicità delle soluzioni lisce di Navier-Stokes potrebbe essere falsa anche per dati iniziali con energia finita.
• Roadmap per la Matematica: Gli autori propongono una congettura basata sui loro risultati: Le equazioni di Navier-Stokes ammettono soluzioni globali lisce distinte per tutte $t > 0$ da condizioni iniziali quasi identiche. Questo offre una direzione di ricerca per i matematici che cercano di provare o confutare l'univocità in modo analitico.
• Comprensione della Turbolenza: Sottolinea la natura intrinsecamente caotica della turbolenza, dove la prevedibilità a lungo termine è limitata non solo dalla fisica, ma anche dalla impossibilità pratica di definire condizioni iniziali con precisione infinita, rendendo la distinzione tra soluzioni diverse una questione fondamentale.

In sintesi, il paper utilizza la potenza computazionale e la precisione matematica della CNS per sfidare l'assunzione convenzionale di unicità nelle equazioni di Navier-Stokes, fornendo evidenze numeriche robuste che le soluzioni possono divergere drasticamente anche da perturbazioni iniziali infinitesimali.