[论文解读] 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求积法处理正则部分。该方法应用于环面上的Laplace方程时，仅使用两个展开项，即可在计算量极小的情况下实现约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.

研究动机与目标

• 为高效离散化在三维曲面边界上具有奇异核的奇异积分算子提供解决方案。
• 通过采用解析奇点相减，减少对自适应求积或预计算表的依赖。
• 为环面区域上的Laplace方程开发一种简单且高阶的Nyström格式。
• 证明仅使用少量解析项（K=1）即可在不增加过多计算成本的前提下实现高精度。

提出的方法

• 该方法将奇异积分算子分解为奇异部分 $ D_{ij}^{K} $（通过高阶解析乘积积分处理）与正则部分 $ D_{ij}^{ullet} $（数值离散化）。
• 对于每个面片 $ \Gamma_{ij} $，利用 $ \Delta_{ij}^{k} $ 展开核函数，实现对主导奇异项的解析积分。
• 奇异项 $ D_{ijk} $ 通过递归计算公式（17）中的积分形式获得，结合多项式展开与Vandermonde线性系统求解。
• 采用三级策略自适应调整离散化：远距离点使用GL10，中距离点使用GL16并插值，近距离点采用奇点相减。
• 该格式采用复合张量积Gauss-Legendre求积，包含 $ p_1, p_2 $ 个面片，并将 $ \mu $ 插值至更精细的256点网格，以在奇点附近实现高精度评估。
• 选择参数 $ K=1 $，即仅使用奇点展开中的两项，以在精度与计算成本之间取得平衡。

实验结果

研究问题

• 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] 提出的混合方法，但提供了一种计算开销极小的高阶简洁替代方案。