[论文解读] Fast solution of boundary integral equations with the generalized Neumann kernel

Mohamed M. S. Nasser|arXiv (Cornell University)|Aug 24, 2013

Numerical methods in engineering参考文献 53被引用 41

一句话总结

该论文提出了一种快速、精确的方法，利用梯形法则的Nyström方法结合快速多极子方法（FMM）与FFT，求解具有广义Neumann核及其伴随核的边界积分方程。该方法对广义Neumann核的复杂度为$O((m+1)nackslashln n)$，对伴随核的复杂度为$O((m+1)n)$，从而可在无需数值微分的情况下，高效求解高连通性和几何复杂的区域问题。

ABSTRACT

A fast method for solving boundary integral equations with the generalized Neumann kernel and the adjoint generalized Neumann kernel is presented. The method is based on discretizing the integral equations by the Nyström method with the trapezoidal rule to obtain $(m+1)n imes(m+1)n$ linear systems where $m+1$ is the multiplicity of the multiply connected domain and $n$ is the number of nodes in the discretization of each boundary component. The obtained linear systems are solved by the generalized minimal residual (GMRES) method. Each iteration of the GMRES method requires a matrix-vector product which can be computed using the Fast Multipole Method (FMM). The complexity of the presented method is $O((m+1)n\ln n)$ for the integral equation with the generalized Neumann kernel and $O((m+1)n)$ for the integral equation with the adjoint generalized Neumann kernel. The presented numerical results illustrate that the presented method gives accurate results even for domains with high connectivity, domains with piecewise smooth boundaries, and domains with close boundaries.

研究动机与目标

• 开发一种用于在多连通区域中求解具有广义Neumann核的边界积分方程的快速且稳定的数值方法。
• 通过重新表述奇异积分算子$\mathbf{M}$的离散化，消除对未知函数进行数值微分的需求。
• 通过结合FMM与FFT，对Nyström离散化产生的稠密线性系统实现最优计算复杂度。
• 为具有高连通性、分段光滑边界及近距离边界的区域提供鲁棒的数值解。
• 通过MATLAB函数实现并验证该方法，用于共形映射与边界值问题。

提出的方法

• 在每个边界分量上使用梯形法则的Nyström方法对广义Neumann核及其伴随核方程进行离散化。
• 将离散化算子$\mathbf{M}$重写为适合FMM的矩阵与适合FFT的循环矩阵块之和。
• 使用FMM在$O((m+1)n)$次操作内计算涉及非循环部分的矩阵-向量乘积。
• 使用FFT在$O((m+1)nackslashln n)$次操作内计算涉及循环部分的矩阵-向量乘积。
• 结合FMM用于矩阵-向量乘积，使用GMRES迭代法求解所得稠密线性系统。
• 实现两个MATLAB函数：FBIE用于广义Neumann核，FBIEad用于伴随核，以及FCAU用于快速Cauchy积分计算。

实验结果

研究问题

• RQ1 能否在不损失精度的前提下，将求解具有广义Neumann核的边界积分方程的计算成本降低至$O((m+1)^2n^2)$以下？
• RQ2 在奇异积分算子$\mathbf{M}$的离散化中，如何避免对未知函数进行数值微分？
• RQ3 边界几何特性（如高连通性、分段光滑性或近距离）对解的收敛性与稳定性有何影响？
• RQ4 为何FBIE的精度依赖于函数$\theta$是否为常数，而FBIEad的精度却不受影响？
• RQ5 FMM与FFT的结合能否为广义Neumann核及其伴随核均实现最优复杂度？

主要发现

• 所提方法对广义Neumann核方程的求解复杂度为$O((m+1)nackslashln n)$，对伴随核的复杂度为$O((m+1)n)$，显著优于标准稠密求解器。
• 该方法在连通数超过1,000、具有分段光滑边界或边界间距小至$10^{-4}$的区域中仍保持高精度。
• 对于FBIE函数，当$\theta$非常数且分离距离$\varepsilon$极小时，精度会下降；而FBIEad的精度则不受$\theta$的影响，保持一致。
• 随着$\varepsilon$减小，GMRES迭代次数、CPU时间和条件数均增加，但因采用奇异性减去法，仍保持可管理状态。
• 系数矩阵的条件数随$\varepsilon$减小而增加，但基于FMM的求解器仍保持鲁棒性。
• 数值证据表明，算子$\mathbf{M}$的离散化是导致FBIE精度依赖$\theta$的原因，而FBIEad中更简单的Fredholm算子则无此依赖。