Uniform inf–sup stability of quartic and quintic Scott–Vogelius elements on Freudenthal meshes: a protected raw edge-star lifting

技术摘要：Freudenthal 网格上四次与五次 Scott–Vogelius 元的一致性 inf–sup 稳定性

问题陈述
本文研究了在均匀三维 Freudenthal 三角剖分上，用于 Stokes 方程的 Scott–Vogelius 有限元的均匀 inf–sup 稳定性。具体针对多项式次数 $k=4$ 和 $k=5$ 的情况。虽然 Zhang 此前已经证明了对于 $k \ge 6$ 的网格一致稳定性，但其构造依赖于插入六次气泡（degree-six bubbles）来修复在构建边迹（edge traces）后产生的单元级散度均值。本文识别出的关键差距在于：对于 $k=4$ 和 $k=5$，边迹函数本身已经具备了 $k \ge 4$ 的次数，然而现有的构造缺乏一种真正局部且与阶数无关的提升机制（lifting mechanism）来处理这些迹，而无需依赖于高阶情形中使用的顺序“借用”（borrowing）逻辑。缺失的部分是一个局部的边阶段提升（edge-stage lifting），该提升结合一个独立的均值修复策略，能够证明对于所有固定的 $k \ge 4$ 的稳定性。

方法论
作者采用了一种基于精确有理算术和整数矩阵恒等式的计算机辅助证明策略，来构造一个“受保护的原始边-星提升”（protected raw edge-star lifting）。该方法通过以下阶段进行：

1. 几何枚举与分类：
   作者对 Freudenthal 网格中所有可能的边配置进行了分类。通过分析相连的四面体和立方体，他们确定了恰好七种几何类和三十七种定向/边界配置。该枚举涵盖了所有边类型（内部边、边界边、面对角线边、体对角线边），适用于每边包含 $N \ge 3$ 个立方体的网格。

2. 精确局部提升构造：
   针对这 37 种配置以及次数 $k=4, 5$，作者构造了一个局部线性映射。

   • 源空间（Source Space）： 目标边上的散度迹空间，受连续性约束（奇异边及边界相容关系）限制。
   • 目标空间（Target Space）： 与该边相连的立方体并集（即“星形区域”，star）上的速度系数。
   • 保护（Protection）： 该映射旨在重现目标边迹，同时确保所有非目标边迹和顶点值保持为零。
   • 验证： 这些映射并非通过非正式的“借用”自由度次数而得出，而是由精确的整数矩阵进行认证。源迹空间被表征为整数矩阵 $H$ 的核，且提升映射是通过满足特定代数恒等式（例如 $HB=0, DB=s_0 I$）的整数矩阵 $B, D, C$ 来构造的。

3. 范数估计：
   利用精确的有理 Bernstein 质量与刚度计算，作者计算了提升映射的 $H^1$ 范数。他们建立了一个统一的参考补丁界限 $C_{ref} < 385$，适用于所有配置。由于仿射缩放性质，该界限与网格尺寸 $h$ 无关。

4. 组装与重叠控制：
   为了组装全局解，作者基于边方向和坐标模 3 的端点坐标定义了一种着色方案。这产生了一个 189 色着色方案，其中相同颜色的补丁是互不相交的。最大支持重叠被限制在 19 以内。组装的依赖深度为零，这意味着所有的边提升可以同时应用，而不存在顺序依赖。

5. 均值修复：
   原始边提升并不一定能保持单个四面体上散度的零均值。作者引入了一个单独的“受保护原始边-星”引理，并结合一个四次（$P_4$）宏补丁算子。该算子通过一个双立方体补丁修复单元级的散度均值，从而确保最终场在每个四面体上具有零均值，且不改变边迹。这种四次修复取代了 Zhang 原本构造中用于均值修正阶段的六次气泡，而新的边阶段提升则处理了 $k=4, 5$ 的迹匹配问题。

核心贡献

• 解决 $k=4, 5$ 的差距： 本文首次证明了 Scott–Vogelius 元在 Freudenthal 网格上对于多项式次数 4 和 5 的均匀 inf–sup 稳定性，完成了对所有 $k \ge 4$ 的结果证明。
• 受保护的边-星提升： 引入了一种“受保护”的提升机制，该机制是局部的、与阶数无关的，并由精确整数矩阵认证，将 Zhang 关于边阶段次数要求的观察转化为对于低阶情形的严格证明，而无需依赖顺序的“借用一个自由度”逻辑。
• 完整的几何普查： 对 37 种边配置和 117 种边界装饰的源包络进行了严格枚举，并验证了其对于 Freudenthal 模式是详尽无遗的。
• 计算机辅助验证： 构造过程依赖于一个包含精确有理证书的可复现软件包。作者利用独立的 Python 验证器和 Lean 形式化工具（使用 Mathlib）来验证代数恒等式、范数界限和组装属性。Lean 代码在进行核心代数检查时，不包含任何 sorry 或未证明的公理。

结果

• 定理 1.1： 对于每一个固定的整数 $k \ge 4$，存在一个与网格尺寸 $h$ 无关的常数 $\beta_k > 0$，使得 Scott–Vogelius 对 $(V^k_h, Q^k_h)$ 在 Freudenthal 网格上满足 inf–sup 条件。
• 统一界限： 构造的右逆 $R_h: Q^k_h \to V^k_h$ 满足 $|R_h q_h|_1 \le C_k \|q_h\|_0$，其中 $C_k$ 与 $h$ 无关。
• 特定常数： 边提升的参考补丁界限为 $C_{ref} < 385$。全局组装涉及 189 色着色，最大重叠为 19。

意义与主张
本文声称提供了一种“真正局部且与阶数无关的提升”，它将关于边阶段次数要求的观察转化为对于低阶情形的严格证明。通过使用有限精确矩阵代替顺序的“借用一个自由度”逻辑进行边阶段处理，并通过使用四次宏补丁算子来修复单元均值（取代了 Zhang 工作中用于该特定目的的六次气泡），作者消除了在 $k=4, 5$ 时进行临时次数提升的需求。该工作证明，仅使用四次和五次元素即可实现稳定性，而无需依赖于之前 $k \ge 6$ 分析中用于均值修复阶段的六次气泡修正。使用精确有理证书和 Lean 形式化确保了稳定性常数和结构属性在没有浮点误差的情况下经过数学验证。该结果统一了 Freudenthal 网格上所有 $k \ge 4$ 的稳定性理论。