On weak solutions for the stationary Cahn-Hillard-Navier-Stokes equations with singular potential

Sintesi Tecnica: Soluzioni Deboli per le Equazioni di Cahn-Hilliard-Navier-Stokes Stazionarie con Potenziale Singolare

Enunciato del Problema
Questo articolo indaga l'esistenza di soluzioni deboli per il sistema stazionario di Navier–Stokes–Cahn–Hilliard che descrive una miscela fluida bifase, comprimibile e tridimensionale con interfaccia diffusa. Il sistema è governato dalla conservazione della massa, della quantità di moto e della concentrazione, accoppiata a un'equazione del potenziale chimico. Una caratteristica distintiva del modello è l'inclusione di una densità di energia libera singolare di tipo Flory–Huggins (logaritmico). Questa scelta garantisce che la frazione di massa $c$ rimanga nell'intervallo fisico $(-1, 1)$ e permette la presenza di stati di vuoto ($\rho = 0$).

Il problema ai valori al contorno specifico coinvolge un dominio limitato $\Omega \subset \mathbb{R}^3$ con condizioni di Dirichlet per la velocità e condizioni di Neumann per la concentrazione e il potenziale chimico. Le principali sfide matematiche derivano da due fonti:

1. Il Contesto Stazionario: A differenza dei problemi evolutivi, il sistema stazionario manca di una disuguaglianza energetica a priori ($E(t) \leq E(0)$) guidata dall'evoluzione temporale, che tipicamente fornisce il controllo necessario sul potenziale singolare e sul campo di concentrazione.
2. Potenziale Singolare e Vuoto: La singolarità logaritmica nell'energia libera impedisce le regolarizzazioni standard tramite sviluppo di Taylor, e la potenziale degenerazione della densità vicino agli stati di vuoto complica gli argomenti di compattezza necessari per il passaggio al limite.

Metodologia
Gli autori stabiliscono l'esistenza attraverso un processo di limite a due livelli applicato a un sistema approssimato accuratamente costruito.

1. Regolarizzazione Specializzata del Potenziale:
   Per affrontare la mancanza di stime $L^2$ su $c$ e il fallimento degli sviluppi di Taylor standard, gli autori introducono una regolarizzazione specializzata $f^\delta_2(c)$ per il termine logaritmico singolare. A differenza di approcci precedenti, questa approssimazione è costruita in modo tale che la differenza $f^\delta_2(c) - \frac{\theta_c}{2}c^2$ cresca solo linearmente quando $|c| \to \infty$. Questa specifica costruzione è progettata per annullare la crescita quadratica indotta dal termine anti-diffusivo, ripristinando così la compattezza necessaria e permettendo stime uniformi.

2. Schema di Approssimazione a Due Livelli:

   • Livello 1 (limite $\epsilon$): Gli autori introducono inizialmente una viscosità artificiale e parametri di regolarizzazione ($\epsilon$) per costruire un sistema regolarizzato. Utilizzando un argomento di punto fisso (Teorema del Punto Fisso di Schaefer), dimostrano l'esistenza di soluzioni forti per questo sistema approssimato. Vengono derivate stime uniformi indipendenti da $\epsilon$, permettendo di passare al limite $\epsilon \to 0$, ottenendo una soluzione debole per un sistema intermedio con energia libera regolarizzata.
   • Livello 2 (limite $\delta$): Il secondo passo consiste nel rimuovere la regolarizzazione del potenziale ($\delta \to 0$). Una difficoltà maggiore qui è controllare l'integrabilità $L^{3/2}$ della densità. Gli autori impiegano una scelta specifica di pressione artificiale e un argomento di interpolazione delicato. Questo argomento sfrutta la struttura del potenziale regolarizzato, in particolare l'osservazione che $\frac{\partial f^\delta_2}{\partial c} - \theta_c c = O(\ln \frac{1}{\delta})$ per $|c| > 1+\delta$. Ciò permette di chiudere le stime per la norma della densità e stabilire limiti uniformi indipendenti da $\delta$.

3. Convergenza Forte e Limiti Fisici:
   Il passaggio al limite $\delta \to 0$ richiede di dimostrare la convergenza forte della densità per gestire il termine di pressione non monotono. Adattando tecniche dalla teoria delle equazioni di Navier–Stokes comprimibili stazionarie, gli autori stabiliscono la convergenza forte della densità. Inoltre, dimostrano che la concentrazione limite $c$ soddisfa i limiti fisici stretti $-1 < c < 1$ quasi ovunque sul supporto della densità ($\{\rho > 0\}$), garantendo che il termine singolare sia ben definito.

Risultati Chiave
Il teorema principale (Teorema 1.1) stabilisce l'esistenza di una soluzione debole $(\rho, u, \mu, c)$ al problema stazionario alle seguenti condizioni:

• Dominio: Un dominio limitato e semplicemente connesso con bordo $C^2$.
• Forze: Le forze esterne $g_1, g_2$ sono limitate.
• Esponente Adiabatico ($\gamma$):
  • $\gamma > 3/2$ se la forza esterna è irrotazionale ($\nabla \times g_1 = 0$).
  • $\gamma > 5/3$ se la forza esterna è rotazionale ($\nabla \times g_1 \neq 0$).
• Vincoli di Massa: La massa totale $m_1 > 0$ e la massa relativa $m_2$ soddisfano $m_2 \in (-m_1, m_1)$.

La soluzione soddisfa l'equazione di continuità renormalizzata, la formulazione debole delle equazioni della quantità di moto e della concentrazione, e la disuguaglianza energetica. Crucialmente, la concentrazione soddisfa $-1 < c(x) < 1$ quasi ovunque dove $\rho(x) > 0$.

Significato e Affermazioni
Gli autori affermano che questo lavoro fornisce il primo risultato di esistenza per soluzioni deboli del sistema stazionario di Navier–Stokes–Cahn–Hilliard comprimibile che incorpora simultaneamente:

1. Un potenziale di energia libera singolare (di tipo Flory–Huggins), che è fisicamente coerente per miscele binarie ma matematicamente impegnativo a causa della mancanza di limiti globali sulla concentrazione.
2. La presenza di stati di vuoto ($\rho = 0$), che introduce degenerazione nelle equazioni.

L'articolo nota che, sebbene l'intervallo degli esponenti adiabatici ammissibili $\gamma$ possa essere estendibile basandosi sui risultati per le equazioni di Navier–Stokes comprimibili stazionarie per fluidi singoli, l'intervallo attuale è limitato dalla difficoltà di ottenere stime pesate per la concentrazione vicino al bordo a causa della singolarità dell'energia libera. Il lavoro rappresenta un passo significativo nel colmare il divario tra modelli evolutivi (dove le disuguaglianze energetiche forniscono controllo) e modelli stazionari con potenziali singolari.