对称翼型低雷诺数小攻角升力系数非线性现象研究

采用Rogers发展的三阶Roe格式,求解非定常不可压N-S方程,时间方向为二阶精度双时间步方法, 数值模拟了对称翼型SD8020低雷诺数(Re=40000,100000)条件下,流场层流分离涡结构和升力系数随攻角的变化．同试验比较证明了数值模拟的正确性．通过对数值模拟时均化流场结果的详细分析,发现对称翼型在小雷诺数0°攻角附近出现的层流分离泡,其内部结构和演化规律都不同于经典层流分离泡模型,从而提出了一种后缘层流分离泡模型．并应用该模型对对称翼型小攻角低雷诺数流场特性以及升力系数非线性效应的形成机理进行了研究和解释．     

低雷诺数下柔性翼型气动性能分析

基于流固耦合方法对吸力面5%至95%弦长处为三段柔性结构的NACA0012翼型绕流进行了数值模拟,研究了不同弹性模量下柔性翼型的气动性能和结构响应.结果表明:在大攻角下,翼面变形影响着翼型表面的非定常流场,起到延缓失速和提高升力的作用;失速后柔性翼的升力系数下降得较为缓慢,且柔性越大,升力系数下降得越平缓;适当减小弹性模量能够提高翼型的气动性能,然而弹性模量过小反而不利于翼型气动性能的提升,并且翼面会产生大幅度的振动.     

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

针对低雷诺数翼型气动性能差的特点, 通过介质阻挡放电(dielectric barrier discharge, DBD)等离子体激励控制的方法, 提高翼型低雷诺数下的气动特性,改善其流场结构. 采用二维准直接数值模拟方法求解非定常不可压Navier-Stokes方程,对具有俯仰运动的NACA0012翼型的低雷诺数流动展开数值模拟.同时将介质阻挡放电激励对流动的作用以彻体力源项的形式加入Navier-Stokes方程,通过数值模拟探究稳态DBD等离子体激励对俯仰振荡NACA0012翼型气动特性和流场特性的影响.为了进行流动控制, 分别在上下表面的前缘和后缘处安装DBD等离子体激励器,并提出四种激励器的开环控制策略,通过对比研究了这些控制策略在不同雷诺数、不同减缩频率以及激励位置下的控制效果.通过流场结构和动态压强分析了等离子体进行流场控制的机理. 结果表明,前缘DBD控制中控制策略B(负攻角时开启上表面激励器,正攻角时开启下表面激励器)效果最好,后缘DBD控制中控制策略C(逆时针旋转时开启上表面激励器,顺时针旋转时开启下表面激励器)效果最好,前缘DBD控制效果会随着减缩频率的增大而下降, 同时会导致阻力增大.而后缘DBD控制可以减小压差阻力, 优于前缘DBD控制,对于计算的所有减缩频率(5.01~11.82)都有较好的增升减阻效果.在不同雷诺数下, DBD控制的增升效果较为稳定, 而减阻效果随着雷诺数的降低而变差,这是由流体黏性效应增强导致的.      

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

针对低雷诺数翼型气动性能差的特点,通过介质阻挡放电(dielectric barrier discharge, DBD)等离子体激励控制的方法,提高翼型低雷诺数下的气动特性,改善其流场结构.采用二维准直接数值模拟方法求解非定常不可压Navier-Stokes方程,对具有俯仰运动的NACA0012翼型的低雷诺数流动展开数值模拟.同时将介质阻挡放电激励对流动的作用以彻体力源项的形式加入Navier-Stokes方程,通过数值模拟探究稳态DBD等离子体激励对俯仰振荡NACA0012翼型气动特性和流场特性的影响.为了进行流动控制,分别在上下表面的前缘和后缘处安装DBD等离子体激励器,并提出四种激励器的开环控制策略,通过对比研究了这些控制策略在不同雷诺数、不同减缩频率以及激励位置下的控制效果.通过流场结构和动态压强分析了等离子体进行流场控制的机理.结果表明,前缘DBD控制中控制策略B(负攻角时开启上表面激励器,正攻角时开启下表面激励器)效果最好,后缘DBD控制中控制策略C(逆时针旋转时开启上表面激励器,顺时针旋转时开启下表面激励器)效果最好,前缘DBD控制效果会随着减缩频率的增大而下降,同时会导致阻力增大.而后缘DBD控制可以减小压差阻力,优于前缘DBD控制,对于计算的所有减缩频率(5.01~11.82)都有较好的增升减阻效果.在不同雷诺数下, DBD控制的增升效果较为稳定,而减阻效果随着雷诺数的降低而变差,这是由流体黏性效应增强导致的.     

顺桨风力机塔叶干涉翼型非稳气动力时频特性

由于风力机叶片与塔筒流场相互干涉,实际气动力与理想情况存在较大差异,这种干涉作用造成的气动力差异给叶片与塔筒结构可靠性带来不可忽视的影响.以翼型DU91-W2-250为研究对象,采用瞬态数值分析与本征正交分解方法,考虑叶片和塔筒流场相互干涉作用,分析顺桨工况翼型非稳气动力时频特性及其影响规律,量化不同雷诺数下塔叶相对位置及几何参数对气动力均值、波动幅度和频率的影响程度,通过流场模态能量分布形态分析,揭示流场干涉对气动力的影响机制.结果表明,翼型气动中心至塔筒几何中心的垂直距离、水平距离以及塔筒直径相对于翼型弦长的无量纲参数y*,x*和D*对气动力均有不同程度影响,其中y*对升阻力系数均值影响最大,对频率无明显影响,y*绝对值越大,C_l均值越接近单翼型C_l值,y*绝对值越小升阻力系数波动幅度越大,y*从-12增大到12,升力系数均值最小值为-0.48,最大值为1.16;x*减小和D*增大,反向阻力均值增大,波动幅度增大,波动频率略有下降,当x*小于临界值5时,带塔翼型阻力均值反向;在计算范围内,带塔翼型升力系数均值相对于单翼型升力系数最大偏差为...     

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

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

雷诺数对运输类飞机气动特性影响的试验研究

开展了某背景飞机的高低速变雷诺数风洞试验,并对试验结果进行了分析;结合国外的一些变雷诺数风洞试验结论,给出了大展弦比运输类飞机雷诺数对升力、阻力以及俯仰力矩特性的影响规律。试验研究表明:升力特性方面,在中小迎角下,随着雷诺数增大,迎角相同时升力增加,大迎角时最大升力系数增大,失速迎角增大;阻力特性方面,随着雷诺数的增大阻力减小,并且最小阻力系数随着雷诺数对数的增加基本呈线性减小;俯仰力矩特性方面,随着雷诺数增大机翼后部载荷增大,低头力矩增大,稳定性增强。除此之外,文中还从雷诺数对气动特性影响机理角度出发,分析了这些规律形成的原因。     

基于涡流发生器的翼型失速流动控制及雷诺数效应影响研究

针对所设计的三角形涡流发生器开展用于翼型失速流动控制的风洞实验研究,重点讨论涡流发生器几何参数、方向角、安装位置及实验雷诺数等因素对翼型失速流动控制的影响。实验结果表明:涡流发生器作用下,在干净翼失速迎角后能够形成一个升力几乎不随迎角变化的相对稳定的高升力状态,抑制了失速流动的发生,与此同时阻力大幅下降;本文所设计的涡流发生器方向角过大时会削弱翼型失速流动控制的效果;同一涡流发生器作用下雷诺数过大其失速流动控制效果会急剧恶化,第一种涡流发生器控制翼型失速的雷诺数有效范围略宽于第二种涡流发生器。     

低雷诺数翼型蒙皮主动振动气动特性及流场结构数值研究

针对低雷诺数(Re)翼型气动性能差的特点,文章通过对翼型柔性蒙皮施加主动振动的方法,提高翼型低Re下的气动特性,改善其流场结构.采用带预处理技术的Roe方法求解非定常可压缩Navier-Stokes方程,对NACA4415翼型低Re流动展开数值模拟.通过时均化和非定常方法对比柔性蒙皮固定和振动两种状态下的升阻力气动特性和层流分离流动结构.初步研究工作表明在低Re下柔性蒙皮采用合适的振幅和频率,时均化升阻力特性显著提高,分离泡结构由后缘层流分离泡转变为近似的经典长层流分离泡,分离点后移,分离区缩小.在此基础上,文章更加细致研究了柔性蒙皮两种状态下单周期内的层流分离结构及壁面压力系数分布非定常特性和演化规律.蒙皮固定状态下分离区前部流场结构和压力分布基本保持稳定,表现为近似定常分离,仅在后缘位置出现类似于卡门涡街的非定常流动现象.柔性蒙皮振动时从分离点附近开始便产生分离涡,并不断向下游移动、脱落,表现为非定常分离并出现大范围的压力脉动.蒙皮振动使流体更加靠近壁面运动,大尺度的层流分离现象得到有效抑制.      

UNSTEADY SEPARATION OF FLOW AROUND AIRFOIL WITH LOCAL ELASTIC STRUCTURE

对在低雷诺数下局部弹性翼型绕流中,局部弹性导致的自激振动所产生的复杂非定常流动分离现象和描述方法进行了分析．采用ALE—CBS方法数值模拟了具有可动边界的绕流流场问题,同时采用Galerkin方法求解局部弹性结构的控制方程．着重研究了翼型的局部弹性对流动分离和翼型性能的影响,并分别从Eulerian和Lagrangian的角度分析了局部弹性结构导致的不同非定常分离现象,其中Lagrangian角度可以方便地揭示出局部弹性翼型大幅度提高升力的机理和流动中的能量迁移．结果表明翼型的局部弹性对非定常分离和分离泡的演化过程有着明显的影响,可以使得流体质点由主流获取动量实现再附,并且在一定的攻角下可以将固定分离转变为移动分离,从而明显地提高了翼型的升力．     

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.     

An experimental study of a bio-inspired corrugated airfoil for micro air vehicle applications

An experimental study was conducted to investigate the aerodynamic characteristics of a bio-inspired corrugated airfoil compared with a smooth-surfaced airfoil and a flat plate at the chord Reynolds number of Re _C  = 58,000–125,000 to explore the potential applications of such bio-inspired corrugated airfoils for micro air vehicle designs. In addition to measuring the aerodynamic lift and drag forces acting on the tested airfoils, a digital particle image velocimetry system was used to conduct detailed flowfield measurements to quantify the transient behavior of vortex and turbulent flow structures around the airfoils. The measurement result revealed clearly that the corrugated airfoil has better performance over the smooth-surfaced airfoil and the flat plate in providing higher lift and preventing large-scale flow separation and airfoil stall at low Reynolds numbers (Re _C  < 100,000). While aerodynamic performance of the smooth-surfaced airfoil and the flat plate would vary considerably with the changing of the chord Reynolds numbers, the aerodynamic performance of the corrugated airfoil was found to be almost insensitive to the Reynolds numbers. The detailed flow field measurements were correlated with the aerodynamic force measurement data to elucidate underlying physics to improve our understanding about how and why the corrugation feature found in dragonfly wings holds aerodynamic advantages for low Reynolds number flight applications.     

Optimal airfoil shapes for low Reynolds number flows

Flow over NACA 0012 airfoil is studied at α = 4^° and 12^° for Re?500. It is seen that the flow is very sensitive to Re. A continuous adjoint based method is formulated and implemented for the design of airfoils at low Reynolds numbers. The airfoil shape is parametrized with a non‐uniform rational B‐splines (NURBS). Optimization studies are carried out using different objective functions namely: (1) minimize drag, (2) maximize lift, (3) maximize lift to drag ratio, (4) minimize drag and maximize lift and (5) minimize drag at constant lift. The effect of Reynolds number and definition of the objective function on the optimization process is investigated. Very interesting shapes are discovered at low Re. It is found that, for the range of Re studied, none of the objective functions considered show a clear preference with respect to the maximum lift that can be achieved. The five objective functions result in fairly diverse geometries. With the addition of an inverse constraint on the volume of the airfoil the range of optimal shapes, produced by different objective functions, is smaller. The non‐monotonic behavior of the objective functions with respect to the design variables is demonstrated. The effect of the number of design parameters on the optimal shapes is studied. As expected, richer design space leads to geometries with better aerodynamic properties. This study demonstrates the need to consider several objective functions to achieve an optimal design when an algorithm that seeks local optima is used. Copyright © 2008 John Wiley & Sons, Ltd.     

Comparative studies of numerical methods for evaluating aerodynamic characteristics of two-dimensional airfoil at low Reynolds numbers

This study investigates the predictability of the aerodynamic performance of some numerical methods at low Reynolds numbers and their dependency on the geometric shape of airfoil. We conducted three-dimensional large-eddy simulations (3-D LES), two-dimensional laminar simulations (2-D Lam), and Reynolds-averaged Navier–Stokes simulations with Baldwin–Lomax (2-D RANS(BL)) and Spalart–Allmaras (2-D RANS(SA)) turbulence models. Although there is little discrepancy between the 3-D LES, 2-D Lam, and 2-D RANS(SA) results in terms of the lift and drag characteristics, significant differences are observed in the predictability of the separation and reattachment points. The predicted lift, separation, and reattachment points of the 2-D Lam are qualitatively similar to those of the 3-D LES, except for high angles of attack at which a massive separation occurs. The 2-D RANS(SA) shows good predictability of the lift and separation points, but it does not estimate reattachment points accurately. The 2-D RANS(BL) fails to predict the precise separation points, which results in a poor lift predictability. These characteristics appear regardless of the airfoil geometry shapes. The results suggest that a 2-D Lam without any turbulence models can be used to estimate qualitative airfoil aerodynamic characteristics at the low Reynolds numbers.     

蜻蜓滑翔时柔性褶皱前翅气动特性分析

蜻蜓是自然界优秀的飞行家,滑翔是其常见且有效的飞行模式.蜻蜓优异的飞行能力来源于其翅膀的巧妙结构,褶皱是蜻蜓翅膀上最为显著的结构之一,不仅提高了翅膀的刚度,还改变了其气动特性,而飞行过程中柔性翅膀会产生变形是蜻蜓翅膀的另一特性.为揭示蜻蜓在滑翔时,柔性褶皱前翅的变形,探究褶皱和柔性的共同作用对其气动特性的影响,基于逆向工程,依据前人的测量数据和研究成果,通过三维建模软件建立了蜻蜓三维褶皱前翅的计算流体力学(computational fluiddynamics,CFD)模型和计算结构力学(computational structuralmechanics,CSD)模型,并通过模态分析验证了此模型有足够的精度.基于CFD方法和CFD/CSD双向流固耦合计算方法分别对蜻蜓滑翔飞行时刚性和柔性褶皱前翅的气动特性进行了数值模拟,结果表明,柔性褶皱前翅受气动载荷后,翅脉和翅膜产生形变,柔性前翅上下表面压力差相较于刚性前翅减小了,从而其升力和阻力也减小了,而在大攻角时,变形后的前缘脉诱导出比刚性前翅更强的前缘涡.因此在攻角小于10$^\circ$时刚性前翅的气动特性优于柔性前翅,继续增大攻角,柔性前翅的气动特性则优于刚性前翅.前翅受载后气动响应时间短,翅尖的变形最大,仅仅产生了垂直于翅膀所在平面方向上的变形,而没有发生扭转,翼根处受到应力最大,褶皱上凸部分承受蜻蜓滑翔时前翅的主要载荷.      

钝后缘风力机翼型的环量控制研究

钝后缘风力机翼型具有结构强度高、对表面污染不敏感等优点,但其较大的阻力系数使得翼型的整体气动特性不够理想. 利用环量控制方法对钝后缘风力机翼型进行了流动控制,以改善钝后缘风力机翼型的气动特性,减弱尾迹区脱体涡强度. 通过对钝后缘风力机翼型环量控制方法进行相关的数值模拟,对比研究了环量控制方法的增升减阻效果, 研究了环量控制下翼型升阻力特性随射流动量系数的变化规律,并对不同射流动量系数下环量控制方法的气动品质因子和控制效率进行了分析. 研究结果表明:环量控制方法能够大幅提升钝后缘风力机翼型的升力系数,同时有效地降低翼型的阻力系数; 翼型的升力系数随射流动量系数的增大而增大,表现出很明显的分离控制阶段和超环量控制阶段的变化规律; 射流能耗的功率系数随射流动量系数的增大而增大,且增长速率逐渐增大;实施环量控制方法后叶片的输出功率同样随射流动量系数增大而增大,但增长速率逐渐降低. 总体来说,环量控制方法可以有效地改善钝后缘风力机翼型的气动特性以及功率输出特性,在大型风力机流动控制中具有很好的应用前景.      

基于低雷诺数串置翼型气动特性研究

基于低雷诺数,应用二维CFD方法,对串置翼型的气动特性进行了数值模拟。比较了串置翼型布局与常规单独翼型的升阻特性。分析了翼差角度对串置翼升阻特性的影响,并且在保持鸭翼±5°偏角不变的同时,研究了两翼之间的距离以及安装的相对高度对整个串置翼型升阻特性的影响。发现上鸭翼与主翼相结合具有相对较高的最大升力系数和临界迎角, 可以显著改善串置翼型的气动特性。     

Low Reynolds unsteady flow simulation around NACA0012 airfoil with active flow control

This research numerically elucidates the effects of suction and blowing on the enhancement of unsteady aerodynamic characteristics of flows and their corresponding impact on stall delay over the well-known NACA0012 airfoil at various angles of attack (\( 12 \le {\text{AOA}} \le 20 \)) under low Reynolds numbers. For this purpose, an in-house solver written in C++ is developed. The numerical code utilizes the Jameson’s cell-centered finite volume numerical method accompanied by a progressive power-law preconditioning approach to suppress the stiffness of the governing equations. Many numerical simulations are performed over the suction-blowing control parameters, namely, the slot location (\( L_{j} \)), suction/blowing amplitudes (\( A_{j} \)), and suction/blowing angle (\( \theta_{j} \)). Most of the analyses are based on the measurements of the unsteady aerodynamic characteristics behaviors (such as lift, drag, moment coefficients, and stall phenomena) over the airfoil. The numerical results confirm that the unsteady behavior of the flow (vortex shedding) is weakened or approximately removed when suction is used, especially near the leading edge. In all of the test cases, the ratio of the average lift coefficient to the average drag coefficient increases with increasing suction and blowing amplitudes, except in the case of perpendicular blowing. Furthermore, the blowing is more sensitive to the blowing angle compared to the suction. From the suction and blowing results, it is concluded that the former has a more positive impact on the lift and drag characteristics, especially in the case of incompressible flow at Low-Reynolds regimes.     

Influence of heat transfer on the aerodynamic performance of a plunging and pitching NACA0012 airfoil at low Reynolds numbers

The unsteady low Reynolds number aerodynamics phenomena around flapping wings are addressed in several investigations. Elsewhere, airfoils at higher Mach numbers and Reynolds numbers have been treated quite comprehensively in the literature. It is duly noted that the influence of heat transfer phenomena on the aerodynamic performance of flapping wings configurations is not well studied. The objective of the present study is to investigate the effect of heat transfer upon the aerodynamic performance of a pitching and plunging NACA0012 airfoil in the low Reynolds number flow regime with particular emphasis upon the airfoil's lift and drag coefficients. The compressible Navier–Stokes equations are solved using a finite volume method. To consider the variation of fluid properties with temperature, the values of dynamic viscosity and thermal diffusivity are evaluated with Sutherland's formula and the Eucken model, respectively. Instantaneous and mean lift and drag coefficients are calculated for several temperature differences between the airfoil surface and freestream within the range 0–100 K. Simulations are performed for a prescribed airfoil motion schedule and flow parameters. It is learnt that the aerodynamic performance in terms of the lift C_L and drag C_D behavior is strongly dependent upon the heat transfer rate from the airfoil to the flow field. In the plunging case, the mean value of C_D tends to increase, whereas the amplitude of C_L tends to decrease with increasing temperature difference. In the pitching case, on the other hand, the mean value and the amplitude of both C_D and C_L decrease. A spectral analysis of C_D and C_L in the pitching case shows that the amplitudes of both C_D and C_L decrease with increasing surface temperature, whereas the harmonic frequencies are not affected.