Sciences & Technologie A – N°48, Décembre 2018, pp 09-17

CONCEPTION ET RÉALISATION D’UN RÉFRIGÉRATEUR SOLAIRE

À ADSORPTION DESTINÉ À LA CONSERVATION DES VACCINS
DANS LES ZONES SOUDANO SAHÉLIENNES DU CAMEROUN.

ARMAND NOËL NGUECHE CHEDOP 1, NOËL DJONGYANG 1, ZAATRI ABDELOUAHAB 2

1 Department of Renewable Energy, The National Advanced School of Engineering of Maroua, University of Maroua, Cameroon.
2 Department of Mechanical Engineering, Faculty of Engineering Science, Frères Mentouri University Constantine 1, Constantine,
Algeria.

Reçu le 30/04/2018 – Accepté le 12/07/2018

Résumé
Cet article présente la conception et la réalisation d’un réfrigérateur solaire à adsorption utilisant le couple silicagel/eau destiné à la
conservation des vaccins dans les zones soudano sahéliennes du Cameroun. Dans un premier temps, une modélisation complète des
différentes parties de la machine à adsorption solaire a été effectuée, ensuite, un programme de calcul sur Matlab est réalisé en utilisant
les conditions climatiques (rayonnement solaire, ensoleillement et température) de la ville de Maroua. Ceci nous a permis d’effectuer
le dimensionnement des éléments essentiels du prototype à savoir : boitier, absorbeur, bahut, condenseur et évaporateur. Dans notre
programme, en fonction de l’ensoleillement et de la température des mois d’Avril et Août 2016 de la ville de Maroua au Cameroun,
nous avons évalué les températures au niveau du condenseur, de l’évaporateur et de la chambre froide. Ce qui nous a permis d’évaluer
pour ses deux mois l’évolution de l’irradiation solaire et de l’énergie utile de désorption à fin de calculer le coefficient de performance
thermique〖COP〗_th, le coefficient de performance solaire 〖COP〗_Sol et le rendement de captation solaire ɳ de la machine. Les
résultats obtenus montrent un fonctionnement satisfaisant du réfrigérateur solaire dans cette ville, avec des coefficients de performances
thermiques ayant des valeurs comprises entre 0,0327 et 0,1937 pour le mois d’Avril contre 0,0325 et 0,1576 pour le mois d’Août, et
les coefficients de performances solaire compris entre 0,0226 et 0,151 pour le mois d’Avril contre 0,0216 et 0,108 pour le mois d’Août.
Mots clés: Capteur solaire ; Programme ; Prototype ; Réfrigérateur Solaire ; Modélisation ; Conditions climatiques.

Abstract
This article presents the design and implementation of a solar adsorption refrigerator using silica gel / water pair for vaccine storage in
Sudano Sahelian areas of Cameroon. At first, a complete modeling of the different parts of the solar adsorption machine was carried
out, then a calculation program on Matlab is carried out using the climatic conditions (solar radiation, sunshine and temperature) of the
city of Maroua. This allowed us to dimension the essential elements of the prototype namely: case, absorber, sideboard, condenser and
evaporator. In our program, according to the sunshine and the temperature of the months of April and August 2016 of the city of
Maroua in Cameroon, we evaluated the temperatures at the condenser, the evaporator and the cold room. This allowed us to evaluate
for its two months the evolution of solar irradiation and the useful energy of desorption in order to calculate the coefficient of thermal
performance 〖COP〗 _th, the coefficient of solar performance 〖COP〗 _Sol and solar capture efficiency ɳ of the machine. The results
obtained show a satisfactory operation of the solar refrigerator in this city, with thermal performance coefficients having values between
0.0327 and 0.1937 for the month of April against 0.0325 and 0.1576 for the month of April. 'August, and the solar performance
coefficients between 0.0226 and 0.151 for the month of April against 0.0216 and 0.108 for the month of August.
Keywords : Solar captor ; Program ; Prototype; Solar refrigerator; Modelization ; Weather conditions.
‫ﻣﻠﺨﺺ‬
‫ زوج اﻟﻤﺎء ﻟﺘﺨزﯿن اﻟﻠﻘﺎﺤﺎت ﻓﻲ ﻤﻨﺎطق ﺴوداﻨو اﻟﺴﺎﺤل‬/ ‫ﺘﻘدم ﻫذﻩ اﻟﻤﻘﺎﻟﺔ ﺘﺼﻤﯿم وﺘﻨﻔﯿذ ﺜﻼﺠﺔ ﻻﻤﺘﺼﺎص اﻟطﺎﻗﺔ اﻟﺸﻤﺴﯿﺔ ﺒﺎﺴﺘﺨدام ﻫﻼم اﻟﺴﯿﻠﯿﻛﺎ‬
‫ ﺒﺎﺴﺘﺨدام‬Matlab ‫ ﺜم ﯿﺘم ﺘﻨﻔﯿذ ﺒرﻨﺎﻤﺞ ﺤﺴﺎب ﻋﻠﻰ‬، ‫ ﺘم وﻀﻊ ﻨﻤوذج ﻛﺎﻤل ﻟﻸﺠزاء اﻟﻤﺨﺘﻠﻔﺔ ﻤن آﻟﺔ اﻤﺘﺼﺎص اﻟطﺎﻗﺔ اﻟﺸﻤﺴﯿﺔ‬، ‫ ﻓﻲ اﻟﺒداﯿﺔ‬.‫ﻓﻲ اﻟﻛﺎﻤﯿرون‬
، ‫ اﻟﺤﺎﻟﺔ‬:‫ ﻫذا ﺴﻤﺢ ﻟﻨﺎ ﺒﻌدا اﻟﻌﻨﺎﺼر اﻷﺴﺎﺴﯿﺔ ﻟﻠﻨﻤوذج اﻷوﻟﻲ وﻫﻲ‬.‫ أﺸﻌﺔ اﻟﺸﻤس ودرﺠﺔ اﻟﺤ اررة( ﻟﻤدﯿﻨﺔ ﻤﺎروا‬، ‫اﻟظروف اﻟﻤﻨﺎﺨﯿﺔ )اﻹﺸﻌﺎع اﻟﺸﻤﺴﻲ‬
‫ ﻗﻤﻨﺎ ﺒﺘﻘﯿﯿم‬، ‫ ﻟﻤدﯿﻨﺔ ﻤﺎروا ﻓﻲ اﻟﻛﺎﻤﯿرون‬2016 ‫ وﻓﻘًﺎ ﻷﺸﻌﺔ اﻟﺸﻤس ودرﺠﺔ ﺤ اررة ﺸﻬري أﺒرﯿل وأﻏﺴطس‬، ‫ ﻓﻲ ﺒرﻨﺎﻤﺠﻨﺎ‬.‫ اﻟﻤﻛﺜف واﻟﻤﺒﺨر‬، ‫ ﺒوﻓﯿﻪ‬، ‫اﻤﺘﺼﺎص‬
‫ ﺴﻤﺢ ﻟﻨﺎ ذﻟك ﺒﺘﻘﯿﯿم ﺘطور اﻹﺸﻌﺎع اﻟﺸﻤﺴﻲ واﻟطﺎﻗﺔ اﻟﻤﻤﺘﺎزة ﻟﻼﻤﺘﺼﺎص ﻟﻤدة ﺸﻬرﯿن ﻤن أﺠل ﺤﺴﺎب‬.‫درﺠﺎت اﻟﺤ اررة ﻓﻲ اﻟﻤﻛﺜف واﻟﻤﺒﺨر واﻟﻐرﻓﺔ اﻟﺒﺎردة‬
‫ أظﻬرت اﻟﻨﺘﺎﺌﺞ اﻟﺘﻲ ﺘم اﻟﺤﺼول‬.‫ ﻤن اﻟﺠﻬﺎز‬ɳ ‫ وﻛﻔﺎءة اﻟﺘﻘﺎط اﻟطﺎﻗﺔ اﻟﺸﻤﺴﯿﺔ‬COP〗 _Sol〗 ‫ ﻤﻌﺎﻤل اﻷداء اﻟﺸﻤﺴﻲ‬، COP 〗_th〗 ‫ﻤﻌﺎﻤل اﻷداء اﻟﺤراري‬
0.0325 ‫ ﻟﺸﻬر أﺒرﯿل ﻤﻘﺎﺒل‬0.1937 ‫ و‬0.0327 ‫ ﻤﻊ وﺠود ﻤﻌﺎﻤﻼت أداء ﺤ اررﯿﺔ ﺘﺘراوح ﻗﯿﻤﺘﻬﺎ ﺒﯿن‬، ‫ﻋﻠﯿﻬﺎ ﻋﻤﻠﯿﺔ ﻤرﻀﯿﺔ ﻟﻠﺜﻼﺠﺔ اﻟﺸﻤﺴﯿﺔ ﻓﻲ ﻫذﻩ اﻟﻤدﯿﻨﺔ‬
.‫ ﻟﺸﻬر أﻏﺴطس‬0.108 ‫ و‬0.0216 ‫ ﻟﺸﻬر أﺒرﯿل ﻤﻘﺎﺒل‬0.151 ‫ و‬0.0226 ‫ أﻏﺴطس وﻤﻌﺎﻤﻼت اﻷداء اﻟﺸﻤﺴﻲ ﺒﯿن‬.‫ ﻟﺸﻬر أﺒرﯿل‬0.1576 ‫و‬

.‫ اﻟﺜﻼﺟﺔ اﻟﺸﻤﺴﯿﺔ ؛ ﺗﺼﻤﯿﻢ ؛ اﻟﺸﺮوط اﻟﻤﻨﺎﺧﯿﺔ‬.‫اﻟﻜﻠﻤﺎت اﻟﻤﻔﺘﺎﺣﯿﺔ ﺟﺎﻣﻊ اﻟﻄﺎﻗﺔ اﻟﺸﻤﺴﯿﺔ اﻟﺒﺮﻧﺎﻣﺞ؛ اﻟﻨﻤﻮذج‬

© Université Frères Mentouri Constantine 1, Algérie, 2018

 ARMAND NOËL NGUECHE CHEDOP, NOËL DJONGYANG et ZAATRI ABDELOUAHAB

NOMENCLATURE INTRODUCTION

Symboles : L’utilisation de l’énergie solaire dans les pays

Mg: masse de glace [Kg] ensoleillés est un moyen efficace pour pallier au manque
d’énergie surtout dans les zones rurales où il est parfois
ΔT: différence de température [K]
difficile et coûteux de les alimenter avec le réseau électrique
q: la fraction adsorbée conventionnel. Parmi les procédés de transformations
P: la pression d'équilibre d'adsorption, thermiques de l’énergie solaire, la réfrigération solaire est
M l M
PS: la pression d' vapeur
de i saturante de l’application la plus indiquée pour le stockage des denrées
l'E:dénergie
b àsolaire
l alimentaires, des produits pharmaceutiques et le confort
utile ou entréeidé
é é
solaire
thermique dans l’habitat puis qu’elle contribue efficacement
(kWh
H / J)
J: rayonnement solaire global reçu dans à la réduction du fort pic de la consommation d’énergie
lA: surface du collecteur
l i (kWh / 2 (m²).
solaire J) électrique en période de canicule responsable de la plupart
T = température dans le capteur solaire des problèmes d’approvisionnement. Associées aux
[K]
I (t): est l'éclairement reçu par la face éventuelles fuites de fluides frigorigènes, ces pics de
d grecques
ll production électrique induisent une augmentation des
Lettres : l i émissions de gaz à effet de serre, accentuant le cercle vicieux
ρ masse volumique, kgm-3 du changement climatique. Par ailleurs, les régions
λ conductivité thermique, W·m-1·K-1 tropicales sèches sont des zones où le potentiel solaire est
très important, l’ensoleillement annuel est toujours supérieur
η : rendement solaire (%) à 20 000 kJ par m² de surface de captation [1-5]. Il est donc
σ = 5.67 x 10-8Stephan-Boltz important d’exploiter cette ressource naturelle notamment
η : Solar Yield (%) dans le domaine de la production de froid. Si les méthodes
σ = 5.67 x 10-8 Stephan-Boltz de production de froid et de climatisation classiques sont un
φ fraction massique frein pour l’efficacité énergétique en particulier dans ses
régions tropicales sèches, pourquoi ne donc pas utiliser
κ constante de Von Karman
comme énergie de substitution l’énergie solaire ?
L’exploitation de l’énergie solaire n’est-elle pas une solution
Indices / Exposants : aux divers problèmes énergétique et climatique que connait
Lfu: chaleur latente [J / Kg] l’Afrique ? Vue le faible produit international brute (PIB) de
Cp: la chaleur spécifique de l'eau la plupart des pays Africains, est-il véritablement possible
d’envisager un développement à grand échelle de
Ta: est la température ambiante [K] l’exploitation de l’énergie solaire à des fins socio-
Levap: la chaleur latente de l'évaporation de environnementaux telles que la climatisation des locaux et
l'Lfu: chaleur
[J / K ]latente [J / Kg] la conservation des denrées alimentaires ? Toutes ces
Bener: les besoins énergétiques (kWh / J) questions imposent donc une étude première de faisabilité.
Plusieurs travaux dans le domaine de la réfrigération par
Ad = gain solaire diffus [W] sorption (absorption et adsorption) ont été couronnés de
Rayon = pertes de rayonnement [W] succès et des réfrigérateurs ont été réalisés. Cependant, ils
Tair = Température ambiante [K] sont parfois peu adaptés puisque rares sont ceux qui ont été
Bener: les besoins énergétiques (kWh / J) testés sous des conditions météorologiques réelles.
Ad = gain solaire diffus [W] Les régions sahéliennes en générales et celles du
COPth: Coefficient de performance Cameroun en particulier sont des zones les plus souvent
h Soli : coefficient de performance solaire dépourvues d’énergies électriques et en période très chaude
COP
(canicule), on observe une augmentation de taux de
Md: Masse de réfrigérant désorbée [Kg] mortalité dans ses zones à cause des maladies telles que la
CP: Chaleur spécifique du réfrigérant [J méningite, le choléra et la poliomyélite. Fort de ce constat il
K
Qe est1 Kle1]
froid produit à l'évaporateur est alors important de se demander comment réduire ses taux
QC: Quantité de chaleur extraite au de mortalité dans ses zones en proposant un réfrigérateur
solaire pour la conservation des vaccins.
UE: dIrradiation nécessaire pour effectuer la
dé Esti la surface avant du capteur
Sav: L’objectif final de ce travail est de concevoir et de
réaliser au laboratoire des Energies Renouvelables de
l’Université de Maroua un réfrigérateur fonctionnant à
l'énergie solaire et permettant de refroidir 5 kg d'eau de 35°C
afin de transporter des vaccins dans des régions reculées. La
glace créée devra se conserver au moins 3 jours en absence
de soleil. La première partie présente la modélisation
complète du réfrigérateur solaire à adsorption. Dans la
deuxième partie, nous présentons le prototype réalisé ainsi

10
 CONCEPTION ET RÉALISATION D’UN RÉFRIGÉRATEUR SOLAIRE À ADSORPTION DESTINÉ À LA CONSERVATION DES VACCINS DANS LES
ZONES SOUDANO SAHÉLIENNES DU CAMEROUN.

que le dimensionnement complet de notre réfrigérateur. • La chauffe portant la température de Tinitiale à

Ensuite une description du programme de calcul développé
Tdésorption [7]
sur le logiciel Matlab est présentée. Enfin la dernière partie
présente les résultats et discussions. Mtot,CPglobal
𝒅𝒅𝒅𝒅
= AD + Ad – Ray (9)
𝒅𝒅𝒅𝒅
2. MODELE MATHEMATIQUE/METHODE
EXPERIMENTALE : Avec [8-9]:

La formule ci-contre nous permet d’obtenir l’énergie Mtot = masse totale [Kg] + massesilicagel + masseacier
nécessaire pour une masse de la glace partir d’une donnée 𝑻𝑻𝟐𝟐
[7-12]: T= 𝒂𝒂𝒂𝒂𝒂𝒂
(10)
𝑻𝑻𝒆𝒆𝒆𝒆

E = 𝑴𝑴𝒈𝒈𝒈𝒈𝒈𝒈𝒈𝒈𝒈𝒈 (Cp*∆T + Lfus) (1) • La masse de l’acier

∆T = (𝑻𝑻𝒂𝒂 - 𝑻𝑻𝒈𝒈 ) (2) Macier = 𝝆𝝆acierVacier= 𝝆𝝆acier (Vtoles + Vgrillage) (11)
2.1. La masse d’eau évaporée • La masse totale
Cette énergie conduit à l’évaporation d’une masse Mtot= Meau + Macier + MSiO2 (12)
d’eau de [7-12]:

Meau =
𝑬𝑬
(3) • La chaleur spécifique globale
𝑳𝑳𝒆𝒆𝒆𝒆𝒆𝒆𝒆𝒆
𝑀𝑀𝑒𝑒𝑒𝑒𝑒𝑒 𝐶𝐶𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃 +𝑀𝑀𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎 𝐶𝐶 𝑀𝑀𝑆𝑆𝑆𝑆𝑆𝑆2 𝐶𝐶
𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃 +
2.2. La masse de l’adsorbant Cpglobale = 𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃2
(13)
𝑀𝑀𝑡𝑡𝑡𝑡𝑡𝑡

La masse de l’adsorbant (silicagel) nécessaire pour une • Le temps de chauffe

masse d’eau donnée correspond à environ 35 à 40 % et est
donnée par [13]: En introduisant ces valeurs dans l’équation de bilan
𝑴𝑴𝒆𝒆𝒆𝒆𝒆𝒆 thermique, et en intégrant de Tinitiale à Tdésorption pour une
𝑴𝑴𝑺𝑺𝑺𝑺𝑺𝑺𝑺𝑺𝑺𝑺𝑺𝑺𝑺𝑺𝑺𝑺𝑺𝑺 = (4) durée allant de 0 à tchauffe, on peut obtenir le temps de chauffe
𝜟𝜟𝜟𝜟
:
La masse de silicagel à introduire est multiplié par un
coefficient de sécurité de 2,5. La quantité ΔM est obtenue à 𝑻𝑻 𝒅𝒅𝒅𝒅 𝟏𝟏 𝒕𝒕 (𝑪𝑪𝑪𝑪𝑪𝑪𝑪𝑪𝑪𝑪𝑪𝑪𝑪𝑪)
∫𝑻𝑻 𝒅𝒅𝒅𝒅𝒅𝒅𝒅𝒅𝒅𝒅𝒅𝒅𝒅𝒅𝒅𝒅𝒅𝒅𝒅𝒅 𝑨𝑨 = = ∫𝟎𝟎 𝒅𝒅𝒅𝒅
partir du diagramme isostérique sur lequel on trace un cycle 𝑰𝑰𝑰𝑰𝑰𝑰𝑰𝑰𝑰𝑰𝑰𝑰𝑰𝑰𝑰𝑰 𝑫𝑫 + 𝑨𝑨𝒅𝒅 −𝑹𝑹𝑹𝑹𝑹𝑹 𝑴𝑴𝒕𝒕𝒕𝒕𝒕𝒕 𝑪𝑪𝑷𝑷𝑷𝑷𝑷𝑷𝑷𝑷𝑷𝑷𝑷𝑷𝑷𝑷
caractérisé par quatre températures (tads, treg, tevap et tcond) ou
calculé à partir de l'équation de Dubinin-astaknov, qui s'écrit (14)
de la façon suivante [14-16]. De cette équation, on obtient le temps de chauffe.
𝒒𝒒 𝑷𝑷𝑷𝑷
ln( ) = -D[𝑻𝑻. 𝐥𝐥𝐥𝐥( )]𝒏𝒏 (5) 2.3.2. La désorption
𝒘𝒘𝟎𝟎 𝝆𝝆 𝑷𝑷
Cette phase nous permet de calculer le temps
2.3. Choix des éléments essentiels de la machine :
nécessaire à désorber toute l’eau du silicagel (tdésorption) :
Les éléments essentiels de la machine sont le capteur,
le condenseur, l’évaporateur, l’enceinte frigorifique et le 𝑨𝑨𝑫𝑫 + 𝑨𝑨𝒅𝒅 − 𝑹𝑹𝑹𝑹𝑹𝑹 = Ldésorption*J (15)
bahut. 𝑴𝑴𝒆𝒆𝒆𝒆𝒆𝒆
Le temps de désorption est alors : Tdésorption =
𝒋𝒋
2.3.1. Capteur :
(16)
• Détermination de la surface du capteur (première
2.3.3. Evaluation des performances de la machine
étape) :
L’évaluation des performances de la machine est
Le capteur, élément essentiel de la machine, est un déterminée à partir des quantités de chaleur Qe , Eu et E :
parallélépipède en aluminium, contenant l’adsorbant, sa
surface extérieure est couverte d’un vitrage. Après avoir • Quantité de chaleur extraite à l’évaporateur [20-22].
déterminé sa surface, nous procéderons à sa fabrication car
il est simple à réaliser. L’énergie produite par les capteurs Qe = 𝚫𝚫𝚫𝚫. MaL – MdCp (TC - Te). 𝚫𝚫𝚫𝚫 (17)
solaires s’exprime selon la formule [17] :
• Quantité de chaleur extraite au niveau du condenseur
E = ɳ*A*Hj (6) [20-22].
𝑩𝑩𝑩𝑩𝑩𝑩𝑩𝑩𝑩𝑩
A= (7)
ɳ∗𝑯𝑯𝒋𝒋 QC = 𝚫𝚫𝚫𝚫. MaL – MdCp (Tg – TC). 𝚫𝚫𝚫𝚫 (18)
• Les besoins énergétiques sont déterminés par • Quantité de chaleur produite par l’adsorbant pendant
l’équation : l’adsorption [23].
Bener = 𝑴𝑴𝒈𝒈𝒈𝒈𝒈𝒈𝒈𝒈𝒈𝒈 (Cp*∆T + Lfus) (8)

11
 ARMAND NOËL NGUECHE CHEDOP, NOËL DJONGYANG et ZAATRI ABDELOUAHAB

𝐓𝐓 𝐓𝐓
QA = ∫𝐓𝐓 𝐚𝐚𝐚𝐚𝐚𝐚 (𝐢𝐢𝐢𝐢𝐢𝐢𝐢𝐢𝐢𝐢𝐢𝐢𝐢𝐢) ∑ 𝐌𝐌𝐢𝐢 𝐂𝐂𝐏𝐏𝐢𝐢 . 𝐝𝐝𝐝𝐝 + ∫𝐓𝐓 𝐚𝐚𝐚𝐚𝐚𝐚 (𝐟𝐟𝐟𝐟𝐟𝐟) ∑ 𝐌𝐌𝐢𝐢 𝐂𝐂𝐏𝐏𝐢𝐢 . 𝐝𝐝𝐝𝐝 + �(AB) de chauffage ;
𝐓𝐓
𝐠𝐠 𝐚𝐚𝐚𝐚𝐚𝐚 (𝐢𝐢𝐢𝐢𝐢𝐢𝐢𝐢𝐢𝐢𝐢𝐢𝐢𝐢) �(BC) de désorption à pression presque constante ;
∫𝐓𝐓 𝐚𝐚𝐚𝐚𝐚𝐚 (𝐟𝐟𝐟𝐟𝐟𝐟) 𝚫𝚫𝐇𝐇. 𝚫𝚫𝚫𝚫. 𝐌𝐌𝐚𝐚 (19) �(CD) du refroidissement ;
𝐚𝐚𝐚𝐚𝐚𝐚 (𝐢𝐢𝐢𝐢𝐢𝐢𝐢𝐢𝐢𝐢𝐢𝐢𝐢𝐢)
�(DA) d’adsorption (production de froid) à pression
• Irradiation nécessaire pour effectuer la désorption : constante.

𝒇𝒇𝒇𝒇𝒇𝒇 𝒅𝒅𝒅𝒅 𝒍𝒍𝒍𝒍 𝒅𝒅é𝒔𝒔𝒔𝒔𝒔𝒔𝒔𝒔𝒔𝒔𝒔𝒔𝒔𝒔𝒔𝒔 On remarque, cependant, au début de la période (DA) une
Eu = Sav ∫𝒍𝒍𝒍𝒍𝒍𝒍é 𝒅𝒅𝒅𝒅 𝒋𝒋𝒋𝒋𝒋𝒋𝒋𝒋 𝑰𝑰(𝒕𝒕)𝒅𝒅𝒅𝒅 (20) élévation de pression à l’intérieur de l’installation
correspondant au début de la production du froid. Cette
• Irradiation globale reçue par le capteur [20-22]. augmentation est due à la chaleur dégagée par le silicagel .
L’effet résultant est un changement de phase : de la phase
𝐓𝐓 𝐓𝐓
E = ∫𝐓𝐓 𝐝𝐝𝐝𝐝𝐝𝐝 ∑ 𝐌𝐌𝐢𝐢 𝐂𝐂𝐏𝐏𝐢𝐢 . 𝐝𝐝𝐝𝐝 + ∫𝐓𝐓 𝐠𝐠 ∑ 𝐌𝐌𝐢𝐢 𝐂𝐂𝐏𝐏𝐢𝐢 . 𝐝𝐝𝐝𝐝 + ΔH. 𝚫𝚫𝚫𝚫. Ma vapeur à la phase absorbée.
𝐚𝐚𝐚𝐚𝐚𝐚 𝐝𝐝𝐝𝐝𝐝𝐝

(21) 50
Ta
T(°C)
∑ 𝐌𝐌𝐢𝐢 𝐂𝐂𝐏𝐏𝐢𝐢 . 𝐝𝐝𝐝𝐝 = ∑ 𝐐𝐐𝐢𝐢 = δ𝑸𝑸𝟏𝟏 + δ𝑸𝑸𝟐𝟐 + δ𝑸𝑸𝟑𝟑 (22)
45
TC
40

35
• Coefficient de performance thermique [8]. 30

25
𝑸𝑸𝒆𝒆
COPth = (23) 20
𝑬𝑬 15

• Coefficient de performance solaire [8]. 10

5
𝑸𝑸𝒆𝒆
COPsol = (24) 0
𝑬𝑬𝒖𝒖 2 4 6 8 10 12 14 16 18 20 22 24
Time(Hours)
• Rendement de captation solaire [8].
𝑬𝑬𝒖𝒖 𝐂𝐂𝐂𝐂𝐂𝐂𝐬𝐬𝐬𝐬𝐬𝐬
ɳ= = (25)
𝑬𝑬 𝐂𝐂𝐂𝐂𝐂𝐂𝐭𝐭𝐭𝐭 140 P(mbar) C
B
120
3. RESULTATS
100

Parmi les résultats obtenus nous allons présenter deux 80

séquences de fonctionnement. La première s’étend du 1er au 60
14 Avril 2016, la seconde est obtenue pour la période du 7
au 19 Août 2016. 40

20
Pour illustrer le fonctionnement du réfrigérateur solaire dans D A D
le climat de Maroua, nous avons choisi quatre jours de la 0
2 4 6 8 10 12 14 16 18 20 22 24
séquence d’Avril de conditions climatiques différentes : Time(Hours)

� un jour clair et bien ensoleillé (11 Avril 2016) ;

� un jour clair le matin avec quelques passages nuageux
l’après-midi (4 Avril 2016) ;
25
T(°C) Te
20 TCF
� un jour d’ensoleillement moyen (5 Avril 2016) ;
� un jour de mauvaises conditions climatiques avec des 15

éclaircies (8 Avril 2016). 10

Les figures 5, 6, 7 et 8, donnent, pour ces jours, l’évolution 5

: 0

-(a) des températures du condenseur et de l’ambiance ; -5

-(b) de la pression dans l’installation ; 2 4 6 8 10 12 14 16 18 20 22 24

Time(Hours)
-(c) des températures de la chambre froide et de
l’évaporateur.
On remarque que l’unité produit toujours du froid Figure 5 : Evolution le 4 Avril de : (a) la température du
même pour la journée de mauvais ensoleillement. condenseur et de l’ambiance, (b) la pression dans
l’installation, (c) la température de l’évaporateur et de la
L’enchaînement des quatre étapes du cycle de chambre froide.
fonctionnement est représenté sur la courbe de pression pour
la journée du 4 Avril (figure 5 (b)). En effet, on retrouve les Les courbes donnant l’évolution de la température du
quatre périodes : condenseur, montrent que ce dernier subit un échauffement
pendant la période correspondant à la phase de désorption

12
 CONCEPTION ET RÉALISATION D’UN RÉFRIGÉRATEUR SOLAIRE À ADSORPTION DESTINÉ À LA CONSERVATION DES VACCINS DANS LES
ZONES SOUDANO SAHÉLIENNES DU CAMEROUN.

(figures 5 (a), 6 (a) et 7 (a)) dû à la chaleur latente de

condensation de l’eau. Plus cette température est grande, 50
plus la quantité de vapeur d’eau désorbée est importante. 45
T(°C) Ta
40
TC
35
50
T(°C) Ta 30
45 TC 25
40
20
35
15
30
10
25
5
20
0
15
2 4 6 8 10 12 14 16 18 20 22 24
10 Time(Hours)
5

0
2 4 6 8 10 12 14 16 18 20 22 24
Time(Hours) 240
P(mbar)
220
200

160 180
P(mbar) 160
140
140
120 120

100 100
80
80
60
60 40
20
40
0
20 2 4 6 8 10 12 14 16 18 20 22 24
Time(Hours)
0
2 4 6 8 10 12 14 16 18 20 22 24
Time(Hours)
25
T(°C) Te
25
T(°C)
Te 20 TCF
20 TCF 15

15
10
10
5
5
0
0

-5
-5

-10 -10
2 4 6 8 10 12 14 16 18 20 22 24 2 4 6 8 10 12 14 16 18 20 22 24
Time(Hours) Time(Hours)

Figure 6 : Evolution le 11 Avril de : (a) la température du Figure 7: Evolution le 5 Avril de : (a) la température du
condenseur et de l’ambiance, (b) la pression dans condenseur et de l’ambiance, (b) la pression dans
l’installation, (c) la température de l’évaporateur et de la l’installation, (c) la température de l’évaporateur et de la
chambre froide. chambre froide.

13
 ARMAND NOËL NGUECHE CHEDOP, NOËL DJONGYANG et ZAATRI ABDELOUAHAB

Il est intéressant de noter que la température ambiante

Pour la journée du 8 Avril caractérisée par un ensoleillement du jour a une influence sur la période de désorption. En effet,
faible, l’évolution de la température du condenseur suit à peu les courbes donnant l’évolution de la pression montrent que
près celle de la température ambiante et donc ne subit aucune pour deux journées d’ensoleillements voisins (4 et 5 Avril),
augmentation (figure 8 (a)). Ce qui indique que la quantité la journée du 5 Avril, caractérisée par une température
désorbée est nulle. La courbe de la figure 8 (c) montre, par ambiante (36.7 °C) plus élevée que celle de la journée du 4
contre, une production de froid à l’évaporateur, ce qui Avril (28°C) connaît un retard à la désorption qui ne
indique que l’unité a adsorbé une quantité d’eau non commence que lorsque la pression est supérieure à 180 mbar
négligeable. (figure 5 (b), figure 7 (b)). Cette remarque rejoint le fait que
la désorption ne commence que lorsque la pression du
capteur atteint celle du condenseur qui dépend de la
50 Ta température ambiante.
T(°C)
45 TC
40 Lorsque les conditions climatiques sont bonnes (journée
35
ensoleillée et claire, température ambiante faible) comme
celle de la journée du 11 Avril, les performances sont
30
optimales.
25

20

15 28000 2
26000
E(Kj/m )
10

5
24000 E
0
22000 EU
2 4 6 8 10 12 14 16 18 20 22 24 20000
Time(Hours) 18000
16000
14000
180
P(mbar) 12000
160 10000
140 8000

120 01/04 02/04 03/04 04/04 05/04 06/04 07/04 08/04 09/04 10/04 11/04 12/04 13/04 14/04
April 2016
100

80
Figure 9: Evolution de l’irradiation solaire (E) et de
60 l’énergie utile de désorption (EU) dans la ville de Maroua
40 au mois d’Avril 2016.
20

0
2 4 6 8 10 12 14 16 18 20 22 24
Time(Hours)

30
T(°C) Te
25 TCF
20

15

10

-5
Figure 10: Evolution de l’irradiation solaire (E) et de
l’énergie utile de désorption (EU) dans la ville de Maroua
-10
2 4 6 8 10 12 14 16 18 20 22 24
au mois d’Août 2016.
Time(Hours)
Les figures 9 et 10 présentent respectivement l’irradiation
solaire (E) et l’énergie utile de désorption (Eu) des séquences
Figure 8: Evolution le 8 Avril de : (a) la température du des mois d’Avril et Août 2016. Nous remarquons que
condenseur et de l’ambiance, (b) la pression dans l’irradiation maximale était obtenue le 11/04 avec une
l’installation, (c) la température de l’évaporateur et de la valeur de 27308 Kj/m² pour la séquence d’Avril contre une
chambre froide. valeur maximale de 25507 Kj/m² obtenue le 09/08 pour la

14
 CONCEPTION ET RÉALISATION D’UN RÉFRIGÉRATEUR SOLAIRE À ADSORPTION DESTINÉ À LA CONSERVATION DES VACCINS DANS LES
ZONES SOUDANO SAHÉLIENNES DU CAMEROUN.

séquence d’Août. Nous remarquons également que deux

journées qui se suivent peuvent avoir des valeurs 0,85

d’irradiation complètement différentes 0,80

Solar Yield
0,75
0,20 COP
0,18 COPth 0,70

0,16 COPSol 0,65

0,14
0,60
0,12

0,10 0,55
0,08
0,50
0,06
0,45
0,04
01/04 02/04 03/04 04/04 05/04 06/04 07/04 08/04 09/04 10/04 11/04 12/04 13/04 14/04
0,02 April 2016
01/04 02/04 03/04 04/04 05/04 06/04 07/04 08/04 09/04 10/04 11/04 12/04 13/04 14/04
April 2016
Figure 13: Evolution du rendement de captation solaire
Figure 11: Evolution des coefficients de performance pour le mois d’Avril 2016
thermique (COPth) et solaire COPsol) au mois d’Avril 2016 La figure 13 montre l’évolution du rendement de captation
La figure 11 présente l’évolution des coefficients de solaire pour les jours de la séquence d’Avril. Nous
performance thermique (COPth) et solaire (COPsol) pour les remarquons qu’il est compris entre 0,49 et 0,81 avec une
jours de la séquence du mois d’Avril 2016. Nous observons valeur maximale obtenue le 14/04. Nous pouvons conclure
que le COPth est entre 0,03 et 0,19 avec un pic obtenu le qu’une grande partie de l’énergie solaire est transformée en
09/04. Le COPsol est compris entre 0,02 et 0,15 ; il est chaleur dans le capteur et que le rendement de captation
supérieur à 0,O9 quand l’irradiation solaire est supérieur à solaire est supérieur à 0,7 quand la journée est claire le
20000 KJ/m². Nous pouvons donc conclure que le mois matin.
d’Avril à des performances satisfaisantes pour le
fonctionnement de notre réfrigérateur solaire. 0,85
Solar Yield
0,80
0,16 COP
COPth 0,75
0,14
COPSol
0,12 0,70

0,10 0,65

0,08 0,60

0,06 0,55
0,04
0,50
0,02
0,45
07/08 08/08 09/08 10/08 11/08 12/08 13/08 14/08 15/08 16/08 17/08 18/08 19/08 07/08 08/08 09/08 10/08 11/08 12/08 13/08 14/08 15/08 16/08 17/08 18/08 19/08
August 2016 August 2016

Figure 12: Evolution des coefficients de performance

thermique (COPth) et solaire COPsol) au mois d’Août 2016. Figure 14: Evolution du rendement de captation solaire
pour le mois d’Août 2016.
La figure 12 présente l’évolution des coefficients de
performance thermique (COPth) et solaire (COPsol) pour les La figure 14 présente l’évolution du rendement de
jours de la séquence du mois d’Août 2016. Nous observons captation solaire pour les jours de la séquence d’Août. Nous
que le COPth est entre 0,03 et 0,15 avec un pic obtenu le remarquons qu’il est compris entre 0,48 et 0,80 avec une
15/08. Le COPsol est compris entre 0,02 et 0,10. Nous valeur maximale obtenue le 11/08. Nous pouvons conclure
remarquons également que pour des journées que le rendement de captation solaire n’est pas tellement
d’ensoleillement faible le système continue à fonctionner. différent pour les deux séquences des mois d’Avril et Août.
Nous pouvons donc conclure que le mois d’Août malgré 4 . CONCLUSIONS
qu’étant le mois le plus défavorable de l’année 2016, à des
performances permettant le fonctionnement de notre Dans ce travail, nous avons effectué dans un premier temps
réfrigérateur solaire. une modélisation complète du réfrigérateur solaire à
adsorption. Ensuite nous avons présenté le prototype réalisé
ainsi que le dimensionnement complet des éléments
essentiels de notre réfrigérateur à savoir : Capteur solaire,
bahuts et portes, évaporateur et bacs, condenseur et

15
 ARMAND NOËL NGUECHE CHEDOP, NOËL DJONGYANG et ZAATRI ABDELOUAHAB

tuyauterie. Enfin le développement d’un programme de d’ingénieurs du Canton de Vaud / HES-SO / LESBAT
calcul sur le logiciel matlab en utilisant les conditions – EIVD.
climatiques (rayonnement solaire et température) des mois [9] Medini N, Marmottant B, El Golli S, Grenier P. Etude
d’Avril et Août 2016 de la ville de Maroua, nous a permis d'une machine solaire autonome à fabriquer de la glace.
dans un premier temps d’évaluer les températures au niveau Int. J. Refrig. Vo1 14 (1991) 363-367.
du condenseur, de l’évaporateur et de la chambre froide,
ensuite d’évaluer pour ses deux mois l’évolution de [8] Pralon A, Leite F, Daguenet M. Performance of a new
l’irradiation solaire et de l’énergie utile de désorption à fin solid adsorption ice maker with solar energy
de calculer le coefficient de performance thermique 𝐶𝐶𝐶𝐶𝐶𝐶𝑡𝑡ℎ , regeneration. Energy Conversion & Management 41
le coefficient de performance solaire 𝐶𝐶𝐶𝐶𝐶𝐶𝑆𝑆𝑆𝑆𝑆𝑆 et le rendement (2000) 1625-1647.
de captation solaire ɳ de la machine. Les résultats obtenus [9] Hildbrand C, Dind P, Pons M, Buchter F. A new solar
montrent un fonctionnement satisfaisant du réfrigérateur powered adsorption refrigerator with high
solaire dans cette ville, nous avons pu atteindre des Performance. HES-SO - Ecole D'Ingénieurs du Canton
températures basses à l’évaporateur allant jusqu’à -8,1°C De Vaud - LESBAT (2001).
pour le mois d’Avril avec des conditions climatiques [10] Lemmini F, Errougani A, Bentayeb F. Réalisation
favorables ainsi que des coefficients de performances expérimentale d’un réfrigérateur solaire a adsorption à
thermiques ayant des valeurs comprises entre 0,0327 et rabat. Laboratoire d’Energie Solaire, Faculté des
0,1937 pour le mois d’Avril contre 0,0325 et 0,1576 pour le Sciences
mois d’Août, et les coefficients de performances solaire
compris entre 0,0226 et 0,151 pour le mois d’Avril contre [11] Lemmini F, Errougani A. Building and experimentation
0,0216 et 0,108 pour le mois d’Août. Le rendement de of a solar powered adsorption refrigerator. Renewable
captation solaire est compris entre 0,49 et 0,81 pour le mois Energy 30 (2005) 1989– 2003.
d’Avril et 0,47 à 0,78 pour le mois d’Août. Dans la suite de [12] Mayor J, Dind P. Construction et test d’un réfrigérateur
nos travaux nous allons finaliser le prototype et effectuer des solaire à adsorption transportable. Laboratoire
tests afin de faire une comparaison des résultats obtenus d'Energétique Solaire et de Physique du bâtiment,
avec les résultats expérimentaux. HES-SO/EIVD (2003). Yverdon-les-Bains. Suisse.
REFERENCES [13] Abu-Hamdeh NH, Alnefaie KA, Almitani KH. Design
and performance characteristics of solar adsorption
[1] [1] R.L. Webb, E.R.G. Eckert and R.J. Goldstein, refrigeration system using parabolic trough collector:
‘Heat Transfer and Friction in Tubes with Repeated Experimental and statistical optimization technique.
Rip Roughness’, International Journal of Heat and Energy Conversion and Management 74 (2013) 162–
Mass Transfer, Vol. 14, N°4, pp.601 – 617, 1971. 170.
[2] K.A. Joudi and A.I. Mohamed, ‘Expérimental [14] Wang J, Wanga LW, Luo WL, Wanga RZ,
Performance of a Solar Air Heater with a ‘V’ Experimental study of a two-stage adsorption freezing
Corrugated Absorber’, Energy Conversion and machine driven by low temperature heat source.
Management, Vol. 26, N°2, pp. 193 – 200, 1986. International Journal of Refrigeration, 36 (2013 )1029-
1036.
[3] M.R. Zugari and J.J. Vullierne, ‘Amélioration des
Performances Thermiques d’un Capteur [15] Aklouche S, Contribution à la conception d’une
chambre froide utilisant le couple zéolithe-eau et
Solaire par l’Utilisation d’une Structure à Lamelles’,
fonctionnant à l’énergie solaire, thèse de magister
Entropie, Vol. 29, N°176, pp. 25 – 30, 1993.
USTHB, spécialité énergie solaire, option conversion
[4] A. N. Ngueche Chedop, N. Djongyang and A. Zaatri. thermodynamique.
‘’Modelling of the performance of a solar electric-
[16] Yahi F, Abbad B, Berdja M, Ouali M, Bouzafour F.
vapor compression refrigeration system in dry
experimental evaluation of a semi pilote cop of solar
tropical regions’’. International Journal of Science
adsorption machine with intermittent cycle,
and Research, Vol. 3, N°11, pp. 1066 – 1076, 2014.
international workshop on new working fluids for
[5] Noël Djongyang, René Tchinda, Donatien Njomo. absorption heat pumps and refrigeration systems, 22-
‘’Estimation of some comfort parameters for 23 july, 2013,at Taragona, Spain.
sleeping environments in dry-tropical sub-Saharan
[17] M. Chahid, contribution à l’étude des machines
Africa region’’. Energy Conversion and
thermiques à adsorption, Thèse de troisième cycle,
Management, 58 (2012), P.110-119.
1985, I. N. P. L, Nancy France.
[6] Chedop ANN, Djongyang N, Tchinda R, Zaatri A, Kana
[18] Tchernev. Closed cycle Zeolite regenerative heat pump.
JD. A field study of indicators of the performance of
Heat transfer enchancement and energy conservation ,
four absorption chillers in the sudano-sahelian region
Hemisphere Publ. Co,1988, pp. 747–55.
of Cameroon. J.
[19] M. Pons and J.J. Guilleminot, Design of an
Fundam. Appl. Sci., 2016, 8(1), 20-40.
experimental solar powered solid- adsorption ice
[7] Mayor J., Dind Ph. (2002). Construction et test d’un maker. ASME J Solar Energy Engng 108, 1986, pp.
réfrigérateur solaire à adsorption transportable. Ecole 332–337.

16
 CONCEPTION ET RÉALISATION D’UN RÉFRIGÉRATEUR SOLAIRE À ADSORPTION DESTINÉ À LA CONSERVATION DES VACCINS DANS LES
ZONES SOUDANO SAHÉLIENNES DU CAMEROUN.

[20] M. Pons and P.H. Grenier, A phenomenological adsorption of ammonia in a cylindrical adsorbent bed:
adsorption equilibrium law extracted from thermal performance study of a combined parabolic
experimental and theoretical considerations applied to solar collector, water heat pipe and adsorber generator
the activated carbon+methanol pair, Carbon, volume assembly, Applied Thermal Engineering, volume 24,
24, 1987, pp. 615–625. 2004, Pages 2537-2555.
[21] Nolwenn Le Pierrés (2005). Procédé solaire de [24] Vichan Tangkengsirisin, Atsushi Kanzawa, Takayuki
production de froid basse température (-28C) par Watanabe, A solar powered adsorption cooling system
sorption solide-gaz, Thèse présentée en vue de using a silica gel-water mixture, Energy, volume 23,
l’obtention de grade de Docteur en Sciences de 1998, pp. 347-353. 123.
l’Ingénieur de l’Université de PERPIGNAN.
[22] R.E. Critoph, Limitations of adsorption cycles for solar
cooling. Solar Energy 41, 1988, pp. 21–31.
[23] F. Aghbalou, A. Mimet, F. Badia, J. Illa, A. El Bouardi
and J. Bougard, Heat and mass transfer during

17