截断圆柱绕射解析解中水线上速度势及其z向偏导数的收敛性研究.pdfVIP

第二十一届全国水动力学研讨会暨第八届全国水动力学学术会议 暨两岸船舶与海洋工程水动力学研讨会文集 截断圆柱绕射解析解中水线上速度势及其z 向偏导数的收敛性研究 耿宝磊 勾莹 滕斌 （大连理工大学 海岸和近海工程国家重点实验室，大连 116024） 摘要： 通过分区离散和边界匹配的方法可以得到波浪关于截断圆柱绕射的解析解，计算结果 的精度取决于解析式中展开项的项数。本研究分别研究了水线上展开项数m、i、j 对φ ，∂φ /∂z ， 2 2 ∂ φ /∂z 的收敛性的影响。并在此基础上，给出了计算中的处理方法和建议取值。 关键词： 截断圆柱；解析解；速度势；z 向偏导数；收敛性 1 引言 相比数值方法而言，通过解析方法可以得到波浪与结构物作用的精确结果，自从Thorne[1] 开始将多极展开方法用于处理表面波浪问题以来，学者们已经得到了很多简单几何形状物体 [2] [3] [4] 的解析解。Black 等计算了二维方箱和直立圆柱的解析解。Havelock 和MacCamy and Fuchs 对无限水深和有限水深中的均匀直立圆柱做了研究，得到了波浪绕射问题的解析解。之后 [5] [6] Garret 和 Yeung 分别得到了三维截断圆柱绕射问题和辐射问题的解析解。解析解的计算精 度取决于展开项项数的选取，虽然只有一些简单的结构形式才存在解析解，但是作为检验数 值方法正确与否的手段，有必要对解析解的精度及其收敛性开展研究。研究了三维截断圆柱 绕射问题解析解的收敛性问题，文中第二部分介绍了通过分区离散和边界匹配建立的波浪关 于截断圆柱绕射的解析方法，第三部分分析水线处解析解的展开项数对速度势、速度势 z 方 向导数及其二阶导数计算结果的影响，第四部分总结得到的结论，给出取值建议。 2 波浪关于截断圆柱的绕射问题解析解的建立 考虑波浪对在水深d 中，半径为a ，吃水为T 的圆柱的绕射问题( 图1)。取柱坐标系orθz ， 使得orθ在静水面上，oz 轴通过柱中心，向上为正。设φ (r,θ,z ）为波浪的速度势，则 满足的 φ 基金项目：国家自然科学基金10772040）和国家高技术研究发展计划（863 计划）（2006AA09A109 －3 ）。作者简介：耿宝磊（1980－），男，河北衡水人，博士生，主要从事波浪对海上建筑物作用方面的研究。 - 93 - 第二十一届全国水动力学研讨会暨第八届全国水动力学学术会议 暨两岸船舶与海洋工程水动力学研讨会文集 边界条件包括： 2 ω φz (r,θ , z ) φ (r,θ , z ) ( z 0) （1） g φz (r,θ , z ) 0 ( z −d ) （2 ） φr (r,θ , z ) 0 (r a, − T z 0) （3 ） φz (r,θ , z ) 0 (r = a, z =−T ) （4 ） φ 分解为入射势φ 和绕射势φ 两部分。φ 向外传播 图1 波浪关于截断圆柱的绕射 i d d 并逐渐衰减，满足无穷远处Sommerfeld 散射条件： 1/ 2  ∂φd  lim r  − ikφd  0