1. Introduction
Flapping airfoils offer the possibility of understanding how both thrust and lift can be generated using solely one surface. Applications of unsteady airfoils span from aircraft elasticity [
1] to marine engineering [
2], but when looking at the lower Reynolds number range, a connection is seen with natural flight.
The first explanation of how flapping airfoils produce thrust and lift came from Knoller [
3] and Betz [
4], who saw that, when an airfoil oscillates, it presents an effective angle of attack that generates an aerodynamic force that has the thrust and lift components. The main limitation of the pioneering theory was that it was only based on the airfoil kinematics, neglecting any real-fluid effects. The work of von Kármán and Burgers [
5] offers a complete explanation based on the vortices of the airfoil’s wake and their orientation. One of the first mathematical theories regarding unsteady airfoils was made by Theodorsen [
6], who predicted the unsteady aerodynamics of pitching and plunging airfoils. Garrick [
7] extended the study to the propulsion capabilities of a movable wing with three degrees of freedom. Nowadays, with the available computational power, numerical methods have evolved to offer flow solutions based on Navier–Stokes equations.
The study of unsteady airfoils encircles flow structures, propulsive performance, energy harvesting, and bionics. Within these topics, flapping is the most studied, followed by pure pitching, and lastly, pure plunging [
8]. Pure plunging studies are generally shadowed by more complex kinematics, but the combined motion of plunging and pitching is not always a guarantee of higher thrust production. However, it often leads to a better propulsive efficiency [
9].
An alternative way to improve the performance of unsteady airfoils is by modifying the leading-edge geometry. This may be accomplished using concepts such as DDLE (Dynamically Deformed Leading-Edge) or VDLE (Variable Droop Leading-Edge), which are commonly used for dynamic stall mitigation [
10]. Innovative concepts, such as rotating cylinders at the leading-edge, have been used as active flow control techniques to prevent flow separation [
11]. Others use bioinspired solutions, for instance, alula deployment at pre-stall conditions to significantly improve lift [
12] or adding an active elastic leading-edge deformation that can improve airfoils with energy extraction purposes [
13]. While these are examples of active systems that deliver aerodynamic and/or propulsive improvements, other investigations have employed passive solutions, which generally offer lower energy consumption and a lower system complexity. Liu et al. [
14] attached a rigid flat plate to the trailing edge of a foil that exhibited passive motion determined by the fluid–structure interaction. The results showed that with an adequate level of flexibility, this modification was responsible for at least half of the thrust production. A similar concept was used by Tian et al. [
15], who studied the power-extracting capabilities of an oscillating airfoil. However, the authors concluded that the passive deformation could not reach the efficacy hoped without some level of active control.
Geissler and van der Wall [
16] identified that flow separation limits the flapping-wing efficiency and proposed a time-dependent deflectable leading-edge that reduces the LEV (Leading-Edge Vortex) size. The authors tested a plunging airfoil at a Mach of
,
, and a reduced frequency equal to
and searched for improvements in the propulsive coefficients. The results showed that the propulsive efficiency of the deformable leading-edge with a droop amplitude of
was roughly doubled compared to the rigid airfoil. Niu et al. [
17] made changes to the leading-edge by using a VDLE to improve the aerodynamic characteristics of a pitching NACA0012 airfoil. By deflecting the leading-edge, adverse pressure gradients were lowered due to the lower local angle of attack near the leading-edge. This effect increased the lift and reduced the drag and pitching moment at higher pitch angles.
The VDLE concept was also used by Chandrasekhara et al. [
18] to control dynamic stall in the compressible regime. They used a VR-12 airfoil with a movable leading-edge as it underwent a sinusoidal pitch oscillation. By drooping the leading-edge, the frontal part of the airfoil presented a lower incidence and the tendency of the airfoil to enter dynamic stall was reduced. The results showed that reductions of up to 75% were achieved, making the concept desirable for rotor dynamics, as it can eliminate torsional instabilities. Similarly, Xing et al. [
19] looked into the role of the leading-edge, proposing an inflatable frontal structure as a way to reduce dynamic stall phenomena.
With a micro air vehicle design in mind, Drost et al. [
20] studied the use of a dynamic leading-edge in a thin, flat, rigid plunging airfoil. By actively moving the leading-edge, the LEV at high angles of attack has its strength reduced. By moving the leading-edge dynamically, the results showed that the boundary layer and flow separation were altered significantly, increasing the lift-to-drag ratio. The authors indicated that this is an adequate alternative when it comes to eliminating the complexity of adding systems, such as blowing or suction.
In their work, Geissler and van der Wall [
16] mentioned that the influence of a dynamic leading-edge deformation in the domain of incompressible flows, where birds and insects live, is of high interest. To our knowledge, no published work has explored such a regime with an innovative geometry. Although dynamically changing the leading-edge has been widely studied in the field of flow separation control and dynamic stall, these works have generally focused on the higher Reynolds spectrum and rotor dynamics. Contrary to those, this study focused on improving the known propulsive characteristics of flapping airfoils, rather than looking at stall or flow separation prevention. Our preliminary studies [
21,
22] explored the low Reynolds regime by proposing a modified version of the NACA0012 airfoil with a movable leading-edge.
The present work focuses on studying the influence of a pitching leading-edge added to the traditional NACA0012 at a lower Reynolds state. The leading-edge pitching amplitude is changed parametrically under different oscillating conditions. The propulsive characteristics are evaluated using propulsive indicators, and associations between force generation, pressure distribution, and flow features are given. The real cause of propulsive augmentation is also exposed.
3. Results and Discussion
The propulsive capabilities of the NACA0012-IK30 were tested numerically at a Reynolds number of with two nondimensional plunging amplitudes, and , and three nondimensional velocities, , and . For each of these conditions, five leading-edge pitching amplitudes were considered ( 0, 5, 10, 15, and 20).
At first, the propulsive performance was evaluated using the three fundamental parameters commonly used for flapping airfoils: the mean propulsive power, the mean required power, and the propulsive efficiency. In
Figure 6, the mean propulsive power coefficient is plotted against the leading-edge pitching amplitude for the three nondimensional velocities considered. Overall, increasing the pitching amplitude was shown to be beneficial up to a certain amplitude, where
then started decreasing, except when
, where the maximum was not reached. The results also show that the slower flapping condition did not benefit from the dynamic leading-edge or presented good propulsive capabilities, being very close to a drag-producing regime. While both nondimensional amplitudes had fairly similar behaviors, with
, improvements in terms of magnitude were much more significant. When it came to the mean propulsive power, improvements of up to
relative to the plunging-only leading-edge were achieved. Additionally, the results indicate that the higher the nondimensional velocity is, the higher the leading-edge pitching amplitude needs to be to achieve the optimal propulsive power.
When it comes to the mean required power coefficient, as shown in
Figure 7, the influence of the dynamically activated leading-edge is beneficial overall but is rather small, except for the case with the lower nondimensional amplitude and highest nondimensional velocity. In this specific case, the required power decreases linearly as
increases, with a
reduction compared to the plunging-only leading-edge.
Combining the two propulsive coefficients results in the propulsive efficiency shown in
Figure 8. The graphs indicate that when the leading-edge is deflected, great improvements are achieved for both nondimensional amplitudes with increases of up to
. Overall, a higher propulsive efficiency is seen for the middle nondimensional velocity, although it is not clear whether
could offer better efficiency at larger pitching amplitudes, since no larger values were considered.
Based on these results, it is clear that adding a pitching leading-edge provides great advantages to the propulsive characteristics of the conventional NACA0012 plunging airfoil. These improvements are inevitably related to changes induced in the flow field by the dynamic leading-edge. In
Figure 9 and
Figure 10, the influence of the leading-edge pitching amplitude on the flow features is analyzed for all conditions presented before. The representation of the numerical field uses the total pressure contour to observe flow separation, since separated flow regions show lower values due to associated losses. Snapshots are given at the end of the oscillation period (
), as they provide a representation of what occurred during the entire period.
As expected, all conditions show some flow separation, as the airfoil is oscillating at nondimensional velocities that surpass its static stall angle (≈ 12). The cases with a plunging-only leading-edge experience a maximum effective angle of attack at between 14 and 45 (), which can be reduced by ≈ 6 when the leading-edge pitches at 20.
Still, even though both nondimensional amplitudes have the same effective angle of attack profile, the phenomena seen around the airfoil are different. When the airfoil oscillates with a larger nondimensional amplitude, flow features are scaled up with the leading-edge vortex occupying the whole upper surface of the airfoil when there is no deflection of the leading-edge.
For both amplitudes, as the nondimensional velocity increases, the presence of the leading-edge vortex becomes more evident and moves upstream. This is expected as the effective angle of attack that triggers separation is reached quicker, and thus, the LEV (leading-edge vortex) formation starts earlier. In a plunging airfoil, the leading-edge vortex only works as a thrust enhancement feature when located upstream of the maximum thickness position, described by Gao et al. [
32] as a thrust booster for plunging airfoils. The NACA0012-IK30 airfoil follows the same principle, even though it should not be considered to be a traditional plunging airfoil.
Focusing on the influence of the leading-edge pitching amplitude, we see that it has a crucial role in dictating the LEV size. The size and growing process of the leading-edge vortex are explained. Images show that the LEV presents its maximum size when there is no leading-edge deflection. With an increase in the pitching amplitude to 20, the vortex reduces its size to approximately 50 percent for both nondimensional amplitudes. This decrease is caused by the reduction in the adverse pressure gradient as a consequence of lowering the effective angle of attack of the leading-edge part.
As seen for the mean propulsive coefficients, the pitching amplitude has a small influence at the lowest nondimensional velocity, which is also seen in the flow fields obtained numerically. However, for the other two, the influence is much more evident. We verify that, when increasing the pitching amplitude, the leading-edge vortex present in becomes smaller and feeds another vortex that appears beyond the maximum thickness position. This feeding process alone is not beneficial for propulsive improvement, since it increases the suction force after the maximum thickness position, which inevitably increases drag.
At the same time, looking back at
Figure 6, we see that its concave-shaped function means that there is an optimal pitching amplitude value that maximizes the mean propulsive power. For instance, at
, that value is somewhere between 10
and 15
, which causes flow separation in the flow field. This result is central to the analysis of the proposed mechanism, because it clearly indicates that we are not pursuing the total elimination of the LEV. Instead, we are looking for a way to exploit its presence. This is contrary to flow separation occurring at a higher Reynolds numbers, which is not desirable at all, meaning that understanding the LEV dynamics (its position, intensity, and residence time on the airfoil) is crucial to optimally use the NACA0012-IK30 as a thrust enhancement mechanism.
Additionally, determining the evolution of the propulsive coefficients over time is crucial to understanding where propulsive improvements come from and, more specifically, the time frames in which they occur. The following graphs are focused on the conditions that show good propulsive characteristics and a clear dependence on the leading-edge pitching amplitude, and .
In
Figure 11, the propulsive power is shown as a function of
for a nondimensional velocity of
. To simplify the analysis, only three leading-edge pitching amplitudes are presented, together with the changes in the effective angle of attack over time. Deflecting the leading-edge improves the propulsive power for the vast majority of the oscillating period with the greatest improvements seen when the airfoil has its maximum effective angle of attack (
and
). However, the results also show that a continuous increase in
may not be beneficial when the pitching amplitude can lead to the degradation of the propulsive power, as observed in
.
At the highest nondimensional velocity, shown in
Figure 12, propulsive power is highly favored by the leading-edge pitching amplitude. Once again, bigger improvements are seen when the airfoil is at its maximum plunging velocity. This coincides with the instant where the leading-edge has its maximum angle of attack, enlarging the frontal area of the airfoil, which leads to an increase in thrust. This is presented in the following text using the pressure distribution around the airfoil. For the range studied, higher pitching amplitudes lead to higher propulsive power values, although a pitching amplitude ceiling similar to the
case is expected, since excessive deflection presents no real benefit. This ceiling was not observed, solely due to the relative importance between the plunging and pitching terms of the effective angle of attack. The reduction in
that the leading-edge can provide is limited to 6
when the pitching amplitude is set at 20
. Thus, increasing the
makes the effective angle of attack induced by the plunging motion more dominant. This is the reason why we need higher pitching amplitudes at higher nondimensional velocities to achieve higher propulsive power values, as seen in
Figure 6.
Concerning the required power, it was seen that its mean value decreases with the leading-edge pitching amplitude but, overall, shows a small dependence on it, except in some cases at the lowest nondimensional amplitude. This is reinforced in
Figure 13, which shows that the deflection of the leading-edge indeed causes a reduction in power consumption during the acceleration phases of the period.
In
Figure 14, the power consumption for the highest nondimensional velocity is shown. Overall, the required power coefficient is similar to the previous case, having its biggest reduction at
and
, where the airfoil experiences a smaller effective angle of attack compared to the baseline case. Although this was also verified for
, for the rest of the period, the leading-edge deflection slightly increases the power consumption, counteracting those gains. This is why some of the cases analyzed in terms of the mean power consumption are almost independent of the pitching amplitude.
To better understand the impact of the dynamically activated leading-edge, some snapshots of the numerical flow field are presented, together with the airfoil pressure distribution. The flow fields are shown only for at where satisfactory propulsive improvements are seen for and . Flow visualization data are also provided for supplemental comparisons.
In
Figure 15, the conditions
and
are presented. The influence of the movable leading-edge on the flow field features is clear, especially regarding the dynamics of the leading-edge vortex and the pressure distribution around the airfoil. By increasing the leading-edge pitching amplitude, the suction peak located at the front decreases. Not only that, a low-pressure zone located beyond the maximum thickness position grows in intensity. For
, the intensity approaches the suction level present at the leading-edge, but interestingly, it does not reduce the propulsion considerably, as seen in
Figure 11 (right).
This is an effect that is worth taking some time to analyze. When looking at and , we verified that a lower suction peak at the leading-edge does not influence the suction force considerably when going from a pitching amplitude of 0 to 10 or even 20. How could this happen if the suction force located at the frontal part of the airfoil is a major contributor to thrust enhancement? This could only be justified by the enlarged frontal area created by the leading-edge deflection. Even if the suction peak becomes smaller, integrating it into a larger frontal area results in a larger or similar suction force, as seen by the propulsive power increase that was verified. This is also verified in the following text for .
In
Figure 16, the same plunging amplitude is considered, but the nondimensional velocity is increased to
. The results are similar to those for
, but the effects are magnified, as expected by the higher effective angle of attack imposed by the larger nondimensional velocity. However, contrary to the previous case, the suction peak is not greatly influenced by the leading-edge deflection, reinforcing the fact that propulsive augmentation comes mainly from the increase in the frontal area created by the pitching of the frontal part. Unlike the previous case, the leading-edge pitching amplitudes considered cannot mitigate LEV formation. For that, a higher
must be considered. As mentioned before, this is a condition where the plunging component of the motion dominates the effective angle of attack.
The presence of the leading-edge vortex is justified by a strong adverse pressure gradient located right at the leading-edge, which is seen in the pressure distribution graphs close to
where there is a rapid increase in pressure. This pressure gradient at the leading-edge can be alleviated using leading-edge deflection, as verified in
Figure 15 when
. However, it comes at the expense of enlarging the suction at the back part of the airfoil. Meanwhile, the size of the vortex its not controlled by the existing pressure gradient. Instead, it is governed by the angular deviation between the shear layer that appears following separation and the leading-edge chord direction, as illustrated in
Figure 17.
Figure 18 summarizes the influence of the pitching amplitude on the leading-edge vortex’s size and its intensity at the highest nondimensional amplitude and velocity. The frames shown were taken at the end of the oscillation period (
) when the LEV was fully developed, being a cumulative representation of what occurred during the full period.
Although the size of the leading-edge vortex is reduced considerably, by about 50%, as mentioned before, its intensity stays fairly unchanged, judging by the similar pressure distributions. This is also observed during the descending phase, as shown in
Figure 16. However, it is worth mentioning that this is a plunging-dominated condition and the suction does remain somewhat unaffected until a certain pitching amplitude, and from there on, the suction present at the leading-edge starts shrinking, as seen in
Figure 15.