12ème Congrès de Mécanique 21-24 Avril 2015 Casablanca (Maroc)

SIMULATION DES PERFORMANCES 4402 en comparant les résultats de calcul CFD avec les
mesures pour différents angles d’attaque, le travail de Yilei
AERODYNAMIQUES D’UNE He et al.([3], 2014) qui s’intéresse à l’optimisation du profil
EOLIENNE A AXE HORIZONTAL S809[4], utilisant un algorithme multi-objectifs MOGA ([5],
2006), ainsi que le travail de MEMON et al. ([6], 2012), qui
EL KHCHINE Y., SRITI M., OUARIAGLI M.
utilise le modèle de surface active pour la modélisation et la
Laboratoire des Sciences de l’Ingénieur - LSI, Université Sidi simulation des écoulements autour du profil éolien S809 [4],
Mohamed Ben Abdellah Faculté Polydisciplinaire de Taza, en comparant les résultats de simulation à ceux obtenus par
BP. 1223, Taza, Maroc des essais expérimentaux à l’aide d’une soufflerie. En outre,
Email: youneselkhchine13@gmail.com Bekhti et al. ([7], 2010) ont travaillé sur la simulation de
l’écoulement de l’air autour d’un profil NREL S809 pour
Résumé: Dans ce travail, nous étudions la simulation
prédire le phénomène de décrochage dynamique et de bien
numérique des performances aérodynamiques d’une éolienne comprendre le comportement physique de l’écoulement
à axe horizontal, qui est basée sur la résolution des équations autour d’un rotor éolien. Ces travaux sont basés sur la
de la mécanique des fluides avec un modèle de turbulence résolution des équations de Navier-Stockes moyennées en
avancé (k-ω). Cette simulation a été faite à l’aide du logiciel
utilisant le modèle RANS, pour la turbulence.
Ansys-Fluent CFD, qui permet de donner des prédictions
fines sur les profils des champs de pression, de vitesse et de Dans le présent travail, l’étude a été portée sur la
turbulence développés autour d’un élément de pale. De simulation numérique à l’aide du logiciel Ansys-Fluent pour
même, il fournit des grandeurs caractérisant les performances comprendre et analyser la structure fine de l’écoulement
aérodynamiques de l’éolienne (coefficients de portance, autour d’un profil éolien du type S809 et de déterminer aussi
trainée, moment, pression,…). Les résultats de simulation les caractéristiques aérodynamiques pour l’optimisation et
comparés avec ceux disponibles dans la littérature sont en l’amélioration de la puissance nette. En plus, la simulation
bonne concordance. Finalement, une analyse détaillée des CFD sert à calculer les efforts aérodynamiques et à dégager
résultats obtenus et des phénomènes développés autour de la des corrélations nécessaires qui seront intégrées dans un
pale a été faite afin de dégager de conclusions. modèle global pour le calcul des performances
Mots clés: Eolienne à axe horizontal, Profil S809, Simulation aérodynamiques des éoliennes équivalentes. Le choix d’un
CFD, Modèle k-ω, Performance aérodynamique, rotor éolien tripale du type S809 a été fait en raison de la
Optimisation. disponibilité des données géométriques et fonctionnelles qui
constituent une référence d’étude.
1. Introduction
2. Modélisation Mathématique du Problème
L’étude des performances aérodynamiques des éoliennes à
Dans cette étude, on s’intéresse à la simulation numérique
axe horizontal a beaucoup suscité les chercheurs pendant ces
de l’écoulement autour d’un rotor éolien à profil du type S809
dernières années à cause des applications potentielles dans la
[4] soumis exposé à la vitesse du vent sous un à angle
production de l’électricité à base de l’énergie cinétique du
d’incidence donné α (Fig.1). L’écoulement est supposé
vent qui rentre dans la catégorie des énergies renouvelables
bidimensionnel axisymétrique, instationnaire et pleinement
propres à l’environnement.
turbulent.
Cette étude est basée sur la simulation CFD (Computation
L’écoulement est gouverné donc par les équations de la
Fluid Dynamic) à l’aide du logiciel Ansys/Fluent [1] qui
conservation de masse et de la quantité de mouvement, qui
permet de résoudre les équations de la mécanique des fluides
s’écrivent respectivement:
pour des écoulements 2D et 3D, intégrant des modèles de
turbulence très avancés. Ce logiciel utilise la méthode des +∇ . =0
éléments finis, des techniques numériques appropriées et
( )
dispose d’un mailleur très adapté à des configurations +∇ . = −∇ + ∇ ̿ +
complexes.
Où t représente le temps (en s), est la masse volumique du
De nombreux travaux de recherche ont été effectués pour fluide, = ( , , ) est la vitesse d'une particule fluide, p
déterminer la structure fine de l’écoulement à travers les pales est la pression statique, ̿ est le tenseur des contraintes
des éoliennes à axe horizontal. A titre d’exemples, nous visqueuses, est la résultante des forces massiques s'exerçant
citons quelques travaux effectués dans ces dernières années, sur une particule de fluide.
tels que celui de Chandrala et al. ([2], 2012), qui porte sur
l’analyse aérodynamique d’un rotor éolien du type NACA
 12ème Congrès de Mécanique 21-24 Avril 2015 Casablanca (Maroc)

L’équation d’état, pour un gaz parfait, s’ajoute également différents angles d'incidence α= 0.0, 2.1°, 4.1°, 6.1°, 8.2°,
pour exprimer la masse volumique en fonction de la pression 10.1°, 11.2°. La pression de référence est la pression
statique. atmosphérique (Po=Patm=101 325 #$). la vitesse de rotation
Pour modéliser le phénomène de la turbulence, notre angulaire de la pale est % = 72 '/ ().
choix est porté sur le modèle K-ω SST [8], qui est connu par
Pour les conditions aux limites du domaine d’étude, nous
son efficacité pour des écoulements fortement cisaillés ainsi
avons imposé à l’entrée la vitesse du vent dont les
que des configurations complexes.
composantes = . cos(-) et = . sin(-), à la sortie la
La simulation est faite sur une section d’une corde de 1m
pression atmosphérique (P=Patm) et les conditions
du profil S809. Le domaine de calcul repose sur un
d’adhérences sur la paroi solide du profil (Vx=0, Vy=0 et Vz=
compromis qui suppose que les frontières soient
ωr), où r est la distance entre l’axe de rotation et la section de
suffisamment éloignées pour que les conditions aux limites
profil.
n’influencent pas sur l’évolution des phénomènes physiques
autour du profil. 3. Résultats et Discussions
Dans cette étude, nous avons calculé les coefficients de
portance et de traînée du profil S809 [4] par la simulation
numérique CFD, pour des conditions de fonctionnement
prédites. Les résultats de calcul sont présentés,
respectivement, sur les figures 3 et 4 illustrant
respectivement, le coefficient de portance 01 et celui de
traînée 02 .en fonction de l’angle d’incidence de l’écoulement
initial α variant de 0 à 11.2°. Ces figures montrent que les
résultats obtenus sont assez proches des résultats de Yilei He
Figure 1 : Modèle géométrique du problème [3] pour les deux coefficients, de portance 01 et de traînée02 .
Dans la figure 3 on constate que la courbe est une droite
Le maillage du domaine d’étude est effectué en adoptant
ascendante jusqu’au α=8.2°, puis subit une inflexion au delà
un raffinement des mailles au voisinage de la paroi du profil
de cette valeur. Ce point d’inflextion se traduit par un début
et desserré ailleurs. Ce maillage est choisit après une série des
de décrochage de l’écoulement au niveau de l’extrados au
tests de façon à satisfaire la stabilité des résultats obtenus. Il
bord de fuite du profil.
composé d’une grille de 18964 éléments finis quadrilatères
(Fig.2). 1,2
Calcul CFD
1
Yilei He
0,8
Cl 0,6
0,4
0,2
0
0 2,1 4,1 6,1 8,2 10,1 11,2
Angle d'attaque α

Figure3 : Distribution du coefficient de portance 01 en

fonction de l’angle d’attaque α avec 34 = 105

Figure 2: Maillage du domaine d’étude

Dans la figure 4, le coefficient de traînée 02 augmente
avec l’angle d’incidence. Ce résultat ne permet pas de se
Dans cette simulation, le fluide considéré est de l'air
prononcer sur le phénomène de décrochage, il serait donc
ambiant à une température To=288.16 K (ou 15 C°), la
important de mesurer l’évolution du rapport CL/CD en
viscosité dynamique µ=1.810-5Kg/m.s, la masse volumique fonction de l’angle d’incidence.
ρ=1.225 / ! , la vitesse relative du vent V=14.6 m/s ce
qui correspond à un nombre de Reynolds Re=106 avec les
 12ème Congrès de Mécanique 21-24 Avril 2015 Casablanca (Maroc)

0,04 Calcul
0,035 CFD Par ailleurs, ces résultats satisfaisants, vont nous permettre
0,03 Yilei He de calculer les forces de pression exercés sur le rotor qui
0,025
seront injectées ensuite sur un modèle de calcul global des
Cd 0,02
0,015 performances aérodynamiques de l’éolienne et de corriger
0,01 certaines corrélations faite dans ce modèle afin d’améliorer sa
0,005 précision et sa performance.
0
0 2,1 4,1 6,1 8,2 10,1 11,2
Angle d'attaque α References
[1] Ansys Fluent CFD, Release 13.0, 2011.
[2]M. Chandrala, A.Choubey and B. Gupta, International
Journal of Engineering Research and Applications (IJERA)
Figure 4 : distribution du coefficient de trainée en fonction de
ISSN: 2248-9622 Vol. 2, Issue6, pp.1244-1248, 2012
l’angle d’attaque α avec 34 = 105
[3] Y. He and Ramesh K. Agarwal, Shape optimization of
NREL S809 airfoil f or wind turbine blades using a
multiobjective genetic algorithm, International Journal of
Aerospace Engineering, Volume 2014, Article ID 864210,
13 pages, 2014
[4] M. JOURIEH, Développement d’un modèle représentatif
d’une éolienne afin d’étudier l’implantation de plusieurs
machines sur un parc éolien, Laboratoire de Simulation
Numérique en Mécaniques des Fluides, N˚ : 2007 ENAM
0032, 2007
[5] H. Locteau, R. Raveaux, S. Adam, Y. Lecourtier,
P.Héroux, E.Trupin, Approximation de courbes par
Figure 5 : Distribution du coefficient de pression 06 le long Algorithme Génétique Multi-Objectif (MOGA), Laboratoire
du profil pour - = 8,2° et 39 = 105 LITIS, Université de Rouen, CIFED, SDN06, pp.151-156,
2006
La figure 5 montre le coefficient de pression sur l’intrados [6] A. A. Memon, H.U. Abbasi, Modeling and Simulation of
et l’extrados le long du profil S809 avec un angle d’attaque flow around Wind Turbine Blade by using Actuator Surface
optimal de 8,2°. Tout d’abord, on remarque qu'il existe une Model, Sindh University research Journal, Vol.44 (4), pp.
grande zone de basse pression proche du bord d’attaque sur 695-698, 2012
l’extrados de profil, la différence de pression entre l’intrados [7] Ahmed Bekhti et Ouahiba Guerri, Simulation de
et l’extrados permette de créer une force de portance l’écoulement autour d’un profil de pale d’éolienne, Revue
maximale, afin d’obtenir une puissance élevée. des Energies Renouvelables, SMEE’10, pp. 49-56, 2010.
[8] F.R. Menter, Two-Equation Eddy-Viscosity Turbulence
4. Conclusion Models for Engineering Applications. AIAA Journal, AIAA
Dans ce travail, nous avons étudié l’écoulement autour Fluid Dynamics Plasma dynamics, and Lasers Conference,
d’un profil S809 pour une éolienne à axe horizontal. La N°23, Orlando Fl, USA, Vol. 32, N°8, pp. 1598-1605, 1994
simulation numérique est faite à l’aide du logiciel
Ansys/Fluent CFD, en utilisant un modèle de turbulence
: − %. Les résultats de simulation sont en bonne concordance
avec ceux de Yeili He [3], concernant les coefficients de
portance et de trainée. Les résultats de calcul ont permis de
soulever un angle d’attaque optimal α=8.2°, pour des
conditions de fonctionnement correspondant à Re=106.
Cependant, il serait important de d’étendre l’angle
d’incidence à des valeurs supérieurs pour mieux évaluer le
phénomène de décrochage de l’écoulement et d’étudier
l’influence des autres paramètres sur le comportement
aérodynamique sur ce type de profil.