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.     

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%.     

Reynolds number, thickness and camber effects on flapping airfoil propulsion

The effect of varying airfoil thickness and camber on plunging and combined pitching and plunging airfoil propulsion at Reynolds number Re=200, 2000, 20 000 and 2×10⁶ was studied by numerical simulations for fully laminar and fully turbulent flow regimes. The thickness study was performed on 2-D NACA symmetric airfoils with 6-50% thick sections undergoing pure plunging motion at reduced frequency k=2 and amplitudes h=0.25 and 0.5, and for combined pitching and plunging motion at k=2, h=0.5, phase ?=90°, pitch angle θ_o=15° and 30° and the pitch axis was located at 1/3 of chord from leading edge. At Re=200 for motions where positive thrust is generated, thin airfoils outperform thick airfoils. At higher Re significant gains could be achieved both in thrust generation and propulsive efficiency by using a thicker airfoil section for plunging and combined motion with low pitch amplitude. The camber study was performed on 2-D NACA airfoils with varying camber locations undergoing pure plunging motion at k=2, h=0.5 and Re=20 000. Little variation in thrust performance was found with camber. The underlying physics behind the alteration in propulsive performance between low and high Reynolds numbers has been explored by comparing viscous Navier-Stokes and inviscid panel method results. The role of leading edge vortices was found to be key to the observed performance variation.     

Thrust augmentation of flapping airfoils in low Reynolds number flow using a flexible membrane

The unsteady aerodynamic thrust and aeroelastic response of a two-dimensional membrane airfoil under prescribed harmonic motion are investigated computationally with a high-order Navier–Stokes solver coupled to a nonlinear membrane structural model. The effects of membrane prestress and elasticity are examined parametrically for selected plunge and pitch–plunge motions at a chord-based Reynolds number of 2500. The importance of inertial membrane loads resulting from the prescribed flapping is also assessed for pure plunging motions. This study compares the period-averaged aerodynamic loads of flexible versus rigid membrane airfoils and highlights the vortex structures and salient fluid–membrane interactions that enable more efficient flapping thrust production in low Reynolds number flows.     

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.     

An experimental study of the unsteady vortex structures in the wake of a root-fixed flapping wing

An experimental study was conducted to characterize the evolution of the unsteady vortex structures in the wake of a root-fixed flapping wing with the wing size, stroke amplitude, and flapping frequency within the range of insect characteristics for the development of novel insect-sized nano-air-vehicles (NAVs). The experiments were conducted in a low-speed wing tunnel with a miniaturized piezoelectric wing (i.e., chord length, C = 12.7 mm) flapping at a frequency of 60 Hz (i.e., f = 60 Hz). The non-dimensional parameters of the flapping wing are chord Reynolds number of Re = 1,200, reduced frequency of k = 3.5, and non-dimensional flapping amplitude at wingtip h = A/C = 1.35. The corresponding Strouhal number (Str) is 0.33, which is well within the optimal range of 0.2 < Str < 0.4 used by flying insects and birds and swimming fishes for locomotion. A digital particle image velocimetry (PIV) system was used to achieve phased-locked and time-averaged flow field measurements to quantify the transient behavior of the wake vortices in relation to the positions of the flapping wing during the upstroke and down stroke flapping cycles. The characteristics of the wake vortex structures in the chordwise cross planes at different wingspan locations were compared quantitatively to elucidate underlying physics for a better understanding of the unsteady aerodynamics of flapping flight and to explore/optimize design paradigms for the development of novel insect-sized, flapping-wing-based NAVs.     

Influence of heat transfer on the aerodynamic performance of a plunging and pitching NACA0012 airfoil at low Reynolds numbers

The unsteady low Reynolds number aerodynamics phenomena around flapping wings are addressed in several investigations. Elsewhere, airfoils at higher Mach numbers and Reynolds numbers have been treated quite comprehensively in the literature. It is duly noted that the influence of heat transfer phenomena on the aerodynamic performance of flapping wings configurations is not well studied. The objective of the present study is to investigate the effect of heat transfer upon the aerodynamic performance of a pitching and plunging NACA0012 airfoil in the low Reynolds number flow regime with particular emphasis upon the airfoil's lift and drag coefficients. The compressible Navier–Stokes equations are solved using a finite volume method. To consider the variation of fluid properties with temperature, the values of dynamic viscosity and thermal diffusivity are evaluated with Sutherland's formula and the Eucken model, respectively. Instantaneous and mean lift and drag coefficients are calculated for several temperature differences between the airfoil surface and freestream within the range 0–100 K. Simulations are performed for a prescribed airfoil motion schedule and flow parameters. It is learnt that the aerodynamic performance in terms of the lift C_L and drag C_D behavior is strongly dependent upon the heat transfer rate from the airfoil to the flow field. In the plunging case, the mean value of C_D tends to increase, whereas the amplitude of C_L tends to decrease with increasing temperature difference. In the pitching case, on the other hand, the mean value and the amplitude of both C_D and C_L decrease. A spectral analysis of C_D and C_L in the pitching case shows that the amplitudes of both C_D and C_L decrease with increasing surface temperature, whereas the harmonic frequencies are not affected.     

An investigation into the effects of unsteady parameters on the aerodynamics of a low Reynolds number pitching airfoil

The growing applications of low Reynolds number (LRN) operating vehicles impose the need for accurate LRN flow solutions. These applications usually involve complex unsteady phenomena, which depend on the kinematics of the vehicle such as pitching, plunging, and flapping of a wing. The objective of the present study is to address the issues related to LRN aerodynamics of a harmonically pitching NACA0012 airfoil. To this end, the influence of unsteady parameters, namely, amplitude of oscillation, d, reduced frequency, k, and Reynolds number, Re, on the aerodynamic performance of the model is investigated. Computational fluid dynamics (CFD) is utilized to solve Navier–Stokes (N–S) equations discretized based on the Finite Volume Method (FVM). The resulting instantaneous lift coefficients are compared with analytical data from Theodorsen’s method. The simulation results reveal that d, k, and Re are of great importance in the aerodynamic performance of the system, as they affect the maximum lift coefficients, hysteresis loops, strength, and number of the generated vortices within the harmonic motion, and the extent of the so-called figure-of-eight phenomenon region. Thus, achieving the optimum lift coefficients demands a careful selection of these parameters.     

Investigation of steady plasma actuation effect on aerodynamic coefficients of oscillating airfoil at low Reynolds number

In this work, numerical study of two dimensional laminar incompressible flow around an oscillating NACA0012 airfoil is proceeded using the open source code Open FOAM. Oscillatory motion types including pitching and flapping are considered. Reynolds number for these motions is assumed to be 12000 and effects of these motions and also different unsteady parameters such as amplitude and reduced frequency on aerodynamic coefficients are studied. For flow control on airfoil, dielectric barrier discharge plasma actuator is used in two different positions on airfoil and its effect is compared for the two types of considered oscillating motions. It is observed that in pitching motion, imposing plasma leads to an improvement in aerodynamic coefficients, but it does not have any positive effect on flapping motion.Also, for the amplitudes and frequencies investigated in this paper, the trailing edge plasma had a more desirable effect than other positions.     

Aerodynamics of laminar separation flutter at a transitional Reynolds number

Experimental observations of self-sustained pitch oscillations of a NACA 0012 airfoil at transitional Reynolds numbers were recently reported. The aeroelastic limit cycle oscillations, herein labelled as laminar separation flutter, occur in the range 5.0×10⁴≤Re_c≤1.3×10⁵. They are well behaved, have a small amplitude and oscillate about θ=0°. It has been speculated that laminar separation leading to the formation of a laminar separation bubble, occurring at these Reynolds numbers, plays an essential role in these oscillations. This paper focuses on the Re_c=7.7×10⁴ case, with the elastic axis located at 18.6% chord. Considering that the experimental rig acts as a dynamic balance, the aerodynamic moment is derived and is empirically modelled as a generalized Duffing–van-der-Pol nonlinearity. As expected, it behaves nonlinearly with pitch displacement and rate. It also indicates a dynamically unstable equilibrium point, i.e. negative aerodynamic damping. In addition, large eddy simulations of the flow around the airfoil undergoing prescribed simple harmonic motion, using the same amplitude and frequency as the aeroelastic oscillations, are performed. The comparison between the experiment and simulations is conclusive. Both approaches show that the work done by the airflow on the airfoil is positive and both have the same magnitude. The large eddy simulation (LES) computations indicate that at θ=0°, the pitching motion induces a lag in the separation point on both surfaces of the airfoil resulting in negative pitching moment when pitching down, and positive moment when pitching up, thus feeding the LCO.     

Camber effects in the dynamic aeroelasticity of compliant airfoils

This paper numerically investigates the effect of chordwise flexibility on the dynamic stability of compliant airfoils. A classical two-dimensional aeroelastic model is expanded with an additional degree of freedom to capture time-varying camber deformations, defined by a parabolic bending profile of the mean aerodynamic chord. Aerodynamic forces are obtained from unsteady thin airfoil theory and the corresponding compliant-airfoil inertia and stiffness from finite-element analysis. V–g and state-space stability methods have been implemented in order to compute flutter speeds. The study looks at physical realizations with an increasing number of degrees of freedom, starting with a camber-alone system. It is shown that single camber leads to flutter, which occurs at a constant reduced frequency and is due to the lock in between the shed wake and the camber motion. The different combinations of camber deformations with pitch and plunge motions are also studied, including parametric analyses of their aeroelastic stability characteristics. A number of situations are identified in which the flutter boundary of the compliant airfoil exhibits a significant dip with respect to the rigid airfoil models. These results can be used as a first estimation of the aeroelastic stability boundaries of membrane-wing micro air vehicles.     

Simulation of incompressible viscous flows around moving objects by a variant of immersed boundary‐lattice Boltzmann method

A variant of immersed boundary‐lattice Boltzmann method (IB‐LBM) is presented in this paper to simulate incompressible viscous flows around moving objects. As compared with the conventional IB‐LBM where the force density is computed explicitly by Hook's law or the direct forcing method and the non‐slip condition is only approximately satisfied, in the present work, the force density term is considered as the velocity correction which is determined by enforcing the non‐slip condition at the boundary. The lift and drag forces on the moving object can be easily calculated via the velocity correction on the boundary points. The capability of the present method for moving objects is well demonstrated through its application to simulate flows around a moving circular cylinder, a rotationally oscillating cylinder, and an elliptic flapping wing. Furthermore, the simulation of flows around a flapping flexible airfoil is carried out to exhibit the ability of the present method for implementing the elastic boundary condition. It was found that under certain conditions, the flapping flexible airfoil can generate larger propulsive force than the flapping rigid airfoil. Copyright © 2009 John Wiley & Sons, Ltd.     

Correlation between vortex structures and unsteady loads for flapping motion in hover

During the past decade, efforts were made to develop a new generation of unmanned aircrafts, qualified as Micro-Air Vehicles. The particularity of these systems resides in their maximum dimension limited to 15 cm, which, in terms of aerodynamics, corresponds to low Reynolds number flows (Re ≈ 10² to 10⁴). At low Reynolds number, the concept of flapping wings seems to be an interesting alternative to the conventional fixed and rotary wings. Despite the fact that this concept may lead to enhanced lift forces and efficiency ratios, it allows hovering coupled with a low-noise generation. Previous studies (Dickinson et al. in Science 284:1954–1960, 1999) revealed that the flow engendered by flapping wings is highly vortical and unsteady, inducing significant temporal variations of the loads experienced by the airfoil. In order to enhance the aerodynamic performance of such flapping wings, it is essential to give further insight into the loads generating mechanisms by correlating the spatial and temporal evolution of the vortical structures together with the time-dependent lift and drag. In this paper, Time Resolved Particle Image Velocimetry is used as a basis to evaluate both unsteady forces and vortical structures generated by an airfoil undergoing complex motion (i.e. asymmetric flapping flight), through the momentum equation approach and a multidimensional wavelet-like vortex parameterization method, respectively. The momentum equation approach relies on the integration of flow variables inside and around a control volume surrounding the airfoil (Noca et al. in J Fluids Struct 11:345–350, 1997; Unal et al. in J Fluids Struct 11:965–971, 1997). Besides the direct link performed between the flow behavior and the force mechanisms, the load characterization is here non-intrusive and specifically convenient for flapping flight studies thanks to its low Reynolds flows’ sensitivity and adaptability to moving bodies. Results are supported by a vortex parameterization which evaluates the circulation of the multiple vortices generated in such complex flows. The temporal evolution of the loads matches the flow behavior and hence reveals the preponderant inertial force component and that due to vortical structures.     

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.     

Parameter dependence of vortex interactions on a two-dimensional plunging plate

The structure and dynamics of the flow field created by a plunging flat-plate airfoil are investigated at a chord Reynolds number of 10,000 while varying plunge amplitude and Strouhal number. Digital particle image velocimetry measurements are used to characterize the shedding patterns and the interactions between the leading- and trailing-edge vortex structures (LEV and TEV), resulting in the development of a wake classification system based on the nature and timing of interactions between the leading- and trailing-edge vortices. The streamwise advancement of the LEV during a plunge cycle and its resulting interaction with the TEV is primarily dependent on reduced frequency; however, for Strouhal numbers above approximately 0.4, significant changes are observed in the formation of vortices shed from the leading and trailing edges, as well as the circulation of the leading-edge vortex. The functional form of the relationship between leading-edge vortex circulation and Strouhal number suggests that the Strouhal number dependence is more specifically a manifestation of the effective angle of attack. Comparison with low-Reynolds-number studies of plunging airfoil aerodynamics reveals a high degree of consistency and suggests applicability of the classification system beyond the range examined in the present work.     

Numerical study of large amplitude,nonsinusoidal motion and camber effects on pitching airfoil propulsion

The effects of large amplitude and nonsinusoidal motion on pitching airfoil aerodynamics for thrust generation were numerically studied with a 2-D NACA0012 airfoil used, and various 2-D NACA asymmetric airfoils were applied for camber effect study. The large amplitude effect study has been undertaken over a wide range of reduced frequency k (from 6 to 18) and pitching amplitude θ (from 5° to 30°) at Re=1.35×10⁴ with sinusoidal pitching profile used. It is shown that the large pitching amplitude results in much more thrust generated than that at low pitching amplitude and the increase of thrust with amplitude becomes slow when the amplitude reaches some degree. However, the propulsive efficiency noticeably decreases with the increase of θ at a fixed k.An adjustable parameter K was employed to realize various nonsinusoidal motions and the effect of nonsinusoidal motion was investigated with various unsteady parameters (θ, k) applied. The results reveal that nonsinusoidal motion has a noticeable effect on the aerodynamic performance, as it affects the instantaneous force coefficients, maximum thrust coefficients and flow structures. An increase in K results in a better thrust generation performance at fixed θ and k, especially for K>0. It is also shown that the larger K noticeably influences the wake pattern and induces a stronger reverse von Karman vortex street in the wake, which in turn leads to the increased thrust. The camber study was performed on various 2-D NACA airfoils with different cambers and camber locations undergoing sinusoidal pitching motion at θ=5° and Re=1.35×10⁴. It is found that varying camber offers little improvement in thrust generation performance.     

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.

     

Schooling behavior of rigid and flexible heaving airfoils

The schooling behavior of rigid and flexible NACA0017 airfoils undergoing a heaving motion was experimentally explored using a merry-go-round configuration. Each airfoil was attached to the end of a horizontal support bar whose other end was connected to a freely rotating vertical axis. The axis was forced to undergo a sinusoidal motion in the vertical direction to generate a pure heaving motion of the airfoil in the frequency range of 0.4 to 4.8 Hz. The propulsion due to the heaving airfoil was expressed as the horizontal rotational speed of the support bar. This experimental setup simulates an infinite schooling of airfoils separated by a streamwise distance d undergoing in-phase heaving motions. The ratio of the distance to the chord length, d/c, was determined by the number of airfoils (2 ≤ n ≤ 8). The variation in rotational frequency F as a function of heaving frequency f was determined using different experimental parameters. The schooling number S = f /(nF), which represents the number of heaving oscillations between each pair of successive airfoils, was introduced to explain the schooling behavior of the airfoils. The effects of airfoil flexibility, d/c and f on the propulsive performance were examined in the context of the schooling behavior of the airfoils.     

Power extraction efficiency improvement of a fully-activated flapping foil: With the help of an auxiliary rotating foil

The power extraction efficiency improvement of a fully-activated flapping foil with the help of an auxiliary rotating foil is numerically examined in this work. A NACA0015 airfoil is placed in a two-dimensional laminar flow and synchronously performs the imposed pitching and plunging motions. An auxiliary smaller foil, which rotates about its center, is arranged below the flapping foil. As a consequence, the vortex interaction between the flapping foil and the rotating foil occurs. At a Reynolds number of 1100 and the position of the pitching axis at one-third chord, the effects of the distance between the flapping foil and the auxiliary foil, the phase difference between the rotating motion and the flapping motion as well as the frequency of flapping motion on the power extraction performance are systematically investigated. It is found that compared to the single flapping foil, the efficiency of power extraction for the flapping foil with an auxiliary device can be improved. Based on the numerical analysis, it is indicated that the enhanced plunging contribution, which is caused by the increased lift force owing to the vortex interaction, directly helps the efficiency improvement.     

How a flexible tail improves the power extraction efficiency of a semi-activated flapping foil system: A numerical study

The improvement of power extraction of a semi-activated flapping foil system via the use of a flexible tail is numerically investigated in this work. A NACA0015 airfoil arranged in a two-dimensional laminar flow synchronously executes a forced pitching motion and an induced plunging motion. A flat plate attached to the trailing edge of the foil is utilized to model a tail, and thereby they are considered as a unit for the purpose of power extraction. The tail is either rigid or deformable. At a Reynolds number of 1100 and the position of the pitching axis at third chord, the effects of the mass and flexibility of the tail as well as the frequency of pitching motion on the net power extraction are systematically examined. It is found that compared to the foil with a rigid tail, the efficiency of net power extraction for the foil with a deformable tail can be improved. Based on the numerical analysis, it is indicated that the enhanced power extraction, which is caused by the increased lift force, directly contributes to the net efficiency improvement. Moreover, owing to high enhancement of power extraction, a flexible tail with high flexibility is recommended in the semi-activated flapping foil based power extraction system.