Propulsive performance and vortex shedding of a foil in flapping flight

Summary. The propulsive performance and vortex shedding of an oscillating foil, which mimics biological locomotion, are investigated based on a computational fluid dynamics analysis. The objectives of this study are to investigate unsteady forces, in particular a thrust force, for the foil in pitching and plunging motion, and to deal with the relations of the propulsive performance with leading-edge vortex structure and vortex shedding in the near wake. The two-dimensional incompressible Navier–Stokes equations in the vorticity and stream-function formulation are solved by fourth-order essentially compact finite difference schemes for the space derivatives and a fourth-order Runge-Kutta scheme for the time advancement. To reveal the mechanism of the propulsive performance, the unsteady forces and the shedding of the leading- and trailing- edge vortices of the foil in the pitching and plunging motion are analyzed. Based on our calculated results, three types of the leading-edge vortex shedding, which have an effective influence on the vortex structures in the wake of the oscillating foil, are identified. The effects of some typical factors, such as the frequency and amplitude of the oscillation, the phase difference between the pitching and plunging motions, and the thickness ratio of the foil, on the vortex shedding and the unsteady forces are discussed.     

Three-dimensional flow and lift characteristics of a hovering ruby-throated hummingbird

A three-dimensional computational fluid dynamics simulation is performed for a ruby-throated hummingbird (Archilochus colubris) in hovering flight. Realistic wing kinematics are adopted in the numerical model by reconstructing the wing motion from high-speed imaging data of the bird. Lift history and the three-dimensional flow pattern around the wing in full stroke cycles are captured in the simulation. Significant asymmetry is observed for lift production within a stroke cycle. In particular, the downstroke generates about 2.5 times as much vertical force as the upstroke, a result that confirms the estimate based on the measurement of the circulation in a previous experimental study. Associated with lift production is the similar power imbalance between the two half strokes. Further analysis shows that in addition to the angle of attack, wing velocity and surface area, drag-based force and wing–wake interaction also contribute significantly to the lift asymmetry. Though the wing–wake interaction could be beneficial for lift enhancement, the isolated stroke simulation shows that this benefit is buried by other opposing effects, e.g. presence of downwash. The leading-edge vortex is stable during the downstroke but may shed during the upstroke. Finally, the full-body simulation result shows that the effects of wing–wing interaction and wing–body interaction are small.     

The wake of hovering flight in bats

Hovering means stationary flight at zero net forward speed, which can be achieved by animals through muscle powered flapping flight. Small bats capable of hovering typically do so with a downstroke in an inclined stroke plane, and with an aerodynamically active outer wing during the upstroke. The magnitude and time history of aerodynamic forces should be reflected by vorticity shed into the wake. We thus expect hovering bats to generate a characteristic wake, but this has until now never been studied. Here we trained nectar-feeding bats, Leptonycteris yerbabuenae, to hover at a feeder and using time-resolved stereoscopic particle image velocimetry in conjunction with high-speed kinematic analysis we show that hovering nectar-feeding bats produce a series of bilateral stacked vortex loops. Vortex visualizations suggest that the downstroke produces the majority of the weight support, but that the upstroke contributes positively to the lift production. However, the relative contributions from downstroke and upstroke could not be determined on the basis of the wake, because wake elements from down- and upstroke mix and interact. We also use a modified actuator disc model to estimate lift force, power and flap efficiency. Based on our quantitative wake-induced velocities, the model accounts for weight support well (108%). Estimates of aerodynamic efficiency suggest hovering flight is less efficient than forward flapping flight, while the overall energy conversion efficiency (mechanical power output/metabolic power) was estimated at 13%.     

Computational investigation of cicada aerodynamics in forward flight

Free forward flight of cicadas is investigated through high-speed photogrammetry, three-dimensional surface reconstruction and computational fluid dynamics simulations. We report two new vortices generated by the cicada''s wide body. One is the thorax-generated vortex, which helps the downwash flow, indicating a new phenomenon of lift enhancement. Another is the cicada posterior body vortex, which entangles with the vortex ring composed of wing tip, trailing edge and wing root vortices. Some other vortex features include: independently developed left- and right-hand side leading edge vortex (LEV), dual-core LEV structure at the mid-wing region and near-wake two-vortex-ring structure. In the cicada forward flight, approximately 79% of the total lift is generated during the downstroke. Cicada wings experience drag in the downstroke, and generate thrust during the upstroke. Energetics study shows that the cicada in free forward flight consumes much more power in the downstroke than in the upstroke, to provide enough lift to support the weight and to overcome drag to move forward.     

Insect-like flapping wing mechanism based on a double spherical Scotch yoke

We describe the rationale, concept, design and implementation of a fixed-motion (non-adjustable) mechanism for insect-like flapping wing micro air vehicles in hover, inspired by two-winged flies (Diptera). This spatial (as opposed to planar) mechanism is based on the novel idea of a double spherical Scotch yoke. The mechanism was constructed for two main purposes: (i) as a test bed for aeromechanical research on hover in flapping flight, and (ii) as a precursor design for a future flapping wing micro air vehicle. Insects fly by oscillating (plunging) and rotating (pitching) their wings through large angles, while sweeping them forwards and backwards. During this motion the wing tip approximately traces a "figure-of-eight" or a "banana" and the wing changes the angle of attack (pitching) significantly. The kinematic and aerodynamic data from free-flying insects are sparse and uncertain, and it is not clear what aerodynamic consequences different wing motions have. Since acquiring the necessary kinematic and dynamic data from biological experiments remains a challenge, a synthetic, controlled study of insect-like flapping is not only of engineering value, but also of biological relevance. Micro air vehicles are defined as flying vehicles approximately 150 mm in size (hand-held), weighing 50-100g, and are developed to reconnoitre in confined spaces (inside buildings, tunnels, etc.). For this application, insect-like flapping wings are an attractive solution and hence the need to realize the functionality of insect flight by engineering means. Since the semi-span of the insect wing is constant, the kinematics are spatial; in fact, an approximate figure-of-eight/banana is traced on a sphere. Hence a natural mechanism implementing such kinematics should be (i) spherical and (ii) generate mathematically convenient curves expressing the figure-of-eight/banana shape. The double spherical Scotch yoke design has property (i) by definition and achieves (ii) by tracing spherical Lissajous curves.     

Prediction of leading-edge vortex behaviour to supplement the suction analogy

Summary An analytical treatment is presented which permits the prediction of the strength and path of leading-edge vortices on thin wings of delta-like planform from a knowledge of the behaviour of the linearized approximation to the attached flow past the wing. This supplements the leading-edge suction analogy, which predicts the forces and moments acting on the wing from the same input. The fundamental assumption is that the vortex lies close to the leading edge. It is represented by a single line-vortex and the presentation is confined to low-speed flow and plane wings. The treatment is independent of the way the basic attached flow is calculated. Results are shown for two simple planforms, using inputs from three-dimensional lifting-surface theory and from slender-body theory.     

Vortex wake,downwash distribution,aerodynamic performance and wingbeat kinematics in slow-flying pied flycatchers

Many small passerines regularly fly slowly when catching prey, flying in cluttered environments or landing on a perch or nest. While flying slowly, passerines generate most of the flight forces during the downstroke, and have a ‘feathered upstroke’ during which they make their wing inactive by retracting it close to the body and by spreading the primary wing feathers. How this flight mode relates aerodynamically to the cruising flight and so-called ‘normal hovering’ as used in hummingbirds is not yet known. Here, we present time-resolved fluid dynamics data in combination with wingbeat kinematics data for three pied flycatchers flying across a range of speeds from near hovering to their calculated minimum power speed. Flycatchers are adapted to low speed flight, which they habitually use when catching insects on the wing. From the wake dynamics data, we constructed average wingbeat wakes and determined the time-resolved flight forces, the time-resolved downwash distributions and the resulting lift-to-drag ratios, span efficiencies and flap efficiencies. During the downstroke, slow-flying flycatchers generate a single-vortex loop wake, which is much more similar to that generated by birds at cruising flight speeds than it is to the double loop vortex wake in hovering hummingbirds. This wake structure results in a relatively high downwash behind the body, which can be explained by the relatively active tail in flycatchers. As a result of this, slow-flying flycatchers have a span efficiency which is similar to that of the birds in cruising flight and which can be assumed to be higher than in hovering hummingbirds. During the upstroke, the wings of slowly flying flycatchers generated no significant forces, but the body–tail configuration added 23 per cent to weight support. This is strikingly similar to the 25 per cent weight support generated by the wing upstroke in hovering hummingbirds. Thus, for slow-flying passerines, the upstroke cannot be regarded as inactive, and the tail may be of importance for flight efficiency and possibly manoeuvrability.     

Hummingbird wing efficacy depends on aspect ratio and compares with helicopter rotors

Hummingbirds are the only birds that can sustain hovering. This unique flight behaviour comes, however, at high energetic cost. Based on helicopter and aeroplane design theory, we expect that hummingbird wing aspect ratio (AR), which ranges from about 3.0 to 4.5, determines aerodynamic efficacy. Previous quasi-steady experiments with a wing spinner set-up provide no support for this prediction. To test this more carefully, we compare the quasi-steady hover performance of 26 wings, from 12 hummingbird taxa. We spun the wings at angular velocities and angles of attack that are representative for every species and measured lift and torque more precisely. The power (aerodynamic torque × angular velocity) required to lift weight depends on aerodynamic efficacy, which is measured by the power factor. Our comparative analysis shows that AR has a modest influence on lift and drag forces, as reported earlier, but interspecific differences in power factor are large. During the downstroke, the power required to hover decreases for larger AR wings at the angles of attack at which hummingbirds flap their wings (p < 0.05). Quantitative flow visualization demonstrates that variation in hover power among hummingbird wings is driven by similar stable leading edge vortices that delay stall during the down- and upstroke. A side-by-side aerodynamic performance comparison of hummingbird wings and an advanced micro helicopter rotor shows that they are remarkably similar.     

Wake structure and wingbeat kinematics of a house-martin Delichon urbica.

The wingbeat kinematics and wake structure of a trained house martin in free, steady flight in a wind tunnel have been studied over a range of flight speeds, and compared and contrasted with similar measurements for a thrush nightingale and a pair of robins. The house martin has a higher aspect ratio (more slender) wing, and is a more obviously agile and aerobatic flyer, catching insects on the wing. The wingbeat is notable for the presence at higher flight speeds of a characteristic pause in the upstroke. The essential characteristics of the wing motions can be reconstructed with a simple two-frequency model derived from Fourier analysis. At slow speeds, the distribution of wake vorticity is more simple than for the other previously measured birds, and the upstroke does not contribute to weight support. The upstroke becomes gradually more significant as the flight speed increases, and although the vortex wake shows a signature of the pause phase, the global circulation measurements are otherwise in good agreement with surprisingly simple aerodynamic models, and with predictions across the different species, implying quite similar aerodynamic performance of the wing sections. The local Reynolds numbers of the wing sections are sufficiently low that the well-known instabilities of attached laminar flows over lifting surfaces, which are known to occur at two to three times this value, may not develop.     

Time-varying span efficiency through the wingbeat of desert locusts

The flight performance of animals depends greatly on the efficacy with which they generate aerodynamic forces. Accordingly, maximum range, load-lifting capacity and peak accelerations during manoeuvres are all constrained by the efficiency of momentum transfer to the wake. Here, we use high-speed particle image velocimetry (1 kHz) to record flow velocities in the near wake of desert locusts (Schistocerca gregaria, Forskål). We use the measured flow fields to calculate time-varying span efficiency throughout the wing stroke cycle. The locusts are found to operate at a maximum span efficiency of 79 per cent, typically at a plateau of about 60 per cent for the majority of the downstroke, but at lower values during the upstroke. Moreover, the calculated span efficiencies are highest when the largest lift forces are being generated (90% of the total lift is generated during the plateau of span efficiency) suggesting that the combination of wing kinematics and morphology in locust flight perform most efficiently when doing the most work.     

Parallel simulation of unsteady hovering rotor wakes

Numerical simulation using low diffusion schemes, for example free‐vortex or vorticity transport methods, and theoretical stability analyses have shown the wakes of rotors in hover to be unsteady. This has also been observed in experiments, although the instabilities are not always repeatable. Hovering rotor wake stability is considered here using a finite‐volume compressible CFD code. An implicit unsteady, multiblock, multigrid, upwind solver, and structured multiblock grid generator are presented, and applied to lifting rotors in hover. To allow the use of very fine meshes and, hence, better representation of the flow physics, a parallel version of the code has been developed, and parallel performance using upto 1024 CPUs is presented. A four‐bladed rotor is considered, and it is demonstrated that once the grid density is sufficient to capture enough turns of the tip vortices, hover exhibits oscillatory behaviour of the wake, even using a steady formulation. An unsteady simulation is then performed, and also shows an unsteady wake. Detailed analysis of the time‐accurate wake history shows that three dominant unsteady modes are captured, for this four‐bladed case, with frequencies of one, four, and eight times the rotational frequency. A comparison with theoretical stability analysis is also presented. Copyright © 2006 John Wiley & Sons, Ltd.     

Three-dimensional vortex wake structure of flapping wings in hovering flight

Flapping wings continuously create and send vortices into their wake, while imparting downward momentum into the surrounding fluid. However, experimental studies concerning the details of the three-dimensional vorticity distribution and evolution in the far wake are limited. In this study, the three-dimensional vortex wake structure in both the near and far field of a dynamically scaled flapping wing was investigated experimentally, using volumetric three-component velocimetry. A single wing, with shape and kinematics similar to those of a fruitfly, was examined. The overall result of the wing action is to create an integrated vortex structure consisting of a tip vortex (TV), trailing-edge shear layer (TESL) and leading-edge vortex. The TESL rolls up into a root vortex (RV) as it is shed from the wing, and together with the TV, contracts radially and stretches tangentially in the downstream wake. The downwash is distributed in an arc-shaped region enclosed by the stretched tangential vorticity of the TVs and the RVs. A closed vortex ring structure is not observed in the current study owing to the lack of well-established starting and stopping vortex structures that smoothly connect the TV and RV. An evaluation of the vorticity transport equation shows that both the TV and the RV undergo vortex stretching while convecting downwards: a three-dimensional phenomenon in rotating flows. It also confirms that convection and secondary tilting and stretching effects dominate the evolution of vorticity.     

Boundary element simulation of inviscid flow around an oscillatory foil with vortex sheet

The dynamic performance of a rigid foil with harmonic vertical and rotational motions in fluid flow has been studied through velocity potential theory. A boundary element based time stepping scheme is introduced to simulate the flow around the foil and the vortex wake. The body surface condition is satisfied on the exact foil surface and the motion and deformation of the wake sheet shed at the trailing edge is tracked. Kelvin condition is satisfied and a Kutta condition for the unsteady motion is proposed to circumvent the singularity at the trailing edge. Point vortex, which is reduced from wake vortex dipole, is introduced to approximate the vorticity. The performance of foil NACA0012 with harmonic vertical and rotational motions are studied extensively; the propulsion/swimming mode, energy harvesting mode and the flying mode are analysed in detail.     

Deformable wing kinematics in free-flying hoverflies

Here, we present a detailed analysis of the deforming wing kinematics of free-flying hoverflies (Eristalis tenax, Linnaeus) during hovering flight. We used four high-speed digital video cameras to reconstruct the motion of approximately 22 points on each wing using photogrammetric techniques. While the root-flapping motion of the wing is similar in both the downstroke and upstroke, and is well modelled as a simple harmonic motion, other wing kinematic parameters show substantial variation between the downstroke and upstroke. Whereas the magnitude of the angle of incidence varies considerably within and between different hoverflies, the twist distribution along the wing is highly stereotyped. The angle of incidence and camber both show a recoil effect as they change abruptly at stroke reversal. Pronation occurs consistently after stroke reversal, which is perhaps surprising, because this has been found to reduce lift production in modelling studies. We find that the alula, a hinged flap near the base of the wing, operates in two discrete states: either in plane with the wing, or flipped approximately normal to it. We hypothesize that the alula may be acting as a flow-control device.     

Deformable wing kinematics in the desert locust: how and why do camber,twist and topography vary through the stroke?

Here, we present a detailed analysis of the wing kinematics and wing deformations of desert locusts (Schistocerca gregaria, Forskål) flying tethered in a wind tunnel. We filmed them using four high-speed digital video cameras, and used photogrammetry to reconstruct the motion of more than 100 identified points. Whereas the hindwing motions were highly stereotyped, the forewing motions showed considerable variation, consistent with a role in flight control. Both wings were positively cambered on the downstroke. The hindwing was cambered through an ‘umbrella effect’ whereby the trailing edge tension compressed the radial veins during the downstroke. Hindwing camber was reversed on the upstroke as the wing fan corrugated, reducing the projected area by 30 per cent, and releasing the tension in the trailing edge. Both the wings were strongly twisted from the root to the tip. The linear decrease in incidence along the hindwing on the downstroke precisely counteracts the linear increase in the angle of attack that would otherwise occur in root flapping for an untwisted wing. The consequent near-constant angle of attack is reminiscent of the optimum for a propeller of constant aerofoil section, wherein a linear twist distribution allows each section to operate at the unique angle of attack maximizing the lift to drag ratio. This implies tuning of the structural, morphological and kinematic parameters of the hindwing for efficient aerodynamic force production.     

Elytra boost lift,but reduce aerodynamic efficiency in flying beetles

Flying insects typically possess two pairs of wings. In beetles, the front pair has evolved into short, hardened structures, the elytra, which protect the second pair of wings and the abdomen. This allows beetles to exploit habitats that would otherwise cause damage to the wings and body. Many beetles fly with the elytra extended, suggesting that they influence aerodynamic performance, but little is known about their role in flight. Using quantitative measurements of the beetle''s wake, we show that the presence of the elytra increases vertical force production by approximately 40 per cent, indicating that they contribute to weight support. The wing-elytra combination creates a complex wake compared with previously studied animal wakes. At mid-downstroke, multiple vortices are visible behind each wing. These include a wingtip and an elytron vortex with the same sense of rotation, a body vortex and an additional vortex of the opposite sense of rotation. This latter vortex reflects a negative interaction between the wing and the elytron, resulting in a single wing span efficiency of approximately 0.77 at mid downstroke. This is lower than that found in birds and bats, suggesting that the extra weight support of the elytra comes at the price of reduced efficiency.     

Numerical simulation of incompressible Navier-Stokes and Euler equations to the vortical flow about a delta wing

Summary An algorithm based on the method of artificial compressibility is used to solve the incompressible Navier-Stokes and Euler equations in 3D generalized curvilinear coordinates. The choice of the pseudocompressibility parameter is discussed and a function form of which accelerates the convergence effectively is presented. Numerical simulation of the vortical flow about a 70° delta wing is made by both Euler and Navier-Stokes methods. Comparison of the computed results and the experimental data indicates that the Euler method can capture the global flow feature, i.e., separation from the sharp leading edge and the formation of the leading-edge vortex; the Navier-Stokes method depicts not only the flow features captured by the Euler method, but also the boundary layer effects and specifically the secondary separation induced by the leading-edge vortex on the wing upper surface. Thus, the Navier-Stokes method shows a much improved correlation with the experimental data.     

Effect of flapping kinematics on aerodynamic force of a flapping two-dimensional flat plate

Potential applications of flapping-wing micro-aerial vehicles (MAVs) have prompted enthusiasm among the engineers and researchers to understand the flow physics associated with flapping flight. An incompressible Navier–Stokes solver that is capable of handling flapping flight kind of moving boundary problem is developed. Arbitrary Lagrangian–Eulerian (ALE) method is used to handle the moving boundaries of the problem. The solver is validated with the results of problems like inline oscillation of a circular cylinder in still fluid and a flat plate rapidly accelerating at constant angle of attack. Numerical simulations of flapping flat plate mimicking the kinematics of those like insect wings are simulated, and the unsteady fluid dynamic phenomena that enhance the aerodynamic force are studied. The solution methodology provides the velocity field and pressure field details, which are used to derive the force coefficients and the vorticity field. Time history of force coefficients and vortical structures gives insight into the unsteady mechanism associated with the unsteady aerodynamic force production. The scope of the work is to develop a computational fluid dynamic (CFD) solver with the ALE method that is capable of handling moving boundary problems, and to understand the flow physics associated with the flapping-wing aerofoil kinematics and flow parameters on aerodynamic forces. Results show that delayed stall, wing–wake interaction and rotational effect are the important unsteady mechanisms that enhance the aerodynamic forces. Major contribution to the lift force is due to the presence of leading edge vortex in delayed stall mechanism.     

柔性扑翼气动性能的数值研究

国内外对扑翼飞行的气动特性进行了大量研究,这些研究大多基于简谐扑动的刚性翼,然而大量观察发现鸟或昆虫飞行时,翅膀存在明显的柔性变形,这种变形对其气动性能具有显著的影响。该文针对一简化的二维柔性扑翼模型,采用数值求解N-S方程并耦合扑翼柔性变形方程的计算方法,研究了扑翼柔性变形对其气动性能的影响。结果显示扑翼的柔性变形改变了扑翼周围的涡结构,从而影响扑翼的气动性能;适当的柔性变形能延迟前缘涡的脱落,从而有效地改善扑翼的推进效率,但同时减弱了扑翼在低雷诺数环境中产生高升力的尾迹捕捉机制。     

Power reduction and the radial limit of stall delay in revolving wings of different aspect ratio

Airplanes and helicopters use high aspect ratio wings to reduce the power required to fly, but must operate at low angle of attack to prevent flow separation and stall. Animals capable of slow sustained flight, such as hummingbirds, have low aspect ratio wings and flap their wings at high angle of attack without stalling. Instead, they generate an attached vortex along the leading edge of the wing that elevates lift. Previous studies have demonstrated that this vortex and high lift can be reproduced by revolving the animal wing at the same angle of attack. How do flapping and revolving animal wings delay stall and reduce power? It has been hypothesized that stall delay derives from having a short radial distance between the shoulder joint and wing tip, measured in chord lengths. This non-dimensional measure of wing length represents the relative magnitude of inertial forces versus rotational accelerations operating in the boundary layer of revolving and flapping wings. Here we show for a suite of aspect ratios, which represent both animal and aircraft wings, that the attachment of the leading edge vortex on a revolving wing is determined by wing aspect ratio, defined with respect to the centre of revolution. At high angle of attack, the vortex remains attached when the local radius is shorter than four chord lengths and separates outboard on higher aspect ratio wings. This radial stall limit explains why revolving high aspect ratio wings (of helicopters) require less power compared with low aspect ratio wings (of hummingbirds) at low angle of attack and vice versa at high angle of attack.