混合式CRP推进器是一种组合传统轴系螺旋桨和吊舱推进器实现的新型对转桨推进方式,具有推进效率高、操纵性能好、激振力和废气排放低等优点,在绿色船舶的背景下具有广阔的应用前景。
2001年,韩国三星船模试验池[1]对超大型集装箱船的单桨、双尾鳍双桨和混合式CRP推进器3种推进系统进行了对比试验。结果表明,混合式CRP推进器总效率比双尾鳍双桨高9%,比常规单桨高5%。Sasaki N等[2]对采用混合式CRP推进器的船舶模型试验规程进行了探索性研究。他们认为采用该型推进器的船模试验需要涉及不带吊舱推进器的阻力试验、带混合式CRP推进器的自航试验、混合式CRP推进器敞水试验等。Sasaki N等[3]进一步研究了混合式CRP推进器的设计,并改进了敞水试验方法。认为混合式CRP敞水试验需包括前桨敞水试验,前桨带吊舱舱体的敞水试验、吊舱推进器敞水试验以及混合式CRP推进器敞水试验等。Black S等[4]在英国纽卡斯尔大学水洞进行了混合式CRP推进器的敞水试验。试验并未考虑前后桨之间的间隙,但在数值计算时计入了桨毂的影响。他们还进一步进行了空泡试验。Shimamoto H等[5]评估了混合式CRP推进器的整体性能。研究发现与传统的单桨推进相比,采用混合式CRP推进器的集装箱船舶不但可以节省功率而且可以提高船舶的操纵性能。Bong J C等[6]采用试验方法研究了采用混合式CRP推进器的集装箱船的实船推进性能。研究发现:在某航速下,前后桨的转速比与前后桨的功率比是唯一相关的;在设计工况下,后桨对前桨的影响是可以忽略不计的。Schez A等[7]介绍了欧盟TRIPOD计划,并从水动力节能考虑,指出混合式CRP推进器的前桨可以考虑采用CLT型桨,进一步介绍了如何将混合式CRP推进器水动力性能模型试验数据外推至实尺度推进器,并给出了如何采用CFD方法准确预报其自航性能。在国内,盛立[8] 、张可[9]率先采用RANS方法结合SST k-ω湍流模型研究了混合式CRP推进器的敞水性能,通过与试验数据的比较验证了数值计算方法的准确性。汪小翔[10]采用RANS方法计算了混合式CRP推进器的定常水动力性能和非定常水动力性能。第26届ITTC推进委员会[11]认为,关于混合式CRP推进器的试验还比较少,相关研究还处于发展阶段,建议下届ITTC推进委员会起草关于混合式CRP推进器的试验规程。
由于混合式CRP推进器的吊舱推进器同时肩负舵的作用,在船舶操纵时必须进行左右回转,因此前后桨必须保持一定的距离才能正常工作。另一方面,混合式CRP推进器采用相互独立的功率输入,前后桨的转速比可以任意调节。鉴于此,本文重点研究前后桨的间距和转速比对混合式CRP推进器水动力性能的影响。
1 数值计算模型本文采用RANS方法结合SST 湍流模型和滑移网格技术计算混合式CRP推进器的敞水性能,该方法在混合式CRP推进器水动力性能预报方面具有较高的精度, 具体可参考文献[8,9]。以海军工程大学设计的4000TEU集装箱船混合式CRP推进器为研究对象,吊舱和前后桨的主要参数如表1和表2所示,混合式CRP推进器的外形轮廓如图1所示。
计算域的划分综合考虑对转桨和吊舱推进器的特点。流域为长方体,入口离前桨盘面中心3Lpod,出口离后桨盘面中心10Lpod,左、右、上、下远场离盘面中心各5Lpod(Lpod为舱体的长度)。整个流域分为4个部分,前方的来流域、前桨旋转域、后桨旋转域和后方的吊舱域。前、后桨旋转域采用圆柱体形状,划分周期性结构化网格,其直径为1.2DF(DF为前桨的直径)。来流域与吊舱域均为静止域,采用纯六面体单元离散,4个流域通过交界面传递通量。螺旋桨和吊舱壁面生成边界层网格,第一层网格的y+为80左右,边界层的增长率为1.05。整个计算域共400万左右的单元。考虑真实的桨毂形状,前后桨的桨叶、吊舱的表面网格分布如图2所示。
名称 | 前桨 | 后桨 |
直径/mm | 240 | 203.636 |
叶数 | 4 | 5 |
(P/D)0.7R | 1.162 2 | 1.302 7 |
毂径比 | 0.227 | 0.216 |
剖面类型 | NACA66mod/a=0.8 | NACA66mod/a=0.8 |
旋向 | 左旋 | 右旋 |
入口设为速度入口,给定均匀来流的速度值;出口设为压力出口;外域边界设为对称面;前桨子域按照滑移网格的要求设为绕x轴以角速度20 r/s旋转;后桨子域速度按照前后桨的转速比进行设置,方向与前桨子域相反;桨叶和桨毂相对子域的旋转速度为0,定义无滑移、不可穿透的壁面边界条件,吊舱亦定义为无滑移、不可穿透的壁面边界条件。进速系数的改变通过来流速度的改变实现。采用有限体积法离散控制方程和湍流模型,对流项和扩散项均采用二阶迎风格式进行离散,压力速度耦合迭代采用SIMPLEC方法。先采用定常方法计算,待收敛后再开启非定常模式,计算时间步为前桨旋转1°所对应的时间。
2 结果与分析
采用前、后桨各自的直径和转速对前、后桨的推力和扭矩系数进行无量纲化。进速系数J、前桨的推力系数、前桨的扭矩系数、后桨的推力系数、后桨的扭矩系数、吊舱单元的推力系数、混合式CRP推进器的推力系数、混合式CRP推进器的扭矩系数、混合式CRP推进器敞水效率 的定义如下:
\[\left\{ \begin{align}
& J=\frac{{{V}_{A}}}{{{\eta }_{F}}{{D}_{F}}},{{K}_{TF}}=\frac{{{T}_{F}}}{pn_{n}^{2}D_{F}^{A}} \\
& {{K}_{QF}}=\frac{QF}{pn_{F}^{2}D_{F}^{5}},{{K}_{TA}}=\frac{{{T}_{A}}}{pn_{A}^{2}D_{A}^{4}} \\
& {{K}_{QA}}=\frac{QA}{pn_{F}^{2}D_{F}^{5}},{{K}_{TA}}=\frac{{{T}_{A}}+{{R}_{pod}}}{pn_{A}^{2}D_{A}^{4}} \\
& {{K}_{T}}=\frac{{{T}_{F}}+{{T}_{A}}+{{R}_{pod}}}{pn_{F}^{2}D_{F}^{4}},{{K}_{Q}}=\frac{{{n}_{F}}{{Q}_{F}}+{{n}_{A}}{{Q}_{A}}}{pn_{F}^{3}D_{F}^{5}} \\
& {{\eta }_{0}}=\frac{\left( {{T}_{F}}+{{T}_{u}} \right){{V}_{A}}}{2\pi \left( {{n}_{F}}{{Q}_{F}}+{{n}_{A}}{{Q}_{A}} \right)} \\
\end{align} \right.\]
(1)
为了研究前桨、后桨和吊舱之间的相互干扰影响,同时计算了单独前桨和单独吊舱推进器的敞水性能。螺旋桨的转速和来流速度与混合式CRP推进器完全一致。单独前桨和混合式CRP推进器的前桨敞水性能的对比如图3所示,单独吊舱推进器和混合式CRP推进器的吊舱推进器敞水性能的对比如图4所示。
从图3可以看出:单独前桨与混合式CRP推进器的前桨敞水性能基本一致,在工作点附近,推力系数的差别不超过3%,而扭矩系数的差别不超过2%。其原因是前桨处于吊舱推进器的上游,后桨对前桨的抽吸作用和吊舱对流体的阻塞作用相互抵消。
从图4可以看出:前桨诱导速度对吊舱推进器的水动力性能产生了重要的影响。由于前桨尾流的加速作用,混合式CRP推进器的后桨和吊舱单元的推力系数和扭矩系数比单独吊舱推进器小很多,这种差别随着进速系数的减小而增大。
J=0.781时,单独前桨、单独吊舱推进器和混合式CRP推进器在y=0剖面处无量纲轴向速度分布的对比如图5所示。
从图5可以看出:单独前桨和混合式CRP推进器的前桨前方来流速度分布大小及形态基本一致;而单独吊舱推进器与混合式CRP推进器的后桨前方来流速度分布差异较大。相对于单独的吊舱推进器,混合式CRP推进器的后桨前方来流极其不均匀,速度值较大。
单独前桨单独吊舱推进器和混合式CRP推进器桨后方的流线形态如图6所示。
从图6可以看出:单独前桨时,后方的流体质点具有很强的旋转能量,桨盘面正后方很长距离内,流线呈螺旋状。对于单独吊舱推进器,吊舱的整流作用使桨后方螺旋状流线分布面积减小;进一步观察发现,流线经过吊舱包时,产生一定的爬升。而对于混合式CRP推进器,经过后桨对转吸收前桨尾流旋转能量以及吊舱的整流作用,桨后方基本上不含螺旋状流线。
2.2 前后桨间距的影响设计了4种前后桨盘面中心间距方案,分别为0.394DF、0.455DF、0.53DF和0.61DF。J=0.781时,前桨、吊舱单元和混合式CRP推进器的敞水性能随间距的变化如图7所示。
从图中可以看出:
1)前桨的推力系数和扭矩系数受间距的影响不大。主要是由于前桨处于后桨的上游,一定范围内,后桨的抽吸作用使前桨来流速度的增加变化不大;
2)后桨和吊舱单元的推力系数随着间距的减小而增大,后桨的扭矩系数亦随间距的减小而增大,增大的幅度随着间距的增加有变缓的趋势,同时低进速系数的增大幅度比高进速系数大。这主要是由于间距减小,后桨回收前桨尾流能量增强引起的;
3)混合式CRP推进器总的推力系数和扭矩系数随着间距的减小而略有增加,敞水效率随着间距的减小而增大。从节能的角度上看,应当在保证吊舱推进器正常回转的前提下尽量减小前后桨的间距。
2.3 前后桨转速比的影响设计了5种转速比方案,分别为nA/nF=1.05、1.075、1.104、1.125和1.15。J=0.781时,前桨、吊舱单元和混合式CRP推进器的敞水性能随转速比的变化如图8所示。后桨占整个混合式CRP推进器的功率百分比随转速比的变化如图9所示。
后桨占整个混合式CRP推进器的功率百分比定义如下:
\[\frac{{{P}_{A}}}{\left( {{P}_{F}}+{{P}_{A}} \right)}=\frac{2\pi {{n}_{A}}{{Q}_{A}}}{\left( 2\pi {{n}_{F}}{{Q}_{F}}+2\pi {{n}_{A}}{{Q}_{A}} \right)}\]
(2)
从图9中可以看出:
1)前桨的推力系数和扭矩系数随着转速比的增大而略有减小,其原因是后桨转速增加,对前方流体的抽吸作用增强,进而导致前桨来流速度增大而引起的,但这种影响极其有限;
2)后桨的推力系数、扭矩系数以及吊舱单元的推力系数随转速比的增大而显著增大;一定进速系数下,吊舱单元的敞水性能与转速比呈线性关系;
3)混合式CRP推进器的推力系数和扭矩系数随转速比的增大而增大,但扭矩系数增大的幅值比推力系数大,因而其敞水效率随转速比的增大而减小。一定进速系数下,混合式CRP推进器的敞水性能与转速比呈线性关系。
4)后桨占整个混合式CRP推进器的功率百分比随转速比和进速系数的增大而增大。
3 结论
采用RANS方法结合SST 湍流模型计算了混合式CRP推进器的敞水性能,详细分析了前后桨与吊舱之间的相互干扰作用,研究了主要设计参数对混合式CRP推进器水动力性能的影响,为混合式CRP推进器的尺度效应研究及非设计工况下水动力性能预报奠定了基础。结果表明:
1)处于后方的吊舱推进器对前桨的干扰较小,设计工况附近可以忽略不计,而前桨尾流对吊舱推进器水动力性能产生较大的影响;
2)混合式CRP推进器总的推力系数随着间距的增大而减小,而扭矩系数则变化不大,敞水效率随着间距的增大而减小,故从节能的角度上看,应当在保证吊舱推进器正常回转的前提下尽量减小前后桨的间距;
3)混合式CRP推进器的推力系数和扭矩系数随转速比的增大而增大,但敞水效率随转速比的增大而减小;相同进速系数下,混合式CRP推进器的敞水性能与转速比呈线性关系;后桨占整个混合式CRP推进器的功率百分比随转速比和进速系数的增大而增大。
[1] | KIM S E, CHOI S H. Model tests on propulsion systems for ultra large container vessel[C]//Proceedings of the 12th International Offshore and Polar Engineering Conference. Osaka, 2002:520-524. |
[2] | SASAKI N, KAWANAMI Y, UKON Y, et al. Model test procedure and analysis of hybrid CRP POD system[C]//Proceedings of the 2nd International Conference on Technical Advances in Podded Propulsion(T-POD). Nantes, France, 2006. |
[3] | SASAKI N, KURODA M, FUJISAWA J, et al. On the model tests and design method of hybrid CRP podded propulsion system of a feeder container ship[C]//Proceedings of the 1st International Symposium on Marine Propulsors smp'09. Trondheim, Norway,2009. |
[4] | BLACK S D, CUSANELLI D S. Design and testing of a hybrid shaft-pod propulsor for a high speed sealift ship[C]//Proceedings of Propellers/Shafting 2003 Symposium. Virginia, 2009. |
[5] | SHIMAMOTO H, TAKEDA A, MIYAKE S. Tandem hybrid CRP (Contra-Rotating Propeller) System[J]. Journal of the Marine Engineering Society in Japan, 2011, 46(3):326-329. |
[6] | CHANG B J, GO S. Study on a procedure for propulsive performance prediction for CRP-POD systems[J]. Journal of Marine Science and Technology, 2011, 16(1):1-7. |
[7] | SÁNCHEZ-CAJA A, PÉREZ-SOBRINO M, QUEREDA R, et al. Combination of Pod, CLT and CRP propulsion for improving ship efficiency:the TRIPOD project[C]//Proceedings of the 3rd International Symposium on Marine Propulsors smp'13. Tasmania, Australia, 2013. |
[8] | 盛立, 熊鹰. 混合式CRP吊舱推进器水动力性能数值模拟及试验[J]. 南京航空航天大学学报, 2012, 44(2):184-190. SHENG Li, XIONG Ying. Numerical simulation and experimental investigation on hydrodynamics performance of hybrid CRP podded propulsion[J]. Journal of Nanjing University of Aeronautics & Astronautics, 2012, 44(2):184-190. |
[9] | 张可. 混合式CRP推进器设计与水动力性能研究[D]. 武汉:海军工程大学, 2013:60-90. ZHANG Ke. Research on hybrid CRP propulsor design and Hydrodynamic Performance[D]. Wuhan:Naval University of Engineering, 2013:60-90. |
[10] | 汪小翔. 混合式CRP吊舱推进器的水动力性能研究[D]. 哈尔滨:哈尔滨工程大学, 2013:47-70.WANG Xiaoxiang. Research on the hydrodynamic performance of hybrid CRP podded propulsion[D]. Harbin:Harbin Engineering University, 2013:47-70. |
[11] | ANZBÖCK R, LEER-ANDERSEN M, ATLAR M, et al. Final report and recommendations to the 26th ITTC[C]//Proceedings of the 26th ITTC. Rio de Janeiro, Brazil, 2011. |