[論文レビュー] Corrected Trapezoidal Rules for Boundary Integral Equations in Three Dimensions

Bowei Wu, Per‐Gunnar Martinsson|arXiv (Cornell University)|Jul 6, 2020

Electromagnetic Scattering and Analysis参考文献 30被引用数 6

ひとこと要約

本稿では、一様グリッドでパラメータ化された3次元表面における境界積分方程式を解くために、補正付き台形則を導入し、係数行列の対角成分のみを修正することで $O(h^3)$ 収束を達成するとともに、9点ステンシルを用いることで $O(h^5)$ 収束を実現する。この手法は、標準的な台形則の単純さを保ちながら、特異核に対する高次精度を可能にする。

ABSTRACT

The manuscript describes a quadrature rule that is designed for the high order discretization of boundary integral equations (BIEs) using the Nyström method. The technique is designed for surfaces that can naturally be parameterized using a uniform grid on a rectangle, such as deformed tori, or channels with periodic boundary conditions. When a BIE on such a geometry is discretized using the Nyström method based on the Trapezoidal quadrature rule, the resulting scheme tends to converge only slowly, due to the singularity in the kernel function. The key finding of the manuscript is that the convergence order can be greatly improved by modifying only a very small number of elements in the coefficient matrix. Specifically, it is demonstrated that by correcting only the diagonal entries in the coefficient matrix, $O(h^{3})$ convergence can be attained for the single and double layer potentials associated with both the Laplace and the Helmholtz kernels. A nine-point correction stencil leads to an $O(h^5)$ scheme. The method proposed can be viewed as a generalization of the quadrature rule of Duan and Rokhlin, which was designed for the 2D Lippmann-Schwinger equation in the plane. The techniques proposed are supported by a rigorous error analysis that relies on Wigner-type limits involving the Epstein zeta function and its parametric derivatives.

研究の動機と目的

• 弱い特異核による3次元表面における境界積分方程式のNyström離散化における標準台形則の収束遅さを解消する。
• 台形則の単純さを保ちつつ、核の特異性を補正する高次 quadrature スキームを開発する。
• 対角成分を含む少数の行列要素のみを修正することで、$O(h^3)$ および $O(h^5)$ の高次収束を達成する。
• 2次元のDuan-Rokhlin quadrature則を周期的またはトーラス型幾何の3次元表面へ一般化する。
• 特異積分の漸近的挙動を解析するため、Wigner型極限およびEpsteinゼータ関数理論に基づく厳密な誤差解析を提供する。

提案手法

• 3次元パラメータ化表面（例：変形トーラスや周期的チャネル）上の一様グリッドに標準台形則を適用し、Nyström法を適用する。
• 非対角成分を $\mathbf{K}(i,j) = \mathcal{K}(\mathbf{x}_i,\mathbf{x}_j)w_j$ のまま維持し、核の特異性を補正するために行列の対角成分のみを修正する。
• 9点補正ステンシルを用いることで、$O(h^3)$ から $O(h^5)$ 収束に向上させ、計算コストの増加を最小限に抑える。
• 特異積分の漸近的挙動を解析するため、Epsteinゼータ関数およびそのパラメトリック導関数を含むWigner型極限を用いて誤差境界を導出する。
• 滑らかなカットオフ関数 $\phi$ と $\eta$ を用いた正則化技術を適用し、核を正則部と特異部に分解することで、誤差の精密な推定を可能にする。
• Fourierドメインで誤差を推定し、$S^{(\beta-j)}(\mathbf{u},s)$ を含む積分をバウンドし、$\eta^{(\ell)}(\mathbf{u}h)$ のテイラー展開を用いて剰余項を制御する。

実験結果

リサーチクエスチョン

• RQ1 標準台形則を用いたNyström法の収束速度を、弱い特異核を持つ3次元境界積分方程式に対して著しく改善できるか？
• RQ2 対角成分を含む少数の行列要素のみを修正することで、どの程度の高次精度が達成可能か？
• RQ3 提案された補正スキームは、Duan-Rokhlin法などの既存の2次元 quadrature 原則を3次元表面へどのように一般化できるか？
• RQ4 補正付き quadrature の理論的誤差境界は何か？また、メッシュ間隔 $h$ および核の滑らかさにどのように依存するか？
• RQ5 特異積分の文脈において、Epsteinゼータ関数およびその導関数といった特殊関数を用いて、誤差解析を厳密に正当化できるか？

主な発見

• Nyström行列の対角成分のみを補正することで、ラプラス核およびヘルムホルツ核の単層および双層ポテンシャルに対して $O(h^3)$ 収束が達成される。
• 9点補正ステンシルを用いることで収束速度が $O(h^5)$ に向上し、システム行列への修正を最小限に抑えつつ、著しい精度向上が達成される。
• 非対角成分については標準台形則の効率性を維持するため、FMM や高速直接ソルバーとの互換性が保たれる。
• 誤差解析はWigner型極限およびEpsteinゼータ関数に基づき、特異積分の収束挙動に対する厳密な基礎を提供する。
• $K$ 次補正に対して理論的誤差境界は $O(h^{2K+1})$ であり、カットオフ関数の導関数のテイラー展開により剰余項を明示的に制御できる。
• 強いかつ特異核に対しても、$\left|\sum_{|\mathbf{k}|\neq 0}\hat{f}_{\phi}(\mathbf{k},h) - Z_1(s)\right| \leq C h^{2K+1}$ という境界により、高次収束が実証されている。ここで $Z_1(s)$ はEpsteinゼータ関数の和である。