当今世界,计算机硬件水平显著提高,数值计算方法飞速发展,船舶CFD凭借其独特的优势越来越引起造船界的关注。通过数值计算研究船桨干扰的方法主要有2种:一是用体积力来模拟代替螺旋桨真实旋转的力场模拟方法[1–7];二是将船桨整体建模计算、考虑螺旋桨真实形状和旋转的整体计算方法。Visonneau M等[8]考虑全附体的影响,研究了模型尺度的不同对于船桨干扰产生的影响。ZHANG Nan等[9]计入自由液面影响,采用滑移网格模拟螺旋桨旋转,计算研究了潜艇艇体与螺旋桨相互干扰对其水动力性能及自航因子的影响。Choi等[10]计算研究了船桨舵之间的相互影响,所得结果与试验结果吻合较好。Han K J等[11]研究了桨舵性能受到舵的位置变化的影响,通过船桨舵干扰的计算结果对舵的布置位置提出了一些参考建议。Roberto M等[12]运用重叠动网格技术研究了船桨舵的非定常干扰问题,所得结果相比定常计算结果更为准确,与实验值能够很好地吻合。
我国对于船桨干扰问题的研究也很多。张志荣等[13]以混合面方法实现了船桨的整体计算,通过与试验结果的对比,表明该方法可以有效模拟船桨之间的相互影响。沈兴荣等[14]研究了粘性流场中桨舵的干扰,并尝试预报了舵的空泡性能,与试验观察结果吻合较好。王金宝等[15]考虑了自由液面效应和非定常流场下的螺旋桨旋转效应,对KCS整体模型进行了计算研究。沈海龙等[16]用DES(Detached eddy simulation)湍流模型,基于滑移网格技术,对孤立螺旋桨、孤立船体和船桨干扰的非定常流场分别进行了计算,所得计算结果与实验结果一致。沈海龙等[17]分别在均匀伴流场和非均匀伴流场中对不同螺旋桨水动力性能进行了预报,体现出非定常方法在研究船桨干扰问题上的优越性。傅慧萍[18]分别运用混合面方法、多参考系方法和滑移网格方法进行了船桨整体计算,对这3种方法的优缺点进行了系统的比较。杨春蕾等[19]通过数值计算表明MRF法、动量源法和滑移网格法这3种方法都可以运用到研究船桨干扰问题的计算中。
以上研究的都是单桨船的尾部干扰问题,多桨船的船桨干扰问题并未得到深入的研究。覃新川等[20]、王展智等[21]分别通过面元法和数值计算方法分析研究了四桨船舶桨的布局对螺旋桨水动力性能的影响。王展智等[22]建立了船桨舵整体计算模型,基于滑移网格探讨了舵的位置对桨性能的影响。然而上述研究并未考虑附体的影响。
由于多桨船尾部桨、舵、附体较多,尾部流场更为复杂,这对于螺旋桨之间的相互干扰产生一定影响。因此,在船体尾部附体的影响下,对多桨船的螺旋桨相互干扰问题进行研究,合理布置船后螺旋桨的相互位置,对于提高螺旋桨效率、改善螺旋桨盘面处伴流场起到一定的作用,能够为多桨船的设计提供指导和建议。本文以某四桨船舶为研究对象,采用整体建模方法,对船体、附体、螺旋桨进行整体建模并计算,探讨四桨布局对螺旋桨水动力性能的影响。
1 数学模型目前,对于N-S方程的统计平均方法以及补充反映湍流特性的其他方程,如湍动能方程和湍流耗散率方程,是目前的基本方法。
连续性方程是质量守恒定律在流体运动中的具体表现形式,流体连续性方程为:
$ \frac{{\partial \rho }}{{\partial t}} + \frac{\partial }{{\partial {x_i}}}(\rho {u_i}) = 0{\text{,}} $ |
雷诺平均N-S方程为:
$ \begin{array}{l} \displaystyle\frac{\partial }{{\partial t}}(\rho {u_i}) + \displaystyle\frac{\partial }{{\partial {x_j}}}(\rho {u_i}{u_j}) = - \displaystyle\frac{{\partial p}}{{\partial x{}_i}} + \\ \displaystyle\frac{\partial }{{\partial {x_j}}}\left[ \begin{array}{l} {\mu _0}(\displaystyle\frac{{\partial {u_i}}}{{\partial {x_j}}} + \displaystyle\frac{{\partial {u_j}}}{{\partial {x_i}}}) - \\ \displaystyle\frac{2}{3}{\mu _0}\displaystyle\frac{{\partial {u_l}}}{{\partial {x_l}}}{\delta _{ij}} \end{array} \right] + \displaystyle\frac{\partial }{{\partial {x_j}}}( - \rho \overline {u_i^{'}u_j^{'}} ) + \rho {f_i}{\text{,}} \end{array} $ |
式中的变量都是时历平均值,其中:
为了求解N-S方程中增加的雷诺应力,需要额外增加方程,以使雷诺应力中的脉动值与时均值联系起来,使方程组封闭,即湍流模型。本文采用的湍流模型为SST
本文所选取的船舶模型长度为8.56 m,采用四桨推进,左右舷各一对,左、右两侧的螺旋桨对称布置,每一侧的两桨前后布置,分别称为外前桨和内后桨,四桨的形状完全相同,均为5叶侧斜桨,直径为0.173 m,螺旋桨模型如图1所示。附体包括轴支架、轴包套、呆木和舵。螺旋桨、附体布置如图2所示,图3为轴支架和轴包套的局部放大图。
由于研究对象和计算流域的对称性,在建模和计算时均只考虑船舶右舷一半。计算坐标系x轴正向指向船尾,z轴正向垂直向上,y轴正向指向船模右舷,原点位于船体中纵剖线与船底水线的交点,如图4所示。计算模型的来流速度为2.572 m/s,螺旋桨转速为1 200 r/min。
图5为该模型前后桨的相对位置示意图。A,B分别表示内后桨和外前桨在X,Y方向上的相对距离,即纵向间距和横向间距。原始的布局方案为:A=2.545D,B=1.264D,其中D表示螺旋桨直径。
按照如下2种方式改变前后桨的相对位置:
1)改变尺寸A的大小,如表1所示。
2)改变尺寸B的大小,如表2所示。
计算域入口离船首1倍船长,出口离船尾2倍船长,外侧面离船中1倍船长,如图6所示。根据MRF模型的原理,将螺旋桨旋转子域单独划分出来,采用四面体的非结构化网格进行划分,每个旋转子域的网格约70万。由于轴包套、轴支架和螺旋桨的相对位置固定,因此螺旋桨位置的改变必然引起轴支架和轴包套位置的改变,考虑到各方案中螺旋桨位置的不断变化,将尾部附体与船体包在一起,尾部包体范围如图7所示。在此包体内采用四面体的非结构化网格,尾部包体区域网格约630万。其他区域则采用六面体的结构化网格,约700万。总共约1 470万网格。网格划分如图8~图10所示。
边界条件设置如下:上游入口为速度入口,输入均匀来流速度值;下游入口为压力出口,设定表压为0,即与参考点静压相等;外域各边界都设置为对称面;螺旋桨旋转子域的流体按照MRF方法绕X轴旋转,旋转速度为1 200 r/min;螺旋桨桨叶和桨毂相对子域的旋转速度为0,定义无滑移、不可穿透的边界条件;湍流模型选取SST
图11为外前桨的KT与纵向间距A的关系曲线,图12为外前桨的10KQ与纵向间距A的关系曲线,图13为内后桨的KT与纵向间距A的关系曲线,图14为内后桨的10KQ与纵向间距A的关系曲线。
由图可以直观地看到在A=2D时,外前桨的
$ {\eta _{\rm{B}}} = \frac{{T{V_0}}}{{2{\text{π}} nQ}}\text{。} $ |
通过表3可以看出,各方案中内后桨的推力系数、扭矩系数和船后效率均大于外前桨。通过表4的数据对比可知,随着外前桨向后移动,外前桨的
图15为外前桨的
由图可以看出,外前桨的水动力系数随基本不受其位置影响,而内后桨的推力系数和扭矩系数都随外前桨横向远离内后桨而减小。表5为5种方案的螺旋桨水动力性能参数,表6为方案B1,B2,B3,B4与原始方案螺旋桨水动力性能参数的比较。
通过表5可以看出,各方案中内后桨的推力系数、扭矩系数和船后效率均大于外前桨。通过表6可以看出,随着外前桨向内后桨横向移动过程中,内后桨的推力系数和扭矩系数均有较大幅度的提高,船后效率也逐渐增加,在B=1D时,推力系数增加8.9%,扭矩系数增加5.9%,船后效率提高2.8%;而外前桨的推力系数、扭矩系数受到其位置变化影响较小,幅度保持在0.8%以内,其船后效率基本保持不变,变化幅度保持在0.3%以内。
5 结 语1)在全附体影响下,随着外前桨逐渐靠近内后桨或者随着前后桨纵向距离变小,内后桨水动力性能基本不变,各参数变化幅度都在0.5%以内,外前桨的船后效率也基本保持不变,这说明外前桨的纵向移动对于前后桨的效率影响不大。
2)在全附体影响下,外前桨向内后桨横向移动的过程中(逐渐靠近内后桨),外前桨的水动力性能变化不大,而内后桨的
[1] |
SIMONSEN C D. Rudder propeller and hull interaction by RANS[D]. Denmark: Technical University of Denmark, 2000: 12–19.
|
[2] |
CHAO K Y. Numerical propulsion simulation for the KCS container ship[C]//Proceedings of CFD Workshop. Tokyo, 2005: 85–105.
|
[3] |
PRAKASH M S, SUBRAMANIANV A. Simulation of propeller-hull interaction using ranse solver[J]. The International Journal of Ocean and Climate Systems, 2010, 1(3): 189-208. |
[4] |
KIM J. RANS computations for KRISO container ship and VLCC tanker using the WAVIS code[C]//Proceedings of CFD Workshop. Tokyo, 2005: 105–121.
|
[5] |
吴召华. 基于体积力法的船桨舵粘性流场的数值研究[D]. 上海: 上海交通大学, 2013.
|
[6] |
TAHAPAI Y. Comparison of free-surface capturing and tracking approaches in application to modern container ship and prognoses for extension to self-propulsion simulator[C]//Proceedings of CFD Workshop. Tokyo, 2005: 145–167.
|
[7] |
CHOU S K, CHAU S W, CHEN W C, et al. Computations of ship flow around commercial hull forms with free surface or propeller effect[C]//A Workshop on Numberical Ship Hydrodynamics Proceedings, Gothenburg, 2000.
|
[8] |
VISONNEAU M, DENG G B, QUEUTEY P. Computation of model and full scale flows around fully-appended ships with an unstructured RANSE solver[C]//26th Symposium on Naval Hydrodynamics. Rome: 2006: 119–132.
|
[9] |
ZHANG N, ZHANG S L. Numerical simulation of hull/propeller interaction of submarine in submergence and near surface conditions[J]. Journal of Hydrodynamics, Ser.B, 2014, 26(1): 50-56. DOI:10.1016/S1001-6058(14)60006-8 |
[10] |
CHOI J E, MIN K S, KIM J H, et al. Resistance and propulsion characteristics of various commercial ships based on CFD results[J]. Ocean Engineering, 2010, 37(7): 549-566. DOI:10.1016/j.oceaneng.2010.02.007 |
[11] |
HAN K J, LARSSON L, REGNSTRÖM B. A numerical study of hull/propeller/rudder interaction[C]//Proceedings of 27th Symposium on Naval Hydrodynamics. Korea: 2008: 147–153.
|
[12] |
ROBERTO M, ANDREA D M. Simulation of the viscous flow around a propeller using a dynamic overlapping grid approach[C]//Proceedings of First International Symposium on Marine Propulsors. Norway: 2009: 32–49.
|
[13] |
张志荣, 李百齐, 赵峰. 螺旋桨/船体粘性流场的整体数值求解[J]. 船舶力学, 2004, 8(5): 19-26. DOI:10.3969/j.issn.1007-7294.2004.05.003 |
[14] |
沈兴荣, 冯学梅, 蔡荣泉. 大型集装箱船桨舵干扰粘性流场的数值计算研究[C]//2007年船舶力学学术会议暨《船舶力学》创刊十周年纪念学术会议论文集. 银川, 2007: 152–161.
|
[15] |
王金宝, 蔡荣泉, 冯学梅. 计及自由面兴波和螺旋桨非定常旋转效应的集装箱船舶绕流场计算[J]. 水动力学研究与进展(A辑), 2007, 22(4): 491-500. |
[16] |
沈海龙, 郑丰勇, 苏玉民. 基于滑动网格技术的船体和螺旋桨非定常干扰研究[C]//第十四届中国海洋(岸)工程学术讨论会论文集. 呼和浩特, 2009: 357–363.
|
[17] |
沈海龙, 苏玉民. 船体黏性非均匀伴流场中螺旋桨非定常水动力性能预报研究[J]. 水动力学研究与进展(A辑), 2009, 24(2): 232-241. |
[18] |
付慧萍. 船桨整体及螺旋桨诱导的船体表面脉动压力计算[J]. 哈尔滨工程大学学报, 2009, 30(7): 728-734. DOI:10.3969/j.issn.1006-7043.2009.07.002 |
[19] |
杨春蕾, 朱仁传, 缪国平, 等. 基于CFD方法的船/桨/舵干扰数值模拟[J]. 水动力学研究与进展(A辑), 2011, 26(6): 667-673. |
[20] |
覃新川, 黄胜, 常欣. 双桨两舵推进系统的水动力干扰研究[J]. 中国造船, 2008, 49(1): 112-116. DOI:10.3969/j.issn.1000-4882.2008.01.018 |
[21] |
王展智, 熊鹰, 齐万江, 等. 船后桨的布局对螺旋桨水动力性能的影响[J]. 哈尔滨工程大学学报, 2012, 33(4): 728-734. |
[22] |
王展智, 熊鹰, 姜治芳. 舵的布置对螺旋桨水动力性能的影响[J]. 华中科技大学学报, 2012, 40(8): 53-56. |