[論文レビュー] Integral equations requiring small numbers of Krylov-subspace iterations for two-dimensional penetrable scattering problems

Yassine Boubendir, Oscar P. Bruno+2 more|arXiv (Cornell University)|Oct 4, 2013

Electromagnetic Scattering and Analysis参考文献 21被引用数 4

ひとこと要約

本稿では、滑らかな境界を持つ2次元の電磁および音響散乱問題における透過性散乱体に対する、新しい正則化型結合境界積分方程式のクラスを提案する。正則化作用素（特に複素波数を用いた層状ポテンシャルまたはフーリエに基づく作用素）を適用することにより、変換問題が一意に解ける第二種フリードホルム方程式に変換され、与えられた精度においてKrylov部分空間反復回数と全体的な計算コストが顕著に削減される。

ABSTRACT

This paper presents a class of boundary integral equations for the solution of problems of electromagnetic and acoustic scattering by two dimensional homogeneous penetrable scatterers with smooth boundaries. The new integral equations, which, as is established in this paper, are uniquely solvable Fredholm equations of the second kind, result from representations of fields as combinations of single and double layer potentials acting on appropriately chosen regularizing operators. As demonstrated in this text by means of a variety of numerical examples (that resulted from a high-order Nystrom computational implementation of the new equations), these "regularized combined equations" can give rise to important reductions in computational costs, for a given accuracy, over those resulting from previous boundary integral formulations for transmission problems.

研究の動機と目的

• 透過性散乱問題における境界積分方程式のKrylov部分空間解法で一般的に見られる高反復回数の問題に取り組む。
• 一意に解ける第二種フリードホルム方程式としての境界積分式を構築し、スペクトルの条件数を改善する。
• 高精度なシミュレーションにおける収束に必要なKrylov反復回数を最小化することで、計算コストを低減する。
• 複素波数を用いた層状ポテンシャルおよびフーリエ記号に基づく正則化作用素を導入し、特異的（ハイパースィングュラー）作用素を安定化・正則化する。
• 滑らかな2次元透過性散乱体を対象とした高次Nyström離散化および数値実験を通じて、新規式の有効性を実証する。

提案手法

• 正則化作用素に作用する単一層および二重層ポテンシャルの組み合わせによる場の表現を構築し、積分方程式の安定化を図る。
• 複素波数（例：$ \tilde{S}_{\tilde{\nu}} $）を用いた正則化作用素を適用し、ハイパースィングュラー核をコンパクトで良好な条件数を持つ作用素に変換する。
• 複素波数層状ポテンシャルの高周波数漸近挙動と一致するフーリエ記号を持つフーリエに基づく正則化作用素を用いる。
• 適切なソボレフ空間における第二種フリードホルム方程式系として問題を定式化し、一意解の存在を保証する。
• グローバルな三角多項式近似、特異部と滑らか部へのカーネル分割、対数特異性に対する明示的 quadrature を用いた高次Nyström法を実装する。
• ベッセル関数およびハイパー記号関数の既知の漸近展開を活用し、正則化作用素の滑らかさ特性を分析する。

実験結果

リサーチクエスチョン

• RQ1 正則化作用素を用いて、透過問題におけるハイパースィングュラー積分作用素を、コンパクトで良好な条件数を持つ第二種方程式に変換できるか？
• RQ2 新規の正則化型結合積分方程式は、古典的式と比較して、Krylov部分空間反復回数を顕著に削減するか？
• RQ3 新規積分作用素のスペクトル的性質は何か？ そして、反復解法の収束速度にどのように影響するか？
• RQ4 複素波数層状ポテンシャルまたはフーリエ記号に基づく正則化作用素が、条件数が改善された第二種方程式を厳密に導出できるか？
• RQ5 滑らかな散乱体に対して、新規方程式の高次Nyström離散化は、最適収束率を達成し、計算コストを低減できるか？

主な発見

• 提案された正則化型結合積分方程式は、適切なソボレフ空間において一意に解ける第二種フリードホルム方程式であり、安定かつ収束する反復解法の実装を保証する。
• 複素波数層状ポテンシャルを正則化作用素として用いることで、システム行列のスペクトル的性質が改善され、Krylov部分空間反復回数が顕著に削減される。
• 数値実験により、新規式は古典的式と比較して、同じ精度においてしばしば2〜5倍の反復回数削減が達成され、全体的な計算コストが低減することが示された。
• パラメトリック単一層ポテンシャルとそのフーリエに基づく正則化形との差は、三階の滑らかさ作用素に相当し、正則化がスペクトル条件数の改善に有効であることを裏付けている。
• 高次Nyström法により、滑らかな境界に対してスペクトル精度が達成され、対数特異性に対する明示的 quadrature 法則により高精度が実現可能である。
• 正則化作用素が誘導する有利なスペクトル的性質のおかげで、特に高周波数散乱問題において、従来の手法と比較して優れた計算効率を達成している。