Analysis of non-symmetrical flapping airfoils

Simulations have been done to assess the lift, thrust and propulsive efficiency of different types of non-symmetrical airfoils under different flapping configurations. The variables involved are reduced frequency, Strouhal number, pitch amplitude and phase angle. In order to analyze the variables more efficiently, the design of experiments using the response surface methodology is applied. Results show that both the variables and shape of the airfoil have a profound effect on the lift, thrust, and efficiency. By using non- symmetrical airfoils, average lift coefficient as high as 2.23 can be obtained. The average thrust coefficient and efficiency also reach high values of 2.53 and 0.61, respectively. The lift production is highly dependent on the airfoil＇s shape while thrust production is influenced more heavily by the variables. Efficiency falls somewhere in between. Two-factor interac- tions are found to exist among the variables. This shows that it is not sufficient to analyze each variable individually. Vorticity diagrams are analyzed to explain the results obtained. Overall, the S1020 airfoil is able to provide relatively good efficiency and at the same time generate high thrust and lift force. These results aid in the design of a better ornithopter＇s wing.     

Numerical analysis of the s1020 airfoils in tandem under different flapping configurations

The objective of this project is to improve the performance of the efficiency, thrust and lift of flapping wings in tandem arrangement. This research investigates the effect of the arrangement of the airfoils in tandem on the performance of the airfoils by varying the phase difference and distance between the airfoils. Three flapping configurations from an earlier phase of a research which gives high efficiency, thrust and lift are used in the tandem simulation. It is found all the different flapping configurations show improvement in the efficiency, thrust or lift when the distance between the two airfoils and the phase angle between the heaving positions of the two airfoils are optimal. The average thrust coefficient of the tandem arrangement managed to attain more than twice that of the single one （4.84 vs. 2.05）. On the other hand, the average lift coefficient of the tandem arrangement also increased to 4.59, as compared to the original single airfoil value of 3.04. All these results obtained will aid in the design of a better ornithopter with tandem wing arrangement.     

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.     

褶皱幅度对仿生翼型滑翔气动特性的影响

蜻蜓翅膀具有独特的褶皱状形貌．研究者们致力于利用仿生学原理,设计在低雷诺数条件下具有更优气动性能的褶皱翼型．本文采用计算流体力学方法,求解二维不可压Navier-Stokes方程组,探讨了四种翼型（平板翼型、流线翼型、小幅度褶皱翼型和大幅度褶皱翼型）的气动表现．在低雷诺数条件下得到以下结果：(1) 较小幅度的褶皱结构有利于增加升力和减小阻力．(2) 雷诺数变化时褶皱翼型的升力系数呈非线性变化;在特定雷诺数区间,幅度相近的褶皱翼型会发生相对气动优势的转变．(3) 褶皱结构内的回流区通过减小粘性阻力,使得翼型总阻力下降．(4) 翼型前缘的极小区域会产生脉冲高升力,对升力表现产生较大影响．这些结果表明,调整褶皱幅度是实现褶皱翼型气动优化的有效方案．     

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

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

On the Design of Lifting Airfoils with High Critical Mach Number Using Full Potential Theory

 We wish to construct airfoils that have the highest free-stream Mach number for a given set of geometric constraints for which the flow is nowhere supersonic. Nonlifting airfoils that maximize the critical Mach number for a given cross-sectional area are known to possess long sonic segments at their critical speed. To construct lifting airfoils, we proceed under the conjecture that an airfoil with a high value of has the longest possible arc length of sonic velocity over its upper and lower surface. In Kropinski et al. (1995) the lifting problem was tackled in transonic small-disturbance theory. In this paper we numerically construct lifting airfoils with high using the full potential theory and we show that these airfoils have significantly higher than some standard airfoils. We also construct airfoils with higher values of the lift coefficient, by relaxing the speed constraint on the lower surface of the airfoil to have a value less than sonic. Received 13 May 1996 accepted 12 September 1996)     

A new direct design method for the medium thickness wind turbine airfoil

The newly developed integral function of airfoil profiles based on Trajkovski conformal transform theory could be used to optimize the profiles for the thin thickness airfoil. However, it is hard to adjust the coefficients of the integral function for the medium thickness airfoil. B-spline curve has an advantage of local adjustment, which makes it to effectively control the airfoil profiles at the trailing edge. Therefore, a new direct design method for the medium thickness wind turbine airfoil based on airfoil integral expression and B-spline curve is presented in this paper. An optimal mathematical model of an airfoil is built. Two new airfoils with similar thickness, based on the new designed method and the original integral method, are designed. According to the comparative analysis, the CQU-A25 airfoil designed based on the new method exhibits better results than that of the CQU-I25 airfoil which is designed based on the original method. It is demonstrated that the new method is feasible to design wind turbine airfoils. Meanwhile, the comparison of the aerodynamic performance for the CQU-A25 airfoil and for the DU91-W2-250 airfoil is studied. Results show that the maximum lift coefficient and the maximum lift/drag ratio of the CQU-A25 airfoil are higher than the ones of DU91-W2-250 airfoil in the same condition. This new airfoil design method would make it possible to design other airfoils with different thicknesses.     

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

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

Flows around two airfoils performing fling and subsequent translation and translation and subsequent clap

The aerodynamic forces and flow structures of two airfoils performing “fling and subsequent translation“ and “translation and subsequent clap“ are studied by numerically solving the Navier-Stokes equations in moving overset grids. These motions are relevant to the flight of very small insects. The Reynolds number, based on the airfoil chord length c and the translation velocity U, is 17. It is shown that: (1) For two airfoils performing fling and subsequent translation, a large lift is generated both in the fling phase and in the early part of the translation phase. During the fling phase,a pair of leading edge vortices of large strength is generated; the generation of the vortex pair in a short period results in a large time rate of change of fluid impulse, which explains the large lift in this period. During the early part of the translation, the two leading edge vortices move with the airfoils;the relative movement of the vortices also results in a large time rate of change of fluid impulse, which explains the large lift in this part of motion. (In the later part of the translation, the vorticity in the vortices is diffused and convected into the wake.) The time averaged lift coefficient is approximately 2.4 times as large as that of a single airfoil performing a similar motion. (2) For two airfoils performing translation and subsequent clap, a large lift is generated in the clap phase. During the clap, a pair of trailing edge vortices of large strength are generated; again, the generation of the vortex pair in a short period (which results in a large time rate of change of fluid impulse) is responsible for the large lift in this period. The time averaged lift coefficient is approximately 1.6 times as large as that of a single airfoil performing a similar motion. (3) When the initial distance between the airfoils (in the case of clap, the final distance between the airfoils) varies from 0.1 to 0.2c, the lift on an airfoil decreases only slightly but the torque decreases greatly. When the distance is about lc， the interference effects between the two airfoils become very small.     

Experimental study of pitching and plunging airfoils at low Reynolds numbers

Measurements of the unsteady flow structure and force time history of pitching and plunging SD7003 and flat plate airfoils at low Reynolds numbers are presented. The airfoils were pitched and plunged in the effective angle of attack range of 2.4°–13.6° (shallow-stall kinematics) and ?6° to 22° (deep-stall kinematics). The shallow-stall kinematics results for the SD7003 airfoil show attached flow and laminar-to-turbulent transition at low effective angle of attack during the down stroke motion, while the flat plate model exhibits leading edge separation. Strong Re-number effects were found for the SD7003 airfoil which produced approximately 25 % increase in the peak lift coefficient at Re = 10,000 compared to higher Re flows. The flat plate airfoil showed reduced Re effects due to leading edge separation at the sharper leading edge, and the measured peak lift coefficient was higher than that predicted by unsteady potential flow theory. The deep-stall kinematics resulted in leading edge separation that led to formation of a large leading edge vortex (LEV) and a small trailing edge vortex (TEV) for both airfoils. The measured peak lift coefficient was significantly higher (~50 %) than that for the shallow-stall kinematics. The effect of airfoil shape on lift force was greater than the Re effect. Turbulence statistics were measured as a function of phase using ensemble averages. The results show anisotropic turbulence for the LEV and isotropic turbulence for the TEV. Comparison of unsteady potential flow theory with the experimental data showed better agreement by using the quasi-steady approximation, or setting C(k) = 1 in Theodorsen theory, for leading edge–separated flows.     

前后双扑翼水平距离及相位差对其气动性能的影响

通过在动态网格上求解Navier-Stokes方程,对前后双扑翼的非定常粘性流场进行了数值模拟和气动干扰分析,考察了前后翼不同水平距离和不同相位差对其气动力和气动效率的影响。结果表明,扑动前翼和静止后翼间的气动干扰在各种不同水平距离下都有利于气动特性的改善,但气动干扰的作用随着前后翼水平距离的增大而减弱;前后双翼扑动的相位差是影响气动性能的重要参数,两翼间的气动干扰是否有利则与相位差和水平距离有直接联系。     

Reynolds number, thickness and camber effects on flapping airfoil propulsion

The effect of varying airfoil thickness and camber on plunging and combined pitching and plunging airfoil propulsion at Reynolds number Re=200, 2000, 20 000 and 2×10⁶ was studied by numerical simulations for fully laminar and fully turbulent flow regimes. The thickness study was performed on 2-D NACA symmetric airfoils with 6-50% thick sections undergoing pure plunging motion at reduced frequency k=2 and amplitudes h=0.25 and 0.5, and for combined pitching and plunging motion at k=2, h=0.5, phase ?=90°, pitch angle θ_o=15° and 30° and the pitch axis was located at 1/3 of chord from leading edge. At Re=200 for motions where positive thrust is generated, thin airfoils outperform thick airfoils. At higher Re significant gains could be achieved both in thrust generation and propulsive efficiency by using a thicker airfoil section for plunging and combined motion with low pitch amplitude. The camber study was performed on 2-D NACA airfoils with varying camber locations undergoing pure plunging motion at k=2, h=0.5 and Re=20 000. Little variation in thrust performance was found with camber. The underlying physics behind the alteration in propulsive performance between low and high Reynolds numbers has been explored by comparing viscous Navier-Stokes and inviscid panel method results. The role of leading edge vortices was found to be key to the observed performance variation.     

Lift enhancement of airfoil and tip flow control for wind turbine

Two techniques that improve the aerodynamic performance of wind turbine airfoils are described. The airfoil S809, designed specially for wind turbine blades, and the airfoil FX60-100, having a higher lift-drag ratio, are selected to verify the flow control techniques. The flow deflector, fixed at the leading edge, is employed to control the boundary layer separation on the airfoil at a high angle of attack. The multi-island genetic algorithm is used to optimize the parameters of the flow deflector. The results indicate that the flow deflector can suppress the flow separation, delay the stall, and enhance the lift. The characteristics of the blade tip vortex, the wake vortex, and the surface pressure distributions of the blades are analyzed. The vortex diffuser, set up at the blade tip, is employed to control the blade tip vortex. The results show that the vortex diffuser can increase the total pressure coefficient of the core of the vortex, decrease the strength of the blade tip vortex, lower the noise, and improve the efficiency of the blade.     

平板大攻角绕流升力和阻力系数的计算

二维平板或二维对称薄翼型大攻角绕流升力和阻力系数与攻角之间存在的函数关系一般用数据表格的形式给出。本文根据垂直平板绕流阻力实验数据和对称薄翼型全攻角绕流实验数据,分析得到了平板大攻角绕流总压力及其升力分量和阻力分量系数的近似计算公式。结果表明：平板总压力系数约等于攻角正弦值的2倍;总压力的阻力分量系数约等于攻角正弦值平方的2倍;升力分量系数约为攻角2倍的正弦值。计算结果与两组试验数据具有较好的一致性。     

Investigation of horizontal-axis wind turbine (HAWT) blade threedimensional rotational effect based on field experiments

Field experiments are performed on a two-bladed 33 kW horizontal-axis wind turbine (HAWT). The pressures are measured with 191 pressure sensors positioned around the surfaces of seven spanwise section airfoils on one of the two blades. Three-dimensional (3D) and two-dimensional (2D) numerical simulations are performed, respectively, on the rotor and the seven airfoils of the blade. The results are compared with the experimental results of the pressure distribution on the seven airfoils and the lift coefficients. The 3D rotational effect on the blade aerodynamic characteristics is then studied with a numerical approach. Finally, some conclusions are drawn as follows. From the tip to the root of the blade, the experimental differential pressure of the blade section airfoil increases at first and then decreases gradually. The calculated 3D result of the pressure distribution on the blade surface is closer to that of the experiment than the 2D result. The 3D rotational effect has a significant impact on the blade surface flow and the aerodynamic load, leading to an increase of the differential pressure on the airfoils and their lift coefficient than that with the 2D one because of the stall delay. The influence of the 3D rotational effect on the wind turbine blade especially takes place on the sections with flow separation.     

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.     

两种椭圆翼型高速气动特性实验研究

新概念旋转机翼飞机的主机翼既能高速旋转作为旋翼,又可锁定作为固定翼,所以只能使用特殊的前后对称翼型。针对主机翼翼型的这一特殊要求,对16％相对厚度,相对弯度分别为0％和3％的两种椭圆翼型的高速气动特性进行了风洞实验研究,试验分别在中国空气动力研究发展中心FL-21风洞和荷兰代尔夫特大学TST-27风洞进行,采用表面测压和尾排型阻测量技术。试验结果的对比分析表明,有弯度椭圆翼型的升力和力矩特性优于无弯度椭圆翼型,而阻力特性和最大升阻比劣于无弯度椭圆翼型。试验结果为旋转机翼飞机主机翼翼型的选取提供了参考。     

一种可用于机翼优化设计中的气动噪声预测模型

发展了一种可用于翼型/机翼外形设计中的气动噪声快速预测方法。相较于传统的半经验噪声预测方法,该方法以两方程非线性k-ε湍流模型模化雷诺应力的雷诺平均方程为背景,考虑了升力系数、三维流动效应以及机翼几何参数等因素对后缘噪声的影响。而相对于直接数值模拟或声类比拟方法,该方法虽不能准确预测噪声强度,但其计算量小,能给出不同翼型/机翼的相对总声压级,以及总声压级随升力系数的变化情况,易于应用于翼型/机翼气动外形优化设计中。通过计算分析二维NACA0012翼型几何参数或来流状态的改变所带来的气动噪声差异,与ANOPP软件及Brooks等计算结果进行对比,验证了该模型的可靠性。最后,计算分析二维、三维翼型/机翼气动噪声,凸显该方法在翼型/机翼气动外形优化设计中的应用价值。     

Design and development considerations for biologically inspired flapping-wing micro air vehicles

In this paper, the decade of numerical and experimental investigations leading to the development of the authors’ unique flapping-wing micro air vehicle is summarized. Early investigations included the study of boundary layer energization by means of a small flapping foil embedded in a flat-plate boundary layer, the reduction of the recirculatory flow region behind a backward-facing step by means of a small flapping foil, and the reduction or suppression of flow separation behind blunt or cusped airfoil trailing edges by flapping a small foil located in the wake flow region. These studies were followed by systematic investigations of the aerodynamic characteristics of single flapping airfoils and airfoil combinations. These unsteady flows were described using flow visualization, laser-Doppler velocimetry in addition to panel and Navier–Stokes computations. It is then shown how this flapping-wing database was used to conceive, design and develop a micro air vehicle which has a fixed wing for lift and two flapping wings for thrust generation. While animal flight is characterized by a coupled force generation, the present design appears to separate lift and thrust. However, in fact, the performance of one surface is closely coupled to the other surfaces.     

Hodograph design of lifting airfoils with high critical mach numbers

We wish to construct airfoils that have the highest free-stream Mach number M for a given set of geometric constraints for which the flow is nowhere supersonic. Nonlifting airfoils which maximize M for a given thickness ratio are known to possess long sonic segments at their critical speed. To construct lifting airfoils, we proceed under the conjecture that the optimal airfoil satisfying a given set of constraints is the one possessing the longest possible arc length of sonic velocity. A boundary-value problem is formulated in the hodograph plane using transonic small-disturbance theory whose solution determines an airfoil with long sonic arcs. For small lift coefficients, the hodograph domain covers two Riemann sheets and a finite-difference method is used to solve the boundary-value problem on this domain. A numerical integration of the solution around the boundary yields an airfoil shape, and three examples are discussed. The performance of these airfoils is compared with standard airfoils having the same lift coefficient and , and it is shown that the calculated airfoils have a 6%–10% increase in critical M .This research was supported by the Air Force Office of Scientific Research under Grant F49620-93-1-0022DEF and by the National Science Foundation under Grant DMS-9157546. Research support for M.C.A.K. was given by the IBM Corporation under a Graduate Research Fellowship.