[论文解读] Consistency result for a non monotone scheme for anisotropic mean curvature flow

Élie Bretin, Éric Bonnetier+1 more|arXiv (Cornell University)|May 26, 2010

Geometric Analysis and Curvature Flows参考文献 37被引用 28

一句话总结

本文提出了一种在 ℝᵈ 中各向异性平均曲率流的新型非单调数值格式，通过在傅里叶空间中线性化各向异性 Allen-Cahn 方程的非线性扩散项，得到一个与核 $ K_{\rho,t}(x) = \mathcal{F}^{-1}[e^{-4\pi^2 t \phi^o(\xi)}] $ 的卷积。尽管该核不具备正性且其二阶矩不可积，该格式仍被证明与各向异性平均曲率流一致，在适当条件下建立了收敛性。

ABSTRACT

In this paper, we propose a new scheme for anisotropic motion by mean curvature in $\\R^d$. The scheme consists of a phase-field approximation of the motion, where the nonlinear diffusive terms in the corresponding anisotropic Allen-Cahn equation are linearized in the Fourier space. In real space, this corresponds to the convolution with a kernel of the form \\[ K_{\\phi,t}(x) = \\F^{-1}\\left[ e^{-4\\pi^2 t \\phi^o(\\xi)} \ ight](x). \\] We analyse the resulting scheme, following the work of Ishii-Pires-Souganidis on the convergence of the Bence-Merriman-Osher algorithm for isotropic motion by mean curvature. The main difficulty here, is that the kernel $K_{\\phi,t}$ is not positive and that its moments of order 2 are not in $L^1(\\R^d)$. Still, we can show that in one sense the scheme is consistent with the anisotropic mean curvature flow.

研究动机与目标

• 开发一种在 ℝᵈ 中各向异性平均曲率流的稳定、高效数值格式，该格式将各向同性平均曲率运动推广至各向异性的设定。
• 通过引入相场近似并在线性化非线性扩散项时采用傅里叶空间方法，克服现有格式的局限性。
• 在粘性解框架与相场近似下，分析一种非单调格式（其核非正且其二阶矩不属于 L¹）与各向异性平均曲率流的一致性。
• 通过二维和三维模拟对格式进行数值验证，展示其收敛于 Wulff 集并准确模拟界面演化。

提出的方法

• 该格式源自各向异性 Allen-Cahn 方程，通过傅里叶空间中的核 $ K_{\rho,t}(x) = \mathcal{F}^{-1}[e^{-4\pi^2 t \phi^o(\xi)}] $ 对非线性扩散项进行线性化，从而在实空间中形成基于卷积的演化。
• 该方法采用分裂策略，将 Allen-Cahn 演化分解为扩散步（通过傅里叶变换）与反应步（通过势能 $ W(s) = \frac{1}{2}s^2(1-s)^2 $），从而实现无条件稳定的时间积分。
• 分析基于 Ishii-Pires-Souganidis 的粘性解框架，经适配以适用于非单调核，从而建立与各向异性平均曲率流的一致性。
• 该格式采用谱方法实现，使用 $ P = 2^8 $ 或 $ 2^7 $ 个傅里叶模态，时间步长 $ \delta_t = 1/P^2 $，且令 $ \epsilon = 1/P $ 以控制扩散界面的宽度。
• 初始条件采用符号距离函数 $ q $，其通过最小化一维能量泛函获得，以确保初始界面清晰。
• 数值验证包括从 Wulff 集和在体积约束下从圆形初始形状出发的演化，通过特征函数的 $ L^1 $ 误差衡量收敛性。

实验结果

研究问题

• RQ1 即使核不具备正性且其二阶矩不属于 $ L^1(\mathbb{R}^d) $，基于傅里叶空间线性化的非单调格式能否与各向异性平均曲率流保持一致？
• RQ2 在体积约束下，该格式是否收敛于正确的 Wulff 集？收敛速率如何？
• RQ3 该格式在模拟非光滑或非 Wulff 初始形状（如环面或圆形）的界面演化时表现如何？
• RQ4 该格式是否如其谱形式所暗示的那样，具有无条件稳定性且比有限差分或有限元相场格式更精确？
• RQ5 在 $ L^1 $ 范数下，扩散界面逼近精确 Wulff 集的定量收敛速率是多少？

主要发现

• 该格式在粘性解意义下与各向异性平均曲率流一致，即使核 $ K_{\rho,t} $ 非正且其二阶矩不属于 $ L^1(\mathbb{R}^d) $，这构成了重大的分析挑战。
• 数值模拟表明，界面从 Wulff 集演化时，半径按 $ R(t) = \sqrt{R_0^2 - 2t} $ 减小，与理论预测一致。
• 在体积约束下从圆形初始形状演化时，最终形状收敛于 Wulff 集，计算结果与精确 Wulff 集之间的 $ L^1 $ 误差为 $ \epsilon $ 阶，表明一阶收敛。
• 该格式表现出无条件稳定行为，并且在二维和三维模拟中，其精度高于标准的有限差分或有限元相场方法。
• 在三维模拟中，该格式成功捕捉了在 $ \phi^o_4 $ 和 $ \phi^o_5 $ 作用下环面的演化，生成了具有预期对称性的 Wulff 集。
• 特征函数误差在 $ L^1 $ 范数下观测到 $ O(\epsilon) $ 的收敛速率，证实了该格式在逼近尖锐界面极限时的鲁棒性与精确性。