泵喷推进器通常由导管、定子(导叶)和转子(叶轮)构成[1]。由于导管的存在其工作流场相比喷泵工作内流场和螺旋桨工作开式流场更为复杂,是一个复杂多连通区域,各部件之间以及推进器与推进载体之间的相互作用会直接影响到推进器的性能[2-3]。近年来,泵喷推进器(简称泵喷)以推进效率高、辐射噪声低、抗空化能力强和临界航速高[4]等优点广泛应用于水下航行器,尤其是潜艇和鱼雷推进。如美国的“Mk48”鱼雷、海狼级弗吉尼亚级攻击核潜艇,英国的机敏级攻击核潜艇等。
当今电力电子技术的飞速发展随之推动了一种新型水下推进装置的诞生——集成电机推进器(integrated motor pump-jet,IMP)。这种IMP泵喷主要有以下优点:1) 原动机与推进器之间不再有机械传动装置,水下辐射噪声的主要贡献源之一减速齿轮箱被省去。电机转子与叶片无缝连接,消除了螺旋桨和机械泵喷都有可能产生的稍涡空泡噪声,提了高声隐身性;2) 避免了输出轴处动密封等摩擦耗功损失,提高推进系统效率。此外IMP始终工作在海水环境中电机产生热量可通过海水自然冷却;3) IMP可以作为独立的推进模块直接安装在推进主体上提高了推进系统的可靠性。目前,国内外对IMP已有相关研究但可查阅文献较少,主要集中在CFD仿真方面。1993年Luciano Verinesi对水下航行器的吊舱式IMP结构做了详尽地介绍[5]。杨琼方运用计算流体力学(CFD)手段分析了“机电一体化”导管桨的水力性能并根据计算结果提出改进优化思路[6]。潘光采用雷诺时均方法对某泵喷推进器进行了定常水动力性能模拟[7]。何东林初步探索了集成电机泵喷推进器技术并进行数值模拟[8]。
经过近三十年的发展,CFD已经普遍应用于研究机构、生产厂家的新型推进器研发。2010年美国海军利用计算流体力学研究了浸没式喷水推进器的推进性能,分析了流道前不同进流面的选取对推进效率的影响[9]。Park等用RANS方法研究了带导叶的导管桨性能和流场[10]。刘业宝建立了环流理论与泵相似理论相结合的泵喷推进器设计方法,并应用CFD软件对设计泵喷推进器性能进行验证[11]。本文借助计算流体力学方法对相同工作条件下常规泵喷与无轴IMP泵喷的工作流场进行定常数值模拟,对比分析了两种泵喷推进器的水力性能差异,为后续IMP无轴泵喷设计人员的水力设计提供参考。
1 数值计算方法验证IMP结构示意图如图 1,电机转子与叶轮叶片叶稍面无缝连接,定子安装在导管内,定、转子之间存在一定的间隙(水隙)来保证静止部件与旋转部件之间的润滑及冷却。
泵喷推进器主要做功部件为叶轮(转子),在叶轮的后部会安置导叶(定子)对叶轮出流进行整流,或将导叶置于叶轮前部为叶轮进流提供一定的预旋。喷水推进泵的几何结构与泵喷类似(图 2),工作原理与泵喷相同都是通过叶轮对水做功将机械能转换为水流的动能和压能高速喷出,产生推力推动推进载体前进。本文以某混流式喷泵为对象进行水力性能分析,验证数值计算方法的可信性。
该混流式喷水推进泵标称直径710 mm,叶轮6叶,导叶11叶。叶轮单通道网格数20万,导叶10万,泵内流场网格总数约300万,均为六面体网格。叶轮采用J型拓扑结构,导叶采用H拓扑,叶片周围为O型网格,叶顶间隙用独立的H型网格嵌入到周围的O型网格之中[12]。文中所述计算模拟均以控制域内流动为定常,不可压缩流动为假设前提,动静交界面采用稳态多参考系方法来处理。喷泵壁面Y+值保持在185以下,满足湍流模型求解要求。边界条件设置为总压进口,静压出口,参考压力为1个大气压。该泵几何以及叶轮、导叶网格如图 3所示。
采用SST模型、Realizable k-ε模型和标准k-ε模型3种湍流模型预报了在不同转速下混流泵轴功率值并与厂家提供实测值比较(见表 1),结果表明SST湍流模型与实测值最为接近。这是因为切应力输运(shear stress transport,SST)模型结合了k-ε和k-ω两种湍流模型的优点,能够较好地模拟近壁面区域流动[13]。3种湍流模型的功率计算值都小于试验值的原因是由于数值计算没有考虑机械损失。此案例间接证明了本文选用湍流模型的适用性和所用数值方法的可信性。
转速/(r·min-1) | 厂家提供值/kW | SST误差/% | Realizable误差 | 标准k-ε误差 |
0.80 | 0.505 | -0.61 | -1.28 | -2.41 |
0.91 | 0.757 | -0.72 | -1.76 | -2.03 |
1 | 1.007 | -0.66 | -1.27 | -1.97 |
1.02 | 1.075 | -1.27 | -0.73 | -1.09 |
1.14 | 1.479 | -1.97 | -1.09 | -1.07 |
研究对象常规泵喷为作者所在研究团队研发的某型泵喷。泵喷轴向长度1.7 m,导管进口直径D =2.6 m,叶轮为9片叶片,导叶叶片数为11,叶顶间隙为5 mm。常规泵喷几何见图 4,为前置导叶后置叶轮结构。图 4中导管内部深色部分为旋转部件叶轮,其余浅色部分为静止部件包括导叶和导管。泵喷导管是以荷兰船模试验池公开导管No.19A为基础进行了部分改进优化的19A改进型。
采用计算流体力学(CFD)方法,对泵喷施加均匀进流流场并进行数值模拟得到其敞水性能,计算域见图 5。以圆柱体包裹泵喷推进器,泵喷进口距进流边界为3D;泵喷喷口距出流边界为5D;直径为3D,轴线与泵喷对称轴重合。
根据模拟流动现象的物理性质,将计算域分为叶轮、导叶、喷口和外域4部分。其中,叶轮为旋转域,导叶、喷口和外部流场计算域为静止域。采用分块六面体结构化网格对各部件进行网格划分。由网格无关性计算曲线(图 6)可知当泵喷内部计算域网格达到330万时泵喷的基本水力性能指标杨程和效率不再有较大变化,即达到流场求解精度要求。此时叶轮单通道网格18.5万,导叶12.5万,壁面第一层网格厚度保持在0.02 mm,控制域总网格节点数690万。
图 7是泵喷整体网格和叶轮叶稍与导管内壁间隙附近网格。计算过程中同样采用SST湍流模型封闭方程组,得到壁面Y+值小于100(图 8)满足SST湍流模型对壁面流动的模拟要求。
图 9显示6 kn均匀来流常规泵喷内外流场流线,可以看出流线基本光顺无漩涡,说明导叶和叶轮匹配良好。常规泵喷在6节均匀进流时的水力性能见表 2。本文以流量Q、扬程H、功率P和效率η等泵的主要参数指标来表征泵喷的水力性能,其中流量Q为单位时间内流经泵喷流体的质量流量;扬程为导管出口与进口流面流体总压之差;力矩N为泵喷内所有旋转部件产生的力矩;轴向分量F表示泵喷射流在出口处水流轴向能量占总能量的百分比;功率P为叶轮、轮毂等旋转部件消耗轴功率:
$P = \frac{{n \times 2\pi \times N}}{{60 \times 1000}}\left( {{\rm{kW}}} \right)$ | (1) |
式中:N为力矩,N·M;n为泵喷转速,r/min;效率η为泵喷吸收轴功率对流体做功使其获得能量增加值占轴功率的比例,效率越高表明泵喷做功过程中的水力损失越少:
$\eta = \frac{{Q \times g \times H}}{{P \times 1000 \times 1000}}$ | (2) |
式中:Q为质量流量,kg/s;g为重力加速度;H为扬程,m;P为轴功率,kW。
转速(r·min-1) | 流量/(kg·s-1) | 扬程/m | 力矩/(N·M) | 功率/kW | 效率/% | 推力/N | 导管/N | 轴向分量/% | 叶轮效率/% |
65 | 12 498 | 0.52 | 10 757.6 | 73.19 | 86.99 | 16 149.1 | 3 707.5 | 99.41 | 91.82 |
由表 2可知该常规泵喷导管产生正向推力为推力导管,在设计转速设计航速下常规泵喷效率超过85%,导叶、导管、喷口等产生的摩擦损失约为5个百分点。图 10所示为叶轮叶稍附近流线和稍涡,可见计算很好地扑捉到流场中的流动细节。在间隙处,流线从压力面被吸到吸力面,稍涡内大量的旋转能量造成叶稍附近流动损失较大。加速导管机械式泵喷推进器中叶轮叶片的梢涡空泡经常成为泵喷推进器噪声的贡献源之一同时也会导致泵喷水力性能的显著下降,本研究对象经Rayleigh-Plesset空化模型计算验证设计航速转速下泵喷内部并未出现空泡。
2.2 IMP泵喷水力性能预报IMP泵喷与常规泵喷相比结构上最大的区别在于其导管内加装电机,叶轮叶片与电机转子连接为一体消除了叶顶间隙,避免了叶稍间隙流损失,但也因此增加了气隙的存在。为确保两种泵喷推进器对比分析条件的可比性即叶轮、导叶和导管的叶型、几何形状等均保持相同,作者在原常规泵喷导管内增加电机转子设置气隙来模拟电机的存在如图 11,气隙宽度15 mm。需要说明的是实际IMP泵喷推进器加装电机的导管厚度较常规泵喷是有很大增加的,这会大大增加导管的迎水面积增大导管阻力,文中对此导管厚度变化带来的影响暂没有考虑。
IMP泵喷水力性能预报的计算域大小和边界条件设置与2.1节所述常规泵喷相同,主要不同体现在旋转部件除叶轮轮毂、叶片外还增加了电机转子。为精确捕捉气隙内流动,采用结构化网格对气隙内流通区域进行网格划分并加密处理,气隙内壁面第一层网格0.08 mm,厚度增长率1.2,全通道网格节点数78万(图 12)。泵喷壁面Y+值最大值45.8满足湍流模型适用要求。
对CFD计算结果进行后处理,取得间隙流场剖面速度矢量如图 13所示(只显示叶轮出流区域)。图中深色箭头部分为贴近电机转子壁面位置,因水流粘性力作用其速度值较大,方向与电机转子外壁面旋转方向切向矢量相同。气隙内流动为逆向流动与泵喷整体工作流场相反,此结论与文献[6]一致。表 3为6 kn均匀来流情况下IMP泵喷水力性能。图 14显示为流经泵喷内外流场和气隙的流线图,对比表 2和表 3轴向分量数值和内部流线图可以看出,与常规泵喷相比由于气隙的存在泵喷出流旋转分量增加。
转速/(r·min-1) | 流量/(kg·s-1) | 扬程/m | 力矩/(N·M) | 功率/kW | 效率 | 推力/N | 导管力/N | 轴向分量/% | 叶轮效率/% |
65 | 12 443.9 | 0.51 | 11 066.5 | 75.29 | 81.98 | 15 673 | 3 287.7 | 99.22 | 87.37 |
3 性能对比分析
为分析对比两种泵喷水力性能,扩展计算了不同进流速度下泵喷敞水性能,对比曲线如图 15。可以看出,相同进流速度条件下IMP泵喷相比常规泵喷流量减小、消耗功率增大,扬程和泵效率亦有所降低。这是由于相同进流情况下,电机定子与转子之间气隙的回流造成了IMP泵喷实际流量的减小、扬程降低,而电机转子作为旋转部件,其较大的旋转半径也大大增加了转矩,进而功率增大效率降低。
图 16中对比了常规机械式泵喷与IMP泵喷在6节进速时导管内外壁的压力分布,其中横坐标z为轴向进口到出口空间位置,纵坐标为压力数Cp,
${C_p} = \left( {p - {p_a}} \right)/\left( {0.5\varepsilon {v^2}} \right)$ | (3) |
式中:pa为参考压力,v为自由来流速度。常规泵喷与IMP泵喷外壁面压力分布几乎重合,IMP内壁面存在两个压力阶跃,即对应叶轮进口处气隙位置和叶轮出口处气隙位置,叶轮出口处压力较高进口处压力较低。这与图 13间隙内流场在逆向压力作用下为逆向流动相吻合。因为IMP泵喷导管内壁面在叶轮对应纵向位置为电机定子内壁面,所以没有出现像常规泵喷壁面压力持续增加而是在气隙进出口出现阶跃,导管内壁面其他位置压力分布大致相同。
4 结论本文通过经验证的数值计算方法对常规机械式泵喷推进器和IMP泵喷推进器设计转速下的水力性能进行预报,成功捕捉到IMP电机转子与定子之间气隙内逆向流动现象,并对比不同进速下其存在的差异,分析得出了噪声指标等其他性能优越的IMP泵喷水力性能下降的主要原因和变化趋势,为新型推进器的性能设计及优化提供参考。本文结论主要有以下几点:
1) 由于气隙的存叶轮出口与叶轮进口形成局部“短路”以至进出口压力下降,泵喷叶轮做功能力减弱,相同进速转速下IMP泵喷敞水效率比常规机械式泵喷低约5%。
2) IMP泵喷电机转子与叶轮叶片直接相连同为旋转部件,所以旋转部件摩擦表面积约为常规机械式泵喷的2.2倍,力矩增大摩擦损失增加,相同进速转速下消耗功率增大。
3) 相同转速不同进速下,随着进速的增加IMP泵喷与常规泵喷流量、推力差值逐渐减小,功率差值逐步增大。
4) IMP泵喷从安装结构和工作原理上分析其噪声性能将比机械式泵喷有显著提高,但会牺牲部分水力性能指标。IMP噪声性能相比常规机械式泵喷降低的量值有待进一步对比分析。
[1] |
李锡群, 王志华. 电机/推进器一体化装置(IMP)介绍[J].
船电技术, 2003, 23(2): 5-6–31.
LI Xiqun, WANG Zhihua. Introduction to integrated motor/propulsor (IMP)[J]. Marine electric & electronic technology, 2003, 23(2): 5-6–31. |
[2] | 王天奎, 唐登海. 泵喷推进器-低噪声的核潜艇推进方式[J]. 现代军事, 2006(7): 52–54. |
[3] |
王涛, 周连第, 张鑫. 轴对称体与导管推进器组合体的三维复杂流场的计算与分析[J].
船舶力学, 2003, 7(2): 21–32.
WANG Tao, ZHOU Liandi, ZHANG Xin. Numerical simulation of 3-D integrative viscous complicated flow field around axisymmetric body with ducted propulsion[J]. Journal of ship mechanic, 2003, 7(2): 21–32. |
[4] | 杨龙塾. 潜艇声隐身技术[J]. 现代舰船, 2009(11): 23–25. |
[5] | VERONESI L. Integral motor propulsor unit for water vehicles: US, 5220231[P]. 1993-06-15. |
[6] |
杨琼方, 王永生, 范露. CFD在"机电一体化"导管桨性能分析及优化设计中的应用[J].
机械工程学报, 2010, 46(1): 162–168.
YANG Qiongfang, WANG Yongsheng, FAN Lu. Hydrodynamic performance analysis and optimization design of mechatronic ducted Propeller with CFD[J]. Journal of mechanical engineering, 2010, 46(1): 162–168. |
[7] |
潘光, 胡斌, 王鹏, 等. 泵喷推进器定常水动力性能数值模拟[J].
上海交通大学学报, 2013, 47(6): 932–937.
PAN Guang, HU Bin, WANG Peng, et al. Numerical simulation of steady hydrodynamic performance of a pump-jet propulsor[J]. Journal of Shanghai jiaotong university, 2013, 47(6): 932–937. |
[8] | 何东林. 集成电机泵喷推进器技术研究[D]. 西安: 西北工业大学, 2005. |
[9] | DINHAM-PEREN T. CFD and submerged waterjets[J]. The Naval Architect, 2010(7/8): 34–35. |
[10] | PARK W G, JANG J H. Numerical simulation of flow field of ducted marine propeller with guide vane[C]//Proceedings of the 4th international conference on pumps and fans. Beijing, China, 2002: 307-314. |
[11] | 刘业宝. 水下航行器泵喷推进器设计方法研究[D]. 哈尔滨: 哈尔滨工程大学, 2013. LIU Yebao. Study on design method of pump jet thruster for underwater vehicles[D]. Harbin: Harbin Engineering University, 2013. |
[12] |
彭云龙, 王永生, 靳栓宝. 轴流式喷水推进泵的三元设计[J].
中南大学学报:自然科学版, 2014, 45(6): 1812–1818.
PENG Yunlong, WANG Yongsheng, JIN Shuanbao. Three-dimensional inverse design method applied to waterjet axial-flow pump[J]. Journal of central south university:science and technology, 2014, 45(6): 1812–1818. |
[13] | LLOYD G, ESPANOLES A. Best practice guidelines for marine applications of computational fluid dynamics[R]. WS Atkins Consultants and Members of the NSC, MARNET-CFD Thernatic Network, 2002. |