[论文解读] A $L^\infty$ bound for the Cahn--Hilliard equation with relaxed non-smooth free energy

Christian Kahle|arXiv (Cornell University)|Nov 9, 2015

Solidification and crystal growth phenomena参考文献 1被引用 3

一句话总结

该论文通过使用 Moreau–Yosida 正则化处理具有松弛双势阱自由能的 Cahn–Hilliard 方程，建立了对物理约束 $|\varphi| \leq 1$ 违反情况的 $L^\infty$ 有界性。结合最优控制理论中的分析技术与 Hölder 正则性估计，作者推导出约束违反程度关于松弛参数的收敛速率，确保相场在多相流模型中保持物理一致性所必需的统一有界性。

ABSTRACT

Phase field models are widely used to describe multiphase systems. Here a smooth indicator function, called phase field, is used to describe the spatial distribution of the phases under investigation. Material properties like density or viscosity are introduced as given functions of the phase field. These parameters typically have physical bounds to fulfil, e.g. positivity of the density. To guarantee these properties, uniform bounds on the phase field are of interest. In this work we derive a uniform bound on the solution of the Cahn--Hilliard system, where we use the double-obstacle free energy, that is relaxed by Moreau--Yosida relaxation.

研究动机与目标

• 建立具有松弛双势阱自由能的 Cahn–Hilliard 方程解在 $L^\infty(\Omega)$ 范围内的统一有界性，以防止相场出现非物理取值。
• 通过保证密度和粘度等材料属性为正，确保相场模型的物理一致性。
• 将最优控制理论中的收敛速率结果推广至具有 Neumann 边界条件的 Cahn–Hilliard 系统。
• 在二维和三维中通过自适应有限元方法与牛顿法对理论有界性进行数值验证。
• 比较松弛系统中罚项的不同数值处理方式（精确积分、质量集中法、插值法）对收敛性与约束违反的影响。

提出的方法

• 构建具有双势阱自由能的 Cahn–Hilliard 系统，并应用 Moreau–Yosida 正则化以处理非光滑势能。
• 将松弛系统解释为 $H^1(\Omega)$ 空间中带约束的最优控制问题的一阶最优性条件。
• 借鉴 [HSW14] 中关于 $C(\bar{\Omega})$-约束最优控制的技术，推导出约束违反的 $L^1(\Omega)$ 有界性，进而获得 $L^\infty(\Omega)$ 有界性。
• 对解施加结构假设，以在初始估计基础上进一步改善 $L^\infty$ 有界性。
• 采用基于残差的误差估计器与 Dörffer 标记策略的自适应有限元方法进行网格加密。
• 在 FEniCS 中使用 PETSc 与 MUMPS 求解器，通过牛顿法求解非线性系统，并对罚项采用不同的积分方案（精确积分、质量集中法、插值法）。

实验结果

研究问题

• RQ1 能否为具有松弛双势阱自由能的 Cahn–Hilliard 方程相场解建立统一的 $L^\infty(\Omega)$ 有界性？
• RQ2 约束违反程度 $\|\lambda(\varphi^s)\|_{L^\infty(\Omega)}$ 随 Moreau–Yosida 松弛参数 $s$ 变化的收敛速率如何？
• RQ3 不同的数值积分方案（精确积分、质量集中法、插值法）在实际中如何影响 $L^\infty$ 有界性的收敛性与精度？
• RQ4 理论 $L^\infty$ 有界性在自适应网格加密与牛顿法求解非线性系统条件下是否依然成立？
• RQ5 所推导的有界性能否在二维与三维空间中，针对不同罚项形式通过数值方法验证？

主要发现

• 论文建立了对物理约束 $|\varphi| \leq 1$ 违反情况的理论 $L^\infty(\Omega)$ 有界性，确保相场始终处于物理上有意义的范围内。
• 随着松弛参数 $s$ 增大，该有界性以 $\mathcal{O}(s^{-1/2})$ 的速率衰减，与最优控制理论结果一致。
• 数值实验表明，理论有界性在所有测试的二维与三维配置中均被准确实现。
• 当罚项 $\lambda(\varphi^s)$ 采用插值处理时，牛顿法在 $s$ 较大时无法收敛，而精确积分与质量集中法则保持鲁棒性。
• 对于更高阶罚项 $\lambda_k$（$k>2$），需增大 $s$ 才能达到相近的约束违反水平，但理论有界性依然成立。
• 自适应有限元策略能有效解析界面层，并在不同空间维度与松弛参数下保持计算精度。