项 松,刘远强
(沈阳航空航天大学 辽宁省通用航空重点实验室,沈阳 110136)
某型三叶螺旋桨的气动特性数值模拟及试验
项 松,刘远强
(沈阳航空航天大学 辽宁省通用航空重点实验室,沈阳 110136)
以某型三叶螺旋桨为研究对象,基于RANS(Reynolds-averaged Naiver-Stokes)方程和SST(Shear Stress Transport)湍流模型的多重参考坐标系MRF(Multiple Reference Frames)方法对该三叶螺旋桨进行准定常数值模拟和性能计算。通过与试验结果对比,对三叶螺旋桨不同转速下的静态拉力、扭矩和效率进行了验证分析,得到拉力偏差值在2%左右,扭矩偏差值在10%左右。经过比较发现,计算结果与试验结果吻合良好,可为通航飞机螺旋桨的模拟和设计提供参考。
三叶螺旋桨;气动性能;数值模拟;风洞试验;通用航空
目前国内外文献主要是采用经验公式的方法进行性能的计算,虽然计算速度快,但是精度不高,而且许多方法对于零风速下的螺旋桨性能预测偏差范围更大。本文将对三叶通航飞机螺旋桨采用基于雷诺平均NS方程的多重参考坐标系MRF方法对其滑流进行准定常数值模拟,分析螺旋桨在零风速下的流场特性和相关性能,并按照等前进比原则进行了螺旋桨的静态拉力试验,验证计算值与试验值的吻合程度较好,可为今后通航飞机螺旋桨气动计算提供参考。
1.1 控制方程
(1)
式(1)中t为时间项,V为流体微元单元控制体,S为围绕单元控制体的封闭曲面,Q为守恒变量,H为无黏通量,Hv为黏性通量项和坐标系转换的添加源项G,其具体求解公式如下:
QT=[ρ,ρu,ρv,ρw,ρE]
(2)
(3)
(4)
(5)
在式(2)~(5)中,Ix,Iy,Iz分别是绝对坐标系下坐标轴方向的单位向量,τij为黏性应力量,ρ,q,qb,u,v,w,E,H,p,ω分别为流体的密度、流体绝对速度、网格速度、绝对速度在旋转坐标系下的3个分量、总能、总焓、压力和旋翼的旋转角速度矢量。f5,g5,h5的表达式如(6),其中k和T分别为热传导系数和温度。
(6)
1.2 湍流模型
虽然前列腺等离子电切技术安全性高于普通单极电切,但是绝大多数前列腺增生患者为高龄人群,且多合并各种基础疾病。因此,安全有效地使用该技术对确保医疗安全、提高前列腺增生患者的生活质量具有重要意义。
本文计算采用基于k-ω的SST湍流模型来使平均运动方程封闭以便求解。基于k-ω的SST湍流模型解决了湍流剪切应力的传输,同时又在湍流开端和在负压梯度下产生的气流分离进行了高度准确的预测,具有较好的普适性。
两方程SST模型方程可以表示为
(7)
涡粘性系数可以由下式确定
(8)
湍流生成项Pk,Pω分别定义为
(9)
模型控制方程中的常数Φ={Cω,σk,σω,β}由下式求得:
Φ=F1Φ1+(1-F1)Φ2
(10)
以下为计算中选取的各拟合常数的数值,其中d为到物面的最小距离。
1.3 计算方法及模型
采用体积法对控制方程进行离散化,对流项的离散格式采用高分辨率格式,湍流数值方程的离散格式也采用高分辨率格式,物理时间步长为0.001s。当残差小于1×10-7时认为迭代收敛,停止计算。
模型采用直径1.6m的电动飞机固定桨距的两叶螺旋桨,采用四面体网格进行划分旋转的流场区域,在物面生成十分之一弦长高度的附面层。网格数量为530万,网格如图1所示。
图1 螺旋桨旋转域网格
静止的流场域采用六面体的结构网格进行模拟,网格数量为280万。两个流场域的交界面采用流体与流体耦合的形式进行计算。图2为边界条件的设置情况,分别为速度入口、压力出口、壁面。由于已知来流速度,将入口边界定义为速度入口条件,添加来流速度变量。螺旋桨表面定义为无滑移壁面,出口和四周定义为环境压力101 kPa,无压力梯度的边界条件。
图2 流场域及边界条件
在各计算状态收敛后,采用CFD-POST进行数据的后处理。分析转速为2380RPM时的流场特性。图3和图4桨叶分别为桨叶迎、背风表面压力云图和桨叶表面极限流线图。
图3 桨叶表面压力云图
从图3可以看出螺旋桨迎风面压力最小值出现在r=0.8~0.9 R的叶素前缘,压力总体上从桨叶尖部前缘向桨根后缘逐步增加。背风面压力在r=0.85~0.95 R桨叶后缘较大,整体变化趋势较为缓和。
图4 桨叶表面极限流线
从图4可以看出来流在经过螺旋桨表面平滑的流向后缘,整个螺旋桨上下表面的流线分布比较均匀,仅桨尖位置有很小范围的流线不是相互平行,可以考虑为桨叶尖部形成的桨尖涡诱导了流线的偏转,由于整个区域的流线都保持的很好,没有出现分离,证明高效率螺旋桨的设计方法是可行的。
为了验证本文数值模拟方法的准确性,根据等前进比相似准则[16],在西北工业大学NF-3风洞进行了风洞试验,模型直径为0.96 m,风洞试验段宽3.5m、高2.5 m、长12 m,截面为切角矩形,切角为0.6 m。湍流度为0.078%。试验使用西工大700型六分量盒式天平(编号TP0701)。六分量天平的电压信号采集由VXI数据采集系统完成。该系统有64个通道,采集速度为100 k/s,采集速度不小于100 kHz。试验对螺旋桨零风速情况下的拉力、力矩进行了测量。螺旋桨在风洞中的安装如图5所示。
图5 螺旋桨在风洞中
螺旋桨的流场数值计算比较复杂,影响因素较多,为了更好的对比实验结果,本文对拉力、力矩与实验数据进行了对比分析。如图6所示,随转速的增加,螺旋桨的拉力系数拉力的吻合的差距越小,两者的差距维持在2%以内。如图7所示,力矩随着转速的增加,差距逐步变大,两者的差距维持在10%以内。两个数据与试验基本一致,吻合度较好。
图6 计算拉力值与试验对比
图7 计算扭矩值与试验对比
本文通过基于雷诺平均Navier-Stokes方程的多重参考坐标系MRF方法对螺旋桨滑流进行准定常数值模拟计算与地面风洞试验,分析了某型三叶螺旋桨的气动特性,得出以下结论。
(1)总体来看,该准定常方法与试验的偏差在10%范围以内,且差距趋势基本保持不变,可以为工程提供较为准确的预测和评估。
(2)从流场的云图及流线等信息可以看出,该方法可以较为准确的评估流场的实际状况,对于螺旋桨的优化设计提供参考。
(3)旋转域采用非结构网格,静止域采用结构网格,这种操作可以大大节省网格的生成时间,对于精度的影响有待进一步的研究。
[1]项松,佟刚,吴江,等.某型三叶螺旋桨的设计及性能试验[J].航空动力学报,2016,30(8):1793-1798.
[2]HANSON D B.Compressible lifting surface theory for propeller performance calculation[J].Journal of Aircraft,1985,22(1):19-27.
[3]SCHULTEN J B H M.Advanced propeller performance calculation by a lifting surface method[J].Journal of Propulsion and Power,1996,12(3):477-485.
[4]LIESER J A,LOHMANN D,ROHARDT C H.Aeroacoustic design of a 6-bladed propeller[J].Aerospace Science and Technology,1997,1(6):381-389.
[5]ANGELO S D,BERARDI F,MINISCI E.Aerodynamic performances of propellers with parametric considerations on the optimal design[J].The Aeronautical Journal,2002,106(1060):313-320
[6]SLAVIK S.Preliminary determination of propeller aerodynamic characteristics for small aeroplanes[J].Acta Polytechnica,2004,44(2):103-108.
[7]SABZEHPARVAR M.In-flight thrust measurements of propeller-driven airplanes[J].Journal of Aircraft,2005,42(6):1543-1547.
[8]GUI O,ROSEN A.Propeller performance at low advance ratio[J].Journal of Aircraft,2005,42(2):435-441.
[9]WALD Q R.The aerodynamics of propellers[J].Progress in Aerospace Sciences,2006,42(2):85-128.
[10]ROMEO G,CESTINO E,PACINO M,et al.Design and testing of a propeller for a two-seater aircraft powered by fuel cells[J].Journal of Aerospace Engineering,2012,226(7):804-816.
[11]夏贞锋,罗淞,杨永.基于激励盘理论的螺旋桨滑流数值模拟研究[J].空气动力学学报,2012,30(2):219-225.
[12]刘远强,郭金锁,项松,等.基于片条理论的螺旋桨性能计算[J].沈阳航空航天大学学报,2013,30(1):43-46.
[13]MORGADO J,ABDOLLAHZADEH M,SILVESTRE MAR,et al.High altitude propeller design and analysis[J].Aerospace Science and Technology,2015,45(6):398-407.
[14]CHEN YAXI,LIU PEIQING,TANG ZHIHAO,et al.Wind tunnel tests of stratospheric airship counter rotating propellers[J].Theoretical and Applied Mechanics Letters,2015,5(1):58-61.
[15]项松,王吉,张利,等.一种高效率螺旋桨设计方法[J].航空动力学报,2015,30(1):136-141.
[9]王建明,刘炜,艾延廷,等.圆柱-平板角区湍流流动控制的风洞试验研究[J].实验流体力学,2011,25(5):10- 15.
[16]刘沛清,马蓉,段中喆,等.平流层飞艇螺旋桨地面风洞试验[J].航空动力学报,2011,26(8):1775-1781.
(责任编辑:吴萍 英文审校:赵欢)
Aerodynamic numerical and test study of a three-blade propeller
XIANG Song,LIU Yuan-qiang
(Liaoning Key Laboratory of General Aviation,Shenyang Aerospace University,Shenyang 110136,China)
Propeller slipstream and performance of a three-blade propeller were simulated using Multiple Reference Frames (MRF) quasi-steady method based on the Reynolds-averaged Naiver-Stokes(RANS)equation and Shear Stress Transport (SST) turbulence model.Thrust,torque and efficiency of the three-blade propellers caused by simulation were analyzed under different rotational speed compared to the results of tests.Deviation value of the tension is around 2%,and deviation value of the torque is about 10%.It is found by comparison that calculation results are in good agreement with those obtained from the test results,which can provide the guidance for propeller design applications.
three-blade propeller;aerodynamic performance;numerical simulation;wind tunnel test;general aviation
2016-12-29
辽宁省教育厅科研项目(项目编号:L201622);沈阳市科技计划项目(项目编号:F16-205-1-07);辽宁省科学技术计划(项目编号:2015020167);辽宁省百千万人才工程项目(项目编号:2014921048)
项松(1978-),男,辽宁灯塔人,高级工程师,主要研究方向:高效率螺旋桨设计,E-mail:xs74342@sina.com。
2095-1248(2017)03-0032-05
V211.44
A
10.3969/j.issn.2095-1248.2017.03.004