Rayonnement Solaire Etudiants
Rayonnement Solaire Etudiants
Chapitre 1:
Le rayonnement solaire
Introduction générale
Le rayonnement solaire qui arrive
au niveau du sol peut être
décomposé en deux parties, la
première provient directement du
Rayons directs
disque solaire, la deuxième a été particules
diffusée par les éléments existants
dans l'atmosphère.
Le soleil est une sphère formée de gaz chaud fortement comprimé, ayant un diamètre de 1.39*109m. La distance moyenne entre la
terre et le soleil est de 1.5 * 1011 m.
Le soleil peut être assimilé a un corps noir ayant une température de 5762 K. Le centre du soleil a une température qui varie entre
8 * 106 et 40 * 106 K. La densité au centre du soleil étant de 100 fois celle de l'eau. L'énergie du soleil provient d'une fusion
thermonucléaire.
Chapitre 1:
La structure du soleil
Zone convective
0.7 ≤ r/rS ≤ 1 Photosphère
T ≈ 2 ∙ 106 K 1 ≤ r/rS ≤ 1.001
ρ ≈ 0.15 g/cm3 T ≈ 5800 K
Transport d‘énergie par convection ρ ≈ 2 ∙ 10-7 g/cm3
Emission du spectre
Zone radiative solaire en continu
0.23 ≤ r/rS ≤ 0.7
T ≈ 7 ∙ 106 K
ρ ≈ 4 g/cm3
Absorption et Chromosphère
emission 1.001 ≤ r/rS ≤ 1.01
T ≈ 4500 K
ρ ≈ 5 ∙ 10-9 g/cm3
Noyeau
(Coeur) Couronne
0 ≤ r/rS ≤ 0.23 1.01 ≤ r/rS ≤ 20
T ≈ 1.5 ∙ 107 K T ≈ 107 K
ρ ≈ 100 g/cm3 ρ ≈ 10-15 g/cm3
Fusion Nucléaire
Chapitre 1:
Soleil Solstice
Solstice 147 d’hiver
d’été 152 millions de millions (21
(22 Juin) km de km Décembre)
= +23.45 = -23.45
Equinoxe d’automne
(21 Septembre)
= -0
La terre décrit autours du soleil un mouvement de rotation sur une trajectoire elliptique dont le soleil occupe l'un des foyers. Une
rotation complète prend 365 j 1/4 (une Figureannée), ce mouvement
II-1-a: est
Trajectoire de la à l'origine
terre autourdes saisons (Fig. II-1-a). Le plan qui contiens la
du soleil.
trajectoire de la terre est appelé plan de l'écliptique
Chapitre 1:
Animation
2327 ’
Axe de
Normale au plan de
rotation
l’écliptique
de la terre
Tropique du cancer
équateur
Tropique du
capricorne
La terre effectue une rotation autours de son axe en 24 h. cette rotation est à l'origine du jour et de la nuit. L'axe de rotation de la
terre fais un angle de 23° 27' avec la perpendiculaire au plan de l'écliptique (Fig. II-1-b).
Chapitre 1:
Plan de l’équateur
360
= 23.45 * sin (n + 284) (1)
365
Avec :
la terre: ( ) , l'angle entre le plan de l'équateur et la direction terre soleil. La déclinaison peut être
On appelle déclinaison360
• JD = JD0 + :de ( ( ) ( ) ( )
. 0.007133. sin JD0 + 0.03268. cos JD0 − 0.000318. sin 2.JD0 + 0.000145. cos 2.JD0 ( ))
calculée par l'équation 2
(1)
•n : •
étant JD ( )
0 = n − 81 *
le quantième
Décembre : n=365.
360
de l'année, en
c'estdegrés.
à dire le rang du jour dans l'année, le premier Janvier étant le premier jour (n = 1), le 31
365
360
•
Si l'on veut calculer avec une meilleure précision, on utilise l'équation :
: pour équilibrer la rotation de la terre : 360° en 365 jours.
365
= 23.45 sin JD.
Chapitre 1:
Au cours d'une journée, JD peut être considéré comme une constante. La figure II-2 donne la valeur de JD en fonction des jours de
l’année.
Chapitre 1:
= 1.27*107m
= 1.39*109m
32
°
En dehors de
l'atmosphère, la densité
de flux qui parvient au
1.495 * 1011 m 1.7 %
niveau de la terre est
donnée par : Fig II-3
L
E=
4..a.2 (1)
Le soleil rayonne un flux de L = 4 * 1026 W dans tout l'espace. En dehors de l'atmosphère, la densité de flux qui parvient au niveau
de la terre est donnée par :[1]
Chapitre 1:
Eclairement extraterrestre
Éclairement extra-atmosphérique
1440
Eclairement [w/m2]
1420
1400 Figure II-4 :
Variation de
1380
l’éclairement extra-
1360 atmosphérique au
cours de l’année.
1340
1320
0 50 100 150 200 250 300 350
jours de l'annee
360.n
E = E 0 .1 + 0.033. cos
➢n étant le quantième.
La valeur approximative de l’éclairement au cours de l'année 365
[1]
peut être exprimée par la fonction :[1]
Les mesures effectuées grâce aux satellites ont montré que la distribution spectrale du rayonnement solaire se fait selon la
Figure II-5 : Distribution spectrale du rayonnement solaire en dehors de l’atmosphère et au niveau de la terre.
figure II-5. On constate que les différents gaz de l’atmosphère absorbent certaines longueurs d’ondes.
La fraction la plus importante de l’éclairement solaire se trouve dans le spectre visible.
Chapitre 1:
réflection
Diffusion de rayleigh molécules
Et absorption (ca. 15%) de l‘air
06/10/2023
Chapitre 1:
AM = 2.5 AM = 2 AM = 1.5 AM = 1
Plus exactement:
θz
h … la hauteur du soleil
Rayonnement direct, Diffus et réfléchi
Rayonnement direct
Rayonnement diffus
G=Gb+Gd+Gr
06/10/2023
Rayonnement Direct : turbidité
Relation Empirique:
06/10/2023
Eclairement sur un plan incliné
Plan perpendiculaire Plan Horizontal Plan incliné
au rayonnement
direct
direct (beam)
radiation
θ
z … Angle Zénital
… Angle d‘incidence sur le plan incliné θz θz
L‘éclairement est
inversement
proportionnel à l‘aire A
de la surface irradiée
Eclairement
Gbn: éclairement direct (beam) normal, Gb: éclairement direct horizontal, Gbt: Ecl. direct sur un plan incliné
06/10/2023
Rayonnement indirect: diffus et Réfléchi
Estimation de l‘éclairement Estimation de l‘éclairement
diffus sur une surface réfléchi sur une surface
inclinée (inclinaison: ): inclinée (inclinaison: angle ):
06/10/2023
Bilan radiatif de la terre
La distribution de Planck pour le soleil (5780K, max. à 0.5 m) et de la
terre (288K, max. à 10 m) (les courbes ne sont pas à l‘échelle)
06/10/2023
Chapitre 1:
Rayonnement de la matière
(Vent solaire)
rayonnement électromagnétique
Trajectoire de la
areas under curves
Longueur d‘onde
indicate P(T)
Max
06/10/2023
Le Soleil: Un radiateur thermique
wave length nm
06/10/2023
2.2: Les coordonnées du soleil
Le rayonnement solaire
Introduction
06/10/2023
Les coordonnées du soleil
Introduction
Afin de localiser la
position du soleil en
chaque instant, nous
devons adopter un
système de coordonnées
dans lequel nous
décrirons la position du
soleil dans le ciel. On
distingue deux systèmes
de coordonnées
06/10/2023
Les coordonnées Azimutales ou horizontales
Le système de coordonnées horizontales
On définit le système de
coordonnées Azimutales en un
point donné de la surface de la
terre comme suit :
06/10/2023
Les coordonnées Azimutales ou horizontales
Les coordonnées du soleil dans un système de coordonnées horizontales
• Le soleil (ou n’importe quel point dans le ciel)
z
est alors décrit par sa hauteur et son azimut.
• La hauteur h : c'est l'angle entre l'horizontale
et le direction terre-soleil OS
a
• L'axe de rotation de la terre fais j avec OX et
se trouve dans le plan ( OX , OZ ).
x
y
06/10/2023
Coordonnées azimutales ou horizontales
Les coordonnées du soleil dans un système de coordonnées horizontales
système d'axes ( ,
OX OY , OZ ),
on a :
x = cosh . cos a
OS = y = cosh . sin a
=
z sinh
Figure II-7-a : Les coordonnées du soleil dans
un système de coordonnées horizontale
06/10/2023
Coordonnées équatoriales ou horaires
06/10/2023
Coordonnées Equatoriales ou Horaires
Rappel
06/10/2023
Coordonnées Equatoriales ou Horaires
Rappel
•L'angle horaire : = C'est l'angle entre les plans (OZ’ , OS’) et ( OZ’, OX’ )
• à midi solaire vrai (TSV), = 0
• augmente de 15° par heure, il a une valeur négative
Z’
06/10/2023
Coordonnées équatoriales ou Horaires
Rappel
06/10/2023
Coordonnées Horizontales
Rappel
On déduit donc :
• La hauteur maximale qui peut être obtenue dans un lieu donné arrive à midi TSV le
jour du solstice d'été lorsque = s = ’
• Le soleil est à la verticale du lieu (h = ) lorsque j = . cette condition ne peut avoir
2
lieu que lorsque s = , c'est à dire dans les zones situées entre le tropique
du cancer et le tropique du capricorne.
06/10/2023
Calcul de l’angle d’incidence sur un plan donné
06/10/2023
Calcul de l’angle d’incidence sur un plan donné
on trouve :
l = sin cos z
U = m = sin . sin
S
n = cos
On déduit alors :
n
Cos = U • OS θ
u x
n β α
u
β m
On se propose dans ce qui suit de calculerl l'angle d'incidence du rayonnement direct à un instant donné sur un plan ayant
l'orientation suivante : α y
m
Inclinaison : angle entre le plan en question et le plan horizontal.
Azimut : c'est l'angle entre la normale au plan et le plan méridien.
Si on cherche les cosinus directeurs de la normale au plan, on trouve:
06/10/2023
Calcul de l’angle d’incidence sur un plan donné
06/10/2023
Calcul de l’angle d’incidence sur un plan donné
Remarques
06/10/2023
Le temps solaire Vrai
Introduction
06/10/2023
Le temps solaire Vrai
La Première correction:
La première correction est due à la différence de temps entre le méridien où est situé l'observateur et le méridien pris comme
méridien de référence pour le calcul du temps local. En effet, le soleil prend 4 minutes pour traverser 1° de longitude.
Le temps local est le temps du méridien qui passe au centre du fuseau horaire dans lequel est situé le lieu d'observation. Un fuseau
horaire est la surface terrestre entre deux méridiens séparés de 15°. Le fuseau horaire de référence est celui qui s'étend entre 7.5°
Est et 7.5° Ouest du méridien de Greenwich.
06/10/2023
Temps solaire Vrai
La deuxième correction est due au fait que la vitesse de rotation de la terre autours du soleil n'est pas constante. Cette perturbation
est causée par le fait que la trajectoire de la terre est elliptique et non circulaire. L'inclinaison de l'axe de rotation de la terre par
rapport à la normale au plan de l'écliptique interviens aussi pour perturber la longueur de la journée. Supposons que la trajectoire de
la terre est circulaire. Du fait que la terre effectue deux mouvements, le premier autours d'elle-même et le second autours du soleil,
le soleil ne repasse au zénith qu'après avoir décrit une rotation de (360° + ). L'angle peut être calculé par l'équation :
ts 1
= V •
avec : d
•V : la vitesse de la terre sur son orbite;
•ts : l'intervalle de temps entre deux passages du soleil au zénith;
•d : la distance terre-soleil.
06/10/2023
Correction de l’équation du temps
Loi de kepler
La vitesse de la terre ( V ) et la distance terre-soleil ( a ), varient au cours de l'année (selon la 2ième loi de Kepler), le produit (V . a)
est une constante.
Pour tenir compte de l'inclinaison de l'axe des pôles par rapport au plan de l'écliptique, il suffit de remplacer V par sa projection
Loi de Kepler
orthogonale sur le plan de l'équateur, ainsi une correction appelée correction de l'équation du temps (ET) doit être effectuée pour
tenir compte de ce phénomène.
ET = 9.87 • Sin(2.JD) − 7.53 • Cos ( JD) − 1.5 • Sin( JD)
ET : est exprimée en minutes [mn] et fraction décimale de mn. Elle varie au cours de l'année entre (-14.3) mn et (+16.4) mn.
360
JD = JD0 + 2 [0.007133 Sin JD0 + 0.032680 Cos JD0 - 0.000318 Sin 2JD0 + 0.000145 Cos 2JD0].
360
JD0 = (n − 81).
365 [degrés].
06/10/2023
Correction de l’équation du temps
L'allure de la courbe de l'équation du temps est donnée par la courbe de la figure II-11 suivante :
Ainsi, le temps solaire vrai est exprimé comme suit :
TSV = TL – 4 . (L – Lref) + ET
avec :
•TL : le temps légal : est le temps lu au niveau des montres et adapté officiellement au niveau d'un pays. D'habitude c'est le
temps du fuseau horaire (sauf exceptions).
•L = longitude du lieu.
•Lref = longitude du méridien de référence pris pour le calcul du temps légal.
•ET : correction de l'équation du temps.
•4 : [mn /°].
06/10/2023
Temps solaire vrai
Le tableau II-1 donne les valeurs de la déclinaison, de l’équation du temps ainsi que du décalage (en minutes) entre le Temps
Solaire Vrai et le Temps légal pour la ville de TUNIS pour les 365 jours de l'année.
06/10/2023
Les diagrammes solaires
Objectifs de la leçon
06/10/2023
Les diagrammes solaires
Introduction
• Pour faciliter le travail des Ingénieurs et
architectes, et leur éviter les calculs fastidieux,
des diagrammes solaires et des abaques sont
proposés. L'abaque le plus utilisé étant celui qui
indique la hauteur et l'azimut en fonction des
heures de la journée et ceci mois par mois. Cet
abaque est utilisée par les architectes pour
calculer les ombres crées par le voisinage.
06/10/2023
Chapitre 1:
90°
80°
70°
60°
50°
40°
30°
20°
10°
00°
Figure II-12-a : Construction de l’abaque : la hauteur du soleil
Mesurée entre l'horizon et la position du soleil au dessus de ce dernier. Les lignes horizontales de l'abaque représentent la hauteur
du soleil avec un incrément de 10°.
Chapitre 1:
Construction du diagramme
L’azimut du soleil
les lignes verticales représentent l'azimut du soleil (par 5°) à l'Est et à l'Ouest du Sud vrai.
Chapitre 1:
Construction du diagramme
Projection de la voûte céleste
Diagramme solaire
Voûte céleste
Horizon
En traçant les lignes horizontales et verticales, on obtiens finalement une grille représentant les coordonnées dans le ciel. La
projection effectuée est une projection polaire.
Chapitre 1:
Construction du diagramme
La trajectoire du soleil
E S O
Une fois la hauteur et l'azimut du soleil connus, la position du soleil est localisée en chaque instant dans le ciel, et la trajectoire est
Figure
tracée pour n'importe quel jour II-12-d
de l'année. Sur:l'abaque,
Construction de l’abaque
sont tracées : la trajectoire
les trajectoires du soleil.
du soleil tout les 21 de chaque mois.
Construction du diagramme
L’abaque solaire de la ville de Tunis :
Lorsqu'un capteur est installé, les bâtiments, arbres et obstacles qui l'entourent peuvent créer des masques sur sa surface. Il est
important de connaître à quel instant ceci arrive, et de calculer l'importance de ces masques et leur effet sur l'énergie reçue à la
surface du capteur.
Chapitre 1:
Le capteur
Pour localiser l'effetest
deàces
l'ombre lorsque
obstacles, h < ho
il faut les reporter sur un diagramme solaire en transformant leurs altitude en hauteur (angle)
au dessus de l'horizon et leur largeur en azimut. La référence étant la position du capteur. Prenons le cas simple d'un bâtiment
ho étant défini
parallélépipédique : par : tg ho = tg hm • cos
Le capteur est à l'ombre lorsque h < ho ; ho étant défini par :
étant l'azimut du soleil
tg ho = tg hm • cos
étant l'azimut
Cette du soleil.
condition estCette condition
vérifiée lorsqueest vérifiée
2 lorsque
< < 1 , < < 1 , en dehors de cet intervalle, le soleil atteint le capteur
sans être caché. On trace les diagrammes d'ombre et on les reporte sur l'abaque d'ensoleillement. On déduit les heures ou le soleil
n'atteint pas le capteur.
Chapitre 1:
• Ciel isotrope
• ………..
Chapitre 1:
Id ID
On a donc : R= Rd + RD
I I
Un capteur solaire plan absorbe les deux composantes de l'éclairement solaire : le diffus et le direct. En général, on dispose de
l'éclairement sur le plan horizontal. Pour estimer l'éclairement sur un plan incliné il faut d'abord connaître la valeur du rapport R.
Chapitre 1:
L’éclairement direct
Pour le cas de l’éclairement direct, le rapport Rd peut être calculé avec précision, vu que ce
rayonnement est directionnel.
L'éclairement direct sur une surface inclinée étant :
cos
I di = I d cos = I dh
sinh
Avec :
• = angle entre la normale au plan et le rayon solaire (Angle d’incidence du rayonnement direct).
cos
d'où Rd =
sinh
Le meilleur angle azimutal pour un capteur solaire plan étant égal à zéro : (a = 0), on a donc :
L’éclairement diffus
Introduction
L’éclairement diffus
Modèle du ciel clair
On peut supposer que le ciel est clair, et que la grande partie du rayonnement
diffus vient d'un cercle proche du soleil (le ciel circumsolaire). Le diffus est à
ce moment très anisotrope. Dans ce cas, on peut traiter l'éclairement diffus de
la même manière que l'éclairement direct : R = Rd
Exemple :
Estimer la valeur de R pour une surface inclinée de 60° pour j= 60° à 9h30 du matin. Le
20 Février. On supposera que le rayonnement diffus est concentré autour du soleil.
L’éclairement diffus
Le modèle de Hottel et Woertz
Hypothèses:
→ Ciel nuageux, donc le rayonnement diffus est isotrope
→ Le sol réfléchit uniformément le rayonnement sur la surface inclinée
dans ce cas: RD = 1
IG = Id . Rd + ID
L’éclairement diffus
Le modèle de Liu et Jordan
Hypothèses:
Une surface inclinée reçoit du rayonnement:
• En provenance du disque solaire (direct)
• Diffus en provenance du ciel
• Réfléchi par le sol
éclairement = f(orientation)
180
160
Irradiation journalière
140
sud
120
moyenne
sud-ouest
100
ouest
80
nord-ouest
60
nord
40
20
0
0 1 2 3 4 5 6 7 8 9 10 11 12
mois
Chapitre 1:
Slide 62
www.dlr.d
Chapitre 1:
Mesures au sol
Slide 63 www.dlr.de/enerMENA
Les pyranomètres
• Plage de mesures: 0.3 μm ≤ λ ≤ 3 μm
• Principe de fonctionnement
• L’absorbeur noir sous les domes de verre est chauffé grace au rayonnement solaire
• Un contact thermique entre la surface noire de l’absorbeur et le corps extérieur du
pyranomètre
• La différence de température entre le centre de l’absorbeur noir et le bord est
mesurée grace à une thermopile.
Anneau pour
cacher le disque
solaire
Sphère pour
À ajuster cacher le disque
régulièrement solaire
(2 fois par semaine)
Avec système de
suivi
Slide 65 www.dlr.de/enerMENA
Pyrheliometer
- Même principe de
fonctionnement
qu’un Pyranomètre
(thermopile)
- Un système de suivi
du soleil est
nécessaire
Slide 67 www.dlr.de/enerMENA
Chapitre 1:
Slide 68 www.dlr.de/enerMENA
Chapitre 1:
Des mesures sur du long terme sont requises pour obtenir des données fiables.
20%
15%
10%
longue période
5%
0%
-5%
-10%
-15%
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
La mesure au sol demande de longues périodes de mesures et nécessite des moyens coûteux
source: Quaschning
Slide 69 www.dlr.de/enerMENA
Chapitre 1:
Données atmosphériques
additionnelles:
• distribution de l’ozone
• contenu en vapeur d’eau
• contenu en Aérosols
Slide 70 www.dlr.de/enerMENA
Chapitre 1:
Satellites géostationnaires
Meteosat 7
0° latitude; 0° longitude
Meteosat 5
0° latitude; 63° longitude
source: DLR
Slide 71 www.dlr.de/enerMENA
Chapitre 1:
Slide 72 www.dlr.de/enerMENA
MERCI DE VOTRE
ATTENTION
Chapitre 2: LE RAYONNEMENT SOLAIRE
06/10/2023 Retour
Chapitre 2: LE RAYONNEMENT SOLAIRE
06/10/2023 Retour