1、第45卷第2 3期2023年1 2 月舰船科学技术SHIP SCIENCEANDTECHNOLOGYVol.45,No.23Dec.,2023定子叶片对推进器水动力性能影响分析张敏革1 2，韩光哲，武天龙，王新锐，仇景!(1.天津科技大学机械工程学院，天津30 0 2 2 2；2.天津瀚海蓝帆海洋科技有限公司，天津30 0 457)摘要：为探究定子叶片对于推进器性能影响，对不同结构的推进器进行水动力分析。设计配备加速型导流罩的7 叶螺旋桨推进器，通过增加不同定子叶片以及设置对转推进器得到5种推进器。基于计算流体力学（CFD）对推进器的流场进行仿真模拟，考察不同结构推进器的水动力特性变化。结果表
2、明：在研究的转速范围内前置定子最高降低扭矩百分比仅为48.7 1%，后置定子可完美平衡到1 0 0.9 9%的扭矩；前置定子降低推力幅度最大可达36.16%，后置定子提高整体推力最高为3.32%；对转推进器在推力方面远高于泵喷推进器2 4倍，且具有扭矩自平衡的效果，流场稳定性远低于泵喷；后置定子泵喷相比无定子推进器效率提高6.2 4%，且具有更加稳定的尾流场。关键词：定子叶片；扭矩平衡；泵喷推进器；CFD中图分类号：U661.31+3文章编号：1 6 7 2-7 6 49(2 0 2 3)2 3-0 0 31-0 8Influence analysis of stator blade on h
3、ydrodynamic performance of propellerZHANG Minge2,HAN Guang-zhe,WU Tian-long,WANG Xin-rui,QIU Jing(1.School of Mechanical Engineering,Tianjin University of Science Technology,Tianjin 300222,China;2.Tianjin Hanhai LanfanAbstract:In order to explore the influence of stator blades on propeller performan
4、ce,the hydrodynamic analysis ofpropellers with different structures was carried out.A seven-bladed propeller equipped with an accelerating hood is designed.Five kinds of propellers are obtained by adding different stator blades and setting counterrotating thrusters.Based on compu-tational fluid dyna
5、mics(CFD),the hydrodynamic characteristics of propellers with different structures were investigated.Theresults show that the highest torque reduction percentage of the front stator is only 48.71%,and the rear stator can be per-fectly balanced to 100.99%torque in the speed range studied.The thrust r
6、eduction of the front stator is up to 36.16%,and theoverall thrust of the rear stator is up to 3.32%.The thrust of counterrotating thruster is 24 times higher than that of pump-jetthruster,and it has the effect of torque self-balance,and the flow field stability is much lower than that of pump-jet t
7、hruster.Compared with the stator-free propeller,the efficiency of the rear-mounted stator pump jet is 6.24%higher and the wakeflow field is more stable.Key words:stator blade;torque balance;pump jet propeller;CFD0引言泵喷推进器具有高效和低噪等特性，主要应用于军事领域。目前，除了中国、印度之外，美、英、法、俄已将其用于核潜艇推进。美国在1 9 44年提出了通过泵喷射来减少流过叶片的速度
8、以提高螺旋桨的性能。英国在2 0 世纪8 0 年代率先在“特拉法尔加”级核潜艇上安装了泵喷结构。泵喷推进器包括转子结构、导管结构、定子结构。与普通导管推进器相收稿日期：2 0 2 2-1 0-1 7作者简介：张敏革（1 9 8 0-），女，博士，高级工程师，研究方向为计算流体力学、水下推进器。文献标识码：AMarine Technology Co.,Ltd.,Tianjin 300457,China)doi:10.3404/j.issn.1672-7649.2023.23.006比，具有前置定子的推进器可改善转子进流状态，称为前置定子泵喷；具有后置定子的推进器可回收转子周向尾流所损失的动能，称
9、为后置定子泵喷 2。泵喷推进装置中的导流罩结构可分为加速式和减速式2种，加速式导管增加了流体流速，提高了总体效率；减速式导管延缓空泡的发生，改善噪声特性 3。不同的定子结构对泵喷推进器的水动力性能具有不同影响。韩瑞德等 4最早提出泵喷推进器的设计方32法及准则，利用不同翼型对转子定子进行设计。刘占一等 5用CFD对泵喷水动力性能进行数值仿真，成功预报了泵喷推进器水动力性能。胡雷俊等 6 基于数值模拟对推进器性能进行了研究，设计了对转泵喷结构，其结构是串联并且反向旋转的前后转子。谷浪 7 对鱼雷上的泵喷推进器进行分析，认为采用定子后置式可提高推进性能，同时验证了CFD方法在泵喷推进器水动力性能预
10、报方面的可行性及准确性。彭云龙等 8 对前置与后置定子泵喷的敞水性能及空化特性进行了数值模拟。Haiting等 9 通过数值模拟，分析不同定子参数在泵喷推进器推进性能和非定常力中的作用。Han等 1 0 研究前置定子的预旋作用，结果表明定子的预涡流效应改善了转子的迎面流体，使转子叶片表面波动强度较低。Negin等II对泵喷推进系统的流体动力性能进行了数值研究，结果显示在恒定进速比下随着定子厚度的增加，效率提高。武建国等 1 2 通过改变后置定子结构参数，优化泵喷推进器的效率。本文基于Ka型7 叶螺旋桨及加速性导流罩，设计多种推进器结构，系统分析了增加前置定子、后置定子、前后定子均增加以及对转推
11、进器等5种推进器的推力、扭矩、效率等特性，对比分析各结构部件及推进器整体的推力、扭矩的变化幅值，以及不同推进器尾流场变化，为泵喷推进器的研究、设计和开发提供参考。1推进器结构及尺寸泵喷推进器整体结构包括螺旋桨、前置定子、后置定子以及导流罩。其中，螺旋桨直径D为2 0 0 mm，叶片数Z为7，螺距比P/D为0.7，盘面比A/A为1；定子翼型采取NACA6412，叶数Zr均为9；导流罩采取加速性结构，直径Dh为2 40 mm。结构部件三维视图如图1 所示，2 种定子安装位置如图2 所示。通过固定螺旋桨与导管的相对位置，增加不同结构定子的泵喷推进器以及对转推进器的三维视图如图3所示。其中，对转推进器
12、是由2 个旋转方向相反的螺旋桨构成，后桨直径是前桨的0.9 倍，研究工况下，前桨转速是后桨的0.9 倍。后桨能够充分利用前桨转动产生的能量，将旋转涡流转变为轴向动能，提高推进器效率。2计算模型及模型验证2.1数学模型所有推进器的流场均为圆柱形的单相流水域，螺舰船科学技术(a)导流罩(c)前置定子图1 结构部件三维视图Fig.1 3D view of structural components前置定子(a)前置定子V1转子(b)后置定子图2 两型定子的位置Fig.2Position of two types of stator旋桨转动时流体流动不涉及热量交换，对于三维不可压流体，流体的运动遵循于
13、质量、动量守恒方程，本文使用RANS控制方程，表达式如下：apui=0,8xi第45卷(b)螺旋浆(d)后置定子转子V4后置定子(1)(4)第45卷(a)无定子推进器(c)前置定子泵喷(e)前后定子泵喷图3推进器三维视图Fig.33Dviewofpropeller(pui)(puiui)otxjxiaxilQujpu.uxi式中：t为时间；p为流体密度；u为动力粘度；p为压力；u,和u,为速度分量；pu,u,为雷诺应力；x和为笛卡尔坐标分量；S为动量方程广义源项。为了模拟复杂的流运动，同时保持计算的可靠性和计算精度，采用集合k-和k-模型优点的SSTk-作为研究水动力性能的流模型。能够准确的预
14、测出流体运动的起始和分离，适合处理边界层问题。基于紊流动能传输方程，k-模型引人了紊流脉动频率的传输方程，并对其进行了剪应力修正得到SSTk-0 湍流模型 1 3。OpkOpkuOpkyOpkwot8yokyu+kyz从+张敏革，等：定子叶片对推进器水动力性能影响分析从+oykoy(b)对转推进器式中，k=2.0，。=2.0，=5/9,=0.0 7 5，=0.0 9。端流动能生产项为：QuiPk=uiS2,S=Sii2(xjxi(6)(d)后置定子泵喷2.2物理模型利用Ansys软件对推进器的几何模型进行网格划分，采用多重参考系法（MRF）处理螺旋桨的转动问题。模型中的流域分为2 部分：包含螺
15、旋桨的动域和包含其他结构的静域，所有计算域均为圆柱形结构，转子中心到人口处的距离是3D，距离出口为7 D，静域直径为5D，减弱静域边界对于推进器的影响，所以可将推进器看作在一个无限宽的范围内运行。其流场分布如图4所示。Quixj一k从+xakz+Pk-p*kw,33Opkapwu+ot3D(2)速度入口在边界条件的设定上，推进器朝向静区域的一侧为速度入口，另一侧设置为压力出口，而在静区域的周围以及推进器表面则是没有滑动的刚性壁面，其余各面默认为内部面interior。计算所设置的压力-速度耦合算法采用SIMPLEC，此算法稳定性较好，在计算中可将亚松弛因子适当放大加速计算速度。流域采用非结构网
16、格划分，为提高推进器计算精度，设计加密块区域用来加密推进器。其中，动域网格尺寸与转子表面网格相同，加密块采用尺寸增长函数，增长比例为1.1，外表面最大值为1 0 mm，静域同样采用增长函数，增长比例为1.1，外表面最大尺寸为(3)30mm。在距离推进器较近的地方，尽可能地精细网Opwwykzkpk=一无滑移壁面2.5D静域动域图4流场分布Fig.4The flowfield distributiont+kW+-Pui7D压力入口(5)34格确保计算的准确性，在距离较远的地方设置尺寸更大的网格，以节省计算时间。网格分布如图5所示。(a)轴向截面网格分布(b)推进器表面网格分布图5网格分布图Fig
17、.5Grid distribution diagram2.3标模验证通过对Ka4-70+19A导管浆标模进行CFD仿真计算，验证模拟仿真的准确性。采取SSTk-流模型，螺旋桨导流罩边界层第一层采取0.0 0 5mm，设置速度进口及压力出口，流域边界及推进器表面为无滑移壁面，采用SIMPLEC算法进行压力速度耦合，动量方程采用二阶迎风模式，选用MRF对流场进行处理，设置各变量的收敛残差为1 0-5且监视曲线不再变化时认为收敛，标模三维视图如图6 所示。图6 标准模型Fig.6 The standard module将由Fluent计算所得数据与查找到的标模实验数据对比，判断实验值EXP与仿真值C
18、FD的误差范围是否满足要求，数据结果如表1 所示。由图7 可以看出，通过CFD计算的推力、扭矩系数与通过实验得出结果变化趋势一致，仿真值结果与实验值结果误差仅在3%以内，仿真结果数据可靠。3结果分析与讨论总推力是推进器所有结构共同产生的使潜器前行方向的力。在数据分析时，各部件的推力取值均选取前行方向一致的力，若某部件产生推力为负值，则表舰船科学技术Tab.1Data comparison resultsJKTCFD0.10.47880.20.41350.41400.50.2468 0.24900.500.450.400.350.300.25图7 仿真与实验数据对比图Fig.7 Comparis
19、on of simulation data and experimental data明该部件降低推力。扭矩统一为从推进器后面向推进器前方看去，顺时针转动方向为正，逆时针为负。若两部件产生的扭矩值符号相反，则表明2 种部件在相互平衡扭矩。所有推进器采取统一工况，进口流速为7 kn，转速范围1 0 0 0 2 0 0 0 r/min。有关的计算如下式：(7)10Kgr=100,/(on?D5),(8)KT,=T,/(on?D),(9)KT=T/(on?D4),(10)KT.Jn=Kor2元式中：V为推进器进口的速度；n为螺旋桨转速；T,为螺旋桨推力；Q,为螺旋桨扭矩；T为整个系统的推力；p为流体
20、密度；Kr为转子、导管和定子整个系统的推力系数；Kor为转子扭矩系数；KTr为转子推力系数；J为进速系数；n为推进器效率。3.1推力性能变化图8 为各推进器推力性能变化。可以看出：前置定子一直产生阻力，其阻力大小的绝对值低于螺旋桨及导流罩产生的推力绝对值，但大于后置定子产生推力的绝对值，且随着转速的增加占总推力的比例慢慢第45卷表1 数据对比结果KTEXP10KoCFD10KoEXP0.46800.44120.42290.35490.10.2进速系数/JJ=V/(nD),K T0.4387-2.31%0.42790.12%0.35060.88%KTCFDKTEXP10KOCFD10KOEX0.
21、30.4AKo-0.56%1.17%-1.23%0.5(11)第45卷-200100012001400160018002000转速/rmin!(a)无定子推进器800转子推力前置定子推力600一总推力一导管推力400力2000-200100012001400160018002000转速/rmin-1(b)前置定子泵喷800一转子推力后置定子推力600F总推力导管推力4002000-200100012001400160018002000转速/rmin-!(c）后置定子泵喷1000+转子推力800后置定子推力一总推力一导管推力600前置定子推力N/4400200010001200140016001
22、8002000转速/rmin(d)前后定子泵喷1800一前转子推力1600后转子推力1400一总推力1200力10008006004002000100012001400160018002000前转子转速/rmin!(e)对转推进器图8 推力性能变化点线图Fig.8 Point plot of thrust performance change张敏革，等：定子叶片对推进器水动力性能影响分析800一一导管推力转子推力一总推力600F400200035减小，最大可达36.1 6%；后置定子在低转速下产生阻力，随转速升高变成了推力，但占总推力比例都在3%左右，影响非常小；在高速运转状态下，转子及导流罩
23、提供主要推力，且都随转速增加而增大，转子提供的推力占比幅度逐渐减小，导流罩推力占比逐渐增加，最大可达2 4.33%。对转推进器推力由前后转子共同提供，后转子产生的推力时刻高于前转子1.2 1.44倍。相比于泵喷推进器，对转推进器在推力方面远高于泵喷推进器，最大可达4倍左右；随着转速的提高，对转推进器相比泵喷推进器增加推力的百分比随之下降。3.2扭矩性能变化图9 为各推进器扭矩性能变化。可以看出：扭矩主要由转子产生，且随转速增加而增大；导管产生的扭矩占总扭矩1%左右，微小可忽略不计；前置定子及后置定子均产生与转子相反的扭矩，随着转速增加前置定子改变扭矩的百分比逐渐下降，最大仅为48.7 1%；后
24、置定子变化趋势与转子大小接近，可达完美100.99%平衡状态。后置定子平衡扭矩的幅度高于前置定子，平衡推进器扭矩整体效率非常高。在前后定子推进器工作时，前后置定子同时产生与转子相反的扭矩，可平衡推进器整体的扭矩，但却使整体的扭矩平衡过量。随着转速增加都逐渐上升，后置定子增加的幅度高于前置定子。对于对转推进器，前后转子产生方向相反的扭矩，可达到平衡扭矩的目的；后转子产生的扭矩绝对值相比前转子大一些，二者平衡扭矩的幅度保持在75%79%之间。3.3水动力性能对比分析图1 0 分别为各泵喷总推力、总扭矩、效率对比图以及后置定子泵喷与对转推进器效率对比图。可知，总推力排名中具有前后定子的泵喷最高，其次
25、是后置定子泵喷、前置定子泵喷，而无定子推进器推力最低。总体扭矩平衡效果是后置定子泵喷强于前后定子泵喷强于前置定子泵喷。总体效率排名中后置定子泵喷为56.8 6%，高于前后定子泵喷的55.40%、前置定子泵喷的51.37%以及无定子推进器的50.6 2%。对转推进器在低转速下，效率要明显高于后置定子泵喷推进器，其中前后转子效率分别可达6 4%、6 2%。当转速达到1 30 0 r/min时，后置定子泵喷推进器效率要高于对转推进器，最高可达56.8 6%3.4流场特性对比流动能是衡量瑞流强度的一种方法，流强度36Fig.9 Point plot of torque performance chan
26、ge舰船科学技术100030前置定子推进器推力328导管扭矩转子扭矩15总扭矩10.N/X50-5-101520-35100012001400160018002000转速/rmin-1(a)无定子推进器20151050T.N/XE-5-10-15-20-25-30-35100012001400160018002000转速/r:min-l(b)前置定子泵喷881510.N/EX50-5-101520-35100012001400160018002000转速/rmin-!(c)后置定子泵喷3020.N/X100-10-20-30100012001400160018002000转速/r-min(d)
27、前后定子泵喷40前转子扭矩30一后转子扭矩一总扭矩201.N/XE100-10-20-3010001200140016001800 2000前转子转速/rmin1(e)对转推进器图9 扭矩性能变化点线图第45卷800一后置定子推进器推力一前后置定子推进器推力600一一无定子推进器推力4002000100012001400160018002000转速/rmin-1(a)总推力对比图2015105.N/X旺0-5-10-15-20转子扭矩-25前置定子扭矩-30100012001400160018002000总扭矩转速/rmin导管扭矩(b)总扭矩对比图70前置定子推进器效率65后置定子推进器效率
28、前后置定子推进器效率60无定子推进器效率555045403530一转子扭矩100012001400160018002000一后置定子扭矩转速/rmin!总扭矩(c)泵喷效率对比图导管扭矩706560%/率x5550454035100012001400160018002000转子扭矩转速/rmin危置定子扭矩(d)泵喷与对转效率对比图等管扣矩前置定子扭矩前置定子推进器扭矩后置定子推进器扭矩前后置定子推进器扭矩无定子推进器扭矩对转无导管推进器前转子效率对转无导管推进器后转子效率后置定子泵喷推进器效率图1 0 推力、扭矩、效率对比图Fig.10Thrust,torque,efficiency com
29、parison chart越大代表的流场稳定性越差。由图1 1 可知：后置定子泵喷尾部流场湍流动能为最低的0.31 6，后置定子具有尾部整流的作用，因此降低了湍流强度；前置定子泵喷尾流场的湍流动能与无定子推进器的几乎无异，可判断前置定子对于尾部湍流动能的影响很微弱；对转推进器的湍流动能最大可达1.1 1，因2 个转子同时做功，尾部流场的稳定性大大降低。由图1 2 可以看出，具有导流罩的推进器尾部呈现明显喷射状态，具有后置定子结构的推进器尾部喷射流更加集中。同时，由于后置定子可以减弱尾流切向第45卷张敏革，等：定子叶片对推进器水动力性能影响分析372:090+011.8B0+0.11.67040
30、11.46e+011.25e+016.560-011.04e4015.92e-01B:35e+005.280-018:.2604004.64e-014.18+004.00e-012:09e+003.36e-010.0084002.72e-01m/s2.08e-011.44e-017.97e-021.57e-02m/6.57e-015.93e-015.28e-014.64e-013.990-013.35e-012.70-012.06e-011.41e-017.67e-021.22e-02m/3.16e-012.86e-012.5e-012.25e-011.950-011.64e-011.34e-
31、011.03e-017.28e-024.23e-021.19e-023.92e-013.54e-013.16e-012.78e-012.40e-012.02e-011.64e-011.26e-018.800-025.00-021.21e-02m/a1.110+001.000+008.92e-017.82e-016.72e-015.63e-014.53e-013.43e-012.34e-011.24e-011.45e-02m图1 1 各推进器尾部端流动能云图Fig.11Turbulent kinetic energy cloud of each propeller tail(a)无定子推进器09
32、0+011.884011:6884011.47e+01(a)无定子推进器1:2604011.050+018:38e+006:28e:004.19e4002:090+000:000-00m/2.090+011.8Be+011:68e+011.4704011.26e+011:05e4018:3884006.280+004.1904002.090+00(b)前置定子泵喷0.00e+00m/s2.080+011.8704011:660+011.45e+011.25e4011.040+018:310+006.23e+004.16e4002:080+000.000+00m/2.090+01(c)后置定子泵
33、喷1.884011.67e+011.468:4011.2584011.04e+018:36:4006260+004.17e4002:09e+000.000+00图1 2 各推进器流场速度云图Fig,12Velocity cloud of each propeller flow field(d)前后定子泵喷诱导速度，从而减少了旋转的能量损耗，使桨毂尾部产生的毂涡降低，减少了推进器的阻力，因此具有后置定子的推进器有着更高效率。对转推进器尾部流场具有更高流速，但由于尾部端流动能过大以及喷射流比较分散，导致流场稳定性很差。4结语1）前置定子可平衡扭矩，但是会降低推力。转速(e)对转推进器越高，降低推力
34、占总推力的比例慢慢减小，平衡扭矩的幅度也逐渐降低。2）后置定子平衡扭矩的作用显著，平衡扭矩的能(b)前置定子泵喷()后置定子泵喷(d)前后定子泵喷(e)对转推进器38力高于前置定子。同时还能在高转速时提高推力，但占总推力的比例略小。产生的扭矩与转子扭矩相近，从而使整个推进器的转矩保持平衡，且随着转速增加变化趋势与转子大小相同，方向相反。3）泵喷模型中，后置定子泵喷效率最高为56.8 6%，平衡扭矩效果最好可达1 0 0.9 9%，同时具有更加稳定的尾流场。相比于泵喷推进器，对转推进器在推力方面远高于泵喷推进器，前者是后者的2 4倍左右。但效率随转速升高一直下降，流场稳定性也远低于泵喷推进器。参
35、考文献：【1】刘业宝.水下航行器泵喷推进器设计方法研究 D.哈尔滨哈尔滨工程大学，2 0 1 3【2 袁建平，王子路，王龙滟.泵喷推进器推进性能及噪声研究综述.舰船科学技术,2 0 2 2,44(6):1-7.YUAN J P,WANG Z L,WANG L Y.A review on the perfor-mance and noise of pump-jet propulsion J.Ship Science andTechnology,2022,44(6):1-7.【3】胡健，王楠，胡洋.加速导管和减速导管的性能比较.北京航空航天大学学报，2 0 1 7,43(2):2 40-2 52.
36、HU J,WANG N,HU Y.Performance comparison of accelerat-ing catheter and decelerating catheterJ.Journal of Beijinguniversity of aeronautics and astronsutics,2017,43(2):240-252.【4韩瑞德,夏长生.泵喷射推进器的设计方法及设计准则 C中国力学学会，1 9 9 9:1 2 2-1 2 7.HAN R D,XIA C S.Design method and criterion of pump jetpropeller CJ/Chin
37、ese Mechanical Society,1999:122-127.5刘占一，宋保维，黄桥高.基于CFD技术的泵喷推进器水动力性能仿真方法 J.西北工业大学学报，2 0 1 0,2 8(5):724-729.LIU Z Y,SONG B W,HUANG Q G.Hydrodynamic perfor-舰船科学技术mance simulation method of pump-jet propeller based on CFDtechnologyJ.Journal of Northwestern Polytechnical Universi-ty,2010,28(5):724729.【6
38、胡雷俊，曹琳琳，车邦祥.基于高空化性能的对转泵喷推进器设计.工程热物理学报,2 0 1 7,38(9):1 8 7 7-1 8 8 1.HU L J,CAO L L,CHE B X.Design of counterrotating pumpjet propulaion based on high altitude perfprmance J.Journalof Engineering Thermophysics,2017,38(9):1877-1881.【7】谷浪.鱼雷泵喷推进器设计及水动力性能预报方法研究 D.哈尔滨：哈尔滨工程大学，2 0 1 9.【8 彭云龙，王永生，刘承江，等.前置与
39、后置定子泵喷推进器的水动力性能对比.哈尔滨工程大学学报，2 0 1 9,40(1):132-140.PENG Y L,WANG Y S,LIU C J,et al.Comparative analysisof the hydrodynamic performance of front-stator and rear-sta-tor pump-jetsJ.Journal of Harbin Engineering University,2019,40(1):132140.9 YU Haiting,DUAN Ningyuan,HUA Hongxing,et al.Propul-sion perfo
40、rmance and unsteady forces of a pump-jet propulsorwith different pre-swirl stator parametersJ.Applied OceanResearch,2020,100:102184.10 LI Han,HUANG Qiaogao,PAN Guang,et al.The transientprediction of a pre-swirl stator pump-jet propulsor and a com-parative study of hybrid RANS/LES simulations on the
41、wakevorticesJJ.Ocean Engineering,2020,203:107224.11 DONYAVIZADEH N,GHADIMI P.Efficacy analysis ofthickness and camber size of cross section of the stator on hy-drodynamic parameters in linear jet propulsion system.Mathematical Problems in Engineering,2020:5861948.12】武建国，武天龙.泵喷推进器扭矩自平衡优化分析.舰船科学技术,2 0 2 2,44(1 2):1 0-1 5.WU J G,WU T L,Torque self-balancing optimization analysisof pump-jet propeller 3.Ship Science and Technology,2022,44(12):1015.13 EPPLER R.A Computer progam for the desig and analysis oflow-sped airfoil J.NASA Technical Memoandum,1980.第45卷
