Airfoil Suction Side Research Articles

Linear stability analysis of turbulent mean flows based on a data-consistent Reynolds-averaged Navier–Stokes model: prediction of three-dimensional stall cells around an airfoil

Stall cells are transverse cellular patterns that often appear on the suction side of airfoils near stalling conditions. Wind-tunnel experiments on a NACA4412 airfoil at Reynolds number ${Re}=3.5 \times 10^5$ show that they appear for angles of attack larger than $\alpha = 11.5^{\circ }\ (\pm 0.5^{\circ })$ . Their onset is further investigated based on global stability analyses of turbulent mean flows computed with the Reynolds-averaged Navier–Stokes (RANS) equations. Using the classical Spalart–Allmaras turbulence model and following Plante et al. (J. Fluid Mech., vol. 908, 2021, A16), we first show that a three-dimensional stationary mode becomes unstable for a critical angle of attack $\alpha = 15.5^{\circ }$ which is much larger than in the experiments. A data-consistent RANS model is then proposed to reinvestigate the onset of these stall cells. Through an adjoint-based data-assimilation approach, several corrections in the turbulence model equation are identified to minimize the differences between assimilated and reference mean-velocity fields, the latter reference field being extracted from direct numerical simulations. Linear stability analysis around the assimilated mean flow obtained with the best correction is performed first using a perturbed eddy-viscosity approach which requires the linearization of both RANS and turbulence model equations. The three-dimensional stationary mode becomes unstable for angle $\alpha = 11^{\circ }$ which is in significantly better agreement with the experimental results. The interest of this perturbed eddy-viscosity approach is demonstrated by comparing with results of two frozen eddy-viscosity approaches that neglect the perturbation of the eddy viscosity. Both approaches predict the primary destabilization of a higher-wavenumber mode which is not experimentally observed. Uncertainties in the stability results are quantified through a sensitivity analysis of the stall cell mode's eigenvalue with respect to residual mean-flow velocity errors. The impact of the correction field on the results of stability analysis is finally assessed.

Embedded shear layers in turbulent boundary layers of a NACA0012 airfoil at high angles of attack

An investigation of turbulent boundary layers (TBLs) is presented for a NACA0012 airfoil at angles of attack 9 and 12 deg. Wall-resolved large eddy simulations (LES) are conducted for a freestream Mach number M=0.2 and chord-based Reynolds number Re=4×105, where the boundary layers are tripped near the airfoil leading edge on the suction side. For the angles of attack analyzed, mild, moderate and strong adverse pressure gradients (APGs) develop over the airfoil. Despite the strong APGs, the mean flow remains attached along the entire airfoil suction side. Similarly to other APG-TBLs investigated in the literature, a secondary peak appears in the Reynolds stress and turbulence production profiles. This secondary peak arises in the outer layer and, for strong APGs, it may overcome the first peak typically observed in the inner layer. The analysis of the turbulence production shows that other components of the production tensor become important in the outer layer besides the shear term. For moderate and strong APGs, the mean velocity profiles depict three inflexion points, the third being unstable under inviscid stability criteria. In this context, an embedded shear layer develops along the outer region of the TBL leading to the formation of two-dimensional rollers typical of a Kelvin–Helmholtz instability which are captured by a spectral proper orthogonal decomposition (SPOD) analysis. The most energetic SPOD spatial modes of the tangential velocity show that streaks form along the airfoil suction side and, as the APG becomes stronger, they grow along the spanwise and wall-normal directions, having a spatial support along the entire boundary layer.

Turbulence Modeling of Iced Wind Turbine Airfoils

Icing is a severe problem faced by wind turbines operating in cold climates. It is affected by various fluctuating parameters. Due to ice accretion, a significant drop in the aerodynamic performance of the blades’ airfoils leads to productivity loss in wind turbines. When ice accretes on airfoils, it leads to a geometry deformation that seriously increases turbulence, particularly on the airfoil suction side at high angles of attack. Modeling and simulation are indispensable tools to estimate the effect of icing on the operation of wind turbines and gain a better understanding of the phenomenon. This paper presents a numerical study to assess the effect of surface roughness distribution, along with the effect of two turbulence models on estimating wind turbine airfoils’ aerodynamic performance losses in the presence of ice. Aerodynamic parameter estimation was performed using ANSYS FLUENT, while ice accretion was simulated using ANSYS FENSAP-ICE. The results using the adopted modeling approaches and the simulation tools were compared with another numerical study and validated against experimental data. The validation process demonstrated the model’s accuracy when considering roughness distribution via the beading model available in ANSYS FENSAP-ICE. The two turbulence models examined (Spalart–Allmaras and k-ω SST) gave comparable results except for the drag at high angles of attack. The k-ω SST model was more efficient in replicating turbulence at high angles of attack, leading to higher accuracy in aerodynamic loss estimation.

Direct numerical simulation of transitional flow past an airfoil with partially covered wavy roughness elements

We perform direct numerical simulation of flow past the NACA (National Advisory Committee for Aeronautics) 0012 airfoil with partially covered wavy roughness elements near the leading edge. The Reynolds number Rec based on the freestream velocity (U0) and airfoil chord length (C), and the angle of attack (AoA) is fixed, i.e., Rec=5×104, AoA=10°. The leading edge roughness elements are characterized by streamwise sinusoidal-wave geometry that is uniform in the spanwise direction with fixed height, whereas varying wave numbers (k) from 0 to 12. Based on the validation of the smooth baseline case (k = 0), the roughness effects on the aerodynamic performance are evaluated in terms of the lift and drag coefficients. The drag coefficient decreases monotonically with k, while the variation of the lift coefficient with k is similar to the “drag crisis” phenomenon observed in cylinder flows. The sharp variations of lift and drag coefficients from k = 6 to 8 indicate that k = 8 is a critical case, beyond which massive separation occurs and almost covers the airfoil's suction side and dominates the airfoil aerodynamic performance. To further reveal the underlying mechanism, the flow structures, pressure, skin friction coefficients, shear layer transition onset, and boundary layer separation are quantified to investigate the roughness effects. The roughness elements strongly modify the separation behavior, whereas they have little effect on the transition onset. The unsteady interactions and convections of separation bubbles downward the trailing edge also change the wake evolution. Based on the scaling of the asymmetric wake profiles, we find that the wake defect is gently decreasing with k, but the increase in the wake width is much stronger. This scaling analysis gains a better insight into the roughness effects on the momentum deficit flow rate, which confirms the drag mechanism with different roughness wave numbers.

Effect of Different Types of Erosion on the Aerodynamic Performance of Wind Turbine Airfoils

Taking the S823 airfoil as the research object, this study investigates the influence of different types of leading-edge erosion on the aerodynamic performance of airfoil by using the computational fluid dynamics method. The effect of leading-edge erosion on the inception of stall vortex is also analysed. The results show that when the angle of attack (AoA) is greater than 5°, the leading-edge erosion results in a significant decrease in the lift coefficient and an increase in the drag coefficient. The deterioration in the drag coefficient of the airfoil caused by leading-edge erosion is much greater than that of the lift coefficient. Moreover, the maximum promotion rate of the drag coefficient can reach 357% at Re = 300,000. The exacerbation of the erosion level leads to a dramatic expansion of the stall vortex on the airfoil suction side at a large AoA and results in a reduction in the pressure difference between the pressure and suction sides of the airfoil. This is also the reason erosion causes the degradation of the aerodynamic performance of the wind turbine airfoil. This work is beneficial to establish the reasonable maintenance cycle of the wind turbine blades working in a sand blown environment.

Switch of tonal noise generation mechanisms in airfoil transitional flows

Large eddy simulations are performed to study tonal noise generation by a NACA0012 airfoil at an angle of attack ${\alpha = 3}$ deg. and freestream Mach number of ${M_{\infty} = 0.3}$. Different Reynolds numbers are analyzed spanning ${0.5 \times 10^5 \le Re \le 4 \times 10^5}$. Results show that the flow patterns responsible for noise generation appear from different laminar separation bubbles, including one observed over the airfoil suction side and another near the trailing edge, on the pressure side. For lower Reynolds numbers, intermittent vortex dynamics on the suction side results in either coherent structures or turbulent packets advected towards the trailing edge. Such flow dynamics also affects the separation bubbles on the pressure side, which become intermittent. Despite the irregular occurrence of laminar-turbulent transition, the noise spectrum depicts a main tone with multiple equidistant secondary tones. Increasing the Reynolds number leads to a permanent turbulent regime on the suction side that reduces the coherence level causing only small scale turbulent eddies to be observed. Furthermore, the laminar separation bubble on the suction side almost vanishes while that on the pressure side becomes more pronounced and permanent. As a consequence, the dominant noise generation mechanism becomes the vortex shedding along the wake.

Characterization Of Shock Buffet On A Supercritical 2D Airfoil In Transonic Flow

The transonic flow over the supercritical OAT15A airfoil is investigated experimentally. Using time-resolved focusing schlieren sequences, we extracted the temporal evolution of the chordwise shock location for a wide range of combinations of Mach number and angle of attack. Based on this data, we identified the conditions of pre-buffet, buffet onset, fully-established buffet, and buffet-offset behavior. Key parameters characterizing most pronounced buffet behavior are presented and quantified in terms of the shock amplitude, its variation over multiple buffet cycles, and the associated frequency of the periodic shock oscillation. Fully-established buffet is identified at a Mach number of 0.72 and angle of attack of 5 deg. It is demonstrated that these operating conditions are characterized by a quasi-sinusoidal shock oscillation at a distinct frequency of 113 Hz, in good agreement with previous studies. This analysis reveals that the harmonic shock oscillation spanning a range of about 16.5% of the airfoil chord is accompanied by a coupled large-scale variation of the global flow topology. Both instantaneous and mean velocity components in the chordwise and vertical directions are discussed for representative phases of the buffet cycle. Instantaneous vector fields of both components provide insight into the flow organization throughout the buffet cycle and allow to identify the streamwise distance traveled by the shock wave within a full cycle. We can thus precisely capture the state of the shock wave and quantify the evolution of the shock-induced separation and the size of the turbulent wake. Severe variations of the flow topology are induced by the large-scale upstream and downstream shock displacement. The results show that incidents of massive flow separation are observed for shock positions being close to the upstream limit, and more general, that the wake influence is more pronounced when induced by an upstream traveling shock wave. Phase-averaged velocity profiles at significant locations along the airfoil suction side complement this analysis. Characteristic sets of velocity profiles are discussed in detail that allow to uniquely describe the development of the flow in each phase of the buffet cycle. Detailed contours and profiles of both velocity and Reynolds shear stress contributions across the wake region highlight the strong downstream influence the shock-induced turbulent wake. Phases of severe flow separation result in reverse flow persisting until the trailing edge which further substantiates the induced momentum deficit.

Adverse-pressure-gradient turbulent boundary layer on convex wall

Direct numerical simulations (DNSs) of an incompressible turbulent boundary layer on an airfoil (suction side) and that on a flat plate are compared to characterize the non-equilibrium turbulence and the effect of wall curvature on the flow. The two simulations effectively impose matching streamwise distributions of adverse pressure gradient (APG) quantified by the acceleration parameter (K). For the airfoil flow, an existing compressible DNS carried out by Wu et al. [“Effects of pressure gradient on the evolution of velocity-gradient tensor invariant dynamics on a controlled-diffusion aerofoil at Rec = 150,000,” J. Fluid Mech. 868, 584–610 (2019)] of the flow around a controlled-diffusion airfoil is used. For the flat-plate flow, a separate simulation is carried out with the aim to reproduce the flow in the region of the airfoil boundary layer with zero to adverse pressure gradients. Comparison between the two cases extracts the effect of a mild convex wall curvature on velocity and wall-pressure statistics in the presence of APG. In the majority part of the boundary layer development, curvature effect on the flow is masked by that of the APG, except for the region with weak pressure gradients or a thick boundary layer where the effect of wall curvature appears to interact with that of APG. High-frequency wall-pressure fluctuations are also augmented by the wall curvature. Overall, the boundary layers are qualitatively similar with and without the wall curvature. This indicates that a flat-plate boundary layer DNS may serve as a low-cost surrogate of a boundary layer over the airfoil or other objects with mild curvatures to capture important flow features to aid modeling efforts.

Power augmentation of Darrieus wind turbine blades using trapped vortex cavity

Enhancing the performance of Darrieus Vertical Axis Wind Turbines (VAWTs) is the key to boost their possible commercialisation in the yet highly competitive wind energy market. This study evaluates an active boundary layer control technique to increase the aerodynamic efficiency of Darrieus VAWTs. A cavity is created on the airfoil's suction side, where suction is applied to create a trapped vortex. High-fidelity Computational Fluid Dynamics (CFD) simulations are used after validation against experimental results. The effect of suction momentum ratio Cμ on the performance of the turbine is investigated. The turbine with suction cavity airfoils (SC) shows a higher power coefficient than the one with baseline airfoil (BL) blades, especially at low tip-speed ratios. For the selected test case, the net maximum predicted power coefficient is 0.435 for the SC turbine at TSR = 2.3, 0.375 for the baseline at TSR = 2.6. At the lower tip-speed ratios, i.e. TSR = 2, the power coefficient of the SC turbine is almost double that of the BL turbine. The feasibility of the active control technique is finally discussed, demonstrating its potential as an effective solution for the performance augmentation of the Darrieus turbine.

Investigation of Blowing and Suction for Turbulent Flow Control on Airfoils

An extensive parametric study of turbulent boundary-layer control on airfoils via uniform blowing or suction is presented. The control is applied on either the suction or pressure side of several four-digit NACA-series airfoils. The considered parameter variations include angle of attack, Reynolds number, control intensity, airfoil camber, and airfoil thickness. Two comprehensive metrics, designed to account for the additional energy required by the control, are introduced to evaluate the net aerodynamic performance enhancements. The study confirms previous findings for suction-side boundary-layer control and demonstrates the interesting potential of blowing on the pressure side under various conditions, which achieves a maximum total net drag saving of 14% within the considered parameter space. The broad parameter space covered by the presented Reynolds-average Navier–Stokes simulations allows for more general conclusions than previous studies and can thus provide guidelines for the design of future detailed experimental or numerical studies on similar boundary-layer control schemes.

Aerodynamic investigation of a linear cascade with tip gap using large-eddy simulation

The flow in a linear compressor cascade with tip gap is simulated using a wall-resolved compressible Large-Eddy Simulation. The cascade is based on the Virginia Tech Low Speed Cascade Wind Tunnel. The Reynolds number based on the chord is 3.88 x 10⁵ and the Mach number is 0.07. The gap considered in this study is 4.0 mm (2.9% of axial chord). An aerodynamic analysis of the tip-leakage flow allow us identifying the main mechanisms responsible for the development and the convection of the tip-leakage vortex downstream of the cascade. A region of high turbulence and vorticity levels is located along an ellipse that borders the top of the tip-leakage vortex. The influence of the airfoil suction side boundary layer development on the tip-leakage vortex is highlighted by tripping the flow. A tripped boundary layer induces a stronger and larger tip-leakage vortex that tends to move further away from the airfoil suction side and from the endwall compared with an untripped flow. The boundary layer turbulent state influences the tip-leakage flow development.

A calculation method for modeling the flow characteristics of the wind turbine airfoil with leading-edge protuberances

An effective flow separation control method is urgently needed to improve the aerodynamic performance of the wind turbine blades. The leading-edge protuberances (LEPs) are the method that can inhibit flow separation and improve the aerodynamic performance of wind turbine airfoil. In order to achieve the optimal design of the LEPs on the wind turbine airfoil, it is necessary to develop a theoretical model of the effect of the LEPs on the flow characteristics. This paper proposes a novel modeling calculation method as achieved by substituting the momentum source into the N–S equation, and establishes the relationship between the geometrical parameters of the LEPs and the flow characteristics. The results show that the aerodynamic and flow characteristics of the modeled and solid airfoil with LEPs maintain good consistency under different stall conditions. Moreover, the high-strength streamwise vortices induced by LEPs increase the area of low pressure in the airfoil suction side, which effectively suppresses the flow separation and improves the aerodynamic performance of the airfoil. The modeling method provides the basis for the optimal design of the wind turbine airfoil with LEPs. The process of model construction explains the flow control mechanism of the LEPs from a new perspective.

Loss Reduction in a 1.5 Stage Axial Turbine by Computer-Driven Stator Hub Contouring

Improvements in stage isentropic efficiency and reductions in total pressure loss are sought in a 1.5 stage axial turbine. This is representative of power generation equipment used in thermal power cycles, which delivers about 80% of the 20 × 1012 kWh world-wide electricity. Component-level improvements are therefore timely and important toward achieving carbon dioxide global emission targets. Secondary flow loss reduction is sought by applying a nonaxisymmetric endwall design to the turbine stator hub. A guide groove directs the pressure side branch of the horseshoe vortex away from the airfoil suction side, using a parametric endwall hub surface, which is defined as to obtain first-order smooth boundary connections to the remainder of the passage geometry. This delays the onset of the passage vortex and reduces its associated loss. The Automatic Process and Optimization Workbench (apow) generates a Kriging surrogate model from a set of Reynolds-averaged Navier–Stokes simulations, which is used to optimize the hub surface. The three-dimensional steady Reynolds-averaged Navier–Stokes model with an axisymmetric hub is validated against reference experimental measurements from the Rheinisch-Westfälische Technische Hochschule (RWTH) Aachen. Comparative computational fluid dynamics (CFD) predictions with an optimized nonaxisymmetric hub show a decrease in the total pressure loss coefficient and an increase in the isentropic stage efficiency at and off design conditions.

On the boundary layer structure on suction side of an airfoil

The structure of boundary layer on an airfoil suction side was studied using stereo time resolved PIV technique experimentally in low Reynolds number, about 33 thousands. The mean velocity field is close to the 2D case, however the instantaneous structure is highly dynamical and 3D. Dynamics of the vortices in the boundary layer has been analysed using POD method. The results show lag of coherence in the streamwise direction, oblique patterns are detected instead.

Unsteady boundary layer development on a wind turbine blade: an experimental study of a surrogate problem

Wind turbines with thick blade profiles experience turbulent, periodic approach flow, leading to unsteady blade loading and large torque fluctuations on the turbine drive shaft. Presented here is an experimental study of a surrogate problem representing some key aspects of the wind turbine unsteady fluid mechanics. This experiment has been designed through joint consideration by experiment and computation, with the ultimate goal of numerical model development for aerodynamics in unsteady and turbulent flows. A cylinder at diameter Reynolds number of 65,000 and Strouhal number of 0.184 is placed 10.67 diameters upstream of a NACA 63215b airfoil with chord Reynolds number of 170,000 and chord-reduced frequency of $$k=2\pi f\frac{c}{2}/V=1.5$$ . Extensive flow field measurements using particle image velocimetry provide a number of insights about this flow, as well as data for model validation and development. Velocity contours on the airfoil suction side in the presence of the upstream cylinder indicate a redistribution of turbulent normal stresses from transverse to streamwise, consistent with rapid distortion theory predictions. A study of the boundary layer over the suction side of the airfoil reveals very low Reynolds number turbulent mean streamwise velocity profiles. The dominance of the high amplitude large eddy passages results in a phase lag in streamwise velocity as a function of distance from the wall. The results and accompanying description provide a new test case incorporating moderate-reduced frequency inflow for computational model validation and development.

Large Eddy Simulations on a pitching airfoil: Analysis of the reduced frequency influence

Large Eddy Simulations (LES) have been performed on a pitching NACA0012 airfoil at Rec=2.104. The influence of the reduced frequency on the dynamic stall phenomenon is investigated. Special attention is given on the aerodynamic coefficients drop incidence delay in comparison with a static case. The analysis is based on vorticity field observations and the introduction of local lift and friction coefficients. The boundary layer events and the presence and the impact of a Leading-Edge Vortex (LEV) are investigated. The total stall delay is decomposed into two main phases: a delay in the boundary layer separation which occurs at higher incidence as the reduced frequency increases and the initiation and the growth of the LEV. For selected reduced frequencies, the influence of the LEV on the lift and drag coefficients is analysed as a function of the airfoil angle of incidence and then as a function of time. The LEV maintains a high level of lift and drag. Its lifetime on the airfoil suction side seems to be shortened by the beginning of the airfoil downstroke motion. Given the fact that the LEV appears at higher incidence as the reduced frequency grows, the LEV lifetime on the airfoil suction side decreases as the reduced frequency increases. However, the airfoil oscillation being faster, the LEV goes through larger incidence ranges during its lifetime and leads to a larger lift and drag drop incidence delay as the reduced frequency increases.

Numerical investigation of transition in a boundary layer subjected to favourable and adverse streamwise pressure gradients and elevated free stream turbulence

Laminar-to-turbulent transition of a boundary layer subjected to streamwise pressure gradients and elevated free stream turbulence is computed through direct numerical simulation. The streamwise pressure distribution and elevated free stream turbulence levels mimic the conditions present on the suction side of highly-cambered airfoils. Longitudinal streamwise streaks form in the laminar boundary layer through the selective inclusion of low-frequency disturbances from the free stream turbulence. The spanwise spacing normalized by local inner variables indicates stabilization of the streaks occurs by the favourable pressure gradient and prevents the development of secondary streak instability modes until downstream of the suction peak. Two distinct processes are found to trigger transition to turbulence in the adverse pressure gradient region of the flow. One involves the development of varicose secondary instability of individual low-speed streaks that results in their breakdown and the formation and growth of discrete turbulent spots. The other involves a rapid amplification of free stream disturbances in the inflectional boundary layer in the adverse pressure gradient region that results in a largely homogeneous breakdown to turbulence across the span. The effect of high-frequency free stream disturbances on the streak secondary instability and on the nonlinear processes within the growing turbulent spot are analysed through the inviscid transport of instantaneous vorticity. The results suggest that free stream turbulence contributes to the growth of the turbulent spot by generating large strain rates that activate vortex-stretching and tilting processes within the spot.

Development of a Fluidic Actuator for Adaptive Flow Control on a Thick Wind Turbine Airfoil

Wind turbines are exposed to unsteady incident flow conditions such as gusts or tower interference. These cause a change in the blades' local angle of attack, which often leads to flow separation at the inner rotor sections. Recirculation areas and dynamic stall may occur, which lead to an uneven load distribution along the blade. In this work, a fluidic actuator is developed that reduces flow separation. The functional principle is adapted from a fluidic amplifier. High pressure air fed by an external supply flows into the interaction region of the actuator. Two control ports, oriented perpendicular to the inlet, allow for a steering of the actuation flow. One of the control ports is connected to the suction side, the other to the pressure side of the airfoil. Depending on the pressure difference that varies with the angle of attack, the actuation air is directed into one of four outlet channels. These guide the air to different chordwise exit locations on the airfoil's suction side. The appropriate actuation location adjusts automatically according to the pressure difference between the control ports and therefore incidence. Suction side flow separation is delayed as the boundary layer is enriched with kinetic energy. Experiments were conducted on a DU97-W-300 airfoil at Re = 2.2 × 105. Compared to the baseline, lift variations due to varying angles of attack were reduced by an order of magnitude. A Fast/Aerodyn simulation of a full wind turbine rotor was performed to show the real world load reduction potential. Additionally, system integration is discussed, which includes suggestions on producibility and operational details.

Comparative research of algorithms of measuring the geometry of complex profiles of gas turbine engine compressor blades

The paper presents three algorithms for the measurement of gasturbine engine (GTE) compressor blades airfoil (suction side, pressure side, leading and trailing edges) geometry using coordinate measuring machines. The first algorithm is a single scan of the heightwise blade section using a test point. The second algorithm is based on the pre-alignment of the nominal profile with the measured actual profile. An iterative closest-point algorithm is used to carry out the procedure of the best possible alignment of the measured ptofile and the nominal one. In the third algorithm, the coordinates of the measured points of the edges are calculated using the spline defined by the coordinates of the test point (TP) center. The research of errors of measuring the geometry of the elements using the algorithms mentioned above was carried out. To perform the research a model was developed simulating the contact of the TP and the surface of the part being measured. The model makes it possible to calculate the coordinates of the contact point, the center coordinates and the coordinates of the MT of the measured points. The search for the coordinates of the contact was performed using the method of sequential quadratic programming. A procedure was also used that makes it possible to simulate the deviation of the profile and surface position from the nominal. In the course of the simulation dependences for determining the average, lower and upper boundaries of the measurement error distribution for any airfoil point being measured were obtained.

An experimental and numerical investigation on the formation of stall-cells on airfoils

Stall Cells (SCs) are large scale three-dimensional structures of separated flow that have been observed on the suction side of airfoils designed for or used on wind turbine blades. SCs are unstable in nature but can be stabilised by means of a localized disturbance; here in the form of a zigzag tape covering 10% of the wing span. Based on extensive tuft flow visualisations, the resulting flow was found macroscopically similar to the undisturbed flow. Next a combined investigation was carried out including pressure recordings, Stereo-PIV measurements and CFD simulations. The investigation parameters were the aspect ratio, the angle of attack and the Re number. Tuft and pressure data were found in good agreement. The 3D CFD simulations reproduced the structure of the SCs in qualitative agreement with the experimental data but had a delay of ~3deg in capturing the first appearance of a SC. The error in Cl max prediction was 7% compared to 19% for the 2D cases. Tests show that SCs grow with Re number and angle of attack. Also analysis of the time averaged computational results indicated the presence of three types of vortices: (a) the trailing edge line vortex (TELV) in the wake, (b) the separation line vortex (SLV) over the wing and (c) the SC vortices. The TELV and SLV run parallel to the trailing edge and are of opposite sign, while the SC vortices start normal to the wing suction surface, then bend towards the SC centre and later extend downstream, with their vorticity parallel to the free stream.