2D CFD Simulation of Dynamic Stall on a Vertical Axis Wind Turbine: Verification and Validation with PIV Measurements

Abstract

After a decrease of interest in the 1990’s, the research on Vertical Axis Wind Turbines (VAWT) a has reappeared in the last years as a result of the its increasing application in the built environment, where VAWTs present several advantages over Horizontal Axis Wind Turbines (HAWT). The VAWT has an inherent unsteady aerodynamic behavior due to the variation of angle of attack with the angle of rotation µ, perceived velocity and consequentially Reynolds number. The phenomenon of dynamic stall is then an inherent efiect of the operation of a Vertical Axis Wind Turbine (VAWT) at low tip speed ratios (‚), having a signiflcant impact in both loads and power. The complexity of the unsteady aerodynamics of the VAWT makes it extremely attractive to be analyzed using Computational Fluid Dynamics (CFD) models, where an approximation of the continuity and momentum equations of the Navier-Stokes equations set is solved. The complexity of the problem and the need for new design approaches for VAWT for the built environment has driven the authors of this work to focus the research of CFD modeling of VAWT not in the perspective of creating one large academic model to test a particular situation, but to develop a work that would: † verify the sensitivity of the model to its grid reflnement (space and time), † evaluate the difierences between the difierent commonly used turbulence models (Laminar, Spalart i Allmaras and k i †), and † be evaluated using Particle Image Velocimetry (PIV) experimental data, thus determining the suitability of this data for model validation. The 2D model created represents the middle section of a single bladed VAWT with inflnite aspect ratio. The model simulates the experimental work b of ∞ow fleld measurement using Particle Image Velocimetry (PIV) by Sim~ao Ferreira et al 1 for a single bladed VAWT, for two difierent reference Reynolds numbers of Re = 52000 and Re = 70000 for three tip speed ratios: ‚ = 2;3;4. The results show the suitability of the PIV data for the validation of the model, the unsuitability of the application of a single turbulent model and the high sensitivity of the model to grid reflnement.