On Lions' density patch problem at a critical level of regularity

1. Il Problema: Il "Density Patch" di Lions

Il paper si concentra sull'evoluzione di un fluido incomprimibile omogeneo che occupa una regione limitata $D \subset \mathbb{R}^2$, circondata dal vuoto. Questo scenario è noto come problema del density patch (o "patch di densità").
La configurazione è descritta dalle equazioni di Navier-Stokes non omogenee:
$$\begin{cases} \partial_t(\rho u) + \text{div}(\rho u \otimes u) - \nu \Delta u + \nabla P = 0, \\ \partial_t \rho + \text{div}(\rho u) = 0, \\ \text{div } u = 0, \end{cases}$$
con dati iniziali $\rho(0, x) = \mathbf{1}_D$ (la funzione caratteristica del dominio $D$) e una velocità iniziale $u_0$.

La questione centrale, sollevata da J.-L. Lions, riguarda la propagazione della regolarità geometrica dell'interfaccia. Se il dominio iniziale $D$ ha un bordo Lipschitziano (o più regolare, come $C^{1,\gamma}$), questa regolarità viene preservata dall'evoluzione temporale?
Il lavoro affronta questo problema in un quadro di regolarità critico, ovvero nello spazio di Besov $\dot{B}^0_{2,1}(\mathbb{R}^2)$ per la velocità iniziale. Questo spazio è critico rispetto alla scalatura naturale delle equazioni e rappresenta il limite inferiore per garantire il controllo della geometria dell'interfaccia, superando le limitazioni dello spazio $L^2$.

2. Metodologia e Strategie Chiave

Gli autori sviluppano una teoria di ben-postezza globale (esistenza, unicità, stabilità) combinando diverse tecniche avanzate di analisi funzionale e PDE:

• Decomposizione Atomica e "Gluing": Per costruire soluzioni con dati iniziali in $L^2$ (e successivamente in $\dot{B}^0_{2,1}$), gli autori utilizzano una procedura di "lifting" per approssimare la densità di vuoto con densità strettamente positive. Successivamente, sfruttano una decomposizione atomica della velocità iniziale $u_0 = \sum u_{0,j}$, dove ogni atomo appartiene a spazi di Sobolev sub-critici ($\dot{H}^\eta \cap \dot{H}^{-\eta}$). Risolvono il sistema linearizzato per ogni atomo e "incollano" le soluzioni per ottenere la soluzione globale. Questo approccio permette di recuperare proprietà strutturali forti (come l'uguaglianza dell'energia) che non sono garantite dalle soluzioni deboli classiche.
• Disuguaglianza di Gagliardo-Nirenberg Dinamica e Localizzata: Una difficoltà principale è la mancanza di integrabilità della velocità fuori dal patch (in $\mathbb{R}^2$). Gli autori dimostrano una disuguaglianza di tipo Gagliardo-Nirenberg adattata al dominio mobile. Essi mostrano che la dinamica è confinata in una regione di raggio $\sqrt{t}$ attorno al dominio iniziale ($D_{R,t} = D + B_{R\sqrt{t}}(0)$). All'interno di questa striscia dinamica, è possibile recuperare stime critiche di integrabilità pesate nel tempo, essenziali per chiudere gli argomenti di unicità.
• Trasformata di Galilei: Per ottenere la regolarità Lipschitziana globale ($\nabla u \in L^1_t L^\infty_x$), gli autori introducono una trasformazione di Galilei che rimuove il moto rigido associato alla quantità di moto media $M$. Questo permette di studiare il decadimento della velocità in un sistema di riferimento in movimento, eliminando il termine costante che altrimenti impedirebbe l'integrabilità temporale del gradiente.
• Stabilità e Unicità Debole-Forte: Viene stabilita una stima di stabilità quantitativa tra una soluzione "forte immediata" (costruita con dati in $\dot{B}^0_{2,1}$) e una soluzione debole generica. Questo porta a un risultato di unicità nel caso in cui i dati iniziali coincidano.

3. Risultati Principali

I teoremi principali del lavoro (Teoremi 1.1, 1.4, 1.5, 1.6, 1.7) possono essere riassunti come segue:

• Esistenza e Unicità Globale (Teorema 1.5): Per dati iniziali $u_0 \in \dot{B}^0_{2,1}(\mathbb{R}^2)$ e $\rho_0 = \mathbf{1}_D$ (con $D$ dominio Lipschitziano limitato), esiste un'unica soluzione globale $(\rho, u)$ che soddisfa l'uguaglianza dell'energia e possiede regolarità sufficiente per garantire $\nabla u \in L^1((0, \infty); L^\infty(\mathbb{R}^2))$.
• Propagazione della Regolarità del Bordo: La soluzione garantisce che la mappa di flusso $X(t, \cdot)$ sia Lipschitziana. Di conseguenza, se il dominio iniziale $D$ è Lipschitziano, allora il dominio evoluto $D_t = X(t, D)$ rimane Lipschitziano per ogni $t > 0$. Inoltre, se $u_0 \in \dot{B}^\gamma_{2,1}$ con $\gamma \in (0,1)$ e $D$ è di classe $C^{1,\gamma}$, allora $D_t$ rimane di classe $C^{1,\gamma}$ per ogni tempo (Teorema 1.6).
• Comportamento Asintotico (Teorema 1.7): Il sistema mostra un comportamento asintotico rigido. La velocità converge esponenzialmente alla sua media spaziale $M$ (quantità di moto per unità di massa) all'interno del patch:
  $$\|u(t) - M\|_{L^2(D_t)} \lesssim \|u_0 - M\|_{L^2(D)} e^{-\lambda t}.$$
  Di conseguenza, il dominio evolve come una traslazione rigida con velocità costante $M$, più una deformazione finita che converge a un dominio limite $D_\infty$. In altre parole, la viscosità smorza le deformazioni interne, lasciando solo il moto rigido.
• Stabilità: Viene dimostrata la stabilità delle soluzioni rispetto alle perturbazioni dei dati iniziali nello spazio critico, confermando la ben-postezza del problema.

4. Significato e Contributi Scientifici

Questo lavoro rappresenta un avanzamento significativo nella teoria dei fluidi incomprimibili non omogenei:

1. Risoluzione del Caso Critico: È il primo risultato che stabilisce la ben-postezza globale e la propagazione della regolarità del bordo per il problema del density patch nel caso critico $\dot{B}^0_{2,1}$ in presenza di vuoto. Risultati precedenti richiedevano dati in $H^1$ o spazi sub-critici.
2. Superamento dell'Ostruzione $L^2$: Dimostra che lo spazio $L^2$ non è sufficiente per preservare la regolarità Lipschitziana del bordo in due dimensioni a causa della mancanza di integrabilità temporale del gradiente di velocità. L'uso dello spazio di Besov critico è essenziale per ottenere il controllo necessario.
3. Caratterizzazione Asintotica: Fornisce una descrizione quantitativa e rigorosa del comportamento a lungo termine del fluido, mostrando come la dissipazione viscosa porti a un moto puramente rigido, risolvendo una questione aperta sulla dinamica a lungo termine delle patch di densità.
4. Nuovi Strumenti Analitici: Le tecniche sviluppate, in particolare la disuguaglianza di Gagliardo-Nirenberg localizzata e l'uso della trasformata di Galilei combinata con la decomposizione atomica, offrono nuovi strumenti per lo studio di problemi di interfaccia libera e fluidi con densità discontinue.

In sintesi, Škondrić e Violini confermano la congettura di Lions per una classe di dati iniziali critica, dimostrando che la regolarità geometrica dell'interfaccia è preservata globalmente nel tempo e che il sistema evolve verso uno stato asintotico di moto rigido.