Navier-Stokes-Cahn-Hilliard system in a $3$D perforated domain with free slip and source term: Existence and homogenization

这是一篇关于多孔介质中 Navier-Stokes-Cahn-Hilliard (NSCH) 系统的数学分析论文。文章主要研究了在三维周期性穿孔域中，具有可变粘度、非保守源项以及混合边界条件（外边界无滑移，固体夹杂表面自由滑移）的扩散界面模型的弱解存在性及其周期均质化（Homogenization）问题。

以下是对该论文的详细技术总结：

1. 研究问题与背景 (Problem Statement)

• 物理背景：研究多孔介质或微结构环境中的多相输运过程（如地下水修复、增强石油采收率、生物组织支架等）。在孔隙尺度上，流动由粘性输运、毛细力、界面能和相分离机制的复杂相互作用控制。
• 数学模型：
  • 考虑一个二元不可压缩混合物，由时间依赖的 NSCH 系统描述。
  • 动量方程：包含非线性对流项、可变粘度张量 $A^\varepsilon$、Korteweg 型毛细力项 $\lambda_\varepsilon \phi_\varepsilon \nabla \mu_\varepsilon$ 以及外力项。
  • 相场方程：Cahn-Hilliard 方程，包含对流项 $u_\varepsilon \cdot \nabla \phi_\varepsilon$、非保守源项 $G(\phi_\varepsilon)$（破坏质量守恒）以及化学势 $\mu_\varepsilon$ 的定义。
  • 边界条件：
    • 外边界 $\partial \Omega$：无滑移条件 ($u_\varepsilon = 0$)。
    • 固体夹杂表面 $\Gamma^\varepsilon_s$：自由滑移条件（Navier slip with zero tangential stress），即法向速度为零 ($u_\varepsilon \cdot n_\varepsilon = 0$)，切向应力为零 ($(A^\varepsilon D(u_\varepsilon)n_\varepsilon) \cdot \tau^\varepsilon = 0$)。
  • 尺度参数：$\varepsilon$ 为微观周期尺度，$\lambda_\varepsilon$ 为毛细强度，满足 $\lim_{\varepsilon \to 0} \lambda_\varepsilon = \lambda \in [0, +\infty)$。

2. 方法论 (Methodology)

论文的分析分为两个主要部分：

第一部分：弱解存在性与先验估计 (Existence and A Priori Estimates)

• 函数空间：定义了适应穿孔几何的 Sobolev 空间，特别是处理了满足混合边界条件的速度空间。
• 关键工具：
  • Bogovski˘ı 算子：在 $\varepsilon$ 依赖的流体域 $\Omega^\varepsilon_p$ 和固体域 $\Omega^\varepsilon_s$ 上构造了具有与 $\varepsilon$ 无关范数界的逆散度算子，用于处理压力项和构造测试函数。
  • 延拓算子：构造了从穿孔域到整体域的延拓算子（包括对偶空间），使得可以在宏观域上应用紧性论证。
  • Korn 不等式与 Poincaré 不等式：建立了与 $\varepsilon$ 无关的均匀估计。
• 存在性证明：采用 Galerkin 逼近法。
  • 利用适应 Navier 边界条件的 Stokes 型特征函数作为速度基。
  • 利用有限维 ODE 系统近似，结合 Gronwall 不等式推导先验估计，证明解在有限时间 $(0, T)$ 内存在且无爆破。
  • 利用紧性论证（Aubin-Lions 引理、Simon 紧性准则）取极限得到弱解。
• 能量估计：推导了关于 $\varepsilon$ 和 $\lambda_\varepsilon$ 一致的能量耗散律。由于源项 $G(\phi_\varepsilon)$ 的存在，质量不守恒，导致浓度均值随时间变化，因此需要特殊的处理来估计化学势和浓度的均值。

第二部分：均质化 (Homogenization)

• 工具：使用周期性展开算子 (Periodic Unfolding Operator) $T^\varepsilon$ 和 $T^*_\varepsilon$，这是处理穿孔域均质化的核心工具，能够将微观振荡转化为双尺度极限。
• 极限过程：令 $\varepsilon \to 0$，根据毛细强度极限 $\lambda = \lim \lambda_\varepsilon$ 的不同，分两种情况讨论：
  1. $\lambda = 0$ (消失毛细力 regime)：外力 $g_\varepsilon$ 的缩放使得惯性项消失。
  2. $\lambda \in (0, +\infty)$ (平衡 regime)：外力缩放使得惯性项保留。
• 细胞问题 (Cell Problems)：通过求解定义在参考单元 $Y_p$ 上的细胞问题，将微观修正项（correctors）表示为宏观变量的函数，从而导出有效张量。

3. 主要结果 (Key Results)

A. 弱解存在性

证明了对于任意固定的 $\varepsilon > 0$，在有限时间 $(0, T)$ 内，NSCH 系统存在弱解 $(u_\varepsilon, \phi_\varepsilon, \mu_\varepsilon)$，并建立了关于 $\varepsilon$ 和 $\lambda_\varepsilon$ 一致的先验估计。

B. 均质化极限模型

根据 $\lambda$ 的值，得到了两个截然不同的宏观有效模型：

1. 情形 1：$\lambda = 0$ (Stokes-Cahn-Hilliard 型)

   • 物理意义：毛细力太弱，无法在宏观尺度上产生对流。
   • 方程特征：
     • 动量方程退化为 Stokes 方程，没有宏观对流项 $(u \cdot \nabla)u$。
     • 相场方程没有宏观平流项 $u \cdot \nabla \phi$。
     • 系统由有效粘度张量 $A^{hom}$、有效毛细耦合张量 $C^{hom}$ 和有效扩散张量 $B^{hom}$ 描述。
   • 形式：
     $$\begin{cases} \partial_t u - \nabla \cdot (A^{hom} D(u)) + \nabla \tilde{p} + \phi C^{hom} \nabla \mu = g \\ \nabla \cdot u = 0 \\ \partial_t \phi + G(\phi) - \nabla \cdot (C^{hom} \nabla \mu) = 0 \\ \mu - F'(\phi) + \nabla \cdot (B^{hom} \nabla \phi) = 0 \end{cases}$$

2. 情形 2：$\lambda \in (0, +\infty)$ (Navier-Stokes-Cahn-Hilliard 型)

   • 物理意义：毛细力与粘性力平衡，界面运动诱导显著的宏观流动。
   • 方程特征：
     • 动量方程保留 Navier-Stokes 形式，包含非线性对流项，但系数为 $\sqrt{\lambda}$。
     • 相场方程保留 平流项 $u \cdot \nabla \phi$，系数同样为 $\sqrt{\lambda}$。
     • 有效张量 $A^{hom}, B^{hom}, C^{hom}$ 与 $\lambda=0$ 情形相同（由细胞问题定义）。
   • 形式：
     $$\begin{cases} \partial_t u + \sqrt{\lambda}(u \cdot \nabla)u - \nabla \cdot (A^{hom} D(u)) + \nabla \tilde{p} + \phi C^{hom} \nabla \mu = g \\ \nabla \cdot u = 0 \\ \partial_t \phi + \sqrt{\lambda} u \cdot \nabla \phi + G(\phi) - \nabla \cdot (C^{hom} \nabla \mu) = 0 \\ \mu - F'(\phi) + \nabla \cdot (B^{hom} \nabla \phi) = 0 \end{cases}$$

C. 能量收敛

证明了微观总能量 $T_\varepsilon$ 在适当缩放后收敛到宏观总能量 $T$，且宏观能量满足与微观系统类似的耗散律。

4. 创新点与贡献 (Contributions & Significance)

1. 混合边界条件的处理：不同于以往文献中常见的无滑移条件，本文处理了固体夹杂表面的**自由滑移（零切向应力）**条件。这在多孔介质流动中更符合物理实际（如超疏水表面或特定孔隙结构），并增加了数学分析的复杂性（需要构造特殊的 Bogovski˘ı 算子和延拓算子）。
2. 非保守源项的引入：模型包含了非保守源项 $G(\phi)$，导致质量不守恒。这使得传统的 Cahn-Hilliard 理论（通常假设质量守恒）无法直接应用，作者通过精细的能量估计和均值控制克服了这一困难。
3. 统一的均质化框架：文章在一个统一的框架下，根据毛细强度参数 $\lambda$ 的极限行为，严格推导出了两种截然不同的宏观模型（Stokes 型 vs. Navier-Stokes 型）。这揭示了毛细力在多孔介质多相流宏观输运中的关键调节作用。
4. 严格性：与以往主要依赖启发式渐近展开（Asymptotic Expansion）的研究不同，本文提供了基于变分方法和紧性论证的严格数学证明（Rigorous Homogenization）。
5. 变量粘度与穿孔域：结合了可变粘度张量和复杂的穿孔几何，扩展了扩散界面模型在多孔介质中的均质化理论。

5. 总结

该论文在数学流体力学和多孔介质理论领域做出了重要贡献。它不仅解决了具有复杂边界条件和非保守源项的 NSCH 系统的适定性问题，还通过严格的均质化分析，揭示了微观毛细力强度如何决定宏观流动是处于粘性主导（Stokes）还是惯性主导（Navier-Stokes）状态。这些结果为理解多孔介质中的多相流、相分离及界面动力学提供了坚实的理论基础。