基于积分核级数展开的多极边界元法及其截断误差分析-燕山大学学报.pdfVIP

37 卷 3 期 燕山大学学报 Vol. 37 No. 3 20 13 5 Journal of Yanshan University May 2013 年 月 文章编号：1007-791X (2013) 03-0254-06 基于积分核级数展开的多极边界元法及其截断误差分析 * 于春肖，苑润浩，于海源，王慧倩 （燕山大学 理学院，河北 秦皇岛 066004） 摘 要：对于二维Helmholtz 方程问题，本文提出一种基于积分核级数展开的多极边界元方法，推导证明了二 维 Helmholtz 方程的多极展开定理，给出了多极边界元法计算公式和计算过程，分析了截断误差，说明截断误 差可由截断项数 控制，并给出一个可广泛应用于实 计算的截断项数 的近似表达式。 关键词：多极边界元法；Helmholtz 方程；误差分析；截断项数 中图分类号：O241; O221.2 文献标识码：A DOI ：10.3969/j .issn.1007-791X .2013.03.013 0 引言 在文献 [4-11] 中所讨论的多极边界元法，按 照Helmholtz 方程基本解展开方式的不同，实 上 多极边界元法[1-3] (FM-BEM) 是近几年发展起 可 分为两类：一类是基于对角形式展开，适于求 来的一种快速数值计算方法，它融合了多极展开法 [5,7-11] [3,6] 解高频率问题 ；一类是基于级数形式展开 ， (FMM) 和边界元法 (BEM) 。这种方法的计算量和 适于求解低频率问题。在高频段使用级数形式展开 存储量都与未知量的个数 成正比，采用 GMRES 的多极边界元法，计算效率随着频率的升高而迅速 ( ) [4] 求解器 ，使得多极边界元方法的计算效率与 下降，这是由于频率的升高需要增加截断项数以保 传统边界元方法相比有显著的提高，在个人计算机 证计算精度；同样，在低频段也不能使用对角形式 上几小时内可实现数十万自由度的复杂工程计算。 多极边界元法，这就使得级数形式展开的多极边界 作为一种高效率的方法，多极边界元法适于求解大 元法值得研究。 规模科学与工程实 问题。 本文对于级数形式展开的多极边界元法进行 在科学与工程计算中，许多专家学者利用多极 了详细的阐述，并对这种方法的截断误差进行了分 边界元法研究二维Helmholtz 方程问题，主要研究 析。首先，简要介绍了Helmholtz 方程及其边界积 集中在声学领域：刘轶军