Equilibria of the Pressureless Euler System in Dimension 1

Sintesi Tecnica: Equilibri del sistema di Euler senza pressione in dimensione 1

Enunciato del problema
Questo articolo investiga il comportamento a lungo termine del sistema di Euler senza pressione (PE) unidimensionale, proposto originariamente da Zeldovich per modellare la formazione delle strutture su larga scala. Il sistema è dato dalle equazioni di continuità e quantità di moto:
$$\partial_t \rho + \partial_y(\rho v) = 0, \quad \partial_t(\rho v) + \partial_y(\rho v^2) = 0$$
dove $\rho$ è una distribuzione di massa (misura di probabilità con secondo momento finito) e $v$ è il campo di velocità. Gli autori si concentrano specificamente sulle Soluzioni a Particelle Appiccicose (SPS), una classe di soluzioni distribuzionali in cui le particelle di massa si uniscono al collidere. Sebbene l'esistenza di tali soluzioni sia stata stabilita tramite una formula di proiezione (Natile-Savaré), le condizioni sotto le quali queste soluzioni convergono a uno stato di equilibrio, la natura di tale equilibrio e il tempo necessario per il collasso sono rimaste questioni aperte. A differenza del sistema Euler-Poisson attrattivo, dove l'equilibrio è tipicamente una singola massa puntiforme, il sistema di Euler senza pressione permette una ricca varietà di distribuzioni di equilibrio.

Metodologia
L'analisi si basa pesantemente sulla formula di proiezione derivata in [10], che traduce il problema dallo "spazio delle posizioni" fisico allo "spazio della massa". Utilizzando la mappa di trasporto ottimale $N_t: (0,1) \to \mathbb{R}$ (dove $N_t \# \mathcal{L} = \rho_t$), l'evoluzione del sistema è caratterizzata dalla proiezione della traiettoria lineare $N_0 + tV_0$ sul cono convesso $K$ delle funzioni non decrescenti in $L^2(0,1)$:
$$N_t = \text{proj}_K(N_0 + tV_0)$$
dove $V_0 = v_0 \circ N_0$. Gli autori utilizzano strumenti di analisi convessa, specificamente le proprietà dell'inviluppo convesso ($\text{conv}$) e della derivata destra, per analizzare il comportamento asintotico di $N_t$. Definiscono una funzione primitiva $F(x) = \int_0^x V_0(w) dw$ e utilizzano l'invarianza galileiana del sistema per normalizzare la quantità di moto totale a zero.

Contributi Chiave e Risultati

1. Condizione Necessaria e Sufficiente per l'Equilibrio:
   L'articolo stabilisce una condizione precisa, indicata come $(\star)$, basata esclusivamente sui dati iniziali $(N_0, V_0)$, che determina se la SPS converga a un equilibrio. La condizione richiede:

   • Quantità di moto totale nulla: $\int_0^1 V_0(w) dw = 0$ (equivalente a $F(1)=0$).
   • Quantità di moto parziale non negativa: $\int_0^x V_0(w) dw \geq 0$ per ogni $x \in (0,1)$ (equivalente a $F(x) \geq 0$).

   Gli autori dimostrano che la SPS converge a un equilibrio nella metrica 1-Wasserstein se e solo se la condizione $(\star)$ è soddisfatta. Se $F(x) < 0$ per qualsiasi $x$, la soluzione diverge (la massa si diffonde all'infinito).

2. Caratterizzazione dello Stato di Equilibrio:
   Quando la condizione $(\star)$ è soddisfatta, l'articolo fornisce una descrizione completa della distribuzione limite $\rho_\infty$. L'equilibrio è determinato dagli zeri della funzione $F$:

   • Regioni Stazionarie: Su qualsiasi intervallo in cui $F$ è identicamente zero, la distribuzione di massa rimane stazionaria per tutto il tempo ($N_t(x) = N_0(x)$).
   • Regioni di Collasso: Su qualsiasi intervallo $(a, b)$ dove $F(a)=F(b)=0$ e $F > 0$ strettamente all'interno, la massa in quell'intervallo collassa in una singola misura di Dirac (massa puntiforme) situata nel centro di massa di quel "pezzo" iniziale di massa.
     Ciò implica che l'equilibrio è una sovrapposizione di massa stazionaria e masse puntiformi collassate.

3. Tempo di Collasso:
   Gli autori derivano una formula esplicita per il tempo di collasso ($T_c$), definito come il tempo minimo $T$ tale che $\rho_t = \rho_\infty$ per ogni $t \geq T$. Se la soluzione converge, $T_c$ è dato dal supremo su tutti gli intervalli $(a,b)$ dove $F$ è positiva:
   $$T_c = \sup_{[a,b] \in I} \sup_{x \in (a,b)} \frac{\int_x^b N_0(w) dw + \frac{x-b}{b-a}\int_a^b N_0(w) dw}{\int_0^x V_0(w) dw}$$
   Se questo valore è infinito, la soluzione si avvicina all'equilibrio asintoticamente ma non lo raggiunge mai in tempo finito.

4. Metriche di Convergenza:
   I risultati sono formulati principalmente in termini di distanza 1-Wasserstein ($W_1$). L'articolo nota che se il supporto della soluzione rimane limitato (il che accade se il supporto iniziale è limitato e $(\star)$ è soddisfatta), la convergenza in $W_p$ per ogni $p$ segue dalla $W_1$. Gli autori discutono l'estensione alla convergenza $W_2$, notando che, sebbene essa sia valida sotto specifiche condizioni di monotonicità vicino ai bordi dell'intervallo di massa, una condizione generale dipendente solo dai dati iniziali per la convergenza $W_2$ rimane una difficoltà aperta a causa della natura globale dell'operazione di inviluppo convesso.

Significatività e Rivendicazioni
L'articolo sostiene di fornire una caratterizzazione completa del comportamento di equilibrio per le soluzioni a particelle appiccicose per il sistema di Euler senza pressione in 1D utilizzando solo le condizioni iniziali. Esso risolve le tre questioni centrali riguardanti l'esistenza dell'equilibrio, la forma specifica di tale equilibrio e il tempo di collasso. Gli autori sottolineano che il loro approccio evita la necessità di assunzioni riguardanti il supporto compatto nel tempo, che erano precedentemente richieste nella letteratura correlata (ad es. [8]). Il lavoro si distingue dal sistema Euler-Poisson repulsivo, dove tali metodi predittivi falliscono a causa di comportamenti patologici, evidenziando come la natura guidata dalla quantità di moto del sistema di Euler senza pressione permetta questa determinazione precisa. I risultati si basano sulla formula di proiezione e sull'analisi convessa, offrendo un collegamento rigoroso tra la geometria della distribuzione di velocità iniziale e lo stato asintotico del sistema.