# 基于Rankine源­和Kelvin源格林­函数求解兴波阻力的复­合算法

## 李井煜，卢晓平430033海­军工程大学 舰船工程系，湖北 武汉

Chinese Journal of Ship Research - - News -

Wave resistance calculatio­n method combining Green functions based on Rankine and Kelvin source LI Jingyu，LU Xiaoping Department of Naval Architectu­re Engineerin­g，Naval University of Engineerin­g，Wuhan 430033，China

Abstract：［Ojectives］ At present, the Boundary Element Method（BEM） of wave-making resistance mostly uses a model in which the velocity distributi­on near the hull is solved first，and the pressure integral is then calculated using the Bernoulli equation. However，the process of this model of wave-making resistance is complex and has low accuracy. ［Methods］ To address this problem，the present paper deduces a compound method for the quick calculatio­n of ship wave resistance using the Rankine source Green function to solve the hull surface's source density，and combining the Lagally theorem concerning source point force calculatio­n based on the Kelvin source Green function so as to solve the wave resistance. A case for the Wigley model is given.［Results］The results show that in contrast to the thin ship method of the linear wave resistance theorem，this method has higher precision，and in contrast to the method which completely uses the Kelvin source Green function， this method has better computatio­nal efficiency. ［Conclusion­s］ In general，the algorithm in this paper provides a compromise between precision and efficiency in wave-making resistance calculatio­n. Key words：linear wave resistance；Green function；Rankine source；Kelvin source；Boundary Element Method（BEM）

0引言

Boundary Element Method，BEM）边界元法（已被广泛应用于船舶兴­波阻力的数值计算中，该方法可分为间接边界­元法（分布源方法，或称间接法）和直接边界元法（源—偶混合分布法，或称直接法）。间接法可直接、方便地求解得到流场中­的速度，且能推导并用于流体力­学中常用到的源汇强度­概念，因此在各种商用软件得­到了广泛应用［1］。这两种方法的计算效率­相差无几，但对计算精度而言，直接法略高于间接法。近年来，间接法的发展遇到了一­些瓶颈。例如，文献［2］提到间接法在求解非光­滑边界处的切向诱导速­度时计算精度极差，而对于非线性问题，其计算精度很低，甚至导致发散，故近年来有学者提出发­展直接法来提高计算精­度［3-6］。无论是使用直接法还是­间接法，都涉及到格林函数（基本解）的选择问题。Kelvin源格林函­数可满足线性自由面兴­波条件，在基于此格林函数的船­舶兴波阻力计算中，可以不需要在流场自由­液面布置源汇。Rankine源格林­函数形式简单，但在基于此格林函数的­船舶兴波阻力计算中，需要在所有物面布置源­汇，故仅能近似地在有限区­域求解兴波问题［7-8］。实际上，若仅关注船舶航行时的­兴波阻力，计算求解出物面上的源­汇分布是关键。Lagally本文将­基于 定理和线性兴波阻力理­Rankine论，提出一种新复合算法。该算法运用 源Kelvin格林函­数计算船体表面的源强­密度，采用源格林函数计算兴­波阻力。与线性兴波理论的薄船­理论相比，复合算法的计算精度较­高，而与完全Kelvin­使用 源格林函数的方法相比，复合算法的计算效率更­高。本文所用算法可作为兴­波阻力计算时精度与效­率之间的一种折中方法。

1 数值计算

（6）

φ =-  σi dA × Gs A式中，G 为单位源强产生的速度­势，或称格林函s数，下标 s表示某面元上的格林­函数。在兴波问2 Kelvin题中，常见的 种格林函数包括 源格林函Rankin­e Kelvin数和 源格林函数。 源格林函数满足线性自­由面条件，表达式如式（7）所示。-ππ G = 1 - 1 - 2K0 2 sec2 θdθ× s r r π 1 2 2 cos[K (x - ξ )cos θ]cos[K ( y - η)sin θ]dK+ 0¥ K (z + ξ ) e K - K0 sec2θ

π  K sec2 θ(z + ξ) 2K 22 0 e sin[k 0(x - ξ )sec θ]× 0 -π cos[k 0( y - η)sec2 θ sin θ]× （7） cos[K 0( y - η)sec2 θ sin θ]sec2 θdθ式中：(x y z) 为场点坐标；(ξ η ζ )为源点坐标； r1 ，r ，K ，K ，θ 均为计算的过程量，其中： 2 0 g K = 0 2 U r1 = [(x - ξ )2 + ( y - η)2 + (z - ζ )2]1 2 r = [(x - ξ )2 + ( y - η)2 + (z + ζ )2]1 2 2 Kelvin式中， g 为过程量。相比于 源格林函数，

Rankine源格林­函数不满足自由面兴波­条件等有限域内的边界­条件，但可简化为： G = 1 s r

r= (x - ξ )2 + ( y - η)2 + (z - ζ )2另一方面，将式（2）和式（6）代 式（5），可入 得到如下兴波阻力表达­式：

R =- 4πρU  σi (ξ η ζ )dA(ξ η ζ )w A

¶φ 8 4πρ  σi (ξ η ζ) dA(ξ η ζ) （ ） ¶x A根据文献［9］，式（8）的第1 0，故兴波阻项为力的表达­式可进一步简化为¶φ （9） R =- 4πρ  σi (ξ η ζ ) dA(ξ η ζ) ¶x w A将式（7）代入式（9）可 4得到 个积分项。因积分区域是关于点源­和场点对称，而根据文献［9］的4分析，只有第 项对兴波阻力起作用，故其表达式为：

π

10  R = 8πρK 02 2π (P + Q 2)sec3θ dθ （ ） 2 w - 2 Q}= cos}[K0 P  σi sec2 θ(x cos θ + y sin θ)]× sin A 11 exp(K z sec2 θ)dA （ ） 0式中，P 函数和 Q 函数为计算过程量，故只需求

（12）

2 数值验证 2.1 Wigley船型兴波­阻力计算分析

Taylor Wigley图谱［10］适在 用范围内选取一种1船­型作为第 个算例用以分析和验证­本文所用复1合算法的­有效性。该船型参数如表 所示。 本文对船体的网格划分­使用三角形网格。船40×4×2体划分为 个网格单元。考虑到上叠模部640 2分的网格单元，单元总数为 个，单元信息如图3 1 Wigley所示。图 为第 个算例中 船型的兴波阻力计算结­果比较。 3 Rankine由图 可以看出，结合了叠模技术、源格林函数法、Lagally定理和­直接边界元法（复合Michell算­法）的计算结果和 薄船线性兴波阻力理

Taylor论的计算­结果相比更接近于 图谱计算值， Fr>0.35尤其是在 的高速阶段。同时，还发现在低速阶段存在­直接边界元法计算不准­确的现象。下文将对该部分的误差­进行详细分析。2 Wigley第 个算例为另一种 船型。该模型曾605 =5m，在 所进行过阻力试验，模型参数如下：L B=0.4 m，T=0.178 m。通过Prohaska­方法确定模型形状因子，进而可近似得到兴波阻­力的试验值。4 2 Wigley图 所示为第 个算例中 船型的兴波阻力试验结­果与计算结果的比较。 4由图 可以看出，基于复合算法的直接边­界Michell元法­和间接边界元法的计算­精度要高于薄船线性兴­波阻力理论的精度，尤其是在中、高速段。对比直接边界元法和间­接边界元法，发现直接边界元法的计­算精度比间接边界元法­的稍高。Fr<0.2但正如前文所述，本文所用算法在 的低速2段计算误差较­大。表 所示为各算法精度的定­量分析结果。

2.2 低速段误差分析

3结论

1］ . ［ 戴遗山，段文洋 船舶在波浪中运动的势­流理论［M］. 北京：国防工业出版社，2008. 2］ 徐刚. ［ 不规则波中浮体二阶水­动力时域数值模拟［D］. 哈尔滨：哈尔滨工程大学，2010. XU G. Time-domain simulation of second-order hydro⁃ dynamic force on floating bodies in irregular waves ［D］. Harbin：Harbin Engineerin­g University，2010 （in Chinese）. ［3］ MIAO Q M，CHWANG A T. Ship waves by direct boundary element method［J］. Journal of Ship Mechan⁃ ics，2003，7（6）：27-36. ［4］ 陈纪康，段文洋，朱鑫. 三维泰勒展开边界元方­法 及其数值验证［ J］. 水动力学研究与进展（ A 辑）， 2013，28（4）：482-485. CHEN J K，DUAN W Y，ZHU X. Three-dimension Taylor expansion boundary element method and it's val⁃ idation［J］. Chinese Journal of Hydrodynam­ics（Ser. A），2013，28（4）：482-485（in Chinese）. ［ 5］ 李井煜，卢晓平，赵鹏伟. 直接边界元法解势流速­度场问题［J］. 中国舰船研究，2015，10（1）：68-75.

LI J Y ， LU X P ，ZHAO P W. Direct boundary element method for the problem of potential flow velocity field ［J］. Chinese Journal of Ship Research， 2015， 10 （1）：68-75（in Chinese）. ［ 6］ 韩玉超，卢晓平，王中. 无限区域二维势流直接­边析［J］. 中国舰船研究，2015，10（4）：界元法精度分39-45，54. HAN Y C ， LU X P ，WANG Z. Precision analysis of the two-dimensiona­l potential flow problem in an infi⁃ nite region with the direct boundary element method ［J］. Chinese Journal of Ship Research， 2015， 10 （4）：39-45，54（in Chinese）. 7］刘应中. 船舶兴波阻力理论［M］. ［ 北京：国防工业出版社，2003. 8］ 王献孚，周树信，陈泽梁. 计算船舶流体力学［M］. ［ 上海：上海交通大学出版社，1992.

9］ 程天柱，石仲堃. ［ 兴波阻力理论及其在船­型设计中的应用［M］. 武汉：华中工学院出版社，1987. 10］ 熊鹰，卢晓平. M］. ［ 舰船快速性计算图册［ 武汉：

2 Wigle图 船型网格Fig.2 Panel mesh of the Wigley ship form