[论文解读] Use of Triangular Elements for Nearly Exact BEM Solutions

Supratik Mukherjee, Nayana Majumdar|ArXiv.org|Apr 12, 2007

Electromagnetic Scattering and Analysis参考文献 24被引用 6

一句话总结

本文提出了一种新型边界元法（BEM）求解器 neBEM，该求解器采用三角形单元上均匀分布奇点的势和通量的精确解析表达式，可在几何奇点（如角点）附近实现高度精确且稳定的解。该方法通过建模分布奇点而非节点奇点，避免了传统BEM中常见的数值边界层和振荡现象，实现了对单位正方形板电容的估计值为 0.3660587 pF/(4πε₀)，其精度与当前最先进的方法相当或更优。

ABSTRACT

A library of C functions yielding exact solutions of potential and flux influences due to uniform surface distribution of singularities on flat triangular and rectangular elements has been developed. This library, ISLES, has been used to develop the neBEM solver that is both precise and fast in solving a wide range of problems of scientific and technological interest. Here we present the exact expressions proposed for computing the influence of uniform singularity distributions on triangular elements and illustrate their accuracy. We also present a study concerning the time taken to evaluate these long and complicated expressions extit{vis a vis} that spent in carrying out simple quadratures. Finally, we solve a classic benchmark problem in electrostatics, namely, estimation of the capacitance of a unit square plate raised to unit volt. For this problem, we present the estimated values of capacitance and compare them successfully with some of the most accurate results available in the literature. In addition, we present the variation of the charge density close to the corner of the plate for various degrees of discretization. The variations are found to be smooth and converging. This is in clear contrast to the criticism commonly leveled against usual BEM solvers.

研究动机与目标

• 开发一种边界元法，以克服传统BEM求解器在几何奇点（如角点和边）附近常见的数值不稳定性和不准确性。
• 建立平面上三角形和矩形单元上均匀奇点分布影响函数的精确闭式表达式库。
• 证明使用分布奇点而非节点奇点的BEM方法在实现物理一致且收敛结果方面的可行性与优越性。
• 在经典静电基准问题——单位正方形板电容问题上验证该方法，实现无需外推或对称性假设的高精度结果。
• 表明使用精确影响函数计算的导体板角点附近电荷密度随网格细化可实现平滑且一致的收敛，与标准BEM中常见的振荡行为形成鲜明对比。

提出的方法

• ISLES库提供了基于平面上三角形和矩形单元上均匀奇点分布的势和通量的精确解析表达式，其推导基于平方反比定律。
• neBEM求解器利用这些精确表达式计算影响矩阵，无需数值积分，从而消除了积分误差，显著提升了近场区域的精度。
• 边界条件在单元质心处施加，但该方法允许在域内任意点（包括边缘和角点）精确计算势和通量。
• 该方法采用三角形与矩形单元的混合网格，其中在必要时使用三角形单元以建模任意三维几何形状，特别是在角点附近。
• 采用高精度数值积分作为基准，验证整个物理域（包括近场区域）中解析表达式的正确性。
• 将求解器应用于单位正方形板问题，通过不同网格细化程度研究角点附近电荷密度的收敛性。

实验结果

研究问题

• RQ1 在三角形单元上使用影响函数的精确解析表达式，是否能显著提升BEM在几何奇点（如角点）附近的精度与稳定性？
• RQ2 基于分布奇点的BEM求解器在边和角附近收敛性与数值振荡方面，与传统基于节点奇点的BEM相比表现如何？
• RQ3 neBEM求解器在无需使用对称性或外推法的前提下，能在单位正方形板电容等基准静电问题上达到多高的精度？
• RQ4 使用精确影响函数计算导体板角点附近的电荷密度时，其行为如何？是否随网格细化实现平滑收敛？
• RQ5 neBEM方法在角奇点问题上的精度是否可与蒙特卡洛或奇异摄动等先进方法相媲美或超越？

主要发现

• neBEM求解器对单位正方形板的电容计算结果为 0.3660587 pF/(4πε₀)，与文献中最精确的结果（如经外推的 0.3667874±1×10⁻⁷）高度一致。
• 单位正方形板角点附近的电荷密度随网格细化表现出平滑且单调的收敛性，无任何振荡，与传统BEM中常见的不规则行为形成鲜明对比。
• 对角点附近电荷密度进行最小二乘拟合后得到的斜率为 0.713567，与旧有及近期研究公认的 0.7034 值非常接近，表明具有高度的物理一致性。
• ISLES库的精确表达式通过与高精度数值积分在整个域内对比得到验证，确认其正确性与鲁棒性，即使在关键的近场区域亦然。
• 该方法表明，当使用精确影响函数时，在单元质心处施加边界条件已足够，且势和通量可在任意点（包括边和角）实现高精度计算。
• 本研究证实，使用分布奇点替代节点奇点可消除数值边界层，实现具有物理真实性的、收敛的解，适用于具有尖锐边缘和角点的问题。