[論文レビュー] High-order accurate Nystrom discretization of integral equations with weakly singular kernels on smooth curves in the plane

Shifeng Hao, Alex H. Barnett+2 more|arXiv (Cornell University)|Dec 29, 2011

Electromagnetic Scattering and Analysis参考文献 20被引用数 29

ひとこと要約

本稿は、平面上の滑らかな1次元曲線上における弱特異積分方程式のNyström離散化のための4つの高次数則を提示し、比較している。主な焦点は対数特異性に向けられている。周期的台形則に基づくグローバル手法（Kapur–Rokhlin、Alpert、Kress）と、カーネル修正を施したパネルベースのガウス＝ルジャンドル則を評価し、Alpertの手法がノードあたりの精度が優れており、修正されたガウス則が適応的リファインメントを可能にすると示している。一方、Kressの手法は収束が最も速いが、FMMとの互換性に欠ける。

ABSTRACT

Boundary integral equations and Nystrom discretization provide a powerful tool for the solution of Laplace and Helmholtz boundary value problems. However, often a weakly-singular kernel arises, in which case specialized quadratures that modify the matrix entries near the diagonal are needed to reach a high accuracy. We describe the construction of four different quadratures which handle logarithmically-singular kernels. Only smooth boundaries are considered, but some of the techniques extend straightforwardly to the case of corners. Three are modifications of the global periodic trapezoid rule, due to Kapur-Rokhlin, to Alpert, and to Kress. The fourth is a modification to a quadrature based on Gauss-Legendre panels due to Kolm-Rokhlin; this formulation allows adaptivity. We compare in numerical experiments the convergence of the four schemes in various settings, including low- and high-frequency planar Helmholtz problems, and 3D axisymmetric Laplace problems. We also find striking differences in performance in an iterative setting. We summarize the relative advantages of the schemes.

研究の動機と目的

• 滑らかな1次元曲線上の対数特異核を持つ積分方程式のNyström離散化のための高次数則の開発と比較。
• 弱特異核における対角特異性に起因する標準的数値積分法の失敗を回避し、高精度を達成する挑戦に応える。
• 収束速度、精度の飽和、反復解法への適合性という観点から、異なる数則の性能を評価。
• 実装の複雑さ、適応的特性、反復解法環境におけるスペクトル特性のトレードオフを評価。
• 境界積分方程式を用いたラプラスおよびヘルムホルツ問題を解く際の数則選定に関する実用的指針を提供。

提案手法

• Kapur–Rokhlin、Alpert、Kressの手法を用いて、対数特異性に対処するための特別な補正を施した周期的台形則の応用。
• パネルベースのガウス＝ルジャンドル則を、カーネル評価のシフトと補間を用いて修正し、近似特異積分に対応させ、適応的リファインメントを可能にした。
• 対角付近の行列要素に対して高次補正を施し、対数特異性が存在してもスペクトル精度を維持。
• カーネルの解析的分解と数値積分の組み合わせにより、高次収束を達成。
• ノード点をコロケーション点として用いるNyströmコロケーション法を採用し、積分方程式から線形方程式系を構築。
• 閉曲線および開曲線の両方、特に高周波数ヘルムホルツ問題および軸対称ラプラス問題を対象とし、手法のロバストネスを評価。

実験結果

リサーチクエスチョン

• RQ1 弱特異積分方程式の滑らかな曲線上における、異なる高次数則の収束速度と精度飽和はどのように比較されるか？
• RQ2 カーネル分解と積分法の種別が、特にGMRES収束に与える影響は何か？
• RQ3 特異性が存在する状況下で、適応的リファインメントがパネルベースとグローバル数則の両方の性能に与える影響は？
• RQ4 修正されたガウス＝ルジャンドル則は、グローバル数則と同等の精度を達成しつつ、局所的リファインメントを可能にするか？
• RQ5 実装の複雑さ、スペクトル条件数、FMMとの互換性という観点から、各数則の相対的利点と欠点は何か？

主な発見

• Kressの数則は超代数的収束を達成し、最小の飽和誤差（10−13 から 10−15）を示すが、解析的カーネル分解を必要とし、FMMと互換性がない。
• Alpertの数則はKapur–Rokhlinを上回り、収束速度とノードあたりの精度の両面で優れている。同じ数の点数で2〜8桁の追加精度を達成。
• 修正されたガウス則は、10次数のAlpertと同等の精度を1.5〜2倍の点数で達成でき、適応的リファインメントを可能にするという主な利点を有する。
• 高周波数領域では、16次数Alpertは10次数Alpertに比べて最大3桁の追加精度を示すが、Kapur–Rokhlinはスペクトル条件数が悪いため、GMRES収束が著しく遅い。
• 修正されたガウス則は、再パrametrizationを必要とせず、開曲線を自然に扱える。これに対して周期的数則は開曲線では次数が著しく低下する。
• Kapur–Rokhlinは実装が最も簡単だが、反復解法の収束を妨げる不自然な大きな固有値を導入し、実用的利用に制限を受ける。