Upper bound of high-order derivatives for Wachspress coordinates on polytopes

技术摘要：多胞体上 Wachspress 坐标的高阶导数上界

问题陈述
广义重心坐标（GBCs）是在多胞体网格上构建有限元方法（FEM）的基础。虽然 GBCs 的梯度界（一阶导数）已得到充分确立，并足以用于二阶椭圆方程的$H^1$范数误差估计，但高阶有限元应用（例如板弯曲等四阶椭圆方程）需要高阶导数（$H^k$范数，$k > 1$）的上界。在此项工作之前，文献中缺乏此类界限，导致基于 GBC 的非协调有限元（如文献 [34] 中提出的 Morley 型单元）的误差分析存在理论空白。

方法论
作者聚焦于定义在简单、非退化、凸$d$维多胞体上的Wachspress 坐标。任意阶导数（$D^\alpha \phi_v$）上界的推导涉及以下技术步骤：

1. 代数分解：Wachspress 坐标$\phi_v$表示为比率$w_v / W$，其中$w_v$是多项式权重函数，$W$是所有权重之和。作者利用广义莱布尼茨公式和多元 Faà di Bruno 公式展开该商式的导数。这将问题简化为估计分子（$w_v$）和分母（$W$）的导数。
2. 几何参数：界限是根据两个关键几何量推导得出的：
   • $h_K$：多胞体$K$的直径。
   • $h^*$：$K$中非关联顶点与面之间的最小距离。
3. 导数估计：
   • 分子（$w_v$）：由于$w_v$是线性距离函数的乘积，其导数通过涉及面数量和直径$h_K$的组合论证进行界定。
   • 分母（$W$）：关键步骤在于建立$W$的严格正下界。作者通过分析面相对于点$x$的分布，利用 GBCs 的线性精度性质和简单多胞体的几何特性，证明了$W(x) \geq C \frac{(h^*)^{|F|}}{h_K^d}$。
4. 综合：结合$w_v$和$W$的导数上界与$W$的下界，作者推导出了$|D^\alpha \phi_v|$的显式上界。

主要贡献

1. 高阶界限的推导：本文提供了简单凸多胞体上 Wachspress 坐标任意阶导数的首个显式上界。主要结果（定理 3.7）指出，当$h^* = O(h_K)$（该条件意味着多胞体在特定意义上不是“细长”或退化的）时，导数满足：
   $$|D^\alpha \phi_v| \leq C h_K^{-|\alpha|}$$
   其中$C$取决于维度$d$、顶点数$|V|$、面数$|F|$以及多重指标$\alpha$。
2. 几何正则性分析：本文系统地比较了凸多胞体的各种形状正则性条件（例如最小边长、最大顶点数、块度参数）。它建立了这些条件之间的逻辑蕴含关系，表明假设$h^* = O(h_K)$（假设 H1）在实践上是合理的，并且蕴含了有限元中使用的其他标准正则性条件。
3. 尖锐性的澄清：作者阐明，虽然当$h^* = O(h_K)$时推导出的界限在渐近阶上是尖锐的，但当$h^* \ll h_K$时，它们可能并不尖锐。数值实验表明，在后一种情况下，导数的行为取决于“短”边相对于导数方向的取向。

数值结果
作者通过二维数值实验验证了其理论发现：

• 网格质量：他们证明，在消除“短”边后，质心 Voronoi 剖分（CVT）自然地满足条件$h^* = O(h_K)$。
• 界限验证：对于满足$h^* = O(h_K)$的随机多边形，计算出的$\Lambda_\alpha = \max \sum |D^\alpha \phi_i|$值与理论预测$O(h_K^{-|\alpha|})$一致。
• 退化情况：在$h^* \ll h_K$的多边形上的实验表明，理论上是保守的。导数的实际增长取决于短边的方向，某些方向导数表现出比一般界限所暗示的更慢的增长率。
• 应用：本文包含了将基于 Morley 型 GBC 的有限元应用于夹支板弯曲问题（四阶方程）的结果。数值收敛率（$H^2$半范数中的$O(h)$和$H^1$半范数中的$O(h^2)$）证实，推导出的导数界限足以证明此类问题的最优收敛性。

意义与主张
本文声称填补了多胞体有限元方法分析中的一个关键理论空白。通过提供必要的高阶导数界限，它使得针对四阶椭圆方程应用的基于 Wachspress 的有限元的最优收敛性得以严格证明。作者强调，他们的方法避免了以往工作（如 [14]）中高度技术化且特定于梯度的技术，而是采用了一种适用于任意阶的不同代数策略。此外，这项工作通过阐明各种形状正则性假设之间的关系，促进了对多胞体几何的理解，并断言条件$h^* = O(h_K)$是生成应用中高质量网格的稳健且实用的假设。