Numerical investigation of transonic flow over deformable airfoil with plunging motion

In this article, the transonic inviscid flow over a deformable airfoil with plunging motion is studied numerically. A finite volume method based on the Roe scheme developed in a generalized coordinate is used along with an arbitrary Lagrangian-Eulerian method and a dynamic mesh algorithm to track the instantaneous position of the airfoil.The effects of different governing parameters such as the phase angle, the deformation amplitude, the initial angle of attack, the flapping frequency, and the Mach number on the unsteady flow field and aerodynamic coefficients are investigated in detail. The results show that maneuverability of the airfoil under various flow conditions is improved by the deformation. In addition, as the oscillation frequency of the airfoil increases, its aerodynamic performance is significantly improved.     

Flow field characteristics study of a flapping airfoil using computational fluid dynamics

The flow field of a flapping airfoil in Low Reynolds Number (LRN) flow regime is associated with complex nonlinear vortex shedding and viscous phenomena. The respective fluid dynamics of such a flow is investigated here through Computational Fluid Dynamics (CFD) based on the Finite Volume Method (FVM). The governing equations are the unsteady, incompressible two-dimensional Navier-Stokes (N-S) equations. The airfoil is a thin ellipsoidal geometry performing a modified figure-of-eight-like flapping pattern. The flow field and vortical patterns around the airfoil are examined in detail, and the effects of several unsteady flow and system parameters on the flow characteristics are explored. The investigated parameters are the amplitude of pitching oscillations, phase angle between pitching and plunging motions, mean angle of attack, Reynolds number (Re), Strouhal number (St) based on the translational amplitudes of oscillations, and the pitching axis location (x/c). It is shown that these parameters change the instantaneous force coefficients quantitatively and qualitatively. It is also observed that the strength, interaction, and convection of the vortical structures surrounding the airfoil are significantly affected by the variations of these parameters.     

Design and development considerations for biologically inspired flapping-wing micro air vehicles

In this paper, the decade of numerical and experimental investigations leading to the development of the authors’ unique flapping-wing micro air vehicle is summarized. Early investigations included the study of boundary layer energization by means of a small flapping foil embedded in a flat-plate boundary layer, the reduction of the recirculatory flow region behind a backward-facing step by means of a small flapping foil, and the reduction or suppression of flow separation behind blunt or cusped airfoil trailing edges by flapping a small foil located in the wake flow region. These studies were followed by systematic investigations of the aerodynamic characteristics of single flapping airfoils and airfoil combinations. These unsteady flows were described using flow visualization, laser-Doppler velocimetry in addition to panel and Navier–Stokes computations. It is then shown how this flapping-wing database was used to conceive, design and develop a micro air vehicle which has a fixed wing for lift and two flapping wings for thrust generation. While animal flight is characterized by a coupled force generation, the present design appears to separate lift and thrust. However, in fact, the performance of one surface is closely coupled to the other surfaces.     

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.     

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.     

微型飞行器低雷诺数空气动力学

微型飞行器(MAVs)设计绝不是常规飞行器在尺度上的简单缩小,面临许多技术难题.其中微型飞行器低雷诺数空气动力学是其最为根本的技术瓶颈之一,也是当前受到广泛关注的热点之一.本文紧密结合微型飞行器技术,对这一领域中所面临的低雷诺数空气动力学问题和近两年来该方向国内一些新的进展进行了较为详细的介绍.按照MAVs飞行方式和结构特性进行分类,简单介绍微型飞行器研究中的低$Re$数空气动力学问题.首先介绍了二维和三维固定翼低雷诺数空气动力学问题:包括层流分离泡,翼型升力系数小攻角非线性效应,静态迟滞效应,以及低$Re$数小展弦比机翼气动特性.第2,介绍了拍动翼低雷诺数空气动力学方面的研究工作.包括前人提出的昆虫低$Re$数下获得高升力的多种非定常拍动翼飞行机制:Wagner效应、Weis-Fogh效应(clap-and-fling)、延迟失速效应(delayedstall)、Kramer效应(rotational forces)、尾迹捕获效应(wakecapture)、附加质量效应(addedmass)等.以及国内学者近几年在拍动翼方面取得的一些研究成果.第3,介绍了柔性翼低雷诺数气动问题.研究表明柔性翼对于固定翼微型飞行器提高抗阵风能力,拍动翼微型飞行器产生足够的升力和推力.最后简单介绍了可变形翼(morphingwing)微型飞行器方面的一些研究工作,指出微型飞行器技术可以通过采用可变形翼设计,突破众多的技术瓶颈.另一方面,可变形翼概念可以通过在低成本,低速的MAVs上进行飞行试验,获得非常好的验证平台.      

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.     

低Reynolds数NACA0012翼型绕流的流动特性分析

在水槽中应用PIV测速技术研究了NACA0012翼型在Reynolds数为8200时的流动特性,重点关注了翼型绕流结构中主频和扰动增长速率随迎角的变化。结果表明,分离剪切层的扰动增长符合指数规律;且随着迎角的增大,转捩过程加速,表现为扰动增长率逐渐增大,转捩的起始位置逐渐向上游移动。在所有实验迎角情况下,流场均由脱落旋涡主导,但其主导作用随着迎角的增大而削弱。     

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.     

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.     

High fidelity numerical simulation of airfoil thickness and kinematics effects on flapping airfoil propulsion

High-fidelity numerical simulations with the spectral difference (SD) method are carried out to investigate the unsteady flow over a series of oscillating NACA 4-digit airfoils. Airfoil thickness and kinematics effects on the flapping airfoil propulsion are highlighted. It is confirmed that the aerodynamic performance of airfoils with different thickness can be very different under the same kinematics. Distinct evolutionary patterns of vortical structures are analyzed to unveil the underlying flow physics behind the diverse flow phenomena associated with different airfoil thickness and kinematics and reveal the synthetic effects of airfoil thickness and kinematics on the propulsive performance. Thickness effects at various reduced frequencies and Strouhal numbers for the same chord length based Reynolds number (=1200) are then discussed in detail. It is found that at relatively small Strouhal number (=0.3), for all types of airfoils with the combined pitching and plunging motion (pitch angle 20°, the pitch axis located at one third of chord length from the leading edge, pitch leading plunge by 75°), low reduced frequency (=1) is conducive for both the thrust production and propulsive efficiency. Moreover, relatively thin airfoils (e.g. NACA0006) can generate larger thrust and maintain higher propulsive efficiency than thick airfoils (e.g. NACA0030). However, with the same kinematics but at relatively large Strouhal number (=0.45), it is found that airfoils with different thickness exhibit diverse trend on thrust production and propulsive efficiency, especially at large reduced frequency (=3.5). Results on effects of airfoil thickness based Reynolds numbers indicate that relative thin airfoils show superior propulsion performance in the tested Reynolds number range. The evolution of leading edge vortices and the interaction between the leading and trailing edge vortices play key roles in flapping airfoil propulsive performance.     

Aerodynamic mechanism of forces generated by twisting model-wing in bat flapping flight

The aerodynamic mechanism of the bat wing membrane Mong the lateral border of its body is studied. The twist-morphing that alters the angle of attack （AOA） along the span-wise direction is observed widely during bat flapping flight. An assumption is made that the linearly distributed AOA is along the span-wise direction. The plate with the aspect ratio of 3 is used to model a bat wing. A three-dimensional （3D） unsteady panel method is used to predict the aerodynamic forces generated by the flapping plate with leading edge separation. It is found that, relative to the rigid wing flapping, twisting motion can increase the averaged lift by as much as 25% and produce thrust instead of drag. Furthermore, the aerodynamic forces （lift/drag） generated by a twisting plate-wing are similar to those of a pitching rigid-wing, meaning that the twisting in bat flight has the same function as the supination/pronation motion in insect flight.     

Dynamics of co-rotating vortices in a flow around a bio-inspired corrugated airfoil

Bio-inspired corrugated airfoils show favourable aerodynamic characteristics such as high coefficient of lift and delayed stall at low Reynolds numbers. Two-dimensional (2D) direct numerical simulation has been performed here on a corrugated airfoil at various angles of attack (0°, +5°, -5°) and Reynolds number of 280 to 6700. The objective is to analyse the pressure variation inside the corrugations and correlate it to the vortex movement across the corrugations and the overall aerodynamic characteristics of the corrugated airfoil. The flow characteristics have been examined based on the local Strouhal numbers in the corrugations of the airfoil. It is observed that the pressure variation in each corrugation is the result of vortex merging and separation in the corrugation which plays a major role in changing the flow characteristics. The Strouhal number of the flow is dictated by the most dominant local Strouhal number. The numerical results are further compared with experimental results obtained using particle image velocimetry, and the two set of results are found to match well. These results are significant because they elucidate the effect of corrugation, angle of attack, and Reynolds number on flow over a corrugated airfoil.     

Bond graph modeling of a typical flapping wing micro-air-vehicle with the elastic articulated wings

In the present research, a new comprehensive model of a flexible articulated flapping wing robot using the bond graph approach is presented. The flapping kinematics of a two-section wing is introduced via the bond graph based approach on a hybrid mechanism providing amplitude and phase characteristics. The aerodynamic quasi-steady approach equipped with stall correlation is utilized according to the reduced flapping frequency and the angle of attack ranges. The local flow velocity and the wing position are calculated in both wing and body coordinates taking into account rotation and translation of the wing different parts. Estimation of the effective angle of attack is performed by calculating the instantaneous torque distribution on both wing sections. Aeroelastic modeling is employed in which the wing structure is assumed as an elastic Euler–Bernoulli beam at the leading edge with three linear torsional modes. In this novel integrated bond graph model, computation of the performance indices including the average lift and thrust, consumed and produced powers by flapping and mechanical efficiency are presented. Due to existence of the numerous geometric and kinematic parameters in articulated flexible flapping wing, such a model is essential for design and optimization. Consequently, an example of a typical parametric study and the results validation are carried out. It is indicated that the sensitivity of the bird performance to relative change in design variables would increase for out of phase flapping, second part stiffness, flapping amplitude, frequency and velocity respectively. It is interesting that by employing the reverse-phase flapping which is possible only via articulated wings, the maximum efficiency could be achieved. In addition, it is shown that adjusting the wing torsional stiffness is a crucial item in design of passive flapping robots. The key advantage of the two-section flapping wing is depicted as the controlling capability of the angle of attack in the outer part of the wing. Finally, the improved version of the bird is being addressed by approximately 15% progress in propulsive efficiency.

     

Design,manufacturing, and flight testing of a fixed wing micro air vehicle with Zimmerman planform

An optimized and comprehensive method is used to design and manufacture a fixed wing micro air vehicle (MAV) with Zimmerman planform. The design process includes four stages which are the specification of the flight mission, determination of the best aspect ratio, identification of the optimum wing loading and thrust loading values, and estimation of the weight of the structural components of the MAV. To do this, various statistical and analytical methods are utilized. Based on an aerodynamic analysis, the results show that an optimum aspect ratio that maximizes the performance of the Zimmerman MAV for a well-defined cruise speed is determined. Considering six possible flights, a constraint analysis is performed and an optimum wing loading value is determined. It is shown that the computational method is beneficial to determine the exact masses for the structural components including the wing, fuselage, and vertical tail. Using the 3D panel method, the determination of the shape of the reflexed airfoil for the MAV is successfully done by minimizing the drag force and the angle of attack to use less powerful motor and avoid any stall effect, respectively. A stability analysis is then performed to check the safe flight of the designed vehicle. During test flight, the results show that the designed Zimmerman MAV satisfies the pre-defined specification. The final characteristics of the manufactured MAV are: wingspan of 44 cm, weight of 450 g, aspect ratio of 1.51, cruise speed of 20 m/s, and flight endurance of 20 min.     

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.     

Investigation of the vortex induced unsteadiness of a separation bubble via time-resolved and scanning PIV measurements

A transitional separation bubble on the suction side of an SD7003 airfoil is considered. The transition process that forces the separated shear layer to reattach seems to be governed by Kelvin–Helmholtz instabilities. Large scale vortices are formed due to this mechanism at the downstream end of the bubble. These vortices possess a three-dimensional structure and detach from the recirculation region, while other vortices are formed within the bubble. This separation of the vortex is a highly unsteady process, which leads to a bubble flapping. The structure of these vortices and the flapping of the separation bubble due to these vortices are temporally and spatially analyzed at angles of attack from 4° to 8° and chord-length based Reynolds numbers Re _c = 20,000–60,000 using time-resolved PIV measurements in a 2D and a 3D set-up, i.e., stereo-scanning PIV measurements are done in the latter case. These measurements complete former studies at a Reynolds number of Re _c = 20,000. The results of the time-resolved PIV measurements in a single light-sheet show the influence of the angle of attack and the Reynolds number. The characteristic parameters of the separation bubble are analyzed focusing on the unsteadiness of the separation bubble, e.g., the varying size of the main recirculation region, which characterizes the bubble flapping, and the corresponding Strouhal number are investigated. Furthermore, the impact of the freestream turbulence is investigated by juxtaposing the current and former results. The stereo-scanning PIV measurements at Reynolds numbers up to 60,000 elucidate the three-dimensional character of the vortical structures, which evolve at the downstream end of the separation bubble. It is shown that the same typical structures are formed, e.g., the c-shape vortex and the screwdriver vortex at each Reynolds number and angle of attack investigated and the occurrence of these patterns in relation to Λ-structures is discussed. To evidence the impact of the freestream turbulence, these results are compared with findings of former measurements.     

Efficiency of an auto-propelled flapping airfoil

The present study deals with an investigation of the flow aerodynamic characteristics and the propulsive velocity of a system equipped with a nature inspired propulsion system. In particular, the study is aimed at studying the effect of the flapping frequency on the flow behavior. We consider a NACA0014 airfoil undergoing a vertical sinusoidal flapping motion. In contrast to nearly all previous studies in the literature, the present work does not impose any velocity on the inlet flow. During each iteration the outer flow velocity is computed after having determined the forces exerted on the airfoil. Forward motion may only be produced by flapping motion of the airfoil. This is more consistent with the physical phenomenon. The non-stationary viscous flow around the flapping airfoil is simulated using Ansys-Fluent 12.0.7. The airfoil movement is achieved using the deformable mesh technique and an in-house developed User Define Function (UDF). Our results show the influence of flapping frequency and amplitude on both the airfoil velocity and the propulsive efficiency. The resulting motion is contrasts to the applied forces. In the present study, the frequency ranges from 0.1 to 20 Hz while the airfoil amplitude values considered are: 10%, 17.5%, 25% and 40%.     

Effect of one-sided unsteady energy supply on aerodynamic characteristics of airfoils in a transonic flow

The possibility of controlling the aerodynamic characteristics of airfoils in transonic flight regimes by means of one-sided pulsed-periodic energy supply is studied. Based on the numerical solution of two-dimensional unsteady gas-dynamic equations, the change in the flow structure in the vicinity of a symmetric airfoil at different angles of attack and the aerodynamic characteristics of the airfoil as functions of the amount of energy supplied asymmetrically (with respect to the airfoil) are determined. The results obtained are compared with the data calculated for the flow past the airfoil at different angles of attack without energy supply. It is found that a given lift force can be obtained with the use of energy supply at a much better lift-to-drag ratio of the airfoil, as compared to the case of the flow past the airfoil at an angle of attack. The moment characteristics of the airfoil are found. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 49, No. 6, pp. 82–87, November–December, 2008.     

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

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