Global existence of weak solutions to incompressible anisotropic Cahn-Hilliard-Navier-Stokes system

Titolo: Esistenza globale di soluzioni deboli per il sistema anisotropo di Cahn-Hilliard-Navier-Stokes incomprimibile con densità variabile

Autori: Azeddine Zaidni, Saad Benjelloun, Radouan Boukharfane.

---

1. Il Problema Studiato

Il lavoro si concentra sull'analisi matematica di un sistema accoppiato che descrive la dinamica di fluidi bifase incomprimibili con densità variabile e energia superficiale anisotropa. Il sistema combina le equazioni di Navier-Stokes (per il flusso del fluido) con l'equazione di Cahn-Hilliard (per la separazione delle fasi), noto come sistema CHNS (Cahn-Hilliard-Navier-Stokes).

Le caratteristiche distintive del modello analizzato sono:

• Densità Variabile: La densità $\rho$ non è costante ma evolve secondo un'equazione di trasporto.
• Anisotropia: A differenza dei modelli classici isotropi, l'energia superficiale è modellata tramite una funzione omogenea di primo grado $\Gamma(\nabla \phi)$, dove $\phi$ è il parametro d'ordine (concentrazione). L'energia libera include un termine di gradiente anisotropo: $F = \int_\Omega \frac{\epsilon}{2} \Gamma^2(\nabla \phi) + \frac{1}{\epsilon} V(\phi)$.
• Dominio: Il sistema è definito in un dominio limitato e regolare $\Omega \subset \mathbb{R}^d$ con $d=2, 3$.
• Condizioni al contorno: Velocità nulla ($u=0$) e condizioni di Neumann omogenee per il potenziale chimico e il parametro d'ordine.

L'obiettivo principale è dimostrare l'esistenza globale di soluzioni deboli per questo sistema non lineare in due e tre dimensioni, un risultato che estende la letteratura precedente limitata principalmente al caso isotropo.

2. Metodologia

La dimostrazione si articola in diverse fasi rigorose, basate su tecniche di approssimazione e analisi funzionale:

1. Potenziale Logaritmico Regolarizzato:
   Poiché il potenziale termodinamico $F(\phi)$ (di tipo logaritmico, definito su $[-1, 1]$) presenta singolarità agli estremi, gli autori introducono una famiglia di potenziali regolarizzati $F_\epsilon(s)$ definiti su $\mathbb{R}$. Questo permette di lavorare con funzioni lisce ($C^2$) durante la fase di approssimazione.

2. Schema di Galerkin:
   Viene costruita una successione di soluzioni approssimate $(\rho_n, u_n, \phi_n, \mu_n)$ proiettando il sistema su sottospazi di dimensione finita generati dagli autovettori degli operatori di Stokes e di Laplace.

3. Esistenza Locale e Prolungamento Globale:

   • Fase Locale: L'esistenza locale delle soluzioni approssimate è garantita dal teorema di Cauchy-Lipschitz applicato a un sistema linearizzato.
   • Prolungamento Globale (Chiave): Per estendere la soluzione da un intervallo locale $[0, T_0)$ a un intervallo globale $[0, T]$, gli autori utilizzano un argomento di punto fisso combinato con la disuguaglianza di Bihari.
   • Viene analizzata una mappa che prende un campo di velocità come input e restituisce la soluzione del sistema linearizzato. Applicando la disuguaglianza di Bihari alle stime energetiche, si dimostra che la soluzione non può esplodere in tempo finito, garantendo l'esistenza globale per il sistema regolarizzato.

4. Stime Uniformi e Passaggio al Limite:

   • Vengono derivate stime energetiche a priori uniformi rispetto ai parametri di regolarizzazione $\epsilon$ e all'indice di approssimazione $n$.
   • Si utilizzano teoremi di compattezza (Ascoli-Arzelà, Aubin-Lions) per estrarre sottosuccessioni convergenti.
   • Un aspetto tecnico cruciale è il trattamento del termine anisotropo $\Gamma\xi(\nabla \phi) \otimes \nabla \phi$. Grazie alla linearità dell'operatore $\Gamma\xi$ (ipotesi H2) e alla forte convergenza di $\phi_n$ in spazi appropriati, si dimostra la convergenza forte di questo termine non lineare.

5. Limite $\epsilon \to 0$:
   Infine, si fa tendere il parametro di regolarizzazione $\epsilon$ a zero, recuperando la soluzione debole per il sistema originale con il potenziale logaritmico singolare.

3. Risultati Principali

Il teorema principale (Teorema 2.1) stabilisce che, sotto ipotesi ragionevoli sui dati iniziali ($\rho_0 \in L^\infty$, $u_0 \in L^2$, $\phi_0 \in H^1 \cap L^\infty$), esiste una quadrupla $(\rho, u, \phi, \mu)$ che soddisfa il sistema in senso debole, con le seguenti regolarità:

• Densità: $\rho \in C([0, T]; L^2) \cap L^\infty \cap W^{1,\infty}(H^{-1})$.
• Velocità: $u \in L^2(0, T; V) \cap B^{1/4}_{2,\infty}(0, T; H)$ (spazi di Besov che catturano la regolarità temporale).
• Parametro d'ordine: $\phi \in L^\infty(0, T; H^1) \cap B^{1/4}_{\infty,\infty}(0, T; L^2)$.
• Potenziale Chimico: $\mu \in L^2(0, T; H^1)$.

Inoltre, viene dimostrata la validità dell'equazione di bilancio energetico globale, che include l'energia cinetica, l'energia interfacciale anisotropa e l'energia libera di bulk.

4. Contributi Chiave

• Estensione all'Anisotropia: Questo è uno dei primi lavori che dimostra l'esistenza globale di soluzioni deboli per il sistema CHNS completo (Navier-Stokes + Cahn-Hilliard) in presenza di energia superficiale anisotropa e densità variabile. La maggior parte dei risultati precedenti si limitava al caso isotropo ($\Gamma(\nabla \phi) = |\nabla \phi|$).
• Gestione della Densità Variabile: Il trattamento della densità variabile in un contesto anisotropo richiede stime energetiche e tecniche di compattezza più sofisticate rispetto al caso a densità costante.
• Uso della Disuguaglianza di Bihari: L'applicazione strategica della disuguaglianza di Bihari per prolungare le soluzioni locali a soluzioni globali rappresenta un contributo metodologico significativo per questo tipo di sistemi non lineari.
• Trattamento del Potenziale Singolare: La costruzione rigorosa della soluzione attraverso la regolarizzazione del potenziale logaritmico e il successivo passaggio al limite.

5. Significato e Implicazioni

Questo lavoro ha un impatto significativo sia nella matematica pura che nella fisica applicata:

• Fondamenti Teorici: Colma una lacuna nella teoria delle equazioni alle derivate parziali non lineari, fornendo una base rigorosa per l'analisi di sistemi multifase complessi con anisotropia.
• Modellazione Fisica: L'anisotropia è fondamentale per descrivere fenomeni reali come la solidificazione di cristalli, la crescita di film sottili e la dinamica di gocce su superfici strutturate, dove l'energia superficiale dipende dalla direzione del gradiente di concentrazione.
• Simulazioni Numeriche: La prova di esistenza globale fornisce la giustificazione teorica necessaria per lo sviluppo e l'affidabilità di schemi numerici destinati a simulare questi sistemi fisici complessi in 2D e 3D.

In sintesi, il paper rappresenta un avanzamento sostanziale nella comprensione matematica dei flussi bifase complessi, unendo rigorosamente la termodinamica dei sistemi anisotropi con la dinamica dei fluidi viscosi.