1

MINISTERE DE L’ENSEIGNEMENT SUPÉRIEUR

ET DE LA RECHERCHE SCIENTIFIQUE
Universitaire de KHEMIS MILIANA
Institut des Sciences et de la Technologie
Département des sciences techniques

Simulation d’un écoulement in-stationnaire autour d’un

rotor d’une éolienne à axe vertical

Présenté par :
Bourorga Amel Proposé et dirigé par:
Ouakrine Fatmazohra
Dr. 2Noura Belkhier
 Objectif du travail

•Une modélisation numérique d’un écoulement in-stationnaire en 2D et en

3D autour d’un rotor éolien à axe vertical de type DARRIEUS à trois pales
avec un profil NACA 0021 en dimension réelles.

•Les équations utilisées sont les équations de Navier Stokes moyennées

(RANS).

•La méthode de résolution est basée sur la méthode des volumes finis, qui
permet la résolution des équations de l’écoulement.

•La rotation du rotor est modélisée par la méthode Moving Références

Frames (MRF). 3
Plan de travail

Introduction

Aérodynamique du
Rotor

Modélisation
numérique

Résultats et
discussion

Conclusion
4
 INTRODUCTION

L’énergie éolienne est une source d’énergie utilisée depuis des siècles. ce type
d’énergie a été exploité sur terre durant au moins les 3000 dernières années. Les
premières machines sont des moulins à vent à axe vertical pour moudre du grain.

moulins
moulins àà vent
vent 5
bateau
bateau àà voile
voile
 LES DIFFÉRENTS TYPES
D ’ÉOLIENNES

Généralement, on classifie les éoliennes suivant l’orientation de leur axe de rotation par
rapport à la direction du vent. On distingue ainsi deux grandes familles : Les
aérogénérateurs à axe horizontal (HAWT) et à axe vertical (VAWT).

6
Eoliennes
Eoliennes àà axe
axe horizontal
horizontal Eoliennes
Eoliennes àà axe
axe vertical
vertical
 Un aérogénérateur, ou éolienne, est un dispositif destiné à convertir l’énergie
cinétique du vent en énergie mécanique, elles sont généralement utilisés pour produire
de l’électricité et catégorie des énergies renouvelables.

Conversion
Conversion de
de l'énergie
l'énergie cinétique
cinétique du
du vent
vent
7
Aérodynamique du Rotor

8
Action
Action du
du vent
vent sur
sur un
un profil
profil de
de pale
pale

V ext
r ado
W s
Axe de référence

sens de rotation
U

α angle
d’incidence

Fonctionnement
Fonctionnement du
du
int darrieus
ra darrieus
d os Vent

le triangle de vitesse
9
pour un profil
 Théorie de Betz et coefficient de
puissance 

variation
variation du
du coefficient
coefficient de
de puissance
puissance en
en fonction
fonction de
de la
la vitesse
vitesse spécifique
spécifique
10
 l'angle d’incidence en
fonction de l’angle azimutal

L'angle d’incidence α est l'angle entre

la vitesse du vent relative et la direction
de la pale en rotation.

En raison de la variation de l'angle

d‘incidence et de la vitesse du vent
relative avec les positions d'azimutal θ ,
les forces varient et peuvent être
considérés en fonction de l'angle
d'azimutal.

sin
1
  tan
cos  
11
La puissance cinétique totale disponible sur la turbine d’une éolienne est donnée
par :

La puissance mécanique de la turbine éolienne est alors:

relation entre le couple et la puissance

12
 Les coefficients de forces peuvent être exprimées comme suit:

Coefficient de portance
L

Coefficient de trainée D

13
MODELISATION NUMERIQUE

14
 EQUATIONS RANS

la résolution numérique des équation de NAVIE -STOKES moyennée (RANS),

considéré uniquement les équation de transport pour les quantités moyennés en
modalisant toues les échèles de la turbulence.

Le modèle Spalart-Allmaras 

Le modèle Spalart-Allmaras (S-A) est un modèle relativement simple qui résout une
équation de transport modélisés pour la viscosité turbulente.

15
Domaine de calcul

16

Géométrie du rotor en 3D
Domaine de calcul

17

Géométrie du rotor en 2D
Génération du maillage 

Maillage
Maillage hybride
hybride

Maillage
Maillage non
non structure
structure
18
Maillage domaine de calcul
 Présentation des conditions
aux limites

wall
wall

Velocity-
Velocity- Pressure-
Pressure-
inlet
inlet otlet
otlet

interface
interface wall
wall

Conditions aux limites

19 en 3D
RESULTATS ET DISCUSSION

20
L’enroulement du
Sillage

21
 sillage du rotor

On observe l’enroulement alterné des nappes cisaillées de part et d’autre du rotor. Il

arrive que l’enroulement des nappes forme plusieurs tourbillons qui vont ensuite se
réunir en un seul tourbillon.

On note aussi que l’apparition d’un vortex

peut être initiée par le croisement de deux
nappes tourbillonnaires, ce qui conditionne
ainsi la position d’un tourbillon.

22 écoulement
écoulement en
en 3D
3D autour
autour du
du rotor
rotor
éolien
éolien
 Interprétation des champs de
pression de l’écoulement 2D

La figure représente la distribution de pression relative de

l’écoulement de l’air autour d’un rotor éolien. On peut voir qu’il y a
un gradient de pression au niveau de bord d’attaque : la portance est
par conséquent inexistante.

C’est donc dans cette position

que la perte d’énergie est la plus
importante une part et une forte
augmentation de la trainée
d’autre part à cause de
découlement de la couche limite.

distribution
distribution de
de Champ
Champ dede23la
la pression
pression àà un
un
angle
angle d’incidence
d’incidence 90°
90°
 Interprétation des champs de vitesse
de l’écoulement 2D

La figure présente une distribution de Champ de la vitesse autour du rotor

On remarque que des couches limites très minces sont observées le long
des parois du profil.

On observe aussi l’écoulement

commence à décoller sur
l’extrados au niveau de bord
d’attaque vers le niveau de bord
de fuite.

la distribution de Champ
24 de la vitesse
a un angle d’incidence 90°
distribution de Champ de la pression
autour du profil
déprissions
déprissions

surprissions
surprissions

Zoom sur borde d’attaque Zoom sur borde de fuite

25
 Interprétation des champs de pression
de l’écoulement 3D

La figure représente la distribution de la pression autour du profil .On observe la

pression est maximal au niveau du bord d’attaque par contre au niveau de bord de fuite
est moins importante.

26
distribution
distribution de
de Champ
Champ de
de la
la pression
pression 3D
3D
 Interprétation des champs de
vitesse de l’écoulement 3D

Cette figure présente le champs de vitesse en 3D dans un plan ZY

On observe l’enroulement
des nappes tourbillonnaire en
3D qui vont ensuite se réunir
en un seul tourbillon.

27
distribution
distribution de
de Champ
Champ de
de la
la vitesse
vitesse 3D
3D
 Courbe de coefficient de
puissance Cp

Ce Graphique donne les variations du coefficient de puissance en fonction du

TSR obtenu expérimentalement .La vitesse réduite agit directement sur le
rendement de la machine.
Coefficient de
Une vitesse réduite trop faible 0.6 puissance Cp
provoque une perte de rendement 0.5

0.4

due au décrochage dynamique, alors 0.3

coefficient de puissance
cp

qu’une vitesse réduite trop élevée 0.2

0.1
résulte aussi en une perte de 0
0.5 1 1.5 2 2.5 3 3.5
rendement à cause de la faible
incidence et de l’influence de la TS
force de traînée. Un bon équilibre R

des effets permet d’obtenir un point Variation

Variation
28
de
de coefficient
coefficient de
de puissance
puissance Cp
Cp en
en
fonction
fonction de
de TSR
TSR
 Conclusion

Dans cette étude ,nous avons modélisé un rotor éolien à axe vertical de type
DARRIEUS-H.

La modélisation numérique permet une meilleure compréhension des phénomènes

physiques à partir des champs de vitesse, de pression et des effets des lâchers
tourbillonnaires ,entrant en jeu dans le mode de fonctionnement des éoliennes à axe
vertical et plus particulièrement le décrochage dynamique.
Les résultats de la puissance obtenue par la simulation sont comparés avec celles
obtenus expérimentaux.

29
 Merci Pour
Votre
Am
el   Attention
30 Fatma- zo hra