[論文レビュー] Sweeping Preconditioner for the Helmholtz Equation: Moving Perfectly Matched Layers

Björn Engquist, Lexing Ying|arXiv (Cornell University)|Jul 25, 2010

Electromagnetic Scattering and Analysis参考文献 24被引用数 9

ひとこと要約

本稿では、2次元および3次元の変数係数ヘルムホルツ方程式に対して、移動する完全一致層（PML）を用いて、$LDL^t$分解におけるスシュール補行列を効率的に近似する、新しいスイーピングプリコンディショナを提案する。この手法は、GMRESと組み合わせることで、周波数や問題サイズに依存しない反復回数で線形の応用コストと収束性を達成し、高周波数ヘルムホルツ問題のスケーラブルな反復的解法を可能にする。

ABSTRACT

This paper introduces a new sweeping preconditioner for the iterative solution of the variable coefficient Helmholtz equation in two and three dimensions. The algorithms follow the general structure of constructing an approximate $LDL^t$ factorization by eliminating the unknowns layer by layer starting from an absorbing layer or boundary condition. The central idea of this paper is to approximate the Schur complement matrices of the factorization using moving perfectly matched layers (PMLs) introduced in the interior of the domain. Applying each Schur complement matrix is equivalent to solving a quasi-1D problem with a banded LU factorization in the 2D case and to solving a quasi-2D problem with a multifrontal method in the 3D case. The resulting preconditioner has linear application cost and the preconditioned iterative solver converges in a number of iterations that is essentially indefinite of the number of unknowns or the frequency. Numerical results are presented in both two and three dimensions to demonstrate the efficiency of this new preconditioner.

研究の動機と目的

• 大きな、高次元で不確定な線形方程式系を生じる高周波数ヘルムホルツ方程式の解法における計算的課題に対処すること。
• 従来のプリコンディショナが3次元で高いセットアップコストと理論的根拠の欠如を抱えていたという限界を克服すること。
• 問題サイズや周波数に依存しない線形の応用コストと強固な収束性を持つプリコンディショナを開発すること。
• PMLの物理的洞察と高速直接解法を活用して、2次元および3次元のヘルムホルツ問題の効率的反復的解法を可能にすること。
• 非一様、適応的、非一様グリッドに適用可能な柔軟なフレームワークを提供し、並列処理を支援すること。

提案手法

• ドメイン内部に導入された移動PMLを用いて、$LDL^t$分解におけるスシュール補行列を近似する。
• 2次元では、最適な再順序付けを経て帯行列LU分解を適用し、PMLを準1次元問題として扱う。
• 3次元では、多スケルタル法を用いて各スシュール補行列を処理し、PMLを準2次元問題として扱う。
• PML層の局所性と狭い幅を活かして、低コストで高速な直接逆行列計算を実現する。
• 吸収境界から始めて、層ごとに順次プリコンディショナを構築し、PMLを用いて層間の相互作用をモデル化する。
• 得られたプリコンディショナをGMRES反復ソルバーと組み合わせることで、強固な収束性を達成する。

実験結果

リサーチクエスチョン

• RQ1 移動するPMLを用いて、ヘルムホルツ方程式のスイーピングプリコンディショナにおけるスシュール補行列を効率的に近似できるか？
• RQ2 提案手法は2次元および3次元の両方で、線形の応用コストと周波数に依存しない収束性を達成するか？
• RQ3 セットアップコストと反復回数の観点から、従来の階層行列ベースの手法と比較して、本手法の性能はどのように異なるか？
• RQ4 非一様および適応的グリッドへも拡張可能か、効率性を維持できるか？
• RQ5 特に多スケルタル解法フェーズにおいて、並列処理に適しているか？

主な発見

• 2次元ではプリコンディショナの応用コストが$O(N)$、3次元では$O(N\log N)$であり、スケーラブルな反復的解法が可能である。
• 2次元および3次元の両方で、未知数の数$N$や周波数$\omega$にほとんど依存しないGMRES反復回数が得られた。
• 3次元問題において$159^3$の未知数と$\omega/(2\pi) = 5$の場合、反復回数は13～15回であり、解法時間は10分未満であった。
• より効率的な多スケルタル実装のおかげで、セットアップ時間は従来の数万秒から数百分間（例：$159^3$で約800秒）に短縮された。
• 異なる速度場やPMLの厚さに対してもロバストであり、$q$および$p$パラメータの変化に関わらず一貫した反復回数が得られた。
• 非一様および適応的グリッドへも拡張可能であり、各PML層内の多スケルタルソルバーは効率的な並列処理を支援する。