Global well-posedness of strong solutions to a bulk-surface Navier-Stokes-Cahn-Hilliard model with non-degenerate mobilities in two dimensions

基本信息

• 论文ID: 2511.06847
• 标题: Global well-posedness of strong solutions to a bulk-surface Navier-Stokes-Cahn-Hilliard model with non-degenerate mobilities in two dimensions
• 作者: Jonas Stange (Universität Regensburg, Germany)
• 分类: math.AP (Analysis of PDEs), math-ph (Mathematical Physics), math.MP (Mathematical Physics)
• 提交时间: 2025年11月10日
• 论文链接: https://arxiv.org/abs/2511.06847v1

摘要

本文研究了一个热力学一致的扩散界面模型，用于描述体-表面粘性流体混合物。该模型由体内的Navier-Stokes-Cahn-Hilliard方程与边界上的表面Navier-Stokes-Cahn-Hilliard系统耦合而成。作者在二维情况下证明了强解的全局适定性，涵盖了物理上有意义的非退化迁移率函数情形。此外，还证明了强解的唯一性及其对初始数据的连续依赖性。该方法基于带非退化迁移率的对流体-表面Cahn-Hilliard方程的新适定性和正则性理论，以及带非常数系数的体-表面Stokes方程。

研究背景与动机

问题背景

1. 物理现象建模: 研究涉及两相流体混合物的扩散界面模型，其中流体具有不同密度，且在体域和边界上都存在粘性动力学效应。这类模型在材料科学、生物膜动力学等领域具有重要应用。
2. 数学挑战: 该模型耦合了：
   • 体域中的Navier-Stokes-Cahn-Hilliard方程（描述流体运动和相分离）
   • 边界上的表面Navier-Stokes-Cahn-Hilliard系统（描述表面动力学）
   • 动态边界条件（允许体-表面之间的质量交换）
3. 现有理论的局限性:
   • 对于带动态边界条件的扩散界面模型，现有文献主要关注弱解的存在性
   • 强解的存在性研究较少，特别是对于非退化迁移率函数的情形
   • 三维情形下的全局强解理论尚不完善

研究动机

1. 物理合理性: 非退化迁移率函数（即迁移率在整个区域内有正的上下界）更符合物理实际
2. 数学严格性: 强解提供了更精细的正则性估计，对理解解的长时间行为至关重要
3. 理论完善: 填补了带动态边界条件的Navier-Stokes-Cahn-Hilliard系统强解理论的空白

核心贡献

本文的主要贡献包括：

1. 全局强解存在性: 在二维情况下，证明了带非退化迁移率的体-表面Navier-Stokes-Cahn-Hilliard系统全局强解的存在性（定理3.1）
2. 唯一性与连续依赖性: 证明了强解的唯一性以及对初始数据的连续依赖性（定理3.4），但需要限制条件 $L \in (0,\infty]$
3. 新的正则性理论: 建立了两个关键的辅助理论：
   • 带非退化迁移率的对流体-表面Cahn-Hilliard方程的适定性和正则性理论（附录A）
   • 带非常数系数的体-表面Stokes方程的正则性理论（附录B）
4. 能量不等式: 证明了解满足能量耗散律： $\frac{d}{dt}E_{tot}(v,w,\phi,\psi) + \text{耗散项} = 0$
5. 分离性质: 证明了相场函数严格保持在物理区间 $(-1,1)$ 内：$|\phi| < 1$ a.e. in $Q$, $|\psi| < 1$ a.e. on $\Sigma$

方法详解

任务定义

考虑二维有界区域 $\Omega \subset \mathbb{R}^2$ 及其边界 $\Gamma = \partial\Omega$。设 $Q = \Omega \times (0,\infty)$, $\Sigma = \Gamma \times (0,\infty)$。

控制方程系统（式1.1）包括：

1. 体域动量方程: $\partial_t(\rho(\phi)v) + \text{div}(v \otimes (\rho(\phi)v + J)) = \text{div} T, \quad \text{div} v = 0 \text{ in } Q$
2. 表面动量方程: $\partial_t(\sigma(\psi)w) + \text{div}_\Gamma(w \otimes (\sigma(\psi)w + K)) = \text{div}_\Gamma T_\Gamma + Z, \quad \text{div}_\Gamma w = 0 \text{ on } \Sigma$
3. 体域Cahn-Hilliard方程: $\partial_t\phi + \text{div}(\phi v) = \text{div}(m_\Omega(\phi)\nabla\mu), \quad \mu = -\Delta\phi + F'(\phi) \text{ in } Q$
4. 表面Cahn-Hilliard方程: $\partial_t\psi + \text{div}_\Gamma(\psi w) = \text{div}_\Gamma(m_\Gamma(\psi)\nabla_\Gamma\theta) - \beta m_\Omega(\phi)\partial_n\mu \text{ on } \Sigma$$\theta = -\Delta_\Gamma\psi + G'(\psi) + \alpha\partial_n\phi \text{ on } \Sigma$
5. 边界条件: $w = v|_\Gamma, \quad v \cdot n = 0 \text{ on } \Sigma$$K\partial_n\phi = \alpha\psi - \phi, \quad Lm_\Omega(\phi)\partial_n\mu = \beta\theta - \mu \text{ on } \Sigma$

其中：

• $v, w$ 为体域和表面速度场
• $\phi, \psi$ 为体域和表面相场变量
• $\mu, \theta$ 为化学势
• $\rho(\phi), \sigma(\psi)$ 为依赖于相场的密度函数
• $m_\Omega, m_\Gamma$ 为迁移率函数（满足 $0 < m_* \leq m_\Omega, m_\Gamma \leq m^*$）

模型架构

证明策略采用半Galerkin逼近方法：

第一步：近似问题构造

1. 引入体-表面Stokes算子 $A$ 的特征函数系 $\{(\tilde{v}_j, \tilde{w}_j)\}_{j\in\mathbb{N}}$
2. 定义有限维子空间 $V_k = \text{span}\{(\tilde{v}_1, \tilde{w}_1), \ldots, (\tilde{v}_k, \tilde{w}_k)\}$
3. 构造近似解 $(v_k, w_k, \phi_k, \psi_k, \mu_k, \theta_k)$ 满足Galerkin形式的弱方程

第二步：不动点论证

对固定的速度场 $(v_*, w_*)$，求解：

• 对流体-表面Cahn-Hilliard方程（4.8），得到 $(\phi_k, \psi_k, \mu_k, \theta_k)$
• Galerkin逼近的Navier-Stokes方程（4.12），得到 $(v_k, w_k)$

通过Schauder不动点定理证明映射 $\Lambda: (v_*, w_*) \mapsto (v_k, w_k)$ 存在不动点。

第三步：先验估计

建立关键的一致估计（独立于 $k$）：

1. 能量估计（式4.35-4.37）: $\|(v_k, w_k)\|_{L^\infty(0,\infty; L^2)} + \|(v_k, w_k)\|_{L^2(0,\infty; H^1)} \leq C$$\|(\phi_k, \psi_k)\|_{L^\infty(0,\infty; H^1)} \leq C$
2. 高阶估计（式4.43-4.45）: $\|(\phi_k, \psi_k)\|_{L^\infty(0,\infty; W^{2,p})} + \|(F'(\phi_k), G'(\psi_k))\|_{L^\infty(0,\infty; L^p)} \leq K_{1,p}$$\|(\mu_k, \theta_k)\|_{L^\infty(0,\infty; H^1)} + \|(\mu_k, \theta_k)\|_{L^2_{uloc}(0,\infty; H^3)} \leq K_2$
3. 速度场高阶估计（式4.61）: $\sup_{t\geq 0} \int_t^{t+1} \|( v_k, w_k)\|_{H^2}^2 + \|(\partial_t v_k, \partial_t w_k)\|_{L^2}^2 + \|(p_k, q_k)\|_{H^1}^2 \, ds \leq C$

第四步：极限过程

利用紧性论证（Aubin-Lions引理）和单调性技巧，在 $k \to \infty$ 时提取收敛子列，得到强解。

技术创新点

1. 处理非退化迁移率: 与常数迁移率情形不同，非退化迁移率引入额外的非线性项。关键技巧是利用椭圆正则性理论（式2.13-2.15）控制这些项。
2. 高阶能量估计: 通过微分能量泛函（式4.46），结合体-表面Stokes正则性理论（定理B.2），建立速度场的 $H^2$ 估计： $\frac{d}{dt}\text{(能量)} + \text{耗散} \leq \varpi\|(v_k, w_k)\|_{H^2}^2 + \text{低阶项}$ 选择适当的 $\varpi$，吸收 $H^2$ 范数。
3. 分离性质的保持: 利用对数势的奇异性和最大值原理的变体，证明相场函数严格保持在 $(-1, 1)$ 内。
4. 唯一性证明: 对两个解的差 $(v, w, \phi, \psi) = (v_1 - v_2, w_1 - w_2, \phi_1 - \phi_2, \psi_1 - \psi_2)$，建立微分不等式（式5.2-5.3）： $\frac{d}{dt}\left(\|(v, w)\|_{L^2}^2 + \|(\phi, \psi)\|_{K,\alpha}^2\right) + \text{耗散} \leq P(t) \times \text{(差的范数)}$ 其中 $P(t) \in L^1(0,T)$，应用Gronwall引理得到唯一性。

实验设置

注意: 本文为纯理论数学论文，不包含数值实验。主要结果是存在性和唯一性的数学证明。

理论验证框架

文章通过以下数学工具验证理论：

1. 函数空间设置:
   • 速度场空间: $H^1_{0,\text{div}} = \{(v,w) \in H^1_0: \text{div} v = 0, \text{div}_\Gamma w = 0\}$
   • 相场空间: $V^1_{L,\beta}$ (依赖于耦合参数 $L, \beta$)
   • 化学势空间: $H^1_{L,\beta}$
2. 正则性要求:
   • 初始数据: $(v_0, w_0) \in H^1_{\text{div}}$, $(\phi_0, \psi_0) \in H^1$ 满足 $\|\phi_0\|_{L^\infty} \leq 1$, $\|\psi_0\|_{L^\infty} \leq 1$
   • 兼容性条件: 存在 $(\mu_0, \theta_0) \in H^1_{L,\beta}$ 满足特定积分关系
3. 参数限制:
   • 空间维度: $d = 2$
   • 耦合参数: $K \in (0, \infty)$, $L \in [0, \infty]$（唯一性需要 $L \in (0, \infty]$）
   • 迁移率: $0 < m_* \leq m_\Omega, m_\Gamma \leq m^* < \infty$

关键假设

(A1-A7) 包括：

• 区域正则性: $\Omega \subset \mathbb{R}^2$ 为 $C^3$ 边界
• 势函数: $F(s) = F_0(s) - \frac{c_F}{2}s^2$, 其中 $F_0$ 强凸且在 $s \to \pm 1$ 时奇异
• 增长条件: $F_0''(s) \geq \Theta > 0$ 或特定增长率限制（(A7.1)-(A7.2)）

实验结果

主要理论结果

定理3.1（全局强解存在性）: 在假设 (A1)-(A7) 下，对 $K \in (0, \infty)$, $L \in [0, \infty]$，存在全局强解 $(v, w, \phi, \psi, \mu, \theta)$ 满足：

1. 正则性:
   • $(v, w) \in BC([0,\infty); H^1_{0,\text{div}}) \cap L^2_{uloc}([0,\infty); H^2) \cap H^1_{uloc}(0,\infty; L^2_{\text{div}})$
   • $(\phi, \psi) \in L^\infty(0,\infty; H^3)$
   • $(F'(\phi), G'(\psi)), (F''(\phi), G''(\psi)) \in L^\infty(0,\infty; L^p)$ 对所有 $2 \leq p < \infty$
2. 分离性质: $|\phi| < 1$ a.e. in $Q$, $|\psi| < 1$ a.e. on $\Sigma$
3. 能量不等式: $E_{tot}(v(t), w(t), \phi(t), \psi(t)) + \int_0^t \text{(耗散)} \, ds = E_{tot}(v_0, w_0, \phi_0, \psi_0)$

定理3.4（唯一性）: 在定理3.1的假设下，若 $L \in (0, \infty]$，则强解唯一。

关键估计

文章建立了以下定量估计（式A.10-A.12）：

$\|(\phi, \psi)\|_{L^\infty(0,\infty; H^1)}^2 + \|(\partial_t\phi, \partial_t\psi)\|_{L^2(0,\infty; (H^1_{L,\beta})')}^2 \leq C(1 + \|(v,w)\|_{L^2(0,\infty; L^2)}^2)$

$\|(\phi, \psi)\|_{L^2_{uloc}([0,\infty); W^{2,p})}^2 + \|(F'(\phi), G'(\psi))\|_{L^2_{uloc}([0,\infty); L^p)}^2 \leq C_p(1 + \|(v,w)\|_{L^2(0,\infty; L^2)}^2)$

$\|(\phi, \psi)\|_{L^4_{uloc}([0,\infty); H^2)}^4 \leq C(1 + \|(v,w)\|_{L^2(0,\infty; L^2)}^2)^2$

技术难点克服

1. $L = 0$ 情形的限制: 当 $L = 0$ 时，边界条件 (1.1h)$_2$ 不允许重新表述法向导数 $\partial_n\mu$，导致唯一性证明失效。作者指出这需要额外的结构假设（备注3.5）。
2. 三维情形的挑战: 文章仅处理二维情形，因为三维中：
   • 失去分离性质的传播
   • 需要更精细的逼近技术控制非线性项
   • 只能期望局部强解而非全局强解

相关工作

主要研究方向

1. 无动态边界条件的模型:
   • Abels-Garcke-Grün (AGG) 模型：2, 4, 16, 17, 21, 22
   • 已有成熟的强解理论
2. 带动态边界条件的模型:
   • 主要关注弱解存在性：9, 14, 15, 19, 30
   • 准强解：Allen-Cahn-Navier-Stokes-Voigt系统 13
3. Cahn-Hilliard方程的正则性:
   • 常数迁移率情形：20, 28, 29
   • 非退化迁移率：11, 39（本文作者的前期工作）

本文的创新性

相比现有工作，本文：

• 首次处理带动态边界条件的Navier-Stokes-Cahn-Hilliard系统的强解
• 首次涵盖非退化迁移率的情形
• 建立了新的体-表面Stokes方程正则性理论（附录B）
• 发展了对流Cahn-Hilliard方程的高阶估计技术（附录A）

结论与讨论

主要结论

1. 在二维情况下，证明了带非退化迁移率的体-表面Navier-Stokes-Cahn-Hilliard系统全局强解的存在性和唯一性（当 $L \in (0, \infty]$）
2. 解满足严格的分离性质和能量耗散律
3. 建立了两个重要的辅助理论框架，可独立应用于其他问题

局限性

1. 维度限制: 结果仅适用于二维情形。三维情形需要：
   • 更精细的非线性估计
   • 可能只能得到局部强解
   • 对常数迁移率情形，弱解唯一性和强解存在性已知
2. 参数限制:
   • $K = 0$ 不可行（非退化迁移率情形）
   • $L = 0$ 的唯一性未解决，除非有额外结构假设
3. 兼容性条件: 当 $L = 0$ 时，需要密度函数满足特殊关系 $\beta(\sigma_2 - \sigma_1) = \rho_2 - \rho_1$（式3.2）

未来方向

1. 三维推广: 发展适用于三维情形的技术，可能需要：
   • 局部存在性理论
   • 小初值或特殊结构下的全局性
   • 常数迁移率假设
2. $L = 0$ 情形: 寻找处理该情形的新方法或确定其不可行性
3. 长时间行为: 研究解的渐近行为、吸引子存在性
4. 数值验证: 发展高效数值格式验证理论预测

深度评价

优点

1. 理论严格性:
   • 证明完整细致，技术处理规范
   • 对各种技术难点给出清晰说明
   • 引理和定理陈述精确
2. 创新性:
   • 首次建立该类系统的强解理论
   • 非退化迁移率的处理具有挑战性
   • 附录中的辅助理论具有独立价值
3. 结构清晰:
   • 主要结果在第3节清晰陈述
   • 证明策略层次分明
   • 技术细节放在附录，不影响主线阅读
4. 物理合理性:
   • 模型热力学一致
   • 非退化迁移率更符合实际
   • 能量耗散律和质量守恒得到保证

不足

1. 维度限制: 二维结果的实际应用受限，三维情形更为重要但未解决
2. 参数限制: $K \in (0, \infty)$ 和 $L \in (0, \infty]$ 的限制排除了某些物理情形
3. 缺少数值验证: 作为纯理论工作，未提供数值算例支持理论结果
4. 技术复杂性: 证明高度技术化，需要大量预备知识，可能限制读者群
5. 常数依赖性: 许多估计中的常数 $C$ 依赖于多个参数，未给出显式表达式

影响力

1. 理论贡献:
   • 填补了带动态边界条件的Navier-Stokes-Cahn-Hilliard系统强解理论的重要空白
   • 为后续研究提供了技术工具和理论框架
2. 方法论价值:
   • 半Galerkin方法与不动点论证的结合可应用于其他耦合系统
   • 高阶估计技术具有一般性
3. 应用潜力:
   • 为数值方法提供理论支撑
   • 可应用于材料科学、生物膜动力学等领域的建模

适用场景

1. 材料科学: 两相流体材料的相分离和界面动力学
2. 生物物理: 细胞膜的流体镶嵌模型，脂质双层动力学
3. 接触线动力学: 固体表面上液滴的运动和润湿现象
4. 理论研究:
   • 非线性偏微分方程的正则性理论
   • 耦合系统的适定性分析
   • 动态边界条件的数学理论

参考文献（精选）

2 H. Abels, H. Garcke, and A. Giorgini. Global regularity and asymptotic stabilization for the incompressible Navier-Stokes-Cahn-Hilliard model with unmatched densities. Math. Ann., 2024.

11 M. Conti, P. Galimberti, S. Gatti, and A. Giorgini. New results for the Cahn-Hilliard equation with non-degenerate mobility. Calc. Var. PDE, 2025.

19 A. Giorgini and P. Knopf. Two-phase flows with bulk-surface interaction: thermodynamically consistent Navier-Stokes-Cahn-Hilliard models. J. Math. Fluid Mech., 2023.

30 P. Knopf and J. Stange. A thermodynamically consistent model for bulk-surface viscous fluid mixtures. Preprint, 2025.

39 J. Stange. Well-posedness and long-time behavior of a bulk-surface Cahn-Hilliard model with non-degenerate mobility. Preprint, 2025.

总体评价: 这是一篇高质量的理论数学论文，在偏微分方程分析领域做出了重要贡献。尽管存在维度和参数限制，但其技术创新和理论深度值得肯定。对于从事相关领域研究的数学家和应用科学家，本文提供了宝贵的理论工具和研究范式。