[論文レビュー] Optimal quadrature rules for isogeometric analysis

Michael Bartoň, Victor M. Calo|arXiv (Cornell University)|Nov 12, 2015

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

ひとこと要約

本稿では、等間隔でないスプライン空間における有限領域上での正確な積分を保証する最小の求積点数を達成するため、ホモトピー連続法に基づく最適ガウス求積則の導出手法を提案する。この手法は、多項式に対する古典的ガウス求積則から出発し、目的のスプライン空間の連続性と knot 配置に合わせて knot ベクトルを変形することで、正の重みと機械精度に近い正確性を持つ求積則を生成する。特に、低連続性スプラインでは、漸近的パターンへの収束が速く、明確な構造を示す。

ABSTRACT

We introduce optimal quadrature rules for spline spaces that are frequently used in Galerkin discretizations to build mass and stiffness matrices. Using the homotopy continuation concept [6] that transforms optimal quadrature rules from source spaces to target spaces, we derive optimal rules for splines generated above finite domains. Starting with the classical Gaussian quadrature for polynomials, which is an optimal rule for a discontinuous odd-degree space, we derive rules for target spaces of higher continuity. We further show how the homotopy methodology handles cases where the source and target rules require different numbers of optimal quadrature points. We demonstrate it by deriving optimal rules for various odd-degree spline spaces, particularly with non-uniform knot sequences and non-uniform multiplicities. We also discuss convergence of our rules to their asymptotic counterparts, that is, the analogues of the midpoint rule of Hughes et al. [34], that are exact and optimal for infinite domains. For spaces of low continuities, we numerically show that the derived rules quickly converge to their asymptotic counterparts as the weights and nodes of a few boundary elements differ from the asymptotic values.

研究の動機と目的

• 有限領域上での正確な積分を保証しつつ、求積点数を最小化するスプライン空間における最適求積則の開発。
• 非一様かつ高連続性を持つスプライン空間に対して、閉形式解または数値的に安定した方法が不足している問題に対処。
• 正の重みと高い数値的正確性を備えた求積則を提供することで、効率的かつ安定したガレルキン離散化を可能に。
• 導出された求積則が、特に無限領域の類似物に対して漸近的対応にどの程度収束するかを調査。
• 既存の $ C^1 $ スプライン用再帰的手法を、一様なホモトピー枠組みにより任意の次数および連続性に一般化。

提案手法

• ホモトピー連続法を用いて、既知のソース空間（例：不連続多項式）から、所望の連続性および knot 配置を持つターゲットスプライン空間への最適求積則の解の経路を追跡する。
• モーメント条件から導かれる多項方程式系のゼロとして求積則を表現し、スプライン空間全体での正確な積分を保証する。
• 既知のソース構成（例：$ C^1 $ 用の重複 knot）からターゲット構成（例：単一 knot）への knot ベクトルの連続的変形を実行し、方程式系の代数的構造を維持する。
• ソース空間とターゲット空間で求積点数が異なる場合に備え、ホモトピー追跡中にシステムを動的に調整する。
• 数値的連続解法ソルバを用いて、変形された方程式系の根（求積ノードおよび重み）を計算し、高い正確性と正の重みを確保する。
• スプライン空間に属する乱数関数を用いたテストにより妥当性を検証し、正確な値と16桁まで一致することを確認。

実験結果

リサーチクエスチョン

• RQ1 任意の次数および連続性を持つスプライン空間に対して、有限領域上での最適求積則を体系的に導出する方法は何か？
• RQ2 ホモトピー連続法は、ソース空間（例：$ C^1 $ 三次スプライン）から、異なる knot 重複度および連続性を持つターゲット空間への最適求積則を効果的に移行できるか？
• RQ3 要素数が増加するに従い、導出された求積則がその漸近的対応（例：無限領域における中点則）にどの程度の速度で収束するか？
• RQ4 導出された求積則は、部分的最適または負の重みを持つ手法と比較して、どの程度正の重みと数値的安定性を維持できるか？
• RQ5 ソース空間とターゲット空間で必要な求積点数が異なる場合に、本手法がその処理を一般化して扱えるか？

主な発見

• ホモトピー連続法は、非一様 knot ベクトルおよび高連続性スプラインを含む、さまざまなスプライン空間に対して、最小の求積点数を持つ最適ガウス求積則を成功裏に導出できた。
• $ C^0 $ および $ C^1 $ スプラインでは、導出された求積則が漸近的対応に急速に収束し、最初の数個の境界要素以外は機械精度を上回る差異を示さない。
• すべての導出された求積則は正の重みを有しており、負の重みを持つ手法と比較して数値的安定性が向上している。
• 本手法は機械精度の正確性を達成した：乱数関数のテスト積分結果が正確な値と16桁まで一致した。
• $ C^1 $ 五次スプラインおよび $ C^1 $ 三次スプラインの求積則には明確な代数的構造が存在し、反復的ソルバを用いずに閉形式で計算可能である。
• 表10は、$ C^0 $ および $ C^1 $ スプラインにおいて、漸近的規則から逸脱する境界要素が1～4個に限られることを示しており、効率的な実装に適した準周期的構造であることを示している。