Next Article in Journal
Finite Element Analysis of Pre-Stressed Ultra High-Performance Concrete (UHPC) Girders
Next Article in Special Issue
A Computationally Time-Efficient Method for Implementing Pressure Load to FE Models with Lagrangian Elements
Previous Article in Journal
Comparative Study on Housing Defect Repair Cost through Linear Regression Model
Previous Article in Special Issue
WebTraceSense—A Framework for the Visualization of User Log Interactions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Improved Lift for Thick Flatback Airfoils in the Inboard Blades of Large Wind Turbines

Department of Energy, Systems, Territory and Constructions Engineering, University of Pisa, Largo Lucio Lazzarino, 56122 Pisa, Italy
*
Author to whom correspondence should be addressed.
Eng 2024, 5(3), 2345-2361; https://doi.org/10.3390/eng5030122
Submission received: 12 August 2024 / Revised: 17 September 2024 / Accepted: 19 September 2024 / Published: 20 September 2024
(This article belongs to the Special Issue Feature Papers in Eng 2024)

Abstract

:
Thick airfoils are often used in the inboard sections of blades in commercial wind turbines. The main reason for this is to give the blade greater structural strength, but it is well known that thick airfoils degrade aerodynamic performance by stalling at relatively small angles of attack. The adoption of flatback airfoils instead of sharp trailing edges allows high lift coefficient to be maintained in thick airfoils. In this paper, we propose a novel airfoil design based on a passive flap to further improve the lift coefficient. This new design was tested by numerical simulation on several airfoils with different maximum thickness and different TE thickness. The improved design for flatback airfoils yields a higher lift coefficient, while the drag behaviour is strictly related to the baseline airfoil shape: some airfoils show a decrease in drag at certain angles of attack, while others exhibit a drag increase. In conclusion, the practical implications of the flap’s utilisation on a state-of-the-art blade designed for a 5 MW wind turbine are analysed. The findings demonstrate that, due to the enhanced lift coefficient, it is feasible to shorten the chord while maintaining the power output, thereby reducing material costs.
Keywords:
flatback; 2D; URANS; NREL 5MW; BEM

1. Introduction

Wind energy is a crucial component of clean energy production and reducing greenhouse gas emissions. To achieve the climate targets, the European Union, for example, aims to cover more than one third of its electricity demand with wind power by 2030 and over 40% by 2050. The trend in recent years has been for wind turbines to increase in size. In fact, to be economically viable, wind turbines should increase their power output, which is proportional to the swept area of the rotor. This requires turbines with increasingly longer blades. At the same time, the mass of the blade increases approximately as the cube of the blade length, resulting in greater gravitational loads, particularly torsional, edgewise bending and flapwise bending loads (as explained in [1]). The latter is the most significant load, and to prevent the turbine from structural problems, it is necessary to adopt thick airfoils in the area of the blade close to the root. A thick airfoil (i.e., airfoil with t / c 0.25 [2], where t is the maximum thickness and c is the airfoil chord) can be utilised to increase the resistance of the blade to load due to the enhanced moment of inertia. The downside is that the increase in thickness reduces the aerodynamic performance of the sharp trailing edge (TE) airfoil, since it stalls at small angles of attack, thereby preventing the achievement of a high lift coefficient C L . As a compromise between structural requirements and aerodynamic performance, blunt trailing edge or flatback (FB) airfoils offer a good solution: they increase the maximum lift coefficient and the lift coefficient slope [3]. This is due to the different pressure distribution along the airfoil compared to the sharp trailing edge airfoil, resulting in a lower adverse pressure gradient on the suction side, which prevents premature boundary layer separation. In addition, flatback airfoils offer some transport advantages [4,5]. They are used in the inboard area where the chord is at its highest, so they allow for smaller dimensions.
On the other hand, flatback airfoils show a higher drag coefficient ( C D ) compared to sharp TE and a wake development characterised by vortex shedding. However, when looking at the turbine as a whole, one of the most important parameters is the torque coefficient. Since the inboard region of the blade has a lower moment arm and lower relative velocity, the contribution to the torque generated by the entire blade is relatively small. Moreover, the kinetic energy available in the wind in a ring of thickness δ r is proportional to the area of the ring, thus exhibiting a decrease towards the root. However, increasing the aerodynamic performance of the inboard region could be beneficial.
Considering Figure 1, a generic blade element located at a radial position r is shown. The infinitesimal lift δ L and drag δ D forces are highlighted as well as the relative velocity to the blade W. The infinitesimal torque and thrust acting on the blade element can be obtained considering the total force along the tangential (t) and normal (n) directions, respectively:
δ Q = ( δ L · s i n ( ϕ ) δ D · c o s ( ϕ ) ) r = 1 2 ρ c W 2 C Q r δ r
δ T = δ L · c o s ( ϕ ) + δ D · s i n ( ϕ ) = 1 2 ρ c W 2 C T δ r
where ϕ is the inflow angle, the sum of the angle of attack α and the twist angle β , and C Q and C T are the torque and thrust coefficients, respectively, defined as follows:
C Q = C L · s i n ( ϕ ) C D · c o s ( ϕ )
C T = C L · c o s ( ϕ ) + C D · s i n ( ϕ )
The main contribution to the torque in the inboard region of the blade is due to the lift coefficient C L : near the root, the inflow angle ϕ is large, so the dominant term in the torque coefficient C Q shown in Equation (3) is that which depends on C L (as observed in [2,6,7]). A high C D is therefore acceptable when combined with a high C L , hence the adoption of flatback airfoils even when they have a high C D .
The aim of this paper is to provide a solution to improve the lift coefficient of a flatback airfoil while keeping the drag coefficient almost the same. The latter, as mentioned, has little influence on the torque at the root of the blade but has a negative effect on the thrust coefficient C T Equation (4). Consequently, it is necessary to identify novel geometries that will not substantially enhance the drag.
A review of the literature reveals numerous studies that have sought to reduce the drag coefficient of flatback airfoils. The primary objective of these studies was to reduce the drag coefficient and subsequently minimise noise emission due to vortex shedding. Various passive devices, including splitters and cavities (as shown in [8,9,10]), reveal a common trend in the aerodynamic performance resulting in a reduction in drag at relatively low angles of attack, accompanied by a concomitant reduction in the maximum lift coefficient achieved. Other works adopt a wavy TE such as in [11], where the drag reduction is accompanied by a lift loss. Another interesting solution is proposed in [4], where the swallow tail concept is studied. The trend shown in these works is that it is possible to attain a drag reduction paying the cost of a lift decrease. In the aforementioned works, it is not clear whether the lift losses are fully compensated for by the drag reduction or whether an increase in chord length may be necessary to maintain the same power output.
Other studies instead aim at increasing the lift coefficient; one example is given in [12] where the Gurney flap and vortex generators are analysed. The Gurney flap is a device that is typically set at a right angle at the TE on the pressure side of the airfoil. Its function is to enhance the camber of the airfoil, which results in a higher maximum lift coefficient and a reduction in the zero-lift angle. Vortex generators, on the other hand, delay the separation of the flow on the suction side of the airfoil, allowing for a higher maximum lift coefficient. These devices are beneficial for the lift, but a significant increase in drag is often linked to the increase in lift. The present paper has a similar objective: to enhance the lift coefficient of thick flatback airfoils. This may result in a slight increase in the drag coefficient, but given our focus on large-scale turbines, specifically those used in offshore applications, we will disregard the noise implications. The goal is to develop more efficient airfoils to achieve higher power production or to attain the same power output with a shorter chord, reducing plant costs. The techniques studied in the literature and listed above are not specifically designed for FB airfoils (except for the swallow tail), resulting in a poor overall performance (considering both the lift and drag coefficient). The present study proposes a flap designed specifically for FB airfoils. It consist of a passive device, which represents an economical solution to control the flow. Furthermore, the device is fixed to minimise maintenance costs.
The paper is structured as follows: Section 2 shows the methodology adopted and describes the geometry of the proposed flap. Section 3 describes the validation process, which encompasses the verification of calculation grids and the setup. To assess the robustness of the methodology, validation was conducted on multiple flatback airfoils, varying the airfoil under analysis. Subsequently, the novel geometry was implemented in a number of airfoils, as detailed in Section 4. In Section 5, a case study is presented for analysis. This case study involves the inboard blade of a turbine from the literature, which was analysed using conventional flatback airfoils. The results of this analysis are then compared with the new geometry proposed in this work. The paper concludes with a summary of its principal findings in Section 6.

2. Methodology

The objective of this study is to assess the potential for enhancing the aerodynamic performance of flatback airfoils. In particular, the use of passive devices has been explored as a means of increasing the C L while minimising the impact of these devices on the C D . Given the pivotal role of the C L in determining the torque generated at the root of the blade, small differences and increments in C D are deemed acceptable. The method used to enhance the lift is a flap of a specific design. The latter is depicted in grey in Figure 2. This kind of device can be added to any flatback airfoil; the slope of the flap is given by the line passing through points 1 and 2 (red in the figure), with point 1 located on the pressure side of the airfoil at a distance of 0.2 % c from the TE, while point 2 is on the TE. The extent of the flap is 5% of the chord length in the x direction, and the thickness is 2% of the chord length.
In order to conduct this study, 2D fluid dynamic simulations were employed. All computational grids were generated with the ANSYS ICEM v19.3 software and were structured multi-block grids. Such grids are a combination of the O-grid and the C-grid concept. This approach allows the advantages inherent to both grid types to be exploited. As documented in the literature [13], C-grids are well-suited for sharp TE airfoils, but they are less effective for finite thickness TE airfoils. Conversely, O-grids are adept at capturing the behaviour of flatback airfoils, but they tend to give poor predictions in the development of the wake and consequently in the C D coefficient. Indeed, O-grid results in an excessively coarse grid in the wake region. Figure 3 shows a generic grid used in this work for flatback airfoils: a C-grid is visible with an O-grid surrounding the airfoil with a thickness of 2 % c , which is enough to encompass the boundary layer thickness (the grid details in the red and violet boxes). The number of grid points on both the suction and pressure sides is 350, and there are 100 points on the TE. The wall distance from the first layer of cells is set to 1 × 10 5 c ; this ensures that y + values are largely less than 1 throughout the airfoil wall. The computational domain extends 33 c in the flow direction and 28 c in the normal direction: the distance from the inlet to the airfoil nose is 14 c as shown in Figure 3 on the left side. Figure 4 illustrates the grid employed in this study for the analysis of all FB airfoils with the flap. It differs from the grid utilised for flatback airfoils (as shown in Figure 3) in that an additional C-grid is incorporated at the TE in order to reproduce more complex geometry at the TE. The entire dominion extents are the same as those of Figure 3 on the left side. The boundary conditions used are a velocity inlet where the x velocity component is specified at the inlet and a pressure outlet at the exit where a free-stream pressure is set. The inlet has a C shape and also encompasses the top and bottom side; therefore, in order to set the free-stream velocity, we have to specify the velocity components: the x-component was set to U x , which is the free-stream velocity, while the y-component was zero. The numerical simulations were carried out with ANSYS Fluent v19.3 using URANS: a pressure-based approach was used with the SIMPLE scheme and default settings (also default controls). The turbulence model is k ω SST, and the convergence criteria were 1 × 10 4 for all residuals. A total of 8000 time steps were used for each simulation, and the simulation time was such that the entire domain was traversed twice from the air flow. Therefore, the time step size was defined as Δ t = 2 L x / ( U x 8000 ) , where L x is the domain length in the x direction.

3. Validation

Since the new flap proposed in this work is intended to be added to any flatback airfoil, the validation of the computational grid and the computational setup was carried out using three different airfoils to ensure a wide range of validity. The first airfoil is the one used in [14] called DU97FB (a flatback airfoil developed by modifying the DU97-W-300 airfoil); the second is called DU97ST and it is the airfoil used in [14] with the “swallow tail” (ST) concept (first appearance in [4]); while the third is the flatback airfoil proposed in [8], called FB-3500-1750. In Figure 5, the x-y coordinates of the profile are shown. The three airfoils are very different: DU97FB and DU97ST have a maximum thickness of 30 % c and a TE thickness of 10 % c , while the FB-3500-1750 airfoil has a maximum thickness of 35 % c and a TE thickness of 17.5 % c . Additional information about these validation cases is listed in Table 1. As previously stated, the simple flatback airfoils were evaluated using a grid type comparable to that illustrated in Figure 3, whereas FB airfoils with modifications at the TE (DU97ST in this case) were examined using a grid type analogous to that depicted in Figure 4.
Figure 6 illustrates the lift and drag coefficients for all the airfoils analysed in the validation task. Comparisons are made between the 2D simulations of the present study and the corresponding experimental measurements. With regard to airfoils DU97FB and DU97ST, the experimental data are available in [14], while for FB-3500-1750, the free-transition measurements available in reference [15] were considered. The general trend of 2D URANS simulations is an overestimation of both the lift and drag coefficient (as also observed in [16]). Furthermore, it is evident that 2D results have difficulty in accurately capturing the occurrence of stall. A similar trend for FB-3500-1750 airfoil is shown in [17]. In particular, the linear trend is accurately represented with regard to the C L , but the simulations overestimate both the maximum C L and the alpha at which it occurs, resulting in a delayed stall. Therefore, when using 2D URANS it is important to operate at angles of attack far enough from the stall condition. For C D , the simulated data consistently exceed the experimental data. These behaviours are consistent with those observed in other 2D numerical studies, as in [11,18,19]. Moreover, 2D URANS simulations have exhibit weakness in drag prediction as shown in [20], where a significant spanwise change in drag was measured.

4. Results

The flap geometry described in Figure 2 was applied to the test airfoils DU97FB and FB-3500-1750. Figure 7 illustrates the lift and drag coefficients obtained with and without the flap. It can be observed that the use of the flap allows for a higher lift coefficient across all angles of attack for both airfoils under consideration. With regard to the drag coefficient, the DU97FB airfoil exhibits a deterioration, manifested as an increase in drag across all considered angles of attack, with the only exception being for α = 0 ° . Conversely, the FB-3500-1750 airfoil benefits from the use of the flap at low angles of attack (up to 12°). Figure 8 illustrates the pressure coefficient C p for both airfoils with and without the flap, at two angles of attack, α = 0 ° and α = 12 ° . The pressure coefficient is defined as follows:
C p = p p 1 2 ρ U 2
where p is the static pressure on the blade, p is the undisturbed pressure and ρ is the air density. It is evident that the incorporation of the flap results in a higher overpressure on the pressure side of the airfoil, accompanied by a reduction in pressure on the suction side. The enhanced pressure difference across the airfoil is responsible for the observed lift enhancement. Conversely, the drag reduction observed at certain angles of attack can be attributed to the interference with the vortex shedding caused by the flap. The pathlines of the static pressure of DU97FB and FB-3500-1750 are plotted in Figure 9 and Figure 10, respectively. The figures illustrate the situation with and without the flap at two different angles of attack, α = 0 ° and α = 12 ° , and due to the oscillating trend of the C L (and C D ) caused by the vortex shedding, the figures represent the instantaneous snapshot at the C L peak. It is evident that the presence of the flap consistently confers an advantage to FB-3500-1750 across all angles of attack. The reduction in flow disturbance in the wake due to vortex shedding is a notable benefit of the flap. In the case of DU97FB, this remains true at α = 0 ° , while the situation with the flap at α = 12 ° exhibits a slight deterioration. This is accompanied by an increase in drag, as illustrated in Figure 7. A similar pattern is evident in Figure 11 and Figure 12, where the root-mean-square error (RMSE) of the mean velocity is plotted. It can be observed that the presence of the flap generally reduces velocity fluctuation in the wake, with the exception of DU97FB at α = 12 ° . The preceding renderings were obtained through the utilisation of the data-sampling option in Fluent, given that vortex shedding is a periodic phenomenon, exhibiting varying frequencies across different scenarios, i.e., different airfoils under consideration. Hence, to ensure a fair comparison, the pressure coefficient of Figure 8 and the velocity fields of Figure 11 and Figure 12 were calculated as the mean value of further 4000 time steps of simulation.
In order to facilitate a comparative analysis between the novel flap geometry proposed in the present work and a widely utilised flap geometry, namely the Gurney flap (GF), the GF was applied to the simple DU97FB flatback airfoil. In this instance, the Gurney flap was scaled to the same dimensions as the proposed flap, with a thickness of 2%c and a vertical length of 5%c. Figure 7(top) presents a comparison of the lift and drag coefficients obtained with the addition of the Gurney flap and with our flap. The GF generates a notable increase in the drag coefficient at all alpha values, while the lift coefficient is comparable to that of our flap up to α = 8 ° . For larger alpha values, the lift coefficient is higher when using the GF, with the exception of very high alpha values (e.g., α = 20 ° ). Figure 13 compares the streamlines colouring by static pressure at α = 10 ° for the simple flatback DU97FB airfoil, the airfoil with the flap and the airfoil with the GF. It is evident that the flap is capable of extending the overpressure region on the pressure side of the airfoil while maintaining relatively low drag values. The wake development in the case of the simple DU97FB (Figure 13a) and in the case with the flap (Figure 13b) are observed to be quite similar. In the case of the GF (Figure 13c), it is evident that the GF enhances the overpressure on the pressure side of the airfoil, but simultaneously increases the drag. The wake displays more pronounced depression regions. Therefore, the flap geometry proposed in this work, in contrast to the GF, allows for improvements in the lift coefficient while limiting the drag coefficient increment in comparison to the simple flatback airfoil case.
The domain extent in the y direction was maintained at 28 c for all simulations presented in this paper. To ensure that no significant boundary effects can occur with such a domain width, we conducted additional simulations, doubling the y extent. As illustrated in Figure 14, the new domain extends 47 c in the flow direction and 56 c in the normal direction. In particular, the new domain contains the old 28 c domain, as illustrated in Figure 14, where the blue line encompasses the 28 c domain. Consequently, the grid refinement in the calculation domain delimited by the blue line remains unchanged. The consideration of a larger domain has resulted in a slight increase in the number of grid cells, from 180 k to 200 k. In order to evaluate the influence of the domain y extent, the DU97FB airfoil was analysed at several angles of attack. The results demonstrate that enlarging the domain has no significant impact on the lift and drag coefficients, as illustrated in Figure 15. This is probably due to the fact that the grids used in this work have high refinement near the walls, so the grid surrounding the airfoil is very fine. Consequently, local effects are of greater consequence than far-field effects.
The flap geometry was then applied to two further airfoils: the DU40FB and the DU35FB airfoils. These airfoils are derived from the sharp trailing edge airfoils called DU40 and DU35, which are described in detail in [21,22] (coordinates available in [23] and in [24]). The DU40 and DU35 airfoils are obtained in turn from DU 99-W3-450 and DU 99-W3-350, respectively. All of the aforementioned airfoils are sharp TE airfoils, and the flatback airfoils were optimised in [25], where the x,y coordinates are provided for analysis. As would be expected, these airfoils have a maximum thickness of 40% and 35%, respectively, and a TE thickness of approximately 11% and 9%. The simulations for both airfoils were performed at R e of 7 × 10 6 , with an undisturbed flow velocity of 34.4 m/s and a chord length of 3 m. All other settings are the same as the previous simulations. Figure 16 illustrates the lift and drag coefficients. For these airfoils, the flap exerts a beneficial influence on lift at all angles of attack, whereas the drag exhibits disparate behaviour. The DU40FB airfoil with the flap demonstrates a general increase in drag. In contrast, the DU35FB airfoil with the flap exhibits a drag reduction up to α = 10 ° .
A uniform conclusion regarding the flap’s influence on drag cannot be drawn from the analysis of the results presented in Figure 7 and Figure 16. The flap confers an advantage on some airfoils at relatively low angles of attack, while this is not the case for others. This behaviour appears to be independent of the thickness of the TE or the ratio between the maximum thickness of the airfoil and the TE thickness. Therefore, the limited number of airfoils tested precludes the formulation of a general rule for leveraging the flap solution in terms of drag reduction. However, the flap proposed in this work appears to offer a promising solution when compared to other devices. For instance, as discussed in reference [12], the adoption of the Gurney flap or vortex generators has been shown to significantly enhance lift, particularly in the case of Gurney flaps. However, this increase in lift is often accompanied by a significant increase in drag. In contrast, the flap proposed in this work has been observed to result in a reduction in drag at certain angles of attack.
Given the inherent limitations of 2D URANS analysis, namely the overestimation of drag and the delay in stall prediction, it is recommended that the proposed flap geometry be subjected to further analysis using a more sophisticated tool, such as 3D simulations or wind tunnel tests. However, without considering quantitative data, it is possible to make a qualitative assessment of the flap geometry proposed in this paper, which suggests that it has the potential to enhance the performance of the airfoil, particularly by increasing the lift coefficient.
In other words, a comprehensive campaign of 2D URANS fast simulations could serve as the initial phase of an airfoil optimisation methodology. However, the final verification should be conducted with more reliable CFD techniques, despite the increased computational cost.

5. Case Study

The DU40FB and DU35FB airfoils were selected for the purpose of analysing a portion of a realistic blade. The turbine in question is the NREL 5MW, for which the distribution of airfoils, chord and twist are available for reference in [21,26]. The focus of this analysis was on the portion of the blade from μ = 0.18 to μ = 0.38 (where μ is the dimensionless radius along the blade defined as μ = r / R and r is the local radius along the blade, while R is the turbine radius that in this case is 63 m). It should be noted that the original blade was constructed with sharp trailing edge airfoils; however, for the purposes of this analysis, these were substituted with FB airfoils. As previously stated, the DU40FB and DU35FB airfoils were optimised in [25], beginning with the sharp TE.
The objective of this section is to undertake a comparative analysis of the same portion of blade utilising both the flatback airfoils and the modified geometry with flaps. The analysis of the blade portion was based on blade element/momentum (BEM) theory. In particular, in order to calculate the torque Q and thrust T generated by the blade, it is necessary to calculate the relative velocity to the blade W (see Figure 1) as follows:
W = U 2 ( 1 a ) 2 + r 2 Ω 2 ( 1 + a ) 2
a = 1 3
a = a ( 1 a ) λ 2 μ 2
λ = R Ω U
where a is the axial induction factor, a is the tangential induction factor, λ is the Tip Speed Ratio (TSR) and Ω is the rotation speed of the turbine. Equation (7) and consequently the a value represent the condition of maximum power extraction. The angle of inflow ϕ and the angle of attack α are calculated as follows:
ϕ = arctan ( 1 a λ μ ( 1 + a ) )
α = ϕ β
In order to evaluate and compare the torque and thrust of the “blade FB” and the “blade flap”, where “blade FB” denotes the blade obtained using simple FB airfoils, while “blade flap” denotes the blade using FB airfoils with the flap, it is first necessary to ascertain the chord distribution. Given the analytical nature of this analysis, it is essential to pay close attention to certain aspects. To ensure a fair comparison in the torque and thrust calculation, it is necessary to use a different chord distribution when changing the airfoil under consideration. The theoretical blade design procedure proposed in [27] is based on the assumption of optimised operating conditions, which are represented by Equations (7) and (8). In this analysis, the fundamental parameter is the geometry parameter, defined as follows:
B c 2 π R λ C L = 4 λ 2 μ 2 a ( 1 a ) 2 + ( λ μ ( 1 + a ) ) 2
Equation (12) demonstrates that the chord distribution and the lift coefficient employed for the design are inversely proportional. Consequently, an airfoil with a higher lift coefficient will result in a shorter chord at each blade section. In light of the aforementioned concepts, modifying the airfoil entails modifying the chord distribution. In Equation (12), the chord distribution is a function of solely the C L . The C D is not considered, which is an acceptable simplification given that in [27] it was demonstrated that an almost equal chord distribution would have been obtained considering also the C D . The drag effect will be evident only on the turbine performance, in particular, on the torque and thrust generated. If we assume that the value of C L remains constant from the tip to the root, and thus the α , then the blade will exhibit a strongly non-linear tapering for μ 0.6 . This is because the chord will result in very high values towards the root. Moreover, if a constant value of C L is assumed, the blade will exhibit a high degree of twist, with a significant difference in the β value between the root and the tip. From an economic standpoint, the design procedure entails a linear tapering from the tip to the root, whereby the chord at the root is reduced and the loss is compensated for by an increase in the value of C L and consequently in α . The left member of Equation (12) ( B c 2 π R λ C L ) necessitates a specific product of c · C L to achieve a given power extraction. Consequently, it is possible to reduce the chord when increasing the value of C L , and vice versa. For the present analysis, we take R, the β distribution, λ and the rated wind speed from the real turbine proposed in [26]. Therefore, using equations from (6) to (11), it is possible to find the α distribution along the blade. The procedure is as follows:
  • Given R, λ , the β -distribution and the U of the NREL 5MW turbine ( R = 63 m, λ = 7 , and U = 11.4 ), we calculate α ;
  • With the found α -distribution, we obtain the C L -distribution;
  • With the C L -distribution, we obtain the new chord distribution through Equation (12);
The ϕ -distribution is shown in Figure 17 (bottom right side) along with the α and β distributions. Using Equation (12), we found the chord distribution shown on the right side of Figure 17 (in the red box). The discontinuity at μ 0.25 is a consequence of the rapid change in the airfoil. For μ 0.25 , the airfoil employed is DU40FB (or DU40FB with the flap), whereas otherwise the airfoil used is DU35FB (or DU35FB with the flap). In the case of a real blade, it would be advisable to incorporate a fitting to soften the impact of the airfoil change (see, for example, the airfoil distribution suggested in [28]). Figure 17 on the left illustrates the optimal chord distribution when utilising the DU35FB airfoil up to the tip of the blade. However, the present study focuses exclusively on the section of the blade between μ 0.18 and μ 0.38, as delineated by the red box.
The alpha value exhibits a notable decline as the μ value rises, which is responsible for the observed increase in the chord for μ 0.25 .
Using Equations (3) and (4), it is possible to assess the torque and thrust coefficients obtained using the new designed blade section. Figure 18 illustrates the torque and thrust coefficients. As anticipated, the flap’s deployment results in an elevated thrust coefficient, attributable to the enhanced lift coefficient and the concomitant increase in drag coefficient at specific angles of attack. With respect to the torque coefficient, the flap enables an increased value across all radial positions, substantiating the aforementioned assertion that at the inboard blade, the dominant influence on the torque coefficient is the lift coefficient, despite the rise in drag coefficient.
By integrating Equations (1) and (2) along the blade section and multiplying by the number of blades, it is now feasible to obtain the torque and thrust, which are shown in Figure 19. The observed trend is consistent with the theoretical prediction, exhibiting a growth pattern from the root to the tip. The integration of the Q and T along the blade section yields an increase in total torque by 4% for the “blade flap” compared to the “blade FB” and at the same time a decrease in total thrust by 1%. It should be noted that all BEM calculations used to determine the Q and T values are based on the airfoil lift and drag coefficients presented in Figure 16, which were calculated at R e = 7 × 10 6 . The Reynolds number along the blade sections varies approximately between 6 × 10 6 and 10 × 10 6 . Despite the simplifications employed in the analysis of this case study, the qualitative benefits of flap adoption are evident. It may be possible to construct the blade section with a shorter chord, reducing material requirements and costs while maintaining power generation and thrust on the structure. This can be achieved compatibly with structural requirements at the root.

6. Conclusions

In this study, we propose a novel flap geometry with the objective of enhancing the lift coefficient of flatback airfoils. The latter are employed at the inboard region of the blade in order to satisfy both structural and aerodynamic requirements. The flap can be added to any FB airfoil, and it typically results in an increase in the lift coefficient. With regard to the drag coefficient, no uniform trend is observed, as it is strictly dependent on the baseline airfoil shape. At specific angles of attack, some airfoils exhibit a reduction in drag when the flap is adopted, while other airfoils show a general increase in drag. The flap was evaluated on a number of airfoils, namely, DU97FB, FB-3500-1750, DU40FB, and DU35FB. However, due to the limited number of airfoils tested, no general conclusions can be drawn. Further analyses are required to assess the flap proposed in this work, as the selected airfoils exhibit significant variation in maximum thickness, thickness of the TE and camber.
The observed behaviour of the flap results in an increase in lift coefficient for all the airfoils, which coincides with the objective of this study. In the inboard region of the blade, where the torque coefficient is mainly influenced by the lift coefficient, enhancing the lift coefficient of the airfoil chosen for the blade design is of paramount importance. As a case study, we have attempted to evaluate the potential of adopting the flap in the inboard region of a realistic blade turbine. The method adopted is analytical, specifically based on BEM evaluations. The blade under consideration is the NREL 5MW turbine. The original blade features a sharp TE; however, a blade with FB airfoils (“blade FB”) and a blade with the same airfoils but with the addition of the flap (“blade flap”) were compared. In order to ensure a fair comparison through analytical evaluations, the chord distribution used for the two blades was modified. In particular, the airfoils with a higher C L , i.e., airfoils with the flap, exhibit a shorter chord. The utilisation of the flap enables the attainment of a higher torque coefficient whilst concurrently a higher thrust coefficient. The total torque and thrust (Q and T) on the inboard region of the blade were evaluated, resulting in a +4% and −1% change, respectively. This is due to the fact that the BEM analysis enables the redesign of the blade, in particular, the chord distribution, as a result of the change in airfoil performance. The chord length exerts an influence on both Q and T; consequently, despite the fact that the thrust coefficient of the “blade flap” case is higher, the total thrust is lower due to the shorter chord. Neglecting quantitative results, the important thing in this case study was to evaluate the potential of adopting the flap in the inboard region of a realistic blade. What we can conclude is that it is possible to obtain approximately the same power output using fewer materials because of a shorter chord.
Future development of this work could be to further analyse the flap behaviour with more sophisticated tools, for instance, LES simulations or testing a couple of airfoil geometries with the flap addition in a wind tunnel. For a realistic blade inboard region, it will be useful to simulate an entire blade with 3D fluid dynamic simulations in order to assess the chord reduction. With 3D simulations, neglecting the BEM analytical approach, it will be possible also to analyse and compare airfoils with and without the flap using the same chord distribution, evaluating whether the flap could lead to higher power output.

Author Contributions

Conceptualisation, M.P. and S.Z.; methodology, M.P. and S.Z.; software, M.P. and S.Z.; validation, M.P. and S.Z.; formal analysis, M.P.; investigation, M.P. and S.Z.; data curation, M.P.; writing—original draft preparation, M.P.; writing—review and editing, M.P.; supervision, S.Z.; project administration, S.Z.; funding acquisition, S.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Recovery and Resilience Plan (NRRP), CN1, Centro Nazionale di Ricerca in High-Performance Computing, Big Data e Quantum Computing, spoke 6: multiscale modelling and engineering applications.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors on request.

Conflicts of Interest

The authors declare no conflicts of interest.

Abbreviations

The following abbreviations are used in this manuscript:
FBFlatback
GFGurney flap
TETrailing edge
TSRTip Speed Ratio

References

  1. Cox, K.; Echtermeyer, A. Structural design and analysis of a 10 MW wind turbine blade. Energy Procedia 2012, 24, 194–201. [Google Scholar] [CrossRef]
  2. Standish, K.; Van Dam, C. Aerodynamic analysis of blunt trailing edge airfoils. J. Sol. Energy Eng. 2003, 125, 479–487. [Google Scholar] [CrossRef]
  3. Baker, J.P.; Mayda, E.A.; van Dam, C.P. Experimental Analysis of Thick Blunt Trailing-Edge Wind Turbine Airfoils. J. Sol. Energy Eng. 2006, 128, 422–431. [Google Scholar] [CrossRef]
  4. Grasso, F.; Ceyhan, O. Non-conventional flat back thick airfoils for very large offshore wind turbines. In Proceedings of the 33rd Wind Energy Symposium, Kissimmee, FL, USA, 5–9 January 2015; p. 0494. [Google Scholar] [CrossRef]
  5. Barone, M.; Berg, D. Aerodynamic and aeroacoustic properties of a flatback airfoil an update. In Proceedings of the 47th AIAA Aerospace Sciences Meeting including The New Horizons Forum and Aerospace Exposition, Orlando, FL, USA, 5–8 January 2009; p. 271. [Google Scholar] [CrossRef]
  6. Jackson, K.; Zuteck, M.V.; Van Dam, C.; Standish, K.; Berry, D. Innovative design approaches for large wind turbine blades. Wind. Energy Int. J. Prog. Appl. Wind. Power Convers. Technol. 2005, 8, 141–171. [Google Scholar] [CrossRef]
  7. Law, S.; Gregorek, G. Wind Tunnel Evaluation of a Truncated NACA 64-621 Airfoil for Wind Turbine Applications; Technical Report; NASA Lewis Research Center: Cleveland, OH, USA, 1987. [Google Scholar] [CrossRef]
  8. Baker, J.P. Drag Reduction of a blUnt Trailing-Edge Airfoil. Ph.D. Thesis, University of California, Davis, CA, USA, 2009. [Google Scholar]
  9. Van Dam, C.; Cooperman, A.; McLennan, A.; Chow, R.; Baker, J. Thick airfoils with blunt trailing edge for wind turbine blades. In Proceedings of the Turbo Expo: Power for Land, Sea, and Air, Glasgow, UK, 14–18 June 2010; Volume 44007, pp. 923–931. [Google Scholar] [CrossRef]
  10. Kahn, D.L.; van Dam, C.; Berg, D.E. Trailing Edge Modifications for Flatback Airfoils; Technical Report; Sandia National Laboratories (SNL): Albuquerque, NM, USA, 2008. [Google Scholar] [CrossRef]
  11. Yang, S.J.; Baeder, J.D. Aerodynamic drag and aeroacoustic noise mitigation of flatback airfoil with spanwise wavy trailing edge. In Proceedings of the 33rd Wind Energy Symposium, Kissimmee, FL, USA, 5–9 January 2015; p. 0993. [Google Scholar] [CrossRef]
  12. Cene, A.; Manolesos, M.; Grasso, F. Aerodynamic and Aeroacoustic Measurements of the Flow Past a Very Thick Flatback Airfoil with Passive Flow Control Devices. In Proceedings of the AIAA SCITECH 2022 Forum, San Diego, CA, USA, 3–7 January 2022; p. 0279. [Google Scholar] [CrossRef]
  13. TPI Composites, Inc. Innovative Design Approaches for Large Wind Turbine Blades; Technical Report, SAND2004-0074; Sandia National Lab.: Albuquerque, NM, USA, 2004. [Google Scholar]
  14. Ceyhan, O.; Timmer, W. Experimental evaluation of a non-conventional flat back thick airfoil concept for large offshore wind turbines. In Proceedings of the 2018 Applied Aerodynamics Conference, Atlanta, GA, USA, 25–29 June 2018; p. 3827. [Google Scholar] [CrossRef]
  15. Mayda, E.A.; van Dam, C.P.; Chao, D.D.; Berg, D.E. Flatback Airfoil Wind Tunnel Experiment; Technical Report; Sandia National Laboratories (SNL): Albuquerque, NM, USA, 2008. [Google Scholar] [CrossRef]
  16. Bangga, G.; Sasongko, H. Dynamic stall prediction of a pitching airfoil using an adjusted two-equation URANS turbulence model. J. Appl. Fluid Mech. 2017, 10, 1–10. [Google Scholar] [CrossRef]
  17. Bangga, G.; Lutz, T.; Krämer, E. Numerical investigation of unsteady aerodynamic effects on thick flatback airfoils. arXiv 2017. [Google Scholar] [CrossRef]
  18. Sarlak, H.; Nishino, T.; Sørensen, J.N. URANS simulations of separated flow with stall cells over an NREL S826 airfoil. AIP Conf. Proc. 2016, 1738, 030039. [Google Scholar]
  19. Szydlowski, J.; Costes, M. Simulation of flow around a static and oscillating in pitch NACA 0015 airfoil using URANS and DES. In Proceedings of the Heat Transfer Summer Conference, Charlotte, NC, USA, 11–15 July 2004; Volume 46911, pp. 891–908. [Google Scholar]
  20. Metzinger, C.N.; Chow, R.; Baker, J.P.; Cooperman, A.M.; Van Dam, C. Experimental and computational investigation of blunt trailing-edge airfoils with splitter plates. AIAA J. 2018, 56, 3229–3239. [Google Scholar] [CrossRef]
  21. Kooijman, H.; Lindenburg, C.; Winkelaar, D.; Van der Hooft, E. DOWEC 6 MW Pre-Design: Aero-Elastic Modeling of the DOWEC 6 MW Pre-Design in PHATASDOWEC Dutch Offshore Wind Energy Converter 1997–2003 Public Reports. 2003. Available online: https://citeseerx.ist.psu.edu/document?repid=rep1&type=pdf&doi=b331d7f80401ca1058058f4d130c3246843a1989 (accessed on 11 August 2024).
  22. Blaylock, M.L.; Maniaci, D.C.; Resor, B.R. Numerical simulations of subscale wind turbine rotor inboard airfoils at low Reynolds number. In Proceedings of the 33rd Wind Energy Symposium, Kissimmee, FL, USA, 5–9 January 2015; p. 0493. [Google Scholar] [CrossRef]
  23. Blaylock, M.L.; Maniaci, D.C.; Resor, B.R. Numerical Simulations of Subscale Wind Turbine Rotor Inboard Airfoils at Low Reynolds Number; Technical Report; Sandia National Laboratories (SNL): Albuquerque, NM, USA, 2015. [Google Scholar] [CrossRef]
  24. Available online: https://github.com/old-NWTC/FAST/tree/master/CertTest/5MW_Baseline/Airfoils (accessed on 11 August 2024).
  25. Kim, S.H.; Bang, H.J.; Shin, H.K.; Jang, M.S. Composite structural analysis of flat-back shaped blade for multi-MW class wind turbine. Appl. Compos. Mater. 2014, 21, 525–539. [Google Scholar] [CrossRef]
  26. Jonkman, J.; Butterfield, S.; Musial, W.; Scott, G. Definition of a 5-MW Reference Wind Turbine for Offshore System Development; Technical Report; National Renewable Energy Lab. (NREL): Golden, CO, USA, 2009. [Google Scholar] [CrossRef]
  27. Burton, T.; Jenkins, N.; Sharpe, D.; Bossanyi, E. Wind Energy Handbook; John Wiley & Sons: Hoboken, NJ, USA, 2011. [Google Scholar]
  28. Griffith, D.; Richards, P.W. The SNL100-03 Blade: Design Studies with Flatback Airfoils for the Sandia 100-meter Blade; Technical Report; Sandia National Lab.(SNL-NM): Albuquerque, NM, USA, 2014. [Google Scholar] [CrossRef]
Figure 1. Qualitative representation of a blade element located at radial position r.
Figure 1. Qualitative representation of a blade element located at radial position r.
Eng 05 00122 g001
Figure 2. Qualitative representation of the flap proposed in this work. The zoom in the green circle highlights the geometric parameters of the flap (this is merely a qualitative representation; the image is not to scale in order to more clearly depict the flap details).
Figure 2. Qualitative representation of the flap proposed in this work. The zoom in the green circle highlights the geometric parameters of the flap (this is merely a qualitative representation; the image is not to scale in order to more clearly depict the flap details).
Eng 05 00122 g002
Figure 3. Example of a computational grid: (left) full domain extents and (right) detail of the grid near the blade (red rectangle) and near the TE (violet rectangle).
Figure 3. Example of a computational grid: (left) full domain extents and (right) detail of the grid near the blade (red rectangle) and near the TE (violet rectangle).
Eng 05 00122 g003
Figure 4. Example of the calculation grid used for FB airfoils with the flap (left), with the detail of the TE in the red rectangle (right).
Figure 4. Example of the calculation grid used for FB airfoils with the flap (left), with the detail of the TE in the red rectangle (right).
Eng 05 00122 g004
Figure 5. x-y coordinates of the three airfoils used for the validation task.
Figure 5. x-y coordinates of the three airfoils used for the validation task.
Eng 05 00122 g005
Figure 6. (Top) Lift and drag coefficient of the DU97FB and the DU97ST airfoils: comparison between the present 2D simulations and experimental measurements of [14]. (Bottom) Lift and drag coefficient of the FB-3500-1750 airfoil: comparison between the present 2D simulations and experimental measurements of [15].
Figure 6. (Top) Lift and drag coefficient of the DU97FB and the DU97ST airfoils: comparison between the present 2D simulations and experimental measurements of [14]. (Bottom) Lift and drag coefficient of the FB-3500-1750 airfoil: comparison between the present 2D simulations and experimental measurements of [15].
Eng 05 00122 g006
Figure 7. Coefficients of lift and drag: comparison between the simple flatback airfoil and the corresponding airfoil with the flap. DU97FB (Top) and FB-3500-1750 (Bottom). In case of DU97FB, the use of the Gurney flap (GF) is also analysed.
Figure 7. Coefficients of lift and drag: comparison between the simple flatback airfoil and the corresponding airfoil with the flap. DU97FB (Top) and FB-3500-1750 (Bottom). In case of DU97FB, the use of the Gurney flap (GF) is also analysed.
Eng 05 00122 g007
Figure 8. Pressure coefficient at α = 0 ° and α = 12 ° (Top) DU97FB with and without the flap (Bottom) FB-3500-1750 with and without the flap.
Figure 8. Pressure coefficient at α = 0 ° and α = 12 ° (Top) DU97FB with and without the flap (Bottom) FB-3500-1750 with and without the flap.
Eng 05 00122 g008
Figure 9. Pathlines coloured by static pressure at α = 0 ° and α = 12 ° for DU97FB with and without the flap.
Figure 9. Pathlines coloured by static pressure at α = 0 ° and α = 12 ° for DU97FB with and without the flap.
Eng 05 00122 g009
Figure 10. Pathlines coloured by static pressure at α = 0 ° and α = 12 ° for FB-3500-1750 with and without the flap.
Figure 10. Pathlines coloured by static pressure at α = 0 ° and α = 12 ° for FB-3500-1750 with and without the flap.
Eng 05 00122 g010
Figure 11. RMSE of the mean velocity magnitude at α = 0 ° and α = 12 ° for DU97FB with and without the flap.
Figure 11. RMSE of the mean velocity magnitude at α = 0 ° and α = 12 ° for DU97FB with and without the flap.
Eng 05 00122 g011
Figure 12. RMSE of the mean velocity magnitude at α = 0 ° and α = 12 ° for FB-3500-1750 with and without the flap.
Figure 12. RMSE of the mean velocity magnitude at α = 0 ° and α = 12 ° for FB-3500-1750 with and without the flap.
Eng 05 00122 g012
Figure 13. Streamlines coloured by static pressure of the DU97FB airfoil (a), DU97FB airfoil with a flap (b), and DU97FB airfoil with a Gurney flap (c) at α = 10 ° .
Figure 13. Streamlines coloured by static pressure of the DU97FB airfoil (a), DU97FB airfoil with a flap (b), and DU97FB airfoil with a Gurney flap (c) at α = 10 ° .
Eng 05 00122 g013
Figure 14. Enlarged calculation domain for sensitivity analysis. The new large domain includes the old domain, indicated by the blue line in the figure. The domain around the airfoil remains unchanged (red rectangle).
Figure 14. Enlarged calculation domain for sensitivity analysis. The new large domain includes the old domain, indicated by the blue line in the figure. The domain around the airfoil remains unchanged (red rectangle).
Eng 05 00122 g014
Figure 15. Coefficients of lift and drag: comparison between the “DU97FB” results (obtained with the 28 c large domain) and the “DU97FB larger domain” results (obtained with the 56 c large domain).
Figure 15. Coefficients of lift and drag: comparison between the “DU97FB” results (obtained with the 28 c large domain) and the “DU97FB larger domain” results (obtained with the 56 c large domain).
Eng 05 00122 g015
Figure 16. Coefficients of lift and drag: comparison between the simple flatback airfoil and the corresponding airfoil with the flap. DU40FB (Top) and DU35FB (Bottom).
Figure 16. Coefficients of lift and drag: comparison between the simple flatback airfoil and the corresponding airfoil with the flap. DU40FB (Top) and DU35FB (Bottom).
Eng 05 00122 g016
Figure 17. (Left) Chord distributions with zoom (red box) on the blade portion of interest. (Right) Chord distribution and ϕ , β and α distribution.
Figure 17. (Left) Chord distributions with zoom (red box) on the blade portion of interest. (Right) Chord distribution and ϕ , β and α distribution.
Eng 05 00122 g017
Figure 18. (Left) Torque coefficient along the inboard region of the blade and (Right) thrust coefficient.
Figure 18. (Left) Torque coefficient along the inboard region of the blade and (Right) thrust coefficient.
Eng 05 00122 g018
Figure 19. (Left) Torque along the inboard region of the blade and (Right) thrust.
Figure 19. (Left) Torque along the inboard region of the blade and (Right) thrust.
Eng 05 00122 g019
Table 1. Validation setup.
Table 1. Validation setup.
DU97FBDU97STFB-3500-1750
chord c [m]0.650.650.66
U x [m/s]252515
R e 1 × 10 6 1 × 10 6 6.66 × 10 5
M a 0.070.070.04
transitionfreefreefree
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Pucci, M.; Zanforlin, S. Improved Lift for Thick Flatback Airfoils in the Inboard Blades of Large Wind Turbines. Eng 2024, 5, 2345-2361. https://doi.org/10.3390/eng5030122

AMA Style

Pucci M, Zanforlin S. Improved Lift for Thick Flatback Airfoils in the Inboard Blades of Large Wind Turbines. Eng. 2024; 5(3):2345-2361. https://doi.org/10.3390/eng5030122

Chicago/Turabian Style

Pucci, Micol, and Stefania Zanforlin. 2024. "Improved Lift for Thick Flatback Airfoils in the Inboard Blades of Large Wind Turbines" Eng 5, no. 3: 2345-2361. https://doi.org/10.3390/eng5030122

APA Style

Pucci, M., & Zanforlin, S. (2024). Improved Lift for Thick Flatback Airfoils in the Inboard Blades of Large Wind Turbines. Eng, 5(3), 2345-2361. https://doi.org/10.3390/eng5030122

Article Metrics

Back to TopTop