[論文レビュー] The Hierarchical Poincare-Steklov (HPS) solver for elliptic PDEs: A tutorial

Per‐Gunnar Martinsson|arXiv (Cornell University)|Jun 3, 2015

Numerical methods for differential equations参考文献 11被引用数 8

ひとこと要約

本稿では、滑らかな解を有する2次元領域における変係数楕円型PDEを解くための直接的で高次スペクトル法である階層的ポアンカレ＝ステクロフ（HPS）スolverを提案する。この手法は木構造の領域分割を用い、葉領域ではスペクトルコロケーションによりディリクレからノイマン（DtN）写像を計算し、連続性およびフラックス条件を用いて階層的に統合する。これにより、前処理の計算量が$O(N^{1.5})$、解法の計算量が$O(N/log N)$に抑えられ、反復法が収束しない高周波数の問題において特に有効である。

ABSTRACT

A numerical method for variable coefficient elliptic problems on two dimensional domains is described. The method is based on high-order spectral approximations and is designed for problems with smooth solutions. The resulting system of linear equations is solved using a direct solver with $O(N^{1.5})$ complexity for the pre-computation and $O(N \log N)$ complexity for the solve. The fact that the solver is direct is a principal feature of the scheme, and makes it particularly well suited to solving problems for which iterative solvers struggle; in particular for problems with highly oscillatory solutions. This note is intended as a tutorial description of the scheme, and draws heavily on previously published material.

研究の動機と目的

• 変係数楕円型PDEの直接スolverを、2次元領域における滑らかな解に対して開発すること。
• 高周波数解を有する問題において反復法の限界を克服すること。
• 木構造の領域分割を用いた階層的スペクトルDtN写像のアセンブリにより、最適な計算量を達成すること。
• 実装のための実用的フレームワークを、既存のスペクトル法およびマルチフロント法に基づいて提供すること。
• 共鳴によってDtN作用素が悪条件化する場合に備え、インピーダンス写像を用いて安定性を確保すること。

提案手法

• 領域は、葉領域がチェビシェフノード上での高次スペクトルコロケーションにより解像される木構造の長方形パッチにタイル化される。
• 各葉領域に対して、辺ごとに$q$個のガウスノードにおける境界補間を用いて、$4q \times 4q$の密行列として離散的ディリクレからノイマン（DtN）写像が計算される。
• 親パッチは、共有界面におけるポテンシャルとフラックスの連続性をブロック行列演算で強制することで、子パッチのDtN写像を統合して形成される。
• 解作用素は二段階のアルゴリズムにより構築される：上向きパスでプロキシ（DtN写像）を構築し、下向きパスで境界データおよび源項から解を再構成する。
• グリーン関数が周期的でない場合には、DtN行列の低ランク構造を活用し、H行列形式によりデータスパース表現を用いることで、最適な$O(N)$の計算量を達成できる。
• 共鳴問題では、DtN写像の代わりに安定な代替手段としてインピーダンス写像が用いられ、ユニタリで常に存在する作用素を保証する。

実験結果

リサーチクエスチョン

• RQ1 変係数楕円型PDEの直接スolverを、近似的に最適な計算量と高次精度で構築可能か？
• RQ2 階層的領域分割とスペクトル法を組み合わせることで、滑らかな解と振動的挙動を効率的に扱えるか？
• RQ3 階層的DtN写像を用いたシステムのアセンブリおよび解法の計算量はどの程度か？
• RQ4 共鳴によって標準的なDtN写像が無効化される状況において、この手法はどのようにして安定性を保てるか？
• RQ5 この手法は体積力および非同次源項を有する問題へ拡張可能か？

主な発見

• HPSスolverは、前処理の計算量が$O(N^{1.5})$、解法フェーズの計算量が$O(N\log N)$に抑えられ、古典的なネストド・ディセクションと一致する。
• 非周期的グリーン関数の場合、H行列形式によるDtN行列の低ランク構造を活用することで、最適な$O(N)$の計算量を達成できる。
• 本手法は、反復法が収束しない高周波数解を有する問題に対して特に効果的である。
• インピーダンス写像の使用により、共鳴状態でも安定性が保証され、標準的なDtN写像の代替として強力な代替手段を提供する。
• アルゴリズムは数値的に安定でスケーラブルであり、プロキシの階層的上向きアセンブリの後、下向きパスで解を効率的に再構成できる。
• 本手法は、純ディリクレ問題および体積力付き問題の両方に対して有効であることが実証されており、特解と同次解を別々に処理する二段階アルゴリズムが用いられている。