低雷诺数俯仰振荡翼型等离子体流动控制

针对低雷诺数翼型气动性能差的特点,通过介质阻挡放电(dielectric barrier discharge,DBD)等离子体激励控制的方法,提高翼型低雷诺数下的气动特性,改善其流场结构.采用二维准直接数值模拟方法求解非定常不可压Navier-Stokes方程,对具有俯仰运动的NACA0012翼型的低雷诺数流动展开数值...     

等离子体环量控制翼型增升的实验研究

利用等离子体激励器发展了新型的环量增升技术,并对二维NACA0012翼型绕流实施控制.由于NACA0012翼型为尖后缘构型,环量增升装置由2个非对称型介质阻挡放电等离子体激励器构成.一个等离子体激励器贴附于翼型吸力面靠近后缘处,其诱导的壁面射流沿来流方向指向下游;另一个等离子体激励器贴附于翼型压力面靠近后缘处,其诱导的壁面射流与来流方向相反指向上游.在风洞中通过时间解析二维PIV系统对翼型绕流流场进行了测量,基于翼型弦长的雷诺数Re=20000.结果表明在等离子体激励器的控制下,翼型压力面靠近后缘处可以形成一个定常回流区,从而起到虚拟气动外形的作用,因此翼型吸力面的流场得到加速,压力面的流场得到减速,使得翼型压力面的吸力以及压力面的压力都得到增加,进而增加了翼型的环量.风洞天平测力实验进一步验证了该环量增升技术的有效性.在整个攻角范围内,施加控制的翼型的升力系数相比没有控制的工况有明显的提高.     

EXPERIMENTAL INVESTIGATION ON LIFT INCREMENT OF A PLASMA CIRCULATION CONTROL AIRFOIL

利用等离子体激励器发展了新型的环量增升技术,并对二维NACA0012翼型绕流实施控制。由于NACA0012翼型为尖后缘构型,环量增升装置由2个非对称型介质阻挡放电等离子体激励器构成。一个等离子体激励器贴附于翼型吸力面靠近后缘处,其诱导的壁面射流沿来流方向指向下游;另一个等离子体激励器贴附于翼型压力面靠近后缘处,其诱导的壁面射流与来流方向相反指向上游。在风洞中通过时间解析二维PIV系统对翼型绕流流场进行了测量,基于翼型弦长的雷诺数Re=20 000。结果表明在等离子体激励器的控制下,翼型压力面靠近后缘处可以形成一个定常回流区,从而起到虚拟气动外形的作用,因此翼型吸力面的流场得到加速,压力面的流场得到减速,使得翼型压力面的吸力以及压力面的压力都得到增加,进而增加了翼型的环量。风洞天平测力实验进一步验证了该环量增升技术的有效性。在整个攻角范围内,施加控制的翼型的升力系数相比没有控制的工况有明显的提高。     

新型等离子体激励器对翼型的减阻效果

论述了在西北工业大学低湍流度风洞中采用新型等离子激励器对NACA0015翼型的减阻实验.实验风速为35m/s,攻角范围取0°~20°.并参照压力分布和总压分布实验结果对减阻效果进行了对比分析.本文还进行了有关等离子体激励抑制翼型流动分离的数值模拟研究,基于等离子体激励器的简化模型将体积力以源项方式引入到N-S方程中求解,得到激励器工作时的流场分布.结果表明在新型等离子体激励器开启后:在小攻角范围内,尾耙的总压分布曲线与坐标轴的纵轴(尾耙高度轴)所围面积变化不大;当攻角≥12°时,尾耙的总压分布曲线与坐标轴的纵轴(尾耙高度轴)所围面积明显减小.从而说明该新型等离子体激励器能够有效地减少翼型的阻力.     

基于后缘小翼的翼型反流动态失速主动控制试验研究

针对直升机旋翼反流区因反流动态失速导致的非定常载荷、阻力激增以及负升力等问题,开展了基于后缘小翼的翼型反流动态失速主动控制试验研究.采用动态压力测量结合翼型表面压力积分的方法,重点分析了后缘小翼不同的振荡相位差、幅值和减缩频率对反流动态失速控制的影响规律,对比了后缘小翼动态偏转和固定偏转的差异,试验雷诺数Re=3.5×10⁵.结果表明,当后缘小翼与翼型以相同的频率正弦振荡运动,且二者的相位差为0°时,能改善反流动态失速过程中钝几何前缘的流动分离,并在反流状态下实现了翼型负升力系数下降21.2%,阻力系数下降37.5%,俯仰力矩系数迟滞环面积下降44.6%的控制效果;动态偏转的后缘小翼对翼型反流动态失速的控制效果随后缘小翼振荡幅值的增加而增加,但进一步增加振荡幅值对于控制效果的提升有限;当减缩频率增加时,动态偏转的后缘小翼对反流状态下翼型阻力的控制效果会更加明显;后缘小翼的动态偏转与固定偏转都能有效改善翼型在反流中的动态气动性能,但是动态偏转对于不同翼型迎角的适应能力优于固定偏转,并取得了更好的非定常载荷控制以及更好的阻力和负升力改善效果.     

圆形坚冰影响翼型气动性能的数值分析

本文用FLUENT软件模拟了结冰后NACA 0012翼型周围流场的变化,并与结冰前NACA 0012翼型的气动性能进行了对比.工作中首先以未结冰的NACA 0012翼型(干净翼型)为标准模型进行了数值验证计算,再以经过检验的方法计算结冰模型,并与结冰风洞试验数据进行了对比.本文计算攻角为0°～20°,温度为250.37K,雷诺数为2,400,000,冰型为圆形坚冰.通过对比升力阻力性能,发现与干净翼型相比,结冰翼型的最大升力系数大约减少了50%,阻力系数增加了约65%,失速攻角降低了4°.结冰后翼型提前失速是造成气动性能恶化的主要原因.     

风力机翼型动态失速等离子体流动控制数值研究

针对动态失速引起的风力机翼型气动性能恶化的问题,本文基于动网格和滑移网格技术, 开展了大涡模拟数值计算研究,探索了非定常脉冲等离子体的动态流动控制机理. 结果表明,等离子体气动激励能够有效控制翼型动态失速, 改善平均和瞬态气动力,减小力矩负峰值和迟滞环面积. 压力分布在等离子体施加范围内出现了负压"凸起",上翼面吸力峰值明显增大.脉冲频率和占空比这两个非定常控制参数对流动控制影响显著,无因次脉冲频率为1.5时等离子体控制效果较好,占空比为0.8时即可接近连续工作模式下的气动收益. 翼型深失速状态,等离子体促使流动分离位置明显向后缘移动, 抵抗了大尺度动态失速涡的发生,分离涡结构破碎耗散、重新附着, 涡流影响范围减小; 浅失速状态,等离子体激励具有较强的剪切层操纵能力, 诱导了翼型边界层提前转捩,促进了与主流的动量掺混. 等离子体气动激励诱导出前缘附近贴体翼面"涡簇",起到了虚拟气动外形的作用.不同尺度、频域的动态涡结构与等离子体气动激励的非线性、强耦合作用导致了气动力/力矩的谐波振荡.      

NUMERICAL STUDY ON DYNAMIC STALL FLOW CONTROL FOR WIND TURBINE AIRFOIL USING PLASMA ACTUATOR 1)

针对动态失速引起的风力机翼型气动性能恶化的问题,本文基于动网格和滑移网格技术, 开展了大涡模拟数值计算研究,探索了非定常脉冲等离子体的动态流动控制机理. 结果表明,等离子体气动激励能够有效控制翼型动态失速, 改善平均和瞬态气动力,减小力矩负峰值和迟滞环面积. 压力分布在等离子体施加范围内出现了负压"凸起",上翼面吸力峰值明显增大.脉冲频率和占空比这两个非定常控制参数对流动控制影响显著,无因次脉冲频率为1.5时等离子体控制效果较好,占空比为0.8时即可接近连续工作模式下的气动收益. 翼型深失速状态,等离子体促使流动分离位置明显向后缘移动, 抵抗了大尺度动态失速涡的发生,分离涡结构破碎耗散、重新附着, 涡流影响范围减小; 浅失速状态,等离子体激励具有较强的剪切层操纵能力, 诱导了翼型边界层提前转捩,促进了与主流的动量掺混. 等离子体气动激励诱导出前缘附近贴体翼面"涡簇",起到了虚拟气动外形的作用.不同尺度、频域的动态涡结构与等离子体气动激励的非线性、强耦合作用导致了气动力/力矩的谐波振荡.     

可渗透壁对翼型声涡相互作用下流场及声场的影响

由仿生学原理构建的可渗透翼型对湍流气动噪声抑制作用已展现良好的应用前景。对NACA 0012可渗透翼型和实体翼型进行了数值计算,得到了声涡相互作用下气动噪声声场和流场,分析了可渗透壁对翼型流场和声场的影响。研究表明,相对实体翼型,可渗透壁通过减小声源强度降低了主纯音噪声声压级幅值和远场总声压级,消除了高阶离散纯音,但对噪声的指向性没有较大改变。进一步的流场分析表明,可渗透壁对翼型气动性能影响不大的情况下能够降低边界层扰动和翼型后缘大尺度涡旋强度,并推迟分离泡转捩和再附位置。     

基于涡粘性模型的翼型湍流流场熵产率计算

利用有限体积法实现了基于非正交同位网格的SIMPLE算法。基于熵分析方法,采用涡粘性模型求解湍流熵产方程,系统研究了湍流模型对二维翼型绕流流场熵产率的影响。通过计算NACA0012翼型在来流雷诺数为2.88×106时,0°攻角~16.5°攻角范围内的翼型表面压力系数分布和升阻力特性,验证了算法及程序的正确性。结果表明,选择不同湍流模型时,翼型流场熵产的计算结果存在差异,湍流耗散是引起流场熵产的主要原因;翼型流场的熵产主要发生在翼型前缘区、壁面边界层和翼型尾流区域,流场熵产率与翼型阻力系数线性相关;当产生分离涡时,粘性耗散引起的熵产下降。     

Flow control over a NACA 0012 airfoil using dielectric-barrier-discharge plasma actuator with a Gurney flap

Flow control study of a NACA 0012 airfoil with a Gurney flap was carried out in a wind tunnel, where it was demonstrated that a dielectric-barrier-discharge (DBD) plasma actuator attached to the flap could increase the lift further, but with a small drag penalty. Time-resolved PIV measurements of the near-wake region indicated that the plasma forcing shifted the wake downwards, reducing its recirculation length. Analysis of wake vortex dynamics suggested that the plasma actuator initially amplified the lower wake shear layer by adding momentum along the downstream surface of the Gurney flap. This enhanced mutual entrainment between the upper and lower wake vortices, leading to an increase in lift on the airfoil.     

基于LBM-LES方法翼型纯音噪声数值研究

基于LBM-LES方法,对中低雷诺数下的NACA0012翼型纯音噪声进行了直接计算,研究了不同迎角和雷诺数对纯音噪声的影响。计算结果表明,翼型的声源主要位于翼型的分离区和后缘处,在不同迎角和雷诺数下的声辐射特征均具有偶极子声场的特点;迎角的增大将引起较大的旋涡尺度和湍流强度,吸力面声源区域前移。声压级频谱分析表明,随着迎角的增大,纯音噪声逐渐消失,噪声谱最终呈现宽频特征;随着雷诺数的增大,后缘压力脉动增大。声压级频谱中,主频频率随着雷诺数的增大而增大,且符合Paterson公式的幂律关系。此外,声压级频谱特性随着雷诺数的增大表现出由离散特性向宽频特性转变的趋势。     

Effect of compressibility on plasma-based transition control for a wing with leading-edge excrescence

Compressibility effects were numerically investigated for use of plasma-based flow control, which was applied to delay transition generated by excrescence on the leading edge of a wing. The wing airfoil section incorporates a geometry that is representative of modern reconnaissance air vehicles, and has an appreciable region of laminar flow at design conditions. Modification of the leading edge can be caused by the accumulation of debris, insect impacts, microscopic ice crystal formation, damage, or structural fatigue, resulting in premature transition and an increase in drag. A dielectric barrier discharge (DBD) plasma actuator, located downstream of the excrescence, was employed to delay transition, mitigate the effects of turbulence, decrease drag, and increase energy efficiency. Solutions were obtained for several Mach numbers, up to the transonic range. The effect of compressibility on transitional behaviour was explored, and the effectiveness of plasma-based control to delay transition with increasing Mach number was determined.     

Investigation of steady plasma actuation effect on aerodynamic coefficients of oscillating airfoil at low Reynolds number

In this work, numerical study of two dimensional laminar incompressible flow around an oscillating NACA0012 airfoil is proceeded using the open source code Open FOAM. Oscillatory motion types including pitching and flapping are considered. Reynolds number for these motions is assumed to be 12000 and effects of these motions and also different unsteady parameters such as amplitude and reduced frequency on aerodynamic coefficients are studied. For flow control on airfoil, dielectric barrier discharge plasma actuator is used in two different positions on airfoil and its effect is compared for the two types of considered oscillating motions. It is observed that in pitching motion, imposing plasma leads to an improvement in aerodynamic coefficients, but it does not have any positive effect on flapping motion.Also, for the amplitudes and frequencies investigated in this paper, the trailing edge plasma had a more desirable effect than other positions.     

LES of transient flows controlled by DBD plasma actuator over a stalled airfoil

Large-eddy simulations (LES) are employed to understand the flow field over a NACA 0015 airfoil controlled by a dielectric barrier discharge (DBD) plasma actuator. The Suzen body force model is utilised to introduce the effect of the DBD plasma actuator. The Reynolds number is fixed at 63,000. Transient processes arising due to non-dimensional excitation frequencies of one and six are discussed. The time required to establish flow authority is between four and six characteristic times, independent of the excitation frequency. If the separation is suppressed, the initial flow conditions do not affect the quasi-steady state, and the lift coefficient of the higher frequency case converges very quickly. The transient states can be categorised into following three stages: (1) the lift and drag decreasing stage, (2) the lift recovery stage, and (3) the lift and drag converging stage. The development of vortices and their influence on control is delineated. The simulations show that in the initial transient state, separation of flow suppression is closely related to the development spanwise vortices while during the later, quasi-steady state, three-dimensional vortices become more important.     

Duty cycle and directional jet effects of a plasma actuator on the flow control around a NACA0015 airfoil

This paper reports on the effects of a series of fluid-dynamic dielectric barrier discharge plasma actuators on a NACA0015 airfoil at high angle of attack. A set of jet actuators able to produce plasma jets with different directions (vectoring effect) and operated at different on/off duty cycle frequencies are used. The experiments are performed in a wind tunnel facility. The vectorized jet and the transient of the flow induced by unsteady duty cycle operation of each actuator are examined and the effectiveness of the actuator to recover stall condition in the range of Reynolds numbers between 1.0 × 10⁵ and 5.0 × 10⁵ (based on airfoil chord), is investigated. The actuator placed on the leading edge of the airfoil presents the most effective stall recovery. No significant effects can be observed for different orientations of the jet. An increase of the stall recovery is detected when the actuator is operated in unsteady operation mode. Moreover, the frequency of the on/off duty cycle that maximizes the stall recovery is found to be a function of the free stream velocity. This frequency seems to scale with the boundary layer thickness at the position of the actuator. A lift coefficient increase at low free stream velocities appears to linearly depend on the supply voltage.     

High fidelity numerical simulation of airfoil thickness and kinematics effects on flapping airfoil propulsion

High-fidelity numerical simulations with the spectral difference (SD) method are carried out to investigate the unsteady flow over a series of oscillating NACA 4-digit airfoils. Airfoil thickness and kinematics effects on the flapping airfoil propulsion are highlighted. It is confirmed that the aerodynamic performance of airfoils with different thickness can be very different under the same kinematics. Distinct evolutionary patterns of vortical structures are analyzed to unveil the underlying flow physics behind the diverse flow phenomena associated with different airfoil thickness and kinematics and reveal the synthetic effects of airfoil thickness and kinematics on the propulsive performance. Thickness effects at various reduced frequencies and Strouhal numbers for the same chord length based Reynolds number (=1200) are then discussed in detail. It is found that at relatively small Strouhal number (=0.3), for all types of airfoils with the combined pitching and plunging motion (pitch angle 20°, the pitch axis located at one third of chord length from the leading edge, pitch leading plunge by 75°), low reduced frequency (=1) is conducive for both the thrust production and propulsive efficiency. Moreover, relatively thin airfoils (e.g. NACA0006) can generate larger thrust and maintain higher propulsive efficiency than thick airfoils (e.g. NACA0030). However, with the same kinematics but at relatively large Strouhal number (=0.45), it is found that airfoils with different thickness exhibit diverse trend on thrust production and propulsive efficiency, especially at large reduced frequency (=3.5). Results on effects of airfoil thickness based Reynolds numbers indicate that relative thin airfoils show superior propulsion performance in the tested Reynolds number range. The evolution of leading edge vortices and the interaction between the leading and trailing edge vortices play key roles in flapping airfoil propulsive performance.     

High-lift airfoil trailing edge separation control using a single dielectric barrier discharge plasma actuator

Control of flow separation from the deflected flap of a high-lift airfoil up to Reynolds numbers of 240,000 (15 m/s) is explored using a single dielectric barrier discharge (DBD) plasma actuator near the flap shoulder. Results show that the plasma discharge can increase or reduce the size of the time-averaged separated region over the flap depending on the frequency of actuation. High-frequency actuation, referred to here as quasi-steady forcing, slightly delays separation while lengthening and flattening the separated region without drastically increasing the measured lift. The actuator is found to be most effective for increasing lift when operated in an unsteady fashion at the natural oscillation frequency of the trailing edge flow field. Results indicate that the primary control mechanism in this configuration is an enhancement of the natural vortex shedding that promotes further momentum transfer between the freestream and separated region. Based on these results, different modulation waveforms for creating unsteady DBD plasma-induced flows are investigated in an effort to improve control authority. Subsequent measurements show that modulation using duty cycles of 50–70% generates stronger velocity perturbations than sinusoidal modulation in quiescent conditions at the expense of an increased power requirement. Investigation of these modulation waveforms for trailing edge separation control similarly shows that additional increases in lift can be obtained. The dependence of these results on the actuator carrier and modulation frequencies is discussed in detail.     

Transient phenomena in separation control over a NACA 0015 airfoil

We present the transient phenomena occurring during the impulsive control of flow separation over a NACA0015 airfoil at an incidence angle of 11° and a chord Reynolds number of 1 million. Actuation is performed via pneumatic vortex generators, impulsively activated in order to analyze the transient phenomena corresponding to the attachment process and, conversely, to transient re-separation occurring when the actuators are switched off. Measurements are performed using a linear array of unsteady pressure transducers and a single traversing crosswire. The pressure transducers are positioned in the separated region of the airfoil, which extends ∼ 0.3c upstream of the trailing edge at the above flow condition. To control the flow, the angled fluidic vortex generators are positioned in a single spanwise array located 0.3c downstream of the leading edge of the airfoil. We establish a statistical relationship between pressure and velocity signals during both the uncontrolled steady state and the transient processes of attachment and separation. The unsteady behavior of the attachment process is also qualitatively analyzed via a 0.3 million Reynold number visualizations. The emission of a “starting vortex” is evidenced. This corresponds to a transient increase of drag.