29-31 Oct 2019 Nantes (France)
Wake deviation of yawed wind turbine by Large-Eddy Simulation
Félix Houtin Mongrolle  1@  , Pierre Benard  1@  , Ghislain Lartigue  1@  , Vincent Moureau  1@  , Laurent Bricteux  2@  , Julien Réveillon  1@  
1 : Complexe de recherche interprofessionnel en aerothermochimie  (CORIA)  -  Website
Centre National de la Recherche Scientifique : UMR6614, Institut national des sciences appliquées Rouen Normandie, Université de Rouen Normandie
Site Universitaire du Madrillet, BP 12, 76801 St Etienne du Rouvray Cedex -  France
2 : University of Mons [Belgium]  (UMONS)  -  Website
20, place du ParcB7000 MonsBelgium -  Belgium

According to the current energetic and environmental challenges, maximizing the electric power generated in windfarms is a societal concern. New strategies such as involving wind turbine yaw angle seem relevant to reduce wake interaction and associated power losses[1]. Therefore, yawed turbine aerodynamics is modified and remains a challenging investigation topic. 

Since experimental data on actual windfarm scales are not affordable and given the constant growth of computational resources, high order numerical simulations tend to be a promising approach[2]. The goal of this study is to evaluate a highly resolved numerical model under yaw condition in a wind tunnel before applying it to actual windfarm. The blade modeling is performed using an Actuator Line Method[3](ALM), coupled to the low Mach-number massively-parallel finite-volume Large-Eddy Simulation (LES) flow solver on unstructured meshes, called YALES2[4,5]. 

The Blind Test 5 experimental configuration led at NTNU[6], gathering numerous experimental data, is reproduced in this study. After the study of a yawed turbine wake interaction with downstream turbine the study of a single yawed turbine (+30°, 0° and -30°) will be presented. The computational domain of these cases will be the NTNU wind tunnel, involving a turbulence grid aiming to create a fully turbulent sheared [6]. The grid will be modeled using multiple Actuator Lines (to mimic the turbine blades) with dedicated polars[7,8].Each computational case is performed on a coarse and a fine mesh with around 20 and 170 millions tetrahedra, respectively.

At first, a comparison of global quantities with experimental data will be led, such as averaged power and thrust coefficient, or yaw moment. This will be followed by local characteristics analysis including velocity and turbulence intensity profiles as well as wake scan downstream of the turbine. As well, a comparison to wake deflection models will be presented. This first validation step will allow to go into a deeper analysis using advanced post processing. A study of the turbulence anisotropy[9,10] in the wake will be provided, showing the effect of yaw angle on the turbulence anisotropy fields and how tip vortices will impact the close and far wake velocity fluctuations.

 

[1] P. M. O. Gebraad et al., Wind Energ.19, 95 (2016)

[2] S-P Breton et al., Phil.Trans.R.Soc.A375(2017)

[3]J. Sørensen t al., J. Fluids Eng124(2), 393 (2002)

[4] P. Benard et al., Comput. Fluids173, 133 (2018)

[5] V. Moureau et al., C.R. Mec.339(2/3), 141, (2011)

[6] F. Mühle et al., Wind Energ. Sci.3, 883, (2018)

[7] M. Mahbub Alam et al., J. Fluid Mech., 669, 432 (2011)

[8] R.J. Martinuzzi et al., Exp. Fluids34, 585 (2013)

[9] M. Emory et al., Center for turbulence research(2014)

[10] N. Hamilton et al., Phys. Fluids,27, 015102-1 (2015)


Online user: 18