[论文解读] A Bernstein-Bezier Sufficient Condition for Invertibility of Polynomial Mapping Functions

Stephen A. Vavasis|ArXiv.org|Aug 11, 2003

Advanced Numerical Analysis Techniques参考文献 13被引用 6

一句话总结

本文提出了一种基于雅可比矩阵列的伯恩斯坦-贝齐尔形式的多项式映射函数的全局与局部可逆性的充分条件。通过确保每个导数分量的控制点凸包可由超平面分离，该方法保证了可逆性，并具备理想的仿射不变性特性，从而能够验证仿射映射的小扰动。

ABSTRACT

We propose a sufficient condition for invertibility of a polynomial mapping function defined on a cube or simplex. This condition is applicable to finite element analysis using curved meshes. The sufficient condition is based on an analysis of the Bernstein-Bézier form of the columns of the derivative.

研究动机与目标

• 为二维和三维中立方体与单形等参考单元上的多项式映射提供全局与局部可逆性的通用充分条件。
• 解决等参数有限元方法中可逆性验证的挑战，其中标准测试（如高斯点评估）不足以满足需求。
• 开发一种既可计算测试又具备仿射不变性特性的条件，以确保对仿射映射及其近似映射的可靠性。
• 将以往仅限于低次多项式或特定单元类型的已有充分条件扩展为适用于任意多项式次数与区域的统一框架。
• 为界定条件数和高阶导数奠定基础，这对于高阶有限元方法中的收敛性至关重要。

提出的方法

• 将多项式映射的雅可比矩阵表示为伯恩斯坦-贝齐尔（BB）形式，将每一列表达为控制点的凸组合。
• 定义集合 $ G_{\xi}, G_{\eta}, G_{\zeta} $ 为偏导数 $ F_\xi, F_\eta, F_\zeta $ 的控制点。
• 基于这些集合由超平面分离建立充分条件：对于每个导数方向，控制点的凸包必须能被具有统一法向量的超平面分离。
• 通过半空间的几何分离证明映射是局部可逆的，确保雅可比矩阵非奇异。
• 通过证明映射在所有方向上的方向导数严格为正，将局部可逆性证明扩展至全局可逆性，从而确保映射不会发生缠绕。
• 引入通过内积 $ \mathbf{h}^T \mathbf{f}_i $ 表征的条件，实现对分离性条件的高效数值检验。

实验结果

研究问题

• RQ1 能否为二维和三维中所有多项式次数和参考域（立方体与单形）的多项式映射可逆性，制定出单一的通用充分条件？
• RQ2 如何利用伯恩斯坦-贝齐尔形式推导出雅可比矩阵非奇异性的几何且可计算测试的条件？
• RQ3 所提出的条件是否具备仿射不变性，以确保其对仿射映射及其小扰动均成立？
• RQ4 该条件相较于现有方法的限制程度如何？是否可进一步强化以界定条件数或高阶导数？
• RQ5 在实际等参数网格生成中，当内部自由度由 Lenoir 公式确定时，该条件的行为表现如何？

主要发现

• 所提出的条件可保证在 $ I^d $ 和 $ \Delta^d $ 上（$ d=2,3 $）任意多项式次数 $ p $ 下，多项式映射的局部与全局可逆性。
• 该条件基于雅可比矩阵列的伯恩斯坦-贝齐尔控制点凸包由超平面分离，从而确保雅可比矩阵非奇异。
• 该方法具有仿射不变性：对非退化仿射映射及其小扰动始终满足，证实了其仿射不变性特性。
• 该条件弱于某些现有专用测试方法（如 Ushakova 针对 $ p=1 $, $ I^3 $ 的方法），但更具通用性，适用于更高次多项式与四面体单元。
• 初步计算证据表明，使用 Lenoir 公式确定内部自由度可减小控制点集合 $ G_{\xi}, G_{\eta}, G_{\zeta} $ 的直径，从而提高条件满足的可能性。
• 该框架可进一步拓展以界定条件数 $ \text{cond}(F) $ 和高阶导数，这对满足 Ciarlet-Raviart 收敛条件至关重要。