3.1. Grid Independence and Flapping Cycle Convergence Tests
A grid independence test was performed for three patterns with grid sizes of
,
, and
for the representative length of the finest grid size in the multiblock to investigate the number of grid sizes with numerical reliability in this analysis. The representative speed of the flapping flight
was set to
, and the resolution was verified when the lift coefficient was stable in the seven flapping cycles of the flapping cycle convergence test.
Figure 8 shows the time history of the lift coefficient in seven cycles, and
Table 3 shows the cycle-averaged lift coefficient per cycle. As shown in
Figure 8, a disturbance was observed at the beginning of the cycle because of the beginning of the wing movement; however, it converged immediately, and little change was found in the time history. In addition, not much difference was observed after the fifth cycle (
Table 3). Therefore, the analysis conducted in this study was performed using the data from the sixth cycle.
For the grid independence test,
Figure 9 shows the time history of the lift coefficient
, the pressure component of the lift coefficient
, and the viscous stress component the of lift coefficient
over one cycle at each resolution, and
Table 4 shows the cycle-averaged values of them
,
, and
. Based on these results, the difference in lift at each resolution decreased as the resolution increased, and not much difference was observed between
and
. Therefore, a resolution of
was selected after considering the computational cost and accuracy.
3.3. Effect of Figure-Eight Motion
A comparison of the vortex structure and fluid forces generated by each motion was conducted to investigate the effect of the figure-eight motion on the aerodynamics. Each motion was represented by varying the elevation motion angle
which is a rotational motion that represents the figure-eight motion shown in Equation (4). Motion with figure-eight motion was represented by the elevation motion angle with the parameters shown in
Section 2.2, and motion without figure-eight motion was represented by
. Equations (41) and (42) show the angles with and without figure-eight motions as follows:
Figure 12 shows a schematic view of each motion represented by the flapping motion angle
and the feathering motion angle
shown in
Section 2.2 in addition to the elevation motion angle.
Figure 13 shows the isosurfaces of the vortex structure
formed via each motion at
,
,
, and
. The isosurfaces were colored on the basis of the normalized helicity density
. As shown in
Figure 13, the normalized helicity density was similar for each motion; however, the lengths of the vortex structures were different and the vortices with a figure-eight motion were longer than those without a figure-eight motion. This is because the vortices with a figure-eight motion stay longer on the top surface of the wings than those without a figure-eight motion, whereas the vortices without a figure-eight motion disappear immediately.
This difference in the vortex structure indicates a difference in the flow field, which indicates a difference in the relative velocity vector between the wings and the fluid. The difference in velocity vectors affects forces such as lift and thrust generated by the wings. Therefore, the time histories of the fluid forces generated via each motion were investigated.
Figure 14 shows the time histories of the lift coefficient
, the thrust coefficient
, and the power coefficient
over one stroke cycle for each motion, and
Table 5 shows the cycle-averaged values of coefficients
,
, and
and the ratio of the lift coefficient to the power coefficient
for each motion.
As shown in
Figure 14a, a comparison between the lift coefficient with and without figure-eight motions showed a different trend throughout the entire cycle and an essentially different trend in the lift coefficient in the first half of each stroke. We examined the vortex structure and pressure coefficient distribution at
, which corresponds to the initial stage of the downward motion of the wings (downstroke), to further discuss the cause of this large difference in the lift coefficient.
Figure 15 shows the isosurfaces of the vortex structure (
) formed by each motion. The isosurfaces were colored on the basis of the normalized pressure coefficient
. As shown in
Figure 15, the motion with a figure-eight motion produced a larger leading-edge vortex and wing tip vortex on the upper surface of the wings than the motion without a figure-eight motion. This difference could be attributed to the fact that the
-directional motion of the figure-eight motion increased the angle relative to the motion direction (angle of attack) of the wings by adding a downward motion in the first half of the downstroke. Consequently, a larger negative pressure was produced on the upper surface of the wings, and the pressure difference between the positive pressure on the lower surface and the negative pressure on the upper surface of the wings generated a larger lift force than the motion without a figure-eight motion. In addition, as shown in
Figure 14a, the relation of the lift coefficient in each motion reversed after the middle of the downstroke. This result was based on the abovementioned fact and was considered to be due to the fact that the figure-eight motion in the
-direction increases the angle of attack of the wings by adding upward motion after the middle of the downstroke. As shown in
Table 5, the average lift coefficient with a figure-eight motion was approximately
larger than that without a figure-eight motion.
As shown in
Figure 14b, a comparison between the thrust coefficient with and without figure-eight motions showed a similar trend throughout the entire cycle, and the cycle-averaged thrust coefficient in each was small (
Table 5). Therefore, the effect of the figure-eight motion on the thrust coefficient was small. This effect could be attributed to the fact that the thrust coefficient was offset by the upward motion of the wings (upstroke) and downstroke.
As shown in
Figure 14c, a comparison between the power coefficient with and without figure-eight motions showed a different trend throughout the entire cycle, and this difference could be attributed to the difference in the lift coefficient for each motion. As shown in
Table 5, a comparison of the cycle-averaged power coefficient in each motion confirmed that the motion with figure-eight motion required more power. This finding could be attributed to the motion in the
-direction.
As shown in
Table 5, the ratio of the lift coefficient to the power coefficient with a figure-eight motion was approximately 17% higher than that without a figure-eight motion. Therefore, these results indicated that a flight with a figure-eight motion consumed more power than a flight without a figure-eight motion but generated more lift, and it may be a more efficient way to fly while hovering.
3.4. Effect of Various Figure-Eight Motions and Reynolds Number
We compared the fluid forces generated via various figure-eight motions in addition to the motions shown in
Section 3.3 to investigate the more efficient figure-eight motion patterns in generating lift during the hovering motion and the relationship between figure-eight motion and Reynolds number. These motions were represented as eight types of figure-eight motions by varying the initial phase
of the elevation motion angle shown in Equation (4) from
to
in increments of
. The flapping motion angle and the feathering motion angle were the same as those shown in
Section 2.2.
Figure 16 shows a schematic view of each motion. Moreover, we investigated the dependence of the Reynolds number on the figure-eight motion.
Table 6 shows the cycle-averaged values of the lift coefficient and the ratio of the lift coefficient to the power coefficient for each motion at each Reynolds number.
Figure 17 shows the relation between the cycle-averaged power coefficient and the cycle-averaged lift coefficient.
Figure 17 shows the cycle-averaged power coefficient on the horizontal axis and the cycle-averaged lift coefficient on the vertical axis. Hence, comparing each type of motion, the motion that generates more lift for a certain amount of power (more efficient motion) is located in the upper left corner. This positional relationship through the comparison of motions can be described similarly by the slope of a line connecting the origin and one point of the results. Thus, the motion with a greater slope of the line is more efficient in generating lift. The dotted line shown in
Figure 17 was used to compare each figure-eight motion and without a figure-eight motion, which was the line connecting the value of without-figure-eight motion and the origin of the graph. In this study, the line was known as “without-8-line”. The motion located to the upper left of “without-8-line” generated lift more efficiently than without a figure-eight motion, whereas the motion located to the lower right of “without-8-line” was less efficient than without a figure-eight motion.
As shown in
Figure 16, the motions with
and
were recognized as U-shaped motions. However, as shown in
Table 6, the most efficient motions were not the U-shaped motions but the figure-eight motions with
at any Reynolds number. Therefore, in this study, the figure-eight motion was defined as the motion with
, without discussing the U-shaped motion independently.
As shown in
Table 6, the motion with a smaller initial phase of the elevation motion angle (
) had a higher ratio of the lift coefficient to the power coefficient, and the motion could generate lift more efficiently. The motion with
at
and
at other Reynolds numbers was the most efficient in generating lift. As shown in
Figure 16, this result could be attributed to the fact that the angle of attack was larger than the original figure-eight motion with
by increasing the vertical motion of the wings in each stroke. On the contrary, the motion with
had a larger vertical motion and a larger angle of attack compared with the motion with
; however, the upward motion was larger. Therefore, the motion was less efficient than the motion with
when it was averaged over the entire stroke. However, horizontal forces dominate rather than vertical forces toward motion at low Reynolds numbers; therefore, increasing the vertical motion of the flapping motion and the angle of attack are important for generating more lift. Therefore, the motion with
at
was the most efficient in generating lift. Moreover, the motion patterns of the most efficient in generating lift at each Reynolds number were investigated in detail.
Figure 18 shows the relation between the initial phase of the elevation angle of the most efficient motion in generating lift and Reynolds number. Note that the most efficient motion in generating lift was calculated from the tangential point where
was the greatest slope on the elliptical approximation shown in
Figure 17. As shown in
Figure 18, the initial phase of the elevation angle of the most efficient motion in generating lift increased as the Reynolds number decreased, whereas the angle decreased and approached zero as the Reynolds number increased.
In addition, as shown in
Table 6 and
Figure 17f, all motions generated more lift as the Reynolds number increased, and more lift was generated compared with the power. The above discussion is based on nondimensionalized values without considering differences in the speed of the wings (flapping frequency) to simplify the comparison of the lift and power generated by each motion. Then, the physical values of the lift and power coefficients were summarized to investigate the effect of the flapping frequency on aerodynamic characteristics.
Table 7 shows the cycle-averaged values of the force in the
direction
, power
, and the ratio of the lift coefficient to the power coefficient
generated via the motion without a figure-eight motion at each Reynolds number. As shown in
Table 7, a higher flapping frequency generated more power than lift.
Moreover, as shown in
Figure 17, the eight types of figure-eight motions showed an elliptical distribution in the relation between the lift coefficient and power coefficient. Therefore, we discuss the relation between the lift coefficient and power coefficient by making an elliptical approximation along the distribution. As shown in
Figure 17a–e, the area surrounded by the ellipse and the W8 line (light blue shown in the figure) decreased as the Reynolds number increased, and the number of motions that generate lift more efficiently than without the figure-eight motion also decreased. This decrease indicates that the effect of figure-eight motion decreased as the Reynolds number increased. The elliptical shape, which is primarily related to the decrease in area, was focused on to discuss this difference in effect in more detail. Therefore, the effect of the change in the Reynolds number on the figure-eight motion was determined by calculating the aspect ratio and inclination angle of the ellipse at each Reynolds number. As shown in
Figure 19, the long and short sides of the ellipse for each Reynolds number were defined as
and
, respectively, and the aspect ratio
can be expressed as
, where
is the ratio of the long side to the short side of the ellipse. In addition, the inclination angle of the ellipse and the location of the center of the ellipse in
Figure 17 were defined as
and
.
Table 8 shows the values of the elliptic approximation, and
Figure 20 shows the relation between the aspect ratio of the ellipse and the Reynolds number.
As shown in
Table 8 and
Figure 20, the inclination angle of the ellipse did not change remarkably; the aspect ratio increased with increasing Reynolds number, and the aspect ratio converged when the Reynolds number increased above a certain value. The increase in the aspect ratio indicated that the elliptical approximation of various figure-eight motions approached the straight “without-8-line” shown in
Figure 17, which indicates that the difference between the motions with and without a figure-eight motion decreased. In addition, the increase in the Reynolds number was synonymous with the increase in insect size. Therefore, varying the elevation motion angle that represents the figure-eight motion was not very effective for insects under high Reynolds number conditions (large size). Particularly, the figure-eight motion may be a vital mechanism for insects to generate lift more efficiently under low Reynolds number conditions (small size).
In this study, insect wings were analyzed as rectangular rigid plates without thickness. Regarding the shape, Kirishna et al. [
38] reported that flight efficiency did not change much depending on the shape of the wings between a rectangular model with a model based on the actual shape of a blowfly. However, the changes in the shape of the wings may affect the efficiency of each figure-eight motion pattern. Therefore, in future studies of insects, it is necessary to evaluate figure-eight motion by considering not only the shape of the wings but also the flexibility of the wings and other parameters.