[論文レビュー] A higher-order singularity subtraction technique for the discretization of singular integral operators on curved surfaces

Johan Helsing|arXiv (Cornell University)|Jan 30, 2013

Electromagnetic Scattering and Analysis参考文献 5被引用数 8

ひとこと要約

本稿では、曲がった表面における特異積分作用素の離散化のための高次の特異性除去技法を提示する。特異部には解析的積分積み上げ法を、正則部には標準的な Nyström 積分法を用いる。トーラス上でのラプラス方程式に適用したところ、最小限の計算負荷でわずか2つの展開項を用いて約10次収束性と10桁未満の精度を達成した。

ABSTRACT

This note is about promoting singularity subtraction as a helpful tool in the discretization of singular integral operators on curved surfaces. Singular and nearly singular kernels are expanded in series whose terms are integrated on parametrically rectangular regions using high-order product integration, thereby reducing the need for spatial adaptivity and precomputed weights. A simple scheme is presented and an application to the interior Dirichlet Laplace problem on some tori gives around ten digit accurate results using only two expansion terms and a modest programming- and computational effort.

研究の動機と目的

• 曲がった3次元表面における特異核を有する特異積分作用素を、効率的に離散化する課題に対処すること。
• 適応的積分や事前計算されたテーブルに依存することを減らすために、解析的特異性除去を用いること。
• トーラス領域におけるラプラス方程式のための単純で高次の Nyström スキームを開発すること。
• わずかな数の解析的項（K=1）が、過度な計算コストを伴わずに高精度を達成できることを示すこと。

提案手法

• 本手法は、特異積分作用素を特異部 $ D_{ij}^{K} $ と正則部 $ D_{ij}^\bullet $ に分割する。特異部は高次の解析的積分積み上げ法で処理し、正則部は数値的に離散化する。
• 各パッチ $ \Gamma_{ij} $ に対して、$ \Delta_{ij}^{k} $ を用いた級数展開によりカーネルを展開し、主な特異項の解析的積分を可能にする。
• 特異項 $ D_{ijk} $ は、式 (17) に示す形の積分を多項式展開とバーデモンド行列を用いて再帰的に評価することで計算する。
• 三段階戦略により離散化を適応させる：遠方の点には GL10、中距離の点には GL16 と補間を用い、近接点には特異性除去を適用する。
• 本手法は、$ p_1, p_2 $ パッチを用いた合成テンソル積のガウス＝レジェンドル積分法を用い、$ \mu $ をより細かい256点グリッドに補間して特異点付近の評価を高精度に保つ。
• パラメータ $ K=1 $ が選ばれ、特異性展開に2つの項しか使用しない。これは精度と計算コストのバランスを取るためである。

実験結果

リサーチクエスチョン

• RQ1 高次の特異性除去技法は、適応的メッシュ細分化や事前計算された重みに依存せずに、曲がった表面における特異積分作用素の離散化で高精度を達成できるか？
• RQ2 メッシュ解像度を高めた場合、とくにほぼ特異核が存在する状況下で、この手法の収束速度はどのように振る舞うか？
• RQ3 特異性除去スキームにおける精度と計算コストのバランスを最適化するための展開項数 $ K $ の最適値は何か？
• RQ4 本手法は、より複雑なハイブリッドソルバーと比較して、精度と効率の面でどの程度競争力を持つのか？

主な発見

• 図2に示すように、解 $ U(r) $ の相対的 $ L^2 $ 誤差において、トーラス中心線に沿って約10次収束性を達成した。
• $ \delta = (0, 0.25) $ の場合、10,000点の特異性除去ステップに要する時間はたったの28秒であり、計算効率が優れていることが示された。
• わずか $ K=1 $ の展開項と、控えめな実装コストで、解の精度が10桁未満に達した。
• メッシュを細分化するに従い、正則部 $ D_{ij}^\circ $ の離散化誤差が徐々に支配的になり、理論的な10次収束を下回るようになる。
• 本手法は、BremerとGimbutas [1] のハイブリッド手法にはまだ及ばないが、最小限の計算オーバーヘッドで単純かつ高次の代替手法を提供する。