A stabilized dual-SAV parametric finite element framework for constrained planar geometric flows with mesh regularization

技术摘要：一种用于约束平面几何流的稳定化双标量辅助变量参数有限元框架

问题陈述
约束几何流的参数有限元离散化面临三个同时存在的挑战：处理高阶几何刚度（其在显式方案中限制时间步长）、管理网格退化（顶点的聚集或扩散）以及施加非线性全局约束（如面积或长度守恒）。现有方法往往难以在单个线性大规模框架内结合动态网格正则化、能量稳定的半隐式时间步进以及非线性约束施加。完全隐式方法虽然稳健，但需要求解大型非线性系统，而显式处理则遭受严重的稳定性限制。此外，将用于网格重分布的人工切向速度纳入保结构方案而不引入虚假的法向驱动力，仍然十分微妙。

方法论
本文提出了一种用于闭合平面曲线的**稳定化双标量辅助变量（SAV）参数有限元（PFEM）**框架。其核心创新在于将物理几何能量与人工网格正则化能量分离为两个独立的标量辅助变量（$r_g$ 和 $r_m$）。

1. 双能结构：

   • 物理几何能量（$E_{geom}$）： 驱动曲线的法向演化。引入辅助变量 $r_g$ 以线性化非线性几何力。
   • 网格正则化能量（$E_{mesh}$）： 一种人工能量（通常为加权狄利克雷能量），仅用于切向重分布网格点。引入独立的辅助变量 $r_m$。关键在于，来自 $E_{mesh}$ 的驱动力仅投影到切向方向。这确保了网格重分布能改善参数化，而不改变预设的法向几何演化。

2. 半隐式冻结度量离散化：
   该方案采用半隐式时间离散化，其中所有依赖于几何的项（度量、力、约束）均在时间步 $n$ 的已知曲线上进行评估（“冻结度量”）。

   • 零阶稳定化： 稳定化双线性形式（拉普拉斯或双调和）作用于新速度本身，而非速度差。这一选择对于维持能量耗散估计所需的代数抵消至关重要，并避免了速度历史的传播。
   • 线性空间响应： 一旦标量辅助变量和拉格朗日乘子确定，法向和切向速度的空间问题即变为线性。

3. 全局约束的代数块约简：
   该方法通过代数块约简处理 $K$ 个全局几何约束（如面积、长度）。

   • 首先求解大规模线性空间响应问题以生成“响应向量”。
   • 剩余的非线性系统仅涉及几何辅助变量（$r_g$）和 $K$ 个拉格朗日乘子。
   • 因此，约简后非线性系统的维度为 $K+1$，独立于网格顶点数 $N$。对于同时的面积和长度约束（$K=2$），这将产生一个三维非线性系统。

主要贡献
本文概述了五项具体贡献：

1. 双 SAV 公式： 引入针对物理能量和网格能量的独立辅助变量，并将网格力严格耦合至切向运动。
2. 冻结度量半隐式方案： 构建一种时间离散化方案，在满足离散修正能量耗散估计的同时，产生线性空间响应问题。
3. 全离散 PFEM： 一种在离散层面保持代数抵消结构的公式，确保修正能量的耗散。
4. 与维度无关的约束处理： 推导了一种代数块约简方法，将大规模线性空间求解与非线性约束施加解耦，将非线性求解的维度降低至 $K+1$。
5. 数值验证： 对二阶和四阶流（曲线缩短、保面积曲线缩短、曲线扩散和 Helfrich 型流）进行了全面测试。

结果
数值实验表明：

• 精度： 曲线缩短流（收缩圆）表现出约一阶的时间精度，这与半隐式后向欧拉结构一致。
• 网格质量： 在保面积曲线缩短中，所提出的切向重分布防止了拉格朗日参考案例中观察到的严重顶点聚集，在严格施加面积约束的同时保持了近乎均匀的网格比率（$Q_{mesh} \approx 1$）。
• 非局部度量： 该框架成功地将非局部 $H^{-1}$ 度量纳入曲线扩散流，而无需通过拉格朗日乘子显式施加面积约束（面积由 $H^{-1}$ 速度的零均值性质自然保持）。
• 高阶约束： 对于具有同时面积和长度约束（$K=2$）的 Helfrich 型流，约简后的非线性系统（三维）被高效求解。约束被强制至机器精度（误差 $\sim 10^{-15}$），且修正的几何 SAV 能量单调耗散。
• 可扩展性： 计算成本分析显示，随着顶点数 $N$ 的增加，约简非线性系统的牛顿迭代次数保持恒定（约 2.2 次），而成本随线性响应求解的次数呈线性增长。

意义与主张
本文主张，这项工作的主要意义不在于各个组件（切向速度、SAV 或约束处理），而在于它们在单个参数有限元框架内的代数解耦。通过将物理几何流与人工网格机制分离，并将非线性约束施加约简为低维系统，该方法实现了：

• 稳定性： 通过零阶稳定化和半隐式步进，稳健处理高阶刚度。
• 效率： 避免大型非线性鞍点求解；非线性成本与网格分辨率无关。
• 保结构性： 维持离散修正能量耗散并准确施加全局几何约束。

作者指出，当前的研究范围仅限于闭合平面曲线。将该框架扩展至 $\mathbb{R}^3$ 中的演化曲面被确定为未来工作，这需要额外的曲面网格正则化机制和稳定性分析。本文还承认，修正的 SAV 能量与原始几何能量之间的关系值得进一步研究，特别是对于较大的时间步长。