Euler vs. Navier-Stokes: qual è la differenza?

Le equazioni di Euler descrivono un fluido con attrito interno nullo. Nessuna viscosità. Le equazioni di Navier-Stokes descrivono lo stesso fluido con la viscosità inclusa.

Matematicamente, tutta la differenza sta in un solo termine: $\nu \Delta u$, il termine di diffusione viscosa. Rimuovilo e Navier-Stokes diventa Euler. Mantienilo e l’equazione acquista un meccanismo di regolarizzazione che cambia sia la fisica sia l’analisi in modi che non ti aspetteresti da un singolo termine aggiuntivo.

Quel termine è il motivo per cui il fumo si dissipa, per cui si formano strati limite lungo le superfici, e per cui il Problema del Millennio di Navier-Stokes ha un carattere completamente diverso dalla corrispondente questione per Euler.

Le equazioni di Euler incomprimibili su $\mathbb{R}^3$ sono

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

Le equazioni di Navier-Stokes incomprimibili sono

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

con $\nu > 0$ la viscosità cinematica. Formalmente, ponendo $\nu = 0$ in Navier-Stokes si recupera Euler. Ma questa sostituzione formale nasconde il fatto che il limite $\nu \to 0$ è singolare: il termine viscoso $\nu \Delta u$ contiene la derivata spaziale di ordine più alto nel sistema, e eliminarlo cambia il tipo di PDE e gli spazi funzionali in cui vivono le soluzioni.

Entrambe le equazioni nella loro forma incomprimibile standard, scritte in modo che il confronto sia evidente:

Euler:

$$\partial_t u + (u \cdot \nabla)u = -\nabla p$$

Navier-Stokes:

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

Stesso lato sinistro: il tasso di variazione della velocità più il termine non lineare di autotrasporto $(u \cdot \nabla)u$. Entrambe impongono l’incomprimibilità tramite $\nabla \cdot u = 0$. L’unica differenza strutturale è il termine viscoso $\nu \Delta u$ sul lato destro di Navier-Stokes.

Il parametro $\nu$ è la viscosità cinematica, una costante fisica del fluido. Il miele ha un $\nu$ grande. L’aria ne ha uno piccolo. Le equazioni di Euler corrispondono a $\nu = 0$: un’idealizzazione perfettamente priva di attrito che può approssimare alcuni flussi ad alto numero di Reynolds lontano dai bordi, ma che non esiste in alcun fluido reale.

Entrambi i sistemi condividono la forma bilineare $(u \cdot \nabla)u$ e la pressione determinata implicitamente dal vincolo di divergenza nulla. Prendendo la divergenza dell’equazione della quantità di moto e usando $\nabla \cdot u = 0$ si ottiene l’equazione di Poisson per la pressione

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

che è identica per entrambi i sistemi. La pressione è un funzionale non locale di $u$ in entrambi i casi.

Il termine viscoso $\nu \Delta u$ è un operatore ellittico lineare del secondo ordine che agisce su ciascuna componente della velocità. È regolarizzazione parabolica: la sua presenza rende Navier-Stokes un sistema parabolico semilineare, mentre Euler incomprimibile è un’equazione di trasporto non lineare del primo ordine con accoppiamento non locale della pressione. Questa differenza nel tipo di PDE determina quasi ogni successiva differenza nella teoria della regolarità.

La viscosità è l’attrito tra strati vicini di fluido. Uno strato veloce accanto a uno lento? La viscosità trasferisce quantità di moto tra i due, smussando la differenza di velocità. Concetto semplice. Le conseguenze sono enormi e separano Navier-Stokes da Euler in tre modi.

• Dissipazione. L’energia cinetica si trasforma in calore. Mescola il caffè, poi fermati. Alla fine torna a riposo perché la viscosità prosciuga il moto convertendolo in energia termica. Euler non può prevederlo affatto, poiché nelle equazioni non c’è alcun meccanismo che dreni energia cinetica in calore.
• Strati limite. I fluidi reali aderiscono alle superfici (la condizione di no-slip), creando sottili strati di rapido cambiamento della velocità vicino alle pareti. Questi generano resistenza sulle ali degli aerei, perdite per attrito nei tubi e l’innesco della turbolenza ad alte velocità. I flussi di Euler soddisfano invece una condizione di scorrimento, quindi mancano completamente la resistenza viscosa di parete.
• Regolarizzazione alle piccole scale. La viscosità elimina i gradienti di velocità più marcati. Senza? Nulla impedisce al flusso di sviluppare strutture infinitamente fini, sempre più affilate. Questa regolarizzazione è esattamente ciò che rende la questione della regolarità per Navier-Stokes una bestia diversa da Euler.

L’identità dell’energia per Navier-Stokes su $\mathbb{R}^3$ (o su un dominio periodico) si legge

$$\frac{d}{dt}\frac{1}{2}\|u\|_{L^2}^2 = -\nu \|\nabla u\|_{L^2}^2,$$

quindi l’energia cinetica viene dissipata monotonicamente. Per Euler ($\nu = 0$), la norma $L^2$ di $u$ è formalmente conservata.

A livello del bordo, Navier-Stokes usa condizioni di no-slip ($u|_{\partial \Omega} = 0$), mentre Euler richiede solo l’impermeabilità ($u \cdot n = 0$). Quando $\nu \to 0$, il disallineamento tra queste condizioni genera lo strato limite di Prandtl. È un fenomeno di perturbazione singolare, e la gente ci lotta fin dall’articolo di Prandtl del 1904.

Fisicamente, la viscosità agisce come un filtro ad alta frequenza: smorza i modi di Fourier al tasso $\nu |k|^2$, sopprimendo preferenzialmente le piccole scale. Questo smorzamento spettrale è il meccanismo dietro la scala di dissipazione di Kolmogorov $\eta \sim (\nu^3 / \varepsilon)^{1/4}$ nella turbolenza. Vedi Numero di Reynolds e turbolenza per il quadro completo degli scalamenti.

Formalmente? Sì. Poni $\nu = 0$ e ottieni Euler. Ma è un pessimo punto in cui smettere di pensarci.

Il limite $\nu \to 0$ è singolare. La viscosità contiene le derivate di ordine più alto nell’equazione, quindi rimuoverla non è una piccola modifica. Cambia completamente il tipo di PDE con cui hai a che fare. Gli strati limite non si assottigliano in modo ordinato. Possono esplodere in turbolenza. Soluzioni perfettamente lisce sotto Navier-Stokes possono sviluppare comportamenti radicalmente diversi sotto Euler.

Sì, le due equazioni condividono il loro DNA matematico. Ma il limite a viscosità nulla è uno dei problemi aperti più profondi di tutta la dinamica dei fluidi, non un calcolo da tovagliolo.

Il limite inviscido $\nu \to 0$ è una perturbazione singolare: $\nu \Delta u$ contiene le derivate spaziali di ordine più alto nel sistema, quindi porre $\nu = 0$ abbassa l’ordine della PDE. Su domini con bordo, il limite è legato alla validità dell’espansione dello strato limite di Prandtl, che può fallire in modo spettacolare (Grenier 2000, Gérard-Varet & Dormy 2010).

Su $\mathbb{R}^3$ o $\mathbb{T}^3$ (senza bordi), le cose diventano più pulite. Se la soluzione di Euler $u^E$ resta liscia su $[0,T]$, allora le soluzioni di Navier-Stokes $u^\nu$ convergono a $u^E$ in $L^2$ quando $\nu \to 0$ (Kato 1972). Più pulite, ma non semplici: se le soluzioni di Euler restino persino lisce globalmente è esso stesso un problema aperto in 3D, ed è questo il punto.

Il limite si scontra anche frontalmente con la teoria della turbolenza. Il quadro di Kolmogorov richiede $\nu > 0$ per definire una scala di dissipazione. Eppure la dissipazione anomala, cioè la persistenza della dissipazione di energia quando $\nu \to 0$, è una congettura aperta da decenni. La congettura di Onsager (ora un teorema: Isett 2018, affinato da Buckmaster, De Lellis, Székelyhidi e Vicol nel 2019) caratterizza esattamente quando le soluzioni di Euler possono dissipare energia senza alcuna viscosità.

Ogni volta che la viscosità è trascurabile rispetto alle altre forze in gioco. Succede più spesso di quanto si pensi:

• Aerodinamica ad alta velocità lontano dalle superfici. Lontano da un’ala, il flusso d’aria è quasi inviscido. Gli ingegneri usano abitualmente solutori di Euler per il flusso principale e inseriscono correzioni di strato limite vicino alla parete.
• Flussi astrofisici. Nubi di gas interstellari, interni stellari, dischi di accrescimento attorno a buchi neri. A quelle scale, la viscosità molecolare è completamente irrilevante (anche se la viscosità effettiva turbolenta può non esserlo).
• Dinamica dei gas comprimibili. Onde d’urto. Detonazioni. Volo supersonico. La fisica dominante è pressione e inerzia, non attrito.
• Teoria pura. Euler merita di essere studiato di per sé, non solo come tappa verso Navier-Stokes. Si collega alla geometria riemanniana, alla dinamica dei vortici e a questioni profonde sulla struttura stessa della turbolenza.

Ma per qualunque cosa in cui contino attrito, resistenza o comportamento al bordo (flusso nei tubi, aerodinamica dei veicoli vicino alle superfici, circolazione sanguigna, tempo meteorologico su scale umane), serve Navier-Stokes. Punto.

Le equazioni di Euler governano il comportamento dominante ad alto numero di Reynolds $\mathrm{Re} = UL/\nu \gg 1$. A questi numeri di Reynolds, gli effetti viscosi sono confinati in sottili strati limite e strati di taglio interni, mentre il flusso principale lontano dalle pareti è ben approssimato da Euler.

Le equazioni di Euler comprimibili, un sistema iperbolico con velocità di propagazione finita, sono il modello standard per la dinamica dei gas, inclusa la formazione di onde d’urto e i problemi di Riemann. Differiscono dalle equazioni di Euler incomprimibili discusse sopra: Euler comprimibile è genuinamente iperbolico, mentre Euler incomprimibile ha accoppiamento non locale della pressione e velocità di propagazione infinita.

Nell’analisi matematica, Euler serve sia come oggetto limite per il problema della viscosità evanescente sia come ricco sistema di PDE con una propria teoria della regolarità, quantità conservate (elicità, Casimir tramite il quadro di Euler-Arnold sul gruppo dei diffeomorfismi) e connessioni con la meccanica geometrica.

È qui che il divario conta di più, ed è qui che le cose diventano davvero interessanti.

Il Problema del Millennio di Navier-Stokes pone una domanda che suona quasi troppo semplice: se si parte da un flusso liscio e ben comportato in tre dimensioni, la soluzione resta liscia per sempre, oppure può esplodere? Nessuno sulla Terra conosce la risposta.

La stessa domanda per Euler è aperta anche in 3D. Ma i due problemi hanno un sapore completamente diverso:

• Navier-Stokes ha la viscosità dalla sua parte. Sempre regolarizzante, sempre dissipativa dell’energia, sempre capace di smorzare i gradienti più ripidi. La vera domanda è se questa regolarizzazione sia abbastanza forte da sopraffare il termine non lineare prima che crei una singolarità.
• Euler non ha nulla. Zero regolarizzazione. Zero dissipazione. Il termine non lineare può amplificare i gradienti di velocità senza alcuna forza contraria, e capire se questo produca davvero una singolarità in tempo finito a partire da dati iniziali lisci in 3D è una delle maggiori questioni aperte nella teoria delle PDE.

In 2D, entrambe le equazioni sono globalmente ben poste per dati iniziali lisci. Risolto. Chiuso. Il mistero vive interamente in tre dimensioni, per entrambe le equazioni, ma per ragioni fondamentalmente diverse.

Il quadro della regolarità:

2D: Esistenza globale e unicità di soluzioni lisce sono note per entrambi i sistemi. Per Euler 2D con dati lisci, Wolibner (1933) dimostrò l’esistenza globale in spazi di Hölder; Yudovich (1963) stabilì l’unicità per dati con vorticità limitata. Per Navier-Stokes 2D, la regolarità globale segue dalla disuguaglianza di Ladyzhenskaya e dal principio del massimo per la vorticità.

Navier-Stokes 3D: Leray (1934) dimostrò l’esistenza globale di soluzioni deboli in $L^2$, ma unicità e regolarità restano aperte. Il teorema di Caffarelli, Kohn e Nirenberg (1982) mostra che l’insieme singolare ha misura di Hausdorff parabolica unidimensionale nulla, quindi qualunque blowup, se si verifica, è estremamente sparso. Il termine viscoso fornisce la stima a priori chiave $\int_0^T \|\nabla u\|_{L^2}^2 \, dt \leq C(u_0)$, ma questo controllo al livello dell'energia è al di sotto della scala critica 3D e non basta a chiudere un argomento di bootstrap. Vedi Perché Navier-Stokes è difficile per il divario di supercriticità.

Euler 3D: Non esiste una teoria globale per dati lisci. La buona posizione locale negli spazi di Sobolev $H^s$, $s > 5/2$, è classica (Kato 1972, Kato e Ponce 1988). Il criterio di Beale, Kato e Majda (1984) riduce il rilevamento del blowup a $\int_0^T \|\omega\|_{L^\infty} \, dt$: la soluzione resta liscia su $[0,T]$ se e solo se questo integrale è finito. Il blowup richiede che la vorticità cresca abbastanza rapidamente da non essere integrabile nel tempo. Elgindi (2021, Annals of Mathematics) ha dimostrato la formazione di singolarità in tempo finito per dati $C^{1,\alpha}$. Una vera svolta, ma al di sotto della soglia liscia ($C^\infty$). Se le soluzioni lisce di Euler esplodano in 3D è ancora aperto.

Turbolenza. È qui che il confronto Euler-vs-Navier-Stokes diventa fisicamente vivido, quasi tangibile.

In un flusso turbolento, l’energia entra alle grandi scale (la dimensione del tubo, dell’ala, della tempesta) e casca verso vortici sempre più piccoli. Questa è la cascata di energia, ed è uno dei fenomeni più impressionanti di tutta la fisica. Al fondo della cascata, la viscosità converte finalmente l’energia cinetica in calore. Fine della corsa.

Euler cattura la dinamica dell’intervallo inerziale: il trasferimento di energia attraverso le scale guidato dalla non linearità. Ma non ha alcun cutoff viscoso. Nessun fondo della cascata. Nessun meccanismo per convertire l’energia cinetica in calore a una scala definita. Se l’energia possa comunque dissiparsi nel limite inviscido, ciò che si chiama dissipazione anomala, resta una profonda questione aperta.

Ecco perché la modellizzazione della turbolenza usa quasi sempre Navier-Stokes. Il numero di Reynolds $\mathrm{Re} = UL/\nu$ dice quanto è ampia la cascata: un $\mathrm{Re}$ alto significa molte decadi di scale che separano l’immissione di energia dal consumo viscoso. La turbolenza reale vive nella tensione tra la cascata inviscida che riversa energia verso il basso e il cutoff viscoso che la distrugge alle scale più piccole.

La teoria di Kolmogorov del 1941 predice uno spettro di energia $E(k) \sim \varepsilon^{2/3} k^{-5/3}$ nell’intervallo inerziale, la regione in cui non dominano né il forzamento a grande scala né la dissipazione viscosa. La scala di dissipazione $\eta = (\nu^3/\varepsilon)^{1/4}$ fissa il fondo di questo intervallo. Sotto $\eta$, vince la viscosità.

Le equazioni di Euler sono il modello formale $\nu = 0$. Ma interpretare le soluzioni di Navier-Stokes come convergenti a soluzioni di Euler nel regime turbolento è davvero sottile, e la questione se la dissipazione di energia persista in questo limite (dissipazione anomala) è il contenuto della congettura di Onsager. Il lato rigido, che mostra assenza di dissipazione per $u \in C^{0,\alpha}$ con $\alpha > 1/3$, è stato dimostrato da Constantin, E e Titi nel 1994. Il lato flessibile, che costruisce soluzioni dissipative di Euler al di sotto di $C^{1/3}$, è stato completato da Isett nel 2018, basandosi sul programma di integrazione convessa di De Lellis e Székelyhidi.

Per Navier-Stokes, l’immagine della cascata è incorporata direttamente nel bilancio energetico: $\varepsilon = \nu \langle |\nabla u|^2 \rangle$. La domanda aperta è netta. Le soluzioni di Navier-Stokes restano lisce abbastanza a lungo perché la teoria statistica di Kolmogorov sia matematicamente giustificata? Il problema di esistenza e regolarità chiede, in parte, esattamente questo: se la dissipazione viscosa sia abbastanza forte da domare la cascata a ogni scala, per tutto il tempo.

La differenza tra Euler e Navier-Stokes è un termine: $\nu \Delta u$. Quel termine cambia tutto.

Euler non è un Navier-Stokes semplificato. È un sistema fondamentalmente diverso che capita condivida gran parte della sua struttura. E la scelta conta nella pratica: scegliere il modello sbagliato (Euler dove la viscosità conta, Navier-Stokes dove non conta) può rovinare completamente una simulazione. Per le equazioni complete, vedi Che cosa sono le equazioni di Navier-Stokes?. Per gli ostacoli, vedi Perché è difficile. Per il premio, vedi Il Problema del Millennio. Per flusso incomprimibile vs. comprimibile, vedi Navier-Stokes incomprimibile vs. comprimibile.

Il termine viscoso $\nu \Delta u$ trasforma un sistema non lineare del primo ordine (Euler) in uno parabolico semilineare (Navier-Stokes). Ciò che si ottiene: dissipazione dell’energia, stime a priori di ordine superiore e la struttura di semigruppo analitico alla base della teoria di Navier-Stokes costruita da Leray nel 1934 ed estesa da Fujita e Kato nel 1964.

Eppure questa regolarizzazione non è sufficiente in 3D per chiudere l’argomento di regolarità globale. Neanche lontanamente. Il Problema del Millennio Clay chiede precisamente se la regolarizzazione parabolica in Navier-Stokes possa controllare il trasferimento non lineare di energia attraverso tutte le scale, per tutto il tempo. Per Euler, la domanda parallela è altrettanto netta: l’assenza totale di regolarizzazione garantisce il blowup in tempo finito da dati lisci? Nessuno lo sa.

Entrambi i problemi si trovano al centro della teoria matematica dei fluidi, e confrontarli rivela esattamente che cosa la viscosità ti dà e che cosa no. La questione della regolarità 3D, per entrambi i sistemi, resta tra i problemi aperti più difficili di tutta l’analisi.