不同工况下螺旋桨空泡的水动力性能研究

2022-06-07吴利红刘铭张文灿封锡盛李一平李硕

哈尔滨工程大学学报 2022年5期

吴利红，刘铭，张文灿，封锡盛，李一平，李硕

(1.大连海事大学船舶与海洋工程学院，辽宁大连 116026； 2.中国科学院沈阳自动化研究所机器人学国家重点实验室，辽宁沈阳 110169)

高负荷螺旋桨容易产生空泡，空泡会带来振动、噪声、侵蚀的负面影响，还会降低螺旋桨的推进效率。工作于高负荷的螺旋桨通常包括3种载体：大负荷水面舰船、水面小型快速游艇和水下机器人。其中，大负荷水面舰船螺旋桨发生空泡现象十分常见[1-3]。通常采用空泡数相等方法，通过降低环境压强，在空泡水筒中进行模型试验或基于两相流和空化模型进行数值模拟。数值研究围绕捕捉空泡图的精度，开展螺旋桨网格生成方法、加密位置、湍流模型选择等研究内容。龚鹏[4]分别采用结构网格和非结构网格对其空泡性能进行验证；叶金铭等[5]用基于面元法，比较不同网格划分计算得到的空泡形状，讨论了网格疏密程度对空泡的影响;朱志峰等[6]认为增加盘面到出口的距离和外流域使用结构网格可以提高片空化预报。胡健等[7]采用叶梢尖端螺旋加密方法捕捉叶梢空泡。Yilmaz等[8]采用自适应网格技术划分叶梢涡滑流区域，结合大涡模拟，研究了船、桨、舵之间耦合对螺旋桨梢涡发展的影响。对大负荷水面舰船船后非均匀伴流对螺旋桨空泡的发生、发展、溃灭的瞬态变化也有相关研究。刘恒等[1]和王超等[9]通过在入口处加载船艉非均匀伴流，研究螺旋桨处于不同强度伴流区域空泡性能，发现桨叶表面空泡范围和大小差异较大，进入高伴流区域空泡加剧，远离高伴流区域空泡减弱。对改善螺旋桨空泡的螺旋桨参数方面也开展了侧斜角、纵斜角、扭曲程度等方面的研究。赵旻晟等[10]采用OpenFOAM对E779A螺旋桨在斜流工况下的水动力特性及桨叶上的非定常空化进行模拟分析，结果表明，来流夹角增加会使桨叶吸力面空泡分布变得不均匀，对桨的水动力性能影响较大。相似结论在文献[11]中也表明，当纵倾向压力面弯曲时，空泡会向内半径延伸，而梢涡空泡减小。Hu等[12]对不同侧斜对梢涡空泡的影响进行了研究。

随着小型桨转速和进速提高，如水面小型快速游艇，水下高速AUV或近冰底航行AUV, 其螺旋桨负荷大大增加，空泡现象不可忽略。这种小型桨对应的临界转速和大负荷水面舰船螺旋桨有很大差别，需要进行预报，有利于避免空泡和提高安全操纵。本文对这3种载体可能发生的空化问题，理论分析了空泡初生转速[13]，结合汽液两相流和空化模型，研究了空泡形成和空泡发生后螺旋桨水动力性能曲线，探讨了纵斜角对改善空泡的作用，讨论了AUV携带螺旋桨自航的非均匀伴流场对空泡动态影响。

1 数学模型

采用多相流模型中的均质混合流模型，求解螺旋桨发生空泡后的汽液两相均质混合流的连续方程：

(1)

ρm=ρl(1-α)+ρvα

(2)

式中：α为气相体积分数；ρm为混合相的密度。

动量方程为：

(3)

式中:下标l、v分别指液相和气相；μm为混合流的动力粘性系数，定义与混合密度类似；μt.m为湍流引起的涡粘性系数，需要对混合流采用湍流模型进行求解。本文采用RNGk-ε两方程模型求解涡粘系数。

关于蒸汽质量分数f的控制方程为：

(4)

式中：Re和Rc分别为由汽化和凝结引起的相变率，f为：

(5)

空化模型采用Zwart-Gerber-Belamri模型：

当p≤pv时

(6)

当p>pv时

(7)

式中：pv为水的饱和蒸气压，与温度相关；气核半径RB默认取值1 μm；αnuc为水中所含气核体积分数。

2 物理模型

试验模型为MAU4-40标准桨，盘面比0.4，螺距比0.9，螺旋桨直径D=0.156 m，如图1为螺旋桨网格正视图。计算流域为一圆柱形流域，入口距桨盘面3D，出口距桨盘面5D，外域直径为6D，如图2所示。整个计算域包括旋转内流域①和静止外流域②，螺旋桨在旋转流域内，数值模拟时，螺旋桨和旋转内部流域一起旋转。内外流域之间用滑移界面连接，界面上的网格采用非一致连接。螺旋桨模型复杂，曲率变化较大，采用非结构网格进行划分，在桨毂和桨连接处、叶梢、导边和随边等处进行网格局部加密，以捕捉局部流场梯度变化，网格总数为1 046 147。

图1 螺旋桨网格模型

图2 螺旋桨空泡模拟的计算域

首先在无空泡的敞水工况下对MAU4-40进行螺旋桨水动力性能模拟。螺旋桨进速不变，通过改变转速，实现不同进速系数下的数值模拟，并与图谱试验值比较，如图3所示[14]。J为进速系数；KT、KQ和η0分别为螺旋桨的推力系数、扭矩系数和敞水效率。数值计算(CFD)和试验结果(EXP)两者曲线基本一致，相对偏差均在较小误差范围内，各曲线均达到了较高的预报精度。

图3 敞水性能曲线的试验结果与数值结果对比

3 空化预报

3.1 临界转速的理论预测

当螺旋桨叶片置于水下一定深度hs,来流V经过叶剖面时，随着螺旋桨高速旋转，叶剖面速度Vb增加，叶剖面上的压强pb就会降低(如图4所示)，列出叶背B点的伯努利方程，见式(8)。对应的最大减压系数ξ如式(9)。相应的最低压强达到该温度的气化压强pv，则对应无因次压强为空泡数σ。当最大减小系数ξ大于空泡数σ，空泡就能发生(p0为环境压强)。

图4 叶剖面来流

(8)

(9)

(10)

ξ≥σ

(11)

叶剖面的最大减压系数，就是叶剖面的最大压力系数Cpmin，可根据式(12)计算获得：

(12)

假设海水温度10 ℃，螺旋桨叶剖面浸没深度hs<1 m，则产生空泡的临界速度Vi为：

(13)

假设螺旋桨进速VA=6 m/s，则临界转速如式(14)所示，叶梢产生空泡的临界转速为n=3 537 r/min。随着进速提高，发生空泡的临界转速则降低，如表1所示。因此水面高速游艇，常采用多对转桨分担推力，避免空泡发生。

表1 近水面临界转速

(14)

3.2 标准大气压下的空泡模型

水面小型快速游艇的螺旋桨轴心压强为近似一个标准大气压。设环境压强为一个标准大气压，模拟MAU 4-40 桨的空泡水动力性能。假定进速为6 m/s，调整转速，实现不同进速系数下的空泡模拟，其空泡在叶背面的发生形状随进速系数的变化如图5所示。当进速系数为J=0.8，梢涡空泡初生(这与刘恒等[1]试验中梢涡空泡首先发生结论一致)，对应的临界转速为2 884 r·min-1，较理论预报值要小18%，主要因为螺旋桨各叶剖面连续作用于水，各剖面的耦合作用，使叶梢处合成速度增大，使空泡发生提前。随着螺旋桨转数的不断提高，空泡现象逐渐明显，主要位于导边至叶梢处(J=0.75,0.7)，形成局部空泡，为空泡第1阶段。当转速进一步提高(J=0.65，0.6，0.5)，片状空泡几乎布满整个叶背，空泡开始延伸到随边外围，形成第2阶段空泡，叶背维持其汽化压力，叶背的压力不再随着转速的增加而降低，则螺旋桨的推力不再随着转速的提高而增加，因而其推力也将比没有汽化压力因素时要低，如图6所示。当空泡完全布满叶背，计入空化影响的螺旋桨推力系数(J=0.5)相对无空泡情况降低28.5%。

图5 不同进速系数下螺旋桨产生空泡情况

图6 标准大气压下螺旋桨空泡后的水动力性能曲线

3.3 不同真空度的空泡模拟

大负荷水面舰艇螺旋桨，其螺旋桨尺寸大，载体航速大，通常通过满足进速系数相等，空泡数相等，在减压空泡水筒中进行模型桨空泡数值模拟和试验。即通过降低环境压强，实现不同真空度下的空泡模拟，数值模拟条件和盛振邦[15]试验条件一致，获得了一致的空泡后螺旋桨性能曲线图。图7～8分别为螺旋桨在进速系数J=0.7,0.6，空泡数分别为0.6、0.4、0.2时的空泡云图。由图7～8可见，相同进速系数下，空泡数越小，空泡越容易发生。空泡数一致时，螺旋桨载荷越大，空泡越严重，即螺旋桨合速度保持不变，提高转速，降低进速，则空泡更容易发生，表明载体螺旋桨高转速启动时，容易产生空泡。由图9的空泡性能曲线可见，低空泡数对应的临界进速相对高空泡数的大些。低空泡数更容易产生空泡；相同进速下，低空泡数的空泡水动力性能下降更大。

图7 J=0.7时不同空泡数下螺旋桨空泡图

图8 J=0.6时不同空泡数下螺旋桨空泡图

图9 不同空泡数下螺旋桨空泡性能曲线

3.4 正压下的空泡模拟

水下机器人工作深度从10～6 000 m变化, 螺旋桨直径较小。其最大电机转速也可以接近5 000 r/min, 遇到螺旋桨高负荷下，是否也会产生空泡，则需要对不同深度的环境压强下进行空泡模拟，即正压下的空泡模拟。如图10所示，深度越大，空泡数越大(这是减压空泡水筒无法模拟的)，理论上越难产生空泡。如果产生了空泡，其空泡后的螺旋桨水动力性能如图11所示。以深度500 m为例，临界空泡对应的进速系数J=0.5, 对应的转速n=4 614 r/min。

图10 深度与空泡数曲线

图11 不同深度螺旋桨空泡性能曲线

3.5 纵斜对螺旋桨空泡的性能的影响

改变纵斜角对空泡的影响结果如图12所示。在相同进速系数J=0.5时，在一个标准大气压下，采用无纵斜桨和纵斜桨(纵斜15°)，空泡图有了很大的改变，前者空泡较为严重，覆盖整个叶背；后者空泡范围减小，局限在中间到随边区域。可见，纵斜有利于改善空泡性能。这主要是因为纵斜可以减小不定常载荷通过轴系传递到载体的轴承力。