不同襟翼偏角梯形翼构型气动特性数值模拟

基于雷诺平均Navier-Stokes(RANS)方程和结构网格技术,采用二阶离散精度的单调迎风格式(MUSCL),结合剪切应力输运(SST)两方程湍流模型和γ-Reθ转捩模型,研究了梯形翼高升力构型襟翼偏角变化对气动特性的影响。主要目的是进一步确认Trisonic Platform(TRIP)软件模拟高升力梯形翼不同襟翼偏角引起微小气动特性变化的能力。首先,简要介绍了采用的计算方法;其次,介绍了两种襟翼偏角的梯形翼模型及风洞试验;最后,在网格收敛性研究的基础上,采用全湍流和转捩两种方式模拟了梯形翼构型不同襟翼偏角对气动特性的影响。与试验数据的对比结果表明,采用全湍流和转捩两种方式均可以较好地模拟不同襟翼偏角对气动特性的影响量,采用γ-Reθ转捩模拟方式可以提高梯形翼构型气动特性的模拟精度,失速迎角附近的气动特性模拟需要进一步研究。     

湍流模型对梯形翼高升力构型的影响

采用"亚跨超CFD软件平台"(TRIP)数值模拟了梯形翼全展长高升力构型,主要目的是考察湍流模型对高升力构型气动特性的影响。相应的风洞试验是1998年在NASA Ames 12英尺增压风洞中完成的。本文采用一方程和两方程湍流模型,数值模拟了全展长高升力构型的气动特性并给出了典型站位的压力分布。研究表明,与修正后的试验数据相比较,在定常可收敛的迎角范围内,两种湍流模型得到的气动力系数和压力分布与试验结果吻合较好,SA一方程湍流模型的计算结果更接近试验值;迎角大于19.19°后,襟翼后缘的较大范围分离是导致采用SST模型不能得到收敛的气动特性的主要原因。     

梯形翼高升力构型的数值模拟技术

基于雷诺平均Navier-Stokes(RANS)方程和结构网格技术,采用亚跨超声速平台(TRIP3.0),数值模拟了美国国家航空航天局(NASA)梯形翼构型。研究了控制方程、网格密度、流动转捩和初始条件等不同影响因素对气动特性的影响。风洞试验是2002年在NASA Langley 14ft×22ft亚声速风洞中完成的,试验结果包括了基本气动力和力矩、表面压力系数和边界层速度型分布。计算结果与试验数据的比较表明:求解完全的RANS方程,提高了翼梢涡的模拟精度;网格密度主要影响翼梢涡的强度;转捩模型提高了边界层的模拟精度,进而提高了升力系数、俯仰力矩系数的模拟精度;最大升力系数及失速迎角对初始条件具有依赖性。     

转捩对风力机翼型和叶片失速特性影响的数值模拟

采用基于k-ωSST湍流模型的Gamma-Theta转捩模型对S809翼型和NREL Phase VI叶片进行了考虑转捩的气动力数值模拟,研究了转捩对其气动特性特别是失速特性的影响。首先对S809翼型在迎角0°～30°范围内开展了数值模拟,比较了转捩模拟和全湍流模拟获得的翼型升力系数和流场特征,发现前缘层流分离泡的存在显著影响了翼型的失速特性;然后对NREL Phase VI叶片开展了类似的数值模拟,结果显示转捩对叶片失速特性和翼型失速特性产生影响的作用方式是相似的。经过对以上数值模拟结果加以分析认为,转捩对翼型和叶片失速特性的影响主要通过前缘层流分离泡的作用体现出来,前缘层流分离泡的存在使翼型和叶片更早地进入深失速。     

半展长襟翼梯形翼构型数值模拟技术研究

基于雷诺平均的 Navier-Stokes 方程和拼接结构网格技术,采用 MUSCL 格式和 SST 湍流模型,研究了网格密度对半展长襟翼梯形翼高升力构型的数值模拟结果的影响。相应的风洞试验是1998年在 NASA Ames 12英尺增压风洞(PWT)中完成的,试验结果包括了总体气动特性、压力分布。研究内容主要包括网格密度对收敛历程、气动力特性、压力分布和表面流线的影响,以及气动力特性随迎角的变化。研究表明,Ma=0.15,α=16.7°时,网格密度对收敛历程、典型站位压力分布和表面流态基本没有影响,气动力特性随网格密度单调变化;采用不同密度的网格,典型剖面的压力分布与试验结果吻合良好;与修正后的试验数据相比较,数值模拟得到的失速迎角前的气动力系数与试验结果吻合良好。     

基于eN转捩预测方法的增升装置失速特性数值模拟研究

基于三维可压缩Navier-Stokes方程耦合e~N转捩预测方法,对30P30N二维多段翼型和DLR-F11三维增升装置进行了数值模拟研究。计算结果表明,e~N方法计算得到的转捩位置与试验值符合较好,增升装置翼面特别是主翼和襟翼下翼面存在大范围层流流动,为提高增升装置绕流计算准度,计算中考虑转捩的影响是非常有必要的。相比全湍计算,e~N转捩计算升力系数与试验值符合更好,提高了增升装置大迎角失速特性计算准度。     

CT-1标模大迎角静态气动特性数值模拟

 采用雷诺平均的Navier-Stokes方程,对接网格技术和TRIP（TRI sonic Platform）软件,数值模拟了CT-1标模大迎角静态气动特性。来流马赫数0.5,迎角0°～100°,侧滑角0°和5°,在与1.2 m跨声速风洞试验结果对比的基础上,重点讨论了对流项离散方法、湍流模型、远场距离和不同尾部处理形式对CT-1标模大迎角静态气动特性影响。研究结果表明,不同尾部处理形式对大迎角静态气动特性数值模拟结果影响显著,模拟试验的尾部支杆是对比大迎角静态测力试验结果与计算结果的基本要求。     

采用TRIP2.0软件计算DLR-F6构型的阻力

采用"亚跨超CFD软件平台"(TRIP2.0)数值模拟了DLR-F6构型,主要目的是通过计算DLR-F6构型的安装阻力考察TRIP2.0软件的数值模拟精度,并为运输机构型的气动特性计算积累经验.本文数值模拟采用的多块对接网格,测压和测力的试验结果均来自AIAA CFD Drag Prediction Workshop II(DPWII),对比计算结果采用了CFL3D的结果.本文详细研究了网格密度、湍流模型对DLR-F6翼身组合体和翼/身/架/舱复杂组合体两种构型的的总体气动特性和压力分布的影响,计算结果与相应的试验结果取得了较好的一致.本文采用SST两方程模型计算两种构型均得到了网格收敛结果;不同的湍流模型对压差阻力影响较小,对摩擦阻力影响较大;不同的网格密度和湍流模型对压力分布影响较小.     

30P-30N多段翼型复杂流场数值模拟技术研究

采用"亚跨超CFD软件平台"(TRIP2.0)数值模拟了30P-30N多段翼型的复杂流场,主要目的是考核湍流模型、转捩位置对多段翼型压力分布和典型站位速度型的影响。本文通过求解任意坐标系下的雷诺平均的N-S方程,采用多块对接结构网格技术,在与相应试验结果对比的基础上,详细研究了SA一方程湍流模型、SST两方程湍流模型、不同的转捩位置对该翼型压力分布和典型站位速度型的影响。本文的研究结果表明,采用全湍流模拟方式可以较好地模拟该多段翼型的压力分布,但对速度型的模拟精度较差;模拟试验的转捩位置可以改善主翼附面层与前缘缝翼边界层尾迹区的模拟精度;采用微吸气技术推迟前缘缝翼的转捩位置,可以进一步提高缝翼尾迹区的数值模拟精度。     

战斗机大迎角/过失速机动下的进气道气动特性

为掌握战斗机在大迎角和过失速机动飞行时进气道的稳、动态气动特性,采用基于动态嵌套网格的非定常雷诺平均Navier-Stokes （URANS）方程和大迎角风洞试验方法对某战斗机进行了研究,并通过大迎角和过失速机动飞行试验进行了验证。结果表明：大迎角稳态下进气道气动性能随迎角增大逐渐降低,天地相关性吻合良好,而计算仿真和飞行试验均捕捉了眼镜蛇机动下进气道的非定常迟滞效应。通过研究获得了战斗机在大迎角和过失速机动下的进气道气动特性,建立了过失速机动下进气道非定常非线性特性问题的研究方法。     

大飞机缝翼滑轨影响研究

通过数值模拟和风洞实验两种手段研究了大飞机缝翼滑轨对飞机气动性能的影响.分析了缝翼滑轨对缝道和机翼表面流动分布的影响,获得了缝翼滑轨参数对飞机气动性能的影响规律.数值模拟结果表明:缝翼滑轨对缝道内的流动形成了阻塞,改变了机翼表面的流动形态,减小了机翼附面层流动速度,降低了飞机的失速性能.实验结果表明:通过减小滑轨宽度、减少滑轨数量、采用圆形截面滑轨和滑轨外弯等能够有效降低滑轨影响,改善飞机失速性能;滑轨参数对小尺度模型实验结果的影响尤为显著.研究结果为3m量级和8m量级风洞缝翼滑轨模型设计提供了参考.     

Delaying stall of morphing wing by periodic trailing-edge deflection

Morphing wings can improve aircraft performance during different flight phases. Recently research has focused on steady aerodynamic characteristics of the morphing wing with a flexible trailing-edge, and the unsteady aerodynamic and stall characteristics in the deflection process of the morphing wing are worthy further investigation. The effects of the angle of attack and deflection rate on aerodynamic characteristics were examined, and based on the aerodynamic characteristics of the morphing wing, a method was developed to delay stall by using the flexible periodic trailing-edge deflection. The numerical results show that the lift coefficients in the deflection process are smaller than those in the static situation at small angles of attack, and that the higher the deflection rate is, the smaller the lift coefficients will be. On the contrary, at large angles of attack, the lift coefficients are higher than those in the static case, and they become larger with the increase of the deflection rate. Further, the periodic deflection of the flexible trailing-edge with a small deflection amplitude and high deflection rate can increase lift coefficients at the critical stall angle.     

飞机增升装置的数值模拟研究

增升装置的设计是飞机设计的一项极其重要的内容.本文通过多块网格技术生成计算网格,运用全湍流SST(Shear stress transport)模型和带有转捩修正的SST湍流模型对飞机三维增升装置的气动力特性和流场特征进行了数值模拟研究,通过对数值计算结果与实验结果的对比分析,研究了两种湍流模型的特点及其在飞机增升装置...     

带有横喷控制的导弹流场数值模拟

从包含多种组分的N-S方程出发，考虑两方程湍流模型，采用NND2M差分格式，对带有横向喷流的双锥旋称体高超声速绕流场进行数值模拟，计算结果与已有的试验数据进行对比，符合较好。在此基础上对带有横向喷流控制系统、型尾翼布局、高超声速飞行的导弹外流场进行数值模拟，研究了迎角、多个喷口、热喷流效应和湍流模型对气动力特性的影响；计算表明在无尾翼情况下有／喷流的气动力差别较小，喷流影响随迎角变化不敏感；对带有尾翼的气动布局，喷口位于背风区时喷流影响较小，喷口位于迎风面时气动力变化较大，压心明显前移；多喷口产生的附加推力和力矩不等于每个单喷口线性相加；湍流模型和热喷流效应引起流场结构改变，但是对总的气动力影响不大。     

Effects of wing flexibility on aerodynamic performance of an aircraft model

The low-speed wind tunnel experiment is carried out on a simplified aircraft model to explore the influence of wing flexibility on the aircraft aerodynamic performance. The investigation involves the measurements of force, membrane deformation and velocity field at Reynolds number of 5.4 × 10⁴–1.1 × 10⁵. In the lift curves, two peaks are observed. The first peak, corresponding to the stall, is sensitive to the wing flexibility much more than the second peak, which nearly keeps constant. For the optimal case, in comparison with the rigid wing model, the delayed stall of nearly 5° is achieved, and the relative lift increment is about 90%. It is revealed that the lift enhanced region corresponds to the larger deformation and stronger vibration, which leads to stronger flow mixing near the flexible wing surface. Thereby, the leading-edge separation is suppressed, and the aerodynamic performance is improved significantly. Furthermore, the effects of sweep angle and Reynolds number on the aerodynamic characteristics of flexible wing are also presented.     

CRM-WB风洞模型高阶精度数值模拟

通过求解雷诺平均Navier-Stokes （RANS）方程,采用5阶空间离散精度的WCNS和多块对接结构网格技术,开展了CRM翼身组合体风洞试验模型的高阶精度数值模拟,计算构型和计算状态来自AIAA第6届阻力预测研讨会（DPW Ⅵ）。主要目的是采用高阶精度方法,评估静气动弹性变形和模型支撑装置对CRM翼身组合体数值模拟结果的影响。通过与刚性外形的计算结果和NASA NTF风洞试验结果的对比,高阶精度数值模拟结果表明：迎角为3.0°时,静气动弹性变形使得机翼上表面激波位置前移并显著降低了外侧机翼上表面激波位置前的负压;迎角为3.0°时,模型支撑装置使得机翼上表面激波位置进一步前移,并导致翼梢处机翼上表面流动结构发生变化;迎角为4.0°时,计算模型中没有包含模型支撑装置是导致升力系数下降的主要原因;计算模型中包含机翼静气动弹性变形和模型支撑装置的数值模拟结果更加接近试验结果。     

鸭翼-前掠翼气动布局纵向气动特性实验研究

前掠翼布局由于其潜在的优势,在未来战斗机的研制中将占有日益重要的地位.本实验通过可变前掠翼和鸭式前翼布局的风洞测力实验,重点分析比较了平板机翼在不同掠角下的纵向气动性能以及鸭翼的影响.实验结果表明,前掠翼在大迎角时能有效提高模型的升力系数,小迎角时其升阻比也略优于后掠翼.前掠翼布局能有效推迟失速,具有良好的失速特性;前掠角较大时,升力系数曲线在失速迎角附近有一个升力系数的"平台",该布局具有"缓失速"特性.距离主机翼较远的鸭式前翼(模型M2)在主机翼前掠和后掠情况下,均可改善整体布局的失速特性,增大失速迎角,增强前掠翼布局缓失速的特点.近距耦合鸭翼(模型M3)显著提高了模型在大迎角下的升力系数.另外,主翼前掠和鸭式前翼布局飞行器具有较好的机动性.     

一种翼身融合飞行器的失速特性研究

翼身融合（BWB）布局飞行器作为下一代商用飞机的主要构型之一，越来越受到重视。对于翼身融合飞行器的研究主要针对其巡航状态的特性，而对其失速特性的研究较少。对一种翼身融合客机构型进行风洞试验研究，采用测力试验方法对其无增升装置的构型，以及具有翼梢小翼、前缘缝翼和机身上部双吊舱的组合部件构型下的纵向特性进行研究，特别是对其失速特性的分析，并通过二维粒子图像测试技术以及油流试验对其失速过程的流动机理进行研究。结果表明，无增升装置的基本构型下，翼身融合飞行器可以保持低速飞行，而各组合构型都具有提高最大升力系数的作用。对失速过程的分析表明，随着迎角的增大，飞机表面流场分离区域从翼梢开始逐渐向翼根以及机身发展，当外翼段完全处于分离区域时，飞机并不会马上失速，因为中心体同样具有提供升力的作用，且中心体的流动分离较外翼的流动分离更晚，所以当外翼在失速迎角出现升力损失时可以通过中心体的升力进行补偿，维持其低速飞行状态，真正的失速发生在中心体出现流动分离之后。     

可变正弦前缘对直机翼气动性能影响的研究

通过风洞测力实验,研究了可变正弦前缘对直机翼气动性能的影响。可变正弦前缘为平板式,安装于NACA0015直机翼的前驻点处,并能自由伸展与缩回。风洞实验结果表明：加装可变正弦前缘后,直机翼呈现出缓失速特性,失速后的升力系数均有不同程度的提高,大迎角下的气动性能得到改善。可变正弦前缘伸展长度对气动性能的影响研究表明,其伸展长度是影响机翼气动性能的重要参数,只有伸展长度较小时,才能获得较好的气动性能。     

Turbulent boundary layer separation control using plasma actuator at Reynolds number 2000000

An experimental investigation was conducted to evaluate the effect of symmetrical plasma actuators on turbulent boundary layer separation control at high Reynolds number. Compared with the traditional control method of plasma actuator, the whole test model was made of aluminum and acted as a covered electrode of the symmetrical plasma actuator. The experimental study of plasma actuators' effect on surrounding air, a canonical zero-pressure gradient turbulent boundary, was carried out using particle image velocimetry(PIV) and laser Doppler velocimetry(LDV) in the 0.75 m × 0.75 m low speed wind tunnel to reveal the symmetrical plasma actuator characterization in an external flow. A half model of wing-body configuration was experimentally investigated in the  3.2 m low speed wind tunnel with a six-component strain gauge balance and PIV. The results show that the turbulent boundary layer separation of wing can be obviously suppressed and the maximum lift coefficient is improved at high Reynolds number with the symmetrical plasma actuator. It turns out that the maximum lift coefficient increased by approximately 8.98% and the stall angle of attack was delayed by approximately 2° at Reynolds number2 ×10~6. The effective mechanism for the turbulent separation control by the symmetrical plasma actuators is to induce the vortex near the wing surface which could create the relatively largescale disturbance and promote momentum mixing between low speed flow and main flow regions.