Perché Navier-Stokes in 2D è più facile che in 3D

Le equazioni di Navier-Stokes descrivono come si muovono i fluidi. Funzionano in 2D (piatto, come l’acqua che si sparge su un tavolo) e in 3D (nella vita reale, come le correnti oceaniche che turbinano attorno a un sottomarino o il vento che sferza un grattacielo). Stesse equazioni. Quasi identiche.

Ecco il punto sorprendente. In 2D, i matematici possono dimostrare che le equazioni si comportano sempre bene, che la matematica non si rompe mai, che le soluzioni restano lisce per sempre. In 3D? Nessuno lo sa. Neppure una persona sulla Terra. Il fluido potrebbe fare qualcosa di così violento e improvviso da far smettere del tutto di funzionare la matematica, e dimostrare se ciò possa accadere è il problema del Millennio del Clay, del valore di un milione di dollari.

Non è semplicemente che “il 3D è più difficile perché c’è più roba”. Un meccanismo specifico presente in 3D non esiste in 2D. E cambia tutto.

Le equazioni di Navier-Stokes incomprimibili su $\mathbb{R}^n$ (o su un dominio periodico $\mathbb{T}^n$) con $\nu > 0$ sono

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

Per $n = 2$, l’esistenza globale e l’unicità di soluzioni classiche per dati iniziali lisci a divergenza nulla è un teorema. I riferimenti chiave sono Ladyzhenskaya (1959), che si basa sul lavoro precedente di Leray (1934). Il risultato si estende a soluzioni lisce su $\mathbb{R}^2$, $\mathbb{T}^2$ e domini limitati con condizioni al bordo standard.

Per $n = 3$, l’esistenza globale di soluzioni classiche a partire da dati lisci arbitrari è aperta. Questo è il contenuto del problema del Millennio del Clay come formulato da Fefferman (2000). Leray (1934) stabilì l’esistenza globale di soluzioni deboli in $L^2(\mathbb{R}^3)$, ma l’unicità e la regolarità di queste soluzioni restano irrisolte.

Il divario tra $n = 2$ e $n = 3$ non è una questione di contabilità. Riflette una differenza strutturale nell’equazione della vorticità, nelle proprietà di scala della non linearità e nelle stime a priori disponibili. Ciascuno di questi aspetti è esaminato nelle sezioni che seguono.

La domanda da un milione di dollari riguarda solo il 3D. Perché? Perché il 2D è risolto. Finito. I matematici hanno dimostrato decenni fa che le soluzioni bidimensionali di Navier-Stokes restano sempre lisce, qualunque siano le condizioni iniziali che si scelgono, per quanto a lungo si aspetti. Nessun premio necessario per un problema già risolto.

Quindi la vera domanda non è “perché il 3D è difficile?”. È “perché il 2D è facile e il 3D difficile?” Che cosa si rompe esattamente quando si aggiunge quella terza dimensione?

Il problema del Clay (Fefferman 2000) considera il problema di Cauchy per Navier-Stokes incomprimibili su $\mathbb{R}^3$ con viscosità $\nu > 0$ e dati iniziali lisci, a divergenza nulla, $u_0 \in C^\infty(\mathbb{R}^3)$ che soddisfano opportune condizioni di decadimento. La domanda: esiste una soluzione liscia $u \in C^\infty(\mathbb{R}^3 \times [0,\infty))$ con energia limitata?

L’enunciato analogo per $\mathbb{R}^2$ è un teorema. Ladyzhenskaya dimostrò l’esistenza globale e l’unicità di soluzioni forti per Navier-Stokes 2D con $u_0 \in H^1(\mathbb{R}^2)$, e un argomento di bootstrap dà la regolarità $C^\infty$ per dati lisci. La dimostrazione si basa su stime a priori specifiche delle due dimensioni che non si estendono a tre.

Il problema del Millennio è dunque puramente tridimensionale. I risultati parziali in 3D (le soluzioni deboli di Leray, 1934; il teorema di regolarità parziale di Caffarelli-Kohn-Nirenberg, 1982; vari criteri di regolarità condizionata) si fermano tutti prima di risolvere la questione completa della regolarità globale. Ogni risultato parziale mette in evidenza una lacuna specifica nel nostro controllo delle soluzioni 3D.

Il 2D ha un’arma segreta. Si chiama vorticità: quanto il fluido ruota in ogni punto.

In 2D, la vorticità è semplicemente un numero. Tutto qui. Senso orario o antiorario, veloce o lenta. Ed ecco che cosa rende le due dimensioni così notevolmente diverse dalle tre: questi piccoli vortici possono spostarsi nel fluido e affievolirsi gradualmente a causa dell’attrito, ma non possono mai, in nessuna circostanza, diventare più forti di quanto fossero all’inizio. Rotazione massima al tempo zero? È la rotazione massima che vedrai mai.

Perché questo conta? Tutto segue da qui. La velocità resta liscia. La pressione resta liscia. La soluzione continua a funzionare per sempre, per quanto assurdamente lontano nel futuro si vada, perché quell’unico vincolo sulla vorticità agisce come la prima tessera di un domino che fa cadere tutte le altre.

In 2D, la vorticità $\omega = \partial_1 u_2 - \partial_2 u_1$ è uno scalare che soddisfa

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

Questa è un’equazione di trasporto-diffusione per la vorticità scalare $\omega$. Il sistema completo resta non lineare perché $u$ si ricava da $\omega$ tramite la legge di Biot-Savart $u = \nabla^\perp (-\Delta)^{-1} \omega$, ma l’equazione non ha alcun termine sorgente di stretching dei vortici: il membro destro contiene solo il Laplaciano, non il termine $(\omega \cdot \nabla)u$ che compare in 3D.

Il principio del massimo scalare si applica direttamente: $\|\omega(t)\|_{L^\infty} \leq \|\omega_0\|_{L^\infty}$ per ogni $t \geq 0$. Allo stesso tempo, le norme $L^p$ di $\omega$ sono non crescenti per ogni $1 \leq p \leq \infty$.

Dal limite $L^\infty$ su $\omega$, si recupera $u \in L^\infty(0,T; W^{1,p})$ per ogni $p < \infty$ via Calderón-Zygmund estimates on the Biot-Savart kernel. Higher regularity follows by differentiating the vorticity equation and applying parabolic bootstrap: each spatial derivative of $\omega$ satisfies a parabolic equation with controllable coefficients, so bounds propagate to all orders.

Il principio del massimo per la vorticità in 2D ha un’origine più antica nel caso inviscido. Wolibner (1933) dimostrò l’esistenza globale per Euler 2D in spazi di Hölder, e Yudovich (1963) stabilì l’unicità per soluzioni di Euler 2D a vorticità limitata. Con viscosità ($\nu > 0$), lo smoothing parabolico rafforza soltanto queste stime. La dimostrazione di Ladyzhenskaya della regolarità globale per Navier-Stokes 2D si basa su questa struttura, combinata con la disuguaglianza di interpolazione di Ladyzhenskaya $\|f\|_{L^4}^2 \leq C \|f\|_{L^2} \|\nabla f\|_{L^2}$ (valida in 2D, con una struttura degli esponenti diversa da quella della controparte 3D).

In 3D, la vorticità non è un numero. È un vettore, che porta con sé sia una direzione sia un’intensità, e dovresti immaginarla come minuscoli tubi di tornado che attraversano il fluido.

Ecco che cosa rovina tutto. Quei tubi possono essere stirati. Tirane uno come caramella mou: si assottiglia e ruota più velocemente. Molto, molto più velocemente. Questo è lo stretching dei vortici, ed è il cattivo dell’intera storia perché significa che il fluido può amplificare la propria rotazione, alimentando scale sempre più piccole finché, forse, la rotazione in un singolo punto diventa infinitamente intensa.

Questo è un blowup. La matematica si rompe.

La viscosità (l’attrito interno del fluido) riesce sempre a frenare bruscamente lo stretching prima che raggiunga l’infinito, oppure a volte lo stretching sopraffà l’attrito e vince? Nessuno lo sa. Questa è, letteralmente, la domanda da un milione di dollari. Questo tiro alla fune tra stretching e attrito è il motivo per cui il problema è così difficile.

In 3D, la vorticità $\omega = \nabla \times u$ soddisfa

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

Il termine $(\omega \cdot \nabla)u$ è il termine di stretching dei vortici, assente in 2D. È quadratico nel senso che $u$ si ricava da $\omega$ tramite la legge di Biot-Savart 3D, quindi $(\omega \cdot \nabla)u$ scala grosso modo come $|\omega|^2$ nel caso peggiore. Questo termine permette la crescita di $\|\omega\|_{L^\infty}$ e distrugge il principio del massimo scalare disponibile in 2D.

Per Euler 3D, il criterio di Beale-Kato-Majda (1984) afferma che una soluzione liscia su $[0, T)$ si estende oltre il tempo $T$ se e solo se

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

Criteri di continuazione analoghi valgono per Navier-Stokes 3D. Il blowup, se si verifica, richiede che $\|\omega\|_{L^\infty}$ diventi non integrabile nel tempo. Il termine di stretching è la sorgente nell’equazione della vorticità che può amplificare la vorticità e impedire un argomento basato sul principio del massimo. In 2D, $\|\omega\|_{L^\infty}$ è limitata dai dati iniziali per ogni tempo; in 3D, controllare questa norma è il problema aperto centrale.

Progresso parziale: Caffarelli-Kohn-Nirenberg (1982) ha mostrato che, per ogni soluzione debole adatta di Navier-Stokes 3D, l’insieme dei punti singolari nello spaziotempo ha misura di Hausdorff parabolica unidimensionale nulla. Le singolarità, se esistono, sono estremamente rarefatte. Ma il teorema non ne esclude l’esistenza.

Per il caso inviscido, Elgindi (2021) ha dimostrato la formazione di singolarità in tempo finito per Euler 3D con dati iniziali $C^{1,\alpha}$ ($\alpha$ piccolo), usando un meccanismo guidato dallo stretching dei vortici lungo un asse di simmetria. Questo non implica direttamente il blowup per Navier-Stokes (la viscosità potrebbe comunque regolarizzare), ma mostra che il meccanismo di stretching è abbastanza forte da produrre singolarità senza smorzamento viscoso.

Lo stretching dei vortici non è l’unico problema. C’è una ragione strutturale più profonda per cui il 3D resiste alla dimostrazione, e si manifesta quando si fa \"zoom\" sul fluido.

Le equazioni di Navier-Stokes hanno un trucco di zoom. Prendi una soluzione qualsiasi, ingrandisci una regione più piccola, accelera il tempo della giusta quantità, e ottieni un’altra soluzione perfettamente valida. Quindi: che cosa succede alla stima di energia quando fai zoom?

• In 2D, ingrandire mantiene invariata la scala dell'energia. I matematici chiamano questo comportamento critico rispetto alla scala. Le stime di energia funzionano a ogni scala. Grande o piccola, non si perde mai il controllo.
• In 3D, il controllo energetico diventa più debole come controllo delle piccole scale. Questo è un comportamento supercritico rispetto alla scala, ed è devastante: le stime note perdono presa proprio alle scale in cui un blowup si concentrerebbe.

Un'analogia. In 2D, la tua torcia è sempre abbastanza luminosa. In 3D, più piccolo è ciò che osservi, più fioca diventa, e il fluido diventa più difficile da risolvere. Finisci al buio.

Non è un semplice inconveniente tecnico che un trucco ingegnoso potrebbe risolvere. È un muro. Gli strumenti matematici standard non riescono a controllare Navier-Stokes 3D alle piccole scale. Serve qualcosa di fondamentalmente nuovo.

Le equazioni di Navier-Stokes sono invarianti rispetto alla riscalatura

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

per ogni $\lambda > 0$. La norma $L^2$ si trasforma come $\|u_\lambda\|_{L^2(\mathbb{R}^n)} = \lambda^{1-n/2} \|u\|_{L^2}$.

• Per $n = 2$: $\|u_\lambda\|_{L^2} = \|u\|_{L^2}$. Invariante per scala. L'equazione è critica per l'energia.
• Per $n = 3$: $\|u_\lambda\|_{L^2} = \lambda^{-1/2} \|u\|_{L^2}$, che diminuisce quando $\lambda \to \infty$ (ingrandendo). L'equazione è supercritica rispetto all'energia: una norma di energia limitata diventa un controllo più debole alle piccole scale.

Lo spazio di Sobolev critico per Navier-Stokes 3D è $\dot{H}^{1/2}(\mathbb{R}^3)$, invariante rispetto alla riscalatura naturale. Ma l'identità di energia controlla solo $u$ in $L^\infty_t L^2_x \cap L^2_t \dot{H}^1_x$. È mezzo derivato sotto il livello critico. Questo è il divario di supercriticità.

In 2D, l'identità di energia $\frac{d}{dt}\|u\|_{L^2}^2 = -2\nu\|\nabla u\|_{L^2}^2$ fornisce esattamente il controllo al livello critico necessario; combinata con il principio del massimo per la vorticità, dà abbastanza regolarità per fare bootstrap fino a $C^\infty$. In 3D, la stessa identità produce una stima più debole del livello critico. Non è nota alcuna ulteriore stima a priori capace di colmare il divario.

Tao ha sottolineato questa barriera. Non ci si aspetta che argomenti basati solo sulla disuguaglianza di energia e sulla riscalatura risolvano la regolarità globale in 3D, quindi qualsiasi dimostrazione riuscita dovrà probabilmente sfruttare una struttura aggiuntiva, qualcosa oltre ciò che la sola analisi di scala può vedere. I metodi “ciechi alla supercriticità” (che trattano l'equazione solo attraverso la sua scala e la sua struttura energetica) non possono avere successo. La specifica struttura algebrica dell'equazione, in particolare la condizione di divergenza nulla e la struttura antisimmetrica della non linearità, dovrebbe giocare un ruolo. Si veda Perché Navier-Stokes è difficile per una trattazione più approfondita.

La dimostrazione in 2D funziona perché la vorticità resta limitata e la riscalatura è critica. Il 3D non ha nessuna delle due cose. Allora che cosa dovrebbe avere una dimostrazione?

Nessuno lo sa. Ma ecco che cosa stanno inseguendo i ricercatori:

• Trovare una nuova “manopola di controllo”. La vorticità è la manopola di controllo del 2D: resta limitata, e tutto il resto segue da questo singolo fatto. In 3D serve una quantità diversa, qualcosa che rimanga domabile qualunque cosa faccia il fluido e che sia abbastanza potente da costringere l'intera soluzione a restare liscia per sempre. Nessuno l'ha trovata. I ricercatori la cercano da decenni, e manca ancora.
• Sfruttare strutture nascoste. I fluidi sono incomprimibili. Non possono essere schiacciati. Questo vincolo limita ciò che lo stiramento dei vortici può fare, e potrebbero esserci schemi geometrici più profondi sepolti nelle equazioni che nessuno ha ancora sfruttato pienamente.
• Dimostrare che si rompe davvero. Forse le soluzioni 3D possono esplodere. Sarebbe altrettanto enorme. Bisognerebbe costruire una specifica condizione iniziale in cui lo stiramento dei vortici sopraffà la viscosità e spinge la soluzione all'infinito in tempo finito; per le equazioni di Euler più semplici (Navier-Stokes senza attrito) la formazione di singolarità è stata dimostrata in contesti correlati, ma il caso viscoso resta completamente aperto.

Per saperne di più su ciò che è stato tentato, si veda Sottoproblemi di Navier-Stokes.

Una dimostrazione della regolarità globale in 3D richiederebbe di colmare il divario di supercriticità. Concretamente, serve una stima a priori della forma $\|u(t)\|_X \leq C(\|u_0\|_Y, t)$ per qualche norma $X$ al livello della scala critica o al di sopra, dove $C$ resta finito per ogni $t$. La stima di energia nota $u \in L^\infty_t L^2_x \cap L^2_t \dot{H}^1_x$ è mezzo derivato sotto il livello critico e non basta.

Diversi programmi di ricerca mirano a questo divario:

• Decomposizione in profili e concentrazione-compattezza. Adattati dal successo delle equazioni dispersive critiche (Kenig-Merle 2006), questi metodi cercano di classificare i profili di blowup. Per Navier-Stokes esistono risultati parziali (per esempio Gallagher-Koch-Planchon 2016), ma la natura supercritica dell'energia rende il programma completo più difficile da realizzare rispetto ai contesti critici per l'energia delle equazioni delle onde o di Schrödinger.
• Estensioni di soluzioni mild. Il quadro di Fujita-Kato (1964) dà buona positura locale in $\dot{H}^{1/2}(\mathbb{R}^3)$ e buona positura globale per dati piccoli in spazi critici ($L^3$, $\dot{H}^{1/2}$, $BMO^{-1}$). La domanda è se le soluzioni con dati grandi possano essere proseguite globalmente, il che richiede di controllare la norma critica.
• Criteri di regolarità. Oltre Beale-Kato-Majda ($\int_0^T \|\omega\|_{\infty} < \infty$), there are Prodi-Serrin conditions ($u \in L^p_t L^q_x$ with $2/p + 3/q = 1$, $q > 3$), Escauriaza-Seregin-Šverák ($u \in L^\infty_t L^3_x$, 2003), e altri criteri di endpoint. Ciascuno riduce la regolarità globale a una singola stima a priori, ma dimostrare quella stima resta aperto.
• Costruire blowup. Tao (2016) ha costruito una soluzione con blowup per un sistema di Navier-Stokes mediato che rispetta l'identità di energia e la riscalatura, ma non la piena struttura a divergenza nulla. Questo ci dice che qualsiasi dimostrazione di regolarità deve usare la specifica struttura geometrica della non linearità, non solo le sue proprietà di scala. Resta aperto se le vere equazioni di Navier-Stokes ammettano blowup.

Per il problema inviscido, il blowup $C^{1,\alpha}$ di Elgindi per Euler 3D (2021) mostra che lo stiramento dei vortici può produrre singolarità sotto la regolarità $C^\infty$. La questione del blowup di Euler liscio ($C^\infty$) resta aperta, così come la domanda se la viscosità possa arrestare tali meccanismi nel contesto di Navier-Stokes.

Tutto quanto sopra, in una tabella:

Non è una questione tecnica. Il divario tra 2D e 3D è un abisso. La strategia dimostrativa che funziona perfettamente in due dimensioni non “ha solo bisogno di un po' più di lavoro” per gestire tre dimensioni; fondamentalmente non può funzionare perché la struttura matematica da cui dipende, il principio del massimo per la vorticità e la criticità energetica che rendono il 2D così trattabile, semplicemente non esiste in 3D.

Per le equazioni complete, si veda Che cosa sono le equazioni di Navier-Stokes? Per l'enunciato preciso del problema aperto, si veda Esistenza e regolarità per Navier-Stokes. Per capire perché è così difficile, si veda Perché Navier-Stokes è difficile.

I seguenti contrasti riassumono la divisione matematica:

Il risultato 2D non è semplicemente un riscaldamento in dimensione inferiore. È un teorema completo il cui meccanismo di dimostrazione (il principio del massimo per la vorticità combinato con la criticità dell'energia) non ha alcuna controparte nota in 3D. Qualsiasi risoluzione del problema 3D, sia essa regolarità o blowup, richiederà idee fondamentalmente nuove. Per lo stato attuale dei risultati parziali e dei programmi di ricerca, si veda Sottoproblemi di Navier-Stokes e Perché Navier-Stokes è difficile.