Thrust augmentation of flapping airfoils in low Reynolds number flow using a flexible membrane

The unsteady aerodynamic thrust and aeroelastic response of a two-dimensional membrane airfoil under prescribed harmonic motion are investigated computationally with a high-order Navier–Stokes solver coupled to a nonlinear membrane structural model. The effects of membrane prestress and elasticity are examined parametrically for selected plunge and pitch–plunge motions at a chord-based Reynolds number of 2500. The importance of inertial membrane loads resulting from the prescribed flapping is also assessed for pure plunging motions. This study compares the period-averaged aerodynamic loads of flexible versus rigid membrane airfoils and highlights the vortex structures and salient fluid–membrane interactions that enable more efficient flapping thrust production in low Reynolds number flows.     

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.     

Flow field characteristics study of a flapping airfoil using computational fluid dynamics

The flow field of a flapping airfoil in Low Reynolds Number (LRN) flow regime is associated with complex nonlinear vortex shedding and viscous phenomena. The respective fluid dynamics of such a flow is investigated here through Computational Fluid Dynamics (CFD) based on the Finite Volume Method (FVM). The governing equations are the unsteady, incompressible two-dimensional Navier-Stokes (N-S) equations. The airfoil is a thin ellipsoidal geometry performing a modified figure-of-eight-like flapping pattern. The flow field and vortical patterns around the airfoil are examined in detail, and the effects of several unsteady flow and system parameters on the flow characteristics are explored. The investigated parameters are the amplitude of pitching oscillations, phase angle between pitching and plunging motions, mean angle of attack, Reynolds number (Re), Strouhal number (St) based on the translational amplitudes of oscillations, and the pitching axis location (x/c). It is shown that these parameters change the instantaneous force coefficients quantitatively and qualitatively. It is also observed that the strength, interaction, and convection of the vortical structures surrounding the airfoil are significantly affected by the variations of these parameters.     

Effect of flexure on aerodynamic propulsive efficiency of flapping flexible airfoil

The aim of present study is to investigate the effect of chord-wise flexure amplitude on unsteady aerodynamic characteristics for a flapping airfoil with various combinations of Reynolds number and reduced frequency. Unsteady, viscous flows over a single flexible airfoil in plunge motion are computed using conformal hybrid meshes. The dynamic mesh technique is applied to illustrate the deformation modes of the flexible flapping airfoil. In order to investigate the influence of the flexure amplitude on the aerodynamic performance of the flapping airfoil, the present study considers eight different flexure amplitudes (a₀) ranging from 0 to 0.7 in intervals of 0.1 under conditions of Re=10⁴, reduced frequency k=2, and dimensionless plunge amplitude h₀=0.4. The computed unsteady flow fields clearly reveal the formation and evolution of a pair of leading edge vortices along the body of the flexible airfoil as it undergoes plunge motion. Thrust-indicative wake structures are generated when the flexure amplitude of the airfoil is less than 0.5 of the chord length. An enhancement in the propulsive efficiency is observed for a flapping airfoil with flexure amplitude of 0.3 of the chord length. This study also calculates the propulsive efficiency and thrust under various Reynolds numbers and reduced frequency conditions. The results indicate that the propulsive efficiency has a strong correlation with the reduced frequency. It is found that the flow conditions which yield the highest propulsive efficiency correspond to Strouhal number St of 0.255.     

Experimental study on interaction of aerodynamics with flexible wings of flapping vehicles in hovering and cruise flight

Flapping wings are promising lift and thrust generators, especially for very low Reynolds numbers. To investigate aeroelastic effects of flexible wings (specifically, wing’s twisting stiffness) on hovering and cruising aerodynamic performance, a flapping-wing system and an experimental setup were designed and built. This system measures the unsteady aerodynamic and inertial forces, power usage, and angular speed of the flapping wing motion for different flapping frequencies and for various wings with different chordwise flexibility. Aerodynamic performance of the vehicle for both no wind (hovering) and cruise condition was investigated. Results show how elastic deformations caused by interaction of inertial and aerodynamic forces with the flexible structure may affect specific power consumption. This information was used here to find a more suitable structural design. The best selected design in our tests performs up to 30% better than others (i.e., less energy consumption for the same lift or thrust generation). This measured aerodynamic information could also be used as a benchmarking data for unsteady flow solvers.     

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.     

Characterization of vortical structures and loads based on time-resolved PIV for asymmetric hovering flapping flight

Flight agility, resistance to gusts, capability to hover coupled with a low noise generation might have been some of the reasons why insects are among the oldest species observed in nature. Biologists and aerodynamicists focused on analyzing such flight performances for diverse purposes: understanding the essence of flapping wings aerodynamics and applying this wing concept to the development of micro-air vehicles (MAVs). In order to put into evidence the fundamentally non-linear unsteady mechanisms responsible for the amount of lift generated by a flapping wing (Dickinson et al. in Science 284:1954–1960, 1999), experimental and numerical studies were carried out on typical insect model wings and kinematics. On the other hand, in the recent context of MAVs development, it is of particular interest to study simplified non-biological flapping configurations which could lead to lift and/or efficiency enhancement. In this paper, we propose a parametrical study of a NACA0012 profile undergoing asymmetric hovering flapping motions at Reynolds 1000. On the contrary to normal hovering, which has been widely studied as being the most common configuration observed in the world of insects, asymmetric hovering is characterized by an inclined stroke plane. Besides the fact that the vertical force is hence a combination of both lift and drag (Wang in J Exp Biol 207:1137–1150, 2004), the specificity of such motions resides in the vortex dynamics which present distinct behaviours, whether the upstroke angle of attack leads to a partially attached or a strong separated flow, giving more or less importance to the wake capture phenomenon. A direct consequence of the previous remarks relies on the enhancement of aerodynamic efficiency with asymmetry. If several studies reported results based on the asymmetric flapping motion of dragonfly, only few works concentrated on parametrizing asymmetric motions (e.g. Wang in Phys Rev Lett 85:2216–2219, 2000). The present study relies on TR-PIV measurements which allow determination of the vorticity fields and provide a basis to evaluate the resulting unsteady forces through the momemtum equation approach.     

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.     

昆虫拍翼方式的非定常流动物理再探讨

基于提出的理论模化方法来探讨昆虫拍翼方式的非定常流动物理. 以悬停飞行为 例,通过对拍翼运动的分析,不仅解释了昆虫利用高频拍翼的方式为何能够克服低雷诺数带 来的气动局限性(St \gg 1/Re),而且还指出高升力产生和调节的3个流动 控制因素：(1) 由于拍翼的变速运动即时引起了流体动力响应,这种附加惯性效应 可产生瞬时的高升力; (2) 保持前缘涡不脱离翼面有助于减少升力的下降; (3) 增大后缘涡的强度并加速其脱离后缘能够有效地提高升力.     

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.     

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

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

Experimental investigation of the effect of chordwise flexibility on the aerodynamics of flapping wings in hovering flight

Ornithopters or mechanical birds produce aerodynamic lift and thrust through the flapping motion of their wings. Here, we use an experimental apparatus to investigate the effects of a wing's twisting stiffness on the generated thrust force and the power required at different flapping frequencies. A flapping wing system and an experimental set-up were designed to measure the unsteady aerodynamic and inertial forces, power usage and angular speed of the flapping wing motion. A data acquisition system was set-up to record important data with the appropriate sampling frequency. The aerodynamic performance of the vehicle under hovering (i.e., no wind) conditions was investigated. The lift and thrust that were produced were measured for different flapping frequencies and for various wings with different chordwise flexibilities. The results show the manner in which the elastic deformation and inertial flapping forces affect the dynamical behavior of the wing. It is shown that the generalization of the actuator disk theory is, at most, only valid for rigid wings, and for flexible wings, the power P varies by a power of about 1.0  of the thrust T. This aerodynamic information can also be used as benchmark data for unsteady flow solvers.     

昆虫飞行的高升力机理

对近年来关于昆虫产生非定常高升力的研究进行了综述和归纳.这方面的工作对生物学研究和微型飞行器等微型机械的仿生设计有重要意义.研究表明:果蝇等昆虫翅膀的拍动运动可产生很大的非定常升力,其平均值是定常值的2~3倍,足够平衡昆虫的重量,并有较大的富余用于机动飞行;产生高升力有三个因素:一是拍动开始阶段翅的快速加速运动,二是拍动中的不失速机制,三是拍动结束阶段翅的快速上仰运动.人们从能耗的角度考察了这些非定常高升力机制的正确性和可行性.当作悬停飞行的果蝇用以上机制产生平衡其重量的升力时,其比功率(支持单位身体质量所需的功率)约为29W/kg, 生化/机械效率约为17%. 这些值与人们基于对昆虫肌肉力学特性的研究所预估的值接近.果蝇前飞时,其比功率随速度变化的曲线是一J形曲线,而不是象飞机或鸟的那样是一U形曲线;这与人们基于昆虫新陈代谢率的测量数据所推断的结果一致.对于蜻蜒等(功能上)有前、后两对翅膀的昆虫,有以下初步结果:翅的下拍主要产生升力,上挥主要产生推力;下拍时的平均升力系数可达2~3,十分大,上挥时的平均推力系数可达1~2, 也很大,它们主要由非定常效应产生;前、后翅的相互干扰并未起增大升力和推力的作用,反而有一定的不利作用.      

Effect of vortex shedding in unsteady aerodynamic forces for a low Reynolds number stationary wing at low angle of attack

This study elucidates the relation between wake vortex shedding and aerodynamic force fluctuations for a low Reynolds number wing from time resolved particle image velocimetry (TR-PIV) experimental measurements. The results reveal a periodic lift and drag variation within the shedding cycle and resolve the frequencies of those fluctuations from a proper orthogonal decomposition (POD) and power spectral density (PSD) analysis. To show the effect of vortex shedding on the body force fluctuations, the evolution of instantaneous aerodynamic forces is compared to the pressure field of the fluid flow and to the vortical structures in the wake of the airfoil. A six step model describing the vortex-force relation is proposed. It shows that changes in lift such as maximum lift and minimum lift are associated with the detachment of a vortex. It also shows that the minimum or local minimum drag value is obtained at the onset formation of a vortex on the airfoil wake. Similarly, the maximum or local maximum drag is obtained at the onset formation of the saddle on the airfoil wake. The model further explains the asymmetry observed in the unsteady drag force evolution. The model can be used to optimize flow control and fluid-structure interaction applications.     

An investigation into the effects of unsteady parameters on the aerodynamics of a low Reynolds number pitching airfoil

The growing applications of low Reynolds number (LRN) operating vehicles impose the need for accurate LRN flow solutions. These applications usually involve complex unsteady phenomena, which depend on the kinematics of the vehicle such as pitching, plunging, and flapping of a wing. The objective of the present study is to address the issues related to LRN aerodynamics of a harmonically pitching NACA0012 airfoil. To this end, the influence of unsteady parameters, namely, amplitude of oscillation, d, reduced frequency, k, and Reynolds number, Re, on the aerodynamic performance of the model is investigated. Computational fluid dynamics (CFD) is utilized to solve Navier–Stokes (N–S) equations discretized based on the Finite Volume Method (FVM). The resulting instantaneous lift coefficients are compared with analytical data from Theodorsen’s method. The simulation results reveal that d, k, and Re are of great importance in the aerodynamic performance of the system, as they affect the maximum lift coefficients, hysteresis loops, strength, and number of the generated vortices within the harmonic motion, and the extent of the so-called figure-of-eight phenomenon region. Thus, achieving the optimum lift coefficients demands a careful selection of these parameters.     

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.     

Effects of a trapped vortex cell on a thick wing airfoil

The effects of a trapped vortex cell (TVC) on the aerodynamic performance of a NACA0024 wing model were investigated experimentally at Re = 10⁶ and 6.67×10⁵6.67\times 10^{5}. The static pressure distributions around the model and the wake velocity profiles were measured to obtain lift and drag coefficients, for both the clean airfoil and the controlled configurations. Suction was applied in the cavity region to stabilize the trapped vortex. For comparison, a classical boundary layer suction configuration was also tested. The drag coefficient curve of the TVC-controlled airfoil showed sharp discontinuities and bifurcative behavior, generating two drag modes. A strong influence of the angle of attack, the suction rate and the Reynolds number on the drag coefficient was observed. With respect to the clean airfoil, the control led to a drag reduction only if the suction was high enough. Compared to the classical boundary layer suction configuration, the drag reduction was higher for the same amount of suction only in a specific range of incidence, i.e., α = −2° to α = 6° and only for the higher Reynolds number. For all the other conditions, the classical boundary layer suction configuration gave better drag performances. Moderate increments of lift were observed for the TVC-controlled airfoil at low incidence, while a 20% lift enhancement was observed in the stall region with respect to the baseline. However, the same lift increments were also observed for the classical boundary layer suction configuration. Pressure fluctuation measurements in the cavity region suggested a very complex interaction of several flow features. The two drag modes were characterized by typical unsteady phenomena observed in rectangular cavity flows, namely the shear layer mode and the wake mode.     

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.     

On the aerodynamic characteristics of hovering rigid and flexible hawkmoth-like wings

Insect wings are subjected to fluid, inertia and gravitational forces during flapping flight. Owing to their limited rigidity, they bent under the influence of these forces. Numerical study by Hamamoto et al. (Adv Robot 21(1–2):1–21, 2007) showed that a flexible wing is able to generate almost as much lift as a rigid wing during flapping. In this paper, we take a closer look at the relationship between wing flexibility (or stiffness) and aerodynamic force generation in flapping hovering flight. The experimental study was conducted in two stages. The first stage consisted of detailed force measurement and flow visualization of a rigid hawkmoth-like wing undergoing hovering hawkmoth flapping motion and simple harmonic flapping motion, with the aim of establishing a benchmark database for the second stage, which involved hawkmoth-like wing of different flexibility performing the same flapping motions. Hawkmoth motion was conducted at Re = 7,254 and reduced frequency of 0.26, while simple harmonic flapping motion at Re = 7,800 and 11,700, and reduced frequency of 0.25. Results show that aerodynamic force generation on the rigid wing is governed primarily by the combined effect of wing acceleration and leading edge vortex generated on the upper surface of the wing, while the remnants of the wake vortices generated from the previous stroke play only a minor role. Our results from the flexible wing study, while generally supportive of the finding by Hamamoto et al. (Adv Robot 21(1–2):1–21, 2007), also reveal the existence of a critical stiffness constant, below which lift coefficient deteriorates significantly. This finding suggests that although using flexible wing in micro air vehicle application may be beneficial in term of lightweight, too much flexibility can lead to deterioration in flapping performance in terms of aerodynamic force generation. The results further show that wings with stiffness constant above the critical value can deliver mean lift coefficient almost the same as a rigid wing when executing hawkmoth motion, but lower than the rigid wing when performing a simple harmonic motion. In all cases studied (7,800 ≤ Re ≤ 11,700), the Reynolds number does not alter the force generation significantly.