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

技术摘要：周期性穿孔域中带有源项的 Cahn–Hilliard 系统的定量同质化

问题陈述
本文研究了在周期性穿孔域中，具有非保守、单调源项的 Cahn–Hilliard 系统的定性与定量同质化。微观模型描述了多孔介质中的相分离过程，其中流体占据着一个连通的孔隙空间 $\Omega^\varepsilon_p \subset \Omega \subset \mathbb{R}^d$（$d \in \{2,3\}$），该空间是通过从有界 Lipschitz 区域中移除周期性阵列的固体夹杂物而获得的。该系统通过以下方程将序参数 $\phi^\varepsilon$ 与化学势 $\mu^\varepsilon$ 耦合：
$$\begin{cases} \partial_t \phi^\varepsilon - \Delta \mu^\varepsilon + G(\phi^\varepsilon) = 0, \\ \mu^\varepsilon = -\Delta \phi^\varepsilon + F'(\phi^\varepsilon), \end{cases}$$
并满足外边界 $\partial \Omega$ 和孔壁 $\Gamma^\varepsilon_s$ 上的齐次 Neumann 边界条件。势函数 $F$ 满足多项式增长和强制性（例如经典的双阱势），而源项 $G$ 是全局 Lipschitz 连续、单调且满足 $G(0)=0$ 的。主要的挑战在于当孔隙尺度 $\varepsilon \to 0$ 时推导有效宏观方程，并建立微观解与同质化解之间的误差锐利收敛率。

方法论
分析采用了适用于穿孔域的周期性展开法（periodic unfolding method），利用展开算子 $T^\varepsilon_*$ 来处理振荡几何结构。该方法分两个主要阶段进行：

1. 定性同质化：

   • 能量估计： 推导了微观解的一致先验估计。一个关键的结构性观察是，单调源项 $G$（且 $G(0)=0$）在能量恒等式中贡献了一个耗散项。具体而言，将化学势方程与 $G(\phi^\varepsilon)$ 进行测试会产生一个正的梯度项，这强化了耗散而非与之竞争。这使得无需使用带有指数时间因子的 Gronwall 论证即可获得独立于 $\varepsilon$ 的一致界限。
   • 两尺度收敛： 利用展开算子和紧致性论证（Aubin–Lions–Simon 引理），作者证明了 $(\phi^\varepsilon, \mu^\varepsilon)$ 向同质化对 $(\phi, \mu)$ 的收敛性。极限系统是宏观域 $\Omega$ 上的 Cahn–Hilliard 方程，具有有效扩散张量 $B_{\text{hom}}$。
   • 单元问题： 有效张量 $B_{\text{hom}}$ 通过在参考孔单元 $Y_p$ 上的标量 Neumann 单元问题来表征。修正子 $\chi_i$ 属于 $H^1_{\text{per},0}(Y_p)$。

2. 定量分析（修正子估计）：

   • 尺度分裂算子： 为了实现定量速率，本文利用 $Q_1$ 尺度分裂算子 $Q_\varepsilon$（来自 Griso [29, 30]）构建了一阶修正近似 $\Phi^\varepsilon$ 和 $M^\varepsilon$。该算子作用于固定的宏观域 $\Omega$，随后被限制在 $\Omega^\varepsilon_p$ 上。至关重要的是，这种方法允许单元修正子保持在 $H^1_{\text{per}}(Y_p)$ 中，避免了经典形式展开所要求的严格 $W^{1,\infty}$ 假设。
   • 变分残差法： 误差分析在负范数框架下进行。序参数误差 $e^\varepsilon_\phi = \phi^\varepsilon - \Phi^\varepsilon$ 与 Neumann 反向拉普拉斯算子 $N_\varepsilon e^\varepsilon_\phi$ 进行测试，而化学势误差则与 $e^\varepsilon_\phi$ 进行测试。这避免了对残差更高正则性的需求。
   • 一致性估计： 证明依赖于推导考虑了靠近 $\partial \Omega$ 的不完整单元边界层效应的变分残差估计（引理 4.6）。

主要结果

• 定性收敛（定理 3.7）： 微观解（经子序列）收敛于同质化系统：
  $$\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}$$
  具有有效的 Neumann 边界条件。张量 $B_{\text{hom}}$ 是对称且一致椭圆的。

• 定量速率（定理 4.7）： 在对同质化解具有附加 $H^2$ 正则性假设以及准备充分的初始数据条件下，本文建立了 $\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^2(0,T; L^2(\Omega^\varepsilon_p))$ 中受控。
  • 修正后的梯度误差在 $L^2((0,T)\times\Omega^\varepsilon_p)$ 中受控。
  • 化学势误差在对偶范数 $L^2(0,T; H^1(\Omega^\varepsilon_p)')$ 中受控。

• 边界层效应： 有界域上 $\varepsilon^{1/2}$ 的速率归因于靠近 $\partial \Omega$ 的不完整单元的边界层。在平坦环面 $\mathbb{T}^d$ 上，不存在此类边界层，速率提升至 $O(\varepsilon)$（备注 4.10）。

意义与贡献
本文对穿孔介质中相场方程的同质化理论做出了两项主要贡献：

1. 定性扩展： 它严谨地推导了带有一般单调非保守源项的同质化 Cahn–Hildard 系统，证明了源项的耗散结构简化了能量估计，相比于一般的非保守强迫项。
2. 定量改进： 它将穿孔域中 Cahn–Hilliard 系统的已知收敛速率从 $\varepsilon^{1/4}$（由 Schmuck–Kalliadasis [49] 建立）提高到了穿孔域中二阶椭圆问题的自然速率 $\varepsilon^{1/2}$。
   • 方法论进步： 不同于以往要求单元修正子在 $W^{1,\infty}$ 中且宏观解在 $C^1(0,T; W^{k,\infty})$ 中（$k \geq 4$）的工作，本方法利用了尺度分裂算子 $Q_\varepsilon$。这降低了对宏观解正则性的要求至 $H^2$，并使单元修正子保持在 $H^1_{\text{per}}$ 中。
   • 变分一致性： 一致性估计完全以变分形式推导，利用了从单元问题中继承的弱无通量条件。这避免了此前导致速率退化至 $\varepsilon^{1/4}$ 的边界项贡献。

作者指出，由于边界层的存在，在当前的变分框架内 $\varepsilon^{1/2}$ 的速率是锐利的；若要在有界域上达到最优的 $O(\varepsilon)$ 速率，可能需要使用对偶论证（如 Kenig–Lin–Shen [40]）或算子估计方法（Birman–Suslina [7–9]），这些工作已留作未来研究。本文并未声称在当前阶数下获得了化学势 $\mu^\varepsilon$ 的梯度速率，因为当前残差的正则性不足以支撑所需的测试。