Research Article  Open Access
Effect of Reynolds Number on Aerodynamics of Airfoil with Gurney Flap
Abstract
Steady state, twodimensional computational investigations performed on NACA 0012 airfoil to analyze the effect of variation in Reynolds number on the aerodynamics of the airfoil without and with a Gurney flap of height of 3% chord are presented in this paper. RANS based oneequation SpalartAllmaras model is used for the computations. Both lift and drag coefficients increase with Gurney flap compared to those without Gurney flap at all Reynolds numbers at all angles of attack. The zero lift angle of attack seems to become more negative as Reynolds number increases due to effective increase of the airfoil camber. However the stall angle of attack decreased by 2° for the airfoil with Gurney flap. Lift coefficient decreases rapidly and drag coefficient increases rapidly when Reynolds number is decreased below critical range. This occurs due to change in flow pattern near Gurney flap at low Reynolds numbers.
1. Introduction
A Gurney flap (GF) is a microtab fitted perpendicular to the airfoil near the trailing edge on its pressure surface which increases the lift by altering the Kutta condition and increasing effective camber. They are extensively used on helicopter stabilizers [1]. Considerable efforts are carried out to study the effects of Gurney flap on the airfoil aerodynamics [2–4]. Many investigations are undertaken to determine the effect of various parameters of Gurney flap such as height, position, mounting angle, and configurations. But most of the investigations are performed at high Reynolds numbers above critical range. Mueller and Batill [5], Selig et al. [6], Brown and Filippone [7], and Traub and Agarwal [8] have studied the performance of airfoils at low Reynolds numbers but a systematic investigation is still not available. For smaller UAVs and MAVs, airfoil chord and velocity of aircraft are lower; hence, they operate at low Reynolds numbers. Reynolds numbers between Re = 3 × 10^{4} and 7 × 10^{4} are of great interest to MAV designers [9]. Most UAVs operate at Re = 5 × 10^{4} and 3 × 10^{5} [10]. Gurney flaps are also used in turbomachinery to improve their performance at low Reynolds number due to their inexpensive method of lift enhancement [11].
A semiempirical formula linking flap height to free stream velocity and the airfoil chord is proposed by Brown and Filippone [7]:where is optimum GF height (mm), Ch is chord length (m), and is free stream velocity (m/s).
Keeping this in mind, computations on the effect of six Reynolds numbers from high Reynolds number to low Reynolds number well below the critical range are carried out to analyze the performance of the airfoil.
2. Computational Methodology
2.1. Geometry and Grid Generation
The airfoil considered in this study is NACA 0012 airfoil with chord length of 5 cm. Ctype domain and grid are created in ICEMCFD with farfield boundaries 12.5 chords away from trailing edge in all directions (Figure 1(a)). Grid is generated in ICEMCFD with at least 300,000 nodes. Flow around the airfoil with GF is highly complex with high intensity vortices. Hence, first layer with cell width of 0.025 mm is generated up to 3% chord length around the airfoil. Second layer with cell width of 0.035 mm up to next 7% of the airfoil chord is generated around the airfoil as illustrated in the blocking strategy shown in Figure 1(b). Thereafter cell growth ratio of 1.05 is set up to farfield. Wall values along the airfoil surface are kept below 1 to account for low Reynolds number flow regime. Boundary layer elements near the leading and trailing edges of the airfoil are shown in Figure 1(c).
(a)
(b)
(c)
2.2. Turbulence Modeling
Previous studies show that out of two equation RANS based models of FLUENT, RNG model is best suited for the present airfoil study [12, 13]. However these computations are carried out at a higher Reynolds number of 2.1 × 10^{6}. RNG theory accounts for low Reynolds number effects. However effective use of this feature depends on appropriate treatment [14]. Hence it was decided to explore oneequation SpalartAllmaras (SA) turbulence model for the present computations at low Reynolds numbers and compare these results with those of RNG model for improved predictions. Validation is done on an airfoil with 1 m chord length.
2.3. Boundary Conditions and Solver Settings
The airfoil boundary is assigned as solidwall with noslip condition while inlet is assigned as velocity inlet and outlet is assigned as pressureoutlet conditions. Density based implicit solving scheme is used with the flow medium being air and Mach number less than 0.3. Hence the fluid is assumed to be incompressible with constant density of 1.225 kg/m^{3} and dynamic viscosity of 1.7894 × 10^{−5 }kg/ms. Under relaxation factors for all the transport variables are set to 0.8. Solution initialization is computed from velocity inlet followed by FMG initialization with solution steering. Equations are solved until a convergence criterion of 10^{−5} for all the residuals is satisfied.
3. Results and Discussions
3.1. Validation
Lift and drag coefficients, and , for the airfoil of 1 m chord are computed with RNG and SA turbulence models and compared at high Reynolds numbers (2.1 and 3.0 × 10^{6}) and at a lower Reynolds number (5.0 × 10^{5}). The results are presented as a function of angle of attack in Figure 2. Comparison is also made with available experimental results at Re = 2.1 × 10^{6} [3].
For high Reynolds number, values obtained by both models are very close (Table 1). Masoud et al. [15] also found good agreement with experimental values of lift and drag coefficients for NACA 23018 airfoil using SpalartAllmaras (SA) turbulence model. At the lower Reynolds number, the values predicted by RNG model decrease rapidly more than SA model. This difference also increases at higher angles of attack as shown in Figure 2. Similar trends are noticed for drag coefficient. Both computational models predict higher drag than experimental results. However the values predicted by the SA model are close to the experimental values (Table 2) and almost 40%–50% lower than RNG model at Re = 5.0 × 10^{5}. Based on the validation studies, it was decided to use SA model for all computations at low Reynolds numbers. It is quite possible that improved predictions particularly drag coefficient may be obtained using advanced turbulence models. This is demonstrated by Arko and McQuilling [16], who used advanced turbulence models of Abe et al. [17, 18] and Kato and Launder [19]. However FLUENT does not have these turbulence models available.


3.2. Effect of Reynolds Number on Aerodynamics of the Airfoil without and with Gurney Flap
Computations are carried out for the airfoil without and with a GF of height of 3% chord at six different Reynolds numbers (Re = 3.0 × 10^{5}, 1.5 × 10^{5}, 1.0 × 10^{5}, 8 × 10^{4}, 5 × 10^{4} and 3 × 10^{4}). Based on the equation for optimum GF given by Traub and Agarwal [8], the optimum GF height varies from 2.7 to 4.2% of airfoil chord as Reynolds number decreases from 3.0 × 10^{5} to 3.0 × 10^{4}. Hence a GF of height of 3% of airfoil chord is chosen for computations at all Reynolds numbers. At low Reynolds numbers, SA turbulence model works better than RNG turbulence model. Hence the former model is used for all the following investigations.
3.2.1. Lift Coefficient
The lift coefficient versus AoA characteristics for the airfoil without and with Gurney flap is presented in Figure 3. At all Reynolds numbers, the lift coefficient for the airfoil with GF is always higher than the lift coefficient for the airfoil without GF at all angles of attack below the stall angle. The lift coefficient obtained at AoA = 8° for the airfoil with GF is higher than maximum lift coefficient obtained for the airfoil without GF. As Reynolds number decreases, the lift coefficient also decreases for both cases: without and with GF. The extent to which Reynolds number affects the performance varies for both the cases and is presented in Table 3. The zero lift angle seems to become more negative as Reynolds number increases (it is estimated to be −2 to −3° as Reynolds number increases from 3.0 × 10^{4} to 3.0 × 10^{5}). This suggests that the effect of GF is to increase the effective camber of the airfoil. However the effective camber decreases as Reynolds number decreases.

For the airfoil without GF, variation in values with Reynolds number is negligible at low angles of attack and increases slightly with increasing angle of attack. When Reynolds number is changed from 1.5 × 10^{5} to 1.0 × 10^{5}, C_{L} decrement for AoA = 4°, 8° and 12° is 1.2%, 1.6%, and 2.5%, respectively. However, for the airfoil with Gurney flap, decrement in is almost the same at lower and intermediate angles of attack. The lift coefficient decreases rapidly near stall. Changing the Reynolds number from 1.5 × 10^{5} to 1.0 × 10^{5}, decrement in values the airfoil with GF at AoA = 4° and AoA = 8° is 2.5%, whereas it at AoA = 12° is 4.2%.
Table 3 represents the percentage decrement in values for every 10,000 decreases in Reynolds number for different ranges of Reynolds number for the airfoil without and with GF. For AoA = 12°, percentage decrease in for the airfoil without GF for every 10,000 decreases in Reynolds number near critical Reynolds number range ( ≈ 1.5 × 10^{5}) is 0.25. Below critical range the percentage decrement increases from 0.6 for the Reynolds number range of 1.0 × 10^{5} to 8.0 × 10^{4} to 1.62 at very low Reynolds number range of 5.0 × 10^{4} to 3.0 × 10^{4}. For AoA = 12°, percentage decrease in for the airfoil with GF for every 10,000 decreases in Reynolds number near critical range is 0.42. Below critical range the percentage decrement increases from 1.17 for the Reynolds number range of 1.0 × 10^{5} to 8.0 × 10^{4} to 2.51 at very low Reynolds number range of 5.0 × 10^{4} to 3.0 × 10^{4}. Hence it may be concluded that Gurney flap has adverse effects on lift coefficient as Reynolds number is decreased. However the lift coefficient for the airfoil with GF is always higher than that for the airfoil without GF at all Reynolds number and at all angles of attack.
The variation of the lift coefficient, , as a function of Reynolds number at representative angles of attack is shown in Figure 4. This figure clearly shows the rapid decrease in below the critical Reynolds number. This figure also shows the lift enhancement capability of Gurney flap at all Reynolds numbers. It is also noted that the performance of the airfoil without GF is less affected by variation in Reynolds number.
3.2.2. Drag Coefficients
Apart from decreasing with decreasing Reynolds number, is also found to increase as shown in Figure 5. For the sake of clarity, drag coefficient is presented at four Reynolds numbers only, namely, 3.0 × 10^{5}, 1.5 × 10^{5}, 8.0 × 10^{4}, and 3.0 × 10^{4}. Magnitude by which increases with decreasing Reynolds number is nearly independent of angle of attack. Increase in is very steep below critical Reynolds number as clearly shown in Figure 6. For higher Reynolds number range, the slope of increment is more for high AoA. Apart from lift enhancement capability of Gurney flap, increased drag penalty is also clearly visible. The obtained at AoA = 12° for the airfoil without GF is obtained at only AoA = 8° with GF. However, the increase in is less for the airfoil with GF when compared to the airfoil without GF.
The percentage increment in values for every 10,000 decreases in Reynolds number for different ranges of Reynolds number for the airfoil without and with GF is presented in Table 4. For AoA = 12°, percentage increase in for the airfoil for every 10,000 decreases in Reynolds number near critical range is 1.39, whereas below critical range the value increases from 3.19 for the Reynolds number range of 1.0 × 10^{5} to 8.0 × 10^{4} to 10.07 at very low Reynolds number range of 5 × 10^{4} to 3 × 10^{4}. The respective values for the airfoil with GF are 1.3, 2.62, and 6.60. Hence it may be concluded that Gurney flap has certain beneficial effects on drag coefficient as Reynolds number is decreased.

3.2.3. LifttoDrag Ratio
The lifttodrag ratio is presented as / at different angles of attack in Figure 7. From the figure it is evident that the lifttodrag ratio for the airfoil with GF is higher than that for the airfoil without GF up to angle of attack of 8°. At this angle of attack, the lifttodrag ratios for the airfoil without and with GF are almost equal at all Reynolds numbers. Beyond this angle of attack, the lifttodrag ratio for the airfoil with GF is lower than that for the airfoil without GF, indicating additional drag penalty due to GF. The difference between the lifttodrag ratios for the airfoil without and with GF seems to be constant at all Reynolds numbers. However the difference between the lifttodrag ratios for the airfoil without and with GF at angle of attack of 14° is substantially higher than that 12°. From Figure 3, it is evident that the airfoil with GF stalls at angle of attack of 12° while the airfoil without GF stalls at angle of attack of 14°. This is the reason for the large difference at angle of attack of 14°. As discussed earlier, the effective camber of the airfoil with GF increases causing earlier stall.
The decrease in / ratio per every 10,000 decreases in Reynolds number for different ranges of Reynolds number is presented in Table 5. This decrement seems to be nearly constant for the airfoil without and with GF at all Reynolds number ranges. However as the Reynolds number range decreases the lift/drag ratio is decreasing at a faster rate. When Reynolds number decreases from 3.0 × 10^{5} to 1.5 × 10^{5}, the decrease in / ratio per every 10,000 decreases in Reynolds number varies from 0.262 to 0.375 for the airfoil without GF as the angle of attack increases from 4° to 12°. The corresponding values for the airfoil with GF are 0.255 and 0.363. When Reynolds number decreases from 3.0 × 10^{5} to 1.5 × 10^{5}, the corresponding values for the airfoil without and with GF are 1.161 and 1.256 and 1.551 and 0.910, respectively.

3.2.4. LifttoDrag Ratio versus Lift Coefficient
The lifttodrag ratio as a function of lift coefficient for the airfoil without and with GF is presented in Figure 8. For the sake of clarity, the lifttodrag ration is presented at four Reynolds numbers only (Re = 3.0 × 10^{4}, 8 × 10^{4}, 1.5 × 10^{5} and 3.0 × 10^{5}). The lifttodrag ratio increases as Reynolds number increases for both the airfoils without and with GF. At lower lift coefficients, there is a drag penalty associated with GF. This drag penalty increases with Reynolds number. At higher lift coefficients, the lifttodrag ratio increased. Hence the lift coefficient increases for a given lifttodrag ratio. However at the highest lift coefficient, the lifttodrag ratio is substantially reduced indicating large drag penalty as the airfoil approaches stall angle. This is true for the airfoil without GF and with GF at all Reynolds numbers.
3.2.5. Static Pressure Distribution
Static pressure distribution in terms of nondimensional coefficient is presented in Figure 9. For the sake of clarity, the static pressure distribution is presented for three Reynolds numbers only (Re = 3.0 × 10^{4}, 8 × 10^{4} and 1.5 × 10^{5}) for the airfoil with and without Gurney flap at AoA = 12°, which is close to the stall angle. For the airfoil without Gurney flap, static pressure distributions are almost overlapping with each other for different Reynolds numbers. The effect of Reynolds number is clearly visible for the airfoil with Gurney flap. As the Reynolds number decreases, suction and pressure are continuously decreasing all along both the respective surfaces.
Maximum pressure buildup is almost the same for the airfoils with and without GF. Value of maximum for variation in Reynolds number from Re = 3.0 × 10^{4} to 3.0 × 10^{5} varies from 1.11 to 1.02 for the airfoil without GF and from 1.09 to 1.02 for the airfoil with GF. The maximum suction decreases as Reynolds number decreases. However, the change is very small for the airfoil without GF. Maximum suction varies from −3.4 to −4.6 for Re = 3.0 × 10^{4} to 3.0 × 10^{5} for the airfoil without GF, whereas, for the same variation in Reynolds number, values vary from −4.5 to −6.8 for the airfoil with GF. On the pressure surface, the static pressure increases for the airfoil with GF, while, on the suction surface, the static pressure decreases. This occurs all along the surface resulting in substantial increase in lift coefficient for the airfoil with GF at all angles of attack and at all Reynolds numbers. Adverse pressure gradient occurs on the pressure surface near the trailing edge of the airfoil with GF. This adverse pressure gradient decreases with Reynolds number. Such adverse pressure gradients were experimentally observed by previous researchers [20]. A recirculating vortex that occurs just upstream of Gurney flap may be the reason for this adverse pressure gradient [2]. Static pressure coefficients are reduced near the Gurney flap as Reynolds number decreases. Unlike the sudden rise in pressure just before Gurney flap on pressure surface at high Reynolds numbers, rate of pressure buildup is not rapid at low Reynolds numbers.
3.2.6. Pathlines
Pathlines are highly affected by the variation in the Reynolds number and the formation of vortices experience drastic change. Pathlines superimposed with the contours of turbulent viscosity at AoA = 12° for the airfoil without and with GF at four different Reynolds numbers presented in Figure 10. Unlike at very high Reynolds number (of the order of Re = 10^{6}), only one vortex is completely formed in the wake of Gurney flap at Re = 3.0 × 10^{5} instead of two vortices. As the Reynolds number further decreases, the only vortex in wake also start to disappear.
A laminar separation bubble starts its formation. As the Reynolds number is decreased, the flow loses its ability to make transition into turbulent flow in the attached boundary layer, hence forming a laminar separation bubble. Below Re = 1.0 × 10^{5}, no vortex is present behind the Gurney flap. The increased pressure on pressure surface is only due to vortex ahead of the Gurney flap which might explain absence of sudden pressure increase near the Gurney flap as shown in static pressure distribution. This laminar separation bubble increases the effective thickness of the airfoil, thereby increasing the pressure drag over the region, which explains the increase in drag at low Reynolds number and degraded performance of the airfoil at low Reynolds numbers. For each Reynolds number, flow is turned towards the Gurney flap whereas, due to absence of any suction, the flow leaves at the airfoil at higher angle without Gurney flap. However the turning of flow towards GF is reduced at lower Reynolds numbers.
4. Conclusions
A computational investigation on the effects of Reynolds number on the aerodynamics of NACA0012 airfoil without and with Gurney flap of height of 3% airfoil chord has been carried out. ANSYS FLUENT commercial CFD code with oneequation SpalartAllmaras turbulence model is used for the six Reynolds numbers varying from 3.0 × 10^{5} to 3.0 × 10^{4}. From this investigation, the following major conclusions are drawn.(1)Reynolds number plays a very major role in the airfoil aerodynamics for the NACA0012 airfoil without and with Gurney flap. Lift decreases and drag increases when Reynolds number is decreased.(2)For the airfoil with GF, Reynolds number has adverse effects on lift coefficient, while drag coefficient of the airfoil with GF has some beneficial effects compared to the airfoil without GF.(3)For high Reynolds number above critical range, decrease in and increase in are negligible.(4)As the Reynolds number is decreased below the critical Reynolds number range, decreases and increases at very steep rate and the performance degrades rapidly.(5)For lower Reynolds numbers, the two vortices behind the Gurney flap vanish. The Gurney flap seems to increase the effective camber of the airfoil, causing negative zero lift angle and reduced stall angle.
Nomenclature
AoA, α:  Angle of attack (deg.) 
Drag coefficient =  
Ch:  Airfoil chord (m) 
Lift coefficient =  
Static pressure coefficient = 2()/  
Drag force (N)  
Optimum Gurney flap (m)  
Gurney flap height as percentage of chord  
Lift force (N)  
Static pressure (Pa)  
Re:  Reynolds number = Ch/ν 
Freestream velocity (m/s)  
Rate of dissipation of turbulent kinetic energy (m^{2}/s^{3})  
Turbulent kinetic energy (m^{2}/s^{2})  
Air density (kgm/m^{3})  
Kinematic viscosity (m^{2}/s). 
2:  Exit. 
CFD:  Computational fluid dynamics 
GF:  Gurney flap 
RNG:  RNG turbulence model 
RANS:  Reynolds averaged Navier Stokes equations 
SA:  SpalartAllmaras turbulence model. 
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
References
 R. W. Prouty, Helicopter Aerodynamics, vol. 2, Philips Publishing, 1985.
 R. H. Liebeck, “Design of subsonic aerofoils for high lift,” Journal of Aircraft, vol. 1, no. 9, pp. 547–561, 1978. View at: Google Scholar
 Y. C. Li, J. J. Wang, and P. F. Zhang, “Effects of Gurney flaps on a NACA0012 airfoil,” Flow, Turbulence and Combustion, vol. 68, no. 1, pp. 27–39, 2002. View at: Publisher Site  Google Scholar
 J. J. Wang, Y. C. Li, and K.S. Choi, “Gurney flapLift enhancement, mechanisms and applications,” Progress in Aerospace Sciences, vol. 44, no. 1, pp. 22–47, 2008. View at: Publisher Site  Google Scholar
 T. J. Mueller and S. M. Batill, “Experimental studies of separation on a twodimensional the airfoil at low Reynolds numbers,” AIAA journal, vol. 20, no. 4, pp. 457–463, 1982. View at: Publisher Site  Google Scholar
 M. S. Selig, J. F. Donovan, and D. B. Fraser, Airfoils at Low Speeds, Soartech 8, Soartech Publications, Virginia Beach, Va, USA, 1989.
 L. Brown and A. Filippone, “Aerofoil at low speeds with Gurney flaps,” The Aeronautical Journal, vol. 107, no. 1075, pp. 539–546, 2003. View at: Google Scholar
 L. W. Traub and G. Agarwal, “Aerodynamic characteristics of a gurney/jet flap at low reynolds numbers,” Journal of Aircraft, vol. 45, no. 2, pp. 424–429, 2008. View at: Publisher Site  Google Scholar
 T. J. Mueller, “Aerodynamic measurements at low Reynolds numbers for fixed wing microair vehicles,” in Proceedings of the Special Course on Development and Operations of UAVs for Military and Civil Applications, RhodeSaintGenèse, Belgium, September 1999. View at: Google Scholar
 G. K. Ananda, P. Sukumar, and M. S. Selig, “Lowtomoderate aspect ratio wings tested at low Reynolds numbers,” AIAA Paper 3026, 2012. View at: Google Scholar
 M. Suresh and N. Sitaram, “Gurney flap applications for aerodynamic flow control,” in Proceedings of the 9th International Conference on Mechanical Engineering, Paper no. ICME 11FL40, Dhaka, Bangladesh, December 2011. View at: Google Scholar
 K. Sriram, S. Jain, and N. Sitaram, “Grid and turbulence model based exhaustive analysis of NACA 0012 airfoil,” Journal of Advanced Research in Applied Mechanics and Computational Fluid Dynamics, vol. 1, no. 1, pp. 13–18, 2014. View at: Google Scholar
 S. Jain, N. Sitaram, and S. Krishnaswamy, “Computational investigations on the effects of gurney flap on airfoil aerodynamics,” International Scholarly Research Notices, vol. 2015, Article ID 402358, 11 pages, 2015. View at: Publisher Site  Google Scholar
 Ansys, ANSYS Fluent Theory Guide, Ansys, Canonsburg, Pa, USA, 2011.
 J. N. Masoud, S. Habibollah, and A. Mekanik, “Numerical investigation of aerodynamic characteristics of NACA 23018 airfoil with a gurney flap,” International Journal of Mechanical Engineering and Robotics Research, vol. 1, no. 3, pp. 341–349, 2012. View at: Google Scholar
 B. M. Arko and M. McQuilling, “Computational study of highlift lowpressure turbine cascade aerodynamics at low reynolds number,” Journal of Propulsion and Power, vol. 29, no. 2, pp. 446–459, 2013. View at: Publisher Site  Google Scholar
 K. Abe, T. Kondoh, and Y. Nagano, “A new turbulence model for predicting fluid flow and heat transfer in separating and reattaching flows—I. Flow field calculations,” International Journal of Heat and Mass Transfer, vol. 37, no. 1, pp. 139–151, 1994. View at: Publisher Site  Google Scholar
 K. Abe, T. Kondoh, and Y. Nagano, “A new turbulence model for predicting fluid flow and heat transfer in separating and reattaching flowsII. Thermal field calculations,” International Journal of Heat and Mass Transfer, vol. 38, no. 8, pp. 1467–1481, 1995. View at: Publisher Site  Google Scholar
 M. Kato and B. E. Launder, “The modeling of turbulent flow around stationary and vibrating square cylinders,” in Proceedings of the 9th Symposium on Turbulent Shear Flows, pp. 10.4.1–10.4.6, Kyoto University, Kyoto, Japan, August 1993. View at: Google Scholar
 R. Myose, I. Heron, and M. Papadakis, “The poststall effect of gurney flaps on a NACA0011 airfoil,” in Proceedings of the Aerospace Atlantic Conference, SAE Paper 961316, Dayton, Ohio, USA, May 1996. View at: Google Scholar
Copyright
Copyright © 2015 Shubham Jain et al. 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.