International Journal of Aerospace Engineering
Volume 2019 (2019), Article ID 5151808, 12 pages
https://doi.org/10.1155/2019/5151808

Research Article

The Effect of Damping Coefficient, Spring Coefficient, and Mass Ratio on the Power Extraction Performance of a Semiactive Flapping Wing

Jianyang Zhu1,2

¹Institute of Robotics and Intelligent Systems, Wuhan University of Science and Technology, Wuhan 430081, China
²Hubei Key Laboratory of Mechanical Transmission and Manufacturing Engineering, Wuhan University of Science and Technology, China

Correspondence should be addressed to Jianyang Zhu; zhujianyang02@163.com

Received 15 January 2019; Accepted 27 February 2019; Published 16 April 2019

Academic Editor: Gustaaf B. Jacobs

Copyright © 2019 Jianyang Zhu. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The effect of varying damping coefficient , spring coefficient , and mass ratio on the semiactive flapping wing power extraction performance was numerically studied in this paper. A numerical code based on Finite Volume method to solve the two-dimensional Navier-Stokes equations and coupled with Finite Center Difference method to solve the passive plunging motion equation is developed. At a Reynolds number of 3400 and the pitching axis at quarter chord from the leading edge of the wing, the power extraction performance of the semiactive flapping wing with different damping coefficient, spring coefficient, and mass ratio is systematically investigated. The optimal set of spring coefficient is found at a value of 1.00. However, the variation of mass ratio cannot increase the maximum mean power coefficient and power efficiency, but it can influence the value of damping coefficient at which the wing achieves the maximum mean power coefficient and power efficiency. Moreover, insensitivity of the mean power coefficient and power efficiency to the variation of damping coefficient is observed for the wing with smaller mass ratio, which indicates the wing with smaller has better working stability.

1. Introduction

The application of flapping wing for energy extraction is inspired by insects and birds, who exhibit excellent aerodynamic performance by extracting wind energy through flapping their wings. Comparing to the conventional rotation energy extraction turbines, flapping wing power generator possesses many advantages, such as simpler design, less construction costs, more efficiency in low stream speeds, and more friendly to the flying creatures in nature [1]. Therefore, more and more recently researches have been focused on this type of power generator [2–5].

Up to now, there are three types of flapping wing power generator, named fully active [6], semiactive, and purely passive system, respectively [7]. Among them, the semiactive flapping wing power generator which has one of the flapping motion (pitching) is actuated and another motion (plunging) is induced by free stream fluid and is considered as a more feasible approach in industry, because it is easy to implement and control [8, 9]. This type of flapping wing power generator is characterized by the interaction coupling of fluid, driven wing pitching, and passive wing plunging, which results the flow around the generator very complex, and the energy extraction performance is still the main research focus of this type of power generator.

To investigate the energy extraction performance of a semiactive flapping wing near solid walls, the mechanical parameters of the wing to achieve power extraction efficiency were optimized by Wu et al. [10]. They concluded that the spring constant at a value of 1.0 and damping coefficient at a value of are an optimal choice to achieve net power extraction efficiency of the wing, when the wing has fixed wing mass at a value of 1.0 and Reynolds number at a value of 1100. Under the same Reynolds number, Zhan et al. [11] performed a similar numerical study to optimize the energy extraction performance of a semiactive flapping wing with fixed wing mass at a value of 1.0 and reduced frequency . It was found that the parameters for the wing that achieved best power extraction are the difference with different pitching amplitude. For the wing with pitching amplitude , the optimal power extraction efficiency is achieved when the wing has the spring constant at a value of 5.0 and damping coefficient at a value of , while for the wing with pitching amplitude , the optimal power extraction efficiency is achieved when the wing has the spring constant at a value of 5.0 and damping coefficient at a value of 0.5π. Moreover, based on the numerical study by Zhu et al. [12], the optimal power extraction performance is the system where there is no spring on the plunge motion and the damping coefficient is at a value of .

On the other hand, to explore an inertial effect on the power extraction performance of a semiactive flapping wing, Deng et al. [13] conducted a strong fluid-structure coupling method to study the power extraction performance of a semiactive flapping wing with different mass ratios (ranging from 0.125 to 100), and it was found that the energy harvesting efficiency decreases monotonically with increasing mass ratio. However, the amount of power extraction changes very little when the wing has mass ratio less than 10. While, according to the experiment study on a flapping wing hydroelectric power generation system by Abiru and Yoshitake [14], the wing with lager mass ratio is needed to excite the hydroelastic response; therefore, the wing with lager mass ratio has better power extraction performance.

As discussed above, the effect of damping coefficient, spring constant, and mass ratio on the power extraction performance of the semiactive flapping wing is still not clearly understood. Therefore, in this paper, a numerical code based on Finite Volume method to solve the two-dimensional Navier-Stokes equations and coupled with Finite Center Difference method to solve passive plunging motion equation is employed. The flow around a semiactive flapping wing is simulated with different damping coefficient, spring constant, and mass ratio. The NACA0015 airfoil is employed to represent the two-dimensional section wing, and the power extraction performance as well as the fluid around the wing are studied details in the following.

2. Problem Definition and Methodology

2.1. Problem Definition

In this paper, NACA0015 airfoil with chord length and mass is employed to model the flapping wing. A damper with damping coefficient and spring constant which is attached to the wing is employed to mimic the power extractor. The schematic view of the semiactive flapping wing is shown in Figure 1, where the profile is defined to drive the pitching motion of the wing; the plunging motion is induced by the lift force of the wing. The pitching point is located at the center line of the wing with a distance 0.25d from the leading edge. To simplify the problem, a cosine pitching motion mode is employed, and the center of mass of the wing is also designed to coincide with the pitching point . Then, the governing equation of the pitching motion and passive plunging motion can be defined as

where is the pitching amplitude, is the pitching frequency, is the lift force in direction as shown in Figure 1. The left and right sides of equation (2) is divided by the mass, then, it can be described as

Figure 1: The schematic view of the semiactive flapping wing.

It is seen from equation (1) to equation (3) that if we fix the pitching amplitude and pitching frequency , the kinematic of the semiactive flapping wing is determined by the following three characteristic parameters:

where is the fluid density, is the wing density, , and S is the area of the wing. Moreover, the Reynolds number Re and reduced frequency can be defined as

where is the stream velocity and is the fluid kinematic viscosity.

When the semiactive flapping wing is working, the power extracted from the fluid can be described as

and the power which is needed to input the system to drive the pitching motion can be described as

where is the fluid torque about the pitching axis and is the pitching velocity. Then, the net power is given by

and the power coefficient is described as

The mean energy coefficient can be defined as

where is the pitching period.

Then, the power extraction efficiency of the generator can be defined as

where is the overall vertical extent of the wing motion (the distance in direction between the highest position and the lowest position reached by the wing, either the leading edge or the trailing edge). It is defined as the maximum value of the following equation:

Note that the drag force has not been mentioned in this work, because it does not contribute to the power production.

2.2. Numerical Method

Based on the free stream velocity and the chord length of the wing, the Reynolds number and Mach number are fixed at a value of 3400 and 0.0015, respectively. According to our previous work [5], at Reynolds number 3400, the flow assuming laminar, the calculating results can match experimental results very well for flapping wing. Therefore, the fluid around the wing is assumed as laminar and incompressible. The governing equation can be described as

where is the velocity vector and is the pressure. Fluent ver. 6.3 proprietary solver is employed to solve equation (13) and equation (14). It has been well validated for simulating flapping wing problems [15, 16]. The second-order upwind algorithm is adopted for the space discretization, and first-order implicit algorithm is selected for time discretization. Meanwhile, SIMPLE algorithm is selected for solving the coupling between the pressure and the velocity. The convergence criterion for the solver is satisfied with residual less than 10⁻³.

Finite Center Difference method is employed to solve passive plunging motion governing equation (3), which is embedded in Fluent ver. 6.3 by using a user-defined function (UDF). Dynamic mesh technique is also employed to update the wing’s position (active pitching and passive plunging) at each time step. For more details about the coupling method to solve interaction of fluid and the semiactive flapping wing, one can refer to our previous work [5, 17].

The computational domain of this work has C type shape. According to the works of Zhu and Peng [18], the distance between the outer boundary and the airfoil is at about 8d lengths; flow field and the aerodynamic performance of the semiactive flapping airfoil are not significant by the far field boundary. Therefore, the radius of the semicircle on the left side of the computational domain is fixed 20d, and the size of the rectangle on the right side of the computational domain is set as , . A triangular mesh system is employed where a C type computational domain (shown in Figure 2) containing an inner flapping domain and an outer fixed domain, and the inner domain which moves according to the wing kinematics, and the first grid is located at 0.001c from the wing surface, the growth rate of the grid size is 1.08.

Figure 2: Schematic of the computational grid and boundary conditions.

Figure 3: The power coefficient of a semiactive flapping wing with different meshes and iteration time steps.

The no-slip wall boundary condition is applied at the semiactive flapping wing. The inlet velocity defined as is imposed on the left, up, and low side of the domain. Pressure-outlet boundary condition which has zero static pressure is imposed on the right side of the domain.

2.3. Validation

To validate the independence of the numerical results on the special-discretization and time-discretization scheme, computations with different meshes and iteration time steps are simulated. These checks are carried out on the representative case of , , , , , and . The power coefficient calculated by the coarse mesh (the first grid is located at 0.005c from the wing surface, the growth rate of the grid size is 1.08, total cells), medium mesh (the first grid is located at 0.001c from the wing surface, the growth rate of the grid size is 1.08, total cells), and fine mesh (the first grid is located at 0.0005c from the wing surface, the growth rate of the grid size is 1.08, total cells) with iteration time steps and is shown in Figure 3. It is seen from this figure that the meshes and time steps have slightly influenced the computation results. Therefore, the medium mesh is sufficiently dense and iteration time steps is sufficiently small to grasp the main flow features over the wing, and they are employed for the next simulations.

To validate the coupling method for the study on the semiactive flapping wing, a typical case of a semiactive flapping wing (wing shape is also NACA0015) which was numerical studied by Teng et al. [19] is simulated. The simulation parameters are set the same as the parameters used by Teng et al. [19]: , and , , , and  kg. The resulting lift coefficient of the wing is plotted against time as shown in Figure 4. Comparing to the result by Teng et al. [19], good agreement is achieved.

Figure 4: The lift coefficient of a semiactive flapping wing.

3. Results and Discussions

As conducted above, the main parameters were related to the response of the semiactive flapping wing including the pitching amplitude , reduced frequency , Reynolds number Re, damping coefficient , spring coefficient , and mass ratio . To simplify the problem, the factor is fixed at 45°, is set at 0.32, and Re keeps a value of . Consequently, the key parameters that affect the power extraction performance of the semiactive flapping wing are the damping coefficient , spring coefficient , and mass ratio . The details of these three parameters studied in this paper are summarized in Table 1. Note that the parameter spaces in this work are almost covering all the parameters studied in the literature as discussed in Introduction, but some controversial conclusions are conducted; therefore, it needs to be studied further.

Table 1: Parameters considered in this paper.

Damping coefficient
Spring coefficient
Mass ratio

0.0625~18.0
0.00~3.00
15~240

3.1. Effect of Damping Coefficients and Spring Coefficients

To examine the effect of damping and spring coefficient, the mass ratio is fixed at 60 first. Figure 5 plots the mean power coefficient and power efficiency against spring coefficient with different damping coefficient . Obviously, in this figure that at a given both mean power coefficient and power efficiency increase first then decrease and achieve a maximum at . However, when is larger than 1.0, both mean power coefficient and power efficiency of the wing decrease sharply, which implies the power extraction performance of the wing is deteriorated. In addition, for , the maximum mean power coefficient of the wing monotonically decrease with , while the maximum power efficiency of the wing increase first then decrease and achieve a maximum at . This conclusion is different with the results obtained by Zhu et al. [12], at which they concluded that the best power extraction performance of the wing is achieved, when the spring coefficient is at 0.00 and the damping coefficient is at a value of . These difference can be due to the potential flow theory used by Zhu et al. [12], for which the viscosity effect and the interaction between wing and vortex are neglected.

(a)

(b)

Figure 5: The mean power coefficient and power efficiency against spring coefficient with different damping coefficient ; (a) mean power coefficient; (b) power efficiency.

To explore the mechanism of how the spring coefficients and damping coefficients influence the power extraction performance of the semiactive flapping wing in detail, we first investigate the performance of the wing with fixed , and three typical cases with , 1.0, and 3.0 are studied. The time variation of power coefficient of the above considered cases is shown in Figure 6. It is clear in this figure that the wing with has the largest power coefficient amplitude, then is the wing with , and the wing with has the smallest power coefficient amplitude, which is consistent with the mean power generation of the wing as shown in Figure 5. Meanwhile, there exists phase difference between the three wings. The wing with has the earliest time to achieve maximum power coefficient, then is the wing with , and the wing with has the latest to achieve maximum power coefficient.

Figure 6: The time variation of power coefficient of the considered cases with and , 1.0, and 3.0.

Besides the power coefficient, the power efficiency of the wing is also influenced by the plunging motion . Figure 7 shows the time variation of of the above three considered cases, where the defined pitching motion is also illustrated for comparison. It is found from this figure that the amplitude of the first increase then decrease with , and the phase difference between pitching and is -0.305T, -0.205T, and 0.475T for the wing with , 1.0, and 3.0, respectively. Moreover, nonharmonic time variation of is observed for the wing with , which indicates the complex lift force induced by the wing with .

Figure 7: The time variation of of the considered cases with , 1.0, and 3.0, .

Figure 8 shows the time variation of lift coefficient of the above three considered cases. It is seen from this figure that different with the time variation of , the amplitude of the first decrease then increase with . The wing with larger amplitude of has smaller amplitude of .

Figure 8: The time variation of lift coefficient of the considered cases with , 1.0, and 3.0, .

To explore the physical mechanism of the effect of variation of spring coefficient, Figure 9 shows the vorticity contours during a semipitching cycle (-0.40T). Due to the symmetrical flapping of the wing, semipitching cycle is considered. It is seen from this figure that the spring coefficient can influence the vortex structure of the wing significantly. Comparing with the wing with , there are less separated vortex in the wake for the wing with , and the separated vortex can reattach to the wing surface during flapping; moreover, there are weaker leading edge vortex generating for the wing with , which results the wing to have smaller lift coefficient as shown in Figure 8. Moreover, the separated vortex in the wake indicates more power disappeared, which is the reason why the wing with has larger mean power coefficient and efficiency. On the other hand, for the wing with , obviously, there are more separated vortex in the wake comparing with the wing with ; however, the stronger leading edge vortex are generated by the wing with , which is the reason why the wing with has larger lift coefficient but smaller mean power coefficient and efficiency.

Figure 9: The vorticity contours of the considered cases with , 1.0, and 3.0, .

After investigating the effect of spring coefficient on the power extraction performance of the semiactive flapping wing, the effect of damping coefficient on the power extraction performance of the semiactive flapping wing is studied. To this end, we study the performance of the wing with fixed , and two typical cases with and 4.0 are investigated in details. The time variation of power coefficient of the above considered cases is shown in Figure 10. It is clear in this figure that the wing with has the larger power coefficient than the wing with almost during the whole flapping cycle, which is consistent with the mean power generation of the wing as shown in Figure 5.

Figure 10: The time variation of power coefficient of the considered cases with , and 4.0.

Figure 11 shows the time variation of of the above two considered cases. The same as the power coefficient, the amplitude of the of the wing with has larger amplitude than the wing with ; however, comparing with the pitching motion, the two considered wings have almost identical phase difference.

Figure 11: The time variation of of the considered cases with and 4.0, .

Figure 12 shows the time variation of lift coefficient of the two considered wings with and 4.0, . It is seen from this figure that contrary with the time variation of , the wing with has larger amplitude of lift coefficient than the wing with .

Figure 12: The time variation of lift coefficient of the considered cases with and 4.0, .

Figure 13 shows the vorticity contours of the two considered wings with and 4.0, during a semipitching cycle. Obviously, in this figure, there exists a vortex separated for the wing with (see ), while for the wing with , no obviously separated vortex is observed, which results that the wing with has smaller mean power coefficient. It is also found that the passive plunging amplitude of the wing with has larger value than the wing with . According to equation (11), when the increased passive plunging amplitude of the wing is larger than the increased mean power coefficient, the efficiency of the wing will fall. This is the reason why the wing with has smaller mean power coefficient but has larger power efficiency.

Figure 13: The vorticity contours of the considered cases with and 4.0, .

3.2. Effect of Mass Ratio

To examine the effect of mass ratio, the spring coefficient is fixed at 1.00, which has the best energy extraction performance under different damping coefficient as section above discussed.

Figure 14 plots the mean power coefficient and power efficiency against damping coefficient with different mass ratio . There are three interesting conclusions that can be concluded from this figure. Firstly, at a given , both mean power coefficient and power efficiency increase first then decrease; moreover, the value of at which the wing achieves maximum mean power coefficient is smaller than the wing at which achieves maximum power efficiency. Secondly, the maximum mean power coefficient and power efficiency, respectively, are almost identical for the wing with different ; however, the value of at which the wing achieves the maximum mean power coefficient and power efficiency decreases with increasing. Thirdly, for the wing with smaller , both the mean power coefficient and power efficiency are insensitive to the variation of damping coefficient , which indicates that the wing with smaller has better working stability that, therefore, leads the wing to have better power extraction performance.

(a)

(b)

Figure 14: The mean power coefficient and power efficiency against damping coefficient with different mass ratio ; (a) mean power coefficient; (b) power efficiency.

To explore the mechanism of how the mass ratio influence the power extraction performance of the semiactive flapping wing in detail, two typical cases with , and , are considered. Figure 15 shows the time variation of power coefficient of the above two considered cases. It is clear in this figure that for the wing with , , it has larger power coefficient than the wing with , almost the whole pitching cycle, and it also has earlier time to achieve maximum power coefficient.

Figure 15: The time variation of power coefficient of the considered cases with , and , .

Figure 16 shows the time variation of of the above two considered cases. The wing with , , has larger amplitude of than the wing with , ; meanwhile, the phase difference between pitching and of the two considered wings is identical; the value is 0.30T.

Figure 16: The time variation of of the considered cases with , and , .

Figure 17 shows the time variation of lift coefficient of the above two relation cases. Different from the time variation of power coefficient and the passive plunging motion , the wing with , , has smaller amplitude of lift coefficient than the wing with , ; meanwhile, comparing to the wing with , , there exists an obvious phase delay for the wing with , , to achieve maximum lift coefficient.

Figure 17: The time variation of lift coefficient of the considered cases with , and , .

Figure 18 shows four instantaneous vorticity contours of the above considered wings during a semipitching cycle. It is seen from this figure that the mass ratio can influence the flow field of the wing significantly. Similar as the wing with different damping coefficient (as shown in Figure 13), there exists a vortex separated for the wing with (see t = 0.10T to 0.30T), while for the wing with , no obviously separated vortex is observed, which results the wing with has smaller mean power coefficient. It is also found that the passive plunging amplitude of the wing with has larger value than the wing with . According to equation (11), when the increased passive plunging amplitude of the wing is larger than the increased mean power coefficient, the efficiency of the wing will fall. This is the reason why the wing with has smaller mean power coefficient but has larger power efficiency.

Figure 18: The vorticity contours of the considered cases with , and , .

4. Conclusion

In this paper, the power extraction performance of the semiactivated flapping wing with different damping coefficient, spring coefficient, and mass ratio is numerically examined. A numerical code based on Finite Volume method to solve the two-dimensional Navier-Stokes equations and coupled with Finite Center Difference method to solve passive plunging motion equation is developed.

We first examine the effect of damping coefficient and spring coefficient and fix the mass ratio at 60. It is found that an optimal set of spring coefficient () is obtained, for which both the high mean power coefficient and power efficiency are achieved no matter what the damping coefficient value of the wing is. Second, the effect of mass ratio is examined, and the spring coefficient is fixed at 1.00. It is concluded that the variation of mass ratio cannot increase the maximum mean power coefficient and power efficiency of the semiactivated flapping wing, but it can influence the value of damping coefficient at which the wing achieves the maximum mean power coefficient and power efficiency; moreover, insensitivity of the mean power coefficient and power efficiency to the variation of damping coefficient is observed for the wing with smaller mass ratio, which indicates that the wing with smaller has better working stability. More specifically, the flow field around the wing is also investigated, and for the wing with appropriate damping coefficient, spring coefficient, and mass ratio (, , and ), the separated vortex can reattach to the wing surface during flapping, which results in no obvious vortex disappeared in the wake. This is the reason for the wing to have better power extraction performance.

Notations

d: | Chord length of wing
: | Mass of the airfoil
: | Damping coefficient
: | Spring constant
: | Pitching motion
: | Passive plunging motion
: | Pitching amplitude
: | Pitching frequency
: | The pitching period
: | The lift force
: | Damping coefficient
: | Spring coefficient
: | Mass ratio
: | Fluid density
: | The wing density
: | Area of the wing
Re: | Reynolds number
: | The stream velocity
: | The fluid kinematic viscosity
: | Reduced frequency
: | The power extracted from the fluid
: | The power needed to input the system
: | The fluid torque
: | The pitching velocity
: | The passive plunging velocity
: | The net power
: | The power coefficient and mean power coefficient
: | The power extraction efficiency
: | The overall vertical extent of the wing motion
: | The velocity vector
: | The pressure.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflict of interest.

Acknowledgments

This study was funded by the National Natural Science Foundation of China (project No. 51505347) and Key Laboratory of Metallurgical Equipment and Control of Education Ministry, Wuhan University of Science and Technology Foundation (2015B07).

References

A. B. Rostami and M. Armandei, “Renewable energy harvesting by vortex-induced motions: review and benchmarking of technologies,” Renewable & Sustainable Energy Reviews, vol. 70, pp. 193–214, 2017.

J. Young, J. C. S. Lai, and M. F. Platzer, “A review of progress and challenges in flapping foil power generation,” Progress in Aerospace Sciences, vol. 67, pp. 2–28, 2014.

M. Pourmahdavi, M. N. Safari, and S. Derakhshan, “Numerical investigation of the power extraction mechanism of flapping foil tidal energy harvesting devices,” Energy & Environment, vol. 30, no. article 0958305X1878732, pp. 193–211, 2019.

A. Javed, M. Jamil, B. Ali, K. Djidjeli, and T. A. Shams, “Performance analysis of flapping foil flow energy harvester subjected to non-sinusoidal pitching motion,” in 2018 Fluid Dynamics Conference, Atlanta, GA, USA, 2018.

J. Zhu and T. Tian, “The time asymmetric pitching effects on the energy extraction performance of a semi-active flapping wing power generator,” European Journal of Mechanics - B/Fluids, vol. 66, pp. 92–101, 2017.

W. Mckinney and J. Delaurier, “Wingmill: an oscillating-wing windmill,” Journal of Energy, vol. 5, no. 2, 115 pages, 1981.

Q. Zhu, “Optimal frequency for flow energy harvesting of a flapping foil,” Journal of Fluid Mechanics, vol. 675, no. 675, pp. 495–517, 2011.

Q. Xiao and Q. Zhu, “A review on flow energy harvesters based on flapping foils,” Journal of Fluids and Structures, vol. 46, no. 2, pp. 174–191, 2014.

J. Wu, Y. L. Chen, and N. Zhao, “Role of induced vortex interaction in a semi-active flapping foil based energy harvester,” Physics of Fluids, vol. 27, no. 9, article 093601, 2015.

J. Wu, Y. L. Qiu, C. Shu, and N. Zhao, “Pitching-motion-activated flapping foil near solid walls for power extraction: a numerical investigation,” Physics of Fluids, vol. 26, no. 8, article 83601, 2014.

J. Zhan, B. Xu, J. Wu, and J. Wu, “Power extraction performance of a semi-activated flapping foil in gusty flow,” Journal of Bionic Engineering, vol. 14, no. 1, pp. 99–110, 2017.

Q. Zhu, M. Haase, and C. H. Wu, “Modeling the capacity of a novel flow-energy harvester,” Applied Mathematical Modelling, vol. 33, no. 5, pp. 2207–2217, 2009.

J. Deng, L. Teng, D. Pan, and X. Shao, “Inertial effects of the semi-passive flapping foil on its energy extraction efficiency,” Physics of Fluids, vol. 27, no. 5, pp. 53103–53191, 2015.

H. Abiru and A. Yoshitake, “Study on a flapping wing hydroelectric power generation system,” Journal of Environment and Engineering, vol. 6, no. 1, pp. 178–186, 2011.

M. A. Ashraf, J. Young, and J. C. S. Lai, “Reynolds number, thickness and camber effects on flapping airfoil propulsion,” Journal of Fluids and Structures, vol. 27, no. 2, pp. 145–160, 2011.

F. M. Bos, D. Lentink, B. W. van Oudheusden, and H. Bijl, “Influence of wing kinematics on aerodynamic performance in hovering insect flight,” Journal of Fluid Mechanics, vol. 594, pp. 341–368, 2008.

Z. Jianyang, J. Lin, and H. Yu, “Aerodynamic performance of the three-dimensional lumped flexibility bionic hovering wing,” International Journal of Aerospace Engineering, vol. 2018, Article ID 1314258, 14 pages, 2018.

Q. Zhu and Z. Peng, “Mode coupling and flow energy harvesting by a flapping foil,” Physics of Fluids, vol. 21, no. 3, article 033601, 2009.

L. Teng, J. Deng, D. Pan, and X. Shao, “Effects of non-sinusoidal pitching motion on energy extraction performance of a semi-active flapping foil,” Renewable Energy, vol. 85, pp. 810–818, 2016.