[論文レビュー] 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.

研究の動機と目的

• 2次元および3次元の立方体や単体などの参考要素における多項式写像の全域的・局所的逆写像可能性に対する一般な十分条件を提供すること。
• 標準的なテスト（例：ガウス点評価）では不十分な、等パラメトリック有限要素法における逆写像可能性の検証の課題に対処すること。
• 計算的にテスト可能であり、アフィン写像および近似的にアフィンな写像に対して信頼性を保証するアフィン性の性質を有する条件を開発すること。
• 従来、低次の多項式や特定の要素タイプに限定されていた十分条件を、任意の多項式次数および領域に適用可能な統一的枠組みへと拡張すること。
• 高次有限要素法における収束性に不可欠な条件数や高階導関数の境界を付ける基盤を築くこと。

提案手法

• 多項式写像のヤコビ行列をベルンシュタイン＝ベジエ（BB）形式で表現し、各列を制御点の凸結合として表現する。
• 各偏導関数 $ F_\xi, F_\eta, F_\zeta $ の制御点をそれぞれ集合 $ G_{\xi}, G_{\eta}, G_{\zeta} $ として定義する。
• 各導関数方向について、制御点の凸包が一様な法線ベクトルを持つ超平面によって分離可能であるという十分条件を確立する。
• 半空間による幾何的分離を用いて、写像が局所的に逆写像可能であることを証明し、ヤコビ行列が非特異のまま保たれることを保証する。
• 方向微分係数がすべての方向で厳密に正であることを示すことにより、局所的逆写像可能性の証明を全域的逆写像可能性へと拡張する。
• 内積 $ \mathbf{h}^T \mathbf{f}_i $ を用いた条件の特徴付けを導入し、分離条件の効率的な数値的テストを可能にする。

実験結果

リサーチクエスチョン

• RQ1 任意の多項式次数および参考領域（2次元および3次元の立方体および単体）にわたる多項式写像の逆写像可能性に対して、一様な十分条件を定式化できるか？
• RQ2 ベルンシュタイン＝ベジエ形式を活用することで、ヤコビ行列の非特異性を保証する幾何的かつ計算的にテスト可能な条件を導出できるか？
• RQ3 提案された条件はアフィン性の性質を有するか？ これにより、アフィン写像およびその小さな摂動に対しても成立するか？
• RQ4 既存の手法と比較して、この条件はどの程度制限が少なく、条件数や高階導関数の境界を付けるために強化可能か？
• RQ5 実用的な等パラメトリックメッシュ生成において、内部自由度がレノワールの公式によって決定される場合、この条件はどのように振る舞うか？

主な発見

• 提案された条件は、$ d=2,3 $ の $ I^d $ および $ \Delta^d $ 上で、多項式次数 $ p $ にかかわらず多項式写像の全域的・局所的逆写像可能性を保証する。
• この条件は、ヤコビ行列の各列のベルンシュタイン＝ベジエ制御点の凸包が超平面によって分離可能であることに基づく。これによりヤコビ行列が非特異のまま保たれる。
• この手法はアフィン不変性を有する：非退化なアフィン写像およびその小さな摂動に対して常に満たされる。これによりアフィン性の性質が確認される。
• 一部の既存の専用テスト（例：$ p=1 $, $ I^3 $ におけるウシャコワの手法）より条件が弱いが、より一般化されており、高次多項式や四面体要素にも適用可能である。
• 初期の計算的証拠から、内部自由度にレノワールの公式を適用することで、$ G_{\xi}, G_{\eta}, G_{\zeta} $ の制御点集合の直径が小さくなる傾向があることが示唆されており、条件を満たしやすくなる。
• この枠組みは、収束性に不可欠な条件数 $ \text{cond}(F) $ や高階導関数の境界を付ける可能性を有し、シアレ・ラヴィアルトの収束条件を満たすために重要である。