TRAVAUX PRATIQUES
GISEMENT SOLAIRE
Travail effectué par :
Soraya LIBRAHIMI
TANAN Abderrahmane
DOUADI Oumaima
Master spécialisé en sciences de l’eau, l’énergie et l’environnement.
Sommaire :
MANIP I- Mesure et estimation du rayonnement solaire……….03
MANIP II- Capteur solaire Plan ……………………………………………….12
MANIP III-Gisement photovoltaïque………………………………………18
2
MANIP I
Mesure et estimation du rayonnement solaire
Objectif :
Le but de cette manipulation est de comparer les mesures, pour une
période déterminée de la journée 07/06/2012, du rayonnement solaire globale
arrivant au sol, sur une surface horizontale, via le Pyranomètre avec les
estimations obtenues basées sur l’utilisation des modèles et des résultats
connues.
Partie théorique :
1. Coefficient de trouble de linke :
C’est un coefficient qui caractérise les phénomènes d’absorption et de
diffusion atmosphérique
Il se déduit de la connaissance des divers paramètres qui les engendrent
TL= 2,5+16 βA+0,5 ln w
TL= 1,6+17 βA +0,5 ln e
Compte tenu des humidités normalement relevés dans l’atmosphère
TL=2,6+16 βA(± 0,3) pour les régions tempérées
TL= 3,2+17 βA(± 0,4) pour les régions intertropicales
Si on ne dispose ni de visibilité ni de tension de vapeur approximative
TL=2,5+0,4 en régions tempérées
TL=3,2+0,5 en régions intertropicales
2. Estimation du rayonnement solaire globale :
Pour estimer l’éclairement solaire global reçu par une surface horizontale
on utilise les trois modèles suivants :
Modèle 1 :
3
Gmax =1150 (sin h) 1,15
Gmoy=1080 (sin h) 1,22
Gmin =990 (sin h) 1,25
Modèle 2 :
G = (1270-56 TL)*(sin h) (TL+36)/33
Modèle 3 :
Imax =1210 exp [-1/ (6 sin (h+1))]
- Gmax = Imax + Dmax
Dmax =87 (sin h) 0, 4
Imoy=1230 exp [-1/(3, 8 sin (h+1, 6))]
- Gmoy= Imoy + Dmoy
Dmoy=125 (sin h) 0, 4
Imin =1260 exp [-1/(2, 3 sin (h+3))]
- Gmin =sin(h) Imin + Dmin
Dmin =187 (sin h) 0, 4
3. Mesure de l’éclairement réellement reçu
Partie pratique :
À partir des modèles présentés dans la partie théorique et des résultats de
mesures obtenus dans la partie pratique, on se propose de modéliser, pour la
journée de 07/06/2015 les éclairements globaux reçus au site de l’ENSET Rabat sur
une surface horizontal.
4
1. Evaluation du coefficient de trouble de linke TL :
Etant donné que cette manipulation est programmée pendant l’été et dans
une zone tempérée, donc :
TL = 2,5+0,4
Donc : TL =2,92
2. Estimation de l’éclairement globale d’une surface horizontale
Pour la journée 07/06/2015, à midi TSV, on estime l’éclairement global d’une
surface horizontale par les trois modèles. Les résultats sont regroupés dans le tableau
ci-dessous :
Modèle 1 Modèle 2 Modèle 3
Imax= 972,305 Dmax= 85,44
Gmax = 1091,66 W/m2 G =1048,9999 W/m2 Imoy= 894,659 Dmoy = 122,758
Gmoy= 1022,0166 W/m2 Imin= 769,19 Dmin= 183,64
Gmin = 935,5367 W/m2
Gmax = 1057,745 W/m2
Gmoy= 1017,417 W/m2
Gmin = 952,83 W/m2
2
Tableau 1 : Eclairement global (W/m ) de la journée sujet d’étude à midi TSV site de l’ENSET Rabat
Pour calculer ces éclairements estimés nous avons besoin de trouver la hauteur de
soleil h calculé comme suit :
J = 31+28+31+30+8 donc J = 128
δ = 23,45 sin (0,986 (J + 284)) donc δ = 16,93°
ω = 15(TSV -12) donc ω = 0°
Rabat (ɸ = 34°2 N, L = 6°48 E), alors :
h =Arc sin [sin(16,93)*sin (34+2/60)+cos( 16,93)*cos( 34+2/60)]
Donc h = 72,8966°
5
3. Les valeurs de l’éclairement globale, reçu à midi TSV, mesurés par
l’enregistreur
L’appareil de mesure utilisé est le Pyranomètre, cet instrument permet d’évaluer l’ensemble
du rayonnement solaire (direct et diffus).
Figure(1) : Pyranomètre
A midi TSV sur un plan incliné, l’enregistreur nous a donné la valeur :
G = 765,7227 w/m2
Pour pouvoir la comparer avec celles des trois modèle on doit la multiplier pas le cos(i) où i
est l’angle formé par la direction du soleil et la normal à la surface réceptrice qui est donné
par la relation suivante :
Cos(i)=cos(h)sin(s)cos(a-ɤ)+sin(h)cos(s)
h est la hauteur du soleil,
a est l’azimut correspondant,
s est l’inclinaison sur l’horizontale,
ɤ est l’orientation de la surface par rapport au sud géographique, elle varie de 0
à - 180° du sud vers le nord dans le sens trigonométrique et, de 0 à 180° dans le
sens direct.
Donc sur un plan horizontal, l’éclairement est égal à :
G = 556,40 w/m2
6
4. Quel est le modèle le plus adapté à l’estimation de l’éclairement
Selon les résultats obtenus, le modèle le plus adapté est le modèle 1
Gmin = 935,5367 W/m2
Qui correspond à:
Gest = Gmin = 990(sin h) 1,25
5. Représentation graphique du Gest(t) et Gmes(t) :
Voir l’annexe (2)
6. Représentation graphique De Gmes(t)=f(Gest(t)) :
Voir l’annexe (2)
7. Ajustement du modèle choisi :
Pour ajuster le modèle choisi afin de le rendre plus proche de la réalité, le paramètre le
paramètre dont il faut agir principalement est le S l’inclinaison sur horizontale, ce dernier
intervient dans la formule cos(i)
Cos(i)=cos(h) sin(s)cos (a-ɤ) +sin(h) cos(s)
8. Conclusion générale
Au final, on conclu que La variation de l’inclinaison au cours de la journée, a pour
conséquence une grande influence sur l’éclairement reçu. D’où le point fort d’un système
suiveur, guidé par la trajectoire du soleil au fils des heures tout au long de la journée,
permettant ainsi d’exploiter le maximum de potentiel d’irradiations solaire.
7
Annexe (1)
Modèle 1 :
Eclairement correction Eclairement
Eclairement
mesuré sur un cos(i) mesuré sur un
enregistreur angle hauteur du estimé sur un
Temps légal TL (h) TSV plan incliné plan horizontal Azimut
(mv) horaire w soleil plan
(w/m2) (w/m2)
horizontal
07:15 7,25 6,7291 0,0898 16,8842 0,9088 15,3444 -79,0630 18,2775 -75,4408 232,348
08:15 8,25 7,7291 0,1737 32,6418 0,6111 19,9463 -64,0630 30,6469 -78,0804 426,410
09:15 9,25 8,7291 0,5715 107,4211 0,9986 107,2683 -49,0630 43,0210 -75,2861 613,872
10:15 10,25 9,7291 2,0073 377,3064 0,8441 318,4727 -34,0630 55,0273 -66,1624 771,830
11:15 11,25 10,7291 3,3409 627,9859 0,8553 537,1364 -19,0630 65,7657 -48,1432 882,173
12:15 12,25 11,7291 3,7969 713,7105 0,6057 432,2735 -4,0630 72,4234 -12,6870 932,567
13:15 13,25 12,7291 4,0736 765,7227 0,6604 505,7194 10,9370 70,2066 31,6320 917,436
14:15 14,25 13,7291 4,0492 761,1203 0,8558 651,3425 25,9370 61,1085 57,9285 838,451
15:15 15,25 14,7291 3,7621 707,1607 0,8061 570,0446 40,9370 49,6072 71,1676 704,382
16:15 16,25 15,7291 1,5913 299,1165 0,8925 266,9667 55,9370 37,3681 77,3473 530,348
17:15 17,25 16,7291 0,5609 105,4342 0,9075 95,6853 70,9370 24,9621 77,3920 336,743
18:15 18,25 17,7291 1,3630 256,1992 0,8895 227,8814 85,9370 12,6783 73,1450 148,720
19:15 19,25 18,7291 0,3903 73,3675 0,9980 73,2234 100,9370 0,7512 66,8001 4,392
12 0 35,37 0 499,858
8
Modèle 2:
Eclairement correction Eclairement
mesuré sur un cos(i) mesuré sur un Eclairement
enregistreur plan incliné plan angle horaire hauteur du éstimé sur un
TL (h) tsv (w/m2) horizontal Azimuth
(mv) w soleil plan
(w/m2) horizontal
7,25 6,7195 0,0898 16,8842 0,9523 16,0792 -79,2073 18,1593 -81,4929 230,535
8,25 7,7195 0,1737 32,6418 0,9257 30,2171 -64,2073 30,5274 -89,4984 424,534
9,25 8,7195 0,5715 107,4211 0,9970 107,0967 -49,2073 42,9029 -81,3974 612,176
10,25 9,7195 2,0073 377,3064 0,8387 316,4397 -34,2073 54,9156 -69,3414 770,513
11,25 10,7195 3,3409 627,9859 0,8364 525,2199 -19,2073 65,6746 -49,8230 881,383
12,25 11,7195 3,7969 713,7105 0,8095 577,7222 -4,2073 72,3961 -13,4189 932,391
13,25 12,7195 4,0736 765,7227 0,7266 556,4024 10,7927 70,2703 32,0498 917,894
14,25 13,7195 4,0492 761,1203 0,8855 673,9925 25,7927 61,2119 59,8115 839,493
15,25 14,7195 3,7621 707,1607 0,9041 639,3348 40,7927 49,7227 75,1900 705,891
16,25 15,7195 1,5913 299,1165 0,7685 229,8688 55,7927 37,4872 85,6158 532,152
17,25 16,7195 0,5609 105,4342 0,9141 96,3735 70,7927 25,0812 85,9163 338,624
18,25 17,7195 1,3630 256,1992 0,8154 208,8971 85,7927 12,7951 78,0637 150,407
19,25 18,7195 0,3903 73,3675 0,7396 54,2596 100,7927 0,8635 70,0257 5,227
9
Annexe 2
Eclairement mesuré et estimé en fonction du TSV
G mesuré G éstimé
1000,0000
900,0000
800,0000
700,0000
600,0000
500,0000
400,0000
300,0000
200,0000 TSV
100,0000
0,0000
10
Eclairement mesuré en fonction de l’éclairement estimé
Gmes
700,0000
600,0000
500,0000
400,0000
300,0000
200,0000
G est
100,0000
0,0000
11
MANIP II
Capteur solaire Plan
Objectif :
Le but de cette manipulation est l’étude expérimentale des pertes thermiques
et l’évaluation du rendement instantané du capteur solaire plan sous les conditions
climatiques de la journée (08/05/2015).
Partie théorique :
1- Bilan énergétique d’un capteur solaire plan :
Pu = Pab – Pp
Pu : c’est la puissance utile du capteur
Pu = ṁ C (Ts-Te)
ṁ : débit volumique de l’eau
C : chaleur massique de l’eau
Ts et Te : les températures d’entrée et sortie du capteur
Pab : c’est la puissance transmise par la couverture de verre et absorbée par
l’absorbeur
Pab = αab*τv*G
αab : coefficient d’absorption
τv : coefficient de transmission de verre
G: éclairement global incident sur la face extérieur
12
Pp : c’est les pertes thermiques totales du capteur
Pp = K(Tab – Ta)
K : coefficient global des pertes
Tab : température de l’absorbeur
Ta: température ambiante
2- Rendement du capteur solaire :
Il est donné par le quotient de l’énergie utile par l’énergie du rayonnement
ɳ= ou encore ɳ = ɳ0 -
Avec
ɳ0 : le rendement optique égal à 0,77
Partie expérimentale :
1- Mode opératoire :
Avant chaque série de mesure, il faut
Ouvrir complètement le robinet de débitmètre
Régler le débit avec le robinet
Chaque série de mesure s’étend sur toute la journée avec un pas d’une heure
Pour chaque série de mesure, on place le capteur dans une position optimale
par rapport au rayonnement solaire
2- Exploitation des résultats :
1) Pour avoir une position optimale il faut que le rayonnement solaire soit perpendiculaire à
la surface du capteur.
2) Voir le tableau en annexe 1 qui donne le rendement instantané par l’utilisation des deux
formules suivantes :
ɳ = ɳ0 - et ɳ=
13
3) Le rendement instantané moyen du capteur déduit du tableau est égale à : 0,418
4) Les températures de départ de panneau sont très proches sur l’intervalle (00 :00 à
09 :00) parce que à cet intervalle du temps l’éclairement est très faible et à ce moment la
température du ballon est plus chaude que celle du capteur et donc la pompe ne
fonctionne et par conséquent l’eau existant entre le panneau et le ballon ne circule pas
et sa température reste relativement constante.
Conclusion générale
D’après la courbe représentant le rendement instantané au long de toute la
journée de 07/11/2015, on peut remarquer que le rendement du capteur n’apparait
que durant la période allant de 10h à 17h où les valeurs intéressants de l’éclairement
permet de mettre la pompe en marche puisque la température du capteur devient
plus grande que celle du ballon.
La deuxième courbe donnant la variation du rendement instantané du
capteur en fonction de la différence de température de l’absorbeur et celle de
l’ambiance sur l’éclairement. On peut remarquer que plus les pertes augmentent
plus le rendement diminue et inversement, d’où l’importance de chercher à les
réduire.
14
(Tc- m
Tc Ta Ta) G :débit(m3/s) Pu(w) rendement
00:00 21,8 22 -0,2 0 0,09 0 0
00:15 21,8 22 -0,2 0 0,09 0 0
00:30 21,8 22,1 -0,3 0 0,09 0 0
00:45 21,8 22 -0,2 0 0,09 0 0
01:00 21,8 21,9 -0,1 0 0,09 0 0
01:15 21,7 21,9 -0,2 0 0,09 0 0
01:30 21,8 22 -0,2 0 0,09 0 0
01:45 21,7 21,9 -0,2 0 0,09 0 0
02:00 21,7 21,9 -0,2 0 0,09 0 0
02:15 21,6 21,8 -0,2 0 0,09 0 0
02:30 21,6 21,8 -0,2 0 0,09 0 0
02:45 21,5 21,8 -0,3 2 0,09 0 0
03:00 21,5 21,7 -0,2 0 0,09 0 0
03:15 21,5 21,8 -0,3 0 0,09 0 0
03:30 21,5 21,7 -0,2 0 0,09 0 0
03:45 21,5 21,7 -0,2 0 0,09 0 0
04:00 21,4 21,6 -0,2 0 0,09 0 0
04:15 21,4 21,6 -0,2 0 0,09 0 0
04:30 21,4 21,6 -0,2 0 0,09 0 0
04:45 21,3 21,6 -0,3 0 0,09 0 0
05:00 21,4 21,6 -0,2 0 0,09 0 0
05:15 21,3 21,5 -0,2 0 0,09 0 0
05:30 21,3 21,5 -0,2 0 0,09 0 0
05:45 21,3 21,5 -0,2 0 0,09 0 0
06:00 21,2 21,4 -0,2 0 0,09 0 0
06:15 21,2 21,4 -0,2 0 0,09 0 0
06:30 21,2 21,4 -0,2 0 0,09 0 0
06:45 21,2 21,4 -0,2 2 0,09 0 0
07:00 21,1 21,3 -0,2 9 0,09 0 0
07:15 21,1 21,4 -0,3 16 0,09 0 0
07:30 21,1 21,3 -0,2 20 0,09 0 0
07:45 21,2 21,3 -0,1 25 0,09 0 0
08:00 21,1 21,3 -0,2 29 0,09 0 0
08:15 21,2 21,4 -0,2 33 0,09 0 0
08:30 21,3 21,5 -0,2 35 0,09 0 0
08:45 21,3 21,5 -0,2 41 0,09 0 0
09:00 21,4 21,6 -0,2 45 0,09 0 0
09:15 21,5 21,7 -0,2 583 0,09 0 0
09:30 24,4 33,4 -9 626 0,09 0 0
09:45 36,4 37,2 -0,8 675 0,09 0 0
10:00 42,2 38,4 3,8 735 0,09 211,27 0,287442177
10:15 44,9 39,3 5,6 815 0,09 351,14 0,430846626
15
10:30 46,7 40,1 6,6 889 0,09 347,12 0,390461192
10:45 47,7 40,7 7 911 0,09 415,32 0,455894621
11:00 48,3 41,4 6,9 948 0,09 456,2 0,481223629
11:15 49 42 7 987 0,09 458,3 0,464336373
11:30 49,6 42,8 6,8 1026 0,09 585 0,570175439
11:45 49,9 43,4 6,5 1059 0,09 601,36 0,567856468
12:00 50,4 44 6,4 1094 0,09 615,02 0,562175503
12:15 50,9 44,6 7,3 1179 0,09 679,08 0,575979644
12:30 51 45,1 5,9 1252 0,09 745 0,595047923
12:45 51,3 45,6 5,7 1259 0,09 754,02 0,598903892
13:00 51,9 46 5,9 1204 0,09 601,5 0,499584718
13:15 52,3 46,5 5,8 1132 0,09 564,2 0,498409894
13:30 52,4 46,7 5,7 1083 0,09 520 0,480147738
13:45 52,9 47,1 5,8 1077 0,09 500,01 0,464261838
14:00 52,7 47,4 5,3 1026 0,09 450 0,438596491
14:15 52,6 47,7 4,9 924 0,09 402,01 0,435075758
14:30 52,2 47,9 4,3 800 0,09 332,1 0,415125
14:45 52,1 48 4,1 705 0,09 281 0,39858156
15:00 51,9 47,9 4 648 0,09 232 0,358024691
15:15 51,6 47,9 3,7 611 0,09 200,2 0,327659574
15:30 51,3 47,9 3,4 584 0,09 165 0,282534247
15:45 50,6 47,6 3 565 0,09 135,25 0,239380531
16:00 49,7 47,7 2 511 0,09 112,36 0,219882583
16:15 48,3 47,4 0,9 421 0,09 90,3 0,214489311
16:30 47,1 47,2 -0,1 314 0,09 65 0,207006369
16:45 44,4 45,8 -1,4 211 0,09 42 0,199052133
17:00 40 41 -1 109 0,09 0 0
17:15 36,7 37,1 -0,4 20 0,09 0 0
17:30 34 34 0 6 0,09 0 0
17:45 31,8 31,7 -0,1 2 0,09 0 0
18:00 30 29,8 0,2 0 0,09 0 0
18:15 28,5 28,4 0,1 0 0,09 0 0
18:30 27,3 27,3 0 0 0,09 0 0
18:45 26,3 26,4 -0,1 2 0,09 0 0
19:00 25,6 25,7 -0,1 0 0,09 0 0
19:15 25,1 25,2 -0,1 0 0,09 0 0
19:30 24,5 24,7 -0,2 0 0,09 0 0
19:45 24,1 24,3 -0,2 0 0,09 0 0
20:00 23,8 24 -0,2 0 0,09 0 0
20:15 23,5 23,8 -0,3 0 0,09 0 0
20:30 23,2 23,5 -0,3 0 0,09 0 0
20:45 23,1 23,3 -0,2 0 0,09 0 0
21:00 23 23,2 -0,2 0 0,09 0 0
21:15 22,9 23,1 -0,2 0 0,09 0 0
16
21:30 22,7 23 -0,3 0 0,09 0 0
21:45 22,6 22,9 -0,3 0 0,09 0 0
22:00 22,6 22,8 -0,2 0 0,09 0 0
22:15 22,6 22,9 -0,3 0 0,09 0 0
22:30 22,5 22,8 -0,3 0 0,09 0 0
22:45 22,4 22,7 -0,3 0 0,09 0 0
23:00 22,4 22,7 -0,3 0 0,09 0 0
23:15 22,3 22,6 -0,3 0 0,09 0 0
23:30 22,3 22,5 -0,2 0 0,09 0 0
23:45 22,3 22,5 -0,2 0 0,09 0 0
le rendement (ȵ)
0,7
0,6
0,5
0,4
Série1
0,3
0,2
0,1
0
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28
Autres calculs et données :
k 4
n0 0,77
Cp j/kg.k 3850
ro kg/m3 1029
17
MANIP III
Gisement photovoltaïque
Objectif :
Le but de ce TP est d’étudier l’influence de l’angle d’inclinaison, la distance et la résistance sur la
puissance produite par la cellule photovoltaïque.
Partie Pratique
Tout d’abord nous avons réalisé le montage de circuit électrique qui contient :
La cellule PV
La résistance variable
L’ampèremètre
Le voltmètre
1. Influence la résistance sur la puissance :
Tout d’abord nous avons modifié la résistance dans les conditions suivantes :
Distance entre la cellule et la lampe : 100 cm
L’angle d’inclinaison : 90°
L’éclairement : 500 w/m²
Les résultats sont représentées dans le tableau suivant :
18
R(Ω) 1 5 10 15 18 19 20
I (A) 120 116 109 97,2 94,3 93 92
U(mV) 1,14 1,34 1,64 1,75 2,28 2,3 2,39
P (w) 136,8 155,44 178,76 170,1 215,004 213,9 219,88
0,25
0,2
0,15
0,1
puissance
0,05
0
4 6 9 10 18 19 20 21 40 Résistance
2. Influence de l’inclinaison sur la puissance :
Pour déterminer l’influence de l’inclinaison on a varié l’inclinaison tout en gardant les paramètres
suivants fixes :
Distance entre la cellule et la lampe : 100 cm
L’éclairement : 500 w/m²
La résistance : 20 Ω
Nous avons obtenu les résultats suivants :
19
Inclinaison 0 30 45 70 90
I 4 31,3 60,1 85,7 92
V 0,09 0,81 1,48 2,19 2,39
P (w) 0,00036 0,025353 0,088948 0,187683 0,21988
0,25
0,2
0,15
0,1 puissance
0,05
0
0 30 45 70 90 inclinaison
3. Influence de la distance sur la puissance :
Pour se faire nous avons varié la distance entre la source et la cellule dans les conditions
suivantes :
Inclinaison : 90°
L’éclairement : 500 w/m²
La résistance : 20 Ω
D (cm) 100 80 70 50
102, 120,
I (A) 90 6 111 9
V (V) 2,8 2,76 2,9 3,51
P (W) 0,252 0,283176 0,3219 0,424359
20
puissance
0,45W
0,4
0,35
0,3
0,25
0,2
0,15
0,1
0,05
0
100 80 70 50
la distance cm
Conclusion
On constate que pour avoir un bon fonctionnement de la cellule photovoltaïque, il faut que :
Le rayonnement solaire soit perpendiculaire à la surface de la cellule ;
La distance entre la source du rayonnement et la surface réceptrice soit minimale.
21