Self-sustained shock oscillations on airfoils at transonic speeds

Self-sustained shock wave oscillations on airfoils at transonic flow conditions are associated with the phenomenon of buffeting. The physical mechanisms of the periodic shock motion are not yet fully understood even though experiments performed over fifty years ago have demonstrated the presence of oscillatory shock waves on the airfoil surfaces at high subsonic speeds. The unsteady pressure fluctuations generated by the low-frequency large-amplitude shock motions are highly undesirable from the structural integrity and aircraft maneuverability point of view. For modern supercritical wing design with thick profiles, the shock-induced fluctuations are particularly severe and methods to reduce the shock wave amplitudes to lower values or even to delay the oscillations to higher Mach numbers or incidence angles will result in expanding the buffet boundary of the airfoil. This review begins with a recapitulation of the classical work on shock-induced bubble separation and trailing edge separation of a turbulent boundary layer. The characteristics of the unsteady pressure fluctuations are used to classify the types of shock-boundary layer interaction. The various modes of shock wave motion for different flow conditions and airfoil configurations are described. The buffet boundaries obtained using the standard trailing edge pressure divergence technique and an alternative approach of measuring the divergence of normal fluctuating forces are compared to show the equivalence. The mechanisms of self-sustained shock oscillations are discussed for symmetrical circular-arc airfoils at zero incidence and for supercritical airfoils at high incidence angles with fully separated flows. The properties of disturbances in the wake are examined from linear stability analysis of two-dimensional compressible flows. The advances in high-speed computing make predictions of buffeting flows possible. Navier–Stokes solvers and approximate boundary layer-inviscid flow interaction methods are shown to give good correlation of frequencies and other unsteady flow characteristics with experiments. Finally, passive and active methods of shock oscillation control show promising results in delaying buffet onset to higher Mach numbers or incidence angles, thus enhancing the transonic performance of airfoils.     

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.     

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.     

超临界翼型尾缘噪声影响因素研究

使用大涡模拟(LES)和FW-H相结合的方法对超临界翼型尾缘噪声进行了研究。针对超临界翼型后缘的不同设计参数(后缘厚度、后缘角和几何形状)对尾缘噪声的影响,使用ANSYS FLUENT的LES湍流模型计算声源,采用FW-H积分方法求解远场噪声总声压级。首先完成二维非定常流圆柱绕流的流场及其噪声的验证计算,计算结果与实验值符合,证明了求解器设置和网格生成的合理性。然后基于此正确的求解器设置和网格生成,使用LES/FWH对比了典型超临界翼型的不同后缘设计参数对远场总声压级的影响,所得结论对低噪声超临界翼型优化设计有参考意义,同时为进一步的噪声控制优化设计提供了基础。     

激波控制鼓包提高翼型跨声速抖振边界

翼型抖振边界是仅次于升阻比的一项重要气动指标.采用定常雷诺平均Navier-Stokes方程,以升力线斜率平缓及激波位置振荡作为基本判据确定了RAE2822翼型在指定跨声速来流条件下的抖振边界.通过大量计算流体力学(CFD)验证,针对RAE2822翼型设计了一种特定外形的激波控制鼓包并确定了其具体安装位置.该激波控制鼓...     

Investigation of upstream moving pressure waves on a supercritical airfoil

The present work is a preliminary experimental and numerical investigation of wave processes taking place in the flow field on a supercritical airfoil in a defined Mach and Reynolds number range. These waves which originate near the trailing edge propagate upstream in the subsonic region and become apparently weak near the leading edge. Time-resolved pressure measurements performed, reveal the unsteady behaviour of these waves. The frequencies measured are in the order of kHz. The wave structures can also be seen on corresponding time-resolved shadowgraphs. Not yet fully understood, however, is the mechanism which generates the waves as well as their interaction with the boundary layer. Early numerical investigations show that the wave formation is coupled with a vortex formation in the boundary layer.     

锯齿尾缘叶片气动特性和绕流流场的数值研究

以基于NACA 0018翼型的锯齿尾缘仿生叶片为研究对象，采用大涡模拟的方法研究锯齿相对齿宽与相对齿高对锯齿尾缘叶片的气动特性和非定常绕流流场的影响规律和机制.研究表明，尾缘锯齿参数对叶片气动性能的影响是复杂的非线性过程，在一定来流攻角范围内能提高升阻比，但失速提前.如在9.4°~14.8°来流攻角范围内，不同相对齿宽系列叶片的升阻比高于原始叶片，升阻比与锯齿相对齿宽之间没有线性关系.研究还表明，锯齿尾缘能延迟边界层分离，加速尾迹的流动掺混和能量扩散，改变非定常涡结构和涡脱落频率.相对齿高的变化对非定常流动特性的影响更为显著.尾缘锯齿诱导的二次湍流射流和吸力面侧反向涡对改变了原始叶片的绕翼环量，进而影响锯齿尾缘叶片的气动特性和绕流流场特性.      

高压涡轮转子间隙泄漏流动的非定常特征研究

利用数值方法求解三维非定常雷诺平均Navier-Stokes方程模拟某跨声速高压涡轮流场,研究了某跨声速高压涡轮流场的非定常特征,通过详细分析动静干涉对间隙泄漏流动的影响,进一步明确了泄漏流周期性变化的规律和成因.研究发现:静子尾缘燕尾波的外侧分支外尾波是间隙内部流动结构变化的主要原因,间隙泄漏涡的周期性变化则受外尾波和尾迹的共同影响.外尾波深入转子通道内部周期性经过间隙,在间隙前缘附近产生很强的逆压梯度,使间隙前部流动方向明显改变而产生大范围分离.外尾波导致间隙泄漏流量明显增加并周期性震荡.在静子尾迹和外尾波共同作用下,泄漏涡强度出现波动且涡的位置前后移动,使泄漏涡呈现明显的非定常性.      

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

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

Mini-TED跨声速气动特性及其增升机理

通过求解二维可压Navier-Stokes方程,研究了NACA0012翼型加装微型后缘增升装置（mini-TED）后的跨声速流场特性,与Gurney flap （GF）对比分析了几何参数对mini-TED后方涡系及翼型气动特性的影响.将mini-TED的几何细节参数定义为弦向长度和有效高度,两者方向正交.在相同迎角下仅改变mini-TED的弦向长度,后缘涡系结构虽发生变化,但翼型气动力几乎没有影响;反之仅改变有效高度则后缘涡系和翼型气动力系数同时发生明显改变,且与同等高度下的GF气动系数相近.结果表明：有效高度是影响翼型气动特性的决定因素.有效高度改变了mini-TED后涡系的发生范围,而相对于整个翼型绕流,后缘涡系的大小是影响翼型流场最重要的因素,而涡系的微观结构和形态的改变影响相对很小.加装mini-TED后上表面激波位置后移、下表面激波强度削弱,从而翼型表面压力分布特性发生了改变.随有效高度增大,mini-TED诱导的涡系发生区域随之增大,引流作用增强,翼型升力系数、阻力系数和低头力矩系数提高,同时相同迎角下翼型的升阻比明显提高.     

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.      

Dynamic stall on a supercritical airfoil

For the present investigations of dynamic stall a supercritical airfoil was chosen. This new airfoil designed by DLR will be used in dynamic stall control research activities (project ADASYS) planned for the near future: the leading edge portion of the airfoil will be drooped down dynamically to improve dynamic stall characteristics on the retreating side during blade motion. The optimised transonic properties of the airfoil, i.e., reduction of shock strength over a Mach number range will improve in addition the performance of the advancing rotor blade. Dynamic stall experiments on the rigid supercritical airfoil have first been carried out in the DNW-TWG transonic wind tunnel with a 1 m × 1 m cross section of the test section and adaptive top and bottom – walls. This tunnel has the advantage to cover the speed range of both retreating and advancing blade. Emphasis has been placed on unsteady pressure measurements along the adaptive walls simultaneously with the unsteady pressure measurements on the pitching model. In addition to the experiments corresponding numerical simulations with a RANS-code have been carried out and their results are compared with the experimental data. Of main concern are the influence of laminar-turbulent boundary-layer transition as well as wind-tunnel-wall interference effects on the unsteady results.     

重气体介质中超临界翼型跨声速流动特性

采用Peng-Robinson非理想气体状态方程模拟重气体介质的热力学特性,并与雷诺平均Navier-Stokes方程结合,形成封闭的重气体介质流动模型。针对超临界翼型流动问题,利用LU-SGS（lower-upper symmetric Gauss-Seidel）隐式时间推进格式和有限体积法,分别求解空气介质和重气体介质下的流动特性。数值模拟结果表明:在跨声速条件下重气体介质中超临界翼型的升阻力增大、超声速区域表面负压增加、边界层位移厚度减小、激波后移、表面摩擦阻力明显增大、后缘流动分离推迟。该研究为后续重气体介质中飞行器颤振特性研究及修正方法的发展提供了基础支持。      

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

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

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

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

对转压气机外伸激波对叶顶泄漏流影响

摘要:为了揭示对转压气机下游转子外伸激波对上游转子泄漏流的影响规律,针对上游转子叶顶间隙分别为0.2、0.5、0.8 mm的对转压气机开展了非定常数值模拟研究。研究发现:受下游转子外伸激波掠扫影响,上游转子尾缘附近压力面会形成弱压缩波,且随上游转子泄漏流增强而逐渐减弱;而该外伸激波在上游转子尾缘附近吸力面,会形成与型线切向相垂直的较强压缩波,且其位置基本不受叶顶间隙大小影响;外伸激波使上游转子尾缘附近吸、压力面压差增大,叶顶泄漏流增强,进而导致其损失增大;随着叶顶间隙增大,上游转子叶尖区弦长前半段压力波动的频率,由通道激波转为叶顶泄漏流主导,且呈现减小的趋势,而弦长后半段压力波动的频率主要由外伸激波主导,且基本不变。      

低雷诺数NACA0012平面叶栅流场直接数值模拟

采用具有7阶精度的weighted essentially non-oscillatory(WENO)差分格式,直接求解可压缩二维非定常N-S方程组,研究了NACA0012翼型平面叶栅低雷诺数流动的特征.直接模拟及与文献对比的结果表明:叶栅尾缘涡脱落的形成过程与圆柱绕流涡脱落的形成过程非常相似.平面叶栅尾迹区的2阶统计量与孤立翼型尾迹区的2阶统计量具有相同的分布特征,但前者的强度显著大于后者.周期性的涡脱落不仅在上下翼面形成非定常分离,也使得尾迹区某点的总压发生准周期性的变化.随着栅距的减小,翼型上的平均分离位置向前缘移动;尾迹区某点的总压变化频率及其幅值均显著地增加;而且栅距越小,速度脉动2阶统计量反而越大.      

NPU翼型的气动力分析和改进设计

 在飞行器设计中用计算方法设计超临面翼型已完全取代了选用现成翼型的设计方法。为考察已设计出的NPU翼型是否满足飞行器设计要求我们对其进行了全面气动分析,发现这些翼型尚有不足之处,有必要进行改进设计。     

仿生学翼型尾缘锯齿降噪机理

采用大涡模拟与声类比的方法研究了尾缘锯齿对翼型自噪声的影响。以SD2030翼型为研究对象,设计的尾缘锯齿幅值为10%弦长,周期为4%弦长。模拟了来流速度为31 m/s、0° 攻角下直尾缘翼型与锯齿尾缘翼型的流场,对应的基于弦长的雷诺数约为310 000。详细分析了尾缘锯齿对翼型尾缘湍流流场的影响,并通过FW-H方程计算大涡模拟提取的声源项,得到直尾缘翼型与锯齿尾缘翼型的声场。研究发现,锯齿尾缘可以明显降低翼型中低频范围内的噪声,在4 000 Hz以下,窄带噪声最多可降低约16 dB。但尾缘锯齿对翼型气动性能有着不利影响。进一步研究表明,该状态下翼型噪声主要由层流边界层引起的涡脱落噪声主导,尾缘锯齿可以抑制层流边界层引起的涡脱落现象,降低翼型升力脉动与尾缘附近的表面压力脉动,减弱尾缘处的低频湍流脉动与涡量,并有效降低尾缘附近涡的展向相关性,这些因素的综合作用使得翼型自噪声降低。     

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.