FILIERE :
Master Spécialisée En Energies
Renouvelables Et Systèmes Energétiques
Module :
SOLAIRE THERMIQUE
Manipulation II : Eclairement direct S*
Sous le thème :
GISEMEMT SOLAIRE
Réalisé par :
AIT BAHA Ibtissam
Encadré par :
[Link] Ayoub
I. Partie théorique
Rayonnement solaire au sol
Notations Comme nous l’avons évoqué précédemment, l’atmosphère ne transmet pas au sol la
totalité du rayonnement solaire qu’elle reçoit :
- Le rayonnement direct est celui qui traverse l’atmosphère sans subir de modifications.
- Le rayonnement diffus est la part du rayonnement solaire diffusé par les particules solides ou
liquides en suspension dans l’atmosphère. Il n’a pas de direction privilégiée.
- Le rayonnement global est la somme du rayonnement direct et diffus.
Les notations utilisées pour les composantes du rayonnement solaire sur une surface horizontale
sont les suivantt
I=Energie reçue pendant une certaine durée( W.m-2. durée-1 ou [Link]ée )
Directe S
Diffuse D
Globale G=S+D
I*=Eclairement solaire Flux instantané W.m-2
Directe S*
Diffuse D*
Globale G*=S*+D*
Rayonnement direct= Eclairement S*
L’éclairement solaire direct S* sur un plan horizontal peut être déterminé de plusieurs manières
en fonction des données disponibles :
Par mesure de G* et D*, on en déduit S* = G* - D*.
A partir de la mesure des irradiations journalières globales G et diffuse D sur un plan
horizontal, on en déduit S = G – D et S* par la fonction de répartition suivante :
Où :
a = 0,409 + 0,502 sin(ωl – 60°)
b = 0,661 - 0,477 sin(ωl – 60°)
A partir de la mesure de l’irradiation journalière globale G, on évalue l’irradiation
journalière diffuse D par la corrélation de Collares-Pereira et Rabl
D = 0,99 G KT ≤ 0,17
D =( 1,188 -2, 272 KT +9, 473 KT2− 21,865 KT3+ 14,648 KT4 )G 0,17 < KT ≤0,75
D = (− 0,54 KT + 0, 632) G 0,75 KT ≤0,80
D=0,2G KT ≥0,80
Où : KT=G/G0
G0 étant l’irradiation journalière sur un plan horizontal placé au-dessus de l’atmosphère
calculable par :
Où ωl est en degré et G0 en kJ.m-2
On calcule ensuite S = G – D et on est ramené au cas précédent.
A partir de la connaissance de la moyenne mensuelle de l’irradiation globale journalière
G, on calcule l’irradiation diffuse journalière moyenne D par la corrélation de Collares-
Pereira et Rabl :
Et on est ramené au cas 2.
A partir de la mesure du taux d’ensoleillement σ, on évalue G par :
G = G0 [0,29 cos(L)+ 0,52 σ] Zone tropicale
G= G0 [ √ 2 σ +1 − 0,72] Franc
Et on est ramené au cas précédent
On ne dispose d’aucune mesure : on peut évaluer le rayonnement direct sur un plan
perpendiculaire au rayonnement solaire par la relation :
Où TL est le facteur de trouble de Linke calculable par]
β est le coefficient de trouble atmosphérique que l’on peut prendre égal à :
β = 0,05 en zone rurale
β = 01 en zone urbaine
β = 0,2 en zone industrielle ou polluée
pv est la pression partielle de vapeur d’eau exprimée en mmHg.
On en déduit S *= I* sin(h)
II. PARTIE PROGRAMMMATION
Dans ce programme on va utiliser la dernière méthode pour calculer l’éclairement solaire direct
S*
% Determiner la hauteur et laZimut du soleil
clear;
close all
clc;
l=input ('Enter la longitude de votre ville ');
L=input ('Enter la la altitude de votre ville ');
jour =input ('enter le jour ');
mois=input ('enter le mois ');
wp=0.984;
k=pi/180;
Iref=0;
j=jour+31*(mois-1);
if(mois<3)
j=j;
else
j=j-floor(0.4*mois+2.3);
end
X=0;
for i=1:24
X=X+1 ;
TL(i)=X;
ET=(0.0001-0.4797*cos(j*wp*k)+3.2265*cos(2*j*wp*k)+0.0903*cos(3*j*wp*k)
+7.3509*sin(wp*j*k)+9.3912*sin(2*wp*j*k)+0.3361*sin(3*wp*k));
TS(i)=TL(i)+(ET/60)+(-1/15);
W(i)=15*(TS(i)-12);
dec=23.45*sin((360/365)*(j-81));
sinH(i)=sin(L*k)*sin(dec*k)+cos(L*k)*cos(dec*k)*cos(W(i)*k);
Hauteur(i)=asin(sinH(i))/k;
sinA=(((cos(dec*k)*sin(W*k))/cos(Hauteur(i)*k)));
Azimut=(asin(sinA))/k;
B=0.1;%(le coefficient de troubl eatmospherique :0.05(zone rurale ),0.1(zone
urbaine ),0,2(zone industriel):')
pv=15.8; %(depend de la temperature )
TRL=2.4+14.6*B+0.4*(1+2*B)*log(pv);
E(i)=1370*exp(-TRL/(0.9+9.4*sinH(i)));
S(i)=E(i)*sinH(i);
if (S(i)<0)
S(i)=0;
else
S(i)=S(i);
end
D(i)=54.8*sqrt(sinH(i))*(TRL-0.5-sqrt(sinH(i)));
if (D(i)<0)
D(i)=0;
else
D(i)=D(i);
end
G(i)=S(i)-D(i);
if (G(i)<0)
G(i)=0;
else
G(i)=G(i);
end
end
plot (TL,S,'b')
ylabel('Eclairement direst (W/m2')
xlabel('le temps legal en heure ')
title ('rayonnement direct ')
grid on
zones climatiques longitude latitude
zone 1 agadir -9.35 30.25
zone2 tanger -5.47 35.46
zone 3 Fes -5 34.01
ZONE 4 ifrane -5.11 33.53
ZONE5 marrakech -7.59 31.38
ZONE 6 errachidia -4.25 31.55
