[论文解读] A novel fast solver for Poisson equation with the Neumann boundary condition

Zu‐Hui Ma, Weng Cho Chew+1 more|arXiv (Cornell University)|Jul 18, 2012

Electromagnetic Scattering and Analysis被引用 3

一句话总结

本文提出了一种针对二维泊松方程与诺伊曼边界条件的新型快速求解器，通过将求解过程分解为两个步骤：首先利用基于环路-树分解的快速树求解器求解电位移场 D，然后通过另一个快速求解器恢复电势 φ。该方法实现了 O(N) 复杂度，并相较于使用迭代求解器的传统有限元法（FEM），显著降低了求解时间，在数值示例中表现出高效率与高精度。

ABSTRACT

In this paper we present a novel fast method to solve Poisson equation in an arbitrary two dimensional region with Neumann boundary condition. The basic idea is to solve the original Poisson problem by a two-step procedure: the first one finds the electric displacement field $\mathbf{D}$ and the second one involves the solution of potential $ϕ$. The first step exploits loop-tree decomposition technique that has been applied widely in integral equations within the computational electromagnetics (CEM) community. We expand the electric displacement field in terms of a tree basis. Then, coefficients of the tree basis can be found by the fast tree solver in O(N) operations. Such obtained solution, however, fails to expand the exact field because the tree basis is not completely curl free. Despite of this, the accurate field could be retrieved by carrying out a procedure of divergence free field removal. Subsequently, potential distribution $ϕ$ can be found rapidly at the second stage with another fast approach of O(N) complexity. As a result, the method dramatically reduces solution time comparing with traditional FEM with iterative method. Numerical examples including electrostatic simulations are presented to demonstrate the efficiency of the proposed method.

研究动机与目标

• 开发一种快速、稳健且高效的数值方法，用于求解任意几何形状下的二维泊松方程诺伊曼边界条件问题。
• 克服传统有限元法（FEM）的局限性，如计算成本高、收敛速度慢，特别是在非规则区域中。
• 利用环路-树分解与快速树求解器，实现在静电学中泊松问题的近线性复杂度 O(N)。
• 通过数值示例验证该方法的准确性和效率，包括解析基准测试与复杂源分布。

提出的方法

• 该方法将求解过程分解为两个顺序步骤：首先求解电位移场 D 的 ∇·D = ρ，然后求解电势 φ 的 ∇φ = -D/ε。
• 电位移场 D 在树基下展开，其系数通过基于环路-树分解的快速树求解器以 O(N) 复杂度计算，该方法源自计算电磁学中的环路-树分解技术。
• 尽管树基并非完全无旋，但该方法通过迭代过程去除无散分量来纠正此问题。
• 电势 φ 通过在修正后的 D 场上应用第二个快速树求解器重建，从而保持整体 O(N) 复杂度。
• 无散场去除步骤涉及求解一个对称正定系统（环路基的格拉姆矩阵），使用共轭梯度法（CG）或广义最小残量法（GMRES）进行迭代求解。
• 通过施加参考电势以确保唯一性，满足诺伊曼问题的要求。

实验结果

研究问题

• RQ1 能否开发一种针对二维泊松方程诺伊曼边界条件的快速直接求解器，使其计算规模随问题规模线性增长？
• RQ2 如何将计算电磁学（CEM）中的环路-树分解技术适配到非规则区域中高效求解泊松方程？
• RQ3 与传统 FEM 搭配迭代求解器相比，该方法在诺伊曼问题上的计算性能提升程度如何？
• RQ4 当树基无法完全满足无旋条件时，该方法在多大程度上仍能保持精度？

主要发现

• 复杂度对比图显示，该方法的求解时间显著低于使用 CG 和 GMRES 求解器的传统 FEM。
• 新方法的计算复杂度趋近于 O(N)，主要开销来自无散场去除步骤，其性能取决于迭代求解器的收敛速度。
• 对于具有解析解的简单诺伊曼问题，模拟得到的电势分布与精确解高度一致，通过沿 y=0.1 的线性比较得到验证。
• 在复杂几何形状下的线源问题中，该方法生成了与 FEM 结果一致的有效电势分布，验证了其在非矩形区域中的鲁棒性。
• 在迭代阶段，GMRES 表现优于 CG，在相同停止准则（0.01）下收敛所需迭代次数更少。
• 该方法通过在某一点固定参考电势，成功处理了诺伊曼问题的非唯一性，确保了解的唯一性。