相对厚度对低雷诺数流动中翼型动态气动力特性的影响

以固定翼微型飞行器为研究背景,研究了相对厚度对低雷诺数流动中翼型动态气动力特性的影响规律.采用Roe迎风差分格式和双时间步迭代方法,数值求解拟压缩性修正不可压Navier-Stokes方程,给出了数值算法与实验数据的对比验证.以翼型弦长为特征长度,在Re=500～50,000情况下,选取不同最大相对厚度和不同最大相对厚度位置的翼型,计算了其等速上仰时的动态气动力,结果表明前者对气动力影响显著,较小最大相对厚度值可获得较大的动态升阻比.     

零质量射流对翼型增升影响的数值模拟研究

针对安装了激振器的NACA0012二维翼型,采用数值模拟的方式对零质量射流技术进行研究.流场解算器采用了中国航空工业空气动力研究院计算所自主研发的ENSMB V1.2计算平台.计算结果给出了湍流模型对零质量射流数值模拟技术的影响,并且分析了激振器参数对翼型气动力特性的影响规律.     

零质量射流对翼型增升影响的数值模拟研究（下）

针对安装了激振器的NACA0012二维翼型,采用数值模拟的方式对零质量射流技术进行研究。流场解算器采用了我院计算所自主研发的ENSMBV1．2计算平台：计算结果给出了湍流模型对零质量射流数值模拟技术的影响,并且分析了激振器参数对翼型气动力特性的影响规律。     

旋翼翼型非定常动态失速特性的CFD模拟及参数分析

构建了一套基于运动嵌套网格技术和可压缩RANS方程的旋翼翼型非定常流动特性模拟的高效、高精度的CFD方法。首先,发展了基于Poisson方程求解的围绕翼型的粘性贴体正交网格生成方法,并提出了基于最小距离法(MDM)改进策略的运动嵌套网格生成方法,克服了弹簧法可能导致网格畸变的不足;其次,为准确模拟由湍流分离和气流再附引起的气动力的迟滞效应,基于RANS方程、双时间方法和高阶插值格式,建立了旋翼翼型非定常气动特性分析的高精度数值方法,并采用能够较好捕捉气流分离现象的S-A湍流模型;再次,针对旋翼后行桨叶动态失速时桨叶剖面来流速度较低、迎角较大的特点,为解决低来流速度时L-B半经验模型在旋翼翼型非定常动态失速计算中的局限性,并克服可压缩方程对低速流场计算收敛困难和精度低的问题,建立了基于Pletcher-Chen低速预处理方法、FAS多重网格法和隐式LU-SGS方法相结合的高效数值方法。应用发展的方法,分别针对NACA0012、SC1095旋翼翼型静态和轻度、深度动态失速进行计算,精确捕捉了气动力迟滞效应以及翼型前缘脱体涡的产生、对流和脱落过程,验证了本文方法的有效性;最后,着重针对NACA0012动态失速状态,开展了振荡参数对旋翼翼型非定常动态失速特性影响的分析,研究结果表明翼型迎角平均值、振幅及减缩频率的变化均能引起迟滞效应的改变并使得气动力峰值发生有规律的前、后移现象等。     

绕翼型分离流结构的数值研究

利用大涡模拟技术,对不同来流迎角下的NACA0012翼型绕流结构进行了数值模拟,详细地给出了翼型绕流分离流结构随迎角的变化特征和翼型在分离绕流中的气动力参数。在数值模拟中,采用了弱压缩流的控制方程,用贴体坐标技术进行了网格生成。     

基于非结构动态网格的翼型积冰过程数值模拟

讨论了积冰气象条件下翼型积冰过程的非结构动态网格数值模拟方法.外部无黏绕流流场采用基于非结构网格的分布式并行计算技术,空间离散采用VanLeer迎风有限体积格式,时间推进为5步显式龙格-库塔方法.采用四阶龙格-库塔方法求解过冷水滴运动轨迹方程,使用了一种合理、简单的物面水积冰量计算方法,提高了计算效率.使用非结构动态网格技术对积冰后的翼型网格进行快速调整.在对NACA0012翼型积冰形状进行校核验证后,对某型飞机在积冰气象条件下飞行时翼型积冰过程进行了数值模拟,得到了与参考文献结果相符的冰形特性.     

S1223翼型俯仰-沉浮运动的非定常气动特性分析

扑翼飞行器是未来航空领域的重要发展方向之一,而鸟类、昆虫等自然界的飞行生物所具有的出色飞行能力,为人造扑翼飞行器的设计工作提供了很好的参照。本文以鸟类的扑翼飞行为研究背景,针对简化的二维S1223翼型的刚性俯仰-沉浮运动,对其非定常气动力特性进行了数值研究。研究采用了动态混合网格技术以及非定常数值计算方法。为了提高动态混合网格的变形能力,采用了基于径向基函数的插值方法求解空间点的位移,翼型整个俯仰-沉浮运动周期内计算网格均维持了较好的质量,没有发生网格重构。非定常算法方面,通过约束单元边界面的法向速度从而满足了运动网格下几何守恒律的要求。空间离散采用了二阶的有限体积格式,时间离散则采用了双时间步和BLU-SGS相结合的隐式时间推进策略。计算得到了不同下拍时间、不同拍动角等条件下的升力、推力以及能耗,对其升力、推力产生机制进行了分析,并通过对气动力以及流场进行对比,分析了各拍动参数的影响。计算结果表明,翼型自身的"静态因素"是其产生升力的主要原因,非定常流动对增加升力起到了促进作用,而下拍时间、拍动角等运动参数对翼型的气动性能影响较大。当下拍时间占到整个拍动周期的约65%-70%时,单位能耗下的时均升力最大,该结论和观测数据较为一致。此外,通过对比分析得到了一组具有较好气动特性的拍动角参数,为后续针对三维问题的研究提供了参考。     

融合体布局超临界翼型设计研究

针对跨声速融合体布局飞机,构建了所需超临界翼型的气动优化设计问题;基于"人在回路"的多轮迭代优化理念建立了翼型多目标多约束优化设计平台,其中:翼型参数化采用基于型函数/类函数变换的参数化方法,气动力解算器采用基于三维雷诺平均Navier-Stokes方程(RANS)的求解程序CCFD-MB,计算网格生成与网格更新采用基于无限插值和椭圆光滑算法的C型结构网格程序,优化算法采用多目标优化遗传算法NSGA-II。设计结果显示,优化翼型具有常规超临界翼型的基本气动力特征,满足阻力发散、抖振边界和低速特性要求;同时具有前加载特征,低头力矩大幅低于常规超临界翼型,满足融合体布局气动力要求。     

非结构网格在积冰数值模拟中的应用

基于非结构网格对翼型表面的积冰进行了数值模拟研究.通过求解气体欧拉方程得到翼型的绕流流场,使用穿透型壁面边界模拟水滴对翼型表面的撞击,采用迭代法对水滴控制方程进行数值求解,得到翼型表面的水滴收集特性;采用不同的积冰冻结模型模拟霜冰和明冰的冻结过程;利用冰层时间推进法模拟积冰过程,基于积冰法向生长假设生成积冰外形.对NA...     

薄翼型阵风响应的数值模拟

 采用有限体积法来数值求解非定常Euler和N-S方程，并在非定常计算中通过引入“网格速度”来计入阵风的影响，实现了对刚性薄翼型阵风响应的有效数值模拟。采用该方法来对薄翼型迎角突然增大的阵风、陡边型阵风这两种典型阵风情况下的翼型气动力响应进行了计算，结果表明：计算结果与现有的一些理论结果都符合较好，与文献中计算结果也吻合。同时通过计算也表明厚翼型的阵风响应需根据工况开展进一步的计算分析。     

多段翼型局部主动变形流动控制的非定常数值模拟

对30P30N三段翼型失速攻角附近的分离流动进行了数值模拟研究。为了抑制大攻角时背风区的流动分离,在主翼段上表面引入了行波壁变形模型和抛物型局部主动振动模型,利用作者以往发展的动态混合网格技术和相应的非定常计算方法,对变形过程中的非定常分离流动进行了数值模拟,分析了各种变形参数对流动分离的影响。计算结果表明,在适当的条件下,局部主动变形能够抑制翼型背风区的分离,由此可以起到增升减阻的作用,改善翼型的气动性能。     

Numerical analysis of broadband noise reduction with wavy leading edge

Large Eddy Simulation(LES) is performed to investigate the airfoil broadband noise reduction with wavy leading edge under anisotropic incoming turbulence. The anisotropic incoming turbulence is generated by a rod with a diameter of 10 mm. The incoming flow velocity is 40 m/s and the corresponding Reynolds numbers based on airfoil chord and rod diameter are about 397000 and 26000, respectively. The far-field acoustic field is predicted using an acoustic analogy method which has been validated by the experiment. A straight leading edge airfoil and a wavy leading edge airfoil are simulated. The results show that wavy leading edge increases the airfoil lift and drag whereas the lift and drag fluctuations are substantially reduced. In addition, wavy leading edge can significantly change the flow pattern around the leading edge and a pair of counter-rotating streamwise vortices stemming from each wavy leading edge peak are observed.An averaged noise reduction of 9.5 dB is observed with the wavy leading edge at the azimuthal angle of 90°. Moreover, the wavy leading edge can mitigate noise radiation at all the azimuthal angles without significantly changing the noise directivity. The underlying noise reduction mechanisms are then analyzed in detail.     

低雷诺数下层流分离的等离子体控制

为有效控制层流分离特性,消除或减弱低雷诺数时小迎角下的升力非线性现象,改善翼型升力特性,并通过翼型的上表面转捩带与油流显示测量对等离子体激励控制机理进行阐述,对厚度为16%椭圆翼型低雷诺数下的气动特性进行了风洞试验研究。在此基础上,在上表面前缘10%弦长处布置激励器,通过压力分布测量观察等离子体激励对层流分离的影响。试验结果表明：当翼型上表面仅发生层流分离时,等离子体激励和转捩带的作用类似,可以有效延迟或者消除后缘层流分离,从而增加升力;当翼型上表面出现层流分离气泡并发生再附现象时,等离子体可以有效减小或者消除层流分离泡的范围,从而减小升力;通过控制层流分离,占空循环等离子体激励可以实现对低雷诺数小迎角下的升力的线性控制。     

多段翼型地面效应数值模拟与分析

通过数值模拟方法研究多段翼型的地面效应,采用有限体积法求解质量加权平均平均Navier-Stokes(N-S)方程,湍流模型采用Spalart-Allmaras(S-A)模型,采用滑移壁面模拟地面的相对运动.计算结果分析表明:随着飞行高度的降低,多段翼型的升力、阻力和低头力矩均减小;升力系数的减小主要是由于地面效应引起...     

低雷诺数下微型飞行器主动变形机翼非定常气动特性数值分析

以微型飞行器的气动外形设计为研究背景,通过数值求解N-S方程分析了主动变形柔性机翼的气动特性。空间离散格式采用中心有限体积法,时间推进采用双时间推进方法,其中子迭代过程由隐式LU-SSOR方法完成。在动网格技术的基础上,模拟分析了薄翼面上作行波运动的流场。计算结果发现,该行波翼型的气动特性相比于刚性翼型有明显改善,有效地抑制了大攻角下大尺度流动分离,同时升阻比提高了约35%,起到了显著的增升减阻的作用。本文的工作对于进一步的柔性机翼优化具有良好的理论研究与数值模拟的参考价值。     

计算谐振机翼气动力的激波膨胀波级数展开法

本文根据激波-膨胀波理论,忽略了反射波对翼面压强分布的影响,把斜激波之后和通过一族主特征线的压强展成级数,给出了沉浮和俯仰谐振翼型表面压强分布的近似解析解。然后应用片条理论导得三维机翼振动导数的闭式解。文中给出三角形机翼俯仰导数的几个实例,数值结果与现有其它方法的结果有很好的一致性。     

用直接测力方法进行软质翼型的气动特性研究

在西北工业大学NF-3低速风洞运用翼型气动力直接测量的方法对软质翼型进行风洞试验研究,对比了软、硬质翼型模型的试验结果。结果表明：软质翼型模型与硬质翼型模型在相同风速下具有不同的气动力特性。在一定风速下,软质翼型模型的表面会发生变化,从而影响了气动力。由于该影响非常复杂,因此在研究软质翼型模型的气动特性时进行风洞试验是必要的。     

抽吸控制对低雷诺数下翼型分离流动的影响

为了系统研究抽吸区域相对于分离点的作用位置、孔间距和孔径对抽吸控制翼型分离流动效果的影响,以NACA0012翼型表面分离流动为基准状态,在其吸力面设计了局部多孔抽吸结构,采用AUSM⁺-up格式、大涡模拟方法和双时间步长（LU-SGS）隐式算法,对低雷诺数下多孔分布式抽吸结构对流动分离的控制效果进行了数值研究。研究结果表明：当抽吸区域位于分离点之后时,抽吸控制效果最好;抽吸系数不仅存在一个下限值以达到快速、有效的控制效果,而且有一个上限值以保证抽吸控制品质因数（FOM）大于1;孔间距和孔径对翼型气动性能的影响较小,但对FOM分布的影响较大。     

抽吸控制对低雷诺数下翼型分离流动的影响简

 为了系统研究抽吸区域相对于分离点的作用位置、孔间距和孔径对抽吸控制翼型分离流动效果的影响,以NACA0012翼型表面分离流动为基准状态,在其吸力面设计了局部多孔抽吸结构,采用AUSM⁺-up格式、大涡模拟方法和双时间步长（LU-SGS）隐式算法,对低雷诺数下多孔分布式抽吸结构对流动分离的控制效果进行了数值研究。研究结果表明：当抽吸区域位于分离点之后时,抽吸控制效果最好;抽吸系数不仅存在一个下限值以达到快速、有效的控制效果,而且有一个上限值以保证抽吸控制品质因数（FOM）大于1;孔间距和孔径对翼型气动性能的影响较小,但对FOM分布的影响较大。     

Effect of flexibility on unsteady aerodynamics forces of a purely plunging airfoil

Introducing flexibility into the design of a vertically flapping wing is an effective way to enhance its aerodynamic performance. As less previous studies on the aerodynamics of vertically flapping flexible wings focused on the lift generated in a wide range of angle of attack·a 2D numerical simulation of a purely plunging flexible airfoil is employed using a loose fluid–structure interaction method. The aerodynamics of a fully flexible airfoil are firstly studied with the flexibility and angle of attack. To verify whether an airfoil could get aerodynamic benefit from the change in structure, partially flexible airfoil with rigid leading edge and flexible trailing edge were further considered. Results show that flexibility could always reduce airfoil drag while lift and lift efficiency both peak at moderate flexibility. When freestream velocity is constant, lift is maximized at a high angle of attack about 40° while this optimal angle of attack reduces to 15° in drag-balanced status. The airfoil drag reduction, lift augmentation as well as efficiency enhancement mainly attribute to the passive pitching other than the camber deformation. Partially deformed airfoil with the longest length of moderate flexible trailing edge can achieve the highest lift. This study may provide some guidance in the wing design of Micro Air Vehicle (MAV).