工作在艇后的推进器会产生周期性脉动的非定常力,周期性脉动非定常力是激发推进器线谱噪声的噪声源,它的预报和分析是优化推进器噪声性能的前提,同时也是预报推进器线谱噪声的基础。对推进器非定常力进行数值预报要求能够较为准确地捕捉到推进器较大的推力总量和微小的脉动量,目前已经能较为准确地对推进器定常下的性能分析,对推力值的数值分析精度已能满足工程需求[1],但是微小的推力脉动量的准确捕捉相对较困难,该方面的研究也相对较少。
随着CFD技术的快速发展,近年来许多学者致力于研究螺旋桨非定常性能数值计算方法。谭廷寿和熊鹰等[2-3]采用速度势面元法对螺旋桨非定常性能进行了数值预报,通过与试验值的对比验证了计算方法的可靠性。胡小菲等[4]采用滑移网格方法对孤立螺旋桨在给定的非均匀伴流场下的非定常力进行数值分析,研究了整体网格尺度对非定常结果的影响,沈海龙和姚震球等[5-6]将螺旋桨与船体作为一个整体,基于滑移网格方法研究了螺旋桨与船体之间的非定常干扰特性,结果表明船体对螺旋桨的影响主要是尾流场的非均匀性改变了螺旋桨的进流速度,螺旋桨对船体的影响主要是改变了船体表明的压力分布。
工作于艇后的推进器与螺旋桨工作环境类似,但是推进器模型与螺旋桨有所差别,推进器模型由导管、前置定子及转子等组成,非定常性能也将发生变化,目前对推进器非定常性能的研究很少。本文基于SST k-ω湍流模型,采用滑移网格方法导管式推进器模型的非定常力数值计算方法,并对不同网格划分方式对推进器非定常性能的影响进行研究。
1 计算方法 1.1 计算模型与网格计算模型采用文献[7]中Hughes设计的带前置定子的导管式推进器,定子采用NACA 63-018翼型,叶片数为9,转子翼型采用NACA a=0.8(mod),叶片数为4,z轴指向出口,计算模型及坐标系设置如图 1所示。整个流场计算域为圆柱体,入口距螺旋桨盘面位置约6倍螺旋桨直径,以保证来流速度到达时均匀,出口距螺旋桨盘面位置约7倍螺旋桨直径,以保证尾流充分发展,圆柱体直径约为4倍螺旋桨直径,以保证计算结果不受边界的影响。
![]() |
图 1 计算模型及坐标系 Fig. 1 Computational model and coordinate system |
整个计算域分为转子域及定子域。整个计算域采用了非结构化网格。转子域与定子域通过interface面连接。为了控制转子域附近的网格尺寸,采用Gambit中的“尺度函数”对转子叶片表面进行加密,尤其是对叶片梢部及导边进行了加密。
1.2 计算方法与边界条件泵喷推进器在艇后非均匀流场中以一定的转速绕轴旋转,属于非定常流动问题。推进器性能分析可采用的数值计算方法有相对参考坐标系、动网格及滑移网格法。相对坐标参考系只适用于推进器的定常流动模拟;动网格方法可进行推进器的非定常数值计算,但在计算过程中需要对网格进行重构,较容易在计算中出现负体积,对网格质量要求高,计算量大,对计算机硬件的要求也更高,所以目前较少使用;目前比较常采用滑移网格方法模拟推进器的非定常性能。
基于SST k-ω湍流模型,本文首先采用相对参考坐标系法对泵喷推进器的定常性能进行数值计算,待定常计算充分收敛后且计算精度较高的情况下,以定常计算结果作为初始条件采用滑移网格方法对其非定常性能进行数值分析。计算时的边界条件为:入口采用速度入口,给定来流速度,本文仅研究均匀来流下的非定常性能。出口采用为压力出口。对流项离散格式为二阶迎风格式,耗散项离散格式为二阶中心差分格式。非定常计算中时间离散格式选用二阶隐式格式。为保证残差收敛精度达到10-4,提高计算精度,压力松弛因子取为0.15,动量矩松弛因子0.3,其余采用默认设置。
非定常计算时的时间步长t按照每个时间步内推进器旋转2°进行设置,即
对导管式推进器进行定常与非定常计算稳定后,可根据式(1)~式(3)得到表征推进器性能的无因次系数,即推力系数KT与转矩系数10KQ。
$ {K_T} = \frac{T}{{\rho {n^2}{D^4}}}, $ | (1) |
$ {K_Q} = \frac{Q}{{\rho {n^2}{D^5}}}, $ | (2) |
$ J = \frac{\upsilon }{{nD}}。 $ | (3) |
式中:T和Q为推进器推力与扭矩;J为进速系数;υ为入口速度;n和D为推进器转速与转子直径;ρ为水的密度。
在采用Fluent对推进器进行非定常计算中,将推进器产生的推力按照每个时间步长t保存1次,得到推力值随时间变化的时域图。然后对其进行傅里叶变换可得到推力值在叶频处的脉动力
$ {L_p} = 20{\log _{10}}(\frac{{\tilde T \cdot \omega }}{{4\pi {c_0}{p_0}}}). $ | (4) |
由式(4)可知,为保证线谱噪声预报达到1 dB的水平,非定常力的数值预报误差不能大于12%。
2 非定常计算结果分析 2.1 计算方法验证采用相对参考坐标系方法与滑移网格方法分别对导管式推进器在进速系数分别为0.6,0.8,1.0下进行定常与非定常的计算。定常计算结果如表 1所示,试验数据来自于参考文献[7]。进速系数为0.6与0.8时,数值计算误差均在5%以下;进速系数为1.0时,计算误差有所增加,最大为6.27%。
![]() |
表 1 定常计算结果 Tab.1 Steady numerical results |
以定常计算结果为初始条件对J=0.6下导管式推进器采用滑移网格方法进行非定常计算,取时域上的平均值计算KT与10KQ,结果如表 2所示。非定常计算结果与试验值相比,推力系数与扭矩系数误差均在10%以下,与定常计算结果相比,推力系数计算精度稍低,但精度已能满足要求,扭矩系数精度得到提高。
![]() |
表 2 非定常计算结果 Tab.2 Unsteady numerical results |
在研究推进器非定常性能时需要对推进器全流道进行数值计算,当转子流道内的网格采用整体划分时,各流道之间的网格无法保证完全一致,各个叶片受力受叶片网格的不均匀性的影响。非定常计算需要同时捕捉推进器总力和微小的脉动值。由于各个叶片受力不同,可能对推进器的总力影响不大,但是会影响非定常力的脉动值。为保证转子各流道之间的网格是完全一致的,采用周期性的网格进行非定常分析,即首先生成转子单流道的计算域网格,再将单流道网格旋转复制得到全流道的计算网格。
为研究各流道网格不均匀性对非定常力的影响,分别对推进器在整体网格与周期性网格下的非定常性能进行数值分析,2套网格保证网格尺度一致,消除网格尺度对计算结果的影响。
2.2.1 计算网格对推进器性能的影响图 2与图 3分别为2套网格下在J=0.6时非定常计算推力系数及扭矩系数的时域图。由结果可知,2套网格下非定常计算时推力系数与扭矩系数随时间变化的规律一致,计算到达稳定的时间基本一致,主要差别在数值大小上。
![]() |
图 2 两套网格下的推力系数KT时域图(J=0.6) Fig. 2 The time domain distribution of thrust coefficient (KT) under two sets of meshes (J=0.6) |
![]() |
图 3 两套网格下的扭矩系数10 KO时域图(J=0.6) Fig. 3 The time domain distribution of torque coefficient (10 KO)under two sets of meshes (J=0.6) |
表 3为周期网格下非定常计算后推力系数与扭矩系数值。与表 1和表 2中的数据相比,周期性网格下定常计算推力系数与扭矩系数计算精度有所降低,然而周期性网格下非定常计算推力系数误差为4.95%,推力系数计算精度得到提高,误差较小。2套网格下定常计算时扭矩计算精度比推力计算精度低,但非定常计算下扭矩计算精度得到大幅提高,推力计算精度降低,使得扭矩计算精度比推力计算精度高。
![]() |
表 3 周期性网格下的计算结果(J=0.6) Tab.3 Numerical results under periodic mesh (J=0.6) |
图 4与图 5分别为2套网格下垂向力Tx与径向力Ty的时域分布图。Tx为推进器X方向的受力,Ty为推进器Y方向的受力。从图中可知,Tx和Ty在2套网格下随时间变化的规律一致,计算稳定所需时间基本一样,主要差别在于脉动幅值。Tx与Ty在周期性网格下的脉动幅值明显比整体网格低。
![]() |
图 4
两套网格下垂向力 |
![]() |
图 5 两套网格下径向力Tv的时域分布图 Fig. 5 The time domain distribution of radial force Tv under two sets of meshes |
基于图 2、图 4与图 5所示的时域分布图,对垂向力Tx、径向力Ty与轴向推力Tz进行傅里叶变换,可分别得到各阶叶频下的非定常脉动力
![]() |
表 4 两套网格下x、y及z方向上前3阶叶频处的非定常力 Tab.4 The unsteady force at the first three order blade frequency in x, y, z direction under two sets of meshes |
$ {{\tilde T}_i} = \sqrt {{{\tilde T}_{i1}}^2 + {{\tilde T}_{i2}}^2 + {{\tilde T}_{i3}}^2} 。 $ | (5) |
式中:下标i为力的方向x、y、z;
由表 4可知,整体网格下轴向平均脉动力为4.23 N,周期性网格下的轴向平均脉动力为4.218 N,2套网格下的脉动强度基本一致。
整体网格下x方向上的平均脉动力为3.549 N,y方向上的平均脉动力为5.249 N;周期性网格下x方向上的平均脉动力为2.07 N,y方向上的平均脉动力为2.974 N。与周期性网格相比,整体网格下x方向与y方向上的脉动较为剧烈,这即是由于整体网格下各个叶片受力不一致导致额外增加了x方向与y方向上的脉动值。
3 结语本文基于滑移网格方法对带前置定子导管式推进器非定常性能数值计算方法进行了研究,主要研究了推进器在不同网格划分方式下的非定常性能。可得到以下结论:
1) 与定常计算结果相比,非定常计算下推进器的推力系数计算精度较低,但是能满足要求;
2) 与整体网格相比,周期性网格下非定常计算后推力系数计算精度得到提高;
3) 与整体网格相比,周期性网格下非定常轴向脉动力基本相同,x方向与y方向上的非定常脉动力明显降低。
[1] | 潘光, 胡斌, 王鹏, 等. 泵喷推进器定常水动力性能数值模拟[J]. 上海交通大学学报 , 2013, 47 (6) :932–937. |
[2] | 谭廷寿, 贺伟. 螺旋桨非定常轴承力计算[J]. 船海工程 , 2006, 35 (2) :42–46. |
[3] | 熊鹰.非均匀流中螺旋桨空泡及脉动压力的数值和试验研究[D].武汉:武汉理工大学, 2002. http://cdmd.cnki.com.cn/article/cdmd-10497-2003095132.htm |
[4] | 胡小菲, 黄振宇, 洪方文. 螺旋桨非定常力的黏性数值分析[J]. 水动力学研究与进展 , 2009, 24 (6) :734–739. |
[5] | 沈海龙, 苏玉民. 船体黏性非均匀伴流场中螺旋桨非定常水动力性能预报研究[J]. 水动力学研究与进展 , 2009, 24 (2) :232–241. |
[6] | 姚震球, 高慧, 杨春蕾. 基于滑移网格的带螺旋桨艇体尾流场数值分析方法[J]. 江苏科技大学学报(自然科学版) , 2008, 22 (2) :15–20. |
[7] | HUGES M J, KINNAS S A. An analysis method for a ducted propeller with pre-swirl stator blades[C]//Proceedings of propeller/shafting's 91 symposium Virginia beach. Virgina, USA:SNAME, 1991. |