[論文レビュー] Isogeometric Boundary Element Method with Hierarchical Matrices

Jürgen Zechner, Benjamin Marussig+2 more|arXiv (Cornell University)|Jun 11, 2014

Advanced Numerical Analysis Techniques参考文献 14被引用数 5

ひとこと要約

本稿では、NURBS幾何学とカーネル補間を用いて、完全に密なシステム行列を圧縮する、等稜的境界要素法（IGA-BEM）と階層行列（H-行列）を組み合わせた手法を提案する。この手法は、近似的に線形の計算複雑性と高い精度を達成し、2次元数値実験により最適収束性と効率的な行列圧縮が確認された。

ABSTRACT

In this work we address the complexity problem of the isogeometric Boundary Element Method by proposing a collocation scheme for practical problems in linear elasticity and the application of hierarchical matrices. For mixed boundary value problems, a block system of matrices similar to Galerkin formulations is constructed allowing an effective application of that matrix format. We introduce a strategy for the geometric bisection of surfaces based on NURBS patches. The approximation of system matrices is carried out by means of kernel interpolation. Numerical results are shown that prove the success of the formulation.

研究の動機と目的

• 完全に密なシステム行列に起因する等稜的境界要素法の高い計算コストを低減すること。
• 線形弾性問題における等稜的配置境界要素法フレームワークに階層行列（H-行列）を統合すること。
• NURBSを用いて幾何的正確性を保持しつつ、大規模問題の効率的行列圧縮と解法を可能にすること。
• NURBSのノット挿入と凸包の性質に基づく幾何的二分法戦略を用いて、行列ブロックのクラスタリングを実現すること。
• ラグランジュ多項式とQR圧縮を用いたカーネル補間により、行列近似の高精度を確保すること。

提案手法

• 境界幾何学を正確に表現するためにNURBSパッチを用いることで、幾何的離散化誤差を排除する。
• 線形弾性問題に適したラメ＝ナヴィエ作用素に基づく境界積分方程式を用いた配置BEMを適用する。
• 混合境界値問題のガレルキン形式に類似したブロック行列系を構築する。
• NURBSのノット挿入と凸包から導かれるバウンディングボックスを用いて、表面の幾何的二分法を実装する。
• 許容可能なクラスタに対して、基本解（ケルビンの解）のカーネル補間を用いて遠方行列ブロックを近似する。
• 低ランク因子分解（A·S·Bᵀ）とQR分解を用いて、低ランク行列ブロックを圧縮し、ストレージと計算コストを削減する。

実験結果

リサーチクエスチョン

• RQ1 階層行列は、線形弾性問題における等稜的BEMの計算複雑性を効果的に低減できるか？
• RQ2 NURBSに基づく幾何学は、H-行列の圧縮に適した効率的な幾何的クラスタリングを可能にするか？
• RQ3 本手法のH-行列近似は、弾性における単層および二重層ポテンシャルにおいて、精度と収束特性はどのように振る舞うか？
• RQ4 カーネル補間とQR分解による行列圧縮は、精度を損なわずにストレージと解法の効率を向上させるか？
• RQ5 実際の計算において、近似的に線形の複雑性を達成しつつ、最適収束率を維持できるか？

主な発見

• 2次元数値実験により、提案手法のIGA-BEMとH-行列は、単層および二重層ポテンシャルの両方で最適収束率を達成した。
• 行列圧縮比（cₕ）は問題サイズにほぼ線形に近づく傾向を示し、実際の計算で近似的に線形の複雑性を達成していることを示している。
• 高い圧縮率下でも高い精度を維持しており、誤差は補間次数とランク選択によって制御可能である。
• ラグランジュ多項式とサポートポイントを用いたカーネル補間により、遠方行列ブロックの有効な低ランク近似が可能になった。
• H-行列ブロックのQR圧縮により、ストレージ要件がさらに削減されたが、精度の著しい損失は生じなかった。
• 形状関数に接続する形で法線方向微分をシフトすることで、実装が双曲線的微分の評価を正確に処理できた。