[论文解读] An Integral Equation Formulation of the $N$-body Dielectric Spheres Problem. Part II: Complexity Analysis

Bérenger Bramas, Muhammad Hassan+1 more|arXiv (Cornell University)|Nov 17, 2019

Electromagnetic Scattering and Analysis参考文献 83被引用 4

一句话总结

本文针对N个介电球体问题，对伽辽金边界积分方程方法进行了复杂度分析，证明了求解线性系统所需的GMRES和CG迭代次数与N无关。结合快速多极子方法（FMM），这确保了计算成本的线性增长，从而确立了该方法在大规模极化介电粒子电势模拟中的线性精度缩放特性。

ABSTRACT

This article is the second in a series of two papers concerning the mathematical study of a boundary integral equation of the second kind that describes the interaction of $N$ dielectric spherical particles undergoing mutual polarisation. The first article presented the numerical analysis of the Galerkin method used to solve this boundary integral equation and derived $N$-independent convergence rates for the induced surface charges and total electrostatic energy. The current article will focus on computational aspects of the algorithm. We provide a convergence analysis of the iterative method used to solve the underlying linear system and show that the number of liner solver iterations required to obtain a solution is independent of $N$. Additionally, we present two linear scaling solution strategies for the computation of the approximate induced surface charges. Finally, we consider a series of numerical experiments designed to validate our theoretical results and explore the dependence of the numerical errors and computational cost of solving the underlying linear system on different system parameters.

研究动机与目标

• 建立N个介电球体的伽辽金边界积分方程方法在计算成本上的线性缩放特性，确保大N下的计算效率。
• 证明求解线性系统所需的GMRES和CG迭代次数与粒子数N无关。
• 通过在不同系统参数下进行的数值实验，验证理论推导的收敛界。
• 证明该方法在实现O(N)计算复杂度的同时，保持了与N无关的误差稳定性。
• 为将该方法扩展至计算静电作用力及其他超越拉普拉斯方程的算子奠定基础。

提出的方法

• 该方法采用针对N个介电球体在均匀介质中的第二类边界积分方程进行伽辽金离散化。
• 求解得到的线性系统采用迭代求解器——GMRES，以及在对称化后采用共轭梯度（CG）方法。
• 使用快速多极子方法（FMM）在O(N)次操作内完成矩阵-向量乘积，从而实现可扩展的线性系统求解。
• 对GMRES和CG的收敛性进行了分析，证明了在给定容差下，迭代次数存在与N无关的上界。
• 对线性系统进行对称化，以允许使用CG方法，这提高了效率并支持并行收敛性分析。
• 在笔记本电脑上对不同N和ℓ_max（5和10）进行了数值实验，求解器容差分别设为10^−6和10^−9，以验证可扩展性。

实验结果

研究问题

• RQ1 求解线性系统所需的GMRES迭代次数是否会随介电球体数量N的增加而增长？
• RQ2 共轭梯度（CG）方法能否有效应用于伽辽金离散化产生的对称化线性系统，且其迭代次数与N无关？
• RQ3 在使用FMM和与N无关的迭代收敛性条件下，该算法的整体计算成本是否与N呈线性关系？
• RQ4 在实际中，离散化误差与FMM误差如何相互作用，且FMM误差是否主导了解的精度？
• RQ5 该方法能否扩展至计算介电粒子间的静电作用力，同时保持相似的可扩展性和精度特性？

主要发现

• 达到给定相对误差容差所需的GMRES迭代次数与N无关，证明了迭代求解器的N无关收敛性。
• 求解对称化线性系统所需的CG迭代次数同样与N无关，其收敛界在定性上与GMRES相似。
• 数值实验确认了GMRES和CG方法在计算时间上均表现出线性缩放行为，且在不同N和ℓ_max下性能几乎一致。
• 在测试的几何构型中，FMM误差并未主导离散化误差，支持了渐近复杂度分析的有效性。
• N-误差稳定的伽辽金离散化与使用FMM的O(N)迭代求解器相结合，确立了该方法在精度上的线性缩放特性。
• 结果支持该框架未来可扩展用于计算静电作用力，且具备相似的可扩展性，如未来工作所计划。