Perché il problema di Navier-Stokes è difficile

Molte delle equazioni che si incontrano per prime in fisica sono lineari: raddoppi l'input e la risposta raddoppia. Navier-Stokes non è così.

Navier-Stokes? Non lineare. La velocità del fluido influenza il proprio tasso di variazione, il che significa che il fluido spinge se stesso. Immagina di cercare di prevedere dove andrà una folla quando il movimento di ogni singola persona dipende da ciò che fanno tutti quelli intorno a lei, e ciò che fanno quelle persone dipende da tutti quelli intorno a loro, in una spirale che si propaga all'infinito. Questa è la situazione che hai davanti.

Il colpevole è il termine di auto-interazione $(u \cdot \nabla)u$. Crea cicli di retroazione in cui piccole perturbazioni si amplificano fino a diventare grandi, ed è il motivo per cui la turbolenza dei fluidi è così incredibilmente complessa (vedi sottoproblemi per saperne di più).

La non linearità convettiva $(u \cdot \nabla)u$ è l'ostacolo fondamentale. Nella formulazione in vorticità $\omega = \nabla \times u$, l'equazione diventa

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

Il termine di stiramento dei vortici $(\omega \cdot \nabla)u$ non ha segno. Può amplificare la vorticità senza limite. In 2D questo termine scompare (poiché $\omega$ è uno scalare perpendicolare al flusso), ed è per questo che la regolarità globale in 2D è nota (Ladyzhenskaya, 1969). In 3D, lo stiramento dei vortici è il principale meccanismo candidato per il blowup in tempo finito.

Ecco il punto cruciale: la non linearità è quadratica in $u$. La stima dell'energia $H^1$ dà $\|\nabla u\|_{L^2}$, ma controllare $(u \cdot \nabla)u$ in $L^2$ richiede tipicamente informazioni più forti come $u \in L^\infty$ insieme a $\nabla u \in L^2$, o un controllo equivalente alla scala critica; la sola classe energetica non lo fornisce.

Le equazioni di Navier-Stokes hanno una simmetria di scala. Ingrandisci una soluzione, rendi tutto più piccolo e più veloce nelle giuste proporzioni, e ottieni un'altra soluzione perfettamente valida. Questa simmetria è matematicamente naturale, ma analiticamente pericolosa.

Perché? L'unica quantità che possiamo controllare in modo affidabile è l'energia totale del fluido, e si trova alla scala completamente sbagliata: ci dice qualcosa sul quadro globale ma assolutamente nulla su ciò che accade alle scale microscopiche dove si formerebbe davvero un blowup.

Pensa a monitorare il consumo totale di elettricità di una città per rilevare una scintilla. Utile? Certo. Abbastanza fine? Neanche lontanamente. Quel divario è l'intero problema.

Sotto la scala naturale $u_\lambda(x,t) = \lambda u(\lambda x, \lambda^2 t)$, lo spazio di Sobolev critico è $\dot{H}^{1/2}(\mathbb{R}^3)$ (equivalentemente $L^3$). Una norma è critica per scala se resta invariata sotto questa trasformazione. Se una norma controllata diventa più piccola quando si fa zoom, allora l'equazione è supercritica rispetto a quel controllo: il vincolo si colloca sotto la scala necessaria per vedere la concentrazione.

La disuguaglianza dell'energia dà il controllo di $u \in L^\infty_t L^2_x \cap L^2_t \dot{H}^1_x$. Al livello dell'energia, il problema 3D è supercritico:

$$\|u_\lambda\|_{L^2} = \lambda^{-1/2} \|u\|_{L^2}, \quad \|u_\lambda\|_{L^2_t \dot{H}^1_x} = \lambda^{-1/2} \|u\|_{L^2_t \dot{H}^1_x}$$

Quando $\lambda \to \infty$ (zoom sulle piccole scale), questi vincoli si indeboliscono rispetto alle quantità invarianti $\|u\|_{L^3}$ e $\|u\|_{\dot{H}^{1/2}}$. I termini non lineare e viscoso scalano insieme; il problema è che la stima a priori nota si colloca sotto la scala critica. Passare dalla classe energetica a una norma critica è la difficoltà centrale.

Osserva un fiume. Il moto del fluido diventa caotico. Turbolento. Grandi vortici si frantumano in vortici più piccoli, che si frantumano in vortici ancora più piccoli, in una cascata che arriva fino alle scale microscopiche dove la viscosità finalmente leviga tutto.

Kolmogorov descrisse questa cascata di energia nel 1941, e le equazioni di Navier-Stokes la catturano magnificamente. Ma ecco ciò che toglie il sonno: e se l'energia si concentrasse in regioni sempre più piccole più velocemente di quanto la viscosità riesca a dissiparla? Questo è un blowup.

Può davvero accadere? Oppure la viscosità vince sempre? Questa è la questione aperta, punto. Per il ponte fisico dal numero di Reynolds a questo quadro delle piccole scale, vedi Numero di Reynolds, turbolenza e perché le piccole scale contano.

La teoria K41 di Kolmogorov predice uno spettro energetico $E(k) \sim \varepsilon^{2/3} k^{-5/3}$ nell'intervallo inerziale $k_f \ll k \ll k_\eta$, dove $k_\eta \sim (\varepsilon/\nu^3)^{1/4}$ è il numero d'onda dissipativo di Kolmogorov. Flusso di energia costante attraverso le scale. Questo è il quadro pulito.

La questione della regolarità chiede qualcosa di più inquietante: questa cascata può degenerare? Può la scala dissipativa $k_\eta^{-1}$ ridursi a zero in tempo finito, con $\|\nabla u\|_{L^2} \to \infty$ mentre l'energia totale resta finita?

La congettura dell'anomalia dissipativa di Onsager (1949) dice di sì, in un certo senso: nel limite di viscosità evanescente $\nu \to 0$, la dissipazione di energia persiste, e soluzioni deboli di Euler possono dissipare energia. Confermato per esponenti di Hölder inferiori a $1/3$ (Isett, 2018; Buckmaster et al., 2018). Ma che cosa questo significhi per la regolarità di Navier-Stokes resta davvero poco chiaro. Per un'intuizione a livello di regime, vedi Numero di Reynolds, turbolenza e perché le piccole scale contano.

La pressione nelle equazioni di Navier-Stokes è strana. Non è affatto una variabile indipendente; la velocità la determina completamente attraverso un unico vincolo: il fluido è incomprimibile, quindi non può essere compresso.

Questo rende la pressione non locale. Nel modello incomprimibile, cambiare la velocità in una regione influenza globalmente il campo di pressione, perché la pressione è determinata dall'intero campo di velocità in quell'istante.

Per l'analisi, questo è devastante. Non puoi studiare ciò che accade in un singolo punto senza tenere conto dell'intero fluido tutto insieme. Ragionamento locale? Scordatelo.

Il vincolo di incomprimibilità $\nabla \cdot u = 0$ determina la pressione tramite l'equazione di Poisson

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

quindi $p = (-\Delta)^{-1} \partial_i \partial_j (u_i u_j)$, coinvolgendo le trasformate di Riesz (operatori integrali singolari). La pressione è una funzione non locale della velocità, e questa non località è l'ostacolo centrale alle stime puntuali o locali nello spazio.

Qui gli argomenti standard basati sul principio del massimo falliscono: anche se il termine viscoso $\nu \Delta u$ è dissipativo, il gradiente di pressione $-\nabla p$ può concentrare energia da regioni lontane. La teoria di Caffarelli-Kohn-Nirenberg gestisce questo tramite disuguaglianze locali dell'energia su cilindri parabolici, ma estrarre da esse la regolarità puntuale resta il passaggio difficile.

Per le equazioni di Navier-Stokes incomprimibili in 2D negli ambienti standard, sono note soluzioni lisce globali; ciò è stato stabilito in lavori classici, incluso quello di Ladyzhenskaya (1969). Vedi perché il 2D è più facile.

Tre dimensioni? Tutto crolla, e la ragione si riduce a un solo meccanismo: stiramento dei vortici. In 2D, i vortici possono ruotare e fondersi, ma non possono stirarsi. In 3D, il fluido può afferrare tubi vorticosi e tirarli rendendoli sempre più sottili, potenzialmente concentrando ogni ultimo frammento di energia in un filamento infinitamente sottile.

Questo stiramento può crescere senza controllo fino all'infinito in tempo finito, oppure la viscosità interviene sempre? Questa è la domanda da un milione di dollari. Letteralmente.

La dicotomia è netta.

2D: La vorticità $\omega$ è uno scalare che soddisfa $\partial_t \omega + u \cdot \nabla \omega = \nu \Delta \omega$. Il principio del massimo dà $\|\omega(t)\|_{L^\infty} \leq \|\omega_0\|_{L^\infty}$, BKM implica la regolarità globale, e il termine di stiramento dei vortici $(\omega \cdot \nabla)u$ è identicamente nullo.

3D: La vorticità $\omega \in \mathbb{R}^3$ soddisfa $\partial_t \omega + (u \cdot \nabla)\omega = (\omega \cdot \nabla)u + \nu \Delta \omega$, dove il termine di stiramento $(\omega \cdot \nabla)u$ può amplificare $|\omega|$ in modo superlineare (formalmente $\sim |\omega|^2$ tramite Biot-Savart). Nessun principio del massimo disponibile.

L'enstrofia $\|\omega\|_{L^2}^2$ soddisfa

$$\frac{d}{dt}\|\omega\|_{L^2}^2 \leq C\|\omega\|_{L^2}^2 \|\nabla u\|_{L^\infty} - 2\nu \|\nabla \omega\|_{L^2}^2$$

Controllare $\|\nabla u\|_{L^\infty}$ richiede $\omega \in L^\infty$, che richiede di controllare $\|\nabla u\|_{L^\infty}$. Circolare. Nessuna tecnica esistente ha spezzato questo ciclo.

Questo articolo fa parte di Il problema.

Questi ostacoli hanno portato i matematici a scomporre il problema in sottoproblemi e a sviluppare approcci specializzati per ciascuno.

Per un contesto sul perché il termine viscoso aiuti ma non sia sufficiente, si veda Euler vs. Navier-Stokes.

Questo articolo fa parte di Il problema.

Il problema della regolarità si scompone in parti trattabili; si veda Sottoproblemi. Gli strumenti analitici costruiti per affrontare questi ostacoli sono trattati in Approcci. Perché il termine viscoso aiuta ma ancora non colma il divario? Euler vs. Navier-Stokes.