双体船粘性阻力数值计算方法.pdfVIP

武汉造船 ’ 年第 ( 期（总第 !(’ 期） 文章编号： （ ） !!#!$% ’ (#!#) 双体船粘性阻力数值计算方法 严 刚 刘启国 马廷淮 （华中科技大学 船舶与海洋工程系，武汉 (+ ） 摘 要：将兰德韦伯回转体速度分布所满足的第 类弗雷德和姆积分方程由来流速度为 推广到一任意 ! ! 值，并由此提出了速度因子的概念；在数值计算过程中，结合实际情况，考虑到多次迭代会带来舍入误差的积 累，采用最小偏差法结束迭代；程序用 , - - 语言编制，在 ./0123 , - - .45) 6 环境下调试成功。 关键词：双体船；粘性阻力；速度因子；等价回转体；数值计算 中图分类号： ； - 文献标识码： 1$+ 6*)! 1$$!6 (! ! 2 ： ， !#$%’$ /7 89:; 54=58 894 ;44 1 :7 ?549=@A 4BCD8:=7 :; E4745D@:F4 G5=A =74 8= D 5D7=A 7CAH45 D7 8947 D I=7# I48:=7 7DA4 ;44#GDI8=5 :; 5D:;4 6 /7 I=7;:45D8:=7 =G DIICAC@D8:=7 =G C7DJ=:DH@4 455=5 C5:7E 894 45:= =G 7CA45:ID@ - - ， ID@IC@D8:=7 894 A489= 7DA4 A:7:AD@ 455=5 :; D@:4 6 K94 5=E5DA :; I=A:@4 D7 4HCEE4 :7 ./0123 , ) 6 ;CI# I4;;GC@@L 6 ： ； ； ； ； ()* +,%-# ,D8DAD5D7 .:;I=C; 54;:;8D7I4 044#GDI8=5 MBC:JD@478 EL5D8=5 NCA45:ID@ ID@IC@D8:=7 直线运动的物体不受阻力。这与我们通常所见的 ! 等价回转体理论 实际情况不符，其原因就是由于粘性作用，在物体 为了计算船舶在运动中的粘性阻力，必须先 表面产生边界层。 求得船体上的势流压力分布，继而计算表面的三 界层内流体特性依然可用 （ ，）描述，由经 % 元层流和紊流边界层特性。但实际上，由于船体 典边界层理论，有如下结论： 复杂的几何形状，加上三元紊流边界层计算方法 ）界层内速度 在 方向分布是随表面压强 ! % ’ 的缺乏，使得该问题的精确求解非常困难。较为 在 方向梯度 变化的，当 时， O ( ( ! 简单地，可采用等价回转体法，即将船体水下部分