# Numerical and Experimental Investigation on the Flow Separation Control of S809 Airfoil with Slot

Authors: Yonghui Xie; Jianhui Chen; Huancheng Qu; Gongnan Xie; Di Zhang; Mohammad Moshfeghi
Journal: Mathematical Problems in Engineering
 (2013)

---

## Abstract

A new flow control approach called split blade is applied on the S809 airfoil in the present study. S809 airfoil was investigated experimentally and numerically with different operating conditions including cascade without control, cascade with slots that generate jets with AOAs of 0 degrees, 10 degrees, 15 degrees and 20 degrees. Good agreement was obtained between the comparison of the experimental and numerical results. The results show that the separation area increases with increase of the AOA and the large separation area appears on the airfoil suction surface at AOA equal to 20 degrees. Numerical results show that the control method has little negative influence on the airfoil performance at small AOAs. Smaller vortices are filled with the large separated area which is divided by the jet generated by split when the AOA is 20 degrees. The analysis on the lift coefficient and drag coefficient shows that the flow is improved with the control. The lift coefficient and drag coefficient do not change in the comparison between the cases before and after control when AOA is 0 degrees and 10 degrees. However, the lift coefficient increases and drag coefficient decreases when AOA is 15 degrees and 20 degrees.


---

## Body

## 1. Introduction

Wind turbine is one of the mechanisms used to accomplish the transformation of wind energy into electric energy and it is the key part of wind turbine generator. The aerodynamic characteristics of wind turbine blade are directly related to the wind turbine generator’s efficiency. The research on wind turbine blade’s aerodynamic characteristics and performance has been an important subject of wind power technique. Wind turbine blade always works at high angle of attack even at deep stall condition. The flow separation phenomena will appear on the suction surface of wind turbine blade at the high angles of attack (AOAs), which affects the aerodynamic performance of wind turbine blade seriously. There is a great demand for investigating and exploiting a simple and useful control approach to suppress the flow separation on the blade surface and enhance the aerodynamic performance of blade.Flow control method is an important research field of fluid mechanics, which has been pushed more than 50 years ago. The energy decreases, the lift of airfoil is reduced, and the resistance grows by the separation. The aim of flow control is to reduce the separation zone and flow resistance.In recent years, the experimental and numerical research about wind turbine aerodynamics covers a wide variety of aspects including the control method and the dynamic wall [1, 2]. Kummer and Dang [1] embedded a crossflow fan near trailing edge into a thick wing for raising lift and the results show that the stall AOA is delayed to 40 degrees and the lift coefficient is enhanced to 6. Dygert and Dang [2] studied the same experimental rig and the same conclusion was obtained. The blowing and suction control method was studied in detail [3–5]. Mitsudharmadi and Cui [3] implemented the coflow jet concept on low Reynolds number airfoil and the results show that the momentum coefficient is found to be in the range of 5% and 13% when CFJ concept was used effectively. Chen et al. [4] used steady blowing and suction for separation and circulation control on an elliptical airfoil, which indicated that blowing is more effective than suction. Pulsed blowing was adopted for the separation control by Coiro et al. [5]. Both numerical simulation and wind tunnel experimental result show that unsteady blowing is one of the most advanced solutions for flow separation. The actuation frequency was optimized in the above study. The synthetic jet actuator (SJA) is still the new technology for wind turbines blades, which needs to be studied. Many scholars have since investigated the aerodynamics, actuation frequency, and the monument range of the SJAs. The experimental results of circular cylinder [6] show that the control effect is the best when the SJA was placed at the separated point. Xia and Luo [7] studied the flow characteristic of SJA using the PIV technology, and the optimal frequency was obtained. Farnsworth et al. [8] simulated the delta wing array SJA numerically. You and Moin [9] simulated the NC0015 airfoil with SJA. Flap is also a kind of flow control method. Zhao and Zheng [10] simulated the wind turbine with Gurney flap. Cole et al. [11] tested the various airfoils with Gurney flaps in the wind tunnel. Feng et al. [12] studied the NACA0012 airfoil with Gurney flap treatment by dielectric-barrier-discharged plasma. Lee and Su [13] investigated the control effects of the wingtip vortex experimentally. Jirásek [14] investigated the vortex generators with various shapes. Delnero et al. [15] investigated the low Reynolds number airfoil with vortex generators experimentally. Seshagiri et al. [16] studied the control effects of steady and unsteady vortex generator on airfoil. Pavel et al. [17] studied the control effects of surface HF plasma actuator on NACA 23012 airfoil. Göksel et al. [18] simulated the airfoil of ornithopter with pulsed plasma actuators. Benard et al. [19] studied the hysteresis effects for closed-loop separation control using plasma actuation.In the present study, a novel passive flow control method is studied as shown in Figure2. The separation on the blade is reduced by this method, which improves the blade aerodynamic efficiency.


## 2. Physical Model

The physical model of this research is the S809 airfoil which is used for the horizontal axis wind turbines (HAWTs). The chord length of S809 airfoil (c) is 300 mm. The CFD models consist of one outer domain and one inner domain. The inner domain is a cylinder with a radius of 300 mm and its origin is located at the leading edge (LE) of the airfoil. The outer domain extends from 600 mm (2c) upstream of the LE to 1650 mm (5.5c) downstream; the downstream length of computational domain is longer than the experimental one with the aim of avoiding the effect of outlet, and the outline of the outer domain is the contour of the experiment test section. The whole domain consists of two subdomains with a span of b=200 mm. The domains and their dimensions are shown in Figure 1.Figure 1
S809 airfoil (mm).Geometry of split airfoil.
(a)
(b)The main advantage of the model of this two-subdomain configuration is that only one set of meshes for the inner and outer domains is required. This strategy retains the same mesh arrangement in both domains at different AOAs. In particular, it removes the effects of the near wall grid sensitivity for the area near the blade surface at different AOAs.A novel passive flow control method was studied to control the flow separation on airfoil surface at different AOAs, which is shown in Figure2. The slot with the width of 3 mm was split though the body of the airfoil. One end of the cross-section of slot was placed at 30% of the chord on the pressure side of the airfoil and the other end was placed at 60% of the chord on the suction side of the airfoil. A jet is generated with the effect of pressure difference between two ends of slot, which disturbs the boundary layer flow and provides energy to the boundary layer.


## 3. Experiment Device

The experimental platform was designed and built to study the separation flow characteristics of airfoil at different AOAs and the control effect of the slot method adopted in this study. The flow configuration is shown in Figure3. The length of the test section is 1300 mm, the width of the test section is 660 mm, and the height of test section is 700 mm. In the experiment, a S809 airfoil with a chord length of 300 mm was installed in a wind tunnel, and the location of the airfoil was placed as shown in Figure 4. The center of rotation for the airfoil is designed at c/2.Figure 3
Boundary condition of simulation model.The geometry of the wing tunnel test section.
(a)
(b)The particle image velocimetry (PIV) technology is used to capture the velocity field in the present experiment. PIV is a whole-flow-field technique for capturing and measuring the instantaneous velocity vector in a cross-section of flow. The typical installation of the PIV system is shown in Figure5. The PIV system includes a charge coupled device (CCD) camera, double pulsed laser, synchronous controller, and computer. PIV measurements were performed using a CCD camera placed normal to the light sheet. The actual image resolution for the digital images at this frame rate is 1028×1028 pixels. The camera was focused onto the area of 50 mm (horizontal) × 50 mm (vertical) and 25 mm (horizontal) × 25 mm (vertical), with the top border of the image coinciding with the lower boundary of the flat plate. The experimental data using the method of averages to be processed, the limiting conditions were chosen to discarded the poor resolution results in the process of interrogation, for example, the min and max velocity, the quality of tracer particles captured by camera and so on.Figure 5
The test section and the PIV system.


## 4. Simulation Model

A commercial CFD soft CFX has been used in this study; the three-dimensional Navier-Stokes equations were solved. SST-k-ω turbulence model was used in this study. SST-k-ω model is one of the most popular RANS turbulence models in computational fluid dynamics. The SST-k-ω is a two-equation eddy-viscosity model. The shear stress transport (SST) formulation has two important features. k-ω formulation is used in the inner part of boundary layer, which solves the equations directly through the viscous sublayer, while the k-ε is used for the free stream zone [20].The equation for the turbulence kinetic energy transport is the first enclosure equation and can be presented as follows:k equation:
(1)∂(ρk)∂t+∂(ρujk)∂xj=P-β*ρωk+∂∂xj[(μ+σkμt)∂k∂xj]ω equation:
(2)∂(ρω)∂t+∂(ρujω)∂xj=γVtP-βρω211111111111111111+∂∂xj[(μ+σωμt)∂ω∂xj]11111111111111111+2(1-F1)σω2ρω∂k∂xj∂ω∂xj,P=τij∂ui∂xj,τij=μt(2Sij-23∂uk∂xkδij)-23δkδij,Sij=12(∂uj∂xi+∂ui∂xj).The turbulence eddy viscosity in the SST model is computed as(3)μt=ρa1kmax(a1ω,SF2).As free stream and inside boundary layer use the different model, all of the constants (Φ) are a blend of inner (Φ1) and outer (Φ2) constants, blended via the following linear equation:
(4)ϕ=F1ϕ1+(1-F1)ϕ2.The additional enclosure equation and constants are also given by(5)F2=tanh{[max(2kβ·ωy,500vy2ω)]2},Pk=min(τij∂ui∂xj,10β*kω),CDkω=max(2ρσω21ω∂k∂xi∂ω∂xi,10-10),F1=tanh{{min[max(kβ·ωy,500vy2ω),4σω2kCDkwy2]}},
where α1=0.05, α2=0.44, β1=3/40, β2=0.0828, β*=0.09, σk1=0.85, σk2=1.0, σω1=0.5, and σω2=0.856.The boundary condition of simulation: inlet is inlet velocity type and outlet is pressure outlet, the lateral side of the domain is wall, and up and down sides of the domain are considered to have symmetric boundary condition, which is shown in Figure3. The ambient pressure is considered to be 1 atm and the air temperature is at 25°.


## 5. Results and Discussion

The research conditions studied in the present experiment and numerical simulation are shown in Table1; the turbulence intensity of the tunnel inlet tested by hot wire anemometer was 0.08%.Table 1
The experiment and simulation condition.
Case
Inlet velocity
Re (105)
Angle of attack (°)
Case  1
30 m/s
10
0
Case  2
10
Case  3
15
Case  4
20The streamline images captured by the PIV camera are shown in Figure6. The white part is the tracer particle which shows the flow near the airfoil. The black part is the airfoil. The flow is smooth as shown by tracer particle when AOA is 0° and 10° and there is no significant separation. However, the separation zone on the airfoil surface appears to be significant which is marked out by red loop when the AOA increases to 15°. The separation point locates at the centre of the chord. The large separation area that appears on airfoil surface at AOA is 20° and the separation point is moved up to the leading edge of airfoil.The experimental results at different AOAs before control.
(a)
(b)
(c)
(d)Figure7 shows the streamlines of separation zone with control captured by CCD camera, in which the jet generated by slot nearly has no effect on the boundary layer flow of slot downstream when AOA is 0°. It mainly results from that the pressure difference between the two sides of the slot is so small that the effect of the jet is faint. The flow near the boundary layer downstream the slot is interrupted by jet when AOA is 10°. The outlet of slot is placed at the location of separation point when AOA is 15° and the jet just flows into the separation zone. The separation point moves down to the downstream of slot and the separation zone is reduced. The outlet of slot is placed at the inner region of separation zone and a jet flows into the separation zone when AOA is 20°.The experimental results at different AOAs after control.
(a)
(b)
(c)
(d)The same physical model was numerically simulated for the research on the flow structure and pressure on the airfoil surface before and after control in detail.The lifts for the sensitivity analysis at AOA is 0 degree are plotted in Table2. From Table 2, it could be seen that when the grid number was larger than 1.8 million, lift coefficient changed a little with the increase of grid number. Therefore, when grid number was up to 1.8 million, computational accuracy could be met.Table 2
The grid-independence.
Item
Grid no./10 K
Lift
Adjacent relative deviation
1
132
1.630
0.6%
2
164
1.640
0.3%
3
181
1.646
0.06%
4
210
1.647
—Figure8 is the contour of Y plus y+ on airfoil surface, and the definition of y+ is y+=Δy(τw/ρ)0.5/v, where Δy is the height of first floor grid on the blade wall, τw is wall shear stress, ρ is the fluid density, and v is kinematic viscosity coefficient. The details of y+ distribution are shown in Figure 8. Boundary layer grid on the blade surface meets the requirement for the turbulent solver.Figure 8
The Y plus distribution on the airfoil surface.Figure9 shows the experimental and CFD results without control in the case 1, case 2, case 3, and case 4 conditions. The CFD results coincide well with the experimental results. The flow is smooth on the airfoil surface when AOA is 0° while a small separation zone appears near the trailing edge when AOA is 10°. The big separation area appears on airfoil suction surface when AOA is 15° and 20°.The experimental (a) and CFD (b) results at different AOAs.
(a)
(b)Figure10 shows the lift and drag coefficient of the airfoil at different AOAs. Both lift coefficient and drag coefficient increase with AOA. The increase of the amplitude of the lift coefficient is reduced when the separation appears on the airfoil surface and the increase of the amplitude of drag coefficient gets larger significantly when AOA is 15° and 20°.Figure 10
Lift and drag coefficient on the airfoil surface.Figure11 illustrates the experimental and numerical results with slot control at different AOAs. The numerical results and experimental results are in a good agreement. It can be seen that the control method has no effect on the velocity field near airfoil when AOA is 0° and 10° by comparing the contours of velocity field for different AOAs. The effects of split in breaking down the separated zone can be observed clearly in both CFD and experiment with a high level of agreement. The separated area is predicted accurately by the SST k-ω computation especially when AOA is 15° and 20°.The experimental (a) and CFD (b) results of the velocity field of airfoil.
(a)
(b)Figure12 shows the lift and drag coefficient of the airfoil before and after control at different AOAs. The lift and drag coefficient did not change after control when AOAs are 0° and 10°. The lift coefficient increases from 61.2% to 62.8% after control when AOA is 15° while it increases from 63.9% to 68.8% when AOA is 20°. However, the drag coefficient increases slightly when AOAs are 10° and 20°.Figure 12
The lift and drag coefficient of the airfoil with and without control.Figure13 shows the reversed velocity contours at different AOAs with and without control. There is no reversed flow zone both before and after control when AOA is 0° and so the control method nearly has no influence on the boundary layer of the airfoil surface. There is a small reversed velocity area near the airfoil trailing edge when AOA is 10°. The outlet of slot is far away from the separation zone and the jet has little influence on separation zone although the reversed velocity area is also reduced after control. When AOA is 15°, the outlet of the slot is located inside the separation zone and near to separation point. When the jet flows into the separation zone, the separation thickness is reduced. The outlet of slot is inside the separation zone when AOA is 20° and the separated zone is broken down completely.Area with reversed velocity before control (a) and after control (b) at different AOAs.
(a)
(b)Figure14 shows the pressure contour before and after control at different AOAs. The pressure distribution on the pressure side is not affected by split. Therefore, the concentration is focused on the suction surface. There is no change of the pressure distribution on the suction surface after control at AOA is 0°, 10° and 15° as shown in Figure 13. The pressure distribution downstream of the outlet of the slot with control is higher than that of the slot without control and there is a high pressure area which can be interpreted as the effect of a vortex core near the outlet of slot.The pressure contour of airfoil before (a) and after control (b) at different AOAs.
(a)
(b)To illustrate the impacts of the slot on the turbulent characteristics of flow, the turbulent structure is shown in Figure15. The isosurface has a constant value of λ2=-1500, where λ2 is a kind of criteria of coherent structure which is the nonrandom composition of turbulent flow. The definition of λ2 is S2+Ω2, where S is the magnitude of the shear strain rate and Ω is the vorticity magnitude. The figures show that when AOA is 0° and 10°, the turbulent structure on the airfoil has no change after control. However, the turbulence structure generated by the separated vortex is reduced when AOA is 15° and 20°. It is noteworthy that the turbulence structure is broken into two small parts after control when AOA is 20°.The contours and isosurfaces ofλ2 before control (a) and after control (b) at different AOAs.
(a)
(b)Axial velocities in the separated area and weak zone are plotted in Figure16. Three different axial locations at x/Cax=(0.7,0.8,0.9) are also plotted, and y coordinate is the distance from the wall. It should be noted that x is the axial location and Cax is the axial chord of airfoil at different AOAs. There is no negative value of the axial velocity when AOA is 0°. The negative axial velocity appears when AOA is 10°, 15°, and 20°. The negative value of the axial velocity is higher before control than after control, which indicates that negative value of the axial velocity thickness is reduced after control.Axial velocity component at different locations at different AOAs.
(a)
(b)
(c)
(d)


## 6. Conclusion

A new passive flow control method was proposed which is opening a slot in the airfoil. The fluid with high pressure flows into the separation zone to apply the energy through the slot and the separation is controlled. The S809 airfoil was studied by experiment and numerical simulation. Different operating conditions were investigated including cases without control and cases with slots that generate jets with AOAs of 0°, 10°, 15° and 20°. The numerical results have good agreement with the experimental results. The result shows that the separation area enlarges with the increase of AOA, and the separation point moves to upstream. The big separation area appears on the airfoil suction surface when AOA is 20°. The detailed analysis of the numerical results shows that the control method has little negative influence on the airfoil performance at small AOAs. The jet generated by slot breaks the separated area into smaller vortices completely when AOA is 20°. The analysis of the lift coefficient and drag coefficient shows that the flow performance is improved with the control. The lift coefficient and drag coefficient have no change after control when AOA is 0° and 10°. However, the lift coefficient increases and drag coefficient decreases when AOA is 15°and 20°.



---

*Source: 301748-2013-12-09.xml*