舰船科学技术  2021, Vol. 43 Issue (10): 25-28    DOI: 10.3404/j.issn.1672-7649.2021.10.006   PDF    
基于CFD的带呆木船舶尾部流场分析
饶志强1,2, 魏波1, 樊涛1, 胡天翔1, 朱伯华1, 吕江3, 郭兴乾4     
1. 江南造船(集团)有限责任公司,上海 200219;
2. 上海交通大学 机械与动力工程学院,上海 200240;
3. 中国船舶科学研究中心,江苏 无锡214082;
4. 中国船舶与海洋工程设计研究院,上海 200011
摘要: 为了研究船舶呆木附近流场特性,本文采用计算流体力学方法研究了一艘带呆木的船舶尾部流场。数值仿真首先分析了直航状态下的尾部呆木压力分布和呆木泄出涡形态,然后分析了回转状态下的泄出涡特征。研究表明呆木末端是全船压力最低点,在高速直航时有明显的泄出涡产生,而回转状态下泄出涡大幅减弱或消失。船舶直航时呆木泄出涡存在空化的可能,在船舶设计中不能忽略。
关键词: 船舶     呆木     尾流     泄出涡     CFD    
CFD based numerical analysis of the wake flow of a ship fitted with a deadwood
RAO Zhi-qiang1,2, WEI Bo1, FAN Tao1, HU Tian-xiang1, ZHU Bo-hua1, LV Jiang3, GUO Xing-qian4     
1. Jiangnan Shipyard (Group) Co. Ltd., Shanghai 200219, China;
2. School of Mechanical Engineering, Shanghai Jiaotong University, Shanghai 200240, China;
3. China Ship Scientific Research Center, Wuxi 214082, China;
4. Marine Design and Research Institute of China, Shanghai 200011, China
Abstract: In this paper CFD based numerical simulations are carried out to study the wake flow characteristics of the deadwood fitted on a ship. Firstly, pressure distribution and shedding vortex shape around deadwood under the condition of straight sailing are simulated. Then the vortex characteristics of the shedding vortex under the turning around condition are compared. The study shows that the deadwood is the region of lowest pressure on the hull surface. And obvious shedding vortex occurs when the ship has a working condition of straight sailing. Under the turning around condition, shedding vortex is much weaker than those of straight sailing. The vortex shedding from deadwood has possibilities of cavitation, which should not be ignored in ship design.
Key words: ship     deadwood     wake flow     shedding vortex     CFD    
0 引 言

随着造船技术和航运业的发展,船舶正向大型化和高航速方向发展。在这个过程中,船舶尾部附体数量和尺寸都会增多或增大。尾部区域附体最多,也最容易发生空化,而空化引起的剥蚀、振动、噪声等问题一直是人们关注的重点。通常认为,船尾容易发生空泡的部件主要是轴支架、螺旋桨和舵叶等附体。对这些附体的空泡现象已经有一些成果供参考。对于轴支架空泡现象,程道明等[1]在减压水池中对某型船开展了试验并观测轴支架空泡。黄洪波等[2]对一艘双桨船的轴支架进行了多安装角度试验分析,得到了不同支架角度对应的空泡形态以及螺旋桨的水动力性能供船舶设计参考。吕江等[3]采用Eppler方法对轴支架剖面进行了设计并开展模型试验,试验证明Eppler设计支架能够提高轴支架的空泡裕度。李亮等[4]采用内窥镜观测了实船支架上的空泡,为实船支架臂空泡形态提供了直观认识。而对于舵叶表面空化问题,于安斌等[5]对扭曲舵的空化特性进行了模型试验,并发现扭曲设计能够有效提高舵叶空泡初生舵角。叶金铭等[6]采用CFD方法对舵叶空化开展了数值研究和实船观测,结果表明舵表面会先于螺旋桨表面发生空化。曹彦涛等[7]同样采用内窥镜对实船舵叶空泡进行了观测,结果表明舵叶表面的牺牲阳极最容易发生空化,且悬空的舵叶下端有涡空泡脱落,舵叶中部在高速时有大面积空泡出现。相对于轴支架和舵叶,螺旋桨空化的研究则丰富得多。熊鹰等[8]、张永坤等[9]、季盛等[10]、刘登成等[11],武坤等[12]采用数值模拟或试验方法开展了螺旋桨空化特性的研究。

近年来,很多大型船舶都在船舶尾部安装了尺寸较大的呆木来提高尾部结构安全性或匹配船坞内建造工艺。一般认为,呆木对船舶性能影响最大的是操纵性。赵占华等[13]采用模型试验的方法证实呆木能提高一艘工程船的拖航稳定性。曹留帅等[14]采用分离型操纵性数学模型研究了呆木对操纵性的影响,结果表明呆木面积与航向稳定性提高成正比。而对于尾部呆木周围的流场及其流动特性研究,则尚未发现相关文献发表。

本文采用计算流体力学方法研究一艘带呆木的船舶尾部流场。数值仿真首先分析了直航状态下的尾部呆木压力分布和呆木泄出涡形态,然后分析了回转状态下的泄出涡特征。研究表明呆木末端是全船压力最低点,在高速直航时有明显的泄出涡产生,而回转状态下泄出涡大幅减弱或消失。船舶直航时呆木泄出涡存在空化的可能,在船舶设计中不能忽略。

1 流场数值分析 1.1 物理模型和数值方法

为了开展流场模拟,在模型尺度船体周围创建计算域。其中,直航工况的计算域为方形,回转工况的计算域为环形。在数值模拟中,采用RANS方法对流场进行计算,湍流模型采用SST k-ω湍流模型。计算域的入口采用速度入口,出口采用压力出口,侧面远场采用对称边界(直航工况)或者滑移壁面边界条件保证流体在远场边界上不穿透(回转工况),消除远场边界条件对船体周围流场的影响。图1图2为CFD数值仿真的计算域及具体边界条件设置,右下角为尾部呆木表面网格。直航工况网格总数约为640万,回转工况网格总数约为1300万,船体和附体表面均划分边界层网格。为了分析流场特性,采用同一个航速进行数值计算。

图 1 直航工况计算域和边界条件设置 Fig. 1 Domain and boundary conditions of straight sailing condition

图 2 回转工况计算域和边界条件设置 Fig. 2 Domain and boundary conditions of turning condition
1.2 直航数值分析

首先对直航工况进行数值分析。在船舶流场中,低压区和涡是形成空化的主要因素。因此可以通过压力分布和涡系特性来定性判断尾流场的基本情况。图3为直航工况的尾部压力分布。从压力分布的结果来看,船尾呆木末端由于形状突扩是一个压力极低的区域,为全船压力最低点(图3虚线圈内)。图4图5为尾部流场的Q准则等值涡面(Q=10000,下同)以及涡面上的压力分布。其中图5还给出了呆木后方的流线。从涡面结果来看,尾部呆木末端形成了2个明显的泄出涡向船尾运动。从流线可以看出,向船尾运动的泄出涡呈明显的螺旋线形式,是名副其实的“涡”。

图 3 直航工况尾部压力分布 Fig. 3 Pressure distribution on stern under straight sailing condition

图 4 直航工况等涡面及其表面压力分布 Fig. 4 ISO-surface of vortex and its pressure distribution under straight sailing condition

图 5 直航工况等涡面及附近流线 Fig. 5 ISO-surface of vortex and streamlines

需要注意的是:2股泄出涡在向后运动过程中一直与船体平板龙骨保持一定距离。呆木末端的低压区和高强度泄出涡使呆木相应位置有发生空化的可能,需要特别注意。

1.3 回转数值分析

对回转工况下的船尾流场进行数值计算和分析。图6为中等舵角回转工况下的尾部压力分布。在回转工况中,呆木末端压力上升明显,这意味着该处发生空化的可能也在降低。图7为尾部流场的Q准则等值涡面及涡面上的压力分布。回转工况下,尾部呆木末端仍然会形成2股泄出涡,但泄出涡体积和强度大幅度减弱。呆木末端压力有所升高,泄出涡的强度减弱,这意味着回转工况下呆木发生空化的概率会降低。

图 6 回转工况尾部压力分布 Fig. 6 Pressure distribution on stern under turning condition

图 7 回转工况等涡面及其表面压力分布 Fig. 7 ISO-surface of vortex and its pressure distribution under turning condition
2 数值分析结果

从数值分析结果来看,直航状态时尾部呆木末端存在明显的低压区和泄出涡。而回转工况下呆木末端的低压区和泄出涡上的压力均有所升高,回转工况发生尾部流场空化的可能性也会有所降低。究其原因是因为直航工况下呆木末端附近的流动以从下到上的垂向流动为主,而回转工况中呆木附近的流动以从左到右(右舵回转时)的侧向流动为主。垂向流动由于船型特征存在突然扩张,这种突扩会引起压力急剧降低[15]

船上发生的空泡很大一部分都是由涡流产生,比如螺旋桨梢涡或毂涡空泡,悬挂舵叶下端的空泡等等。对本船而言,呆木末端为全船压力最低点,在该点又形成了高强度的涡流,这种低压和高强度泄出涡的叠加很可能会形成空泡。这个潜在的问题在船体线型设计时尤其需要注意。

3 结 语

本文采用CFD数值模拟方法对一艘带呆木的船舶尾部流场进行了研究,得出以下结论:

1)船舶直航时,呆木末端是全船压力最低点且呆木末端有明显的泄出涡产生,呆木末端处存在空化风险,应予以特别关注。

2)船舶回转时,呆木末端压力有所升高,且没有明显的泄出涡出现。

3)CFD数值技术有助于发现船舶设计初期的的潜在问题。

船舶尾流场一直是船舶外流场最复杂的区域,其对尾部各附体都有影响。本文开展的尾部呆木流场分析对后续的船舶设计和建造具有一定的参考意义。

参考文献
[1]
程明道, 王莉萍, 朱建良. 高速排水型船支架空化观测[J]. 船舶力学, 2000, 4(2): 19-26.
CHENG M D, WANG L P, ZHU J L. Observation of cavitation for a high-speed displacement ship[J]. Journal of Ship Mechanics, 2000, 4(2): 19-26.
[2]
黄红波, 许晖, 王建芳, 等. 多桨船双臂支架空泡性能优化及其对螺旋桨空泡性能影响研究[J]. 中国造船, 2015, 56(2): 150-158.
HUANG H, XU H, WANG J, et al. Cavitation optimization of v-bracket and its influence on cavitation performance of propellers for multi-propulsor ship[J]. Shipbuilding of China, 2015, 56(2): 150-158. DOI:10.3969/j.issn.1000-4882.2015.02.018
[3]
吕江, 王东涛, 曾志波, 等. 轴系支架空化性能优化研究[C]//第十四届全国水动力学术会议暨第二十八届全国水动力学研讨会, 中国, 吉林长春, 2017.
LV J, WANG D T, ZENG Z B, et al. Optimization research on cavitation performance of the shaft-strut[C]//The 14th National Congress on Hydrodynamics & 28th National Conference on Hydrodynamics, Changchun, Jilin, China, 2017.
[4]
李亮, 张伟, 陆芳, 等. 实船轴支架空泡观测技术研究[C]//第十四届全国水动力学术会议暨第二十八届全国水动力学研讨会, 中国, 吉林长春, 2017.
LI L, ZHANG W, LU F, et al. Research on cavitation observation technique for full Scale ship shaft bracket[C]//The 14th National Congress on Hydrodynamics & 28th National Conference on Hydrodynamics, Changchun, Jilin, China, 2017.
[5]
于安斌, 叶金铭, 王友乾. 扭曲舵抗空化性能模型试验研究[J]. 推进技术, 2019, 40(1): 215-222.
YU A B, YE J M, WANG Y Q. Model test research on anti-cavitation performance of twisted rudder model[J]. Journal of propulsion technology, 2019, 40(1): 215-222.
[6]
叶金铭, 于安斌, 王友乾, 等. 桨后舵片空泡的分离涡模型模拟及实船舵空泡试验[J]. 哈尔滨工程大学学报, 2019, 40(5): 913-919.
YE J, YU A, WANG Y, et al. Detached eddy simulation of sheet cavitation of rudder behind propeller and real ship rudder cavitation test[J]. Journal of Harbin Engineering University, 2019, 40(5): 913-919.
[7]
曹彦涛, 叶金铭, 陆芳, 等. 实船舵空化现象试验观察[C]//纪念《船舶力学》创刊二十周年学术会议, 中国, 浙江舟山, 2017.
CAO Y, YE J, LU F, et al. Observation of Cavitation on a Full Scale Rudder[C]// Commemoration Conference for the 20th Anniversary of Journal of Ship Mechanics, Zhoushan, Zhejiang, China, 2017.
[8]
熊鹰, 叶金铭, 樊晓斌, 等. 螺旋桨空泡及其诱导的脉动压力的数值和试验研究[J]. 船舶力学, 2006, 10(5): 32-40.
XIONG Y, YE J M, FAN X B, et al. Numerical and experimental research on pressure fluctuations induced by propeller cavitation[J]. Journal of Ship Mechanics, 2006, 10(5): 32-40. DOI:10.3969/j.issn.1007-7294.2006.05.005
[9]
张永坤, 熊鹰, 叶金铭. 水中含气量对螺旋桨空泡噪声影响的试验研究[J]. 武汉理工大学学报(交通科学与工程版), 2009, 33(2): 234-237.
Zhang Y, Xiong Y, Ye J. Experimental investigations of influence of gas content in water to propeller cavitation noise[J]. Journal of Wuhan University of Technology (Transportation Science & Engineering), 2009, 33(2): 234-237.
[10]
季盛, 董国祥, 陈建挺, 等. 带节能附体螺旋桨空泡脉动压力的试验研究[J]. 上海船舶运输科学研究所学报, 2012, 35(2): 39-46.
JI S, DONG G, CHEN J, et al. Experimental study of propeller induced cavitation and pressure pulse with energy-saving appendages[J]. Journal of Shanghai Ship and Shipping Research Institute, 2012, 35(2): 39-46. DOI:10.3969/j.issn.1674-5949.2012.02.007
[11]
LIU D C, ZHOU W X. Numerical predictions of the propeller cavitation pressure fluctuation behind ship and comparison with experiment[J]. Journal of Ship Mechanics, 2019, 23(3): 245-254.
[12]
武珅, 芮伟, 曾志波. 冰块阻塞状态的螺旋桨空泡激振力特性分析[J]. 舰船科学技术, 2019, 41(19): 21-26.
WU S, RUI W, ZENG Z B. Investigation on the characteristics of cavitation excited force of propeller in ice block condition[J]. Ship Science and Technology, 2019, 41(19): 21-26.
[13]
赵战华, 匡晓峰, 范亚丽, 等. 呆木对拖航稳定性影响的试验研究[C]//第十三届全国水动力学学术会议暨第二十六届全国水动力学研讨会, 山东青岛, 2015.
ZHAO Z H, KUANG X F, FAN Y L, et al. Experimental study of effect of dead wood on towing stability[C]//The 13th National Congress on Hydrodynamics & 26th National Conference on Hydrodynamics, Qingdao, Shandong, China, 2017.
[14]
曹留帅, 朱军, 黄昆仑. 舰船呆木设计及对操纵性影响研究[J]. 中国舰船研究, 2012, 7(2): 20-23+28.
CAO L S, ZHU J, HUANG K L. Design of the dead wood and its effect on ship maneuverability[J]. Chinese Journal of Ship Research, 2012, 7(2): 20-23+28. DOI:10.3969/j.issn.1673-3185.2012.02.004
[15]
L.普朗特, K.奥斯瓦提奇, K.维格哈特著, 郭永怀, 陆士嘉译. 流体力学概论[M]. 北京: 科学出版社, 1984: 276−277.