Modèle numérique de production annuelle d’une éolienne à axe horizontal

Description des éoliennes à axe horizontal

La production annuelle d’une éolienne repose sur la courbe de puissance de l’éolienne, les vitesses de vent et le nombre d’heures de vent durant une année pour un site proposé. Les détails du calcul de la production annuelle d’une éolienne sont présentés dans la section 2.2 du chapitre 2, ici nous identifions les éléments les plus importants qui influencent sa valeur et leur lien avec le phénomène de givrage. La puissance d’une éolienne est déterminée à partir de son coefficient ~e puissance, de la surface balayée par ses pales, de la densité de l’air et de la vitesse du vent au cube. La puissance dépend principalement de la vitesse du vent, mais est aussi légèrement affectée par la densité de l’air, cette dernière étant déterminée à son tour par la température et la pression barométrique ou par l’altitude [7]. Le coefficient de puissance dépend des caractéristiques des profils aérodynamiques composants la pale, de la vitesse du vent, du nombre de pales, de leur rayon, de leur angle de calage et de leur vitesse de rotation ([8] et [9]). Les vitesses de vent sont représentées par une fonction de distribution qui représente la probabilité d’avoir une certaine vitesse du vent pendant une période donnée. Le vent ayant un caractère aléatoire, son analyse est basée sur des méthodes statistiques. Le régime des vents sur un site se caractérise par la valeur de la vitesse moyenne du vent et les lois statistiques permettant de représenter sa variation temporelle [10]. Parmi ces lois, la loi de Weibull est souvent utilisée en génie éolien, car elle représente bien les vitesses du vent observées à long terme sur un grand nombre de sites [5] et prend en compte l’asymétrie des données empiriques [11].

Un cas particulier de la distribution de Weibull est la distribution de Rayleigh. C’est la distribution de la vitesse de vent la plus simple pour représenter les ressources éoliennes, car elle exige seulement de connaître la vitesse moyenne du vent, l’asymétrie étant fixée à deux. Le calcul de la production annuelle est basé sur une vitesse moyenne du vent à la hauteur du moyeu pour le site étudié. En raison du phénomène de cisaillement du vent, la vitesse du vent à la hauteur du moyeu est habituellement plus élevée que celle mesurée à la hauteur habituelle où sont installés les anémomètres. Pour cette raison une extrapolation verticale des vitesses du vent mesurées est effectuée, au moyen d’une loi de puissance appliquée à la hauteur [12]. Si les lois de Weibull et de Rayleigh permettent de déterminer le comportement statistique de la vitesse du vent dans le temps, la loi de puissance permet de mesurer la variabilité du vent au même endroit, au même moment mais à des hauteurs différentes. Le calcul de la production annuelle, basé sur l’intégration du produit entre la puissance produite par l’éolienne à une vitesse de vent donnée et le nombre d’heures durant l’année pendant lesquelles le vent souffle à cette vitesse est valable pour des conditions standards de température et de pression atmosphérique. La production doit être corrigée en fonction de la pression atmosphérique, des conditions de température et des pertes inhérentes aux éoliennes telles que les pertes par effet de sillage, les pertes par encrassement des pales ou résultantes du verglas, les pertes pour causes d’arrêts et les autres pertes [5].

Exploitation des éoliennes en milieu nordique En milieu nordique, les basses températures augmentent la densité de l’air. Une densité d’air plus élevée augmente l’énergie cinétique du vent et par conséquent la puissance captée par l’éolienne. L’air est par exemple 27 % plus dense à -33°C qu’à +35°C. La puissance fournie par une éolienne, qui est directement proportionnelle à la densité de l’air, est donc 27 % plus importante à -33°C qu’à +35°C [13]. Les milieux nordiques disposent généralement de régions très bien exposées au vent et ces régions semblent idéales pour l’exploitation des éoliennes. Malheureusement, les problèmes liés au givrage et aux très basses températures y sont fréquents et ont de sérieuses conséquences sur la production des éoliennes, leur maintenance et leur durée de vie [14]. La présence de glace sur les pales d’éoliennes diminue considérablement les performances du rotor. Lors de conditions d’évènements givrants, un dépôt de glace se forme sur les pales (Figure 2), celles-ci perdent leurs caractéristiques aérodynamiques [15], le couple rotatif s’affaiblit et la puissance électrique fournie par l’éolienne diminue. Lorsque l’évènement givrant est très sévère [13], la production électrique est interrompue en raison d’un couple rotatif nul qui conduit à l’arrêt de l’éolienne [14].

Le Tableau 1 présente un ordre de grandeur des valeurs moyennes annuelles des pertes de production dues au givrage pour un parc éolien [16]. Les pertes de production dépendent à la fois de l’intensité, de la longueur et de la fréquence des évènements givrants. Lors d’un évènement givrant assez important, la masse de glace accumulée peut augmenter considérablement le chargement radial sur les pales en raIson de la force centrifuge. Le système de fixation des pales au moyeu doit alors être spécialement dimensionné pour supporter le surplus de contraintes et éviter tout bris mécanique. Les éoliennes peuvent aussi être arrêtées en raison de vibrations trop importantes qUI apparaissent lorsque le dépôt de glace n’est pas réparti uniformément [18]. Ces vibrations favorisent l’expulsion dangereuse de morceaux de glace [19]. La présence de la glace sur les appareils de mesure d’une éolienne altère également le fonctionnement de ces appareils. Un anémomètre givré [15] ou une girouette givrée ne sont pas en mesure de fournir des données cohérentes, le contrôle de l’éolienne devient alors difficile ou même impossible. Sur certains sites éoliens à climat très froid (régions nordiques et/ou sites montagneux), les températures peuvent descendre en dessous des seuils de fonctionnement de certaines composantes, notamment électroniques, de l’éolienne.

Pour éviter tout bris mécanique, l’huile des boites de vitesses est chauffée régulièrement. Enfin, à très basse température et pour des vitesses de vent favorable, une éolienne peut paradoxalement se retrouver en situation de surproduction à cause de la densité de l’air trop élevée. La génératrice doit alors être dimensionnée pour supporter ce surplus de puissance sans être endommagée [15]. Certaines adaptations sont obligatoires pour envisager d’exploiter des éoliennes en milieu nordique. Ces adaptations concernent à la fois la conception mécanique de l’éolienne (matériaux adaptés aux faibles températures, dimensionnement du moyeu adapté aux contraintes radiales supplémentaires dues à l’accumulation de glace sur les pales), le choix des équipements de l’éolienne (appareils de mesures météorologiques chauffés, boites électriques chauffées, mise en place de systèmes de dégivrage automatisés), et enfin l’entretien et la maintenance de l’éolienne (dégivrage manuel des pales, remplacement des pièces brisées ou usées par le froid et le givre, chauffage de l’huile des boîtes de vitesses).

Modèles numériques de givrage

Plusieurs modèles numériques ont été développés pour simuler la formation de glace (le givrage) sur un profil aérodynamique en deux dimensions afin de prédire le type de glace, la forme et la location du dépôt, ainsi que leur effet sur l’aérodynamique des avions et des hélicoptères en fonction des conditions atmosphériques mesurées. Ils utilisent comme conditions d’entrée la température et l’ humidité relative de l’air, la pression barométrique, la vitesse de vent, la teneur en eau liquide de l’air (LWC), le diamètre volumétrique médian (MVD) des gouttelettes d’eau en surfusion et la durée de l’évènement, et font des calculs successifs de l’écoulement autour du profil et des trajectoires des gouttelettes d’eau pour prédire l’efficacité de collection, la couche limite et les coefficients de transfert, du bilan de masse et transfert de chaleur à la surface de l’objet pour prédire la masse de glace accumulée, et des changements de géométrie pour prédire la forme de dépôt de glace. Les codes validés pour l’aéronautique, en deux dimensions, sont LEWICE développé aux États- Unis par Wright [23] en 1995, ONERA développé en France par Gent [24] en 1990, TRAJICE2D développé en Grande-Bretagne par Guffond [25] en 1992, CANICE développé au Canada (École polytechnique) par Paraschivoiu [26] en 1994, CIRA développé en Italie par Mingione [27] en 1996, CIRALIMA développé au Québec par G. Fortin [20] en 2003. Le code FENSAPICE 3D, développé au Canada (Université McGill) par Habashi [28] en 2002, est le seul qui peut prédire la formation de glace sur un avion en 3D. Ces différents codes sont décrits et comparés plus en détail dans la thèse de doctorat de G. Fortin [20].

L’Institut technique de la recherche en Finlande (Technical Research Institute of Finland VTT) a développé le code TURBICE pour la modélisation de l’accrétion de glace sur les pales d’éolienne. La version initiale a été établie en 1992 par Finstand et Makkonnen [29] et a été amélioré par la suite [30]. Le code CANICE-WT a été développé à l’École polytechnique de Montréal (Canada) pour la simulation et l’analyse de la formation de glace (givre et verglas) sur \es éoliennes à axe horizontal. Le code est capable de prédire la masse et la forme de glace se déposant le long d’une pale d’éolienne connaissant la géométrie de l’éolienne, les conditions d’opération et les conditions atmosphériques [22]. Étant donné que l’accès du code CIRALIMA était facile et que ses résultats sont fiables et précis comparés aux autres logiciels, ce code est utilisé dans ce travail. CIRALIMA permet de prédire en 20 la forme des profils aérodynamiques après l’accumulation de glace, en régime sec (givre) et humide (verglas), et ce pour différentes conditions atmosphériques et différentes positions sur la pale.

Table des matières

RÉSUMÉ
ABSTRACT
REMERCIEMENTS ET AVANT PROPOS
T ABLE DES MA TIÈRES
LISTE DES FIGURES XVI
LISTE DES TABLEAUX XX
LISTE DES SYMBOLES ET ABRÉVIATIONS
CHAPITRE 1 INTRODUCTION
1. 1. C ONTEXTE DU PROJET
1.1.1. Description des éoliennes à axe horizontal
1.1.2. Production annuelle d ‘une éolienne
1.1.3. Exploitation des éoliennes en milieu nordique
1.2. PROBLÉMATIQUE
1.3. O BJECTIFS
1.3. 1. Objectif principal
1.3.2. Objectifs spécifiques
1.4. L ES CARACTÉRISTIQUES PHYSIQUES DU GIVRAGE ET SA MODÉLISATION
1.4.1. Le givrage
1.4.2. Types et formes de glace
1.4.3. Modèles numériques de givrage
1.4.4. Résultats aérodynamiques des profils givrés
1.4.4.1. Impact de givrage et de la forme du dépôt de glace
1.4.4. 2. Impact de la rugosité
1.4.5. Courbes de puissance des éoliennes givrées
1.4.6. Pertes de production annuelle dues au givrage
1.5. CONTENU DU MÉMOIRE
CHAPITRE 2 MODÈLE NUMÉRIQUE DE PRODUCTION ANNUELLE D’UNE ÉOLIENNE À AXE HORIZONTAL
2.1 .INTRODUCTION
2.2. CALCUL DE LA PRODUCTION ANNUELLE D’UNE ÉOLIENNE
2.2. 1. Calcul de l ‘énergie extraite par une éolienne
2.2.2. Production énergétique non corrigée
2. 2.2.1. Distribution de la vitesse du vent
2. 2.2.2. Courbe d ‘énergie
2.2.2.3. Production énergétique non corrigée
2.2.3. Production énergétique corrigée
2.3. logiciel PROPID
2.4. LOGICIEL XFOIL
2.4.1. Description générale
2.4.2. Modèle physique
2.4.3. Modèle numérique
2.4.4. Inconvénients de XFO/L
2.5. L OGICIEL CIRALIMA
2.6. CONCLUSION
CHAPITRE 3 MODÈLE ÉOLIEN GÉNÉRIQUE REPRÉSENTATIF D’UNE ÉOLIENNE VESTAS V80 1.8 MW
3.1. Introduction
3.2. ÉOLIENNE V ESTAS V80 1.8 MW
3.3. PALE DU MODÈLE ÉOLIEN GÉNÉRIQUE
3.3.1. Passage de la pale 3D au profil aérodynamique
3.3.2. Vitesses et angles sur un élément de pale
3.3.3. Choix du profil. Distribution de l’épaisseur le long de la pale
3.3.4. Distribution de la corde le long de la pale
3.3.5. Distribution de l’angle de vrillage le long de la pale
3.3.6. Discrétisation de la pale
3.3.7. Synthèse
3.4. CALCUL DES COEFFICIENTS DE PORTANCE ET DE TRAÎNÉE
3.5. CALCUL DE LA COURBE DE PUISSANCE AUX CONDITIONS ATMOSPHÉRIQUES STANDARDS
3.6. VALIDATION DU MEG
3.7. C ALCUL DE LA Production ANNUELLE
3.7.1. Calcul de la production annuelle de la Vestas V80 1.8 MW
3.7.2. Calcul de la production annuelle du MEG
3.7.3. Validation de PRO/CET
3.8. CONCLUSION
CHAPITRE 4 SIMULA TION NUMÉRIQUE DE LA COURBE DE PUISSANCE DU MEG ET CALCUL DE SA PRODUCTION ANNUELLE SOUS PRÉCIPITATIONS GIVRANTES
4.1. INTRODUCTION
4.2. CONDITIONS MÉTÉOROLOGIQUES AU SOL
4.3. CONDITIONS ATMOSPHÉRIQUES DE GIVRE
4.4. TEMPÉRATURE DE _ 6 oC
4.4.1. Calcul des formes de glace
4.4.2. Calcul des coefficients de portance et de traînée
4.4.3. Calcul de la courbe de puissance
4.5. TEMPÉRATURE DE – 2° C
4.5.1. Calcul des formes de glace
4.5.2. Uniformisation de la surface de glace
4.5.3. Calcul des coefficients de portance et de traînée
4.5.4. Calcul de la courbe de puissance
4.6. CALCUL DE LA PRODUCfION ANNUELLE DU MEG SOUS PRÉCIPITATIONS GIVRANTES
4.7. CONCLUSION
CHAPITRE 5 CONCLUSION
TRAVAUX FUTURS
ANNEXE A CARACTÉRISTIQUES DE L’ÉOLIENNE VESTAS VSO 1.8 MW 100
ANNEXE B IMPACT DE LA TEMPÉRATURE SUR LA PUISSANCE ET LE COEFFICIENT DE PUISSANCE D’ UNE ÉOLIENNE
B.l . I NTRODUCTION
B.2. I MPACT DE LA TEMPÉRATURE
B.3. CONCLUS ION
ANNEXE C CARACTÉRISTIQUES GÉOMÉTRIQUES D’UN PROFIL. PROFILS DE TYPE NACA
c. l. C ARACTÉRISTIQUES GÉOMÉT RIQ UES D’UN PROFIL
C.2. PROFILS DE TYPE NACA
C.2.1. Famille de profils à quatre chiffres
C.2.2. Famille de profils à cinq chiffres
C.2.3. La série 6 de profils NACA
ANNEXE D CALCUL DE LA PRODUCTION ANNUELLE D’UNE ÉOLIENNE AVEC LE LOGICIEL RETSCREEN
ANNEXE E SIMULA TION NUMÉRIQUE DES COEFFICIENTS DE PORT ANCE ET DE TRAÎNÉE AVEC LE CODE COMMERCIAL XFOIL
INTRODUCTION
E.2. ÉTUDE DE SENSIBILITÉ
E.2.1. Effet du nombre d’itération
E.2.2. Effet du nombre de Mach
E.2.3. Effet du nombre de Reynolds
E.2.4. Validation des résultats
F.3 IMPACT DE LA T EMPÉRATURE
EA. CONCLUS ION
ANNEXE F
SIMULATION NUMÉRIQUE DES COURBES DE PUISSANCE AVEC LE LOGICIEL PROPID
F.l . I NTRODUCTION
F.2. I MPACT DE LA T EMPÉRATURE
F.2. 1. Courbe de puissance
F.2.2. Validation des résultats
F.3. C ONCLUSION
BILIOGRAPHIE

Cours gratuitTélécharger le document complet

Télécharger aussi :

Laisser un commentaire

Votre adresse e-mail ne sera pas publiée. Les champs obligatoires sont indiqués avec *