微型凸起流动控制用于翼型的减阻特性研究

微型凸起作为减小阻力的一种有效措施已经备受关注.开展了等熵压缩"弱激波"干扰鼓包用于RAE2822超临界翼型的减阻作用机制研究以及NACA0012对称翼型表面脊状结构减阻特性的数值模拟研究.结果表明:通过鼓包参数最优匹配,可达到弱化激波、减小波阻、提高升阻比、延缓抖振边界等目的;同时,通过对比不同脊状结构、不同网格密度对计算结果的影响,总结了多个速度下脊状表面的减阻规律.所得结论为进一步开展微型凸起类流动控制用于机翼的减阻特性研究奠定了坚实基础.     

激波控制鼓包对跨声速抖振影响的数值研究

在跨声速飞行时,激波控制鼓包不仅能够减弱机翼上表面的激波强度从而降低波阻,对跨声速抖振也有一定的改善作用。通过URANS方法数值模拟来探究二维激波控制鼓包对OAT15A超临界翼型跨声速抖振性能的影响规律,并研究以巡航设计点减阻与抖振状态减振2种目标设计的鼓包的区别。以巡航设计点减阻优化设计出的鼓包,在抖振条件下,能够推迟了翼型上表面的压力恢复,减弱了激波与边界层的相互干扰作用,达到减弱抖振幅度的效果,然而不能对抖振实现完全抑制。通过改变鼓包相对位置、高度和长度计算得到鼓包参数对抖振的影响规律,分析典型流场得到鼓包抑制抖振现象的工作机理是：鼓包减弱了激波强度的同时,阻碍了鼓包尾部边界层向上游移动与激波相互干扰,从而稳定了激波抖振现象。另外,基于巡航设计点减阻设计的2个鼓包相对参考位置距离分别为0.04c和0.10c(c为翼型弦长),与同等高度鼓包在抖振状态完全抑制抖振且不降低升力的位置范围的[-0.01,0.02]c和[0.01,0.08]c不同,二者位置最小相差0.02c,而鼓包这段距离差异对巡航特性和抖振性能都有着重要影响作用。总而言之,以巡航设计点减阻与抖振状态减振2种目标设计得到...     

超临界翼型低雷诺数流动分析及优化设计

以高空长航时无人机(UAV)翼型研究为背景,对超临界RAE2822翼型在高空高亚声速下的低雷诺数气动特性进行了数值模拟及优化设计研究。采用求解雷诺平均N-S方程的有限体积法,对典型低雷诺数下RAE2822翼型绕流进行数值模拟,验证了SST k-ω湍流模型的可靠性和准确性;基于不同高度不同雷诺数下RAE2822翼型的计算气动力对比分析,研究了高度增大所带来的低雷诺数效应;通过对低雷诺数下超临界翼型表面流场结构及流动机理的详细分析,提出了一种弱化激波的翼型设计思想,并通过优化算例验证了该思想的可行性。     

开缝空腔抑制翼型跨声速抖振的数值模拟

采用非定常雷诺平均Navier-Stokes(URANS)方法计算了18%双圆弧翼型的跨声速抖振特性,分析了翼面激波振荡及流场结构演化的特点,研究了在翼型表面开通气空腔抑制跨声速抖振的可行性,对空腔深度、开缝数目对激波振荡的抑制效果进行了对比分析。计算发现,18%双圆弧翼型的跨声速激波自激振荡只有向前的运动,没有向后的运动,开缝空腔能够抑制翼型跨声速抖振,但对抖振频率影响不大;空腔深度大,抑制效果好,但空腔深度变化对振荡频率影响不大;开2、3、4个槽缝抑制抖振的效果差别不大,开缝数量对抖振频率影响不大。     

翼型激波抖振的无模型自适应控制

由于受到流动环境中的不确定因素影响,针对流动系统设计的闭环控制律需要考虑这些不确定因素,自动适应随机扰动和突发扰动.基于流动模型所设计的闭环最优控制,虽然可以以较小的控制量获得满足目标函数的最优控制效果,但随着流动模型的引入,随之而来的建模精度、未建模动态等问题也会出现,固定的流动模型也会限制控制系统的自适应能力.针对翼型在跨声速流动中遇到的激波抖振问题,为了消除不同来流状态的激波抖振带来的脉动载荷,开展了基于数据驱动方法的无模型自适应控制.流场数值仿真采用URANS方法,作动机构采用尾缘舵面,以升力系数作为反馈信号.当流动状态变化时,无模型自适应控制利用输入输出数据在线将流动系统等价转化为动态线性化数据模型,最小化性能指标得到控制律,使系统自动地工作于最优或接近最优的状态.时域仿真结果显示,无模型自适应控制效果优于开环控制和比例控制,并且当来流状态随时间发生变化时,无模型自适应控制也能够完全消除抖振脉动载荷.     

机翼跨声速抖振研究进展

跨声速飞行中,激波附面层干扰会引起激波周期性自激振荡,这种现象称之为跨声速抖振。跨声速抖振引起的脉动载荷有可能造成结构疲劳甚至引发飞行事故,该问题一直是航空领域的研究难点。本文从风洞试验和数值模拟两个方面总结了跨声速抖振研究的主要方法,重点综述了近几年在跨声速抖振机理研究方面以及结构运动对抖振响应特性的影响方面取得的新进展,简要总结了跨声速抖振的被动控制及主动控制方法。在此基础上,提出了跨声速抖振后续研究的重点方向:抖振机理的进一步探究,结构弹性特征对抖振特性的影响以及闭环主动控制等。     

基于谐振舵面的跨声速抖振抑制探究

跨声速抖振引起的非定常脉动载荷会造成飞行器结构疲劳甚至引发飞行事故,所以跨声速抖振的控制研究逐渐成为航空领域的热点。采用基于Spalart-Allmaras(S-A)湍流模型的非定常雷诺平均方程开展了基于谐振舵面的跨声速抖振抑制研究。首先验证静止NACA0012翼型的抖振边界和频率特性,然后分别从舵偏平衡位置、舵偏幅值、频率以及相角等角度研究了谐振舵面的控制效果。舵偏平衡位置等效于减小了翼型的有效迎角;幅值和频率对抖振抑制效果影响较大,当舵面振荡频率与抖振频率接近时发生共振现象;相角对控制效果有一定影响,在270°相角附近,升力系数幅值减小了60%。在合适的舵偏幅值、频率以及相角组合下,谐振舵面有可能成为跨声速抖振的有效开环控制策略。     

基于涡动力学的翼型虚拟质量的理论与计算模型

从涡量矩定理出发,通过保形变换法推导出翼型在不可压流场中的虚拟质量。利用变换前后的涡层守恒关系、结合边界元方法的原理,革新了翼型外平面到单位圆外平面的保形映射算法。理论推导结果表明,翼型的虚拟质量只与不可压流体的密度及翼型本身的形状有关。算法在变换对称与非对称Joukowski翼型、RAE2822翼型,以及若干NACA系列翼型时,节点位置的偏差在0.2%内,虚拟质量的计算值与理论值吻合良好。     

基于POD和DMD方法的跨声速抖振模态分析

跨声速抖振现象是由于非定常跨声速流动中激波的自激振荡而引起的结构强迫振荡,这种现象在跨声速飞行器中普遍存在,对飞机的结构强度和疲劳寿命有不利影响。基于模态分解的分析方法是进一步发展抖振控制手段的有效工具。本文通过两类典型模态分析方法（本征正交分解（POD）和动态模态分解（DMD））对OAT15A翼型的跨声速抖振现象进行分析,通过对模态频率、翼面压力分布、流场重构误差等方面的研究,将两种模态分解方法进行对比。发现基于频率特征的DMD方法能够准确捕捉抖振的临界稳定特征和抖振主频的典型模态,同时能够更准确反映流场变量在激波间断附近随时间的变化过程;而POD方法尽管在流场重构时具有较小的总体误差,但对激波附近压强随时间的变化历程拟合较差。     

固定转捩和自由转捩流动预测方法验证

采用固定转捩位置和基于SST k-ω湍流模型的γ-Reθt转捩模型自由转捩模拟两种方法研究了RAE2822翼型流动转捩问题。先采用γ-Reθt转捩模型确定转捩位置,再采用所获得的转捩位置实施固定转捩位置计算,并与全湍流计算结果和其他软件计算结果进行了对比。表明两种方法均能有效捕捉超临界翼型流动转捩现象,固定转捩位置模拟结果基本上重现了自由转捩预测结果。     

Transonic buffet control research with two types of shock control bump based on RAE2822 airfoil

Current research shows that the traditional shock control bump (SCB) can weaken the intensity of shock and better the transonic buffet performance.The author finds that when SCB is placed downstream of the shock,it can decrease the adverse pressure gradient.This may prevent the shock foot separation bubble to merge with the trailing edge separation and finally improve the buffet performance.Based on RAE2822 airfoil,two types of SCB are designed according to the two different mechanisms.By using Reynolds-averaged Navier-Stokes (RANS) and unsteady Reynolds-averaged Navier-Stokes (URANS) methods to analyze the properties of RAE2822 airfoil with and without SCB,the results show that the downstream SCB can better the buffet performance under a wide range of freestream Mach number and the steady aerodynamics characteristic is similar to that of RAE2822 airfoil.The traditional SCB can only weaken the intensity of the shock under the design condition.Under the off-design conditions,the SCB does not do much to or even worsen the buffet performance.Indeed,the use of backward bump can flatten the leeward side of the airfoil,and this is similar to the mechanism that supercritical airfoil can weaken the recompression of shock wave.     

Flow characteristics around airfoils near transonic buffet onset conditions

In transonic flow, buffet is a phenomenon of flow instability caused by shock wave/boundary layer interaction and flow separation. The phenomenon is common in transonic flow, and it has serious impact on the structural strength and fatigue life of aircraft. In this paper, three typical airfoils: the supercritical OAT15A, the high-speed symmetrical NACA64A010, and the thin, transonic/supersonic NACA64A204 are selected as the research objects. The flow fields of these airfoils under pre-buffet and buffet onset conditions are simulated by Unsteady Reynolds Averaged Navier-Stokes (URANS) method, and the mode analysis of numerical results is carried out by Dynamic Mode Decomposition (DMD). Qualitative and quantitative analysis of the shock wave motion, shock wave intensity, shock foot bubble and trailing edge separation, and pressure coefficient fluctuation were performed to attain deep insight of transonic buffet flow features of different airfoils near buffet onset conditions. The results of DMD analysis show that the energy proportion of the steady mode of these airfoils decreases dramatically when approaching the buffet onset angle of attack, while the growth rate of the primary mode increases inversely. It was found that at the onset of buffet, there exist different degrees of merging behavior between shock foot bubble and trailing edge separation during one buffet cycle, and the instability of shock wave and separation induced shear layer are closely related to the merging behavior.     

跨声速机翼抖振初始迎角N-S方程定常计算分析

 应用代数法和椭圆型方程优化相结合的方法生成翼身组合体块结构网格。采用有限体积空间离散法和五步Runge-Kutta显式时间推进法求解N-S方程。基于雷诺平均N-S（RANS）方程,选用Spalart-Allmaras（SA）一方程模型模拟紊流。根据N-S方程定常计算结果,采用升力曲线、俯仰力矩曲线、后缘压力发散、跨声速激波位置以及机翼表面极限流线等几种判据,对跨声速机翼的抖振初始迎角进行了合理的预测分析。     

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

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

Shock control on a swept wing

Shock and boundary layer control by contour bumps and local boundary layer suction have been investigated experimentally and numerically on a transonic swept wing. Additional 2-D numerical investigations were performed for the airfoil, corresponding to the wing. The investigations were primarily stimulated by the question concerned with the influence of sweep on the bump effectiveness. This influence has been found to be rather small; the drag reduction by the bump is slightly lower for the swept wing than for the airfoil. A location of the bump in the shock region has shown its effectiveness for reducing shock strength and hence wave drag. A position of the bump downstream of the shock wave has been shown to reduce viscous drag and to postpone buffet-onset to higher lift coefficients. Furthermore, the results indicate that boundary layer suction is a powerful device for drag reduction, but the effectiveness decreases with increasing Reynolds number. Higher effectiveness of suction can be attained, when it is coupled with a contour bump. The parameters height and position (relative to the shock) of the bump, optimized in terms of drag, depend on the shock strength; an influence of the boundary layer thickness upstream of the shock on the optimal bump parameters has not been found. A possibility to control an adaptive bump, mounted on an aircraft wing, is to employ the trailing edge pressure.     

基于混合笛卡尔网格的翼型绕流数值模拟

在对重叠网格技术研究的基础上,一套以笛卡尔网格为主体、贴体的结构网格为特殊固壁边界的混合笛卡尔网格(Hybrid Cartesian Grid,HCG)生成方法得到了发展,同时基于流场特征的笛卡尔网格自适应技术和有限体积方法的雷诺平均Navier-Stokes(RANS)的数值求解方法也得到了开发。通过对经典跨音速RAE2822翼型、双NACA0012翼型和高升力两段NLR-7301翼型进行绕流数值模拟,获得了较好的结果,验证了方法的可行性与精确性,为研究二维复杂外形流动问题提供了一种解决方法。