低雷诺数下二维翼型绕流的流场特性分析

采用高精度有限差分格式,对低雷诺数下二维翼型绕流进行了直接数值模拟,计算了雷诺数为1.0×104,NACA0012翼型0°,4°以及10°攻角下的流场。计算结果表明:在0°和4°攻角条件下,翼型绕流尾迹区的统计特性相似,0°攻角下的统计量值具有很好的对称性;在距翼型尾缘0.3弦长以后的尾迹区,旋涡排列成类似涡街的结构,涡量的极值、压力的极小值都位于旋涡中心,沿着流向,涡量极值的绝对值逐渐减小,压力的极小值逐渐增大。10°攻角下,翼型上表面从前缘开始分离,尾迹区统计分析结果所得图象与0°和4°完全不同,数值上较后者结果大;在翼型尾缘处,涡量的卷吸,高压力区域的形成,是旋涡脱落的条件,正向和反向旋涡的交替脱落,形成了类似涡街的结构。      

圆柱／翼型干涉流场的试验研究

在风洞出口低速流中,以NACA0012尾缘钝化翼型为模型,利用粒子成像测速系统,研究了圆柱/翼型结构干涉流动时翼型前缘、近壁和尾缘区域的流场.试验结果表明,由于上游圆柱引起的卡门涡街和翼型相互干涉,在翼型前缘存在大尺度涡的变形、拉伸和破裂,在翼型表面近壁区域和尾迹流场中仅存在小尺度湍流涡,由此可以推断翼型前缘可能是干涉噪声的主要声源区.     

翼型近尾迹流动的PIV研究—运动学特性

利用在线式PIV系统(ParticleImageVelocimetry),在低速风洞中对NACA0012翼型在雷诺数2.39×105,0°和4°攻角下的近尾迹流动进行了实验研究。实验结果表明,在较高的雷诺数下翼型近尾迹流动是一种以旋涡的运动学和动力学特性为主导的湍流剪切流。在测量范围内,翼型的尾缘处是近尾迹涡街的形成区;尾缘后0.5倍弦长的区域存在类似于卡门涡街的有序结构,是旋涡发展区域,旋涡具有较好的稳定性;距翼型尾缘0.5倍弦长至1倍弦长的区域,是翼型近尾迹流动由有序走向无序区域,旋涡开始破裂。翼型表面边界层对翼型近尾迹湍流剪切流的演化有重要影响。实验结果还给出了近尾迹流动的平均速度、湍流强度和剪切应变变化率,以及速度脉动量的二阶关联量u'u',u'v'和v'v' 的分布。      

翼型近尾迹流动的PIV研究—动力学机制

利用在线式互相关PIV(ParticleImageVelocimetry)系统,在低速风洞中对NACA0012翼型在雷诺数2.39×105,0°和4°攻角下的近尾迹流动进行了详细测量。实验结果表明,翼型近尾迹存在有序的涡街结构,涡街在尾缘处形成后,在向下游的迁移中,会经历一个发展壮大、失稳破碎的演化过程,流动从有序走向无序。翼型的近尾迹是一种以旋涡的运动学特性和动力学机制为主导的流动现象。本文着重探讨了翼型尾缘处的涡街形成机理,尾迹内的流动机制,以及近尾迹的流动稳定性。      

卡门涡街与翼型的相互作用

本文利用保角变换,用离散涡方法求解了卡门涡街与翼型的相互作用。较详细地分析了它作为偶极子声源的特性,并模拟了在涡翼相互作用过程中尾迹的发生和演化。计算结果表明:在相互作用时,尾迹所吸收的声能可以是正值也可以是负值;升力和声偶极子强度出现了周期性脉动,这种脉动与翼型的振动和噪声有直接的关系。     

S1223翼型低雷诺数下分离泡及气动特性

采用表面测压技术,测量了低雷诺数下(Re=6.0×104、1.0×105、2.0×105)S1223翼型的表面压力分布,通过时均化处理及瞬态处理方法,分别获得了翼型稳态和瞬态压力系数、升力系数,分析了流场结构随雷诺数及攻角的变化规律,研究了雷诺数及攻角对翼型升力的影响机理.结果表明,从时均升力系数随攻角的变化规律来看,S1223翼型在低雷诺数下存在"静态滞回"效应.攻角由负逐渐增大至0°时,下翼面由完全分离转变为出现层流分离泡,随后分离泡逐渐减小直至消失,导致升力系数斜率呈现随攻角逐渐增大的非线性现象.当攻角超过临界攻角后,不同雷诺数下翼型流场结构随攻角的变化规律存在本质不同,Re=6.0×104和1.0×105时,翼型周围流场迅速发生大范围流动分离,升力系数迅速减小;而Re=2.0×105时,上翼面周期性生成短泡,引发低频振荡现象,升力系数呈现准周期性变化,α=16°时上翼面时均流场呈现40％弦长的长泡结构.     

垂直交叉双圆柱绕流数值模拟及涡结构分析

用分块耦合方法数值模拟了垂直交叉双圆柱绕流,两圆柱间距为圆柱直径的5倍,所取雷诺数为200.在交叉结构的中心流场比较复杂,用λ2方法准确地描述了其中的三维涡结构,发现上游圆柱规则的卡门涡街被下游圆柱截断成复杂的缠绕结构.与串列双圆柱相比,垂直交叉双圆柱绕流上下圆柱的升阻力特性、Strouhal数等都呈现出许多不同的特点,下游圆柱升力振幅小于上游圆柱,阻力平均值大于上游圆柱.本文旨在模拟出这种特定间距下的流场涡结构以及圆柱升阻力特性.     

Kármán涡街中旋涡三维变形的初步研究

本文用涡动力学及LIA方法作为基本理论模型,数值研究了尾迹中孤立涡和细涡丝的三维演化规律。结果表明,圆柱分离尾迹中的旋涡存在三维不稳定性。卡门涡在尾迹平均流场中演化产生三维的类似于马蹄形-勺子形的流场结构。细涡丝在涡辫区的三维演化形成趋向流场最大拉伸变形方向的流向涡结构。     

高升力下二维构型阻力准确测量研究

翼型风洞试验阻力测量常使用尾迹流场测量积分求取阻力的方法,但各积分公式均建立在一定的假设基础上,有一定适用范围。在多段翼型流场N-S方程数值模拟和风洞试验的基础上,研究高升力情况下低速风洞阻力精确测量技术。通过N-S方程数值模拟求解多段翼型绕流场,分析尾迹流场的特点和常规风洞试验阻力计算公式推导时所作假设,提出新的更为准确的型阻计算公式;利用多段翼型绕流的数值模拟结果,积分表面压力和摩擦力求得翼型的气动特性,并利用计算得到的尾迹流场信息按照常规和新提出的风洞试验型阻计算公式计算阻力,将三者进行比较,检验提出的新型阻计算公式的准确性;通过风洞试验检验数值模拟得到的流场特点和新型阻计算公式。研究表明：在高升力条件下,传统型阻计算公式有很大的局限性,必须进行改进;提出的考虑尾迹区流动特点的新型阻计算公式能够得到更准确的阻力值。     

仿鸟柔性扑翼气动特性与能耗的数值研究

建立了适当的三维仿鸟柔性扑翼模型,并以配平重力和平衡阻力为条件,数值计算了它的低雷诺数非定常流场.研究揭示了翼面初始扭转角度、动态俯仰幅度等重要设计参数与飞行性能的关系,表明扑翼平面的初始扭转程度、扑翼柔性材料的选择以及两者之间的合理搭配对扑翼机的成功飞行至关重要.研究分析了仿鸟扑翼的流场涡结构、升力推力产生原理,下扑过程附着上翼面的前缘涡是升力产生的重要机制.对扑翼气动功率的比较分析也发现,人造扑翼机需要的气动功率明显高出同等大小的鸟类,在效率方面尚不及扑翼飞行生物.     

Individual influence of pitching and plunging motions on flow structures over an airfoil during dynamic stall

The individual influence of pitching and plunging motions on flow structures is studied experimentally by changing the phase lag between the geometrical angle of attack and the plunging angle of attack. Five phase lags are chosen as the experimental parameters, while the Strouhal number, the reduced frequency and the Reynolds number are fixed. During the motion of the airfoil, the leading edge vortex, the reattached vortex and the secondary vortex are observed in the flow field. The leading edge vortex is found to be the main flow structure through the proper orthogonal decomposition. The increase of phase lag results in the increase of the leading edge velocity, which strongly influences the leading edge shear layer and the leading edge vortex. The plunging motion contributes to the development of the leading edge shear layer, while the pitching motion is the key reason for instability of the leading edge shear layer. It is also found that a certain increase of phase lag, around 34.15° in this research, can increase the airfoil lift.     

扑动翼型的低雷诺数气动特性分析

通过求解引入拟压缩项的不可压Navier-Stokes方程,数值模拟了绕扑动翼型的低雷诺数非定常流动。针对厚度在4%-12%之间的NACA对称翼型,分析了翼型厚度等参数对扑动翼型气动特性的影响。在低雷诺数条件下,对于纯俯仰运动,随着翼型厚度的减小,平均阻力系数也变小。而对于纯沉浮运动,发现翼型厚度对气动特性的影响和俯仰运动有很大的差别,平均阻力系数随着翼型厚度的减小而变大。通过对沉浮运动一个周期流线图的分析,认为这是翼型前缘涡的影响造成的。由于前缘涡的影响,翼型厚度增加,平均压差阻力系数变小,甚至会出现负值。雷诺数的影响研究表明,随着雷诺数的增加,扑动翼型的阻力系数减小的趋势越缓慢。     

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

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

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.     

不同雷诺数下翼型气动特性及层流分离现象演化

低雷诺数下空气黏性效应突出,翼型表面普遍存在层流分离现象,相比常规雷诺数情况气动特性显著恶化。采用带预处理的Roe方法求解非定常可压缩Navier-Stokes方程的数值模拟技术和低雷诺数低湍流度风洞油流显示试验技术,对FX63-137翼型不同雷诺数下气动特性和流动结构展开深入研究。通过风洞油流显示试验可以清晰获得低雷诺数层流分离流动的两道油流汇集线。数值模拟结果表明其分别为时均化主分离线和二次分离线,两种结果定性定量均吻合较好,证明了本文的研究方法有效可靠;雷诺数从500 000降至20 000,翼型气动特性和层流分离流动结构均发生显著的变化,伴随阻力系数剧增和升力系数剧降,时均化流动结构从附体至出现经典的长层流分离泡,并最终演化为后缘层流分离泡,相应的两种分离泡的非定常流动结构也存在显著差异;对于阻力系数和升力系数而言,存在不同的临界雷诺数,因为导致阻力系数剧增的机理在于经典长层流分离泡的产生使翼型压差阻力大增,而造成升力系数剧降的主要原因在于后缘层流分离泡使得等效翼型后部弯度减小;非定常结果显示正是由于翼型表面漩涡周期性的生成与脱落,才造成了低雷诺数下升力系数的周期性波动。翼型上表面主分离涡即将脱落时,流线在后缘附近再附,升力系数达到峰值;而当流体从下表面向上卷起二次分离涡时,尾部流线大尺度分离,升力系数降至谷值。     

Computations of unsteady transonic flow over airfoil at low Reynolds number

Transonic flow over a thin airfoil at low Reynolds number was studied numerically by directly solving two-dimensional full Navier-Stokes equations through 5th order weighted essentially non-oscillatory(WENO) scheme without using any turbulence model.A series of distinguished unsteady phenomena for a thin 2-D transonic airfoil flow were presented.Due to continuous adverse pressure gradient in the subsonic flow downstream of the sonic line, the unsteady separated boundary layer with main vortex and secondary vortex was developed at the rear of the airfoil.At the trailing edge,the vortex-shedding was characterized by periodical connection of the main vortex and secondary vortex on the other side of the airfoil.The unsteady separation and vortex-shedding occurred with the same period.On the airfoil surface,the average pulse pressure related to the unsteady supersonic region was obviously smaller than that related to the vortex-shedding at the trailing edge.With the attack angle increasing from 0° to 2°, the frequency of vortex-shedding decreases about 4.2%.At last, the turbulence intensity and many second-order statistics in the wake region were investigated.      

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

通过对非定常N-S方程的数值求解,研究了最大厚度为12％的Karman-Trafftz翼型在Re数为1000时的大迎角俯仰振动。其中着重分析了旋涡结构与表面压强分布的关系。数值研究表明,后缘形状、折合频率等对涡结构演化有重要影响。后缘涡顺利地从后缘脱落时,失速涡在上翼面能诱导出较大的吸力。后缘涡在翼面上驻留时,各涡产生复杂的相互干扰,对失速涡在上翼面产生吸力有不利影响。     

Experimental and numerical investigation of threedimensional vortex structures of a pitching airfoil at a transitional Reynolds number

This research examines the vortex behaviors and aerodynamic forces in dynamic stall phenomena at a transitional Reynolds number(Re = 90000) using experimental and numerical approaches.Periodic sinusoidal pitching motion at two different reduced frequencies is used to achieve the dynamic stall of a NACA 0012 airfoil.Several leading edge vortices form and detach in the dynamic stall stage.The flow then quickly transitions to a full separation zone in the stall stage when the angle of attack starts to decrease.There is discrepancy between the phaseaveraged and instantaneous flow field in that the small flow structures increased with angle of attack, which is a characteristic of the flow field at the transitional Reynolds number.The interaction between the streamwise vortices in the three-dimensional numerical results and the leading edge vortex are the main contribution to the turbulent flow.In addition, the leading edge vortex that supplies vortex lift is more stable at higher reduced frequency, which decreases the lift fluctuation in the dynamic stall stage.The leading edge vortex at higher reduced frequency is strong enough to stabilize the flow, even when the airfoil is in the down-stroke phase.