Numerical investigation of aerodynamic flow actuation produced by surface plasma actuator on 2D oscillating airfoil

Numerical simulation of unsteady flow control over an oscillating NACA0012 airfoil is investigated. Flow actuation of a turbulent flow over the airfoil is provided by low current DC surface glow discharge plasma actuator which is analytically modeled as an ion pressure force produced in the cathode sheath region. The modeled plasma actuator has an induced pressure force of about 2 k Pa under a typical experiment condition and is placed on the airfoil surface at 0% chord length and/or at 10% chord length. The plasma actuator at deep-stall angles(from 5° to 25°) is able to slightly delay a dynamic stall and to weaken a pressure fluctuation in down-stroke motion. As a result, the wake region is reduced. The actuation effect varies with different plasma pulse frequencies, actuator locations and reduced frequencies. A lift coefficient can increase up to 70% by a selective operation of the plasma actuator with various plasma frequencies and locations as the angle of attack changes. Active flow control which is a key advantageous feature of the plasma actuator reveals that a dynamic stall phenomenon can be controlled by the surface plasma actuator with less power consumption if a careful control scheme of the plasma actuator is employed with the optimized plasma pulse frequency and actuator location corresponding to a dynamic change in reduced frequency.     

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方程与滚转运动方程耦合计算方法,比较研究了不同后掠角的双三角翼和翼身组合体的滚转运动特性,分析了机翼前缘后掠角及细长机身对非定常滚转力矩时滞环、动态流场结构和物面瞬时压力分布的影响。研究结果表明：主翼迎风面上的融合涡能量在80°/60°双三角翼上耗散较小,而在76°/40°双三角翼上耗散严重,这是造成两模型滚转力矩稳定性与时滞特性差异的主要因素;机身对气流的扰动作用,大幅增强了滚转力矩的线性分量;机身对气流的上洗作用,增强了边条涡与融合涡吸力及其时滞性,同时加剧了主翼背风面的两涡干扰;大滚转角时机身对横流流动的干扰,使得主翼背风面压力分布的时滞差异显著增加。该研究结果有助于认识后掠角与细长机身影响双三角翼滚转运动特性的物理机理。     

Experiments on unsteady vortex flowfield of typical rotor airfoils under dynamic stall conditions

A new experiment for airfoil dynamic stall is conducted by employing the advanced particle image velocimetry(PIV) technology in an open-return wind tunnel. The aim of this experimental investigation is to demonstrate the influences of different motion parameters on the convection velocity, position and strength of leading edge vortex(LEV) of airfoil under different dynamic stall conditions. Two different typical rotor airfoils, OA209 and SC1095, are measured at different free stream velocities, oscillation frequencies, and angles of attack. It is demonstrated by the measured data that the airfoil with larger leading edge radius could notably decrease the strength of LEV. The angle of attack(Ao A) of airfoil can obviously influence the dynamic stall characteristics of airfoil,and the LEV would be effectively inhibited by decreasing the mean pitch angle. In addition, the convection velocity of LEV is estimated in this measurement, and the results demonstrate that the influence of airfoil shape on convection velocity of LEV is limited, but the convection velocity of LEV would be increased by enlarging the oscillation frequency. Meanwhile, the convection velocity of LEV is a time variant value, and this value would increase as the LEV convects to the trailing edge of airfoil.     

基于后缘小翼的旋翼翼型动态失速控制分析

针对后缘小翼(TEF)的典型运动参数对旋翼翼型动态失速特性的控制进行了研究。发展了一套适用于带有后缘小翼控制的旋翼翼型非定常流动特性模拟的高效、高精度CFD方法。通过求解Poisson方程生成围绕旋翼翼型的黏性贴体和正交网格,为保证后缘小翼附近的网格生成质量,建立了基于翼型点重构的方法来描述后缘小翼的偏转运动;为克服变形网格方法可能导致网格畸变的不足,发展了一套适用于带有后缘小翼控制的旋翼翼型运动嵌套网格方法。基于非定常雷诺平均Navier-Stokes(URANS)方程、双时间法、Spalart-Allmaras(S-A)湍流模型和Roe-Monotone Upwind-centered Scheme for Conservation Laws(Roe-MUSCL)插值格式,发展了旋翼翼型非定常气动特性分析的高精度数值方法,并采用Lower-Upper Symmetric Gauss-Seidel(LU-SGS)隐式时间推进方法及并行技术提高计算效率。以有试验结果验证的HH-02翼型和SC1095翼型为算例,精确捕捉了动态失速状态下的气动力迟滞效应,验证了本文方法的有效性。着重针对SC1095旋翼翼型的动态失速状态开展后缘小翼的控制分析,提出了可以体现翼型升力、阻力及力矩综合特性的关系式Po和Pc,揭示了后缘小翼振荡频率、相位差和偏转幅值对动态失速特性影响的规律。研究结果表明:当后缘小翼偏转的相对运动频率为1.0,且小翼运动规律与翼型振荡规律之间的相位差为0°时,后缘小翼能够更好地抑制翼型动态失速现象;在此状态下,当偏转幅值为10°时,SC1095翼型最大阻力系数和最大力矩系数可以分别降低19%和27%。     

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.     

计入离心力影响的直升机旋翼翼型结冰数值模拟

建立了一套计入离心力影响的直升机旋翼翼型结冰的数值模拟方法.首先生成围绕翼型的贴体正交网格,然后用Navier-Stokes（N-S）方程求解黏性绕流流场.在此基础上,利用拉格朗日法建立水滴运动方程.其中,为提高计算效率,提出了结合位移矢量的水滴所处单元寻觅方法.最后,结合桨叶工作特点,发展了一种计入离心力影响的三维结冰模型.通过与桨叶结冰实验的对比,验证了本文结冰预测方法的可靠性.对比常规结冰模型,桨叶结冰量减少22.3%;若考虑桨叶的挥舞运动影响,桨叶结冰量进一步减少,表明了离心力及桨叶运动在结冰数值模拟中的重要性.通过不同剖面间的结冰量和冰形对比,分析并获得了桨叶结冰特征.结果表明离心力的影响程度随径向位移的增加而增加,下翼面结冰量随挥舞角的增加而减少.      

俯仰振荡翼型推阻力转变滞后机制数值研究

扑翼产生的反卡门涡街被认为是一种推力型尾迹,但已有研究指出,随着斯特劳哈尔数(St)增大,低雷诺数下俯仰振荡翼型的净推力产生明显滞后于反卡门涡街的出现。为探究该现象背后的物理机制,对NACA0012翼型在雷诺数1 000条件下作简谐俯仰运动的流场进行了数值模拟。采用翼型表面积分方法和基于有限控制体的气动力估计方法分别研究了翼面分布力特性和尾迹流场特性变化对阻力-推力转变临界点的影响。翼面分布力积分结果表明,当翼型振荡参数进入对应反卡门涡街尾迹形态区域时,翼面压力分布产生的推力分量无法克服剪切层的黏滞阻力,是造成俯仰翼型推力产生滞后于反卡门涡街出现的主要原因。对尾迹流场及相应的推阻力特性变化的分析表明,尽管反卡门涡街会产生"喷流效应",但在St较小时,其产生的动量诱导推力无法克服反卡门涡自身的涡致阻力和尾迹流场压力损失引入的附加阻力,因此即使存在反卡门涡街也不能产生净推力,进而从流场分析的角度进一步解释了这一滞后现象的发生机制。     

Aerodynamic shape optimization for alleviating dynamic stall characteristics of helicopter rotor airfoil

In order to alleviate the dynamic stall effects in helicopter rotor, the sequential quadratic programming(SQP) method is employed to optimize the characteristics of airfoil under dynamic stall conditions based on the SC1095 airfoil. The geometry of airfoil is parameterized by the class-shape-transformation(CST) method, and the C-topology body-fitted mesh is then automatically generated around the airfoil by solving the Poisson equations. Based on the grid generation technology, the unsteady Reynolds-averaged Navier-Stokes(RANS) equations are chosen as the governing equations for predicting airfoil flow field and the highly-efficient implicit scheme of lower–upper symmetric Gauss–Seidel(LU-SGS) is adopted for temporal discretization. To capture the dynamic stall phenomenon of the rotor more accurately, the Spalart–Allmaras turbulence model is employed to close the RANS equations. The optimized airfoil with a larger leading edge radius and camber is obtained. The leading edge vortex and trailing edge separation of the optimized airfoil under unsteady conditions are obviously weakened, and the dynamic stall characteristics of optimized airfoil at different Mach numbers, reduced frequencies and angles of attack are also obviously improved compared with the baseline SC1095 airfoil. It is demonstrated that the optimized method is effective and the optimized airfoil is suitable as the helicopter rotor airfoil.     

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.     

适用于多迎角结冰试验的混合翼型设计

为突破结冰风洞对翼型模型尺寸的限制,提出了一种新的混合翼型设计方法,可使用一套混合翼型模拟原始翼型在不同迎角下的结冰试验,弥补了以往混合翼型只能在单个设计迎角下使用的缺陷.方法采用多段翼的形式设计混合翼型,以多目标迎角等结冰试验条件作为设计输入,优化设计主翼外形和襟翼的位置、偏转角度,利用襟翼位置和偏转角度的变化实现混...     

Experimental investigation of shock-buffet criteria on a pitching airfoil

An experimental investigation of the shock-buffet phenomenon subject to unsteady pitching supercritical airfoil around its quarter chord has been conducted in a transonic wind tunnel. The model was equipped with pressure taps connected to the fast response pressuretransducers. Measurements were conducted at different free-stream Mach number from 0.61 to0.76. The principle goal of this investigation was to experimentally discuss the shock-buffet criterion over a SC(2)-0410 supercritical pitching ...     

变来流速度下旋翼翼型非定常气动特性分析

为分析变来流速度状态下的旋翼翼型气动特性,提出了利用翼型平移来模拟来流速度变化的数值方法.在此方法基础上,采用基于隐式LU-SGS(lower upper symmetric Gauss-Seidal)方法的非定常雷诺平均N-S(Navier-Stokes)(RANS)方程,模拟了SC1095旋翼翼型在定迎角 变来流速度及变迎角 变来流速度状态下的非定常气动特性.通过对比分析发现:翼型在变速度-定迎角状态下会表现出明显的非定常现象,产生了前缘分离涡,气动特性会出现明显的迟滞效应及波动现象,脉动速度越大,非定常效果越明显.并且基准速度越大,翼型气动特性的峰值越大;翼型迎角越大,非定常涡出现的也越早.考虑直升机旋翼翼型实际工作环境,在变速度-动态失速状态下,翼型最大迎角处的气动力会得到一定程度的削弱,在小迎角下的气动力得到一定程度的增强,且脉动速度越大,翼型的非定常特性也越强.      

Numerical study of separation on the trailing edge of a symmetrical airfoil at a low Reynolds number

This study focuses on the trailing-edge separation of a symmetrical airfoil at a low Rey-nolds number. Finite volume method is adopted to solve the unsteady Reynolds-averaged Navier-Stokes (RANS) equation. Flow of the symmetrical airfoil SD8020 at a low Reynolds number has been simulated. Laminar separation bubble in the flow field of the airfoil is observed and process of unsteady bubble burst and vortex shedding from airfoil surfaces is investigated. The time-dependent lift coefficient is characteristic of periodic fluctuations and the lift curve varies nonlinearly with the attack of angle. Laminar separation occurs on both surfaces of airfoil at small angles of attack. With the increase of angle of attack, laminar separation occurs and then reattaches near the trailing edge on the upper surface of airfoil, which forms laminar separation bubble. When the attack of angle reaches certain value, the laminar separation bubble is unstable and produces two kinds of large scale vortex, i.e. primary vortex and secondary vortex. The periodic processes that include secondary vortex production, motion of secondary vortex and vortex shedding cause fluctuation of the lift coefficient. The periodic time varies with attack of angle. The secondary vortex is relatively stronger than the primary vortex, which means its influence is relatively stronger than the primary vortex.     

基于中弧线曲率控制的压气机叶型优化

为提高压气机叶型优化设计水平,基于中弧线曲率控制方法编写了压气机叶片造型程序,将中弧线曲率控制参数作为优化变量,结合粒子群寻优算法对传统可控扩散叶型(CDA)进行了优化研究。结果表明:基于中弧线曲率控制的叶片造型程序能够对CDA叶型进行较好的拟合,拟合叶型的气动性能与设计要求较符。优化叶型在设计点的总压损失降低了约6.34%,优化叶型总压损失随攻角变化较为平缓。在一定攻角范围内,叶型中弧线曲率峰值的前移能够将吸力面马赫数峰值前移,提高叶型吸力面的扩压能力,降低总压损失。在大攻角工况下,改进的中弧线曲率分布能够显著降低叶型总压损失。将中弧线曲率控制参数作为优化变量进行CDA叶型的优化是可行的。     

带有气动及结构非线性的二元机翼颤振分析

研究了翼型在低马赫数条件下的非定常气动特性,从翼型表面气流运动的角度对Leishman-Beddoes(L-B)模型进行了修正,并在此基础上建立了适合低马赫数颤振研究且带有气动及结构非线性的二元机翼气弹系统分析模型.对比低马赫数翼型气动载荷试验结果表明对L-B模型的修正是有效的,且机翼颤振试验结果亦验证了二元机翼气弹分析模型.研究结果表明:二元机翼气弹系统的失速颤振与初始变距角和来流速度密切相关,且耦合的三次非线性变距和浮沉刚度是造成系统呈现准周期运动的主要原因.      

振荡叶栅气动性能实验研究(英文)

为了提供双凸翼型叶栅在所有叶间相位角及实际折合频率下高亚音和跨音定常和振荡气动性能的基本数据,在 NASA路易斯研究中心的跨音振荡叶栅试验设备上进行了一系列试验。为了进行该试验,研制出并使用了非定常气动影响系数法,即叶栅中一次只有一个翼型振荡,测出该振荡翼型与邻近静止翼型上非定常压力的矢量和,藉此就能确定在特定叶间相位角下相当于全部叶片振荡的叶栅的非定常气动性能。      

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

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

仿生正弦前缘对翼面动态失速的影响

动态失速导致叶片气动载荷急剧变化，造成振动载荷激增，桨叶寿命大幅衰减。针对动态失速问题，从座头鲸胸鳍在动态倾转下取得良好的流动特性获得启示，据此模化出仿生正弦前缘翼面（包含3种波峰和2种波长），旨在实现动态失速控制。借助三维非定常数值模拟方法，采用运动网格技术，基于SC1095旋翼翼型，研究了仿生前缘动态失速流动控制机理及运动参数和来流速度的影响。结果表明：正弦前缘大幅度降低俯仰力矩系数峰值和阻力系数峰值；前缘波峰越大、波长越小，阻力系数峰值与俯仰力矩系数峰值的抑制效果越明显，虽然升力系数峰值减小，但其减小量远小于前两者，例如其中一种仿生翼使俯仰力矩系数峰值减小了47.7%，阻力系数峰值减小了36.4%，升力系数峰值减小14.1%；在最大迎角附近，正弦前缘能够缓和失速特性，使载荷变化更为平缓；在高平均迎角、低俯仰频率、低马赫数下，仿生翼动态失速控制效果更强，相比较而言迎角振幅的影响较小。     

一种鼻锥钝化高超声速轴对称进气道流动特性实验

前缘钝化尺度是高超声速进气道设计中的关键参数。针对一种前体锥加弯曲压缩面的高超声速轴对称进气道,选取最大尺度为3.2mm(5%唇缘半径)的几种典型鼻锥钝化半径,在马赫数Ma=6来流,及模型安装攻角为0°、4°、7°的条件下开展鼻锥钝化尺度对进气道流动性能影响的实验研究。采用纹影拍摄及压力测量记录各来流条件下进气道前体流场结构及壁面压强分布,并在无攻角来流条件下利用微型扰流器进行边界层强制转捩研究。结果表明,对无攻角来流而言,即使是尺度高达3.2mm的钝化半径对进气道前体流场结构及壁面静压分布也基本没有影响。此来流条件下,几种不同鼻锥钝化半径的前体压缩面均出现小范围流动分离,而添加扰流器后该分离区均消失。钝化尺度的影响随着攻角的增加而显现,尽管不同鼻锥钝化尺度下迎风面流场及壁面压强分布几乎没有差别,但背风面随钝化尺度增大表现为边界层明显增厚、流动趋于不稳定。其中最大钝化尺度R=3.2mm的构型在4°攻角来流时背风面即出现明显的分离区,而7°攻角来流时背风面更是出现大范围流动分离、进气道背风侧不起动,并导致进气道内部壁面压强显著下降。