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

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

仿生多孔翼型后缘自噪声数值研究

采用大涡模拟(LES)和声比拟(AA)结合的方法,在低马赫数条件下,对SD7003多孔翼型自噪声声场进行计算,研究多孔材料在不同来流迎角下对翼型自噪声远场声压的影响规律,并阐述噪声的影响机理。结果表明：LES与AA相结合的方法能够较为准确地计算多孔翼型的边界层特征以及远场声压,且数值模拟结果与实验结果非常吻合。来流迎角为0°时,多孔翼型整体声压级随着渗透率的减小先减小后增大,渗透率为5×10^-11 m²时多孔材料的降噪效果最佳。多孔翼型整体声压级随着迎角的增大而增大,但多孔材料的降噪效果随着迎角的增大而减小。多孔材料的存在能够有效地削弱翼型表面的大尺度相干结构,使得翼型表面的RMS压强脉动峰值减小,从而抑制翼型远场噪声。     

基于温敏光纤的边界层转捩测量技术研究

利用LUNA公司的ODi SI-B光纤分布式传感系统,开展基于温敏光纤的高精度转捩测量技术研究。试验在NF-3低速风洞中进行。在钢架木质OA309旋翼翼型模型展向中剖面沿弦向开槽埋入电热丝及温敏光纤,用ODi SI-B光纤分布式传感系统测量翼型表面弦向温度分布,按照模型表面温度变化进行边界层转捩判断。通过牛顿冷却公式推导出换热系数比,并根据换热系数比的阶跃来判断转捩位置。结果表明:电热丝与光纤间距的不均匀会导致加热不均匀,进而影响转捩判断;在一定迎角范围下,无论是层流区域还是湍流区域内,温度下降量线重合,换热系数比曲线也重合;转捩区发生从层流线到湍流线跳跃的规律,可以用于边界层转捩判断;可以实现毫米级空间准度的转捩位置测量,得到的转捩位置结果与已有文献的结果一致,转捩判断结果随迎角和雷诺数的变化也符合气动规律。     

弯度翼型大迎角时分离涡流场的多平衡态现象

基于kω的SST两方程湍流模型,在时间域求解雷诺平均Navier-Stokes方程,模拟弯度翼型大迎角时的分离流动。通过给翼型施加一定形式的扰动,重点关注了翼型弯度对大迎角分离涡流场平衡态转移的影响。研究结果表明:与相同厚度20%以上的对称翼型相比,2%弯度的翼型出现分离涡流场平衡态转移的起始迎角变小2°左右,迎角区间变宽约1°;在厚度相对较小的NACA2416翼型上也发现上述分离涡平衡态转移现象。由此说明翼型弯度在一定程度上促使了分离涡平衡态的转移。     

一种层流翼型附面层流动的数值模拟和实验研究

论文通过数值模拟方法对NACA6412层流翼型的低速绕流特性进行了研究,采用修正的BaldwinLomax湍流模型对来流速度为25m/s、8°迎角时翼型的流场、背风区的速度型分布、不同弦向位置附面层厚度的变化规律以及转捩位置进行了模拟,数值模拟结果与风洞实验结果具有一定的相关性.     

翼型厚度对气动弹性特性的影响

结合基于$k$-$\omega$的SST两方程湍流模型,求解雷诺平均Navier-Stokes方程获得定常和非定常气动力,耦合翼型弹性运动方程,在时间域内模拟了不同厚度对称翼型在不同迎角下的气动弹性动态过程,并重点研究了较大迎角下的不同厚度翼型流场特征和气动弹性的性质,研究结果表明：在论文所涉及的参数情况下,对于迎角从零到大迎角范围,翼型颤振临界速度随迎角的变化不是单调的. 翼型颤振临界速度迅速下降的起始迎角比最大升力系数对应的迎角小很多.     

带涡襟翼翼型流场的数值模拟

以数值计算为手段,分析了带涡襟翼的翼型的流场特性,分别对迎角及扰流板偏角对翼型气动性能的影响做了分析。结果表明,在小迎角来流情况下,保持迎角不变,涡襟翼偏转角度越大,升力越小,阻力越大,呈现较好的线性关系。在大迎角情况下,绕翼型的流动发生分离,通过适当控制涡襟翼的偏转角度,能够有效的改善翼型的失速特性,从而达到流动控制的目的,迎角越大,涡襟翼所需偏转的角度越大。     

雷诺数效应对翼型抖振特性的影响

飞行器抖振是一种非线性气动弹性问题,当飞行器进入抖振阶段时,将会对飞行器的性能产生严重影响。而在跨声速条件下,激波附面层相互作用会诱导机翼抖振。本文开展了跨声速条件下翼型抖振特性雷诺数效应的实验研究,揭示了翼型跨声速抖振起始迎角、激波运动前缘边界、频谱特性、抖振频率与雷诺数变化的基本规律。结论如下:雷诺数变化会导致抖振起始边界的改变,对抖振起始迎角下的功率谱密度峰值有明显影响;随着雷诺数的增大,激波运动的前缘后移。雷诺数变化对抖振频率有明显影响,随着马赫数增大,雷诺数效应增强。     

合成射流对失速状态下翼型大分离流动控制的试验研究

为研究低速状态合成射流在抑制翼型气流分离和推迟失速方面的控制机理, 开展了NACA0021 翼型失速特性射流控制的风洞试验研究. 通过系统性的模型测力、翼型瞬态流场粒子图像测速和边界层速度测定的对比试验, 深入探索了合成射流各参数对翼型失速特性控制效果的影响规律. 试验结果表明射流偏角在翼型升力和失速迎角控制方面的效果对射流动量系数较为敏感: 当动量系数较大时, 近切向射流的控制效果更好. 射流动量系数为0.033 时, 偏角30°的射流使得翼型升力系数峰值提高23.56%, 失速迎角增大5°; 而动量系数较小时, 偏角较大的射流能够获得最佳控制效果. 射流动量系数为0.0026 时, 法向射流对翼型最大升力系数控制效果最好(提高9.2%).      

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

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

The flow over a thin airfoil subjected to elevated levels of freestream turbulence at low Reynolds numbers

Micro Air Vehicles (MAVs) can be difficult to control in the outdoor environment as they fly at relatively low speeds and are of low mass, yet exposed to high levels of freestream turbulence present within the Atmospheric Boundary Layer. In order to examine transient flow phenomena, two turbulence conditions of nominally the same longitudinal integral length scale (Lxx/c?=?1) but with significantly different intensities (Ti?=?7.2?% and 12.3?%) were generated within a wind tunnel; time-varying surface pressure measurements, smoke flow visualization, and wake velocity measurements were made on a thin flat plate airfoil. Rapid changes in oncoming flow pitch angle resulted in the shear layer to separate from the leading edge of the airfoil even at lower geometric angles of attack. At higher geometric angles of attack, massive flow separation occurred at the leading edge followed by enhanced roll up of the shear layer. This lead to the formation of large Leading Edge Vortices (LEVs) that advected at a rate much lower than the mean flow speed while imparting high pressure fluctuations over the airfoil. The rate of LEV formation was dependent on the angle of attack until 10° and it was independent of the turbulence properties tested. The fluctuations in surface pressures and consequently aerodynamic loads were considerably limited on the airfoil bottom surface due to the favorable pressure gradient.     

Dynamics of co-rotating vortices in a flow around a bio-inspired corrugated airfoil

Bio-inspired corrugated airfoils show favourable aerodynamic characteristics such as high coefficient of lift and delayed stall at low Reynolds numbers. Two-dimensional (2D) direct numerical simulation has been performed here on a corrugated airfoil at various angles of attack (0°, +5°, -5°) and Reynolds number of 280 to 6700. The objective is to analyse the pressure variation inside the corrugations and correlate it to the vortex movement across the corrugations and the overall aerodynamic characteristics of the corrugated airfoil. The flow characteristics have been examined based on the local Strouhal numbers in the corrugations of the airfoil. It is observed that the pressure variation in each corrugation is the result of vortex merging and separation in the corrugation which plays a major role in changing the flow characteristics. The Strouhal number of the flow is dictated by the most dominant local Strouhal number. The numerical results are further compared with experimental results obtained using particle image velocimetry, and the two set of results are found to match well. These results are significant because they elucidate the effect of corrugation, angle of attack, and Reynolds number on flow over a corrugated airfoil.     

Parametric study of separation and transition characteristics over an airfoil at low Reynolds numbers

Time-resolved surface pressure measurements are used to experimentally investigate characteristics of separation and transition over a NACA 0018 airfoil for the relatively wide range of chord Reynolds numbers from 50,000 to 250,000 and angles of attack from 0° to 21°. The results provide a comprehensive data set of characteristic parameters for separated shear layer development and reveal important dependencies of these quantities on flow conditions. Mean surface pressure measurements are used to explore the variation in separation bubble position, edge velocity in the separated shear layer, and lift coefficients with angle of attack and Reynolds number. Consistent with previous studies, the separation bubble is found to move upstream and decrease in length as the Reynolds number and angle of attack increase. Above a certain angle of attack, the proximity of the separation bubble to the location of the suction peak results in a reduced lift slope compared to that observed at lower angles. Simultaneous measurements of the time-varying component of surface pressure at various spatial locations on the model are used to estimate the frequency of shear layer instability, maximum root-mean-square (RMS) surface pressure, spatial amplification rates of RMS surface pressure, and convection speeds of the pressure fluctuations in the separation bubble. A power-law correlation between the shear layer instability frequency and Reynolds number is shown to provide an order of magnitude estimate of the central frequency of disturbance amplification for various airfoil geometries at low Reynolds numbers. Maximum RMS surface pressures are found to agree with values measured in separation bubbles over geometries other than airfoils, when normalized by the dynamic pressure based on edge velocity. Spatial amplification rates in the separation bubble increase with both Reynolds number and angle of attack, causing the accompanying decrease in separation bubble length. Values of the convection speed of pressure fluctuations in the separated shear layer are measured to be between 35 and 50% of the edge velocity, consistent with predictions of linear stability theory for separated shear layers.     

Finite element computation of turbulent flow past a multi-element airfoil

Flow past multi-element airfoil is studied via two-dimensional numerical simulations. The incompressible Reynolds averaged Navier–Stokes equations, in primitive variables, are solved using a stabilized finite element formulation. The Spalart–Allmaras and Baldwin–Lomax models are employed for turbulence closure. The implementation of the Spalart–Allmaras model is verified by computing flow over a flat plate with a specified trip location. Good agreement is seen between the results obtained with the two models for flow past a NACA 0012 airfoil at 5° angle of attack. Results for the multi-element airfoil, with the two turbulence models, are compared with experiments for various angles of attack. In general, the pressure distribution, from both the models matches quite well with the experimental results. However, at larger angles of attack, the computational results overpredict the suction peak on the slat. The velocity profiles from the Baldwin–Lomax model are, in general, more diffused compared to those from the Spalart–Allmaras model. The agreement between the computed and experimental results is not too good in the flap region for large angles of attack. Both the models are unable to predict the stall; the flow remains attached even for relatively large angles of attack. Consequently, the lift coefficient is over predicted at large α by the computations. Overall, compared to the Baldwin–Lomax model, the predictions from the Spalart–Allmaras model are closer to experimental measurements.     

Influence of pitching angle of incidence on the dynamic stall behavior of a symmetric airfoil

This study presents the influence of pitch angle of an airfoil on its near-field vortex structure as well as the aerodynamic loads during a dynamic stall process. Dynamic stall behavior in a sinusoidally pitching airfoil is usually analyzed at low to medium reduced frequencies and with the maximum angle of attack of the airfoil not exceeding 25°. In this work, we study dynamic stall of a symmetric airfoil at medium to high reduced frequencies even as the maximum angle of attack goes from 25° to 45°. The evolution and growth of the laminar separation bubble, also known as a dynamic stall vortex, at the leading edge and the trailing edge are studied as the pitch cycle goes from the minimum to the maximum angle of attack. The effect of reduced frequencies on the vortex structure as well as the aerodynamic load coefficients is investigated. The reduced frequency is shown to be a bifurcation parameter triggering period doubling behavior. However, the bifurcation pattern is dependent on the variation of the pitch angle of incidence of the airfoil.     

Effect of flow regime change from subsonic to transonic on the air loads of an oscillating airfoil

In this research, the effect of flow regime change from subsonic to transonic on the air loads of a pitching NACA0012 airfoil is investigated. To do this, the effect of change in flow regime on the lift and pitching moment coefficients hysteresis cycles is studied. The harmonic balance approach is utilized for numerical calculation due to its low computational time. Verifications are also made with previous works and good agreements are observed. The assessment of flow regime change on the aforementioned hysteresis cycles is accomplished in the Mach number range of M=0.65–0.755. The reduced frequency and pitch amplitude also vary from k=0.03 to 0.1 and α₀=1–2.51°, respectively. Results show that the effect of increase in Mach number is to increase and decrease the lift coefficient during downstroke and upstroke, respectively, whereas at low reduced frequencies, the effect of increase in Mach number may lead to a reverse manner when airfoil moves toward its extremum angle of attack. Results also reveal that as the pitch amplitude varies, the shape of lift coefficient hysteresis cycle depends more on the pitch amplitude than on the appearance of shock. It is shown that as the Mach number increases, the incidence angles correspond to the extremum pitching moment, and depending on the reduced frequency, lie between zero and extremum angle of attack. These incidence angles shift toward the extremum angle of attack as the reduced frequency decreases. Results also show that the increase in pitch amplitude at low Mach number, in such a way that leads to the formation of shock around the extremum angle of attack, causes the extremum pitching moment to appear around these angles and at high Mach number, depending on the reduced frequency, the extremum pitching moment incidence angles would be between zero and extremum incidence angle.     

Airfoil flow instabilities induced by background flow oscillations

 The effect of background flow oscillations on transonic airfoil (NACA 0012) flow was investigated experimentally. The oscillations were generated by means of a rotating plate placed downstream of the airfoil. Owing to oscillating chocking of the flow caused by the plate, the airfoil flow periodically accelerated and decelerated. This led to strong variations in the surface pressure and the airfoil loading. The results are presented for two angles of attack, α=4° and α=8.5°, which correspond to the attached and separated steady airfoil flows, respectively. Received: 6 June 2000 / Accepted: 18 October 2001     

Effect of one-sided unsteady energy supply on aerodynamic characteristics of airfoils in a transonic flow

The possibility of controlling the aerodynamic characteristics of airfoils in transonic flight regimes by means of one-sided pulsed-periodic energy supply is studied. Based on the numerical solution of two-dimensional unsteady gas-dynamic equations, the change in the flow structure in the vicinity of a symmetric airfoil at different angles of attack and the aerodynamic characteristics of the airfoil as functions of the amount of energy supplied asymmetrically (with respect to the airfoil) are determined. The results obtained are compared with the data calculated for the flow past the airfoil at different angles of attack without energy supply. It is found that a given lift force can be obtained with the use of energy supply at a much better lift-to-drag ratio of the airfoil, as compared to the case of the flow past the airfoil at an angle of attack. The moment characteristics of the airfoil are found. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 49, No. 6, pp. 82–87, November–December, 2008.     

Assessment of 2D/3D numerical modeling for deep dynamic stall experiments

The results of computational fluid dynamics (CFD) simulations in two and three spatial dimensions are compared to pressure measurements and particle image velocimetry (PIV) flow surveys to assess the suitability of numerical models for the simulation of deep dynamic stall experiments carried out on a pitching NACA 23012 airfoil. A sinusoidal pitching motion with a 10° amplitude and a reduced frequency of 0.1 is imposed around two different mean angles of attack of 10° and 15°. The comparison of the airloads curves and of the pressure distribution over the airfoil surface shows that a three-dimensional numerical model can better reproduce the flow structures and the airfoil performance for the deep dynamic stall regime. Also, the vortical structures observed by PIV in the flow field are better captured by the three-dimensional model. This feature highlighted the relevance of three-dimensional effects on the flow field in deep dynamic stall.