Aerofólio S809 - 1 - s2 0 - s0196890412002634 - diversosperfiss8 - muitobom

Aerofólio S809 - 1 - s2 0 - s0196890412002634 - diversosperfiss8 - muitobom

(Parte 1 de 2)

Aerodynamic analysis of different wind-turbine-blade profiles using finite-volume method

Mohamed A. Sayed, Hamdy A. Kandil⇑, Ahmed Shaltot Mechatronics Engineering Department, Faculty of Engineering and Materials Science, German University in Cairo (GUC) Cairo, Egypt article i nfo

Article history: Available online 26 September 2012

Keywords: Wind energy Wind turbines Blade profiles Computational Fluid Dynamics Finite volume method abstra ct

In order to economically gain the maximum energy from the wind turbine, the performance of the blade profile must be obtained. In this paper, the results of aerodynamic simulations of the steady low-speed flow past two-dimensional S-series wind-turbine-blade profiles, developed by the National Renewable Energy Laboratory (NREL), are presented. The aerodynamic simulations were performed using a Computational Fluid Dynamics (CFD) method based on the finite-volume approach. The governing equations used in the simulations are the Reynolds-Averaged-Navier–Stokes (RANS) equations. The wind conditions during the simulations were developed from the wind speeds over different sites in Egypt. The lift and drag forces are the most important parameters in studying the wind-turbine performance. Therefore, an attempt to study the lift and drag forces on the wind turbine blades at various sections is presented. The maximum sliding ratio (lift/drag ratio) is desired in order to gain the maximum power from the wind turbine. The performance of different blade profiles at different wind speeds was investigated and the optimum blade profile for each wind speed is determined based on the maximum sliding ratio. Moreover, the optimum Angle Of Attack (AOA) for each blade profile is determined at the different wind speeds. The numerical results are benchmarked against wind tunnel measurements. The comparisons show that the CFD code used in this study can accurately predict the wind-turbine blades aerodynamic loads. 2012 Elsevier Ltd. All rights reserved.

1. Introduction

Wind power has the greatest future prospects among all types of renewable and sustainable energy technologies. Moreover, wind energy is becoming more interesting throughout the world because it is found everywhere and it is useful renewable energy while it does not affect green house due to any radiation [1].B y the 11th century, people in the Middle East were using wind mills for food production [2]. The Dutch refined the wind mills and utilized them for draining lakes and swamps. When settlers took this technology to the new world in the late 19th century they began using wind mills to pump water to farms. Later at the last decade of the 19th century, it was improved to generate electricity for homes and industry [2]. However industrialization sparked the development of larger wind mills to generate electricity which is commonly called wind turbines that appeared in Denmark as early as 1890 [3]. During 1940 the largest wind turbine operated on Vermont hill top known as Grandpa’s Knob. Afterwards, wind turbines have been improved gradually at the last decade of the twentieth century with a rapidly hugeness in their size [4,5]. At the end of 1989 a 300-kW wind turbine with a 30-meter-rotor diameter was considered very large, while by early 2004 a 4–5 MW wind turbines were commercially available [4].

The efficiency of wind energy constructions are getting important because wind energy applications have rapidly increased in the world. Moreover, wind energy is a low density source of power. To make wind power economically feasible, it is important to maximize the efficiency of converting wind energy into mechanical energy. Among the different aspects involved is rotor aerodynamics which is a key determinant for achieving this goal. There are three different approaches that may be used to analyze the flow around and downstream of a wind turbine as stated by [6] which are: field testing, which provides accurate results but is highly complex and expensive; analytical and semi-empirical models, which adopt simplifying assumptions and are thus not universally reliable; and CFD, which offers the best alternative to direct measurements.

Nowadays, industrial rotor design codes are still built based on

Blade Element Momentum theory (BEM) [6–9]. Nevertheless, in the last decade the opinion has been reached that aerodynamic modeling of Horizontal-Axis-Wind-Turbine (HAWT) rotors by means of the conventional engineering methods has reached a point where no further improvement can be expected without a

0196-8904/$ - see front matter 2012 Elsevier Ltd. All rights reserved.

⇑ Corresponding author. Address: Faculty of Engineering, Alexandria University, Alexandria, Egypt. Tel.: +20 2 2759 5539x1139, mobile: +20 100 6559603, hotline: 16482; fax: +20 2 2758 1041.

E-mail addresses: (M.A. Sayed), hamdy.kandil@- (H.A. Kandil), (A. Shaltot).

Energy Conversion and Management 64 (2012) 541–550 Contents lists available at SciVerse ScienceDirect

Energy Conversion and Management jo urn al homep ag e: w. elsev ier .co m/l oca te/e nco nma n full understanding of the flow physics [9]. Thus, the recent years have seen the rise of numerical studies on all HAWT aerodynamics features, performed on many different levels, ranging from BEM methods integrated by CFD calculations to full 3D-Navier–Stokes models.

Two-dimensional computations employing the CFD-RANS equations have been carried out at low Reynolds numbers for the wind-turbine blade profiles S809 and S826 by Sayed et al. [10]. The simulations were carried out at a wind speed of 1 m/s which is the average wind speed at Gulf of El-Zayt in Egypt. Also, 2D CFD-RANS simulations have been carried out at high Reynolds number (Re > 106) for the wind turbine blade profiles NACA 0008 and NACA0012 by Hoogedoorn et al. [1]. The unsteady separated turbulent flow around an oscillating NACA 0012 airfoil pitching in a sinusoidal pattern was investigated numerically employing the URANS approach with free stream velocity of V = 14 m/s by Wang et al. [12]. The literature shows that most researchers focus is on the range of the high wind speeds and there is no researchers focus on the low-speed range.

In order to increase the efficiency of the wind turbine, the maximum power developed by the wind turbine is desired. The power is determined in the present study by analyzing the flow around the turbine blades. The wind turbine blade profiles are selected from the profiles developed by the NREL. When comparing different airfoils as presented in this study, lift and drag coefficients for each profile are determined at different angles of attack.

Lift and drag coefficients are dimensionless numbers used to measure the aerodynamic lift and drag forces that vary with the angle of attack (a) and the shape of the airfoil. They can be defined as follows:

Therefore, the sliding ratio can be defined as, e ¼ CL

The purpose of this study is to present the aerodynamic loads of the wind-turbine blade profiles in order to find the suitable profile for the wind conditions in Egypt. In order to determine the aerodynamic loads on the selected blade profiles, ANSYS commercial software is used.

2. Method of analysis

Wind turbine blades cross sections use airfoils profiles in order to develop mechanical power. The dimensions of the blades are determined based on the desired aerodynamic performance; the maximum desired power, the airfoil properties and strength considerations. Modern HAWT blades have been consistently designed using well tested airfoil ‘families’. That is, the blade tip is designed using a thin airfoil, for high sliding ratio, and the root region is designed using a thick version of the same airfoil for structural support. Before the detailed design of wind turbine power production are explained, aerodynamic characteristics related to airfoils need to be determined [2].

The aerodynamic simulations of the steady low-speed flow past 2D S-series wind-turbine blade-profiles, developed by NREL [13] are performed using a CFD technique based on the finite-volume approach. The governing equations used in the simulations are the RANS equations. The blades performances at different wind speeds are conducted and the suitable blade profile for each wind speed will be determined based on the maximum power the blade can produce. The maximum achievable power coefficient, Cp, for turbines with an optimum blade shape but with a finite number of blades and aerodynamic drag has been calculated by Wilson et al. (1976), [2]. Their results fit the available data with accuracy within 0.5% for tip speed ratios from 4 to 20, sliding ratio from 25 to infinity and from 1 to 3 blades (B), where:


CL lift coefficient CD drag coefficient CP pressure coefficient

Cp power coefficient Re Reynolds number c airfoil chord (m) L lift force (N) D drag force (N) q air density (kg/m3)

V wind speed (m/s) A cross-sectional area of the blade (m2) B number of blades a angle of attack ( ) e sliding ratio k tip speed ratio 2D two-dimensional 3D three-dimensional

Fig. 1. Maximum achievable power coefficients of a three-bladed optimum rotor a function of the sliding ratio C /C [2].

Table 1 Selected Reynolds numbers for the present study at different wind speeds.

Wind speed(m/s) Reynolds numbers (Re)

542 M.A. Sayed et al./Energy Conversion and Management 64 (2012) 541–550

eðk þ 1 2BÞ ð4Þ

The effect of the sliding ratio on maximum achievable power coefficients for a three-bladed rotor is shown in Fig. 1. It is clear from the results on the figure that there is a significant reduction in maximum achievable power as the airfoil drag increases.

In this study, for each airfoil, the lift and drag coefficient are calculated at different mean wind speeds, as shown in Table 1.

The blade profiles used in the simulations are the S-series profiles, namely; S809, S814, S815, S817, S818, S819, S820, S821, S822, S823, S825, S826, S827, S828, S829, S830, S831, S832, S833, S834 and S835 that are depicted in Fig. 2. The ANSYS CFX [14] has been used to simulate the flow field on the selected profiles and deter-

Fig. 2. NREL selected profiles.

M.A. Sayed et al./Energy Conversion and Management 64 (2012) 541–550 543

mine the forces acting on the surfaces and their surrounding boundaries.

The computational domain size is optimized to get a domainsize independent solution and the final domain size is shown in Fig. 3. Moreover, the discretized domain is optimized based on the number of cells and the cell shape and the final mesh used in the simulations is shown in Fig. 4.

For the prediction of wind turbine aerodynamics, The Shear

Stress Transport model (SST) method has been chosen to capture the turbulence. The SST model was chosen for accurate boundary layer detection due to its ability to capture the influence of different factors that affect transition such as the free-stream turbulence and pressure gradients. For controlling and reducing the numerical solution errors, the upwind scheme method has been selected. In near wall regions, predicting the velocity gradients produced by boundary layer phenomena need elements with high aspect ratios. Inflation is applied on the airfoil surface in order to capture the effect of the boundary layer on the aerodynamic behavior. The meshing strategy selected in this study is an extruded 2D mesh in order to generate a 2D mesh of one cell thickness.

The boundary conditions applied to all regions in the domain are inlets, outlets, walls or symmetry planes. It should be noted that the flow is always subsonic due to the low range of velocity used in the study. The boundary condition at the airfoil surface has been set to no-slip-solid wall boundary. Inlet turbulence was assumed with a medium intensity (intensity = 5%). Pressure boundary conditions are applied at the domain outlets, and the average static pressure method is used in order to allow the pressure to vary locally on the boundary. Whereas for all selected airfoils under consideration, the target residual for convergence was set at 10 5.

3. Results and discussion

The simulations were carried out at different airflow velocities, as listed in Table 1. The wind speed over different sites in Egypt lies within that range which is the reason for selecting that range of velocity. The simulations are conducted for the whole range of the AOA from 5 to 15 because it is the normal operating range of the wind turbine designs. The objective of the simulation also is to find the optimum operating AOA that produces the maximum power from the wind turbine blades based on the maximum lift to drag forces. The simulations are performed for the NREL blade profiles shown in Fig. 2.

Fig. 5a–f shows the results of the sliding ratios for the selected

NREL profiles which are plotted against the AOA for the proposed velocities.

As shown, as the wind speed increases the sliding ratio increases for all profiles. At a constant wind speed, the results show that the sliding ratio first increases with increasing the AOA to reach a maximum value then decreases to almost the same low value for all profiles. The AOA at which the maximum sliding value exists changes with changing the profile.

The geometric parameters that affect the aerodynamic performance of an airfoil include: the leading edge radius, mean camber line, maximum thickness and thickness distribution of the profile and the trailing edge angle, as shown in Fig. 6. At low AOA the lift coefficient can be increased and drag can often be decreased by using a cambered airfoil [2]. So, the profiles S818, S825, S826, S830, S831 and S832 have the highest sliding ratio because all of them have the highest camber.

Moreover, the profiles S814, S815, S817, S820, S821, S822, S823,

S827, S828, S833, S834 and S835 have the intermediate sliding ratio and the profiles S809, S819 and S829 have the lowest sliding ratio because they are almost symmetric profiles which results in a slight difference in pressure between the upper surface and the lower surface of the aerofoil at zero angle of attack. In order to determine the best operating AOA range of the profiles, the

Fig. 3. Computational domain. Fig. 4. The final mesh.

544 M.A. Sayed et al./Energy Conversion and Management 64 (2012) 541–550

(a) Wind speed of 5 m/sec

(b) Wind speed of 7 m/sec (c) Wind speed of 9 m/sec

(d) Wind speed of 1 m/sec.

(e) Wind speed of 13 m/sec

(f) Wind speed of 15 m/sec

Fig. 5. The sliding ratio of the NREL profiles at different wind speeds. (a) Wind speed of 5 m/s, (b) wind speed of 7 m/s, (c) wind speed of 9 m/s (d) wind speed of 1 m/s,(e) wind speed of 13 m/s and (f) wind speed of 15 m/s.

Fig. 6. Airfoil nomenclature [2].

(Parte 1 de 2)