Introduction
Un bilan énergétique est une évaluation systématique et détaillée des flux d'énergie entrant et
sortant d'un système donné. Que ce soit au niveau d'un bâtiment, d'une ville, d'un pays ou même à
l'échelle mondiale, un bilan énergétique permet de quantifier et de caractériser les différentes sources
d'énergie utilisées, les conversions d'énergie qui se produisent et les pertes énergétiques. Ce processus
d'évaluation offre une compréhension approfondie de la consommation d'énergie d'un système, de son
efficacité énergétique et de son impact sur l'environnement. Il permet d'identifier les secteurs les plus
énergivores, les sources d'énergie les plus utilisées, ainsi que les opportunités d'amélioration et de
transition vers des sources d'énergie plus durables
Figure III.3 : Schéma de principe capteur solaire plan à air.
III.1 Bilan thermique global de capteur plan a air en régime transitoire :
Le but est de prévoir les températures de chaque couche du capteur (vitre, absorbeur, le
fluide caloporteur, isolant.). Pour cela on utilise le principe de conservation de l’énergie pour
chaque élément du système qui peut donc être présenté par une équation de bilan d’énergie qui
s’écrit :
𝑙𝑎 𝑣𝑟i𝑎𝑡i𝑜𝑛 𝑑𝑒 𝑙′é𝑛𝑒𝑟𝑔i𝑒 𝑙′é𝑛𝑒𝑟𝑔i𝑒 𝑝𝑒𝑟𝑑𝑢
reçue
𝑙 é𝑛𝑒𝑟𝑔i𝑒
′
={ −
𝑡𝑒𝑟𝑚𝑒 𝑎𝑐𝑐𝑢𝑚𝑢𝑙𝑎𝑡i𝑜𝑛 𝑝𝑎𝑟 𝑒𝑙𝑒𝑚𝑒𝑛𝑡 (𝑒𝑛𝑡𝑟𝑎𝑛𝑡)
i𝑛𝑡𝑒𝑟𝑛𝑒 𝑝𝑎𝑟 𝑒𝑙𝑒𝑚𝑒𝑛𝑡 (𝑠𝑜𝑟𝑡𝑎𝑛𝑡
Ou encore sous forme mathématique, en considérant une section quelconque du
système à l’instant t. le bilan au nœud « i » s’écrit :
CHAPITRE III : MODELISATION ET SIMULATION D’UN CAPTEUR SOLAIRE PLAN Á AIR
𝑑𝑇i (III. 22)
𝑀𝐶 = ∑𝜑 −∑𝜑
i i 𝑑𝑡 i i
𝑒 𝑠
Oú𝑀iet Ci Sont respectivement la masse et la chaleur spécifique du composant indice "i".
𝑀𝐶 𝑑𝑇iest la puissance thermique stockée par le composant indice i.
𝑑𝑡
∑𝑒 𝜑i est la somme des puissances absorbé parle composant d'indice i par trois modes
d'échange thermique (conduction; convection et rayonnement).
∑𝑠 𝜑i est la somme des énergies perdues par le composant d'indice i par les trois modes
d'échange thermique.
𝝋𝒓𝒄𝒊𝒆𝒍−𝒗𝒆
Vitre
𝝋𝒄𝒅𝒗
Sortie
𝝋𝒄𝒗𝒗𝒆−𝒂𝒎 l'air
𝒗
𝝋𝒓𝒗𝒊−𝒂𝒃𝒔
Absorbeur
𝝋𝒄𝒗𝒇−𝒂𝒃𝒔
𝝋𝒄𝒗𝒗𝒊−𝒇 Isolant
𝝋𝒔𝒂𝒃𝒔
𝝋𝒄𝒗𝒊𝒔𝒐𝒍−𝒂𝒎
Entrée 𝝋𝒓𝒔𝒐𝒍−𝒊𝒔𝒐𝒍
l'air
Figure III.7 : Echanges thermiques dans un capteur solaire plan à air
: Nœud représentant la température du vitrage coté extérieure.
: Nœud représentant la température du vitrage coté intérieure.
: Nœud représentant la température du fluide (air).
: Nœud représentant la température de la couche l’absorbeur.
: Nœud représentant la température de la couche isolant.
2
III.1.1 Expressions analytiques de bilan énergétique :
Le bilan énergétique général est donné par les équations suivantes:
Au niveau du vitrage (coté extérieur) :
𝑑𝑇𝑣𝑒
𝑀 𝐶 =𝜑 −𝜑 −𝜑 −𝜑 −𝜑 (III. 23)
𝑣 𝑃𝑣 𝑑𝑡 𝑠𝑣 𝑟𝑣𝑒−𝑐i𝑒𝑙 𝑐𝑣𝑣𝑒−𝑎𝑚 𝑟𝑣𝑒−𝑠𝑜𝑙 𝑐𝑑𝑣𝑒−𝑣i
Flux solaire absorbé par la vitre:
𝜑𝑠𝑣 = 𝑎𝑣𝑒. 𝐺. 𝑆 (III. 24)
Avec:
𝑎𝑣 : Coefficient d’absorption de vitre.
G : éclairement (densité de flux) solaire incident sur le capteur (w/m²).
S : surface de la captation (m).
Finalement cette équation devient:
𝑀𝑣𝐶𝑃𝑣 𝑑𝑇𝑣𝑒 = 𝑆[𝑎𝑣𝐺−𝑟 (𝑇𝑣𝑒−𝑇𝑐i𝑒𝑙)− (𝑇𝑣𝑒−𝑇𝑎𝑚)−𝑟 (𝑇𝑣𝑒 −
𝑑𝑡 𝑣𝑒−𝑐i𝑒𝑙 𝑐𝑣𝑣𝑒−𝑎𝑚 𝑣𝑒−𝑠𝑜𝑙
𝐾 .
𝑇𝑠𝑜𝑙) − 𝑣 (𝑇𝑣𝑒 − 𝑇𝑣i)] (III. 25)
𝑒𝑣
Au niveau du vitrage (coté intérieur):
𝑑𝑇𝑣i
𝑀𝐶 =𝜑 −𝜑 −𝜑 (III. 26)
𝑣 𝑃𝑣 𝑑𝑡 𝑐𝑑𝑣𝑒−𝑣i 𝑟𝑣i−𝑎𝑏𝑠 𝑐𝑣𝑣i−ƒ
Donc cette équation devient:
𝑑𝑇𝑣i 𝐾 𝑆
𝑀 𝐶 = 𝑣 (𝑇 − 𝑇 ) − S(𝑇 − 𝑇 ) − 𝑆(𝑇 − 𝑇 ) (III. 27)
𝑣 𝑃𝑣 𝑑𝑡 𝑒 𝑣𝑒 𝑣i 𝑟_𝑣i−𝑎𝑏𝑠𝑣i 𝑎𝑏𝑠 𝑐𝑣 _𝑣i−𝑓 𝑣i 𝑣𝑓
𝑣
Pour le fluide caloporteur (air) : [15]
𝑑𝑇
𝑚 𝐶 =𝜑 −𝜑 (III. 28)
ƒ 𝑃ƒ 𝑑𝑥 𝑐𝑣𝑣i−ƒ 𝑐𝑣ƒ−𝑎𝑏𝑠
Avec:
𝑚𝑓: Débit massique de l’air (kg/s)
𝑑𝑜𝑛𝑐 :
𝑑𝑇
𝑚 𝐶 = 𝑆(𝑇 − 𝑇 )− S(𝑇 − 𝑇 ) (III. 29)
𝑓 𝑃𝑓 𝑑𝗑 𝑐𝑣𝑣i−𝑓 𝑣i 𝑓 𝑐𝑣𝑓−𝑎𝑏𝑠 𝑓 𝑎𝑏𝑠
Au niveau de l’absorbeur (couche sélective) :
3
(III. 30)
𝑑𝑇𝑎𝑏𝑠
𝑀 𝐶 =𝜑 +𝜑 −𝜑 −𝜑
𝑎𝑏𝑠 𝑝𝑎𝑏𝑠 𝑑𝑡 𝑠𝑎𝑏𝑠 𝑟𝑣i−𝑎𝑏𝑠 𝑐𝑣ƒ−𝑎𝑏𝑠
𝑐𝑑𝑎𝑏𝑠
Flux solaire absorbé par absorbeur s’écrit :
𝜑𝑠𝑎𝑏𝑠 = 𝑟𝑣. 𝑎𝑎𝑏𝑠. 𝐺. 𝑆
Avec :
𝑟𝑣 : Coefficient de transmission de la vitre par rapport au rayonnement solaire.
𝑎𝑎𝑏𝑠:Coefficient d’absorption de la paroi absorbante par rapport au rayonnement solaire.
S: Surface de la paroi absorbante (m²)
Cette équation devient :
𝑀𝑎𝑏𝑠𝐶𝑝 𝑑𝑇𝑎𝑏𝑠 = 𝑟𝑣𝑎𝑎𝑏𝑠𝐺𝑆+𝑟 𝑆(𝑇𝑣i − 𝑇𝑎𝑏𝑠) − 𝑐𝑣 𝑆(𝑇𝑓 − 𝑇𝑎𝑏𝑠) −
𝐾𝑎𝑏𝑠.𝑆𝑎𝑏𝑠 𝑑𝑡 𝑣i−𝑎𝑏𝑠 𝑓−𝑎𝑏𝑠
(𝑇𝑎𝑏𝑠 − 𝑇i𝑠𝑜𝑙)
𝑒𝑎𝑏𝑠 (III. 32)
Au niveau de l’isolant :
(III. 33)
𝑑𝑇i𝑠𝑜𝑙
𝑀 𝐶 =𝜑 −𝜑 +𝜑
i𝑠𝑜𝑙 𝑝i𝑠𝑜𝑙 𝑟𝑠𝑜𝑙−i𝑠𝑜𝑙 𝑐𝑣i𝑠𝑜𝑙−𝑎𝑚
𝑐𝑑𝑎𝑏𝑠 𝑑𝑡
Donc :
𝑑𝑇i𝑠𝑜𝑙
𝑀i𝑠𝑜𝑙𝐶𝑝 = 𝑟 𝑆(𝑇i𝑠𝑜𝑙 − 𝑇𝑠𝑜𝑙)−𝑐𝑣 𝑆(𝑇i𝑠𝑜𝑙 − 𝑇𝑎𝑚) +
i𝑠𝑜𝑙 𝑑𝑡 𝑠𝑜𝑙−i𝑠𝑜𝑙 i𝑠𝑜𝑙−𝑎𝑚 (III. 34)
𝐾𝑎𝑏𝑠𝑆
(𝑇𝑎𝑏𝑠 − 𝑇i𝑠𝑜𝑙)
𝑒𝑎𝑏𝑠
III.1
III.2 Rendement instantané:
Le rendement instantané [17] [14] du capteur que l’on définit par le rapport suivant :
𝑝𝑢i𝑠𝑠𝑎𝑛𝑐𝑒 𝑡𝑒𝑟𝑚i𝑞𝑢𝑒 𝑢𝑡i𝑙𝑒 𝑝𝑎𝑟 𝑚² 𝑑𝑢 𝑐𝑎𝑝𝑡𝑒𝑢𝑟 O (III. 39)
y= = 𝑢
ƒ𝑙𝑢𝑥 𝑠𝑜𝑙𝑎i𝑟𝑒 i𝑛𝑐i𝑑𝑒𝑛𝑡 𝑠𝑢𝑟 𝑙𝑒 𝑝𝑙𝑎𝑛 𝑑𝑒 𝑐𝑎𝑝𝑡𝑒𝑢𝑟 𝐴 𝑐𝐺
Qu:Puissance utile récupérée par le fluide caloporteur (W), elle est donné par :
Qu = 𝑚 𝑐𝑝(𝑇𝑓,𝑠 − 𝑇𝑓,𝑒) (III. 40)
Puissance utile = puissance captée – pertes
Qu = 𝐴𝑐(𝑎𝑟)𝑒𝑓𝑓𝐺 − 𝑈𝐿𝐴𝑐(𝑇𝑎𝑏𝑠 − 𝑇𝑎𝑚) (III. 41)
4
O𝑢 𝐴𝑎𝑏𝑠[(𝑎𝑟)𝑒𝑓𝑓𝐺 − 𝑈𝐿(𝑇𝑎𝑏𝑠 − 𝑇𝑎𝑚)] (III. 42)
y= =
𝐴𝑐𝐺 𝐴𝑐𝐺
𝑎et𝑟: sont respectivement coefficient d'absorption de l'absorbeur et transparence de vitrage
Figure III.6: schéma équivalence électrique-thermique du capteur solaire à air.
conclusion
- En conclusion , le bilan énergétique du capteur solaire met en évidence les performances et
l’efficacité de cette technologie dans la conversion de l’énergie solaire en énergie utilisable.
- L’objectif principal de cette étude était d’évaluer l’efficacité du capteur solaire dans la collecte
et la transformation de l’énergie solaire en chaleur ou en électricité
5
6
CHAPITRE III : MODELISATION ET SIMULATION D’UN CAPTEUR SOLAIRE PLAN Á AIR