[論文レビュー] A polynomial multigrid smoother for the iterative solution of the heterogeneous Helmholtz problem

Wim Vanroose, Bram Reps+1 more|arXiv (Cornell University)|Dec 24, 2010

Electromagnetic Scattering and Analysis参考文献 18被引用数 3

ひとこと要約

本稿では、複素数に伸びたグリッドを用いて、元の作用素の固有値スペクトルに近いスペクトルを持つプリコンディショナを構築することで、不均一で不定なヘルムホルツ問題を解くためのマルチグリッドプリコンディショニング付きクリロフ部分空間法を提案する。多項式スムージング（具体的には3回のGMRES反復）を用いることで、さまざまなテスト問題において安定かつロバストな収束が達成され、グリッドの伸びるパラメータに関する導出された条件下で安定性が保証される。

ABSTRACT

Abstract. In this paper we develop a robust multigrid preconditioned Krylov subspace method for the solution of heterogeneous indefinite Helmholtz problems. The preconditioning operator is constructed by discretizing the original Helmholtz equation on a complex stretched grid. As this preconditioning operator has the same wavenumber as the original problem, and the only difference stems from the complex stretching of the discretization grid, its spectrum closely approximates the spectrum of the original Helmholtz operator, resulting in a fast converging Krylov subspace method. We have analyzed a multigrid cycle based on polynomial smoothers and have found a condition on the parameters used in the complex stretching of the discretization grid such that the existence of a stable third order polynomial smoother is guaranteed. In practice we use three iterations of GMRES as a smoother, and have found this to be a viable smoother for the multigrid preconditioning process. We apply the method to various test problems and report on the observed convergence rate. 1. Introduction. The

研究の動機と目的

• 不均一で不定なヘルムホルツ問題に対するロバストな反復解法の開発。これは、不定的かつ高周波数の振動を示すため、困難である。
• このような問題のプリコンディショニングの難しさに対処するため、元のヘルムホルツ作用素のスペクトルに近いスペクトルを持つプリコンディショナを構築する。
• 3次多項式スムージングのための複素数に伸びたパrameterに関する条件を導出し、マルチグリッドスムージングの安定性と効率を保証する。
• 多様なテスト問題に対して手法を検証し、一貫した収束特性を示す。

提案手法

• プリコンディショナは、元のヘルムホルツ方程式を複素数に伸びた計算グリッド上で離散化することによって構築され、元の問題と同じ波数を維持する。
• 複素数に伸びたグリッドは、元の作用素のスペクトル近似を改善し、プリコンディショナの有効性を高める。
• 実際には、3回のGMRES反復をスムージングとして用いる多項式スムージングを用いたマルチグリッドサイクルが採用される。
• 複素数に伸びたパラメータのパラメータに基づき、3次多項式スムージングの安定性を保証する理論的条件が導出される。
• マルチグリッド反復とクリロフ部分空間加速を組み合わせ、元の作用素とプリコンディショニング作用素のスペクトル類似性を活用する。
• 多様な不均一なヘルムホルツテスト問題にこの手法を適用し、収束速度とロバスト性を評価する。

実験結果

リサーチクエスチョン

• RQ1 元の不定的ヘルムホルツ作用素のスペクトルに近いスペクトルを持つマルチグリッドプリコンディショナを構築することは可能か？
• RQ2 3次多項式スムージングのマルチグリッドサイクルにおける安定性を保証する複素数に伸びたパラメータの条件は何か？
• RQ3 実際のこのクラスの問題に対して、3回のGMRES反復をスムージングとして用いることはどれほど有効か？
• RQ4 提案手法は、多様な不均一なヘルムホルツテスト問題においてロバストな収束を達成するか？
• RQ5 元の作用素とプリコンディショニング作用素のスペクトル類似性は、クリロフ部分空間の収束を加速するために活用可能か？

主な発見

• 複素数に伸びたグリッドの離散化によって構築されたプリコンディショナは、元のヘルムホルツ作用素のスペクトルに類似しており、クリロフ法の高速収束を可能にする。
• 複素数に伸びたパラメータに基づき、3次多項式スムージングの安定性を保証する条件が理論的に導出され、マルチグリッドサイクルの収束が保証される。
• 実際には、3回のGMRES反復をスムージングとして用いることで、有効かつ安定した性能が得られ、マルチグリッドプリコンディショニングフレームワークの有効な構成要素となる。
• 本手法は、多様な不均一なヘルムホルツテスト問題においてロバストな収束を示し、広範な適用可能性を示している。
• プリコンディショナのスペクトル近似の質が、観察された収束速度に直接寄与しており、コア設計原理の妥当性が裏付けられた。