UNIVERSITE D’ANTANANARIVO

*******************
ECOLE SUPERIEURE POLYTECHNIQUE
D’ANTANANARIVO
*******************
Départements : Génie Mécanique Productique
Génie Electrique

Mémoire de fin d’étude en vue de l’obtention du diplôme d’ingénieur en Génie Industriel

ETUDE ET REALISATION D’UN CAPTEUR SOLAIRE A

ABSORBEUR ALVEOLAIRE

Soutenu par :
RAZAFINDRABE ANDRIANJAKA Manantsoa

Sous l’encadrement de :
Monsieur RANARIJAONA Jean Désiré,

Date de soutenance : 24 Septembre 2011

Promotion 2010
 UNIVERSITE D’ANTANANARIVO
*******************
ECOLE SUPERIEURE POLYTECHNIQUE
D’ANTANANARIVO
*******************
Départements : Génie Mécanique Productique
Génie Electrique

Mémoire de fin d’étude en vue de l’obtention du diplôme d’ingénieur en Génie Industriel

ETUDE ET REALISATION D’UN CAPTEUR SOLAIRE A

ABSORBEUR ALVEOLAIRE

Soutenu par :
RAZAFINDRABE ANDRIANJAKA Manantsoa

Directeur de mémoire : Monsieur RANARIJAONA Jean Désiré,

Enseignant chercheur à l’ESPA
Président de jury: Monsieur RANDRIAMORASATA Josoa Albert,
Enseignant chercheur l’ESPA
Examinateurs : Monsieur JOELIHARITAHAKA Rabeatoandro,
Enseignant chercheur à l’ESPA
Monsieur RANDRIAMORASATA ANDRIANOELY Andy Ravaka,
Enseignant chercheur à l’ESPA
Date de soutenance : 24 Septembre 2011

Promotion 2010
 REMERCIEMENTS

Je remercie Dieu de m’avoir aidé durant toutes mes années d’études et pendant la
préparation de ce mémoire.
Je tiens aussi à remercier les personnes suivantes :
Monsieur ANDRIANARY Philippe, Directeur de l’Ecole Supérieure Polytechnique
d’Antananarivo, qui nous a permis de mener à terme nos études
Messieurs RAKOTOMANANA Charles Rodin et RAKOTONIAINA Solofo Hery,
respectivement Chef de département Génie Mécanique Productique et Chef de Département
Génie électrique, pour leur dévouement à notre égard pendant ces années de formations
Monsieur RANDRIAMORASATA Josoa Albert Enseignant chercheur à l’ESPA qui a
accepté de présider ce mémoire ;
Monsieur RANARIJAONA Jean Désiré, Enseignant chercheur à l’ESPA qui m’a
encadré pendant la réalisation de ce mémoire ;
Et à :
Monsieur JOELIHARITAHAKA Rabeatoandro, Enseignant chercheur à l’ESPA ;
Monsieur RANDRIAMORASATA ANDRIANOELY Andy Ravaka, Enseignant
chercheur à l’ESPA
qui ont bien voulu accepter de faire partie du jury.
Je remercie aussi toute ma famille qui m’a toujours soutenu et encouragé toute au
long de mes études universitaires.
J’exprime aussi mes sincères reconnaissances à mes collègues de la classe Génie
industriel, mes amis et à tous ceux qui de près ou de loin ont contribué à l’élaboration de ce
travail

RAZAFINDRABE ANDRIANJAKA Manantsoa I

 TABLE DES MATIERES
REMERCIEMENTS ................................................................................................................... I

TABLE DES MATIERES .........................................................................................................II

LISTE DES FIGURES ............................................................................................................. IV

LISTE DES TABLEAUX ......................................................................................................... V

LISTE DES ABREVIATIONS ET SYMBOLES .................................................................... VI

INTRODUCTION ...................................................................................................................... 1

PARTIE 1 GENERALITES……………………………………………………………………2

CHAPITRE 1 Les capteurs solaires ................................................................................... 3

1-1 Principe .................................................................................................................... 3

1-2 Types de capteurs ..................................................................................................... 3

1-3 Composants du capteur plan .................................................................................... 5

CHAPITRE 2 Le gisement solaire ...................................................................................... 6

2-1 Description du soleil[2] et [3] .................................................................................. 6

2-2 Mouvement de la Terre autour du soleil .................................................................. 6

2-3 Position du soleil ...................................................................................................... 7

2-4 Le rayonnement solaire .......................................................................................... 11

CHAPITRE 3 Notions de transfert thermiques ................................................................. 14

3-1 Conduction[4] ........................................................................................................ 14

3-2 Convection ............................................................................................................. 15

3-3 Rayonnement ......................................................................................................... 17

PARTIE 2 METHODOLOGIE………………………………………………………………20

CHAPITRE 4 Détermination de la structure alvéolaire .................................................... 21

4-1 La forme de la structure ......................................................................................... 21

4-2 Dimensions de la structure ..................................................................................... 22

CHAPITRE 5 Modélisation du capteur solaire ................................................................. 26

5-1 Hypothèses ............................................................................................................. 26

5-2 Principe de résolution ............................................................................................ 26

RAZAFINDRABE ANDRIANJAKA Manantsoa II

 5-3 Calcul des coefficients d’échanges thermiques [3] et [5] ...................................... 28

5-4 Calcul des facteurs de formes [4]........................................................................... 32

5-5 Flux absorbé par la vitre, l’absorbeur, et la structure alvéolaire ............................ 35

5-6 Mise en équation .................................................................................................... 39

PARTIE 3 REALISATION, MESURES ET INTERPRETATIONS………………………...41

CHAPITRE 6 Réalisation .................................................................................................. 42

6-1 Caractéristiques du capteur .................................................................................... 42

6-2 Coordonnées géographiques .................................................................................. 43

CHAPITRE 7 Mesures et interprétations .......................................................................... 44

7-1 Dispositif expérimental de mesure......................................................................... 44

7-2 Variation du rayonnement global, direct et diffus ................................................. 44

7-3 Comparaison entre les températures théoriques et réelles ..................................... 47

7-4 Comparaison entre les températures mesurées des différents absorbeurs ............. 53

7-5 Comparaison entre les températures des absorbeurs pour des distances vitre-
absorbeur différentes ........................................................................................................ 54

CHAPITRE 8 Etudes- économiques ................................................................................. 56

8-1 Devis estimatif du capteur alvéolaire ..................................................................... 56

8-2 Calcul de la VAN (valeur actuelle nette) et le TRI (taux de rentabilité interne) ... 58

8-3 Le POT (Pay out time) ou le temps de retour de l’investissement ........................ 59

CONCLUSION ........................................................................................................................ 60

BIBLIOGRAPHIE ...................................................................................................................... i

WEBOGRAPHIE ....................................................................................................................... ii

ANNEXE 1 ............................................................................................................................... iv

ANNEXE 2 ............................................................................................................................... vi

ANNEXE 3 ............................................................................................................................. viii

ANNEXE 4 ................................................................................................................................ x

RAZAFINDRABE ANDRIANJAKA Manantsoa III

 LISTE DES FIGURES
Figure 1 : Capteur à tubes sous vide .......................................................................................... 4
Figure 2 : Capteur plan [1] ........................................................................................................ 4
Figure 3 : Course du soleil ......................................................................................................... 6
Figure 4 : Coordonnées horaires ................................................................................................ 7
Figure 5 : Coordonnées azimutales [2] ...................................................................................... 8
Figure 6: Position d'un capteur ................................................................................................... 9
Figure 7 : Rayonnement solaire [3] .......................................................................................... 11
Figure 8:Phénomène de conduction ......................................................................................... 14
Figure 9:Phénomène de convection ......................................................................................... 15
Figure 10:Phénomène de rayonnement .................................................................................... 17
Figure 11: Réception du rayonnement ..................................................................................... 18
Figure 12 : Irradiation reçue par différents plaques ................................................................. 22
Figure 13: Numérotation des plaques d'une alvéole................................................................. 22
Figure 14 : Courbes représentant les flux en fonction du temps solaire vrai .......................... 23
Figure 15 : Courbes représentant les différents angles de repérage en fonction du temps
solaire vrai ................................................................................................................................ 24
Figure 16: Courbe représentant l'angle d'incidence du rayonnement en fonction du temps
solaire vrai ................................................................................................................................ 24
Figure 17:Les échanges thermiques dans une tranche de capteur ............................................ 28
Figure 18:Transmission du rayonnement ................................................................................. 36
Figure 19:Schéma électrique équivalent du capteur plan......................................................... 39
Figure 20: Le cadre .................................................................................................................. 43
Figure 21: Le capteur plan ....................................................................................................... 43
Figure 22 : Courbes des puissances solaires ............................................................................ 45
Figure 23: Courbes de variation des puissances absorbées des absorbeurs ............................. 46
Figure 24 : Courbes des températures de l’absorbeur1 ............................................................ 47
Figure 25: Courbe de température de l’absorbeur2 .................................................................. 48
Figure 26 : Courbe de température de l’absorbeur3 ................................................................. 49
Figure 27 : Courbe de température de la vitre1 ........................................................................ 50
Figure 28: Courbe de température de la vitre2 ......................................................................... 51
Figure 29 : Courbe de température de la vitre3 ........................................................................ 52
Figure 30: Courbes des températures des différents absorbeurs .............................................. 53
Figure 31 : Courbes de températures de l’absorbeur2.............................................................. 54

RAZAFINDRABE ANDRIANJAKA Manantsoa IV

 LISTE DES TABLEAUX

Tableau 1:Coefficient de trouble [3] ........................................................................................ 12

Tableau 2: Valeurs des albédos ................................................................................................ 13
Tableau 3 : Journées représentatives de la moyenne mensuelle .............................................. 13
Tableau 4:Analogie électrique ................................................................................................. 27
Tableau 5 : Caractéristiques thermo physique du capteur........................................................ 42
Tableau 6 : Caractéristiques géométriques du capteur ............................................................. 42
Tableau 7 : Caractéristique de la structure alvéolaire .............................................................. 42
Tableau 8 Devis du capteur ...................................................................................................... 56
Tableau 9 Calcul du cash-flow ................................................................................................. 57

RAZAFINDRABE ANDRIANJAKA Manantsoa V

 LISTE DES ABREVIATIONS ET SYMBOLES
a : arête de l’hexagone

A : coefficient de trouble

alb : albédo ou émissivité du sol

B : coefficient de trouble

C : coefficient de correction de la distance Terre-Soleil

Cn : capacité thermique de l’absorbeur

Cpi : capacité thermique de la plaque i du nid d’abeille

Cv : capacité thermique de la vitre

d : durée d’ensoleillement

Dh : diamètre hydraulique

env : distance entre la vitre et l’absorbeur

ET : équation du temps

F : facteur correctif égal à 6 pour une entrée avec contraction brusque

F12 : facteur de forme géométrique

g : gravité

Gz : nombre de GRAETZ

h : hauteur du soleil

H : hauteur du nid d’abeille

hca : coefficient d’échange thermique par conduction à travers la lame d’air

hcv : coefficient d’échange thermique par conduction à travers la vitre

hcvn : coefficient d’échange thermique par convection entre la vitre et l’absorbeur

hrnni : coefficient d’échange thermique par rayonnement entre l’absorbeur et une plaque i
du nid d’abeille

hrnv : coefficient d’échange thermique par rayonnement entre la vitre et l’absorbeur

hrv : coefficient d’échange thermique par rayonnement entre la vitre et la voûte céleste

hrvni : coefficient d’échange thermique par rayonnement entre la vitre et une plaque i du nid
d’abeille

RAZAFINDRABE ANDRIANJAKA Manantsoa VI

hv : coefficient d’échange thermique par convection entre la vitre et l’air ambiant

i : angle d’incidence du rayonnement

Io : constante solaire

j : numéro du jour de l’année

Kλ : coefficient d’extinction correspondant à la longueur d’onde

L : latitude

Lo : longitude

Mn : masse de l’absorbeur

Mpi : masse de la plaque i du nid d’abeilles

mv : masse de la vitre

n : indice de réfraction

NF : nombre de fuseau horaire

Nu : nombre de NUSSELT

Pdif : flux diffus

Pdir : flux direct

Pdirp : flux direct reçu par une surface

Pds : flux diffus venant du sol

Pdvc : flux émis par la voute céleste

Pr : nombre de PRANDTL

Ra : nombre de RAYLEIGH

Re : nombre de REYNOLDS

Ta : température ambiante

TL : temps local

Tn : température de l’absorbeur

Tpi : température de la plaque i de l’hexagone

TSV : Temps solaire vrai

TSVcouch : heure du couchée du soleil

RAZAFINDRABE ANDRIANJAKA Manantsoa VII

TSVlev : heure de levée du soleil

Tvc : température de la voûte céleste

Tve : température de la vitre externe

Tvi : température de la vitre interne

Vv : vitesse du vent

XJOUR : numéro du jour dans l’année

XMOIS : numéro du mois dans l’année

α : azimut du soleil

αa : facteur d’absorption de l’absorbeur

β : inclinaison du capteur

γ : azimut de surface

δ : déclinaison horaire

δa : épaisseur de la lame d’air statique

δv : épaisseur de la vitre

εv : émissivité de la voûte céleste

θ : angle du rayon réfracté

λ : longueur d’onde

λa : conductivité thermique de l’air

λv : conductivité thermique de la vitre

µe : viscosité dynamique du fluide

µp : viscosité dynamique calculée à la paroi

ρd : facteur de réflexion hémisphérique du vitrage

ζ : constante de BOLTZMAN

η : facteur de transmission de la vitre

ηd : facteur de transmission hémisphérique du vitrage

υ : viscosité cinématique

ω : angle horaire

RAZAFINDRABE ANDRIANJAKA Manantsoa VIII

 INTRODUCTION

L’énergie est un facteur clé pour le développement d’un pays. Actuellement, on se

penche plus sur les énergies renouvelables.

L’énergie solaire est gratuite, écologique. Le potentiel solaire de Madagascar est

estimé à 2000 W/m2. Afin de bien exploiter cette énergie, on doit utiliser une installation
solaire à haute performance.

Pour améliorer les performances d’un capteur, il existe plusieurs méthodes. On peut
recouvrir l’absorbeur par une couche sélective, ou bien augmenter le nombre de vitrage. On
peut aussi mettre une structure au dessus de l’absorbeur. C’est pourquoi on a choisi le thème
« ETUDE ET REALISATION D’UN CAPTEUR SOLAIRE A ABSORBEUR
ALVEOLAIRE »

Ce mémoire se divise en 3 grandes parties.

La première partie parlera des généralités sur les capteurs, le rayonnement solaire et
quelques notions de transferts thermiques.

La deuxième partie sera la méthodologie de la détermination de la structure alvéolaire et

la modélisation du capteur solaire

La dernière partie parlera de la réalisation et des interprétations des mesures obtenues et

des évaluations économiques

RAZAFINDRABE ANDRIANJAKA Manantsoa 1

RAZAFINDRABE ANDRIANJAKA Manantsoa 2
GENERALITES

CHAPITRE 1 Les capteurs solaires

1-1 Principe
Le capteur solaire thermique transforme le rayonnement solaire qu’il reçoit en énergie
calorifique utilisable. L’énergie est utilisée par un fluide caloporteur. Les fluides les plus
utilisés sont l’eau et l’air. L’air est généralement utilisé pour le séchage et l’eau pour la
climatisation.

L’énergie solaire est captée par la surface absorbante du capteur solaire. Le fluide
caloporteur absorbe et transmet une partie de l’énergie captée par l’absorbeur.

1-2 Types de capteurs

Il existe différents types de capteurs solaires :

1-2-1 Capteurs plans sans vitrage

Ils sont faits de plastique polymère noir et sont posés sur un support en bois ou sur un
toit. Ces capteurs absorbent bien l’énergie solaire mais les pertes thermiques vers
l’environnement augmentent rapidement avec la température extérieure. Ils sont la plupart du
temps utilisés pour des applications nécessitant une basse température (piscines, chaleur
industriel,…)

1-2-2 Capteurs plans avec vitrage

L’énergie solaire est emprisonnée dans le capteur à cause du vitrage (effet de serre).
Ils sont utilisés pour des applications à température modérée (chauffage de l’eau sanitaire,
chauffage de locaux,…)

1-2-3 Capteurs à tubes sous vides

Ces capteurs solaires sont composés d’un absorbeur revêtu d’une surface sélective et
enfermé sous vide dans un tube de verre. Le vide d’air dans les tubes réduit les pertes par
convection et par conduction. Pour le reste, ces capteurs fonctionnent comme un capteur plan.
Ils sont utilisés pour des applications nécessitant des températures moyennes ou hautes (eau
chaude domestique, chauffage de locaux et applications de chauffage industriel pour des
températures de 60°C à 80°C selon la température extérieure)

RAZAFINDRABE ANDRIANJAKA Manantsoa 3

GENERALITES

Figure 1 : Capteur à tubes sous vide

Source : Etude et conception d’un convertisseur solaire pour la production d’eau chaude
sanitaire [1]

Dans notre étude, on utilise un capteur plan avec vitrage.

Figure 2 : Capteur plan [1]

Source : Etude et conception d’un convertisseur solaire pour la production d’eau chaude
sanitaire [1]

RAZAFINDRABE ANDRIANJAKA Manantsoa 4

GENERALITES

1-3 Composants du capteur plan

1-3-1 L’absorbeur
L’absorbeur est l’élément principal du capteur. Le choix du matériau et le mode de
construction ont une grande influence sur la qualité du capteur. Il est soit en cuivre, soit en
acier, soit en aluminium. Les absorbeurs en cuivre ou en acier sont généralement revêtus
d’une couche sélective.

1-3-2 La couverture transparente

On utilise du verre sécurisé résistant aux chocs, auxcontraintes mécaniques et aux
chocs thermiques (brusque refroidissement). Le verre est généralement pauvre en oxyde de
fer, de ce fait transparent à l’entrée du rayonnement solaire mais opaque aux rayonnements
infrarouges. On distingue des capteurs à plusieurs vitrages.

1-3-3 L’isolation thermique

Il est constitué soit de laines minérales ou de matières synthétiques. Il doit résister aux
hautes températures à l’intérieur du capteur. L’isolant n’est pas utilisé pour les capteurs à
tubes évacués

1-3-4 Le cadre
Il est fait en aluminium ou en acier mais aussi en bois pour qu’il soit léger et facile à
déplacer.

RAZAFINDRABE ANDRIANJAKA Manantsoa 5

GENERALITES

CHAPITRE 2 Le gisement solaire

2-1 Description du soleil[2] et [3]

Le soleil est une sphère gazeuse formée d’hydrogène H et de hélium He et 70 autres
éléments plus lourds.

Le rayon solaire est composé de :

Lumière
Rayon γ, rayon X
Rayon ultra-violet (UV) λ = 0.05µm (longueur d’onde)
Rayons visibles (violet, bleu, vert, jaune, rouge) λ = 0.4 à 0.7µm
Rayon infrarouge IR λ> 2.5µm

2-2 Mouvement de la Terre autour du soleil

Figure 3 : Course du soleil

Source : Ingénierie des systèmes solaire applications à l’habitat [2]

RAZAFINDRABE ANDRIANJAKA Manantsoa 6

GENERALITES

2-3 Position du soleil

La position du soleil peut être repérée par les coordonnées azimutales et horaires.

2-3-1 Coordonnées horaires

Figure 4 : Coordonnées horaires

Déclinaison horaire δ

C’est l’angle entre la direction Terre-Soleil et le plan équatorial terrestre. Il varie de -

23°24’ au solstice d’hiver (21 Juin) et nul aux équinoxes.

(2.1)

j = n° du jour dans l’année

Angle horaire ω

C’est l’angle le long de l’équateur céleste et le plan de méridien du lieu (plan passant
par la direction sud) et celui contenant la direction des pôles célestes. Cet angle est compté
rétrograde vers l’ouest.

(2.2)

Avec TSV : Temps solaire vrai

Le TSV est défini à partir de la rotation de la Terre sur elle-même.

(2.3)

RAZAFINDRABE ANDRIANJAKA Manantsoa 7

GENERALITES

Avec TL : Temps local

NF : Nombre de fuseau horaire

ET : Equation du temps défini par :

ET=0.0072*cos(j) – 0.0523*cos(2j) - 0.0012*cos(3j) - 0.01229*sin(j) - (2.4)

0.1565*sin(2j)-0.0041*sin(3j)
j= n° du jour de l’année

Lo : Longitude, positive à l’Est

2-3-2 Coordonnées azimutales

Figure 5 : Coordonnées azimutales [2]

Hauteur du Soleil

C’est l’angle formé par le plan horizontal du lieu considéré et de la direction point
local-Soleil.

h= arcsin[sin(L)sin(δ)+cos(L)cos(δ)cos(ω)] (2.5)
Avec L : la Latitude du lieu est positive dans l’hémisphère Nord et négative dans
l’hémisphère Sud

δ : Déclinaison horaire

ω : Angle horaire

RAZAFINDRABE ANDRIANJAKA Manantsoa 8

GENERALITES

Azimut du soleil α

C’est l’angle entre la projection de l’axe OS sur le plan horizontal et le Sud. Il est
compté positivement vers l’Est et négativement vers l’Ouest.

(2.6)

2-3-3 Angle d’incidence du rayonnement solaire sur un plan incliné et

d’orientation quelconque
Une surface est repérée par:

L’inclinaison β :

C’est l’angle entre le plan horizontal et le plan considéré.

Azimut γ :

C’est l’angle entre la normale et le plan méridien.

Figure 6: Position d'un capteur

RAZAFINDRABE ANDRIANJAKA Manantsoa 9

GENERALITES

Le cosinus directeur l, m, n du vecteur u normal au plan est :

Le cosinus directeur X, Y, Z de la direction du soleil est défini par :

L’angle i entre la normale et le plan d’incidence du soleil sur le plan est :

(2.7)

α,h : azimut et hauteur du soleil

2-3-4 Problème d’heure

Heure de lever et de coucher du soleil

On a :

[Heure] (2.8)

[Heure] (2.9)

Durée de l’ensoleillement

Elle se détermine par la formule

(2.10)

RAZAFINDRABE ANDRIANJAKA Manantsoa 10

GENERALITES

2-4 Le rayonnement solaire

Le soleil émet plusieurs rayons dont une partie est captée par la Terre. Ces rayons sont
classés en :

Figure 7 : Rayonnement solaire [3]

2-4-1 Le flux direct

(2.11)

Io : constante solaire ; Io = 1367 W/m2

C : coefficient de correction de la distance Terre – Soleil ;

(2.12)

XMOIS= n° du mois dans l’année

XJOUR= n° du jour dans le mois

A et B étant des coefficients de trouble en fonction des phénomènes météorologiques

RAZAFINDRABE ANDRIANJAKA Manantsoa 11

GENERALITES

Ciel très pur Conditions normales Zone industrielle

A 0.87 0.88 0.91
B 0.17 0.26 0.43
Tableau 1:Coefficient de trouble [3]

Le flux direct reçu par la surface d’un capteur plan est :

[W/m2] (2.13)

i : angle d’incidence du rayonnement

2-4-2 Le flux diffus

Ce sont les rayons solaires absorbés puis diffusés par les différentes particules en suspension
dans l’atmosphère.

[W/m2] (2.14)

Il y a 2 sortes de flux diffus reçu par un capteur plan :

- Flux émis par la voute céleste

(2.15)

β : inclinaison du capteur

- Flux venant du sol

(2.16)

alb : albédo ou émissivité du sol

RAZAFINDRABE ANDRIANJAKA Manantsoa 12

GENERALITES

surface Réflectivité moyenne

Sol enneigé 0.7
Sol recouvert de feuilles mortes 0.3
Herbe verte 0.26
Forêt en automne ou champs dorés 0.26
Galet et pierres blanches 0.20
Herbe sèche 0.20
Sol argileux 0.17
Forêt en hiver 0.07
Plan d’eau 0.07
Tableau 2: Valeurs des albédos

L’irradiation globale varie très peu d’un jour à l’autre dans un même mois. Le calcul
peut donc se faire une fois par mois. La valeur de l’irradiation au cours d’un mois est donc à
peu près celle de la journée représentative du mois. Les journées représentatives sont
présentées dans le tableau suivant :

MOIS Jour du mois Jour de l’année

Janvier 17 17
Février 16 47
Mars 16 75
Avril 15 105
Mai 15 135
Juin 11 162
Juillet 17 198
Aout 16 228
Septembre 15 258
Octobre 15 288
Novembre 14 318
Décembre 10 344
Tableau 3 : Journées représentatives de la moyenne mensuelle

RAZAFINDRABE ANDRIANJAKA Manantsoa 13

GENERALITES

CHAPITRE 3 Notions de transfert thermiques

3-1 Conduction[4]
Il y a conduction si la chaleur se propage à travers le corps sans transfert de matière.

Figure 8:Phénomène de conduction

Pour un matériau homogène, une surface élémentaire S à travers laquelle passe de la
chaleur par conduction, on admet à l’instant t :

Loi de Fourrier

(3.1)

Ф en [W] : flux de chaleur transmis par conduction

x: variable d’espace dans la direction du flux

S : aire de la section de passage du flux de chaleur [m2]

Le coefficient d’échange par conduction est :

e : épaisseur du matériau en [m]

RAZAFINDRABE ANDRIANJAKA Manantsoa 14

GENERALITES

3-2 Convection
Il y a phénomène de convection s’il y a transmission de chaleur à chaque fois qu’un
fluide se déplace par rapport à des éléments fixes (parois, plaques, tubes) à une température
différente.

Figure 9:Phénomène de convection

L’écoulement d’un fluide est laminaire si l’écoulement s’effectue très lentement.

L’écoulement s’effectue par couches pratiquement indépendantes. La transmission de chaleur
par écoulement laminaire est faible due à la mauvaise conductibilité des fluides.

L’écoulement est turbulent quand les particules se déplacent de façon tourbillonnaire et

désordonnée.

Il y a 2 sortes de convection :

3-2-1 Convection naturelle

Tous les mouvements du fluide sont dus aux variations de son poids volumique
résultant des échanges de chaleur. Elle se présente quand un corps est placé dans un milieu
ambiant à une température différente de celle du corps.

3-2-2 Convection forcée

Les mouvements du fluide sont provoqués par une action mécanique extérieure (ex :
ventilateur)

RAZAFINDRABE ANDRIANJAKA Manantsoa 15

GENERALITES

Le flux de chaleur peut s’écrire :

(3.2)

h : coefficient d’échange par convection d’unité [w/m2 degré]

S : surface soumise à la convection [m2]

ΔT : différence de température entre le corps et le fluide

Par analogie à la conduction thermique, la résistance thermique par convection est

Les facteurs conditionnant l’écoulement d’un fluide sont englobés dans les relations
adimensionnelles :

Nombre de REYNOLDS

(3.3)

D : dimension caractéristique de la conduite : diamètre s’il s’agit d’une conduite circulaire,

diamètre hydraulique pour les autres cas : Dh =4.A/P avec A la section de la conduite, P le
périmètre mouillé.

ν: viscosité cinématique en [m2/s]

Ve : vitesse d’écoulement [m/s]

Nombre de NUSSELT

(3.4)

hc : coefficient local de convection thermique

Nombre de PRANDLT

(3.5)

Cp : capacité thermique massique en [J/ (Kg. °K)

µ : viscosité dynamique en [Kg/ms]

RAZAFINDRABE ANDRIANJAKA Manantsoa 16

GENERALITES

3-3 Rayonnement

Figure 10:Phénomène de rayonnement

Le transfert de chaleur par rayonnement est un phénomène électromagnétique.

L’énergie thermique se transmet au moyen des ondes électromagnétiques. Chaque corps
transforme une partie de son énergie en ondes électromagnétiques qui se déplacent dans
l’espace et frappe un autre corps où une partie de leur énergie est absorbée et convertie en
énergie.

Le rayonnement se propage dans l’espace libre à la vitesse de la lumière (vl=300.000 km/s)

Energie d'onde
Energie interne Energie interne
électromagnétique

La longueur d’onde du rayonnement est :

(3.6)

νr : fréquence de rayonnement, dépend de la nature de la source

La quantité totale du rayonnement qu’un corps émet par unité de surface et de temps
est le pouvoir émissif total.

RAZAFINDRABE ANDRIANJAKA Manantsoa 17

GENERALITES

Lorsque le rayonnement tombe sur un corps, une partie est réfléchie, une autre
absorbée, et une transmise.

Rayon Réfléchi
incident

.
Absorbé

Transmi
s
Figure 11: Réception du rayonnement

On note

- α : facteur d’absorption (dans le corps)

- ρ: facteur de réflexion (à la surface du corps)
- η : facteur de transmission (à travers le corps)

Avec

α, ρ, η dépendent du matériau (épaisseur, état de surface, …). Les corps qui ne transmettent
pas de rayonnement sont appelés opaques

-
-

Corps noir

C’est un corps qui absorbe totalement le rayonnement incident sur lui. Il n’y a pas de
réflexion et de transmission.

Corps gris (corps réel)

Ce sont les solides qui émettent une fraction constante de l’intensité du rayonnement
du corps noir à la même longueur d’onde et à la même température.

RAZAFINDRABE ANDRIANJAKA Manantsoa 18

GENERALITES

Loi de KIRCHOFF

Le rapport du pouvoir émissif E d’une surface et son facteur d’absorption α est le

même pour tout le corps.

Facteur d’émission

(3.6)

E : pouvoir émissif d’un corps gris

En : pouvoir émissif d’un corps noir

E<En

Pour un corps noir α=1

Loi de STEFAN-BOLTZMAN

Pour un corps noir,

Pour un corps gris,

Avec

Flux de chaleur échangé par rayonnement

La quantité de chaleur transmise entre 2 corps gris (ε1, T1) et (ε2, T2)

(3.7)

Avec

F12 : facteur de forme géométrique

Le coefficient d’échange par rayonnement est :

(3.8)

RAZAFINDRABE ANDRIANJAKA Manantsoa 19

RAZAFINDRABE ANDRIANJAKA Manantsoa 20
METHODOLOGIE

CHAPITRE 4 Détermination de la structure alvéolaire

Afin de bien exploiter le rayonnement solaire, la surface de captage (absorbeur) devrait

être orientée autant que possible perpendiculairement à ce rayonnement. Comme la position
du soleil par rapport à la terre varie, l’angle d’incidence du rayonnement sur le capteur varie
aussi. On choisit donc l’inclinaison du capteur la mieux adaptée au moment où le soleil brille
le plus, c’est-à-dire aux environs de midi (angle horaire w=0).

L’orientation plein sud permet le captage le plus efficace de l’énergie solaire. En faisant varier
l’inclinaison du capteur pour une orientation 0°, et pour w=0, l’irradiation est maximal pour
i=30°.

4-1 La forme de la structure

Afin d’augmenter l’irradiation absorbée par un capteur tout au long de la journée, on
va mettre une structure au dessus de l’absorbeur. L’irradiation reçue par l’absorbeur sera
augmenter grâce aux diverses réflexions des parois de la structure.

Une plaque orientée vers l’est et l’ouest capte plus l’irradiation pendant la matinée et la
soirée.

Pour fermer la structure en reliant les plaques orientées vers l’est et l’ouest, on peut utiliser
des plaques orientées vers le sud ou des plaques orientées de 30° et -30°.

RAZAFINDRABE ANDRIANJAKA Manantsoa 21

METHODOLOGIE

Figure 12 : Irradiation reçue par différents plaques

Grace aux formules (2.13) et (2.14), on remarque que la somme des irradiations des
plaques orientées de 30° et -30° est supérieure à celle de la plaque orientée vers le sud (0°).

La structure aura une forme d’hexagone.

4-2 Dimensions de la structure

On adoptera les numérotations suivantes :

Figure 13: Numérotation des plaques d'une alvéole

RAZAFINDRABE ANDRIANJAKA Manantsoa 22

METHODOLOGIE

On a :

(4.1)

(4.2)

L’irradiation reçue par la plaque4 est :

Figure 14 : Courbes représentant les flux en fonction du temps solaire vrai

Pour calculer les différents angles de repérage du soleil on a :

 Angle horaire

(4.3)

 Hauteur du soleil

(4.4)

 Azimut du soleil

(4.5)

RAZAFINDRABE ANDRIANJAKA Manantsoa 23

METHODOLOGIE

Figure 15 : Courbes représentant les différents angles de repérage en fonction du temps

solaire vrai
L’angle d’incidence du rayonnement sur le capteur est :

(4.6)

Figure 16: Courbe représentant l'angle d'incidence du rayonnement en fonction du

temps solaire vrai

RAZAFINDRABE ANDRIANJAKA Manantsoa 24

METHODOLOGIE

L’angle d’incidence minimal du rayonnement sur la plaque est de 20°.

Or

a : arête de l’hexagone

H : hauteur de la structure

On a une symétrie pour la plaque1 pour l’angle d’incidence minimal, donc :

(4.7)

RAZAFINDRABE ANDRIANJAKA Manantsoa 25

METHODOLOGIE

CHAPITRE 5 Modélisation du capteur solaire

5-1 Hypothèses
- Le vent souffle parallèlement aux faces du capteur
- La conduction thermique parallèle au verre est négligeable
- Les différents milieux solides ont une température uniforme dans le plan normal
de l’écoulement
- Dans chaque partie du système, l’écoulement est unidimensionnel
- La conduction entre la structure alvéolaire et l’absorbeur est négligeable par
rapport aux autres transferts de chaleur

5-2 Principe de résolution

Le système étudié est relativement complexe à cause du nombre élevé des paramètres.
Comme l’ensemble du système doit être étudié en régime transitoire, nous sommes en
présence d’un problème bidimensionnel dans lequel le temps t et la variable d’espace x joue
un rôle similaire. Le principe de résolution est un découpage spatial en un certain nombre de
section dx selon l’axe des x (direction de l’écoulement) et un découpage temporel fictif Δt
selon l’axe du temps t.

Pour un intervalle de temps Δt, on décrit tout le système de la première à la dernière

section spatiale. L’écriture et la résolution des bilans thermiques et massiques donnent la
répartition de température et des humidités relatives à chaque section. On fait de même pour
les autres intervalles de temps jusqu’à la fin du temps t. Le découpage temporel est le même
pour tout le système.

Pour mieux visualiser les transferts thermiques pour la mise en équation du système,
on fait une analogie électrique.

RAZAFINDRABE ANDRIANJAKA Manantsoa 26

METHODOLOGIE

Grandeurs thermiques Grandeurs électriques

Température T [K] Potentiel électrique V [V]
Flux thermique Ф [W] Courant électrique I [A]
Résistance thermique Rt [K/W] Résistance électrique R [Ω]
Capacité thermique Ct [J/W] Capacité électrique C [farad]
Densité de flux de chaleur θ [w/m2] Densité de courant i [A/m2]
Conductivité thermique λ [W/mK] Conductivité électrique γ [A/mV]
Tableau 4:Analogie électrique
Il ne reste plus qu’à faire la loi d’ohm sur les nœuds. La plupart des résistances Rij
dépendent des potentiels aux bornes.

Soit une section spatiale quelconque du système au temps t, on suppose qu’elle est une
entité indépendante. i est l’indice dans l’un des milieux représentés dans cette section.

Mi [kg] : masse

Ci [J.kg-1.K-1] chaleur spécifique

Ti [K] température

Le bilan thermique ou loi d’ohm aux nœuds i est :

(5.1)

Ji : ensemble des milieux J pour lesquels Tj est un potentiel connectés à Ti

Pour discrétiser les équations de la forme y’=f (xi,y), on utilisera la méthode explicite de
différences finies en arrière

(5.2)

RAZAFINDRABE ANDRIANJAKA Manantsoa 27

METHODOLOGIE

Visualisation schématique des échanges thermiques dans une tranche de capteur

Figure 17:Les échanges thermiques dans une tranche de capteur

Pv : flux absorbé par la vitre

Pn : flux absorbé par l’absorbeur

Pna : flux absorbé par le nid d’abeilles

5-3 Calcul des coefficients d’échanges thermiques [3] et [5]

5-3-1 Par convection

 Entre la vitre et l’extérieur

Cet échange est l’effet du vent et de l’air ambiant sur une plaque plane. Il se calcule par la
corrélation de HOTTEL et WOERTZ

(5.3)
Vv : vitesse du vent

 Entre l’absorbeur et le fluide caloporteur

C’est un écoulement forcé dans une conduite. Le diamètre hydraulique

S : surface droite de la conduite

P : périmètre de la section droite

RAZAFINDRABE ANDRIANJAKA Manantsoa 28

METHODOLOGIE

- Nombre de NUSSELT

(5.4)

- Nombre de REYNOLDS

(5.5)

Dh : diamètre hydraulique

υ: viscosité cinématique [m2/s]

- Nombre de PRANDTL

(5.6)

μe: viscosité dynamique du fluide [kg/ms]

- Nombre de GRAETZ

(5.7)

L : longueur caractéristique

- Nombre de GRASHOF

(5.8)

T1-T2 : écart de température entre 2 plans parallèles

 Ecoulement laminaire Re< 2100

Gz< 100

(5.9)

μa: viscosité dynamique du fluide

μp : viscosité dynamique du fluide à la paroi

RAZAFINDRABE ANDRIANJAKA Manantsoa 29
METHODOLOGIE

Gz> 100

(5.10)

 Zone de transition 2100 <Re< 10 000

HAUSSEN

(5.11)

 Ecoulement turbulent Re> 10 000

Pour (5.12)

Pour , il faut tenir compte des effets de bords

Nu0 : nombre de HAUSSEN pour

F : facteur correctif égal à 6 pour une entrée avec contraction brusque ou un demi-tour

Dans toutes ces corrélations, on néglige le terme

KAYES

Convection naturelle entre 2 plans parallèles

On utilise la corrélation suivante :

(5.13)

RAZAFINDRABE ANDRIANJAKA Manantsoa 30

METHODOLOGIE

Avec

s : inclinaison du capteur

 Entre la vitre et l’absorbeur

(5.14)

env : distance entre l’absorbeur et la vitre

5-3-2 Par conduction

 A travers la couverture

(5.15)

λv: conductivité thermique de la vitre

δv : épaisseur de la vitre

 A travers la lame d’air

(5.16)

5-3-3 Par rayonnement

 Entre la vitre et la voute céleste

(5.17)

εv : émissivité de la voute

ζ : constante de BOLTZMAN

Tvc : température de la voute céleste

 Entre la vitre et l’absorbeur

(5.18)

RAZAFINDRABE ANDRIANJAKA Manantsoa 31

METHODOLOGIE

 Entre la vitre et une plaque i du nid d’abeilles

(5.19)

 Entre une plaque i du nid d’abeilles et l’absorbeur

(5.20)

5-4 Calcul des facteurs de formes [4]

Considérons pour une alvéole, les numérotations des plaques suivant :

On va calculer les facteurs de formes entre la plaque1 et les différentes autres plaques.

Le facteur de forme géométrique Fik d’une surface Si par rapport à Sk est défini par la
relation [4]:

(5.21)

RAZAFINDRABE ANDRIANJAKA Manantsoa 32

METHODOLOGIE

 Entre la plaque1 et la plaque2

(5.22)

Avec θ=2*π/3 rad

 Entre la plaque1 et la plaque4

(5.23)

Avec

b : hauteur

a : arête de l’hexagone

c : distance entre 2 plaques parallèles

 Entre la plaque1 et la plaque7

Considérons la configuration suivante

RAZAFINDRABE ANDRIANJAKA Manantsoa 33

METHODOLOGIE

F1-7’ : facteur de forme entre la surface1 et la surface 7’

F1-7 » : facteur de forme entre la surface 1 et la surface 7 »

Avec

Avec

 Entre la plaque1 et la plaque5

Considérons la configuration suivante :

RAZAFINDRABE ANDRIANJAKA Manantsoa 34

METHODOLOGIE

(5.24)

Avec

5-5 Flux absorbé par la vitre, l’absorbeur, et la structure alvéolaire

5-5-1 Flux absorbé par la vitre

Soit n1 et n2 les indices de réfraction de deux différents milieux. Pour un rayon frappant
l’interface entre les deux milieux avec un angle d’incidence i, il sera partiellement réfléchi et
partiellement réfracté. Le rayon réfracté aura une incidence θ et on a :

RAZAFINDRABE ANDRIANJAKA Manantsoa 35

METHODOLOGIE

(5.25)

Figure 18:Transmission du rayonnement

La réflectivité est caractérisée par deux coefficients ρn et ρp qui sont les directions de
polarisation du rayon dans la direction perpendiculaire et parallèle au plan d’incidence.

(5.26)

(5.27)

Le coefficient de transmission correspondant θ et la longueur d’onde λ est :

(5.28)

Avec :

Kλ : Coefficient d’extinction correspondant à la longueur d’onde

x : distance parcourue par le rayon réfracté

Pour Kλ est considéré comme constant, on ne s’intéresse qu’à une marge limitée de λ, on a :

Pour une lame de verre épaisse, les coefficients de transmission, de réflexion et d’absorption
sont :

RAZAFINDRABE ANDRIANJAKA Manantsoa 36

METHODOLOGIE

(5.29)

(5.30)

Avec :

Pour déterminer ρ et η en fonction de l’angle d’incidence i, on a :

(5.31)

(5.32)

Et

(5.33)

5-5-2 Flux absorbé la structure alvéolaire

Le flux transmis η Pdirp par la vitre sera absorbé par la plaque du nid d’abeille αna.η.Pdirp.
Ensuite, ρna.η.Pdirp sera transmis à l’absorbeur.

(5.34)

5-5-3 Flux absorbé par l’absorbeur

Le flux direct Pdirp incident sur la surface de la vitre sera partiellement absorbé par la vitre
et ηPdirp sera transmis. αaηPdirpsera absorbé par l’absorbeur. ηρaPdirp sera réfléchi par
l’absorbeur et arrivera sur la vitre. ηρaρdPdirpsera réfléchi par la vitre et αaηρaρdPdirp sera
absorbé par l’absorbeur.

Après plusieurs réflexions, on a :

Or

(5.35)

RAZAFINDRABE ANDRIANJAKA Manantsoa 37

METHODOLOGIE

η : facteur de transmission de la vitre

ρd : facteur de réflexion hémisphérique du vitrage qui est souvent pris comme étant égale à
ρ (60°)

αa : facteur d’absorption de l’absorbeur pour les longueurs d’ondes λ correspondant au

rayonnement solaire

Pour un flux diffus Pdif, on a :

(5.36)

ηd : facteur de transmission hémisphérique du vitrage, que l’on peut aussi prendre pour égale à
η(60°)

Le flux total absorbé par l’absorbeur est la somme entre Ia et le flux absorbé réfléchi
par le nid d’abeille αa.ρna.η.Pdirp.

(5.37)

Pour le schéma électrique équivalent du capteur, on négligera la convection naturelle

entre l’air et la structure du nid d’abeilles et celle entre l’air et l’absorbeur (cf Annexes)

RAZAFINDRABE ANDRIANJAKA Manantsoa 38

METHODOLOGIE

Figure 19:Schéma électrique équivalent du capteur plan

5-6 Mise en équation

Au niveau du nid d’abeille, on va considérer chaque plaque dans une tranche du capteur.

Avec n : nombre total d’arêtes des hexagones sur la tranche de capteur.

RAZAFINDRABE ANDRIANJAKA Manantsoa 39

METHODOLOGIE

Dans la vitre transparente au nœud Tve

(5.38)

Au nœud Tvi

(5.39)

Au nœud Tni pour i=1 ,n

(5.40)

Au nœud Tn

(5.41)

Afin de déterminer les températures des différents éléments du capteur, on utilisera la

méthode de GAUSS-SEIDEL.(cf.Annexe)

RAZAFINDRABE ANDRIANJAKA Manantsoa 40

RAZAFINDRABE ANDRIANJAKA Manantsoa 41
REALISATION ET INTERPRETATIONS

CHAPITRE 6 Réalisation
L’objectif est de réaliser un prototype de capteur avec 3 sortes d’absorbeurs.

6-1 Caractéristiques du capteur

6-1-1 Caractéristique thermo physique

Elément de Matériaux Masse Chaleur Conductivité

construction volumique spécifique thermique
Kg/m3 J/Kg.°K W/m.°K
couverture Verre 2700 800 1.15

Absorbeur acier 7833 465 54

Boitier bois 535 272 0.144

Tableau 5 : Caractéristiques thermo physique du capteur

6-1-2 Caractéristiques géométriques

Elément de Longueur Largeur Epaisseur

construction mm mm mm
Couverture 240 1000 5

Absorbeur 240 1000 1

Tableau 6 : Caractéristiques géométriques du capteur

6-1-3 Caractéristique de l’alvéole

Arête Hauteur Rapport H/a

mm mm
Alvéole1 45 15 0.33

Alvéole2 60 15 0.25

Tableau 7 : Caractéristique de la structure alvéolaire

Pour fixer la structure alvéolaire sur les absorbeurs, on utilisera une pate thermique.

RAZAFINDRABE ANDRIANJAKA Manantsoa 42

REALISATION ET INTERPRETATIONS

Le cadre sera de la forme :

Vue d’ensemble vue éclatée

Figure 20: Le cadre

Pour l’isolation, on utilisera du contreplaqué.

6-2 Coordonnées géographiques

Notre étude est portée à Madagascar sur la région de Vontovorona de coordonnées
géographiques : Latitude : 18° Sud, Longitude : 44 Est. Pour les calculs on a choisi le mois
d’août.

L’absorbeur sans la structure alvéolaire sera notée absorbeur1. Celui avec la structure d’arête
45mm sera l’absorbeur2. Et celui avec l’arête 60mm sera l’absorbeur3.

Vue de face vue de haut

Figure 21: Le capteur plan

RAZAFINDRABE ANDRIANJAKA Manantsoa 43

REALISATION ET INTERPRETATIONS

CHAPITRE 7 Mesures et interprétations

7-1 Dispositif expérimental de mesure

Pour faciliter les différents calculs théoriques, on utilisera le logiciel MATLAB 7.0 qui est
très souple et facile à manipuler.

Pour les prises de mesures, on utilisera des thermocouples installés par serrage par vis sur
les absorbeurs et reliés à un HT10 pour la lecture des différentes températures au cours de la
journée.

Pour le pas de temps Δt, on prendra une valeur de 15 mn c’est-à-dire 900 s.

7-2 Variation du rayonnement global, direct et diffus

Pour les données géographiques précédentes, on a les courbes suivantes : pour une inclinaison
de 30° et d’orientation 0° le 10 Aout 2011 grâce aux formules (2.13) et (2.14)

RAZAFINDRABE ANDRIANJAKA Manantsoa 44

REALISATION ET INTERPRETATIONS

Figure 22 : Courbes des puissances solaires

Il y a une symétrie par rapport au midi solaire. Le soleil est alors à son zénith (hauteur
maximale) à midi donc l’irradiation absorbée est aussi à son maximum. Le flux diffus est très
faible par rapport au flux direct et le flux global.

RAZAFINDRABE ANDRIANJAKA Manantsoa 45

REALISATION ET INTERPRETATIONS

Variation de l’irradiation absorbée par les 3 différents absorbeurs

Figure 23: Courbes de variation des puissances absorbées des absorbeurs

Les puissances absorbées par les absorbeurs 2 et 3 d’après la formule (5.37) sont plus
élevées que celle de l’absorbeur1. Cela est due aux diverses réflexions de la structure
alvéolaire donc plus d’irradiation reçues. Le facteur d’absorption d’une tôle noire est de 0.95
tandis que celle d’une tôle galvanisée est de 0.65

L’irradiation reçue par l’absorbeur2 est supérieur à celle de l’absorbeur3 car

l’absorbeur 2 a plus de surface réfléchissante.

RAZAFINDRABE ANDRIANJAKA Manantsoa 46

REALISATION ET INTERPRETATIONS

7-3 Comparaison entre les températures théoriques et réelles

Pour l’absorbeur1

Figure 24 : Courbes des températures de l’absorbeur1

RAZAFINDRABE ANDRIANJAKA Manantsoa 47

REALISATION ET INTERPRETATIONS

Absorbeur 2

Figure 25: Courbe de température de l’absorbeur2

RAZAFINDRABE ANDRIANJAKA Manantsoa 48

REALISATION ET INTERPRETATIONS

Absorbeur3

Figure 26 : Courbe de température de l’absorbeur3

Les écarts entre les valeurs théoriques et les valeurs réelles pour les 3 absorbeurs ne
sont pas uniformes. Ceci est dû à la différence de montage entre les trois absorbeurs.
L’isolation des 3 absorbeurs n’est pas très exactement le même et le contact entre la structure
alvéolaire sur l’absorbeur n’est pas très parfait.

RAZAFINDRABE ANDRIANJAKA Manantsoa 49

REALISATION ET INTERPRETATIONS

Vitre1

Figure 27 : Courbe de température de la vitre1

RAZAFINDRABE ANDRIANJAKA Manantsoa 50

REALISATION ET INTERPRETATIONS

Vitre2

Figure 28: Courbe de température de la vitre2

RAZAFINDRABE ANDRIANJAKA Manantsoa 51

REALISATION ET INTERPRETATIONS

Vitre3

Figure 29 : Courbe de température de la vitre3

L’écart entre les valeurs théoriques et réelles des vitres sont à peu près les mêmes pour les 3
vitres. Elles sont presque soumises aux mêmes conditions (vitesse de l’air extérieur,
température ambiante,…)

Les valeurs mesurées sont plus faibles que celles calculées théoriquement. Plusieurs
facteurs sont les responsables de cet écart. Il y a le passage de certains nuages pendant les
mesures créant la diminution de l’ensoleillement, les pertes dans l’isolation, la température
ambiante de l’air varie aussi avec la journée à cause de la vitesse de l’air extérieur.

RAZAFINDRABE ANDRIANJAKA Manantsoa 52

REALISATION ET INTERPRETATIONS

7-4 Comparaison entre les températures mesurées des différents absorbeurs

Figure 30: Courbes des températures des différents absorbeurs

La température de l’absorbeur2 est plus élevée que celles des 2 autres tout au long de
la journée. Les surfaces réfléchissantes sont plus nombreuses et la conduction entre les
plaques verticales et la plaque horizontale est aussi plus élevé même si elle est relativement
faible par rapport au rayonnement entre les plaques.

RAZAFINDRABE ANDRIANJAKA Manantsoa 53

REALISATION ET INTERPRETATIONS

7-5 Comparaison entre les températures des absorbeurs pour des distances vitre-
absorbeur différentes
On utilisera 2 distances différentes. e1=25mm et e2=50mm.

Pour l’absorbeur2

Figure 31 : Courbes de températures de l’absorbeur2

Les températures des absorbeurs ayant une distance vitre-absorbeur de 25 [mm] sont
plus élevées que celles ayant la distance de 50 [mm]. La masse d’air entre la vitre et
l’absorbeur est plus élevé pour ces derniers donc il y une différence entre les conductions
d’air. La distance parcourue par le rayon solaire aussi devient plus large donc il y a des pertes.

RAZAFINDRABE ANDRIANJAKA Manantsoa 54

REALISATION ET INTERPRETATIONS

L’irradiation journalière reçue par un capteur avec un simple absorbeur est aux
environs de 996 [W], d’après la figure(23) tandis que celle reçue par un capteur avec un
absorbeur alvéolaire d’arête d’hexagone 45mm est 1054.3 [W]. Pour l’absorbeur3 dont l’arête
est de 60mm, l’irradiation absorbée est 1020.2 [W]

Pour l’absorbeur2, on a un gain d’énergie journalière de 6.5% par rapport au capteur

simple, tandis que pour l’absorbeur3, ce gain s’élève à 1,024 %

L’absorbeur2 est plus efficace que les 2 autres absorbeurs avec une épaisseur de lame
d’air de 25mm.

RAZAFINDRABE ANDRIANJAKA Manantsoa 55

REALISATION ET INTERPRETATIONS

CHAPITRE 8 Etudes- économiques

8-1 Devis estimatif du capteur alvéolaire

Les investissements initiaux pour la réalisation d’un capteur thermique solaire avec absorbeur
alvéolaire sont représentés dans le tableau suivant :

matériels Prix [Ar]

tôle noire 10 000

tôle galvanisée 5 000

vitre 10 000

pate thermique 5 000

cadre 10 000

peinture noire mat 4 000

Frais de recherche 20 000

total 64 000

Tableau 8 Devis du capteur

Les couts d’installation et de transport sont de 10%

L’investissement total est donc de 70 400 Ar

Le gain énergétique est de 1056 [Wh]. Supposons que le prix du Wh est de 0.2 Ar.

Le bénéfice journalier est de 1056*0.2=211.3 Ar

Donc le bénéfice annuel.est de 211*365= 77088Ar

On a comme amortissement :

En bref, nous pouvons établir l’échéancier suivant :

RAZAFINDRABE ANDRIANJAKA Manantsoa 56

REALISATION ET INTERPRETATIONS

Bilan

Durée de vie du projet : 5 ans

Investissement total : 70 400 Ar

Valeur de la revente à la 5 ième année est nulle

Exploitation

Bénéfice annuel : 77088 Ar

Charges=bénéfice annuel*20% = 23126.4 Ar

L’amortissement linéaire : 2346.66 Ar

Impôts sur les bénéfices : 0% car le projet est personnel

On a alors le tableau suivant :

années 1 2 5

recettes 70400 70400 70400

charges 23126,4 23126,4 23126,4

dotations aux amortissements 2346,66667 2346,66667 2346,66667

résultat avant amortissement 44926,9333 44926,9333 44926,9333

résultat après amortissement 44926,9333 44926,9333 44926,9333

IBS 0 0 0

Cash flow 47273,6 47273,6 47273,6

Tableau 9 : Calcul du Cash-Flow

IBS : impôts sur les bénéfices

RAZAFINDRABE ANDRIANJAKA Manantsoa 57

REALISATION ET INTERPRETATIONS

8-2 Calcul de la VAN (valeur actuelle nette) et le TRI (taux de rentabilité interne)
La VAN est la somme des revenus annuels actualisés, calculée pendant toute la durée
de vie de l’investissement sur le montant du capital investi (CA)

(8.1)

CA : capital investi

Le TRI est égal à la valeur du taux d’actualisation ‘i’ pour lequel la VAN s’annule.

Pour un taux d’actualisation i de 15% on a :

VAN à i=15%

En posant

(8.2)

(8.3)

VAN=88069.4 Ar

VAN à i=60%

VAN à i=65%

La valeur actuelle nette s’annule pour un taux d’actualisation entre 60% et 65%

Par interpolation on a :

RAZAFINDRABE ANDRIANJAKA Manantsoa 58

REALISATION ET INTERPRETATIONS

8-3 Le POT (Pay out time) ou le temps de retour de l’investissement

C’est le temps correspondant à la durée d’exploitation nécessaire pour récupérer le montant de
l’investissement.

(8.4)

(8.5)

POT=1.48 ans

En conclusion, le temps de retour de l’investissement est de 1.48ans pour une durée de vie de
5 ans pour le prototype. Le TRI atteint 60.9% (> 15%), Le capteur solaire alvéolaire est alors
rentable.

RAZAFINDRABE ANDRIANJAKA Manantsoa 59

 CONCLUSION

L’exploitation de l’énergie solaire est toujours rentable car elle est renouvelable,
gratuite et non polluante. Une amélioration des capteurs est nécessaire pour bien utiliser cette
énergie.

Des études théoriques effectuées dans le cadre de cette recherche ont permis de
montrer que les performances d’un absorbeur alvéolaire sont meilleurs que ceux d’un capteur
plan. La réalisation de trois types d’absorbeurs et les mesures effectuées sur ceux-ci
permettent d’estimer les performances respectives grâce à la différence de température. Les
résultats confirment davantage que le capteur à nid d’abeille est meilleur qu’un capteur simple
et que le gain d’énergie journalière est de 6% pour un capteur alvéolaire dont les arêtes sont
de 45mm et la hauteur de 15mm. Ces mesures montrent également que la quantité d’énergie
reçue est meilleure si la distance qui sépare la vitre et l’absorbeur est de 25mm. Afin de se
rapprocher davantage des performances théoriques, dans la pratique il convient d’améliorer
l’isolation des absorbeurs et le plaquage des alvéoles sur l’absorbeur.

Par ailleurs, l’évaluation économique préliminaire justifie la rentabilité de l’utilisation d’un

capteur à nid d’abeille, étant donné que l’investissement est de la VAN est de 88068 Ar et le
TRI est 60.9 Ar.

Cette étude s’est limitée seulement à l’absorbeur, mais peut s’appliquer favorablement
à l’étude d’un système de séchage ou de chauffe – eau, avec un fluide caloporteur.

Nous espérons que ce travail apportera son aide pour compléter d’autres recherches.

RAZAFINDRABE ANDRIANJAKA Manantsoa 60

 BIBLIOGRAPHIE

[1] ABDELKARIM BOURAS, « Etude et conception d’un convertisseur solaire pour la

production d’eau chaude sanitaire », mémoire en physique énergétique 2007

[2] ,A.-A.SFE et G.GUARRACINO, « Ingénierie des systèmes solaires applications à

l’habitat », Paris 1981

[3]Harinjaka RASAMIMANANA, « Modélisation d’un système froid solaire par adsorption

solide utilisant le couple frigorigène CaCl2-CH3NH2 pour la conservation des produits
pharmaceutiques dans des zones non électrifié », mémoire génie industriel promotion 2003

[4] Yves JANNOT, « Transfert thermique », école des mines Nancy, 2009

[5] RAKOTOVAO José Denis, « Analyse informatique d’un système de séchage solaire
indirect en couches minces pour fruits et légumes », thèse

[6]Boumediene TOUATI, « Etude théorique et expérimentale du séchage solaire des feuilles

de la menthe verte (Menthaviridis) », Thèse ,2008

[7] Magdeleine HUETZ-AUBERT, « Rayonnement thermique des matériaux opaques »,

Technique de l’Ingénieur

[8]RAJOELIMANANA Andryfitiavana, « Estimation de l’irradiation solaire diffuse reçue sur

une surface inclinée à Antananarivo », DEA Physique 2006

[9]J.Brau, « Transfert de chaleur et de masse » , INSA de Lyon, 2006

[10] RASOLONJATOVO NomenjanaharyBerthin, « Modélisation d’un chauffe-eau solaire,

Système toiture EST/OEUST, Mémoire Génie industriel 2006

RAZAFINDRABE ANDRIANJAKA Manantsoa i

 WEBOGRAPHIE
1. http://www.Ledifice.com.materielsolairethermique

2. http://www.cansia.ca

3. http://www.tecsol.fr

4 http:/www.these.recherches.gov.mg

RAZAFINDRABE ANDRIANJAKA Manantsoa ii

RAZAFINDRABE ANDRIANJAKA Manantsoa iii
 ANNEXE 1
Méthode de GAUSS-SEIDEL

Si le nombre n d’équation du système linéaire est élevé, on utilise la méthode

itérative.

La méthode classique est du type linéaire

Principe

Considérons le système A.x=b

On décompose la matrice A sous la forme A=M-N

Avec ;

Supposons que A est inversible aii≠ 0 ; i=1,…, n donc M est aussi inversible

Ou

Calcul des composantes

RAZAFINDRABE ANDRIANJAKA Manantsoa iv

 On peut ainsi calculer la ièmecomposante de x(m+1) en fonction des i-1 premières
composantes de x(m) sans être obligé de garder en mémoire les autres composantes de x(m).

La méthode de GAUSS - SEIDEL converge si :

Ou

RAZAFINDRABE ANDRIANJAKA Manantsoa v

 ANNEXE 2
Organigramme pour la modélisation du capteur

Début

Données caractéristiques

Initialisation des températures

Boucle du temps J

Calcul du rayonnement
Réinitialisation des températures

Boucle d’espace I

Calcul des surfaces d’échanges thermiques

Propriétés de l’air
Nombres adimensionnels
GAUSS-SEIDEL
 Coefficients d’échanges thermiques
 Températures respectives des éléments

NON

Tj(i) - Tj*(i) <ε

OUI

Température de l’absorbeur I=I+1 J=J+1

OUI
I ≤ NTR
NON
Température de l’absorbeur

OUI
J ≤ Jmax
NON

Courbes des températures

Fin

RAZAFINDRABE ANDRIANJAKA Manantsoa vi

NTR : nombre de tranches

Tj(i) : température de la tranche I au temps t

Tj*(i) : température de la tranche I au temps t-Δt

RAZAFINDRABE ANDRIANJAKA Manantsoa vii

 ANNEXE 3

Tableaux des mesures des températures en [°C]

T1 T2 T3 vitre1 vitre2 vitre3

8h50 25 47 44,6 23 24,5 23
9h05 32,5 65,4 41 24,2 24 25
9h20 45 71,4 52,6 24,6 24,4 25,2
9h35 61 78,2 45,2 28 30,7 25,3
9h50 67,6 80,2 68,2 26,4 25 31,7
10h05 64,9 85,7 79 26,1 26,3 26,9
10h20 68,3 90 66 27,6 27,7 27
10h35 70,5 91,4 84,3 27,8 28 27
10h50 74,6 94,4 81,4 31 31,3 30,4
11h05 70,4 101,1 63,2 28 28,8 27,8
11h20 74,2 100,5 89 28 27,8 28
11h35 77,3 102,3 101,9 29 28,2 27,9
11h50 80 97 98 28,5 28 28,3
12h05 60,3 102,9 91,8 32,3 30 31,3
12h20 59 103,2 95,1 32,5 29,7 29,9
12h35 71 104,5 99,2 33,3 29,5 29,5
12h50 60,1 102,2 97,7 25,2 27,4 28,1
13h05 54,5 98,9 85,4 26,3 28,1 27
13h20 70 95,8 95,9 26 27,9 26,7
13h35 77,5 92,3 96 25,8 27,5 26,5
13h50 61,7 93,9 94 25,5 27 27
14h05 65 89,5 84,6 25 26,7 27,1
14h20 67,9 86 77 24,7 26,5 26,8
14h35 49,4 81 82,2 24,2 26,3 26,5
14h50 55,6 81 79,5 24,3 26 26,2
15h05 44,5 75 76,7 24 25,9 26,4
15h20 43,3 70,5 68,9 24,1 25,8 26
15h35 46,9 63,9 65 24,2 25,5 25,9
15h50 41,9 56,9 61 24 25 24,8
16h05 41,3 53 51,8 23,9 24,2 24
16h20 38,7 43 41 23,1 23,9 23,9
16h35 35,8 36 35,1 23 23,6 23,5
16h50 33,3 34 33 23 23 23

RAZAFINDRABE ANDRIANJAKA Manantsoa viii

E=50mm

T1 T2 T3
8h50 30,1 33,4 24,3
9h05 35,2 49,3 37,5
9h20 43,3 54 46,2
9h35 46,2 59,6 52,1
9h50 53,7 64,1 59,4
10h05 60 73,4 66
10h20 68,2 84,5 69,8
10h35 70,5 84,6 71,3
10h50 72,1 84,9 74,7
11h05 76,3 84,9 79,5
11h20 64,7 72,9 75,6
11h35 65 84,5 74,5
11h50 67,8 90,6 86,5
12h05 69,7 96,8 91
12h20 65 80 81
12h35 64,4 88,2 83,4
12h50 64,8 76,9 77,4
13h05 69,4 81,3 73,3
13h20 60,6 89,4 85,5
13h35 64,2 80,8 80
13h50 62,6 76,1 75,6
14h05 56,8 72,9 74,2
14h20 53,1 72 69,6
14h35 52,5 71,7 66,4
14h50 51,6 72 67,1
15h05 44,2 65,2 52,8
15h20 44,4 63,9 50,4
15h35 44,6 60,8 48,1
15h50 40 48 40
16h05 38,1 54,2 44,4
16h20 32,7 42 37
16h35 30,2 36,5 34,3
16h50 25 26,1 25,3

RAZAFINDRABE ANDRIANJAKA Manantsoa ix

 ANNEXE 4
Ordre de grandeur du coefficient de transmission de chaleur par convection

configuration h(W m-2.°C-1)

Convection naturelle
Dans un gaz 2-10 (546 W.m-2°K-1)
Dans un liquide 100-1000
Convection forcée
Avec un gaz 10-200
Avec un liquide 100-5000
Ebullition de l’eau
Dans un récipient 2500-35000
En écoulement dans un tube 5000-100000
Condensation de l’eau sous 1atm
Sur une surface verticale 1000-11000
A l’extérieur des tubes hotrizentaux 10000-25000

Valeur des coefficients d’échange par rayonnement entre la vitre et la plaque1 du nid’ abeilles
et entre la plaque1 du nid d’abeilles et l’absorbeur

hrvp1 hrp1n
(W.m-2.°K) (W.m-2.°K)
6h30 0.0089 0.1539
6h45 0.0089 0.1542
7h00 0.0095 0.1647
7h15 0.0112 0.1939
7h30 0.0113 0.1956
7h45 0.0102 0.1932
8h00 0.0138 0.2404
8h15 0.0116 0.2395
8h30 0.0140 0.2432
8h45 0.0117 0.2409
9h00 0.0117 0.2421
9h15 0.0165 0.2861
9h30 0.0131 0.2869
9h45 0.0131 0.2882

RAZAFINDRABE ANDRIANJAKA Manantsoa x

10h00 0.0131 0.2894
10h15 0.0167 0.2903
10h30 0.0132 0.2913
10h45 0.0132 0.2922
11h00 0.0169 0.2934
11h15 0.0170 0.2955
11h30 0.0171 0.2974
11h45 0.0172 0.2995
12h00 0.0174 0.3016
12h15 0.0174 0.3024
12h30 0.0174 0.3029
12h45 0.0174 0.3031
13h00 0.0174 0.3031
13h15 0.0135 0.3045
13h30 0.0135 0.3041
13h45 0.0134 0.3035
14h00 0.0134 0.3027
14h15 0.0134 0.3016
14h30 0.0134 0.3001
14h45 0.0145 0.2517
15h00 0.0145 0.2517
15h15 0.0119 0.2515
15h30 0.0119 0.2508
15h45 0.0130 0.2264
16h00 0.0111 0.2259
16h15 0.0117 0.2029
16h30 0.0116 0.2017
16h45 0.0103 0,2013

RAZAFINDRABE ANDRIANJAKA Manantsoa xi

Auteur : RAZAFINDRABE ANRIANJAKA Manantsoa

Adresse : Lot I H20 Isotry Antananarivo 101

Télephone:0341938294/0330947727

Email:andrianjakamanantsoa@yahoo.fr

Thème de mémoire:

« ETUDE ET REALISATION D’UN CAPTEUR SOLAIRE A ABSORBEUR

ALVEOLAIRE »

Encadreur : Monsieur RANARIJAONA Jean Désiré

Nombre de pages : 59

Nombre de figures : 31

Nombre de tableaux : 9

Résumé

Le présent mémoire a pour but d’améliorer les apports énergétiques d’un

capteur grâce à la mise en place d’une structure alvéolaire ayant une forme d’hexagone. Pour
pouvoir mener à bien cette étude, nous avons estimé le rayonnement solaire incident sur le
capteur, puis on a sorti le modèle mathématique régissant le capteur. On a réalisé ensuite un
prototype avec différents types d’absorbeurs afin de comparer les valeurs théoriques et réelles.

Les résultats obtenus ont bien démontré que l’’absorbeur ayant une structure alvéolaire est
plus avantageux du point de vue énergétique.et rentable économiquement.

Abstract

The purpose of the present report is to improve the energy contributions of a sensor
thanks to the installation of a honeycomb structure having hexagon form. To be able to
conclude this study, we considered the incidental radiation solar on the sensor then one left
the mathematical model governing the sensor. One then produced a prototype with various
types of absorbers in order to compare the theoretical and real values. The results obtained
showed well that the ' absorber with the honeycomb structure is more advantageous from the
energy and profitable economically.

Mots clés : irradiation, capteur, alvéole, température;rayonnement