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.     

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.     

Reduced coupled flapping wing-fluid computational model with unsteady vortex wake

Insect flight research is propelled by their unmatched flight capabilities. However, complex underlying aerodynamic phenomena make computational modeling of insect-type flapping flight a challenging task, limiting our ability in understanding insect flight and producing aerial vehicles exploiting same aerodynamic phenomena. To this end, novel mid-fidelity approach to modeling insect-type flapping vehicles is proposed. The approach is computationally efficient enough to be used within optimal design and optimal control loops, while not requiring experimental data for fitting model parameters, as opposed to widely used quasi-steady aerodynamic models. The proposed algorithm is based on Helmholtz–Hodge decomposition of fluid velocity into curl-free and divergence-free parts. Curl-free flow is used to accurately model added inertia effects (in almost exact manner), while expressing system dynamics by using wing variables only, after employing symplectic reduction of the coupled wing-fluid system at zero level of vorticity (thus reducing out fluid variables in the process). To this end, all terms in the coupled body-fluid system equations of motion are taken into account, including often neglected terms related to the changing nature of the added inertia matrix (opposed to the constant nature of rigid body mass and inertia matrix). On the other hand—in order to model flapping wing system vorticity effects—divergence-free part of the flow is modeled by a wake of point vortices shed from both leading (characteristic for insect flight) and trailing wing edges. The approach is evaluated for a numerical case involving fruit fly hovering, while quasi-steady aerodynamic model is used as benchmark tool with experimentally validated parameters for the selected test case. The results indicate that the proposed approach is capable of mid-fidelity accurate calculation of aerodynamic loads on the insect-type flapping wings.

     

An analytic approach to theoretical modeling of highly unsteady viscous flow excited by wing flapping in small insects

Numerous studies on the aerodynamics of insect wing flapping were carried out on different approaches of flight investigations, model experiments, and numerical simulations, but the theoretical modeling remains to be explored. In the present paper, an analytic approach is presented to model the flow interactions of wing flapping in air for small insects with the surrounding flow fields being highly unsteady and highly viscous. The model of wing flapping is a 2-D flat plate, which makes plunging and pitching oscillations as well as quick rotations reversing its positions of leading and trailing edges, respectively, during stroke reversals. It contains three simplified aerodynamic assumptions: (i) unsteady potential flow; (ii) discrete vortices shed from both leading and trailing edges of the wing; (iii) Kutta conditions applied at both edges. Then the problem is reduced to the solution of the unsteady Laplace equation, by using distributed singularities, i.e., sources/sinks, and vortices in the field. To validate the present physical model and analytic method proposed via benchmark examples, two elemental motions in wing flapping and a case of whole flapping cycles are analyzed, and the predicted results agree well with available experimental and numerical data. This verifies that the present analytical approach may give qualitatively correct and quantitatively reasonable results. Furthermore, the total fluid-dynamic force in the present method can be decomposed into three parts: one due to the added inertial (or mass) effect, the other and the third due to the induction of vortices shed from the leading-and the trailing-edge and their images respectively, and this helps to reveal the flow control mechanisms in insect wing flapping. The project supported by the National Natural Science Foundation of China (10072066) and the Chinese Academy of Sciences (KJCX-SW-LO4, KJCX2-SW-L2)     

Numerical simulation of flapping‐wing insect hovering flight at unsteady flow

A computational fluid dynamics (CFD) analysis was conducted to study the unsteady aerodynamics of a virtual flying bumblebee during hovering flight. The integrated geometry of bumblebee was established to define the shape of a three‐dimensional virtual bumblebee model with beating its wings, accurately mimicking the three‐dimensional movements of wings during hovering flight. The kinematics data of wings documented from the measurement to the bumblebee in normal hovering flight aided by the high‐speed video. The Navier–Stokes equations are solved numerically. The solution provides the flow and pressure fields, from which the aerodynamic forces and vorticity wake structure are obtained. Insights into the unsteady aerodynamic force generation process are gained from the force and flow‐structure information. The CFD analysis has established an overall understanding of the viscous and unsteady flow around the virtual flying bumblebee and of the time course of instantaneous force production, which reveals that hovering flight is dominated by the unsteady aerodynamics of both the instantaneous dynamics and also the past history of the wing. A coherent leading‐edge vortex with axial flow and the attached wingtip vortex and trailing edge vortex were detected. The leading edge vortex, wing tip vortex and trailing edge vortex, which caused by the pressure difference between the upper and the lower surface of wings. The axial flow, which include the spanwise flow and chordwise flow, is derived from the spanwise pressure gradient and chordwise pressure gradient, will stabilize the vortex and gives it a characteristic spiral conical shape. Copyright © 2006 John Wiley & Sons, Ltd.     

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.     

Forward flight of a model butterfly: Simulation by equations of motion coupled with the Navier-Stokes equations

The forward flight of a model butterfly was studied by simulation using the equations of motion coupled with the Navier-Stokes equations. The model butterfly moved under the action of aerodynamic and gravitational forces, where the aerodynamic forces were generated by flapping wings which moved with the body, allowing the body oscillations of the model butterfly to be simulated. The main results are as follows: (1) The aerodynamic force produced by the wings is approximately perpendicular to the long-axis of body and is much larger in the downstroke than in the upstroke. In the downstroke the body pitch angle is small and the large aerodynamic force points up and slightly backward, giving the weight-supporting vertical force and a small negative horizontal force, whilst in the upstroke, the body angle is large and the relatively small aerodynamic force points forward and slightly downward, giving a positive horizontal force which overcomes the body drag and the negative horizontal force generated in the downstroke. (2) Pitching oscillation of the butterfly body plays an equivalent role of the wing-rotation of many other insects. (3) The body-massspecific power of the model butterfly is 33.3 W/kg, not very different from that of many other insects, e.g., fruitflies and dragonflies.     

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.     

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.     

昆虫飞行的高升力机理

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

High-speed stereo DPIV measurement of wakes of two bat species flying freely in a wind tunnel

Previous studies on wake flow visualization of live animals using DPIV have typically used low repetition rate lasers and 2D imaging. Repetition rates of around 10 Hz allow ~1 image per wingbeat in small birds and bats, and even fewer in insects. To accumulate data representing an entire wingbeat therefore requires the stitching-together of images captured from different wingbeats, and at different locations along the wing span for 3D-construction of wake topologies. A 200 Hz stereo DPIV system has recently been installed in the Lund University wind tunnel facility and the high-frame rate can be used to calculate all three velocity components in a cube, whose third dimension is constructed using the Taylor hypothesis. We studied two bat species differing in body size, Glossophaga soricina and Leptonycteris curasoa. Both species shed a tip vortex during the downstroke that was present well into the upstroke, and a vortex of opposite sign to the tip vortex was shed from the wing root. At the transition between upstroke/downstroke, a vortex loop was shed from each wing, inducing an upwash. Vorticity iso-surfaces confirmed the overall wake topology derived in a previous study. The measured dimensionless circulation, Γ/Uc, which is proportional to a wing section lift coefficient, suggests that unsteady phenomena play a role in the aerodynamics of both species.     

昆虫飞行的空气动力学

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

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.     

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

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

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.     

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.     

Enclosure enhancement of flight performance

We use a potential flow solver to investigate the aerodynamic aspects of flapping flights in enclosed spaces. The enclosure effects are simulated by the method of images. Our study complements previous aerodynamic analyses which considered only the near-ground flight. The present results show that flying in the proximity of an enclosure affects the aerodynamic performance of flapping wings in terms of lift and thrust generation and power consumption. It leads to higher flight efficiency and more than 5% increase of the generation of lift and thrust.     

Robustness of the Leading Edge Vortex on Rotating Wings to Unsteady Perturbations

The vast majority of research works on low aspect ratio rotating wings report that, at high angle of attack, the leading edge vortex that forms on the upper surface of the wing is stable. This ‘trick’ is used by insects and auto-rotating seeds, for example, to achieve the desirable amount of lift. Yet, a few experimental studies suggest that leading edge vortices might be unstable under similar, low Rossby number, conditions. While it is unclear what causes vortex shedding in these studies, the present communication explores the sensitivity of leading edge vortex attachment to perturbations of the rotating speed and demonstrates that shedding can be triggered even for very small perturbations, corresponding to wing tip displacements lower than 1% of the wing chord.

     

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.