低速轴流压气机旋转失速边界的模型描述

从低速轴流压气机完整旋转失速过程所具有的突变性和迟滞性等典型动力学特征出发,通过对其拓扑性质的对比分析,证明压气机旋转失速与突变理论尖点突变模型在拓扑结构上具有相似性;以拓扑不变性定理为基础,通过发展有效的拓扑映射方法,建立了尖点突变模型的突变点集与压气机旋转失速突变点集之间的拓扑映射关系.在两台低速轴流压气机上的模型实际应用效果表明:尖点突变模型能够定性地描述具有突变型失速特点的低速压气机不同转速工况下完整旋转失速过程的突变性、迟滞性以及分岔特性,通过该发展的拓扑映射方法,可以将实验获取的失速点和恢复点完全映射到模型线上,这表明尖点突变模型能较为准确地描述旋转失速的失速边界.      

Hysteresis behaviors of compressor rotating stall with cusp catastrophic model

Rotating stall is a complex nonlinear dynamic phenomenon which is always characterized by catastrophe and hysteresis in high aerodynamic-loading compressor. Exploring the key contributing factors and characteristic rules of hysteresis is very important for compressor design and flow instability control. In this paper, a novel model method is proposed to analyze the hysteresis behaviors to extend the understanding of compressor rotating stall. The equilibrium states of compressor system under different conditions are first described based on Moore-Greitzer model. Then, through assessing the stability of the equilibrium points by Liapunov's theorem, the ratio of shutoff head to compressor characteristic semi-height is found to affect the stall hysteresis: the size of hysteresis loop will gradually decrease, even disappear with the increase of the ratio. Combing the effects of both the ratio and throttle coefficient, the hysteresis behaviors of compressor stall under multi-parameters can be found to be consistent with the topological properties of cusp catastrophic model by Thom’s catastrophe theory. Finally, according to topological invariant rules, from the perspective of potential function, the equilibrium surface equation of compressor system is developed by standard cusp catastrophic model to describe the various hysteresis behaviors of compressor rotating stall along different control routes.     

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

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

翼型动态失速等离子体流动控制试验

针对动态失速引起的翼型气动性能恶化的问题,利用小型化的激励电源和介质阻挡放电等离子体激励器,借助动态压力测量和外触发式粒子图像测速（PIV）等手段开展了翼型动态失速等离子体流动控制试验研究。结果表明,等离子体气动激励能够有效控制翼型动态失速,改善平均气动力,提高翼型气动效率,减小气动力随迎角变化的迟滞区域。等离子体诱导出前缘附近的贴体翼面涡,促进分离流再附;增加了上翼面0.2~0.4弦长区域的吸力,减小了升力系数功率谱密度（PSD）分布的二、三、四阶能量幅值,在研究工况下实现了平均升力系数增加7.1%、失速迎角推迟1.3°和迟滞区域减小4.5%的明显控制效果;4°~9°迎角段,等离子体使得翼型平均阻力系数减小40%。此外,振荡频率增加使翼型绕流的非定常性增强,较高雷诺数下的翼型动态分离涡更加难以被抑制,均需要增加等离子体激励强度才能达到较好的控制效果。     

RANS analysis of the low-Reynolds number flow around the SD7003 airfoil

This paper presents a numerical analysis of the incompressible flow at Reynolds number 6.0×¹⁰⁴ around the Selig–Donovan 7003 airfoil. The airfoil performances have been computed by the Reynolds averaged Navier–Stokes equations and large eddy simulations. The airfoil stall and preliminary post-stall have been obtained by both the methods. Some limitations of the RANS turbulence models for low-Reynolds number flows have been overcome by the κ–ω SST-LR model, a recent modification of the well-known SST model. Large-eddy simulations have also been performed for a more detailed analysis of the results. The relevance in the stall mechanism of the laminar separation bubble arising on the airfoil is highlighted. The stall occurs when the laminar bubble present in the leading edge zone and a separated region forming on the central part of the airfoil join together. The κ–ω SST-LR model returns the same stall mechanism as the large eddy simulation. Flows at low-Reynolds numbers can be simulated by the RANS methods, but the choice of the turbulence model is crucial. The κ–ω SST-LR model has provided results in good agreement with the large eddy simulation and the available experimental data.     

基于IDDES方法的翼型结冰失速分离流动数值模拟

应用基于k-ωSST湍流模型的IDDES(Improved Delayed Detached Eddy Simulation)方法,就失速点附近翼型前缘典型双角状积冰导致的复杂分离流动进行了数值模拟研究.通过与风洞试验结果进行对比,表明对于此类分离流动问题,IDDES方法能够在壁面附近取得良好的速度预测结果,有效解析分离区域内的中小尺度湍流结构,较为准确地描述大尺度时均分离泡的再附位置和形态特征,适用于翼型结冰后复杂流动的精细分析.同时计算结果显示当此带冰翼型位于失速点附近时,角状冰后方脱落剪切层内部的旋涡不稳定析出和输运过程促进了外部流动与回流区域流动间的掺混,将导致流动发生非定常再附现象.     

Experimental study of ice accretion effects on aerodynamic performance of an NACA 23012 airfoil

In this paper,the effects of icing on an NACA 23012 airfoil have been studied.Experiments were applied on the clean airfoil,runback ice,horn ice,and spanwise ridge ice at a Reynolds number of 0.6 106 over angles of attack from 8° to 20°,and then results are compared.Generally,it is found that ice accretion on the airfoil can contribute to formation of a flow separation bubble on the upper surface downstream from the leading edge.In addition,it is made clear that spanwise ridge ice provides the greatest negative effect on the aerodynamic performance of the airfoil.In this case,the stall angle drops about 10° and the maximum lift coefficient reduces about50% which is hazardous for an airplane.While horn ice leads to a stall angle drop of about 4° and a maximum lift coefficient reduction to 21%,runback ice has the least effect on the flow pattern around the airfoil and the aerodynamic coefficients so as the stall angle decreases 2° and the maximum lift reduces about 8%.     

前缘下垂结合内吹式襟翼失速特性研究

内吹式襟翼具有高效的增升能力,但失速迎角在较高的吹气动量系数下下降明显,为改善其失速特性,研究内吹式襟翼加装前缘下垂后的失速特性。对前缘下垂结合无缝襟翼的亚声速翼型在环量控制作用下的流场进行数值模拟,研究吹气动量系数对失速特性的影响规律,前缘刚性偏转、弯度变化和厚度变化对失速特性的改善作用,以及改变襟翼偏角研究前缘下垂的作用效果。结果表明：随着吹气动量系数的增加,失速迎角先迅速下降再略微增加;前缘下垂装置减小了翼型上表面逆压梯度,延缓了翼型边界层动量厚度随迎角增加而增加的趋势,能有效提高失速迎角;通过逐渐改变前缘表面曲率,实现了前缘下垂设计对失速特性改善的最好效果。     

基于联合射流的翼型动态失速流动控制研究

动态失速控制对于提高翼型气动特性具有重要意义。采用联合射流方法对翼型俯仰动态失速控制进行数值模拟,完成两方面的研究:一是射流关闭时射流通道对动态失速特性的影响,二是射流开启时不同射流动量系数对动态失速控制的影响和分析。结果表明:射流关闭时,射流通道的存在对翼型上仰过程中附着流阶段的气动特性产生不利影响,使得升力系数明显下降,但是对翼型下俯过程中失速分离流阶段的气动特性影响不明显;射流开启后,动态失速特性得到极大改善,迟滞环面积显著减小,升力增加,阻力减小,且阻力和力矩的峰值显著减小,原基准翼型力矩曲线的负阻尼区域消失。     

飞虫粘附翼型对翼型气动特性影响的实验研究

在西北工业大学NF-3低速风洞进行了飞虫粘附翼型对翼型气动性能影响的风洞实验研究。结果表明：模型表面粘附的飞虫数量累计到一定程度时,将会导致翼型表面50%区域以上的面积发生分离,引起翼型失速,并且实验前模型表面的飞虫数量会改变翼型的失速迎角,因此实验前必须将模型擦拭干净。鉴于飞虫粘着数量的不确定性,对于翼型在复杂环境下使用时其气动性能的变化需要加以关注。     

基于CST方法的高空低雷诺数吸附式叶型耦合优化设计

研究了一种基于类别形状函数变换（CST）方法的吸附式叶型优化设计方法,该方法可以在高空低雷诺数条件下对叶型和抽吸方案耦合优化.结果表明:优化之后在20km高空低雷诺数条件下总压损失降低了65%,静压升提高了0.02,气动性能得到较大提升.而且由于优化过程中罚函数的引入使得优化后吸附式叶型在地面条件下性能也有所提高.对于高空低雷诺数条件下吸附式叶型在抽吸位置之前适当的增加叶型负荷,再通过抽吸来控制附面层,效果最优.并且最佳抽吸位置位于层流分离泡作用区域内.在层流分离泡作用区域内抽吸可以完全消除层流分离泡对叶型性能的影响,并且可以较好控制附面层位移厚度和动量厚度的增加,有效地减小附面层内的动量损失.      

翼伞弧面下反角、翼型和前缘切口对翼伞气动性能的影响

为了研究翼伞弧面下反角、翼型和前缘切口对翼伞气动性能的影响,对带气室的、展弦比为3的不同特征几何参数翼伞模型的流场进行了三维、定常数值模拟。运用有限体积法对三维坐标系下不可压雷诺时均Navier-Stokes (RANS)方程进行了直接求解,采用剪切应力输运(SST)k-ω两方程湍流模型对湍流进行模拟。数值模拟得出的原始翼伞的气动性能参数与试验数据在总趋势上符合很好,不同几何参数翼伞模型计算结果表明:翼伞弧面下反角越大,升力及诱导阻力越小,升阻比变化不大;前缘半径、厚度小的翼伞翼型,阻力更小,升阻比大;前缘切口对翼伞影响区域限于前缘附近,压力分布同干净翼类似,降低了其失速迎角,对升力影响不大,但明显增大阻力。该数值方法可为进一步研究更多不同几何参数的翼伞模型提供参考。     

任意翼型有小分离气泡跨音速绕流的迭代算法

 本文给出一种带小分离气泡的任意翼型粘性跨音速绕流的计算方法,采用有粘-无粘干扰迭代的概念。无粘流的全速势方程用AF差分格式在保角变换法生成的计算网格中求解,粘流附面层方程用C-S盒式法求解,用逆算法消除分离点处的奇性。本文对Ma_∞=0.8,Re_∞=2×10~6,迎角α=3.5°和4°的NACA64A010翼型粘性绕流进行了计算,结果与实验相比较,吻合良好。     

大振幅下薄翼和流线型箱梁的气动迟滞研究

回顾了机翼和箱梁的气动迟滞效应研究现状,介绍了大振幅下获取薄翼和流线型箱梁自激气动力的风洞试验。研究结果表明：薄翼在超过失速角的振幅条件下,升力并没有出现明显的失速现象,但气动力的高阶谐波分量显著增加;大振幅条件下,流线型箱梁的气动力高阶谐波分量也比较显著,并以第2和第3阶谐波分量为主。此外,在大振幅条件下,流线型箱梁的气动力矩迟滞曲线可出现“8字环”,即气动力在一个振动周期内既做了正功也做了负功;其中,迟滞曲线“8字环”的气动正功部分随着振幅和折算风速的增大而增加。该现象可导致桥梁也出现类似于失速机翼的极限环震荡。最后,基于不同振幅下流线型箱梁的力矩迟滞曲线,简要讨论了大跨度桥梁在颤振后状态可能出现的振动形式和气动稳定性。     

表面凹凸充气机翼的气动特性研究

充气机翼在变体飞机和飞艇中具有潜在的应用前景。充气机翼的结构特征与传统硬质机翼显著不同,其外形与传统机翼相比最大的差异在于表面的片条状鼓包,这种外形带来的气动特性、气弹行为等越来越受到人们的关注。以NACA0015翼型为原形,设计制作充气机翼模型,并利用CFX对传统光滑的NACA0015翼型和凹凸表面的0015F2翼型进行定常和非定常气动行为分析。结果表明:充气机翼的凹凸表面外形增加翼型的失速攻角,但其升力线斜率及升阻比都较光滑翼型要小;0015F2翼型的速度梯度过度区大于NACA0015翼型;0015F2翼型在每一个凹槽区生成驻涡,驻涡的存在使得充气机翼的附面层呈现紊流附面层的特性,驻涡的外移改变了机翼后缘的尾涡形成,推迟了分离,使得失速攻角增大。     

某高负荷子午加速风机气动设计与数值模拟

兼顾民用风机的气动性能与制造成本,采用三种方案(等厚度串列静子方案、叶型串列静子方案和等厚度单排静子方案)对某高负荷子午加速风机进行了气动设计和数值模拟研究.结果表明采用叶型叶片与等厚度叶片相比对风机气动性能的改善很小,近设计点时叶型串列静子方案的压比和效率仅比等厚度串列静子方案高0.18%和0.8%.串列静子抑制了单排静子通道内的角区分离,减小了通道内的堵塞,使等厚度串列静子方案的失速裕度比等厚度单排静子方案提高了11%.综合考虑气动性能与制造成本选定等厚度串列静子方案为风机的最终方案,其不同转速的特性表明该方案的性能完全满足设计指标的要求.      

薄翼失速翼型前缘分离泡对失速特性的影响

DES方法结合了RANS(Reynolds-averaged Navier-Stokes)和LES(Large Eddy Simulation approaches)的优点。在近壁面它体现为RANS模型的特点而在远离物面处又起到LES的亚格子模型的特性。论文应用DES(Detached-EddySimulation)方法讨论了影响薄翼失速的分离泡对翼型的升力特性影响。     

低雷诺数翼型局部振动非定常气动特性

针对低雷诺数翼型特殊的气动特性,采用基于动网格的非定常数值模拟方法,研究翼型表面不同弦向位置的局部蒙皮以不同频率及振幅振动时对低雷诺数翼型气动特性及流场结构的影响,揭示蒙皮振动增升减阻的机理。研究表明,在低雷诺数条件下局部蒙皮振动可有效提高翼型气动特性,与刚性翼型相比蒙皮局部振动可使翼型升力系数提高,阻力系数降低,升阻比提高。振动位置对翼型气动特性及流场结构有显著的影响,振动表面位于翼型前缘附近或位于层流分离泡中心时可有效控制翼型层流分离,从而提高翼型气动特性。振动频率对翼型表面层流分离及转捩位置均有显著的影响,随着振动频率增加,翼型气动特性出现最优值。与刚性翼型相比,表面振动使翼型转捩位置略向上游移动,摩擦阻力增加,但振动使等效翼型相对厚度减小,压差阻力明显减小。在小幅振动范围内,随着振幅增加,流场非定常特性更加显著,翼型升阻比增加。     

低速翼型分离流动的等离子体主动控制研究

为了研究等离子体激励器的放电形式及其诱导气流的规律,以及翼型迎角、自由来流速度分别对翼型流动分离抑制效果的影响。在低速、低雷诺数条件下利用介质阻挡放电等离子体激励器对NACA0015翼型进行了主动流动控制研究。结果表明：介质阻挡放电的形式为丝状放电;等离子体激励器诱导气流的方向由裸露电极指向覆盖电极,由电极的布置方式决定,与接线方式无关;当来流速度为25m／s,雷诺数为2．03×10^5时,等离子体气动激励可以有效地抑制翼型吸力面的流动分离,翼型最大升力系数增大约为9．7％,翼型l临界失速迎角由17．5°增大到20．5°;翼型失速延迟的真正原因并非单纯的气流加速;等离子体激励器的作用效果随着来流速度的提高而减弱,研究非定常激励或等离子体激励器与流场之间的耦合效应,也许更加具有潜力。     

Parametric analyses on dynamic stall control of rotor airfoil via synthetic jet

The effects of synthetic jet control on unsteady dynamic stall over rotor airfoil are investigated numerically. A moving-embedded grid method and an Unsteady Reynolds Averaged Navier-Stokes(URANS) solver coupled with k-x Shear Stress Transport(SST) turbulence model are established for predicting the complex flowfields of oscillatory airfoil under jet control. Additionally, a velocity boundary condition modeled by sinusoidal function has been developed to fulfill the perturbation effect of periodic jet. The validity of present CFD method is evaluated by comparisons of the calculated results of baseline dynamic stall case for rotor airfoil and jet control case for VR-7 B airfoil with experimental data. Then, parametric analyses are conducted emphatically for an OA212 rotor airfoil to investigate the effects of jet control parameters(jet location, dimensionless frequency, momentum coefficient, jet angle, jet type and dual-jet) on dynamic stall characteristics of rotor airfoil. It is demonstrated by the calculated results that efficiency of jet control could be improved with specific momentum coefficient and jet angle when the jet is located near separation point of rotor airfoil. Furthermore, the dual-jet could improve control efficiency more obviously on dynamic stall of rotor airfoil with respect to the unique jet, and the influence laws of dual-jet's angles and momentum coefficients on control effects are similar to those of the unique jet. Finally,unsteady aerodynamic characteristics of rotor via synthetic jet which is located on the upper surface of rotor blade in forward flight are calculated, and as a result, the aerodynamic characteristics of rotor are improved compared with the baseline. The results indicate that synthetic jet has the capability in improving aerodynamic characteristics of rotor.