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

一句话总结

本文提出并比较了四种用于平面内光滑1D曲线上弱奇异积分方程Nyström离散化的高阶积分方法，重点针对对数奇异性。评估了全局周期性梯形法（Kapur–Rokhlin、Alpert、Kress）与基于分段的Gauss-Legendre方法（结合核函数修正），结果表明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.

研究动机与目标

• 开发并比较用于光滑1D曲线上对数奇异核积分方程Nyström离散化的高阶积分方法。
• 解决在弱奇异核的对角奇异性附近标准积分方法失效时实现高精度的挑战。
• 从收敛速率、精度饱和度以及与迭代求解器的适用性等方面评估不同方法的性能。
• 评估在实现复杂度、自适应性与迭代设置中谱性质之间的权衡。
• 为通过边界积分方程求解Laplace与Helmholtz问题提供实际的方案选择指导。

提出的方法

• 通过Kapur–Rokhlin、Alpert与Kress方法对周期性梯形规则进行专门校正，以处理对数奇异性。
• 通过核函数评估偏移与插值对Gauss-Legendre分段积分进行修正，以处理近似奇异积分，从而支持自适应加密。
• 对靠近对角线的矩阵条目应用高阶校正，以在存在对数奇异性时保持谱精度。
• 结合核函数的解析分解与数值积分，实现高阶收敛。
• 采用Nyström配点法，以节点作为配点位置，从积分方程构建线性系统。
• 在闭合与开放曲线上测试各方法，包括高频Helmholtz与轴对称Laplace问题，以评估其鲁棒性。

实验结果

研究问题

• RQ1 不同高阶积分方法在光滑曲线上弱奇异积分方程的收敛速率与精度饱和度方面如何比较？
• RQ2 核函数分解与积分类型对迭代求解器性能（尤其是GMRES收敛性）有何影响？
• RQ3 在奇异性存在时，自适应性如何影响基于分段与全局积分方法的性能？
• RQ4 修正后的Gauss-Legendre积分能否在支持局部加密的同时实现与全局方法相当的精度？
• RQ5 在实现复杂度、谱条件数与FMM兼容性方面，各方法的相对优缺点是什么？

主要发现

• Kress方法实现超代数收敛且饱和误差最低（10−13至10−15），但需要解析核函数分解，且不兼容快速多极方法。
• Alpert方法在收敛速度与单位节点精度上均优于Kapur–Rokhlin，相同点数下可多获得2至8位有效数字的精度。
• 修正后的高斯积分方法在1.5至2倍点数下可达到与10阶Alpert方法相当的精度，其主要优势在于支持自适应加密。
• 在高频条件下，16阶Alpert方法相比10阶Alpert可额外提升约3位有效数字精度，而Kapur–Rokhlin因谱条件数差，GMRES收敛显著变慢。
• 修正后的高斯方法天然适用于开放曲线，无需重参数化，而周期性方法在开放曲线上会退化至低阶精度。
• Kapur–Rokhlin方法虽实现最简单，但引入了虚假的大特征值，阻碍迭代收敛，尽管其精度为高阶，实际应用仍受限。