两种格林函数快速多极子求解方法.PDFVIP

第2卷 第1期 16 中国科技论文在线 SCIENCEPAPER ONLINE 2007年 1月 两种格林函数的快速多极子求解方法 宁德志，滕 斌，勾 莹 (大连理工大学海岸和近海工程国家重点实验室，大连 116023) 摘 要：高阶边界元法在很多工程计算中得到了广泛的应用，但由于需要多次求解格林函数及其导数，高昂的 计算量和存储量使其很难应用于大型工程问题。本文采用快速多极子方法，对于形如 1/r 的简单格林函数采用 球坐标系下的双谐展开，对于满足自由水面条件的复杂格林函数则在柱坐标系下展开，使其计算量和存储量都 由未知量的平方量级降为未知量的线性量级。对于高阶边界元法中的固角系数和柯西主值积分本文采用了直接 求解的方法。通过对无限区域中水流对圆球的绕射和波浪与漂浮方箱的数值计算，表明对于大中型计算问题 FMM 算法更加有效。 关键词：高阶边界元法；快速多极子方法；格林函数 中图分类号：O353.2 文献标识码：A 文章编号：1673－7180(2007)01 －0016 －5 0 引言 解波浪与结构物相互作用问题，但格林函数采用 边界积分方程法也叫边界元法（BEM ），具 的是Rankine源，多极子展开在球坐标系下进行。 有降低问题的维数和计算精度较高这两大显著特 此外在波浪与结构物的频域分析中，满足自由水 点，作为有限元法、有限差分法等域内方法的发 面和无穷远边界条件的格林函数应用也十分广 展与补充，在处理开区域、非线性、非定常、耦 泛，但其形式较为复杂，无法像Rankine源那样在 合等复杂问题中，有其独到的长处。边界元法虽 球 坐 标 系 中 做 多 极 子 展 开 。 Utsunomiya 和 然将问题降维，但系数矩阵为满阵存储，因此存 Watanabe[5]提出了一种在柱坐标系下做多极子展 储量级和计算量级都为未知量 N 的平方量级，当 开的方法，利用贝塞尔函数的Graf加法定理对满 未知量个数较多时，计算量和存储量以平方量级 足自由水面条件的复杂格林函数做多极子展开， 递增，这成了边界元法应用的一个瓶颈。 并应用于高阶边界元方法中。 本文采用多极子展开方法来降低边界元法的 在高阶边界元方法中，一个主要的难点是固 计算量和存储量，多极子展开法（Fast Multipole 角系数和柯西主值奇异积分的计算，前人多采用 Method ）是最近发展起来的一种快速的数值计算 增加一个内部积分的间接方法消除[6] ，这对多极 方法，最早由Rokhlin[1]提出，之后Greengard[2]应 子展开方法带来了许多不便，降低了多极子展开 用到二维和三维多物体的库仑势问题以及其他领 方法的高效性。 域中。多极子展开法和迭代法联合应用，不形成 本文对于固角系数和柯西主值积分采用直接 系数矩阵A ，减少了矩阵与向量相乘计算的复杂 方法计算，并对前面提到的两种格林函数的快速 性，因此降低了计算量和存储量。当未知量个数 多极子方法都进行了分析计算，建立了低存储、 很大时该方法就显示了它的优越性，详细理论可 高效率的波浪与结构物作用的分析程序。最后， 见Greengard[2] ，Ken-ichi[3]等的文章。 分别以水流对无界区域中的固定圆球绕射和波浪 在过去不到二十年中，多级子方法被应用到 与漂浮方箱作用为例，通过与解析解和传统高