Approcci al problema di Navier-Stokes

L'approccio più antico parte dall'energia. Un fluido in movimento trasporta energia cinetica, e la viscosità la dissipa, come l'attrito che porta le cose a fermarsi. L'energia totale può solo diminuire nel tempo, supponendo che nulla immetta energia dall'esterno.

Leray lo capì nel 1934 e fece una mossa decisiva: usare il vincolo energetico per dimostrare che deve esistere una soluzione debole globale con energia cinetica finita. Si costruiscono soluzioni approssimate, artificialmente regolarizzate. Si mostra che soddisfano tutte il vincolo energetico. Si passa al limite. Qualcosa deve sopravvivere in quel limite, e così accade.

Ma c'è un problema. I vincoli energetici sono strumenti grossolani. Garantiscono che il fluido abbia energia totale finita, certo, ma non possono dirti che la velocità resta finita in ogni singolo punto dello spazio e del tempo. Questo divario tra "energia finita" e "regolare ovunque" è esattamente il problema della regolarità, ed è aperto da novant'anni.

Link agli articoli: Leray, Sur le mouvement d'un liquide visqueux emplissant l'espace (1934); Hopf, Über die Anfangswertaufgabe für die hydrodynamischen Grundgleichungen (1951).

La costruzione di Leray-Hopf procede in due varianti principali: approssimazione di Galerkin (proiezione su sottospazi finito-dimensionali) oppure mollificazione di Friedrichs (regolarizzazione della non linearità tramite convoluzione). Entrambe condividono lo stesso schema in quattro passaggi. La strategia standard ha quattro passi:

1. Approssimazione: Risolvere il sistema mollificato $\partial_t u_\varepsilon + (J_\varepsilon u_\varepsilon \cdot \nabla) u_\varepsilon = \nu \Delta u_\varepsilon - \nabla p_\varepsilon$ su sottospazi finito-dimensionali.
2. Stima energetica: La stima a priori $\|u_\varepsilon(t)\|_{L^2}^2 + 2\nu \int_0^t \|\nabla u_\varepsilon\|_{L^2}^2 \leq \|u_0\|_{L^2}^2$ vale uniformemente in $\varepsilon$.
3. Compattezza: Estrarre una sottosuccessione debolmente convergente $u_\varepsilon \rightharpoonup u$ in $L^2_t \dot{H}^1_x$ usando il lemma di Aubin-Lions.
4. Passaggio al limite: Il termine non lineare converge grazie alla convergenza forte in $L^2_{\text{loc}}$ di $u_\varepsilon$.

La soluzione debole risultante soddisfa la disuguaglianza energetica, non l'uguaglianza. L'energia può essere persa in tempi irregolari. Ed è questo l'intero problema: il divario tra la classe energetica $L^\infty_t L^2_x \cap L^2_t \dot{H}^1_x$ e la reale regolarità è precisamente ciò che non riusciamo a colmare.

Link agli articoli: Leray, Sur le mouvement d'un liquide visqueux emplissant l'espace (1934); Hopf, Über die Anfangswertaufgabe für die hydrodynamischen Grundgleichungen (1951).

L'approccio di Caffarelli-Kohn-Nirenberg (1982) non tenta di dimostrare la regolarità completa. Pone una domanda del tutto diversa: quanto possono essere davvero gravi le singolarità?

Quasi per nulla. Il loro teorema di $\varepsilon$-regolarità dice che, se certe quantità locali invarianti per scala sono abbastanza piccole in una piccola regione spazio-temporale, la soluzione è automaticamente regolare lì. E poiché l'energia totale è finita, semplicemente non c'è abbastanza "budget" perché molti punti singolari coesistano.

Pensala così. Un muro può avere crepe. Ma la lunghezza totale di tutte quelle crepe messe insieme è zero, il che significa che l'insieme singolare è estremamente piccolo nel senso parabolico della teoria della misura (misura di Hausdorff parabolica unidimensionale nulla).

Link agli articoli: Caffarelli-Kohn-Nirenberg, Partial regularity of suitable weak solutions of the Navier-Stokes equations (1982); Albritton-Barker-Prange, Epsilon regularity for the Navier-Stokes equations via weak-strong uniqueness.

Il teorema di $\varepsilon$-regolarità di CKN: per soluzioni deboli adatte, la piccolezza di certe quantità locali invarianti per scala (che coinvolgono sia $u$ sia $p$ su un cilindro parabolico $Q_r(z_0) = B_r(x_0) \times (t_0 - r^2, t_0)$) implica la regolarità (continuità Hölderiana) vicino al punto $z_0 = (x_0, t_0)$.

La dimostrazione combina la disuguaglianza energetica locale con un'iterazione di tipo Campanato: se l'energia invariante per scala è piccola, un argomento di bootstrap mostra che $u$ è limitata, poi Hölderiana, quindi liscia per la teoria classica di Schauder.

La stima dimensionale $\mathcal{P}^1(\Sigma) = 0$ segue da un argomento di ricoprimento di Vitali. I punti singolari impongono limiti inferiori su quantità localizzate di dissipazione/energia invarianti per scala; il ricoprimento controlla poi la misura parabolica unidimensionale dell'insieme singolare tramite una misura locale dell'energia, che si annulla alle piccole scale.

Link agli articoli: Caffarelli-Kohn-Nirenberg, Partial regularity of suitable weak solutions of the Navier-Stokes equations (1982); Albritton-Barker-Prange, Epsilon regularity for the Navier-Stokes equations via weak-strong uniqueness.

Ecco una riduzione netta dell'intero problema. Beale, Kato e Majda dimostrarono nel 1984 che, per le equazioni di Euler 3D, il blow-up può avvenire solo se si perde il controllo della vorticità. Criteri analoghi furono poi stabiliti per Navier-Stokes. Tutto qui. Una sola condizione.

La vorticità misura la rotazione locale. Il criterio BKM dice: mantieni limitata la massima rotazione nella norma giusta, e la soluzione resta regolare. Tutto il resto segue automaticamente.

Una famiglia di quantità da controllare. Sfortunatamente, controllarle davvero si è rivelato esattamente difficile quanto il problema originale. La riduzione è pulita. L'esecuzione resta fuori portata.

Link agli articoli: Beale-Kato-Majda, Remarks on the breakdown of smooth solutions for the 3-D Euler equations (1984); Kozono-Taniuchi, Bilinear estimates in BMO and the Navier-Stokes equations (2000).

Il teorema originale di Beale-Kato-Majda (1984) riguarda le equazioni di Euler 3D. Per Navier-Stokes, criteri di continuazione analoghi implicano che una soluzione liscia $u$ su $[0, T^*)$ si estende oltre $T^*$ ogni volta che

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

dove $\omega = \nabla \times u$ è la vorticità. I raffinamenti includono:

• Kozono-Taniuchi (2000): $\|\omega\|_{L^\infty}$ può essere sostituito da $\|\omega\|_{\mathrm{BMO}}$
• Varianti in spazi di Besov: anche un controllo critico o borderline in spazi di Besov può servire come criterio di continuazione
• Criteri con restrizione direzionale: condizioni di tipo Serrin su componenti di $\nabla u$ possono servire anch'esse come criteri di continuazione (si veda ad es. Beirão da Veiga, 1995)

Questi criteri si collegano al quadro dello stiramento dei vortici: qualunque singolarità in tempo finito deve costringere la vorticità ad accumularsi troppo rapidamente perché l'integrale temporale sopra resti finito.

Link agli articoli: Beale-Kato-Majda, Remarks on the breakdown of smooth solutions for the 3-D Euler equations (1984); Kozono-Taniuchi, Bilinear estimates in BMO and the Navier-Stokes equations (2000); Chemin-Planchon, Self-improving bounds for the Navier-Stokes equations (2012).

Un punto di vista più moderno lavora con spazi funzionali (come $L^3$ o $\dot{H}^{1/2}$) che si collocano esattamente al confine di ciò che la simmetria di scala consente. Questi sono spazi critici, ed è lì che si trovano i risultati di regolarità più fini.

La logica è chiara: se si riesce a mostrare che una soluzione rimane entro certi limiti in spazi critici, la regolarità segue automaticamente. Diversi gruppi hanno dimostrato risultati di questo tipo, costruendo un intero repertorio di criteri di regolarità (condizioni che garantiscono la regolarità se si riesce a verificarle).

Il problema è il divario. I metodi energetici danno un controllo supercritico rispetto all'energia, quindi al di sotto della scala critica. Ci servono stime critiche. Quel divario è stretto, talvolta appena una derivata di regolarità, ma ha resistito a ogni tentativo di colmarlo.

Link agli articoli: Koch-Tataru, Well-posedness for the Navier-Stokes equations (2001); Kenig-Koch, An alternative approach to regularity for the Navier-Stokes equations in critical spaces; Gallagher-Koch-Planchon, A profile decomposition approach to the $L^\infty_t(L^3_x)$ Navier-Stokes regularity criterion. Per un confronto dettagliato sul perché la criticità energetica funziona in 2D ma fallisce in 3D, si veda Perché il 2D è più facile del 3D.

Principali programmi sulla regolarità in spazi critici:

• Koch-Tataru (2001): Buona positura globale per dati piccoli in $\text{BMO}^{-1}$. La stima chiave è un vincolo di punto fisso per la mappa bilineare di Duhamel nello spazio di soluzioni di Koch-Tataru costruito sul flusso di calore con dati in $\text{BMO}^{-1}$, non una semplice stima statica di prodotto in $\text{BMO}^{-1}$. È tra i risultati più fini noti per i metodi perturbativi.
• Gallagher-Koch-Planchon (2013): Approccio tramite decomposizione in profili al criterio di regolarità $L^\infty_t L^3_x$ per Navier-Stokes. Ogni successione di soluzioni con norma critica limitata ammette una sottosuccessione che si decompone in profili asintoticamente disaccoppiati.

L’ostruzione fondamentale: non si conosce alcun funzionale coercivo che sia al tempo stesso controllato dall’evoluzione e critico rispetto alla scala di Navier-Stokes.

Link agli articoli: Koch-Tataru, Well-posedness for the Navier-Stokes equations (2001); Kenig-Koch, An alternative approach to regularity for the Navier-Stokes equations in critical spaces; Gallagher-Koch-Planchon, A profile decomposition approach to the $L^\infty_t(L^3_x)$ Navier-Stokes regularity criterion.

La teoria moderna delle PDE attinge ampiamente dall’analisi armonica. L’idea di fondo: scomporre una funzione in onde a frequenze diverse, come si separerebbe un accordo musicale nelle sue singole note. Solo che qui le “note” sono oscillazioni spaziali della velocità del fluido a scale enormemente diverse.

La decomposizione di Littlewood-Paley fa esattamente questo. Scompone il campo di velocità in componenti scala per scala. Traccia come l’energia fluisce tra di esse. All’improvviso l’intuizione fisica informale della “cascata di energia” diventa qualcosa su cui si possono effettivamente dimostrare teoremi, e teoremi precisi. Questi metodi hanno prodotto molti dei risultati più fini sui criteri di regolarità e sui tassi di blowup.

Link agli articoli: Cannone-Meyer, Littlewood-Paley decomposition and Navier-Stokes equations (1995); Gallagher-Koch-Planchon, A profile decomposition approach to the $L^\infty_t(L^3_x)$ Navier-Stokes regularity criterion.

La teoria di Littlewood-Paley decompone $u = \sum_j \Delta_j u$, dove $\Delta_j$ localizza alle frequenze $|\xi| \sim 2^j$. Applicata alla non linearità di trasporto $(u \cdot \nabla)u$:

La decomposizione in paraprodotti della non linearità $(u \cdot \nabla)u$ si scinde in interazioni di frequenza basso-alto, alto-basso e alto-alto:

$$(u \cdot \nabla)u = T_u \nabla u + T_{\nabla u} u + R(u, \nabla u)$$

dove $T$ è il paraprodotto e $R$ il resto. Ciascun termine ha proprietà di regolarità diverse negli spazi di Besov $\dot{B}^s_{p,q}$.

Risultati chiave che usano questa macchina:

• Spazi di Chemin-Lerner: $\widetilde{L}^\rho_T \dot{B}^s_{p,q}$ forniscono il quadro naturale per la buona positura critica: la forma bilineare di Navier-Stokes manda $\widetilde{L}^\infty_T \dot{B}^{-1+3/p}_{p,q} \times \widetilde{L}^1_T \dot{B}^{1+3/p}_{p,q} \to \widetilde{L}^1_T \dot{B}^{-1+3/p}_{p,q}$.
• Cannone-Meyer: I metodi di Littlewood-Paley danno una formulazione wavelet/Besov pulita della teoria per dati piccoli.

Link agli articoli: Cannone-Meyer, Littlewood-Paley decomposition and Navier-Stokes equations (1995); Gallagher-Koch-Planchon, A profile decomposition approach to the $L^\infty_t(L^3_x)$ Navier-Stokes regularity criterion.

Ecco un istinto completamente diverso. Invece di seguire numeri (norme, energie), questi metodi studiano la forma della soluzione: come si piegano i tubi vorticosi, come le regioni di intensa rotazione si dispongono nello spazio.

L’intuizione chiave è che il blowup non riguarda soltanto qualcosa che diventa grande. Riguarda il fluido che si organizza in una configurazione geometrica molto specifica. Se si può mostrare che quella configurazione è impossibile (perché contraddice la struttura energia-dissipazione, o l’incomprimibilità, o entrambe), si è escluso il blowup senza mai calcolare una norma.

Questo punto di vista geometrico è diventato una prospettiva che ha ispirato diversi criteri rigorosi di regolarità accanto ai metodi puramente analitici. E ha un sapore diverso. Chiede che forma prende il disastro? invece di quanto può diventare grande questo numero?

Link agli articoli: Constantin-Fefferman, Geometric constraints on potentially singular solutions for the 3-D Euler equations (1993); Albritton-Barker-Prange, Localized smoothing and concentration for the Navier-Stokes equations in the half space.

Gli approcci geometrico-topologici sfruttano vincoli strutturali invisibili ai metodi puramente analitici:

• Geometria delle linee vorticose: Constantin e Fefferman (1993) hanno mostrato che, se il campo delle direzioni della vorticità $\hat{\omega} = \omega/|\omega|$ è Lipschitz nelle regioni di alta vorticità, la soluzione è regolare. Il blowup richiede che la direzione della vorticità sviluppi una singolarità simultaneamente alla sua intensità.
• Argomenti di incompatibilità: se una configurazione di blowup è geometricamente vincolata (ad esempio tramite limiti di impacchettamento sul numero di regioni di concentrazione indipendenti che possono rientrare nei bilanci di energia e dissipazione), si può derivare una contraddizione senza stimare direttamente norme critiche.
• Partizione in casi (speculativa/programmatica): una strategia proposta classificherebbe ogni regione spaziale come appartenente a uno fra un numero finito di scenari (ad esempio localmente regolare, di tipo I, di tipo II, densamente impacchettata) e tenterebbe di mostrare che ciascuno scenario o dà regolarità oppure trasferisce il problema a un argomento di conteggio limitato. Questo rimane un programma di ricerca più che un risultato consolidato.

La pagina sulla dimostrazione di questo sito esplora argomenti di questo tipo; nulla qui è presentato come una dimostrazione formale completa del problema del Millennio.

Link agli articoli: Constantin-Fefferman, Geometric constraints on potentially singular solutions for the 3-D Euler equations (1993); Albritton-Barker-Prange, Localized smoothing and concentration for the Navier-Stokes equations in the half space.

Questo ha colto molti di sorpresa. Le soluzioni deboli ottenute con il metodo di Leray (Sezione 1) risultano essere non uniche, almeno quando è presente una forzante esterna.

L’arma è l’integrazione convessa, una tecnica originariamente sviluppata per problemi di geometria e adattata alle equazioni dei fluidi da De Lellis e Székelyhidi a partire circa dal 2009. L’idea: costruire soluzioni “selvagge” aggiungendo iterativamente correzioni ad alta frequenza che, nel loro insieme, soddisfano l’equazione ma si comportano in modo irregolare.

Per le Navier-Stokes 3D, Buckmaster e Vicol (2019) hanno dimostrato la non unicità di soluzioni deboli a energia finita al di sotto della classe di Leray-Hopf. La non unicità per Euler nasce invece dal programma di integrazione convessa di De Lellis-Székelyhidi, Isett e Buckmaster. Poi, nel 2022, Albritton, Brué e Colombo hanno dimostrato che persino le soluzioni di Leray-Hopf delle Navier-Stokes 3D non sono uniche quando è presente una forza esterna. Resta aperto se la non unicità persista per le equazioni di Navier-Stokes non forzate.

Perché è importante? Perché “esiste una soluzione debole” è stato il risultato principale fin dal 1934. Ora sappiamo che non determina una risposta unica. La domanda diventa più precisa: quale soluzione, se ce n’è una, è quella fisicamente corretta?

Link agli articoli: De Lellis-Székelyhidi, Dissipative continuous Euler flows (2013); Buckmaster-Vicol, Nonuniqueness of weak solutions to the Navier-Stokes equation (2019); Albritton-Brué-Colombo, Non-uniqueness of Leray solutions of the forced Navier-Stokes equations (2022).

L’integrazione convessa per le equazioni dei fluidi nasce nel programma di De Lellis-Székelyhidi (2009–2013), adattando la tecnica di immersione isometrica Nash-Kuiper $C^1$ per costruire soluzioni deboli delle equazioni di Euler che dissipano energia. Tappe chiave:

• De Lellis-Székelyhidi (2013): esistenza di flussi di Euler dissipativi continui ($C^0$) su $\mathbb{T}^3$ (regolarità di Hölder sotto $1/5$ ottenuta nel successivo lavoro Buckmaster-De Lellis-Isett-Székelyhidi del 2015; in seguito migliorata a $< 1/3$ da Isett, nel 2018, risolvendo il lato flessibile della congettura di Onsager).
• Buckmaster-Vicol (2019): non unicità di soluzioni deboli a energia finita per Navier-Stokes 3D, con soluzioni in $C^0_t H^\beta_x$ per qualche $\beta > 0$. La costruzione usa flussi di Beltrami intermittenti come mattoni di base, aggiungendo correzioni oscillanti a ogni passo iterativo mantenendo il controllo dello stress di Reynolds. Questo è al di sotto della classe energetica di Leray-Hopf, quindi non contraddice direttamente l’unicità di Leray.
• Albritton-Brué-Colombo (2022): non unicità delle soluzioni di Leray-Hopf per le equazioni di Navier-Stokes 3D forzate. La dimostrazione costruisce una soluzione di fondo autosimilare instabile e usa un meccanismo di instabilità per diramarsi in soluzioni di Leray-Hopf distinte a partire dagli stessi dati iniziali. Ciò mostra che la sola disuguaglianza di energia non seleziona una soluzione unica quando è presente una forzante.

La questione aperta centrale è se la non unicità persista per le equazioni di Navier-Stokes non forzate nella classe di Leray-Hopf. Il risultato forzato mostra che la disuguaglianza di energia non è un principio di selezione sufficiente, ma non chiarisce se l’equazione non forzata abbia una struttura aggiuntiva che ripristini l’unicità.

Link agli articoli: De Lellis-Székelyhidi, Dissipative continuous Euler flows (2013); Buckmaster-Vicol, Nonuniqueness of weak solutions to the Navier-Stokes equation (2019); Albritton-Brué-Colombo, Non-uniqueness of Leray solutions of the forced Navier-Stokes equations (2022).

Possiamo almeno escludere certe strategie di dimostrazione? Terence Tao ha mostrato nel 2016 che sì, possiamo. E il risultato è istruttivo.

Tao ha costruito una versione modificata delle equazioni di Navier-Stokes, un sistema “mediato”, che conserva molte caratteristiche strutturali chiave delle equazioni reali: l’identità di energia, il modo in cui cresce l’enstrofia (una misura dell’intensità della vorticità), la simmetria di scala. Ma in questo sistema modificato, le soluzioni esplodono in tempo finito.

L’implicazione esclude ampie famiglie di strategie di dimostrazione. Qualunque dimostrazione della regolarità globale per le equazioni reali deve usare qualche proprietà strutturale specifica della vera non linearità che il sistema mediato non possiede. Non si può dimostrare la regolarità usando soltanto stime di energia, scaling e crescita dell’enstrofia. Questi strumenti, da soli, sono compatibili con il blow-up.

Questo non dice che le equazioni reali esplodano. Dice che intere famiglie di strategie dimostrative sono vicoli ciechi. La dimostrazione finale (se la regolarità vale) dovrà essere più precisa di un argomento energetico generico. Molto più precisa.

Link agli articoli: Tao, Finite time blowup for an averaged three-dimensional Navier-Stokes equation (2016).

Tao (2016) considera un sistema della forma

$$\partial_t u + \tilde{B}(u, u) = \nu \Delta u - \nabla p, \quad \nabla \cdot u = 0,$$

dove $\tilde{B}$ è un operatore bilineare che coincide con la vera non linearità di Navier-Stokes $(u \cdot \nabla)u$ nei seguenti sensi:

• è un moltiplicatore di Fourier di ordine 1, che preserva lo scaling $u(x,t) \mapsto \lambda u(\lambda x, \lambda^2 t)$
• soddisfa la stessa identità di energia: $\langle \tilde{B}(u,u), u \rangle = 0$
• riproduce la struttura di crescita dell’enstrofia

Per questo sistema mediato, Tao costruisce dati iniziali lisci la cui soluzione esplode in tempo finito. Il meccanismo di blow-up programma una sequenza di concentrazioni di enstrofia sempre più acute a scale sempre più piccole, con ogni stadio che raddoppia l’enstrofia in una cascata controllata.

L’ostruzione che ne deriva: qualunque quantità supercritica che sia (a) controllata dall’evoluzione nella classe energetica e (b) invariante sotto lo scaling di Navier-Stokes non può, da sola, escludere il blow-up, perché sarebbe controllata anche nel sistema mediato, che invece esplode. Una dimostrazione di regolarità deve sfruttare la specifica struttura di cancellazione algebrica del vero termine di trasporto $(u \cdot \nabla)u$ che l’operatore mediato $\tilde{B}$ non condivide.

Link agli articoli: Tao, Finite time blowup for an averaged three-dimensional Navier-Stokes equation (2016).

Questo articolo fa parte di Progressi.

Dalla dimostrazione di esistenza di Leray del 1934 fino all’integrazione convessa e alle barriere dimostrative di Tao, queste sono le principali strategie che sono state tentate sul problema di Navier-Stokes 3D. Nessuna ha risolto il problema completo della regolarità 3D. Per il contesto su come la viscosità modelli la matematica rispetto alle equazioni inviscide di Euler, vedi Euler vs. Navier-Stokes. Per lo stato attuale, vedi Il problema di Navier-Stokes è risolto? Per l’enunciato formale esatto, torna a Il problema del Millennio.

Questo articolo fa parte di Progressi.

Gli approcci sopra descritti rappresentano i principali filoni rigorosi nella letteratura sulla regolarità e l'unicità fino al 2022. Nessuna combinazione ha risolto il problema completo in 3D. Il campo continua a evolversi. I metodi di dimostrazione assistita dal computer sono un'area attiva trattata separatamente in questo sito.

Per il confronto tra sistemi viscosi e inviscidi (e sul perché la viscosità aiuti ma non abbastanza), si veda Euler vs. Navier-Stokes. Per i sottoproblemi a cui mirano questi approcci, si veda Sottoproblemi. Per gli ostacoli legati allo scaling, si veda Perché è difficile.