Simulation du procédé de dessalement de l’eau de mer

par pompe à chaleur avec capteur solaire intégré

Mustapha Douani, Djilali Bouzina, Fatima Zahra Rahma, N. Tahraoui Douma

To cite this version:

Mustapha Douani, Djilali Bouzina, Fatima Zahra Rahma, N. Tahraoui Douma. Simulation
du procédé de dessalement de l’eau de mer par pompe à chaleur avec capteur solaire intégré.
Jean-Jacques BEZIAN. JITH 2007, Aug 2007, Albi, France. ENSTIMAC, 5p., 2007. <hal-
00164644>

HAL Id: hal-00164644

[Link]
Submitted on 29 Aug 2007

HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est

archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents
entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non,
lished or not. The documents may come from émanant des établissements d’enseignement et de
teaching and research institutions in France or recherche français ou étrangers, des laboratoires
abroad, or from public or private research centers. publics ou privés.
 13èmes Journées Internationales de Thermique

SIMULATION DU PROCEDE DE DESSALEMENT DE L’EAU

DE MER PAR POMPE A CHALEUR AVEC CAPTEUR SOLAIRE
INTEGRE.
Mustapha DOUANI, N. Djilali BOUZINA., [Link]., et [Link] DOUMA
Département de génie des procédés, Univ. Hassiba Ben Bouali, Chlef. Algérie
douani_mustapha@[Link].

Résumé : Le dessalement de l’eau de mer nécessite une dépense importante d’énergie afin
d’éliminer le soluté (eau). La pompe à chaleur assistée de capteurs solaires est en mesure de
répondre favorablement en termes d’économie d’énergie afin d’éliminer le soluté (eau) tout en
palliant les difficultés inhérentes aux conditions de leur fonctionnement (souplesse et adaptation
aux fluctuations saisonnières, etc.). Le fluide R-113 est retenu comme fluide caloporteur pour le
fonctionnement d’une pompe à chaleur à compression. La simulation par la méthode de Newton-
Raphson montre que le cycle conduit à des valeurs de COP dépassant la valeur de six (06) avec
un capteur solaire de rendement avoisinant 50%. L’influence des variables opératoires de
conception (Tcond, G, (Tpm-Tpa), Me) sur la performance de l’installation est également présentée.

Mots clés : Pompe à chaleur à compression, capteur solaire, simulation.

1. INTRODUCTION
Pour répondre à ses besoins hydriques, l’homme s’est concentré autour des rivières, des
réservoirs d’eau souterraine. L’eau va devenir un formidable enjeu durant le siècle à venir [1].
L’Algérie est amenée en plus soit à traiter et régénérer ses eaux usées et résiduaires industrielles
soit à dessaler de l’eau de mer [2,3]. De par sa nature, l’évaporation nécessite de grandes
quantités d'énergie et toute valorisation de sources locales d’énergie disponible est à exploiter
[4,5]. Le couplage solaire-pompe à chaleur est une des idées attrayantes. En effet, pour
l’exploiter, certains procédés physiques ou chimiques et mécaniques sont à mettre en œuvre [6].

2. BASES THEORIQUES ET TECHNIQUES

On montre que l’objectif de Compresseur de vapeur
l’implantation de la pompe à chaleur d’eau.
(PAC) est la valorisation de la quantité Eau salée +
de chaleur rejetée à bas niveau cristaux de sel S Eau salée
thermique TEv. En effet, le condenseur 6 5
Cri
de la pompe à chaleur est utilisé sta 4 Rayonnement
comme un cristalliseur tout en Détendeur incident.
exploitant la puissance fournie à la
source chaude sous forme de chaleur SC Flux
7
latente de condensation. Parmi les perdu.
configurations possibles, l’évaporateur
est conçu de sorte que l’énergie solaire
Evaporateur
incidente est à l’origine de la Compresseur
vaporisation du fluide réfrigérant: c’est Condensât
le capteur solaire (figure 1). Figure 1. Schéma de couplage pompe
à chaleur- capteur solaire.
Albi, France du 28 au 30 Août 2007 1
 13èmes Journées Internationales de Thermique

Le choix du fluide caloporteur est dicté par les critères d’Unterleutner [7] .Le concept du
coefficient de performance (COP) est introduit pour analyser la performance du cycle ; il est
défini par le rapport:
Effet utile de chauffage.
C.O.P = (1)
Energie nette fournie d ' un milieu exterieur.
où l’effet utile de cristallisation (Qcond) est la chaleur récupérée au condenseur. Le calcul des
propriétés thermophysiques de R113 (enthalpie, entropie, fugacité, volumes molaires liquide et
vapeur, etc.) nécessite une analyse thermodynamique faisant appel à l’équation d'état de Martin et
Hou. [8] Les sources de pertes d’exergie sont:
- Transfert de chaleur dû à un pincement fini de température δT [9].
- Pertes de charges au niveau des éléments constitutifs de l’installation.

2.1. Bilans de transfert et concept de degrés de liberté des différents éléments

Compte tenu du mode général du fonctionnement d’un capteur plan, le bilan thermique sur la
paroi absorbante s'écrit: ϕ sa = ϕ p + ϕ u + ϕ st (2)

En régime stationnaire, ϕ st est nul. Par ailleurs, la littérature abonde en matière de définitions de
la performance du capteur solaire et selon l'objectif, on définit le rendement global par la relation
ϕ
suivante [9] : η= *
u
(3)
G ( i ,γ ) *S
La puissance absorbée par le capteur s'écrit :
ϕ sa = τ cs × α ps × G(*i ,γ ) × S
(4)
Dans le cas où le fluide réfrigérant subit uniquement un changement d'état, le flux utile est
proportionnel à son enthalpie de vaporisation. Pour un capteur solaire de type II, les pertes
thermiques sont la somme des flux dissipés au niveau du haut et du bas, d’où:
 
 
 Tpm − Ta Tpm − Ta 
ϕp = ϕp + ϕp = 
↑ ↓
+ S (5)
1 1 ei 1 
 + + 
 hc, p −c + hr , p −c hvent + hr , a −c λi hvent 
L’application de la méthode de calcul de degrés de liberté pour chaque élément permet de
déterminer toutes les variables de conception, contrôlables extérieurement [10]. Si C désigne le
nombre de constituants d’un courant, alors on peut lui associer (C+3) variables qui peuvent être,
la composition (C), le débit L, la pression P et la température T. On note qu’une analyse détaillée
montre que le nombre de degré de liberté de l’unité est :
(ND )unité = C2 + 13 (6)

2.2. Simulation du fonctionnement du système de dessalement

Albi, France du 28 au 30 Août 2007 2

 13èmes Journées Internationales de Thermique

Après identification des paramètres opératoires par l’application du concept de degrés de

liberté, les étapes de la simulation utilisent un algorithme afin de déterminer:
• Le débit d’eau à traiter.
• Les températures Tpm et Tcm du capteur.
• La composition du mélange liquide-vapeur dans le circuit du réfrigérant.
• Le COP de la [Link] le rendement du capteur solaire ainsi que sa surface.

3. RESULTATS ET ANALYSE.

3.1. Influence de Tcond de R113 sur la performance technique du capteur.

Ce paramètre est particulièrement important, Qcond, dépend exclusivement de la température

au condenseur Tcond et que des conditions de la cinétique d’évaporation. On note une
décroissance du COP et du rendement du capteur avec Tcond laquelle s’explique par le fait que
l’augmentation de Tcond conduit inévitablement à l’augmentation des irréversibilités dans le cycle.
Cette influence fait l’objet de la figure n°2. On note que Tcond est inversement proportionnelle à
Qcond. Comparativement, on note que la pompe à chaleur à absorption donne des COP
relativement plus faibles (~1.37) [11].

3.2. Influence du débit du réfrigérant sur l'efficacité du cristalliseur.

On remarque que la quantité d'eau à traiter par unité du temps dépend tout naturellement du
débit de fluide réfrigérant. A travers la figure n°3, on voit clairement que cette dépendance est
quasiment linéaire; or cette croissance affecte la taille de l'installation et par conséquent, le coût
de l'investissement. Il ressort que la définition de la rentabilité économique impose la quantité
d'eau salée à traiter. Il y a lieu de préciser que cette dernière dépend des conditions de
fonctionnement et du flux utile (Qév) fourni par le capteur pour des COP constants.
2

20
Ames(Kg/m jour)*10

Rend [%-] 550

18 COP[--]
525
16
2

500
14
475
12
10 450

8 425
6 400
4 375
2 350
355 360 365 370 375 380 385 390 4,5 5 5,5 6 6,5 7
Tcond (K) mR(Kg/sec)*10
3

Figure 2. Influence de la température Tcond Figure 3. Variation du débit de l'eau en

sur la performance d'un système. fonction du débit du fluide réfrigérant.
T ev = 50 °C , T pa = 30 °C m R = 18 kg / h , T cm = 40 °C ,
x = 36 . 3 gNaCl / 100 g d ' eau , θ = 45 °, G = 1000 W / m 2 .

3.3. Influence du débit du réfrigérant sur le rendement et le débit de l'eau traitée.

Compte tenu de l'interaction énergétique entre les différents éléments constitutifs de

l'installation, nous avons jugé utile de présenter graphiquement la variation de la quantité d'eau

Albi, France du 28 au 30 Août 2007 3

 13èmes Journées Internationales de Thermique

traitée ainsi que le rendement du capteur solaire en fonction du débit du fluide réfrigérant,
laquelle fait l'objet de la figure n°4. On constate que le débit d'eau salée et le rendement sont
directement proportionnels au débit du caloporteur. Un tel résultat est parfaitement confirmé par
l'analyse fonctionnelle des relations de bilans thermiques sur le condenseur (cristalliseur ) et sur
le capteur solaire.
η[%]
2,50 60
2,25 50
2,00 40
1,75 30
1,50 20
Ames*1E+2(Kg
+1
1,25 η*10 [%] 10
1,00 0
4,5 5,0 5,5 6,0 6,5 7,0
300 400 500 600 700 800 900 1000 1100
+3
mR*10 (Kg/sec) G (watt/m2)

Figure 4. Influence du débit du réfrigérant Figure 5. Variation du rendement du capteur

sur la performance du système. en fonction de la densité de flux solaire.
T ev = 50 ° C , T pa = 30 ° C T cond = 85 ° C , Tcm = 40 ° C ,
x = 36 . 3 gNaCl / 100 g d ' eau , θ = 45 ° , G = 1000 W / m 2 .

3.4. Influence de l'intensité du flux solaire incident sur le rendement du capteur

L'énergie solaire G transmise à la surface terrestre par rayonnement est d'une intensité
maximale de l’ordre du kilowatt par mètre carré. Compte tenu de l'intérêt énergétique de
l’évolution diurne du flux solaire incident [12], nous avons estimé suffisant de présenter
l'influence de la variable G sur la performance du capteur à travers la figure n°5. Si le flux utile
est exclusivement dicté par les conditions de fonctionnement de la P.A.C, alors on a:
ϕ P = f (TPm − Tcm ) (7)
Alors lorsque ϕ ab passe du simple au triple 270
727
( )
. On note que ϕ p subit une croissance
exponentielle pour aller se quintupler. Il en découle que pour maintenir l'impact positif de G sur
la performance du capteur, il y a lieu d’augmenter sa surface. On précise que pour un pincement
(Tpm-Tcm)=10°C, le rendement est de ordre de 52%.

4. CONCLUSION:
Il ressort de cette analyse que le procédé de dessalement par pompe à chaleur assisté de
capteur solaire est particulièrement rentable sur le plan économique compte tenu du fait qu’on a
pu obtenir des COP de l’ordre de 6, avec un gain d’énergie qui atteindraient 500%. Par ailleurs,
certaines imperfections dues à l’échange thermique sont à l’origine de la création de flux
d’exergie. Les résultats essentiels obtenus mettent en évidence l'influence positive du débit du
fluide réfrigérant sur le rendement du capteur, la quantité d’eau salée à traiter. On doit noter que
la différence de température ( Tpm- Tcm) affecte négativement la surface du capteur. Cependant,
l’influence de la température à l’évaporateur joue en faveur d’une augmentation conséquente du
COP et du rendement du capteur solaire. Néanmoins, le rendement du capteur est assez sensible
lorsque l’angle d’inclinaison du capteur est de l’ordre de 40 degrés par rapport à l’horizontale.

Albi, France du 28 au 30 Août 2007 4

 13èmes Journées Internationales de Thermique

Nomenclature :

C2 : Nombre de constituants de l’eau salée (=2)

COP : Coefficient de performance Symboles grecs
* 2
G (i,γ) : Densité de flux solaire incident,W/m φst: Variation temporel de l'enthalpie de la masse
H: Enthalpie du réfrigérant, kcal/kg du fluide contenu dans le capteur
0 τcz : Coefficient de transmission de la couverture
m cf : Débit massique de l’eau salée, kg/sec transparente par rapport au rayonnement solaire
Mr : Débit du fluide frigorigène, kmol/h αpz : Coefficient d'absorption de la paroi
Q : Puissances calorifiques échangées, kcal/h
η: Rendement du capteur solaire
T : Température, K
x : Titre massique en NaCl , g. NaCl/ 100g.H2O

Références :
[1] [Link], et [Link], "Procédés thermiques de dessalement à haut rendement ",
Groupe pour la Recherche sur les Echangeurs Thermiques ([Link]); Commissariat à l'Energie
Atomique (C.E.A)/ Grenoble; n°44 - pp, 91, Hiver 2000-2001.
[2] J.V. Anderson, J.W. Mitchell et W.A. Beckman., "A design method for parallel solar heat
pump systems" , Solar Energy. 25 (2): pp.155-163, 1980.
[3] A.R. Day, T.G. Karayiannis, "Solar-assisted heat pump research and development", Building
Ser. [Link]. Technol. 15(2): pp. 71-80, 1994.
[4] E. Dietrich., P. Le Goff, "Pompe à chaleur à absorption avec stockage incorporé de l’exergie,
Evaluations énergétique et économique". Revue Générale de Thermique, n° 293, pp.273-287, mai
1986.
[5] M.N.A. Hawlader, K. D. Prasanta, Diab Sufyan, Chan. Ying Chung, " Solar assisted heat
pump desalination system",Elsevier, Desalination 168, pp.49-54, 2004.
[6] R. Tchinda, [Link], " Simulation numérique des performances d'un distillateur solaire
fonctionnant en mode indirect ", African Journal of Science and Technology, (AJST), Vol. n° 5,
n°1, pp. 79-91, Juin 2004.
[7] Unterleutner., H. Huemer et F. Moser, "Search for optimal working fluids in high
temperature heat pumps", International Congress of Refrigeration, Paris, 1983.
[8] A. Robert, “ASHRAE Handbook.” Chapitre 19. « Thermo physical Properties of Refrigerents
» I-P. pp. 3-4, Edit. 1997.
[9] Duffie., J.A., Beckman. W.A, " Solar engineering of thermal processes", J. Wiley and Sons,
New York, 1980.
[10] Mooson Kwauk , "Specification of design variables", A. I. Ch. E. Journal, pp. 240-249, Juin
1956.
[11] M. Douani, ‘’ Contribution à l’étude du couplage séchoir-pompe à chaleur à absorption
fonctionnant avec le système eau-triéthylène glycol’’. Thèse de Doctorat., [Link]. Oct.
(1989).
[12] B. Bouchekima., et al. " Etude théorique et application pratique du distillateur solaire à film
capillaire «. Int. J. Therm. Sci., 39, pp. 442-459, 2000.

Albi, France du 28 au 30 Août 2007 5