Wing kinematics measurement and aerodynamics of free-flight maneuvers in drone-flies

The time courses of wing and body kinematics of two free-flying drone-flies, as they performed saccades, were measured using 3D high-speed video, and the morpho- logical parameters of the wings and body of the insects were also measured. The measured wing kinematics was used in a Navier-Stokes solver to compute the aerodynamic forces and moments acting on the insects. The main results are as following. （1） The turn is mainly a 90° change of heading. It is made in about 10 wingbeats （about 55 ms）. It is of interest to note that the number of wingbeats taken to make the turn is approximately the same as and the turning time is only a little different from that of fruitflies measured recently by the same approach, even if the weight of the droneflies is more than 100 times larger than that of the fruitflies. The long axis of body is about 40° from the horizontal during the maneuver. （2） Although the body rotation is mainly about a vertical axis, a relatively large moment around the yaw axis （axis perpendicular to the long axis of body）, called as yaw moment, is mainly needed for the turn, because moment of inertial of the body about the yaw axis is much larger than that about the long axis. （3） The yaw moment is mainly pro- duced by changes in wing angles of attack： in a right turn, for example, the dronefly lets its right wing to have a rather large angle of attack in the downstroke （generally larger than 50°） and a small one in the upstroke to start the turn, and lets its left wing to do so to stop the turn, unlike the fruitflies who generate the yaw moment mainly by changes in the stroke plane and stroke amplitude.     

Aerodynamics and mechanisms of elementary morphing models for flapping wing in forward flight of bat

Large active wing deformation is a significant way to generate high aerodynamic forces required in bat's flapping flight. Besides the twisting, elementary morphing models of a bat wing are proposed, including wing-bending in the spanwise direction, wing-cambering in the chordwise direction, and wing area-changing. A plate of aspect ratio 3 is used to model a bat wing, and a three-dimensional unsteady panel method is used to predict the aerodynamic forces. It is found that the cambering model has great positive influence on the lift, followed by the area-changing model and then the bending model. Further study indicates that the vortex control is a main mechanism to produce high aerodynamic forces. The mechanisms of aerodynamic force enhancement are asymmetry of the cambered wing and amplification effects of wing area-changing and wing bending. Lift and thrust are generated mainly during downstroke, and they are almost negligible during upstroke by the integrated morphing model-wing.     

The aerodynamic forces and pressure distribution of a revolving pigeon wing

The aerodynamic forces acting on a revolving dried pigeon wing and a flat card replica were measured with a propeller rig, effectively simulating a wing in continual downstroke. Two methods were adopted: direct measurement of the reaction vertical force and torque via a forceplate, and a map of the pressures along and across the wing measured with differential pressure sensors. Wings were tested at Reynolds numbers up to 108,000, typical for slow-flying pigeons, and considerably above previous similar measurements applied to insect and hummingbird wing and wing models. The pigeon wing out-performed the flat card replica, reaching lift coefficients of 1.64 compared with 1.44. Both real and model wings achieved much higher maximum lift coefficients, and at much higher geometric angles of attack (43°), than would be expected from wings tested in a windtunnel simulating translating flight. It therefore appears that some high-lift mechanisms, possibly analogous to those of slow-flying insects, may be available for birds flapping with wings at high angles of attack. The net magnitude and orientation of aerodynamic forces acting on a revolving pigeon wing can be determined from the differential pressure maps with a moderate degree of precision. With increasing angle of attack, variability in the pressure signals suddenly increases at an angle of attack between 33° and 38°, close to the angle of highest vertical force coefficient or lift coefficient; stall appears to be delayed compared with measurements from wings in windtunnels.     

A flow control mechanism in wing flapping with stroke asymmetry during insect forward flight

A theoretical modeling approach as well as an unsteady analytical method is used to study aerodynamic characteristics of wing flapping with asymmetric stroke-cycles in connection with an oblique stroke plane during insect forward flight. It is revealed that the aerodynamic asymmetry between the downstroke and the upstroke due to stroke-asymmetrical flapping is a key to understand the flow physics of generation and modulation of the lift and the thrust. Predicted results for examples of given kinematics validate more specifically some viewpoints that the wing lift is more easily produced when the forward speed is higher and the thrust is harder, and the lift and the thrust are generated mainly during downstroke and upstroke, respectively. The effects of three controlling parameters, i.e. the angles of tilted stroke plane, the different downstroke duration ratios, and the different angles of attack in both down- and up-stroke, are further discussed. It is found that larger oblique angles of stroke planes generate larger thrust but smaller lift; larger downstroke duration ratios lead to larger thrust, while making little change in lift and input aerodynamic power; and again, a small increase of the angle of attack in downstroke or upstroke may cause remarkable changes in aerodynamic performance in the relevant stroke.The project supported by the National Natural Science Foundation of China (10072066, 90305009) and the Chinese Academy of Sciences (KJCX-SW-L04, KJCX2-SW-L2)The English text was polished by Ron Marshall.     

Aerodynamic force and flow structures of two airfoils in flapping motions

Aerodynamic force and flow structures of two airfoils in a tandem configuration in flapping motions are studied, by solving the Navier-Stokes equations in moving overset grids. Three typical phase differences between the fore- and aftairfoil flapping cycles are considered. It is shown that: (1) in the case of no interaction (single airfoil), the time average of the vertical force coefficient over the downstroke is 2.74, which is about 3 times as large as the maximum steady-state lift coefficient of a dragonfly wing; the time average of the horizontal force coefficient is 1.97, which is also large. The reasons for the large force coefficients are the acceleration at the beginning of a stroke, the delayed stall and the “pitching-up” motion near the end of the stroke. (2) In the cases of two-airfoils, the time-variations of the force and moment coefficients on each airfoil are broadly similar to that of the single airfoil in that the vertical force is mainly produced in downstroke and the horizontal force in upstroke, but very large differences exist due to the interaction. (3) For in-phase stroking, the major differences caused by the interaction are that the vertical force on FA in downstroke is increased and the horizontal force on FA in upstroke decreased. As a result, the magnitude of the resultant force is almost unchanged but it inclines less forward. (4) For counter stroking, the major differences are that the vertical force on AA in downstroke and the horizontal force on FA in upstroke are decreased. As a result, the magnitude of the resultant force is decreased by about 20 percent but its direction is almost unchanged. (5) For 90°-phase-difference stroking, the major differences are that the vertical force on AA in downstroke and the horizontal force on FA in upstroke are decreased greatly and the horizontal force on AA in upstroke increased. As a result, the magnitude of the resultant force is decreased by about 28% and it inclines more forward. (6) Among the three cases of phase angles, inphase flapping produces the largest vertical force (also the largest resultant force); the 90°-phase-difference flapping results in the largest horizontal force, but the smallest resultant force.     

An unconventional mechanism of lift production during the downstroke in a hovering bird (<Emphasis Type="Italic">Zosterops japonicus</Emphasis>)

An unconventional mechanism of ventral clap is exploited by hovering passerines to produce lift. Quantitative visualization of the wake flow, analysis of kinematics and evaluation of the transient lift force was conducted to dissect the biomechanical role of the ventral clap in the asymmetrical hovering flight of passerines. The ventral clap can first abate and then augment lift production during the downstroke; the net effect of the ventral clap on lift production is, however, positive because the extent of lift augmentation is greater than the extent of lift abatement. Moreover, the ventral clap is inferred to compensate for the zero lift production of the upstroke because the clapping wings induce a substantial elevation of the lift force at the end of the downstroke. Overall, our observations shed light on the aerodynamic function of the ventral clap and offer biomechanical insight into how a bird hovers without kinematically mimicking hovering hummingbirds.     

膜扑翼飞行器的变形研究

最近昆虫翼的变形成了研究热点,而扑翼飞行器的变形力学研究却几乎无人问津.然而,无论昆虫、鸟类还是扑翼飞行器在飞行时,翼的变形都是存在的,要精确计算翼扑动产生的气动力,必须考虑其变形.本文比较了导致变形产生的膜扑翼飞行器的惯性力和气动力在一个周期中的变化情况,发现它们的峰值比值为2左右,然后提出了在随体坐标系中的固支边界条件,采用有限元法计算了惯性力和气动力分别对变形的影响,发现扑翼飞行器的气动力对变形的影响是不可忽略的重要因素,而惯性力与气动力的合力引起的最大正变形发生在下扑初始阶段,最大负变形发生在上扑初始阶段.本文为扑翼飞行器的设计提供了力学分析基础.     

Force production and asymmetric deformation of a flexible flapping wing in forward flight

Insect wings usually are flexible and deform significantly under the combined inertial and aerodynamic load. To study the effect of wing flexibility on both lift and thrust production in forward flight, a two-dimensional numerical simulation is employed to compute the fluid–structure interaction of an elastic wing section translating in an inclined stroke plane while pitching around its leading ledge. The effects of the wing stiffness, mass ratio, stroke plane angle, and flight speed are considered. The results show that the passive pitching due to wing deformation can significantly increase thrust while either maintaining lift at the same level or increasing it simultaneously. Another important finding is that even though the wing structure and actuation kinematics are symmetric, chordwise deformation of the wing shows a larger magnitude during upstroke than during downstroke. The asymmetry is more pronounced when the wing has a low mass ratio so that the fluid-induced deformation is significant. Such an aerodynamic cause may serve as an additional mechanism for the asymmetric deformation pattern observed in real insects.     

A simulation-based study on longitudinal gust response of flexible flapping wings

Winged animals such as insects are capable of flying and surviving in an unsteady and unpredictable aerial environment. They generate and control aerodynamic forces by flapping their flexible wings. While the dynamic shape changes of their flapping wings are known to enhance the efficiency of their flight, they can also affect the stability of a flapping wing flyer under unpredictable disturbances by responding to the sudden changes of aerodynamic forces on the wing. In order to test the hypothesis, the gust response of flexible flapping wings is investigated numerically with a specific focus on the passive maintenance of aerodynamic forces by the wing flexibility. The computational model is based on a dynamic flight simulator that can incorporate the realistic morphology, the kinematics, the structural dynamics, the aerodynamics and the fluid–structure interactions of a hovering hawkmoth. The longitudinal gusts are imposed against the tethered model of a hovering hawkmoth with flexible flapping wings. It is found that the aerodynamic forces on the flapping wings are affected by the gust, because of the increase or decrease in relative wingtip velocity or kinematic angle of attack. The passive shape change of flexible wings can, however, reduce the changes in the magnitude and direction of aerodynamic forces by the gusts from various directions, except for the downward gust. Such adaptive response of the flexible structure to stabilise the attitude can be classified into the mechanical feedback, which works passively with minimal delay, and is of great importance to the design of bio-inspired flapping wings for micro-air vehicles.     

Effects of unsteady deformation of flapping wing on its aerodynamic forces

Effects of unsteady deformation of a flapping model insect wing on its aerodynamic force production are studied by solving the Navier-Stokes equations on a dynamically deforming grid.Aerodynamic forces on the flapping wing are not much affected by considerable twist,but affected by camber deformation.The effect of combined camber and twist deformation is similar to that of camber deformation.With a deformation of 6% camber and 20°twist(typical values observed for wings of many insects),lift is increased bv 10%～20%and lift-to-drag ratio by around 10%compared with the case of a rigid flat-plate wing.As a result.the deformation can increase the maximum lift coefficient of an insect.and reduce its power requirement for flight.For example,for a hovering bumblebee with dynamically deforming wings(6?mber and 20°twist),aerodynamic power required is reduced by about 16%compared with the case of rigid wings.     

A computational study of the wing–wing and wing–body interactions of a model insect

The aerodynamic interaction between the contralateral wings and between the body and wings of a model insect are studied, by using the method of numerically solving the Navier-Stokes equations over moving overset grids, under typical hovering and forward flight conditions. Both the interaction between the contralateral wings and the interaction between the body and wings are very weak, e.g. at hovering, changes in aerodynamic forces of a wing due to the present of the other wing are less than 3% and changes in aerodynamic forces of the wings due to presence of the body are less than 2%. The reason for this is as following. During each down- or up-stroke, a wing produces a vortex ring, which induces a relatively large jet-like flow inside the ring but very small flow outside the ring. The vortex rings of the left and right wings are on the two sides of the body. Thus one wing is outside vortex ring of the other wing and the body is outside the vortex rings of the left and right wings, resulting in the weak interactions.     

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.     

Structural Analysis of a Dragonfly Wing

Dragonfly wings are highly corrugated, which increases the stiffness and strength of the wing significantly, and results in a lightweight structure with good aerodynamic performance. How insect wings carry aerodynamic and inertial loads, and how the resonant frequency of the flapping wings is tuned for carrying these loads, is however not fully understood. To study this we made a three-dimensional scan of a dragonfly (Sympetrum vulgatum) fore- and hindwing with a micro-CT scanner. The scans contain the complete venation pattern including thickness variations throughout both wings. We subsequently approximated the forewing architecture with an efficient three-dimensional beam and shell model. We then determined the wing’s natural vibration modes and the wing deformation resulting from analytical estimates of 8 load cases containing aerodynamic and inertial loads (using the finite element solver Abaqus). Based on our computations we find that the inertial loads are 1.5 to 3 times higher than aerodynamic pressure loads. We further find that wing deformation is smaller during the downstroke than during the upstroke, due to structural asymmetry. The natural vibration mode analysis revealed that the structural natural frequency of a dragonfly wing in vacuum is 154 Hz, which is approximately 4.8 times higher than the natural flapping frequency of dragonflies in hovering flight (32.3 Hz). This insight in the structural properties of dragonfly wings could inspire the design of more effective wings for insect-sized flapping micro air vehicles: The passive shape of aeroelastically tailored wings inspired by dragonflies can in principle be designed more precisely compared to sail like wings —which can make the dragonfly-like wings more aerodynamically effective.     

The effects of corrugation and wing planform on the aerodynamic force production of sweeping model insect wings

The effects of corrugation and wing planform (shape and aspect ratio) on the aerodynamic force production of model insect wings in sweeping (rotating after an initial start) motion at Reynolds number 200 and 3500 at angle of attack 40° are investigated, using the method of computational fluid dynamics. A representative wing corrugation is considered. Wing-shape and aspect ratio (AR) of ten representative insect wings are considered; they are the wings of fruit fly, cranefly, dronefly, hoverfly, ladybird, bumblebee, honeybee, lacewing (forewing), hawkmoth and dragonfly (forewing), respectively (AR of these wings varies greatly, from 2.84 to 5.45). The following facts are shown. (1) The corrugated and flat-plate wings produce approximately the same aerodynamic forces. This is because for a sweeping wing at large angle of attack, the length scale of the corrugation is much smaller than the size of the separated flow region or the size of the leading edge vortex (LEV). (2) The variation in wing shape can have considerable effects on the aerodynamic force; but it has only minor effects on the force coefficients when the velocity at r ₂ (the radius of the second moment of wing area) is used as the reference velocity; i.e. the force coefficients are almost unaffected by the variation in wing shape. (3) The effects of AR are remarkably small: when AR increases from 2.8 to 5.5, the force coefficients vary only slightly; flowfield results show that when AR is relatively large, the part of the LEV on the outer part of the wings sheds during the sweeping motion. As AR is increased, on one hand, the force coefficients will be increased due to the reduction of 3-dimensional flow effects; on the other hand, they will be decreased due to the shedding of part of the LEV; these two effects approximately cancel each other, resulting in only minor change of the force coefficients.     

昆虫飞行的空气动力学

昆虫是最早出现、数量最多和体积最小的飞行者. 它们能悬停、跃升、急停、快速加速和转弯, 飞行技巧十分高超. 由于尺寸小, 因而翅膀的相对速度很小, 从而进行上述飞行所需的升力系数很大. 但昆虫翅膀的雷诺数又很低. 它们是如何在低雷诺数下产生高升力的, 是流体力学和生物学工作者都十分关心的问题. 近年来这一领域有了许多研究进展. 该文对这些进展进行综述, 并对今后工作提一些建议. 因2005 年前的工作已在几篇综述文章有了详细介绍, 该文主要介绍2005 年以来的工作. 首先简述昆虫翅的拍动运动及昆虫绕流的基本方程和相似参数; 然后对2005 年之前的工作做一简要回顾. 之后介绍2005 年后的进展, 依次为: 运动学观测; 前缘涡; 翅膀柔性变形及皱褶的影响; 拍动翅的尾涡结构; 翼/身、左右翅气动干扰及地面效应; 微小昆虫; 蝴蝶与蜻蜓; 机动飞行. 最后为对今后工作的建议.      

The effects of corrugation and wing planform on the aerodynamic force production of sweeping model insect wings

The effects of corrugation and wing planform (shape and aspect ratio) on the aerodynamic force production of model insect wings in sweeping (rotating after an initial start) motion at Reynolds number 200 and 3500 at angle of attack 40° are investigated, using the method of computational fluid dynamics. A representative wing corrugation is considered. Wing-shape and aspect ratio (AR) of ten representative insect wings are considered; they are the wings of fruit fly, cranefly, dronefly, hoverfly, ladybird, bumblebee, honeybee, lacewing (forewing), hawkmoth and dragonfly (forewing), respectively (AR of these wings varies greatly, from 2.84 to 5.45). The following facts are shown. (1) The corrugated and flat-plate wings produce approximately the same aerodynamic forces. This is because for a sweeping wing at large angle of attack, the length scale of the corrugation is much smaller than the size of the separated flow region or the size of the leading edge vortex (LEV). (2) The variation in wing shape can have considerable effects on the aerodynamic force; but it has only minor effects on the force coefficients when the velocity at r ₂ (the radius of the second moment of wing area) is used as the reference velocity; i.e. the force coefficients are almost unaffected by the variation in wing shape. (3) The effects of AR are remarkably small: when AR increases from 2.8 to 5.5, the force coefficients vary only slightly; flowfield results show that when AR is relatively large, the part of the LEV on the outer part of the wings sheds during the sweeping motion. As AR is increased, on one hand, the force coefficients will be increased due to the reduction of 3-dimensional flow effects; on the other hand, they will be decreased due to the shedding of part of the LEV; these two effects approximately cancel each other, resulting in only minor change of the force coefficients. The project supported by the National Natural Science Foundation of China (10232010 and 10472008) and Ph. D. Student Foundation of Chinese Ministry of Education (20030006022) The English text was polished by Keren Wang.     

Aerodynamic and functional consequences of wing compliance

A growing body of evidence indicates that a majority of insects experience some degree of wing deformation during flight. With no musculature distal to the wing base, the instantaneous shape of an insect wing is dictated by the interaction of aerodynamic forces with the inertial and elastic forces that arise from periodic accelerations of the wing. Passive wing deformation is an unavoidable feature of flapping flight for many insects due to the inertial loads that accompany rapid stroke reversals—loads that well exceed the mean aerodynamic force. Although wing compliance has been implicated in a few lift-enhancing mechanisms (e.g., favorable camber), the direct aerodynamic consequences of wing deformation remain generally unresolved. In this paper, we present new experimental data on how wing compliance may affect the overall induced flow in the hawkmoth, Manduca sexta. Real moth wings were subjected to robotic actuation in their dominant plane of rotation at a natural wing beat frequency of 25 Hz. We used digital particle image velocimetry at exceptionally high temporal resolution (2,100 fps) to assess the influence of wing compliance on the mean advective flows, relying on a natural variation in wing stiffness to alter the amount of emergent deformation (freshly extracted wings are flexible and exhibit greater compliance than those that are desiccated). We find that flexible wings yield mean advective flows with substantially greater magnitudes and orientations more beneficial to lift than those of stiff wings. Our results confirm that wing compliance plays a critical role in the production of flight forces. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Dynamic flight stability of hovering model insects:theory versus simulation using equations of motion coupled with Navier-Stokes equations

In the present paper,the longitudinal dynamic flight stability properties of two model insects are predicted by an approximate theory and computed by numerical simulation.The theory is based on the averaged model(which assumes that the frequency of wingbeat is sufficiently higher than that of the body motion,so that the flapping wings' degrees of freedom relative to the body can be dropped and the wings can be replaced by wingbeat-cycle-average forces and moments);the simulation solves the complete equation...     

The influence of the wake of a flapping wing on the production of aerodynamic forces

The effect of the wake of previous strokes on the aerodynamic forces of a flapping model insect wing is studied using the method of computational fluid dynamics. The wake effect is isolated by comparing the forces and flows of the starting stroke (when the wake has not developed) with those of a later stroke (when the wake has developed). The following has been shown. (1) The wake effect may increase or decrease the lift and drag at the beginning of a half-stroke (downstroke or upstroke), depending on the wing kinematics at stroke reversal. The reason for this is that at the beginning of the half-stroke, the wing ``impinges' on the spanwise vorticity generated by the wing during stroke reversal and the distribution of the vorticity is sensitive to the wing kinematics at stroke reversal. (2) The wake effect decreases the lift and increases the drag in the rest part of the half-stroke. This is because the wing moves in a downwash field induced by previous half-stroke's starting vortex, tip vortices and attached leading edge vortex (these vortices form a downwash producing vortex ring). (3) The wake effect decreases the mean lift by 6%–18% (depending on wing kinematics at stroke reversal) and slightly increases the mean drag. Therefore, it is detrimental to the aerodynamic performance of the flapping wing. The project supported by the National Natural Science Foundation of China (10232010) and the National Aeronautic Science Fund of China(03A51049) The English text was polished by Xing Zhang     

Application of digital particle image velocimetry to insect aerodynamics: measurement of the leading-edge vortex and near wake of a Hawkmoth

Some insects use leading-edge vortices to generate high lift forces, as has been inferred from qualitative smoke visualisations of the flow around their wings. Here we present the first Digital Particle Image Velocimetry (DPIV) data and quantitative analysis of an insect’s leading-edge vortex and near wake at two flight speeds. This allows us to describe objectively 2D slices through the flow field of a tethered Tobacco Hawkmoth (Manduca sexta). The near-field vortex wake appears to braodly resemble elliptical vortex loops. The presence of a leading-edge vortex towards the end of the downstroke is found to coincide with peak upward force production measured by a six-component force–moment balance. The topology of Manduca’s leading-edge vortex differs from that previously described because late in the downstroke, the structure extends continuously from wingtip across the thorax to the other wingtip.