汲国瑞，刘雪琴，尹晓辉，刘建国

(1.中国船舶及海洋工程设计研究院，上海 200011；2.喷水推进技术重点实验室，上海 200011)

0 引 言

喷水推进是一种广泛应用在中高速船舶、具有时代特征的先进推进方式[1]。与传统螺旋桨推进相比具有很多优点，但是在喷水推进船设计阶段，只能通过相似换算将模型实验数据换算为实尺度喷水推进船的推力和阻力，从而实现快速性预报。在此过程中尺度效应对快速性预报的准确性具有重大影响。因此研究尺度效应对喷水推进器进出口流场的影响对喷水推进船的快速性预报具有重要意义。

1 研究对象

以某喷水推进船为研究对象，该喷水推进船安装两台混流式喷水推进泵，设计航速Vs为32 kn，喷泵设计转速ns为372 r/min，设计流量Qs为30 m3/s。

图1 喷水推进船Fig.1 Water jet propulsion ship

2 相似换算

船模航速可以根据傅汝德相似得出：

在满足进速系数相等且假定伴流无尺度作用时，叶轮转速满足以下相似关系：

喷泵流量和转速满足的相似公式如下：

根据以上相似关系式可以求出不同缩尺比的船速及喷泵参数如表1 所示。

表1 不同缩尺比喷泵参数Tab.1 Parameters of pump with different scale ratio

由此可以得到不同缩尺比下喷泵的参数，进而可以对不同缩尺比下的喷泵进行数值模拟。

3 叶轮和导叶体模拟方法

叶轮和导叶体是喷水推进系统重要水动力部件，在数值模拟喷水推进船自航状态时对于叶轮和导叶的数值模拟通常采用滑移网格法、旋转坐标系法、虚拟盘、力源法4 种方法[5]，其中滑移网格法和旋转坐标系法需要对叶轮和导叶体全部划分网格，对计算机性能要求较高，且耗时较长。而力源法没有考虑水流的旋转，与实际流场相差较大，因此决定采用虚拟盘方法对叶轮和导叶体进行数值模拟。

喷水推进器叶轮的扭矩系数和的推力系数定义如下：

转矩系数

式中：Q为转子转矩，N·m；ρ为水的密度，kg/m3；n为转速，r/min；D为喷泵特征直径，m。

根据相似关系，不同缩尺比的喷泵的转矩系数相等。

推力系数

式中：T为喷水推进器产生的推力。

与螺旋桨推进船不同，喷水推进船喷水推进器产生的推力要大于叶轮旋转产生的轴向力，因此在计算喷水推进器的推力系数时，推力T取喷水推进器产生的推力，也就是船舶的阻力值，因此不同缩尺比喷泵的KT是不同的，需要具体计算如表2 所示。

表2 不同缩尺比转子扭矩系数和推力系数Tab.2 Torque and thrust coefficients of propeller with different scale ratios

4 数值计算

数值计算时只计算一半船体和1台喷泵，采用两相流计算，计算域的长度为6L（L为船长，船首之前计算域长度为2L），宽度为L，高度为1.3L（水线以上部分高度为0.3L）。

图2 计算域Fig.2 Calculation domain

4.1 网格划分

在Star CCM 中划分网格，并且在划分网格时将其分为船体域、叶轮域、导叶体域3 个部分，叶轮域和船体域采用Trim 网格，导叶体域采用多面体网格。在水线面附近和喷口以后区域网格进行加密。

图3 计算模型划分网格Fig.3 Mesh of computing model

4.2 边界条件

计算模型入口设为速度入口，速度为船速；计算域出口设为压力出口；计算域2 个侧面设为对称边界；计算域底部设为壁面；计算域顶部设为对称边界；在叶轮域添加虚拟盘，并设定好不同缩尺比下叶轮的KT，KQ和转速。

图4 转子域虚拟盘Fig.4 Virtual disk of propeller domain

4.3 计算结果

计算得到不同缩尺比喷口处流量与相似换算得到的流量误差如表3 所示。

表3 喷口处流量计算误差Tab.3 Calculation error of flow rate at nozzle

4.3.1 进口

进口形状采用椭圆形，其位置取在进口流道与船底板切点前方1D处，椭圆的长半轴由流线追踪法确定，椭圆的短半轴由质量守恒确定，从而得到不同缩尺比下椭圆的形状，进而可以计算出边界层影响系数α。对不同缩尺比下椭圆的长半轴和短半轴无因次化如表4所示。

表4 进口计算结果Tab.4 Results at inlet

对椭圆的长半轴和短半轴无因次化后可以看出，不同缩尺比下椭圆的长半轴基本在1 附近，而椭圆的短半轴会随着缩尺比的增加而不断增加，无因次化后的椭圆面积也会随着缩尺比的增加而不断增加，即缩尺比越大，获流区也会相应增加。不同缩尺比的边界层影响系数相差不大。

定义喷水推进船实泵产生的推力如下：

式中：Fs为实泵产生的推力；ρs为实船航行时水的密度；Vjs为实泵喷口速度；αs为实船边界层影响系数。

定义缩尺比为λ 的模型泵产生的推力如下：

式中：Fλ为模型泵产生的推力；ρλ为实船航行时水的密度；Vjλ为模型泵喷口速度；αλ为实船边界层影响系数；V0λ为模型船航速。

由以上2 个公式可得：

由喷泵相似换算式（3），上式可以转化为：

定义无因次系数

则式（9）可以转化为：

式中：Cs为实船C值；Cλ为缩尺比为λ 时C值。

由上式可以看出，在将实验测量的缩尺比为λ 的喷水推进泵产生的推力换算为实泵产生的推力时，其不仅与功率和缩尺比有关，与无因次系数C也有关系。根据以上计算结果，不同缩尺比的喷水推进泵无因次系数C值如表5 所示。

表5 不同缩尺比下C 值Tab.5 C of different scale ratio

将以上数据绘制成图5。可以看出，随着缩尺比的增加，C值也不断增加，但是增加的速度会逐渐变缓，C值与缩尺比λ 的关系可以用四次多项式拟合。通过拟合出的关系式就可以计算出缩尺比为1 时的C值，进而换算出实泵产生的推力。

图5 不同缩尺比C 值Fig.5 C of different scale ratio

将缩尺比为λ 喷水推进器进口速度Vin无因次化，不同缩尺比下Vin/V0云图如图6 所示。

图6 不同缩尺比Vin/V0 云图Fig.6 Vin/V0 contour of different scale ratio

可以看出，进口Vin/V0云图分布大部分区域都小于1，即进口的速度分布中大部分区域都小于船速。缩尺比对进口的速度分布影响不大。

4.3.2 出口

将缩尺比为λ 喷水推进器出口速度Vjλ无因次化，不同缩尺比Vjλ/V0λ云图如图7 所示。

图7 不同缩尺比Vjλ/V0λ 云图Fig.7 Vjλ/V0λ contour of different scale ratio

受导流帽影响，喷口Vjλ/V0λ云图中间区域数值较小，而边侧区域速度相对较高。缩尺比对出口的速度影响不明显。

定义出口不均匀度系数如下式：

式中：Q为截面流量；u为截面上任一点的速度；U为截面平均速度。

不同缩尺比喷水推进出口不均匀度系数如表6 所示。可以看出，不同缩尺比下喷水推进出口不均匀度基本相同。

表6 不同缩尺比下出口不均匀度Tab.6 Non-uniformity at nozzle of different scale ratio

5 结 语

通过以上计算分析得出以下结论：

1）将模型尺度喷水推进装置推力换算到实泵时，引入了无因次系数C，该系数会随着的缩尺比的增加而不断增加。可以通过计算不同缩尺比时C 值，拟合出一个C 值与缩尺比的曲线，进而求出实泵的C 值，然后通过相似换算即可得出实泵产生推力。

2）在缩尺比较大时，缩尺比对喷水推进进口和出口速度分布以及出口不均匀度影响不大。