The effect of phase angle and wing spacing on tandem flapping wings

In a tandem wing configuration, the hindwing often operates in the wake of the forewing and, hence, its performance is affected by the vortices shed by the forewing. Changes in the phase angle between the flapping motions of the fore and the hind wings, as well as the spacing between them, can affect the resulting vortex/wing and vortex/vortex interactions. This study uses 2D numerical simulations to investigate how these changes affect the leading dege vortexes (LEV) generated by the hindwing and the resulting effect on the lift and thrust coefficients as well as the efficiencies. The tandem wing configuration was simulated using an incompressible Navier-Stokes solver at a chord-based Reynolds number of 5 000. A harmonic single frequency sinusoidal oscillation consisting of a combined pitch and plunge motion was used for the flapping wing kinematics at a Strouhal number of 0.3. Four different spacings ranging from 0.1 chords to 1 chord were tested at three different phase angles, 0°, 90° and 180°. It was found that changes in the spacing and phase angle affected the timing of the interaction between the vortex shed from the forewing and the hindwing. Such an interaction affects the LEV formation on the hindwing and results in changes in aerodynamic force production and efficiencies of the hindwing. It is also observed that changing the phase angle has a similar effect as changing the spacing. The results further show that at different spacings the peak force generation occurs at different phase angles, as do the peak efficiencies.     

Vortical structures on a flapping wing

A wing in the form of a rectangular flat plate is subjected to periodic flapping motion. Space–time imaging provides quantitative representations of the flow structure along the wing. Regions of spanwise flow exist along the wing surface; and depending on the location along the span, the flow is either toward or away from the tip of the wing. Onset and development of large-scale, streamwise-oriented vortical structures occur at locations inboard of the tip of the wing, and they can attain values of circulation of the order of one-half the circulation of the tip vortex. Time-shifted images indicate that these streamwise vortical structures persist over a major share of the wing chord. Space–time volume constructions define the form and duration of these structures, relative to the tip vortex.     

Experimental investigation of a flapping wing model

The main objective of this research study was to investigate the aerodynamic forces of an avian flapping wing model system. The model size and the flow conditions were chosen to approximate the flight of a goose. Direct force measurements, using a three-component balance, and PIV flow field measurements parallel and perpendicular to the oncoming flow, were performed in a wind tunnel at Reynolds numbers between 28,000 and 141,000 (3–15 m/s), throughout a range of reduced frequencies between 0.04 and 0.20. The appropriateness of quasi-steady assumptions used to compare 2D, time-averaged particle image velocimetry (PIV) measurements in the wake with direct force measurements was evaluated. The vertical force coefficient for flapping wings was typically significantly higher than the maximum coefficient of the fixed wing, implying the influence of unsteady effects, such as delayed stall, even at low reduced frequencies. This puts the validity of the quasi-steady assumption into question. The (local) change in circulation over the wing beat cycle and the circulation distribution along the wingspan were obtained from the measurements in the tip and transverse vortex planes. Flow separation could be observed in the distribution of the circulation, and while the circulation derived from the wake measurements failed to agree exactly with the absolute value of the circulation, the change in circulation over the wing beat cycle was in excellent agreement for low and moderate reduced frequencies. The comparison between the PIV measurements in the two perpendicular planes and the direct force balance measurements, show that within certain limitations the wake visualization is a powerful tool to gain insight into force generation and the flow behavior on flapping wings over the wing beat cycle.     

Exploration of the rotational power consumption of a rigid flapping wing

The development of Micro Air Vehicles with flapping wings is inspired from the observation and study of natural flyers such as insects and birds. This article explores the rotational power consumption of a flapping wing using a mechanical flapper at Re ≃ 4,500. This mechanical flapper is simplified to a 2D translation and a rotation in a water tank. Moreover, the wing kinematics are reduced to a linear translation and a rotation for the purpose of our study. We introduce the notion of non-ideal flapper and associated non-ideal rotational power. Such non-ideal devices are defined as consuming power for adding and removing mechanical power to and from the flow, respectively. First, we use a traditional symmetrical wing kinematic which is a simplified kinematic inspired from natural flyers. The lift coefficient of this flapping is about C _L ≃ 1.5. This symmetrical wing kinematic is chosen as a reference. Further, wing kinematics with asymmetric rotations are then compared with this one. These new kinematics are built using a differential velocity defined according to the translational kinematics, a time lag and a distance, r _kp. The analogy of this distance is discussed as a key point to follow along the chord. First, the wing kinematics are varied keeping a similar shape for the profiles of the angular velocity. It is shown that when compared to the reference wing kinematic, a 10% reduction in the rotational power is obtained whilst the lift is reduced by 9%. Second, we release the limitation to a similar shape for the profiles of the angular velocity leading to a novel shape for the angular velocity profile named here as “double bump” profile. With these new wing kinematics, we show that a 60% reduction in the non-ideal rotational power can be achieved whilst the lift coefficient is only reduced by 1.7%. Such “double bump kinematics” could then be of interest to increase the endurance of Micro Air Vehicles.     

Propulsion velocity of a flapping wing at low Reynolds number

This paper presents a computational fluid–structure interaction analysis for free movements with a flapping wing in a quiescent fluid. We demonstrated the moving velocity of a flapping wing according to the phase difference between the angle of attack and the positional angle in the case of a fruit fly with a Reynolds number of 136. If we considered the moving velocity of the flapping wing, the physics were different from that of hovering flight of previous studies, which did not consider the propulsive velocity and presented the advanced rotation of the angle of attack as the best mechanism for propulsion force, as compared to symmetric rotation and delayed rotation. We found that symmetric rotation produced a better propulsion velocity with less fluctuation in other direction than the advanced rotation. The hairpin vortex generated at the end of a stroke did not clearly contribute to the enhancement of propulsion; the wake capture is considered to be one of the main enhancements of the advanced rotation in a previous studies. We studied the effects of the angle of attack to determine why the fruit fly uses a large angle of attack during a constant angle of attack period. Larger angles of attack produced greater propulsion velocities. Further, larger angles of attack did not generate greater peak force during the rotation of the angle of attack at the reversal of stroke, but they produced less fluctuation at the reversal of the stroke and greater force during the constant angle of attack period.     

Numerical investigation on aerodynamic performance of a bionic flapping wing

This paper numerically studies the aerodynamic performance of a bird-like bionic flapping wing. The geometry and kinematics are designed based on a seagull wing,in which flapping, folding, swaying, and twisting are considered. An in-house unsteady flow solver based on hybrid moving grids is adopted for unsteady flow simulations. We focus on two main issues in this study, i.e., the influence of the proportion of down-stroke and the effect of span-wise twisting. Numerical results show that the proportion of downstroke is closely related to the efficiency of the flapping process. The preferable proportion is about 0.7 by using the present geometry and kinematic model, which is very close to the observed data. Another finding is that the drag and the power consumption can be greatly reduced by the proper span-wise twisting. Two cases with different reduced frequencies are simulated and compared with each other. The numerical results show that the power consumption reduces by more than 20%, and the drag coefficient reduces by more than 60% through a proper twisting motion for both cases. The flow mechanism is mainly due to controlling of unsteady flow separation by adjusting the local effective angle of attack. These conclusions will be helpful for the high-performance micro air vehicle(MAV) design.     

Analogy between a flapping wing and a wind turbine with a vertical axis of revolution

Based on an analysis of available experimental data, the hypothesis about an analogy between a flapping wing and a wind turbine of the Darrieus rotor type is justified. It is demonstrated that the torque on the shaft of the Darrieus rotor is generated by thrust forces acting on the blades in a pulsed flow. A conclusion is drawn that it is necessary to perform aerodynamic calculations of blades on the basis of the nonlinear theory of the wing in an unsteady flow with allowance for the airfoil thickness. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 50, No. 2, pp. 152–155, March–April, 2009.     

Towards simulation of flapping wings using immersed boundary method

In this work the immersed boundary method is applied to simulate incompressible turbulent flows around stationary and moving objects. The goal is to demonstrate that the immersed boundary technique along with a large eddy simulation approach is capable of simulating the effect of the so‐called leading edge vortex (LEV), which can be found in flapping wing aerodynamics. A Lagrangian method is used to approximate the solutions in the freshly cleared cells that lay within solid objects at one time step and emerge into fluid domain at the next time step. Flow around a stationary cylinder at Re_D = 20, 40, and 3900 (based on cylinder diameter D) is first studied to validate the immersed boundary solver based on the finite volume scheme using a staggered grid. Then, a harmonically oscillating cylinder at Re_D = 10 000 is considered to test the solver after the Lagrangian method is implemented to interpolate the solution in the freshly cleared cells. Finally, this approach is used to study flows around a stationary flat‐plate at several angles of attack and fast pitching flat‐plate. The rapidly pitching plate creates a dynamic LEV that can be used to improve the efficiency of flapping wings of micro air vehicle and to determine the optimum flapping frequency. Copyright © 2012 John Wiley & Sons, Ltd.     

Time-resolved reconstruction of the full velocity field around a dynamically-scaled flapping wing

The understanding of the physics of flapping flight has long been limited due to the obvious experimental difficulties in studying the flow field around real insects. In this study the time-dependent three-dimensional velocity field around a flapping wing was measured quantitatively for the first time. This was done using a dynamically-scaled wing moving in mineral oil in a pattern based on the kinematics obtained from real insects. The periodic flow is very reproducible, due to the relatively low Reynolds number and precise control of the wing. This repeatability was used to reconstruct the full evolving flow field around the wing from separate stereoscopic particle image velocimetry measurements for a number of spanwise planes and time steps. Typical results for two cases (an impulsive start and a simplified flapping pattern) are reported. Visualizations of the obtained data confirm the general picture of the leading-edge vortex that has been reported in recent publications, but allow a refinement of the detailed structure: rather than a single strand of vorticity, we find a stable pair of counter-rotating structures. We show that the data can also be used for quantitative studies, such as lift and drag prediction.

---

Influence of wing flexibility on the aerodynamic performance of a tethered flapping bumblebee

The sophisticated structures of flapping insect wings make it challenging to study the role of wing flexibility in insect flight. In this study, a mass-spring system is used to model wing structural dynamics as a thin, flexible membrane supported by a network of veins. The vein mechanical properties can be estimated based on their diameters and the Young's modulus of cuticle. In order to analyze the effect of wing flexibility, the Young's modulus is varied to make a comparison between two different wing models that we refer to as flexible and highly flexible. The wing models are coupled with a pseudo-spectral code solving the incompressible Navier–Stokes equations, allowing us to investigate the influence of wing deformation on the aerodynamic efficiency of a tethered flapping bumblebee. Compared to the bumblebee model with rigid wings, the one with flexible wings flies more efficiently, characterized by a larger lift-to-power ratio.     

Experimental study of wings undergoing active root flapping and pitching

This paper presents the results of experiments carried out on mechanical wings undergoing active root flapping and pitching in the wind tunnel. The objective of the work is to investigate the effect of the pitch angle oscillations and wing profile on the aerodynamic forces generated by the wings. The experiments were repeated for a different reduced frequency, airspeed, flapping and pitching kinematics, geometric angle of attack and wing sections (one symmetric and two cambered airfoils). A specially designed mechanical flapper was used, modelled on large migrating birds. It is shown that, under pitch leading conditions, good thrust generation can be obtained at a wide range of Strouhal numbers if the pitch angle oscillation is adjusted accordingly. Consequently, high thrust was measured at both the lowest and highest tested Strouhal numbers. Furthermore, the work demonstrates that the aerodynamic forces can be sensitive to the Reynolds number, depending on the camber of the wings. Under pitch lagging conditions, where the effective angle of attack amplitude is highest, the symmetric wing was affected by the Reynolds number, generating less thrust at the lowest tested Reynolds value. In contrast, under pure flapping conditions, where the effective angle of attack amplitude was lower but still significant, it was the cambered wings that demonstrated Reynolds sensitivity.     

Experimental study of the flow field induced by a resonating piezoelectric flapping wing

Flexible plate structures with integrated piezoelectric patches offer interesting possibilities when considered as actuation mechanisms for energy harvesting devices, cooling devices and propulsion devices of micro-aerial vehicles. Most of the studies reported in literature are based on the assumption of a 2D aerodynamic flow. However, the flow behind a finite span wing is significantly more complex than that of an infinite span wing. In order to corroborate this statement, the present experimental study contains high-speed particle image velocimetry measurements performed on a piezoelectric finite span wing oscillating in air, at 84.8 Hz. The paper focuses on the situation of low Keulegan–Carpenter numbers (KC < 3). The dimensionless KC number describes the relative importance of the drag forces over inertia forces for objects that oscillate in a fluid flow at rest. The evolution of the unsteady vortex structures near the plate is characterized for different conditions. This allows a better understanding of the unsteady aerodynamics of flapping flight. The accomplished experimental data analysis has shown that the flow phenomena are strongly dependent on the KC values.     

Flow visualization and wall shear stress of a flapping model hummingbird wing

The unsteady low Reynolds number aerodynamics of flapping flight was investigated experimentally through flow visualization by suspended particle imagery and wall shear stress measurement from micro-array hot-film anemometry. In conjunction, a mechanism was developed to create a flapping motion with three degrees of freedom and adjustable flapping frequency. The flapping kinematics and wing shape were selected for dynamic similarity to a hummingbird during hovering flight. Flow visualization was used to validate the anemometry observations of leading edge vortex (LEV) characteristics and to investigate the necessity of spanwise flow in LEV stability. The shear sensors determined LEV characteristics throughout the translation section of the stroke period for various wing speeds. It was observed that a minimum frequency between 2 and 3.5 Hz is required for the formation and stabilization of a LEV. The vortex strength peaked around 30% of the flapping cycle (corresponding to just past the translation midpoint), which agrees with results from previous studies conducted by others. The shear sensors also indicated a mild growth in LEV size during translation sections of the wing’s motion. This growth magnitude was nearly constant through a range of operating frequencies.     

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     

Aeroelastic study of flexible flapping wings by a coupled lattice Boltzmann-finite element approach with immersed boundary method

In this paper, the behavior of two-dimensional symmetric flapping wings moving in a viscous fluid is investigated. Harmonic motion is applied to idealize flying organisms with flexible wings and extensive testing is carried out to investigate the resultant flight behavior related to the ability to take-off or accelerate the flapping wing system away from a starting location. Special attention is paid to analyze the effect of the main mechanical parameters, as well as the effect of lateral wind on flight performances. Moreover, aiming to investigate the possible benefits of flying in flocks, a couple of synchronously flapping wings is considered in addition to the single arrangement. The numerical simulations are performed by solving the fluid–structure interaction problem through a strongly coupled partitioned approach. Fluid dynamics are modeled at the mesoscopic scale by the lattice Boltzmann method. The resulting macroscopic quantities are derived, as usual, based on the statistical molecular-level interpretation.Wings are modeled by geometrically nonlinear, elastic beam finite elements and structure dynamics is solved by the time discontinuous Galerkin method. Fluid–structure interface conditions are handled using the immersed boundary method. The resultant numerical approach combines simplicity and high computational efficiency. A Monte Carlo simulation strategy is employed to characterize the flight behavior subjected to lateral wind. Various scenarios are discussed.     

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)     

Experimental investigation on aerodynamic performance of a flapping wing vehicle in forward flight

The aerodynamic performance of a flexible membrane flapping wing has been investigated here. For this purpose, a flapping-wing system and an experimental set-up were designed to measure the unsteady aerodynamic forces of the flapping wing motion. A one-component force balance was set up to record the temporal variations of aerodynamic forces. The flapping wing was studied in a large low-speed wind tunnel. The lift and thrust of this mechanism were measured for different flapping frequencies, angles of attack and for various wind tunnel velocities. Results indicate that the thrust increases with the flapping frequency. An increase in the wind tunnel speed and flow angle of attack leads to reduction in the thrust value and increases the lift component. The aerodynamic and performance parameters were nondimensionalized. Appropriate models were introduced which show its aerodynamic performance and may be used in the design process and also optimization of the flapping wing.     

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.