[论文解读] 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

一句话总结

本文提出了一种校正的梯形求积规则，用于在由均匀网格参数化的三维表面上求解边界积分方程，通过仅修改系数矩阵的对角线元素，实现 $O(h^3)$ 收敛；若采用九点模板，则可达到 $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.

研究动机与目标

• 解决标准梯形求积在三维表面边界积分方程 Nyström 离散化中因弱奇性核导致的收敛缓慢问题。
• 开发一种高阶求积格式，在保持梯形规则简单性的同时校正核的奇异性。
• 通过仅修改少量矩阵元素（特别是对角线项）实现高阶收敛（$O(h^3)$ 和 $O(h^5)$）。
• 将二维 Duan-Rokhlin 求积规则推广至具有周期性或环面几何结构的三维表面。
• 基于 Wigner 型极限和 Epstein zeta 函数理论，为奇异积分提供严格的误差分析。

提出的方法

• 在参数化三维表面（如变形环面或周期性通道）的均匀网格上应用 Nyström 方法与标准梯形求积。
• 通过仅修改对角线项来校正系数矩阵，以考虑核的奇异性，同时保持非对角线项为 $\mathbf{K}(i,j) = \mathcal{K}(\mathbf{x}_i,\mathbf{x}_j)w_j$。
• 采用九点校正模板以实现 $O(h^5)$ 收敛，而非 $O(h^3)$，在计算开销几乎不变的前提下显著提升精度。
• 利用涉及 Epstein zeta 函数及其参数导数的 Wigner 型极限推导误差界，分析奇异积分的渐近行为。
• 采用平滑截断函数 $\phi$ 和 $\eta$ 的正则化技术，将核分解为正则部分与奇异部分，以实现精确的误差估计。
• 通过有界涉及 $S^{(\beta-j)}(\mathbf{u},s)$ 的积分并在傅里叶域中估计误差，利用 $\eta^{(\ell)}(\mathbf{u}h)$ 的泰勒展开控制余项。

实验结果

研究问题

• RQ1 能否显著提升标准梯形求积在三维边界积分方程中因弱奇性核导致的 Nyström 方法收敛速率？
• RQ2 通过仅修改系数矩阵中极少数元素（特别是对角线项）能在多大程度上实现高阶精度？
• RQ3 所提出的校正方案在多大程度上可将现有二维求积规则（如 Duan-Rokhlin 方法）推广至三维曲面？
• RQ4 校正求积规则的理论误差界是什么？其如何依赖于网格间距 $h$ 和核的光滑性？
• RQ5 能否基于特殊函数（如 Epstein zeta 函数及其导数）在奇异积分背景下严格证明误差分析的合理性？

主要发现

• 仅校正 Nyström 矩阵的对角线元素，即可实现拉普拉斯和赫姆霍兹核的单层与双层势函数的 $O(h^3)$ 收敛。
• 采用九点校正模板可将收敛率提升至 $O(h^5)$，在对系统矩阵的修改极小的情况下显著提高精度。
• 该方法保持了标准梯形求积对非对角线元素的高效性，从而与快速多极方法（FMM）和快速直接求解器兼容。
• 误差分析基于 Wigner 型极限和 Epstein zeta 函数，为奇异积分的收敛行为提供了严格的理论基础。
• 对于 $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 zeta 函数和。