Navier-Stokes incomprimibile vs comprimibile

"Incomprimibile vs comprimibile" si riduce alla densità. Resta costante, oppure cambia?

Prova. Riempi una siringa d'acqua e spingi lo stantuffo. L'acqua si muove, ma in condizioni quotidiane non si comprime in modo apprezzabile. L'acqua resiste alla compressione così fortemente che trattarla come incomprimibile è un'approssimazione eccellente. Incomprimibile. Ora riempi quella siringa d'aria e sigilla l'estremità. Spingi lo stantuffo e sentirai l'aria cedere, la stessa massa d'aria ora stipata in un volume minore mentre si comprime sotto il tuo pollice. Questo è flusso comprimibile.

Nel flusso incomprimibile, la densità $\rho$ è costante in tutto il fluido, e ogni minuscola particella mantiene il proprio volume mentre si muove nello spazio. Il flusso comprimibile è diverso. La densità diventa una variabile, libera di cambiare da luogo a luogo e da istante a istante. Aria intorno a un motore a reazione, gas in un'esplosione, l'atmosfera su grandi scale: tutti comprimibili, tutti guidati da variazioni di densità.

Perché importa? Perché questa distinzione rimodella l'aspetto delle equazioni di Navier-Stokes, ciò che predicono e quanto sono difficili da analizzare e risolvere.

L'ipotesi di incomprimibilità impone che il campo di densità $\rho$ sia una costante positiva in tutto il dominio del flusso. Fisicamente, questo significa che ogni elemento materiale di volume conserva il proprio volume sotto la mappa di flusso. Nel caso comprimibile, $\rho(x,t)$ diventa un'incognita a tutti gli effetti, governata da una propria equazione di evoluzione.

Queste non sono due notazioni per lo stesso sistema. Differiscono nel numero di incognite, nella struttura delle equazioni di vincolo e nel carattere della pressione. Il sistema incomprimibile ha $d+1$ incognite scalari ($u$ e $p$ in $\mathbb{R}^d$); il sistema comprimibile ha tipicamente $d+2$ campi primari (velocità, densità ed energia interna o temperatura), con la pressione determinata tramite un'equazione di stato.

La distinzione cambia anche il tipo PDE della pressione: ellittico nel modello incomprimibile, dove la pressione è determinata da una risoluzione spaziale istantanea, contro una struttura mista iperbolico-parabolica nel caso comprimibile viscoso. I modi acustici nella parte inviscida o iperbolica si propagano alla velocità finita del suono, mentre viscosità e conduzione termica introducono diffusione parabolica dal punto di vista matematico.

Le equazioni di Navier-Stokes incomprimibili descrivono fluidi la cui densità è costante. Sono la versione che compare nel Problema del Millennio del Clay e la versione su cui si concentra questo sito.

Il sistema ha due parti. L'equazione della quantità di moto:

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

e il vincolo di incomprimibilità:

$$\nabla \cdot u = 0$$

Il vincolo $\nabla \cdot u = 0$ dice che il campo di velocità è a divergenza nulla: il fluido non si accumula né si rarefa in alcun punto. Qualunque cosa entri in una minuscola regione deve uscirne allo stesso tasso. Questa singola condizione sostituisce l'intera equazione della densità. La densità non cambia, quindi non serve un'equazione per seguirla.

Qui la pressione svolge un ruolo speciale. Non è determinata da una legge termodinamica (come la legge dei gas ideali). Invece, si aggiusta istantaneamente ovunque per mantenere il flusso a divergenza nulla. Matematicamente, $p$ risolve un'equazione di Poisson derivata dal vincolo. Le variazioni di pressione si propagano a velocità infinita. Non c'è "velocità del suono" nel flusso incomprimibile.

Il sistema di Navier-Stokes incomprimibile ha due campi incogniti: la velocità $u$ e la pressione $p$. Questa semplicità è ingannevole. Il termine non lineare $(u \cdot \nabla)u$ rende comunque il sistema estremamente difficile in tre dimensioni.

Il sistema di Navier-Stokes incomprimibile su $\mathbb{R}^3$ con viscosità cinematica $\nu > 0$:

$$\partial_t u + (u \cdot \nabla)u = -\nabla p + \nu \Delta u + f, \qquad x \in \mathbb{R}^3,\; t > 0,$$

$$\nabla \cdot u = 0, \qquad u(x,0) = u_0(x).$$

La condizione di divergenza nulla $\nabla \cdot u = 0$ è un vincolo puntuale, non un'equazione di evoluzione. Codifica la conservazione locale del volume: la mappa di flusso $\Phi_t$ soddisfa $\det(D\Phi_t) = 1$ per ogni $t$.

Applicando l'operatore divergenza all'equazione della quantità di moto e usando l'incomprimibilità si ottiene l'equazione di Poisson per la pressione:

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

Questa è un'equazione ellittica per $p$ a ogni istante fissato. La pressione non è una variabile termodinamica indipendente; nell'interpretazione PDE standard, agisce come moltiplicatore di Lagrange che impone il vincolo di divergenza nulla, determinato globalmente e istantaneamente dal campo di velocità. L'informazione si propaga a velocità infinita attraverso la pressione, una differenza strutturale rispetto al sistema comprimibile che non può essere mascherata.

Le incognite sono $u : \mathbb{R}^3 \times [0,T) \to \mathbb{R}^3$ e $p : \mathbb{R}^3 \times [0,T) \to \mathbb{R}$. Per la formulazione del Clay (Fefferman, 2000), $u_0 \in C^\infty(\mathbb{R}^3)$ è a divergenza nulla e la domanda è se $u$ rimanga in $C^\infty(\mathbb{R}^3 \times [0,\infty))$ con energia limitata.

Le equazioni di Navier-Stokes comprimibili governano flussi in cui la densità varia. Sistema più grande. Più incognite. Più equazioni.

Hai ancora un'equazione della quantità di moto, ma ora la densità $\rho$ compare esplicitamente:

$$\partial_t (\rho u) + \nabla \cdot (\rho u \otimes u) = -\nabla p + \nabla \cdot \tau + \rho f$$

Il vincolo $\nabla \cdot u = 0$ scompare. Al suo posto, ottieni un'equazione di continuità che tiene traccia dell'evoluzione della densità:

$$\partial_t \rho + \nabla \cdot (\rho u) = 0$$

Questo dice che la massa si conserva: la densità cambia perché il flusso comprime o espande le particelle fluide.

Il sistema ha anche bisogno di un'equazione dell'energia e di un'equazione di stato, una relazione termodinamica come $p = \rho R T$ (la legge dei gas ideali) che lega la pressione a densità e temperatura. La pressione non è più un esecutore passivo di un vincolo. La sua parte acustica viaggia alla velocità del suono, mentre il sistema viscoso completo contiene anche effetti diffusivi.

Il sistema comprimibile è essenziale per l'aerodinamica ad alte velocità, la gasdinamica astrofisica, la combustione e qualunque flusso in cui le variazioni di densità contano. Ma è un oggetto matematico genuinamente diverso dalle equazioni incomprimibili. Più incognite, più equazioni, una struttura PDE del tutto diversa.

Il sistema di Navier-Stokes comprimibile accoppia la velocità $u(x,t)$, la densità $\rho(x,t)$, la pressione $p(x,t)$ e l'energia interna specifica $e(x,t)$ (o temperatura $\theta$). In forma conservativa:

Continuità: $$\partial_t \rho + \nabla \cdot (\rho u) = 0$$

Quantità di moto: $$\partial_t (\rho u) + \nabla \cdot (\rho u \otimes u) + \nabla p = \nabla \cdot \tau + \rho f$$

Energia: $$\partial_t (\rho E) + \nabla \cdot ((\rho E + p)u) = \nabla \cdot (\tau \cdot u) + \nabla \cdot (\kappa \nabla \theta) + \rho f \cdot u$$

dove $E = e + \tfrac{1}{2}|u|^2$ è l'energia specifica totale, $\tau$ è il tensore degli sforzi viscosi (per un fluido newtoniano, $\tau = \mu(\nabla u + \nabla u^T) + \lambda (\nabla \cdot u)I$ con viscosità di volume $\lambda$), e $\kappa$ è la conducibilità termica.

La chiusura richiede un'equazione di stato, ad esempio $p = (\gamma - 1)\rho e$ per un gas ideale con indice adiabatico $\gamma$.

Nella parte inviscida o iperbolica del sistema comprimibile, le perturbazioni acustiche si propagano a velocità finita, con velocità del suono caratteristica $c = \sqrt{\partial p / \partial \rho |_s}$. Il sistema completo Navier-Stokes-Fourier viscoso è però misto iperbolico-parabolico, quindi la diffusione viscosa e termica introduce propagazione infinita dal punto di vista matematico. Le equazioni comprimibili supportano onde d'urto, rarefazioni e discontinuità di contatto che non hanno analogo nel flusso incomprimibile.

Quando conta la comprimibilità? Decide un solo numero: il numero di Mach.

$$\text{Ma} = \frac{|u|}{c}$$

$|u|$ è la velocità del flusso. $c$ è la velocità del suono. Il loro rapporto ti dice quanto rapidamente si muove il flusso rispetto alla velocità con cui le perturbazioni di pressione possono propagarsi nel mezzo, e quel confronto determina se puoi ignorare in sicurezza le variazioni di densità o se esse domineranno la fisica.

Quando $\text{Ma} < 0.3$, density changes by less than about 5%. Incompressible equations work. Air in a room, water in a pipe, wind around a building: all low-Mach flows where pressure disturbances travel so much faster than the flow itself that density barely budges.

Sopra $\text{Ma} \approx 0.3$, la comprimibilità comincia a farsi sentire, e intorno a $\text{Ma} \approx 1$ si entra nel regime transonico, dove compaiono sacche supersoniche locali e si formano onde d'urto. Caccia militari. Ugelli di razzi. Veicoli spaziali in rientro.

Non è un interruttore binario. La maggior parte dei flussi fluidi quotidiani, e il Problema del Millennio del Clay, si collocano saldamente nel regime a basso Mach in cui si applicano le equazioni incomprimibili.

Il numero di Mach $\text{Ma} = |u|/c$ parametrizza l'importanza della comprimibilità, dove $c = \sqrt{\partial p / \partial \rho |_s}$ è la velocità del suono isentropica. Formalmente, le equazioni incomprimibili emergono come limite a basso Mach del sistema comprimibile.

L'espansione asintotica in potenze di $\text{Ma}^2$ (si vedano Klainerman & Majda, 1981, 1982; Schochet, 1986) mostra che, quando $\text{Ma} \to 0$ con dati iniziali opportuni, le soluzioni comprimibili convergono alla soluzione incomprimibile. Nelle adimensionalizzazioni standard a basso Mach con dati ben preparati, si scrive spesso la pressione come un fondo termodinamico quasi spazialmente uniforme più una correzione dinamica più piccola che impone il vincolo di incomprimibilità nel limite.

Questo è un limite singolare: la velocità del suono $c \to \infty$ e il carattere iperbolico del sistema comprimibile degenera nell'equazione ellittica della pressione del flusso incomprimibile. I modi acustici diventano infinitamente veloci e si disaccoppiano dalla dinamica vorticosa.

Il regime $\text{Ma} < 0.3$ is an engineering heuristic. It reflects the empirical observation that relative density variation $\delta\rho / \rho \sim \text{Ma}^2 / 2$ stays below $\sim$5% in this range. The mathematical justification is the convergence theorem for the low-Mach limit, which requires well-prepared initial data and compatible boundary conditions.

Il Problema del Millennio del Clay pone una domanda precisa: data una velocità iniziale liscia e a divergenza nulla su $\mathbb{R}^3$, il sistema di Navier-Stokes incomprimibile produce sempre una soluzione liscia che esiste per ogni tempo?

Perché proprio incomprimibile? Tre ragioni.

Primo, è già abbastanza difficile. Le equazioni incomprimibili 3D hanno resistito alla dimostrazione della regolarità globale fin dal lavoro fondamentale di Leray nel 1934. Aggiungere densità variabile, termodinamica e onde d’urto renderebbe il problema enormemente più difficile, non più trattabile.

Secondo, la difficoltà è pura meccanica dei fluidi. Il sistema incomprimibile isola la sfida matematica centrale, la competizione tra l’avvezione non lineare $(u \cdot \nabla)u$ e la dissipazione viscosa $\nu \Delta u$, senza complicazioni termodinamiche o acustiche. È l’arena più pulita in cui porre la questione della regolarità.

Terzo, la fisica è pulita. Le equazioni incomprimibili modellano i flussi quotidiani più comuni. Se possano produrre singolarità a partire da dati lisci è una questione fondamentale sulla coerenza matematica della meccanica dei fluidi classica.

Il sistema comprimibile ha i suoi profondi problemi aperti (esistenza di soluzioni globali con dati grandi, formazione e interazione di onde d’urto), ma sono problemi diversi con strutture diverse. Il premio Clay riguarda il caso incomprimibile perché questa è la specifica questione di regolarità formulata da Fefferman per Navier-Stokes 3D.

La formulazione ufficiale del Clay (Fefferman, 2000) specifica il sistema incomprimibile su $\mathbb{R}^3$:

$$\partial_t u + (u \cdot \nabla)u = -\nabla p + \nu \Delta u, \qquad \nabla \cdot u = 0, \qquad u|_{t=0} = u_0,$$

con $u_0 \in C^\infty(\mathbb{R}^3)$ a divergenza nulla, e la questione è se $u \in C^\infty(\mathbb{R}^3 \times [0,\infty))$ con $\int_{\mathbb{R}^3} |u(x,t)|^2\,dx$ limitata per ogni $t \geq 0$.

La scelta del sistema incomprimibile è motivata matematicamente. Il problema aperto chiave, il divario tra le soluzioni deboli di Leray-Hopf (che esistono globalmente ma possono non essere uniche o lisce) e le soluzioni classiche lisce (che esistono localmente ma possono esplodere), è specifico delle equazioni incomprimibili 3D. In 2D, la regolarità globale per soluzioni lisce incomprimibili di Navier-Stokes è nota; il caso 3D rimane aperto.

Il sistema comprimibile introduce difficoltà qualitativamente diverse: formazione di onde d’urto (che si verifica anche per le equazioni di Eulero con dati lisci), stati di vuoto ($\rho \to 0$) e l’accoppiamento tra vorticità e modi acustici. Sono problemi aperti importanti, ma strutturalmente distinti dalla questione della regolarità incomprimibile.

Il problema incomprimibile isola la competizione tra la non linearità supercritica rispetto all’energia e la dissipazione viscosa. La stima naturale dell’energia dà $u \in L^\infty_t L^2_x \cap L^2_t \dot{H}^1_x$, che in 3D manca di mezza derivata rispetto allo spazio critico per scaling $L^\infty_t \dot{H}^{1/2}_x$. Colmare questo divario, o dimostrare che non può essere colmato, è il cuore del Problema del Millennio.

Comincia da qui. Vuoi che ogni termine del sistema incomprimibile sia scomposto, con il significato fisico e il ruolo matematico di ciascuna parte spiegati da zero? Che cosa sono le equazioni di Navier-Stokes?

Da dove viene questo sistema? Derivazione delle equazioni di Navier-Stokes.

Elimina la viscosità e ottieni le equazioni di Eulero, più vecchie di un secolo, più semplici sulla pagina e, per certi versi, persino più difficili da comprendere matematicamente perché si perde l’effetto regolarizzante del termine diffusivo. Eulero vs. Navier-Stokes.

Il premio. Il problema di esistenza e regolarità per Navier-Stokes.

Percorsi suggeriti da qui:

• Che cosa sono le equazioni di Navier-Stokes?, il sistema incomprimibile completo con analisi termine per termine
• Derivazione delle equazioni di Navier-Stokes, dagli assiomi della meccanica dei continui al sistema di PDE
• Eulero vs. Navier-Stokes, il limite non viscoso $\nu \to 0$ e il ruolo della viscosità nella regolarità
• Il problema di esistenza e regolarità per Navier-Stokes, la formulazione precisa del Clay, la teoria di Leray-Hopf e il divario tra soluzioni deboli e forti