[论文解读] Development of a BEM Solver using Exact Expressions for Computing the Influence of Singularity Distributions

Supratik Mukhopadhyay, Nayana Majumdar|ArXiv.org|Apr 5, 2006

Electromagnetic Scattering and Analysis被引用 4

一句话总结

本文提出了一种快速且精确的边界元法（BEM）求解器，该求解器使用平直3D表面上均匀分布源的势场和力场的精确解析表达式。通过消除数值积分误差并精确捕捉近场奇异性，该求解器实现了对单位正方形板（0.3667869）和单位立方体（0.6606746）的高精度电容计算，其结果与现有最先进文献值高度一致，验证了其准确性。

ABSTRACT

Closed form expressions for three-dimensional potential and force field due to singularities distributed on a finite flat surface have been used to develop a fast and accurate BEM solver. The exact expressions have been investigated in detail and found to be valid throughout the complete physical domain. Thus, it has been possible to precisely simulate the complicated characteristics of the potential and force field in the near-field domain. The BEM solver has been used to compute the capacitance of a unit square plate and a unit cube. The obtained values have compared well with very accurate values of the capacitances available in the literature.

研究动机与目标

• 开发一种边界元法求解器，以消除边界元法中离散源近似带来的数值误差。
• 准确模拟平直表面附近具有奇点的复杂近场势场和力场。
• 计算经典静电问题——单位正方形板和单位立方体——的高精度电容值，这些问题长期以来被视为挑战。
• 在整个物理域内（包括存在不连续性和陡峭梯度的区域）验证精确表达式的有效性。

提出的方法

• 求解器使用矩形平直表面均匀源分布产生的势和场分量的闭式解析表达式。
• 通过包含几何项的对数和反双曲正切函数的精确积分表达式计算势。
• 场分量（E_x, E_y, E_z）通过势的精确偏导数求得。
• 由于影响函数具有解析性质，该方法可无缝处理所有边界条件（狄利克雷、诺伊曼、罗宾）。
• 在边和角附近采用非均匀的网格加密策略，以提高高电荷密度区域的分辨率。
• 通过与高度精确的文献值及不同离散化水平下的数值基准进行对比，对求解器进行了验证。

实验结果

研究问题

• RQ1 平直表面上均匀分布源的影响是否可通过精确解析表达式完全消除BEM模拟中的数值误差？
• RQ2 该新型BEM求解器在近场区域对边缘和角附近势场和场强的奇异性行为的再现精度如何？
• RQ3 采用此精确BEM方法计算单位正方形板和单位立方体电容的可实现精度是多少？
• RQ4 随着网格细化程度的提高，该精确求解器与传统BEM的性能表现相比如何？

主要发现

• 势和场分量的精确表达式在整个物理域内均有效，包括源表面本身及其附近区域。
• 求解器成功再现了边缘和角附近的尖锐场不连续性和奇异性，而点源近似则难以准确捕捉这些特征。
• 单位正方形板的计算电容为0.3667869，与最精确的文献值（例如，Read等人给出的0.3667874 ± 1×10⁻⁷）高度一致。
• 单位立方体的计算电容为0.6606746，与最佳文献结果非常接近（例如，Read等人给出的0.6606785 ± 6×10⁻⁷，Mascagni给出的0.6606780 ± 2.7×10⁻⁷）。
• 即使采用高达100×100个单元的数值离散化，E_x和表面附近的势仍表现出显著振荡，而精确求解器则无此类伪影。
• 单位立方体顶面的电荷密度在靠近边缘和角附近急剧增加，该求解器通过边界附近高分辨率的网格划分正确捕捉了这一现象。