飞机着陆过程中提高气动性能的一种新方法

为了提高飞机在着陆过程中的气动性能,提出了一种新方法:将翼型上翼面的一段表面设计为活动部分。当飞机进入着陆阶段的较大迎角时,通过活动部分在上翼面形成一个台阶产生稳定的驻涡,再联合Gurney襟翼,达到同时提高翼型的升力、失速迎角及增加翼型阻力的目的。在NACA2415翼型上对上述方法进行了验证。结果表明,翼型最大升力系数从原始翼型的1.548 232提高到2.160 687,最大升力系数所对应的迎角可以从原始翼型的17°提高到20°。可见,所提出的新方法对提高飞机的着陆性能是有效的。     

涡襟翼在不同雷诺数下的控制分离特性研究

受鸟类抬起羽毛控制分离流的启发,涡襟翼成为翼型大迎角分离流的控制措施之一。采用数值模拟方法研究不同雷诺数下涡襟翼在控制翼型大迎角分离流动时的气动特性及其物理机制。结果表明：涡襟翼在低雷诺数下能够极大地改善翼型的大迎角升力特性,其物理机理是涡襟翼将翼型主分离涡的涡心位置控制在离翼型更近的区域,且涡心位置的涡量得到大幅提升,使得涡心附近的低压特性影响到翼型上表面,而且涡襟翼能够将翼型上方前区的低压与下游的高压隔开;但是在高雷诺数（对应常规飞机雷诺数）下涡襟翼改善翼型大迎角气动特性的效果远不如低雷诺数情况,由此解释了为什么鸟类能够通过羽毛抬起提高升力特性,而常规飞机的涡襟翼只能作为阻力板使用的原因。     

翼型等速上仰绕流结构的观测

 翼型在不同转速下等速上仰时,其绕流结构不一样,转速越高,前缘涡开始形成的迎角越大。给定转速时,初始迎角大前缘涡形成得早。终止迎角对流动结构的影响与转速有关。平板翼型上仰时,前缘涡在较小迎角开始形成,涡的尺寸较大。     

适用于多迎角结冰试验的混合翼型设计

为突破结冰风洞对翼型模型尺寸的限制,提出了一种新的混合翼型设计方法,可使用一套混合翼型模拟原始翼型在不同迎角下的结冰试验,弥补了以往混合翼型只能在单个设计迎角下使用的缺陷。方法采用多段翼的形式设计混合翼型,以多目标迎角等结冰试验条件作为设计输入,优化设计主翼外形和襟翼的位置、偏转角度,利用襟翼位置和偏转角度的变化实现混合翼型和原始翼型前缘表面的压力系数在不同迎角下能够保持一致,继而保证前缘结冰外形的一致性。设计的混合翼型较原始翼型弦长减小接近40%,在冰风洞中对混合翼型和原始翼型在目标结冰条件下进行试验,对比结果显示,在选取的多个目标迎角下,混合翼型和原始翼型二维截面的结冰外形基本一致,提出的混合翼型设计方法能够有效的模拟原始翼型在不同迎角下的结冰外形。     

变来流速度下旋翼翼型非定常气动特性分析

为分析变来流速度状态下的旋翼翼型气动特性,提出了利用翼型平移来模拟来流速度变化的数值方法.在此方法基础上,采用基于隐式LU-SGS(lower upper symmetric Gauss-Seidal)方法的非定常雷诺平均N-S(Navier-Stokes)(RANS)方程,模拟了SC1095旋翼翼型在定迎角 变来流速度及变迎角 变来流速度状态下的非定常气动特性.通过对比分析发现:翼型在变速度-定迎角状态下会表现出明显的非定常现象,产生了前缘分离涡,气动特性会出现明显的迟滞效应及波动现象,脉动速度越大,非定常效果越明显.并且基准速度越大,翼型气动特性的峰值越大；翼型迎角越大,非定常涡出现的也越早.考虑直升机旋翼翼型实际工作环境,在变速度-动态失速状态下,翼型最大迎角处的气动力会得到一定程度的削弱,在小迎角下的气动力得到一定程度的增强,且脉动速度越大,翼型的非定常特性也越强.      

大迎角水平振动翼型绕流的数值模拟

从涡量流函数形式的N-S方程出发,数值计算了大迎角水平振动翼型的粘性绕流问题,得到近场涡结构及非定常演化过程。探讨了振动频率、振幅及迎角对流场涡结构及翼型升、阻特性的影响。结果表明,振动频率及振幅增加有助于提高翼型的平均升力;过大的迎角不利于升力提高。     

Mini-TED跨声速气动特性及其增升机理

通过求解二维可压Navier-Stokes方程,研究了NACA0012翼型加装微型后缘增升装置（mini-TED）后的跨声速流场特性,与Gurney flap （GF）对比分析了几何参数对mini-TED后方涡系及翼型气动特性的影响.将mini-TED的几何细节参数定义为弦向长度和有效高度,两者方向正交.在相同迎角下仅改变mini-TED的弦向长度,后缘涡系结构虽发生变化,但翼型气动力几乎没有影响;反之仅改变有效高度则后缘涡系和翼型气动力系数同时发生明显改变,且与同等高度下的GF气动系数相近.结果表明：有效高度是影响翼型气动特性的决定因素.有效高度改变了mini-TED后涡系的发生范围,而相对于整个翼型绕流,后缘涡系的大小是影响翼型流场最重要的因素,而涡系的微观结构和形态的改变影响相对很小.加装mini-TED后上表面激波位置后移、下表面激波强度削弱,从而翼型表面压力分布特性发生了改变.随有效高度增大,mini-TED诱导的涡系发生区域随之增大,引流作用增强,翼型升力系数、阻力系数和低头力矩系数提高,同时相同迎角下翼型的升阻比明显提高.     

NACA0012翼型低雷诺数绕流的实验研究

通过水槽氧气泡流动显示和PIV测速实验研究了NACA0012翼型在雷诺数为8200时的流动特性,重点炎注了翼型绕流结构随迎角的变化。研究发脱：分离点和分离翦切层形成旋涡的位置随迎角的增大而向上游移动,同时翼型上表面流动分离后形成的回流区尺寸随着翼利迎角的增加而增大。当流动再附于翼型上表面时,在再附点附近能够观测到展向涡的三维演化过程,并能观测到展向涡的局部配对现象。     

高速自然层流翼型的设计与风洞实验研究(英文)

为了发展商务客机用的低阻力翼型,本文在荷兰科技大学TST27风洞开展了高速自然层流翼型的实验研究工作。为了减小测量技术对模型流场的干扰,提高测量精度,分别采用了IR,PIV以及阴影法等非接触测量技术。通过这些非接触测量技术的使用,对翼型的转捩位置、后缘分离以及升阻力等气动特性进行了测量研究。从实验结果可以看出,在小迎角范围内,翼型可以在其表面保持较长的层流段,并且后缘分离涡不会随迎角增大而向前发展,使翼型具有高升力、低阻力的特性。从实验结果和计算结果的差别可以看出,实验数据是可靠的,设计的翼型基本满足了设计要求。     

大厚度翼型分离流场的稳定性探讨

基于k-ω的SST两方程湍流模型,求解了雷诺平均Navier-Stokes方程获得定常和非定常气动力,耦合翼型弹性振动方程,在时间域内模拟了弹性支撑下的大厚度对称翼型流场及其翼型动态过程,研究结果发现：在所给定的雷诺数条件下,弹性支撑下的大厚度翼型在一定大迎角范围内存在分离流场稳态的转化现象,此结果的意义在于：1）同一基本迎角下的翼型分离涡存在多个稳态流场;2）弹性支撑的翼型在大迎角分离流条件下气动力在物理上本身就存在非唯一性,从一个新的角度解释了大迎角气动力测量中的分散性。     

翼型极大迎角风洞试验技术研究

在NF-3风洞的二元试验段开展了翼型极大迎角（±180°）条件下气动特性的试验技术研究。针对翼型极大迎角风洞试验的洞壁干扰,提出了风洞壁压信息洞壁干扰修正的改进方法。试验结果表明,发展的试验技术和提出的洞壁干扰修正方法适合于翼型极大迎角试验。     

基于高斯过程回归和遗传算法的翼型优化设计

针对高升阻比风力机翼型前缘曲率半径较大的问题,传统的翼型参数化方法前缘控制能力不足,且基于面元法XFOIL预测精度差的问题,利用增强类函数/形函数转换（CST）参数化方法控制翼型的形状变化、拉丁超立方实验设计、计算流体力学（CFD）流场计算模块、高斯过程回归模型和遗传算法,提出了基于高可信度Reynolds average Navier-Stocks（RANS）和高斯回归模型辅助遗传算法的翼型优化设计方法。结果表明:基于高斯回归模型的翼型优化方法,可以将优化所用CFD计算次数降低一阶,从而大幅度提升优化设计效率。由标准算例超临界翼型RAE2822的降阻设计表明,在百次量级的CFD次数阻力降低43.16%,激波被削弱且升力、力矩和面积严格满足约束。由风力机翼型NACA64618的最大化升阻比优化设计表明,所设计翼型不仅在设计攻角和副设计攻角处升阻比大大增加,在整个小攻角范围内其气动性能都得到了提升,且两个主设计点,无不良阻力的产生。      

模拟厚度翼型非定常运动的涡方法

本文研究二维厚度翼型的大攻角非定常运动。用沿翼型表面分布的源和面涡来模拟翼型的厚度效应和升力效应,应用离散涡方法模拟前缘分离涡层和尾涡层。计算了大攻角翼型作俯仰运动的情形,与实验结果符合较好。     

应用协同射流原理的旋翼翼型增升减阻试验研究

协同射流技术作为一种新型主动流动控制技术,是突破旋翼翼型高增升减阻设计的最有潜力的发展方向之一。以 OA312 旋翼翼型作为基准翼型,研制微型涵道风扇组为驱动的旋翼翼型 CFJ 风洞测力模型,开展基于前缘高负压零质量内循环协同射流原理的旋翼翼型高增升减阻低速风洞试验,研究吹气口大小、吸气口大小和上翼面下沉量等基础参数对增升减阻的影响规律,探讨 CFJ 旋翼翼型关键参数最佳取值。结果表明：与OA312 基准翼型相比,小攻角状态时,CFJ 旋翼翼型可显著降低阻力系数,甚至出现“负阻力”现象,实现了零升俯仰力矩基本不变;大攻角状态时,CFJ 旋翼翼型可显著提升最大升力系数和失速迎角,其中,最大升力系数可提升约 67.5%,失速迎角推迟了近 14.8°。     

Wake structure and similar behavior of wake profiles downstream of a plunging airfoil

Very limited attention has already been paid to the velocity behavior in the wake region in unsteady aerodynamic problems.A series of tests has been performed on a flapping airfoil in a subsonic wind tunnel to study the wake structure for different sets of mean angle of attack,plunging amplitude and reduced frequency.In this study,the velocity profiles in the wake for various oscillation parameters have been measured using a wide shoulder rake,especially designed for the present experiments.The airfoil under consideration was a critical section of a 660 kW wind turbine.The results show that for a flapping airfoil the wake structure can be of drag producing type,thrust producing or neutral,depending on the mean angle of attack,oscillation amplitude and reduced frequency.In a thrust producing wake,a high-momentum high-velocity jet flow is formed in the core region of the wake instead of the conventional low-momentum flow.As a result,the drag force normally experienced by the body due to the momentum deficit would be replaced by a thrust force.According to the results,the momentum loss in the wake decreases as the reduced frequency increases.The thrust producing wake pattern for the flapping airfoil has been observed for suffi ciently low angles of attack in the absence of the viscous effects.This phenomenon has also been observed for either high oscillation amplitudes or high reduced frequencies.According to the results,for different reduced frequencies and plunging amplitudes,such that the product of them be a constant,the velocity profiles exhibit similar behavior and coalesce on each other.This simi larity parameter works excellently at small angles of attack.However,at near stall boundaries,the similarity is not as evident as before.     

Numerical study of separation on the trailing edge of a symmetrical airfoil at a low Reynolds number

This study focuses on the trailing-edge separation of a symmetrical airfoil at a low Rey-nolds number. Finite volume method is adopted to solve the unsteady Reynolds-averaged Navier-Stokes (RANS) equation. Flow of the symmetrical airfoil SD8020 at a low Reynolds number has been simulated. Laminar separation bubble in the flow field of the airfoil is observed and process of unsteady bubble burst and vortex shedding from airfoil surfaces is investigated. The time-dependent lift coefficient is characteristic of periodic fluctuations and the lift curve varies nonlinearly with the attack of angle. Laminar separation occurs on both surfaces of airfoil at small angles of attack. With the increase of angle of attack, laminar separation occurs and then reattaches near the trailing edge on the upper surface of airfoil, which forms laminar separation bubble. When the attack of angle reaches certain value, the laminar separation bubble is unstable and produces two kinds of large scale vortex, i.e. primary vortex and secondary vortex. The periodic processes that include secondary vortex production, motion of secondary vortex and vortex shedding cause fluctuation of the lift coefficient. The periodic time varies with attack of angle. The secondary vortex is relatively stronger than the primary vortex, which means its influence is relatively stronger than the primary vortex.     

翼型大攻角绕流的数值模拟

以求解二维N-S方程数值模拟NACA0012翼型大攻角状态的可压绕流特性;N-S方程是在贴体坐标系中给出的,以代数方法生成C型网格系统。采用LU-ADI格式计算,为提高格式的稳定性在隐式和显式部分分别添加了2阶和4阶人工粘性项。应用BaldwinLomax湍流二层代数模型模拟了大攻角时前缘分离涡的形成,旋涡对流及其非定常现象。在某些Mach数和攻角下NACA0012翼型的湍流解具有周期性。通过比较,本文数值计算结果同实验及国外相应的数值计算结果基本吻合。     

Experiments on unsteady vortex flowfield of typical rotor airfoils under dynamic stall conditions

A new experiment for airfoil dynamic stall is conducted by employing the advanced particle image velocimetry(PIV) technology in an open-return wind tunnel. The aim of this experimental investigation is to demonstrate the influences of different motion parameters on the convection velocity, position and strength of leading edge vortex(LEV) of airfoil under different dynamic stall conditions. Two different typical rotor airfoils, OA209 and SC1095, are measured at different free stream velocities, oscillation frequencies, and angles of attack. It is demonstrated by the measured data that the airfoil with larger leading edge radius could notably decrease the strength of LEV. The angle of attack(Ao A) of airfoil can obviously influence the dynamic stall characteristics of airfoil,and the LEV would be effectively inhibited by decreasing the mean pitch angle. In addition, the convection velocity of LEV is estimated in this measurement, and the results demonstrate that the influence of airfoil shape on convection velocity of LEV is limited, but the convection velocity of LEV would be increased by enlarging the oscillation frequency. Meanwhile, the convection velocity of LEV is a time variant value, and this value would increase as the LEV convects to the trailing edge of airfoil.     

扇翼飞行器翼型附面层控制数值模拟

基于扇翼飞行器翼型特殊的几何形状及流场特性,在原有翼型的弧形槽下方和后缘加装控制阀门,通过调节阀门开启及开启尺寸的大小,利用弧形槽低压涡所产生的吸力对翼型后缘的附面层进行一定的控制,达到增升减阻的效果。通过采用计算流体力学的方法对其机理及阀门开启尺寸的影响进行了详细计算和分析,研究表明当阀门开启的尺寸为10 mm时,修改翼型的最大升力系数、失速迎角及相同迎角下的升力系数和推力系数均大于基本翼型;随着阀门开启尺寸的增大,修改翼型的最大升力系数和失速迎角均减小,但是在失速前,修改翼型在相同迎角下的升力系数大于基本翼型。此方法可以改变先前通过增大横流风扇的转速来提高其气动性能的做法,减小了能量的消耗,增大了整个飞行器的航程,为扇翼飞行器能够早日投入实际运用奠定了一定的理论基础。