Numerical simulation of a flapping four-wing micro-aerial vehicle

A three-dimensional numerical simulation of a four-wing (two wings on each side, one on top of another) flapping micro-aerial vehicle (FMAV), known as the Delfly micro, is performed using an immersed boundary method Navier–Stokes finite volume solver at Reynolds numbers of 5500 (forward flight condition). The objective of the present investigation is to gain an insight to the aerodynamics of flapping wing biplane configuration, by making an analysis on a geometry that is simplified, yet captures the major aspects of the wing behavior. The fractional step method is used to solve the Navier–Stokes equations. Results show that in comparison to the Delfly II flapping kinematics (a similar FMAV configuration but smaller flapping stroke angles), the Delfly-Micro flapping kinematics provides more thrust while maintaining the same efficiency. The Delfly-Micro biplane configuration generates more lift than expected when the inclination angle increases, due to the formation of a uniform leading edge vortex. Estimates of the lift produced in the forward flight conditions confirm that in the current design, the MAV is able to sustain forward flight. The potential effect of wing flexibility on the aerodynamic performance in the biplane configuration context is investigated through prescribed flexibility in the simulations. Increasing the wing׳ spanwise flexibility increases thrust but increasing chordwise flexibility causes thrust to first increase and then decrease. Moreover, combining both spanwise and chordwise flexibility outperforms cases with only either spanwise or chordwise flexibility.     

Lift enhancement through flexibility of plunging wings at low Reynolds numbers

In this paper the combined effect of two mechanisms for lift enhancement at low Reynolds numbers are considered, wing oscillations and wing flexibility. The force, deformation and flow fields of rigid and flexible low aspect ratio (AR=3) and high aspect ratio (AR=6) wings oscillating at a fixed post-stall angle of attack of 15° and amplitude of 15% of chord are measured. The force measurements show that flexibility can increase the time-averaged lift coefficient significantly. For low aspect ratio wings the maximum lift coefficient across all Strouhal numbers was C_l=1.38 for the rigid wing as opposed to C_l=2.77 for the flexible wing. Very similar trends were observed for the high aspect ratio wings. This increase is associated with significant deformation of the wing. The root is sinusoidally plunged with small amplitude but this motion is amplified along the span resulting in a larger tip motion but with a phase lag. The amount it is amplified strongly depends on Strouhal number. A Strouhal number of Sr_c=1.5 was selected for detailed flow field measurements due to it being central to the high-lift region of the flexible wings, producing approximately double the lift of the rigid wing. For this Strouhal number the rigid wings exhibit a Leading Edge Vortex (LEV) ring. This is where the clockwise upper-surface LEV pairs with the counter-clockwise lower-surface LEV to form a vortex ring that self-advects upstream and away from the wing's upper surface. Conversely the deformation of the flexible wings inhibits the formation of the LEV ring. Instead a strong upper-surface LEV forms during the downward motion and convects close to the airfoil upper surface thus explaining the significantly higher lift. These measurements demonstrate the significant gains that can be achieved through the combination of unsteady aerodynamics with flexible structures.     

Analysis of non-symmetrical flapping airfoils

Simulations have been done to assess the lift, thrust and propulsive efficiency of different types of non-symmetrical airfoils under different flapping configurations. The variables involved are reduced frequency, Strouhal number, pitch amplitude and phase angle. In order to analyze the variables more efficiently, the design of experiments using the response surface methodology is applied. Results show that both the variables and shape of the airfoil have a profound effect on the lift, thrust, and efficiency. By using non- symmetrical airfoils, average lift coefficient as high as 2.23 can be obtained. The average thrust coefficient and efficiency also reach high values of 2.53 and 0.61, respectively. The lift production is highly dependent on the airfoil＇s shape while thrust production is influenced more heavily by the variables. Efficiency falls somewhere in between. Two-factor interac- tions are found to exist among the variables. This shows that it is not sufficient to analyze each variable individually. Vorticity diagrams are analyzed to explain the results obtained. Overall, the S1020 airfoil is able to provide relatively good efficiency and at the same time generate high thrust and lift force. These results aid in the design of a better ornithopter＇s wing.     

Investigation of vortex development on accelerating spanwise-flexible wings

The effect of spanwise flexibility on the development of leading-edge vortices for impulsively started flat plates at low Reynolds numbers has been investigated. A theoretical model is proposed, based on the Euler–Bernoulli beam theory, coupled with a vortex growth model based on vorticity flux through a leading-edge shear layer. The model was validated for rigid and flexible flat plates undergoing a towing motion at an angle-of-attack of 45°. It is shown that a phase-delay in lift and drag generation occurs between rigid and flexible cases. The model indicates decreasing vorticity along the span as the wing is accelerated and begins to bend. Particle image velocimetry measurements conducted at three different spanwise planes showed a delay in vortex growth along the span, despite a uniform spanwise circulation. This uniform spanwise distribution of circulation is in contrast to the quasi-two-dimensional model, which predicted a reduced circulation near the profile tip where plate motion was delayed. It is therefore concluded that circulation must be dynamically redistributed through vorticity convection during the impulsive motion.     

Energy harvesting of a heaving and forward pitching wing with a passively actuated trailing edge

This study experimentally investigates the energy harvesting capabilities of an oscillating wing with a passively actuated trailing edge. The oscillation kinematics are composed of a combined heaving and forward pitching motions, where the pitching axis is well behind the wing center of mass. Passive actuation is attained by connecting the trailing edge with the wing body using a torsion rod. The degree of flexibility of the trailing edge is represented by the Strouhal number based on the trailing edge natural frequency. The trailing edge passive response is studied for oscillation Strouhal numbers of 0.017, 0.025 and 0.033. Instantaneous aerodynamic forces are measured in a closed loop wind tunnel at a Reynolds number of 40 000, based on the free stream velocity and the wing chord length. Measured results include the effective angle of attack induced by the trailing edge actuation as well as the lift and moment during the oscillation cycle. For the imposed kinematics in this study, the pitching motion has a positive contribution to the mean power output whereas the heaving motion has a relatively small but negative contribution. Additionally, by decreasing the natural frequency of the trailing edge closer to that of the imposed oscillation frequency, the magnitude of the lift and moment forces and hence the mean power output, increases. It is found that there exists a strong correlation between mean power output and the effective angle of attack, shown through the passive trailing edge response, resulting in an increase in energy harvesting potential.     

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 investigation on tip vortices and aerodynamics

The tip vortices and aerodynamics of a NACA0012 wing in the vicinity of the ground were studied in a wind tunnel. The wing tip vortex structures and lift/drag forces were measured by a seven-hole probe and a force balance, respectively. The evolution of the flow structures and aerodynamics with a ground height were analyzed. The vorticity of tip vortices was found to reduce with the decreasing of the ground height, and the position of vortex-core moved gradually to the outboard of the wing tip. Therefore, the down-wash flow induced by the tip vortices was weakened. However, vortex breakdown occurred as the wing lowered to the ground. From the experimental results of aerodynamics, the maximum lift-to-drag ratio was observed when the angle of attack was 2.5° and the ground clearance was 0.2.     

Near-field tip vortex behind a swept wing model

The near-field flow structure of a tip vortex behind a sweptback and tapered NACA 0015 wing was investigated and compared with a rectangular wing at the same lift force and Re=1.81×10⁵. The tangential velocity decreased with the downstream distance while increased with the airfoil incidence. The core radius was about 3% of the root chord c _r, regardless of the downstream distance and α for α<8°. The core axial velocity was always wake-like. The core Γ_c and total Γ_o circulation of the tip vortex remained nearly constant up to x/c _r=3.5 and had a Γ_c/Γ_o ratio of 0.63. The total circulation of the tip vortex accounted for only about 40% of the bound root circulation Γ_b. For a rectangular wing, the axial flow exhibited islands of wake- and jet-like velocity distributions with Γ_c/Γ_o=0.75 and Γ_o/Γ_b=0.70. For the sweptback and tapered wing tested, the inner region of the tip vortex flow exhibited a self-similar behavior for x/c _r≥1.0. The lift force computed from the spanwise circulation distributions agreed well with the force-balance data. A large difference in the lift-induced drag was, however, observed between the wake integral method and the inviscid lifting-line theory.     

Large aerodynamic forces on a sweeping wing at low Reynolds number

The aerodynamic forces and flow structure of a model insect wing is studied by solving the Navier-Stokes equations numerically. After an initial start from rest, the wing is made to execute an azimuthal rotation (sweeping) at a large angle of attack and constant angular velocity. The Reynolds number (Re) considered in the present note is 480 (Re is based on the mean chord length of the wing and the speed at 60% wing length from the wing root). During the constant-speed sweeping motion, the stall is absent and large and approximately constant lift and drag coefficients can be maintained. The mechanism for the absence of the stall or the maintenance of large aerodynamic force coefficients is as follows. Soon after the initial start, a vortex ring, which consists of the leading-edge vortex (LEV), the starting vortex, and the two wing-tip vortices, is formed in the wake of the wing. During the subsequent motion of the wing, a base-to-tip spanwise flow converts the vorticity in the LEV to the wing tip and the LEV keeps an approximately constant strength. This prevents the LEV from shedding. As a result, the size of the vortex ring increases approximately linearly with time, resulting in an approximately constant time rate of the first moment of vorticity, or approximately constant lift and drag coefficients. The variation of the relative velocity along the wing span causes a pressure gradient along the wingspan. The base-to-tip spanwise flow is mainly maintained by the pressure-gradient force. The project supported by the National Natural Science Foundation of China (10232010)     

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.     

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

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

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.     

Role of wing morphing in thrust generation

In this paper, we investigate the role of morphing on flight dynamics of two birds by simulating the flow over rigid and morphing wings that have the characteristics of two different birds, namely the Giant Petrel and Dove Prion. The simulation of a flapping rigid wing shows that the root of the wing should be placed at a specific angle of attack in order to generate enough lift to balance the weight of the bird. However, in this case the generated thrust is either very small, or even negative, depending on the wing shape. Further, results show that morphing of the wing enables a significant increase in the thrust and propulsive efficiency. This indicates that the birds actually utilize some sort of active wing twisting and bending to produce enough thrust. This study should facilitate better guidance for the design of flapping air vehicles.     

Forces on oscillating foils for propulsion and maneuvering

Experiments were performed on an oscillating foil to assess its performance in producing large forces for propulsion and effective maneuvering. First, experiments on a harmonically heaving and pitching foil were performed to determine its propulsive efficiency under conditions of significant thrust production, as function of the principal parameters: the heave amplitude, Strouhal number, angle of attack, and phase angle between heave and pitch. Planform area thrust coefficients of 2.4 were recorded for 35° maximum angle of attack and efficiencies of up to 71.5% were recorded for 15° maximum angle of attack. A plateau of good efficiency, in the range of 50–60%, is noted. A phase angle of 90–100° between pitch and heave is found to produce the best thrust performance. Also, the introduction of higher harmonics in the heave motion, so as to ensure a sinusoidal variation in the angle of attack produced much higher thrust coefficient at high Strouhal numbers. Second, experiments on a harmonically oscillating foil with a superposed pitch bias, as well as experiments on impulsively moving foils in still water, were conducted to assess the capability of the foil to produce large lateral forces for maneuvering. Mean side force coefficients of up to 5.5, and instantaneous lift coefficients of up to 15 were recorded, demonstrating an outstanding capability for maneuvering force production.     

A parametric investigation of the propulsion of 2D chordwise-flexible flapping wings at low Reynolds number using numerical simulations

This paper presents a numerical investigation of the effects of chordwise flexibility on flapping wings at low Reynolds number. The numerical simulations are performed with a partitioned fluid–structure interaction algorithm using artificial compressibility stabilization. The choice of the structural dimensionless parameters is based on scaling arguments and is compared against parameters used by other authors. The different regimes, namely inertia-driven and pressure-driven wing deformations, are presented along with their effects on the topology of the flow and on the performance of a heaving and pitching flapping wing in propulsion regime. It is found that pressure-driven deformations can significantly increase the thrust efficiency if a suitable amount of flexibility is used. Significant thrust increases are also observed in zero pitching amplitude cases. The effects of the second and third deformation modes on the performances of pressure-driven deformation cases are discussed. On the other hand, inertia-driven deformations generally deteriorate aerodynamic performances of flapping wings unless the behavior of the wing deformation is modified by the presence of sustainable superharmonics in a way that produces slight improvements. It is also shown that wing flexibility can act as an efficient passive pitching mechanism that allows fair thrust and better efficiency to be achieved when compared to a rigid pitching–heaving wing.     

Reynolds number effects on leading edge vortex development on a waving wing

The waving wing experiment is a fully three-dimensional simplification of the flapping wing motion observed in nature. The spanwise velocity gradient and wing starting and stopping acceleration that exist on an insect-like flapping wing are generated by rotational motion of a finite span wing. The flow development around a waving wing at Reynolds number between 10,000 and 60,000 has been studied using flow visualization and high-speed PIV to capture the unsteady velocity field. Lift and drag forces have been measured over a range of angles of attack, and the lift curve shape was similar in all cases. A transient high-lift peak approximately 1.5 times the quasi-steady value occurred in the first chord length of travel, caused by the formation of a strong attached leading edge vortex. This vortex appears to develop and shed more quickly at lower Reynolds numbers. The circulation of the leading edge vortex has been measured and agrees well with force data.     

Effects of mass and chordwise flexibility on 2D self-propelled flapping wings

A self-propelled flexible flapping wing 2D numerical model undergoing a combined pitching and heaving motion is presented. Since such freely moving foil experiences zero net thrust, a definition of efficiency for this kind of problem is proposed and discussed against other formulations found in the literature. It is also shown that the deviation motion of wings such as that found in natural flyers is likely a consequence of the fluid–structure dynamics of the wings. The passive deviation motion observed in numerical simulations is either a consequence of a feathering mechanism referred to as rigid feathering or of the inertial displacement caused by the wing deformation. The effects of flexibility on the performance of the wing are also presented. It is found that flexibility may significantly enhance the efficiency in pressure-driven deformation cases. The rigid feathering mechanism is found to have an effect similar to that of the feathering caused by wing flexibility on the performances of pressure-driven deformation cases.     

Effect of the domain spanwise periodic length on the flow around a circular cylinder

We assess the effect of the choice of spanwise periodic length on simulations of the flow around a fixed circular cylinder. The Reynolds number is set to 400 because, at this value, both lift coefficient and shedding frequency show significant drop due to three-dimensional flow structures. From the analysis of the three-dimensionalities of the wake and of the integral quantities such as Strouhal number, RMS of lift coefficient and energy contained in the three-dimensional portion of the flow we obtain an estimate of the minimum spanwise length to satisfactorily represent the flow. Furthermore, we observe a distinct wake behavior when the spanwise length is approximately the mode B instability wavelength.     

扑翼飞行器带翼刀机翼气动特性实验研究

通过进行微型扑翼飞行器低速风洞试验,研究了带弯度机翼下翼面翼刀对扑翼飞行器升阻特性的影响。文中进行了带翼刀机翼和不带翼刀机翼在不同迎角下的风洞吹风试验。试验结果表明,带翼刀机翼升力系数大于不带翼刀机翼升力系数,从而证明了翼刀可以阻止机翼下表面气流展向流动,起到增加机翼升力的作用。当扑翼在小迎角飞行时,带翼刀机翼可以有效地提高扑翼的气动效率,改善扑翼的飞行性能。研究结果可为带翼刀机翼在扑翼飞行器上的应用提供技术支持。