Quantitative Homogenization of a Cahn--Hilliard System with Source Term in Periodically Perforated Domains

Sintesi Tecnica: Omogeneizzazione Quantitativa di un Sistema di Cahn–Hilliard con Termine Sorgente in Domini Periodicamente Perforati

Enunciato del Problema
Questo articolo investiga l'omogeneizzazione qualitativa e quantitativa di un sistema di Cahn–Hilliard caratterizzato da un termine sorgente monotono non conservativo all'interno di un dominio periodicamente perforato. Il modello microscopico descrive la separazione di fase in un mezzo poroso in cui il fluido occupa uno spazio dei pori connesso $\Omega^\varepsilon_p \subset \Omega \subset \mathbb{R}^d$ ($d \in \{2,3\}$), ottenuto rimuovendo un array periodico di inclusioni solide da un dominio Lipschitz limitato. Il sistema accoppia il parametro d'ordine $\phi^\varepsilon$ e il potenziale chimico $\mu^\varepsilon$ tramite:
$$\begin{cases} \partial_t \phi^\varepsilon - \Delta \mu^\varepsilon + G(\phi^\varepsilon) = 0, \\ \mu^\varepsilon = -\Delta \phi^\varepsilon + F'(\phi^\varepsilon), \end{cases}$$
soggetto a condizioni al contorno di Neumann omogenee sia sul bordo esterno $\partial \Omega$ che sulle pareti dei pori $\Gamma^\varepsilon_s$. Il potenziale $F$ soddisfa una crescita polinomiale e coercitività (ad esempio, il classico potenziale a doppia buca), mentre il termine sorgente $G$ è globalmente Lipschitziano, monotono e soddisfa $G(0)=0$. La sfida principale consiste nel derivare le equazioni macroscopiche efficaci al limite della scala dei pori $\varepsilon \to 0$ ed estendere i tassi di convergenza precisi per l'errore tra le soluzioni microscopiche e quelle omogeneizzate.

Metodologia
L'analisi impiega il metodo dell'unfolding periodico adattato ai domini perforati, utilizzando l'operatore di unfolding $T^\varepsilon_*$ per gestire la geometria oscillante. La metodologia procede in due fasi principali:

1. Omogeneizzazione Qualitativa:

   • Stime di Energia: Vengono derivate stime a priori uniformi per le soluzioni microscopiche. Un'osservazione strutturale chiave è che il termine sorgente monotono $G$ (con $G(0)=0$) contribuisce con un termine dissipativo nell'identità dell'energia. Nello specifico, testando l'equazione del potenziale chimico con $G(\phi^\varepsilon)$ si ottiene un termine di gradiente positivo, che rinforza la dissipazione anziché competere con essa. Ciò consente di ottenere limiti uniformi indipendenti da $\varepsilon$ senza richiedere un argomento di Gronwall con fattori esponenziali temporali.
   • Convergenza a Due Scale: Utilizzando l'operatore di unfolding e argomenti di compattezza (lemma di Aubin–Lions–Simon), gli autori dimostrano la convergenza di $(\phi^\varepsilon, \mu^\varepsilon)$ verso una coppia omogeneizzata $(\phi, \mu)$. Il sistema limite è un'equazione di Cahn–Hilliard sul dominio macroscopico $\Omega$ con un tensore di diffusione efficace $B_{\text{hom}}$.
   • Problemi di Cella: Il tensore efficace $B_{\text{hom}}$ è caratterizzato da problemi di cella di Neumann scalari sulla cella di riferimento $Y_p$. I correttori $\chi_i$ appartengono a $H^1_{\text{per},0}(Y_p)$.

2. Analisi Quantitativa (Stime dei Correttori):

   • Operatore di Scissione di Scala: Per ottenere tassi quantitativi, l'articolo costruisce approssimazioni correttrici del primo ordine $\Phi^\varepsilon$ e $M^\varepsilon$ utilizzando l'operatore di scissione di scala $Q_1$ $Q_\varepsilon$ (da Griso [29, 30]). Questo operatore agisce sul dominio macroscopico fisso $\Omega$ ed è successivamente ristretto a $\Omega^\varepsilon_p$. Fondamentalmente, questo approccio permette ai correttori di cella di rimanere in $H^1_{\text{per}}(Y_p)$, evitando le restrittive assunzioni $W^{1,\infty}$ richieste dalle classiche espansioni formali.
   • Metodo del Residuo Variazionale: L'analisi dell'errore è condotta in un framework di norma negativa. L'errore del parametro di fase $e^\varepsilon_\phi = \phi^\varepsilon - \Phi^\varepsilon$ viene testato contro l'inverso del Laplaciano di Neumann $N_\varepsilon e^\varepsilon_\phi$, mentre l'errore del potenziale chimico viene testato contro $e^\varepsilon_\phi$. Ciò evita la necessità di una maggiore regolarità dei residui.
   • Stime di Consistenza: La dimostrazione si basa sulla derivazione di stime del residuo variazionale (Lemma 4.6) che tengono conto dello strato limite delle celle incomplete vicino a $\partial \Omega$.

Risultati Chiave

• Convergenza Qualitativa (Teorema 3.7): Le soluzioni microscopiche convergono (fino a sottosequenze) al sistema omogeneizzato:
  $$\begin{cases} \partial_t \phi - \text{div}(B_{\text{hom}}\nabla \mu) + G(\phi) = 0, \\ \mu = -\text{div}(B_{\text{hom}}\nabla \phi) + F'(\phi), \end{cases}$$
  con condizioni al contorno di Neumann efficaci. Il tensore $B_{\text{hom}}$ è simmetrico ed uniformemente ellittico.

• Tassi Quantitativi (Teorema 4.7): Sotto un'ulteriore assunzione di regolarità $H^2$ sulla soluzione omogeneizzata e dati iniziali ben preparati, l'articolo stabilisce una stima dell'errore di ordine $\varepsilon^{1/2}$:
  $$\|\phi^\varepsilon - \phi\|_{L^2(0,T; L^2(\Omega^\varepsilon_p))} + \|\nabla \phi^\varepsilon - (\nabla \phi + \sum Q_\varepsilon(\partial_{x_i}\phi)\nabla_y \chi_i(\tfrac{x}{\varepsilon}))\|_{L^2((0,T)\times\Omega^\varepsilon_p)} + \|\mu^\varepsilon - \mu\|_{L^2(0,T; H^1(\Omega^\varepsilon_p)')} \leq C\varepsilon^{1/2}.$$

  • L'errore non corretto del parametro d'ordine è controllato in $L^2(0,T; L^2(\Omega^\varepsilon_p))$.
  • L'errore del gradiente corretto è controllato in $L^2((0,T)\times\Omega^\varepsilon_p)$.
  • L'errore del potenziale chimico è controllato nella norma duale $L^2(0,T; H^1(\Omega^\varepsilon_p)')$.

• Effetti dello Strato Limite: Il tasso $\varepsilon^{1/2}$ su domini limitati è attribuito allo strato limite delle celle incomplete vicino a $\partial \Omega$. Sul toro piatto $\mathbb{T}^d$, dove non esiste tale strato limite, il tasso migliora a $O(\varepsilon)$ (Nota 4.10).

Significatività e Contributi
L'articolo apporta due contributi primari alla teoria dell'omogeneizzazione per equazioni di campo di fase in mezzi perforati:

1. Estensione Qualitativa: Deriva rigorosamente il sistema di Cahn–Hilliard omogeneizzato con un termine sorgente non conservativo monotono generale, dimostrando che la struttura dissipativa del termine sorgente semplifica le stime di energia rispetto a una forza non conservativa generica.
2. Miglioramento Quantitativo: Migliora il tasso di convergenza noto per i sistemi di Cahn–Hilliard in domini perforati da $\varepsilon^{1/4}$ (stabilito da Schmuck–Kalliadasis [49]) al tasso naturale di $\varepsilon^{1/2}$ per problemi ellittici del secondo ordine in domini perforati.
   • Avanzamento Metodologico: A differenza dei lavori precedenti che richiedevano correttori di cella in $W^{1,\infty}$ e soluzioni macroscopiche in $C^1(0,T; W^{k,\infty})$ con $k \geq 4$, questo approccio utilizza l'operatore di scissione di scala $Q_\varepsilon$. Ciò riduce la regolarità richiesta della soluzione macroscopica a $H^2$ e mantiene i correttori di cella in $H^1_{\text{per}}$.
   • Consistenza Variazionale: Le stime di consistenza sono derivate interamente in forma variazionale, utilizzando la debole condizione di non flusso ereditata dal problema di cella. Ciò evita i contributi ai termini di bordo che precedentemente causavano la degradazione del tasso a $\varepsilon^{1/4}$.

Gli autori osservano che, sebbene il tasso $\varepsilon^{1/2}$ sia stretto all'interno dell'attuale framework variazionale a causa dello strato limite, ottenere il tasso ottimale $O(\varepsilon)$ su domini limitati richiederebbe probabilmente un argomento di dualità (come in Kenig–Lin–Shen [40]) o un approccio basato su stime di operatore (Birman–Suslina [7–9]), che sono lasciati come lavoro futuro. L'articolo non pretende di ottenere un tasso per il gradiente del potenziale chimico $\mu^\varepsilon$ a questo ordine, poiché l'attuale regolarità dei residui non giustifica i test necessari.