-
0 引言
-
民用飞机前缘装置主要作用于增升装置打开时,起到延迟大攻角时的分离,从而推后失速攻角,提高最大升力系数的作用。空客A350XWB-900飞机前缘增升装置为了满足起飞升阻比及二阶段爬升需求使用了内侧前缘下垂,外侧前缘缝翼的形式。该形式起飞构型阻力较低,着陆延迟失速提高升力系数的目的[1-3]。国内外对缝翼这类形式的前缘增升装置研究较多,但对前缘下垂(DND)这种较新型的前缘形式研究较少。Strüber等人[4-8]对A350XWB-900的前缘内外侧增升装置进行了介绍,将DND应用在机翼内侧,对气动特性和整体设计都能够带来一定的好处,内侧较大的弦长导致较低的当地升力,应用前缘下垂降低当地攻角可以进一步延迟失速,DND经典压力分布可见图1,头部吸力峰与后侧第二个吸力峰压力分布相近,与缝翼相比,DND形式在减阻降噪上都有优势。前缘下垂不仅利于飞机气动性能提升,而且更适合与飞机所装配的罗·罗挂机与发动机紧耦合发动机组合,内侧相对较长的弦长导致当地升力不足,前缘下偏可以有效降低当地攻角达到延迟失速的作用。但对影响前缘下垂的设计参数研究较少,本文从设计出发,针对前缘下垂典型二维多段翼进行数值模拟,对影响多段翼型气动性能的主要4个参数,前缘下垂弦长、前缘下垂偏度、转轴高度以及尾缘夹角开展了研究。
-
图1 DND翼型典型压力分布
-
1 数值方法及相关气动验证
-
1.1 控制方程
-
基于雷诺平均的Navier-Stokes(RANS)方程的流体求解方法因其在计算精度与计算资源消耗上达到一个较好的平衡而在工业领域广泛应用。RANS控制方程如公式(1)所示:
-
式中:U为解向量,F、G、H为通量,J为源项。基于FLUENT17.0求解。湍流模型S-A。
-
1.2 标模验证
-
数值方法采用标模30P30N加以验证。30P30N为1993年NASA Langley举行的CFD Challenge Workshop算例,其风洞试验主要是90年代在NASA Langley的低湍流度增压风洞中完成的[9]。弦长Cref =1.9m,雷诺数Re=9×109,速度Ma=0.175。网格使用混合网格,总网格量为25万。第一层网格厚度为10-6倍弦长,采用速度基+COUPLE耦合,对流项采用二阶迎风格式,扩散项为中心差分,计算收敛精度为10-6。计算结果如图2所示。从图中可以看出,升力与力矩对试验对标较好;阻力计算结果偏大。
-
图2 30P30N计算与试验力系数对比
-
2 数值模拟
-
2.1 模型及网格说明
-
模拟所使用的带DND的多段翼基于超临界翼型标模SC(2)-0410[10]设计,其中襟翼弦长为0.2c,扰流板弦长为0.126c。前缘下垂部件外形如图3所示。
-
图3 前缘下垂主要参数
-
各部件偏度及缝道参数如表1~表5所示。表1为各状态后缘的偏度及缝道参数;表2为前缘下垂弦长设置参数;表3为前缘下垂高度影响设置;表4为前缘下垂前缘偏度影响设置;表5为前缘下垂尾缘夹角影响设置。
-
气动模拟网格量22万,使用混合网格,如图4所示。
-
图4 前缘下垂网格示意
-
2.2 网格无关性验证
-
针对弦长11%c,前缘下垂偏度26°,转轴高度2mm,尾缘夹角8°的DND构型进行网格无关性验证,对比网格量为9万粗网格,22万中网格及30万细网格的三种网格策略,使用流动参数弦长Cref=1m,雷诺数Re=9×109,速度Ma=0.2,本文后续所有计算均使用了该工况设置。表6为三种粗细网格的升力系数差量,在线性段及近失速升力系数表现相近,最差量都在0.1%左右以内,可以认为在此网格量范围内对计算结果影响不大。结合计算量和精度综合考虑,选取22万中间粗细度网格作为本文中后续计算工作的基准网格。
-
3 前缘下垂影响分析
-
表7~表9及图5展示了四个参数对升力系数影响。由表可知:1)在线性段,以4°为例,弦长增加与偏度增大在线性升力均有负面影响,但转轴高度与尾缘夹角对线性升力几乎没有影响;2)在近失速段,以12°为例,弦长增加与尾缘夹角的提升,可明显提升升力;下垂偏度增加,以26°为界,小于该角度可提升升力,但大于该角度后,影响不再明显;随转轴高度下降,升力出现一定提升,但到2mm后反而有所下降;前缘下垂尾缘夹角增大,可提升近失速段升力;3)最大升力系数规律与近失速段规律相似。弦长对升力特性的影响最明显,弦长为13%c时最大升力系数最高,可达到3.134,而弦长为6%c时,最大升力系数最小,为3.008。
-
图5 多段翼升力规律
-
图6为表2中弦长11%c构型在AOA=12°下的速度云图,图中画出Ma >0.45的区域以表征加速区域。从图中可知前缘下垂构型在头部存在两个区域存在流动加速,第一个是前缘下垂头部,第二个区域为固定翼与主翼拱起区域。这两处加速区域会产生两个吸力峰,作为下文中研究前缘下垂参数机理的主要因素。
-
图6 前缘下垂速度云图
-
为进一步了解这一现象的机理,我们给出了前缘下垂上翼面的压力分布分析,着重研究参数对前缘下垂翼型头部两个吸力峰的影响。弦长参数增加影响与经典多段翼气动理论相符合,在本文中不再分析。图7为转轴高度对近失速段(AOA=12°)压力分布的影响,从图中可知,随着转轴高度增加,即转动半径减小,前缘下垂面积增量减小,在前缘下垂头部吸力有所提高,前缘下垂尾缘吸力变化不明显,固定翼前半部分升力下降,后半部分升力提高,但整体升力影响不明显。
-
图7 转轴高度对主翼头部吸力影响
-
图8为前缘偏度对近失速段攻角12°时压力分布的影响,从图中可知,随前缘下垂偏度增大,前缘下垂头部吸力下降明显,前缘下垂尾缘吸力下降,固定翼前半部分升力提升,后半部分升力提高。
-
图8 前缘偏度对主翼头部吸力影响(AOA=12°)
-
图9为尾缘夹角对近失速段(AOA=12°)压力分布的影响,从图中可知,随尾缘夹角增大,前缘下垂头部吸力变化不明显,但前缘下垂尾缘吸力下降,固定翼前半部分升力下降,后半部分升力显著提高。
-
图9 尾缘夹角对主翼头部吸力影响
-
4 结论
-
通过数值模拟研究了前缘下垂的四参数对典型增升构型二维多段翼升力影响,得出以下结论,可为前缘增升装置设计提供指导:
-
1) 在线性段,弦长增加与偏度增大在线性升力均有负面影响,但转轴高度与尾缘夹角对线性升力几乎没有影响;
-
2) 在近失速及失速段,弦长增加与尾缘夹角的提升,可明显提升升力;下垂偏度增加,以26°为界,小于该角度可提升升力,但大于该角度后,影响不再明显;随转轴高度下降,升力出现一定提升,但到2mm后反而有所下降;前缘下垂尾缘夹角增大,可提升近失速段升力;
-
3) 四参数对升力影响主要体现在头部两段吸力峰的消长,设计中需综合考虑吸力峰特征,并加以应用。
-
参考文献
-
[1] 《飞机设计手册》总编委会.飞机设计手册:第6册气动设计[M].北京:航空工业出版社,2002.
-
[2] DAM C P V.The aerodynamic design of multi-element high-lift systems for transport airplanes[J].Progress in Aerospace Sciences,2002,38(2):101-144.
-
[3] SMITH A M O.High-lift aerodynamics[J].Journal of Aircraft,1975,12(6):501-530.
-
[4] STRÜBER H.The aerodynamic design of the A350 XWB-900 high-lift system[C]//Proceedings of the 29th Congress of the International Council of the Aeronautical Sciences.St.Petersburg:[s.n.],2014.
-
[5] RECKZEH D.Multifunctional wing movables:design of the A350XWB and the way to future concepts[C]//Proceedings of the 29th Congress of the International Council of the Aeronautical Sciences.St.Petersburg:[s.n.],2014.
-
[6] SUTCLIFFE M,RECKZEH D,FISCHER M.Hicon aerodynamics-high-lift aerodynamic design for the future[C]//Congress of the International Council of the Aeronautical Sciences.Hamburg:ICAS,2006.
-
[7] RECKZEH D.Aerodynamic design of the high-lift wing for a megaliner aircraft[J].Aerospace Science and Technology,2003,7(2),107-119.
-
[8] SCHINDLER K,RECKZEH D,SCHOLZ U,et al.Aerodynamic design of high-lift devices for civil transport aircraft using RANS CFD[C]//28th AIAA Applied Aerodynamics Conference.Chicago:AIAA,2010.
-
[9] KLAUSMEYER S M,LIN J C.Comparative results from a CFD challenge over a 2D three-element high-lift airfoil[J].NASA Langley Technical Report Server,1997.
-
[10] HARRIS C D.NASA super-critical airfoils:a matrix of family-related airfoils:NASA Technical Paper 2969[R].USA:NASA,1990.
-
摘要
采用计算流体力学方法,针对基于典型二维多段翼型NACA0410设计的带前缘下垂构型的多段翼型进行数值模拟,研究了前缘下垂四种参数对多段翼升力特性的影响。在所研究的范围内,结果表明:1)在线性段,弦长增加与偏度增大对线性升力均有负面影响,但转轴高度与尾缘夹角对线性升力几乎没有影响;2)在近失速及失速段,弦长增加与尾缘夹角的提升,可明显提升升力;下垂偏度增加,以26°为界,小于该角度可提升升力,但大于该角度后,影响不再明显;随转轴高度下降,升力出现一定提升,但到2 mm后反而有所下降;前缘下垂尾缘夹角增大,可提升近失速段升力;3)四参数对升力影响主要体现在头部两段吸力峰的消长,设计中需综合考虑吸力峰特征,并加以应用。
Abstract
In this paper, the computational fluid dynamic method was applied to simulate the flow field around 2D multi-element airfoil with droop nose devices (DND) at leading edge. A typical 2D airfoil NACA0410 was used to design the multi-element airfoil. The effect of four factors used to design DND on aerodynamic characteristics of a multi-element airfoil was studied. The result shows that, 1) in the linear segment, increases in chord length and deflection angle both have negative effect on lift, but the angle between the hinge height and the trailing edge has little effect on lift; 2) close to /at stalling point, increases in chord length and trailing edge angle can obviously increase lift; lift increases when deflection angle is less than 26° and lift increases as hinge height decreases, but decreases when hinge height is less than 2 mm; 3) 4 parameters mainly affect the fluctuation of two suction peaks at the leading edge. These four parameters should be comprehensively considered and applied during design.
Keywords
droop nose devices(DND) ; lift ; civil aircraft ; aerodynamic design ; high lift devices design ; CFD