Manuscrit auteur, publié dans "2ème Colloque International Francophone en Energétique et Mécanique, CIFEM 2012,

Ouagadougou : Burkina Faso (2012)"

CIFEM2012 – ART-1-48

ETUDE SIMULEE D’UN SYSTEME DE DESSALEMENT D’EAU DE MER ET

DE PRODUCTION DE FROID PAR THERMOFRIGOPOMPE COUPLEE A
DES PANNEAUX SOLAIRES

a,* a a a
Paul BYRNE , Jacques MIRIEL , Laurent SERRES , Redouane GHOUBALI
a
Université Européenne de Bretagne
Equipe MTRhéo – Laboratoire LGCGM - INSA de Rennes et Université de Rennes1
IUT Génie Civil – 3 rue du Clos Courtel – BP 90422 – 35704 Rennes Cedex 7

RÉSUMÉ

L’accès à l’eau potable est un des enjeux majeurs des prochaines decennies. La recherche sur les technologies de
dessalement à forte efficacité énergétique doit donc être particulièrement active. Cet article présente un système de
hal-00729502, version 1 - 7 Sep 2012

thermofrigopompe produisant simultanément du froid pour une enceinte réfrigérée et de la chaleur utilisée par une unité
de dessalement par distillation membranaire. L’énergie électrique est fournie par des panneaux photovoltaïques pour une
possible implantation en site isolé. Ce système a été modélisé à l’aide des logiciels Trnsys et EES. Les simulations
montrent des résultats de ratios de production d’eau potable assez faibles. Toutefois, les pistes d’amélioration de
l’efficacité du système s’avèrent nombreuses.

Mots Clés : Simulation, thermofrigopompe, froid, dessalement

NOMENCLATURE

Symboles : tfct temps de fonctionnement, s

Lettres latines :
2
A surface, m Lettres grecques :
B besoins thermiques, Wh χ tortuosité, -
C consommation, Wh δ épaisseur, m
COP coefficient de performance frigorifique, - ε porosité, -
-1 -1
λ
Cp chaleur spécifique, [Link] .K -1 -1
conductivité, W.m .K
csel fraction molaire en sel, - µ viscosité dynamique, Pa.s
d diamètre, m ρ masse volumique, kg.m
-3

Dv/a coefficient de diffusion de la vapeur dans l’air, -

ft facteur taille, - Indices / Exposants :
-1
h enthalpie, [Link] a air
-1 -2
Jv flux de vapeur, kg.s .m c chaud
-1 -2 -1
K perméabilité, kg.s .m .Pa cd condensation
Kn nombre de Knudsen, - e entrée
L longueur caractéristique, m el électrique
lpm libre parcours moyen, m ev évaporation
-1
M masse molaire, [Link] f froid
-1
m débit massique, kg.s fct fonctionnement
m masse, kg h hydraulique
Nu nombre de Nusselt, - la lame d’air
p pression, Pa m membrane
P puissance, W mol molécule
Pr nombre de Prandtl, - p pore
3 -1
Qv débit volumique, m .s s sortie
2 -1
R résistance thermique, m .K.W T totale
Re nombre de Reynolds, - v vapeur
-1 -1
Ru constante des gaz parfaits, [Link] .K scm surface chaude membrane
T température, °C sfm surface froide membrane

*
auteur correspondant : Tel: +33 2 23 23 42 97 ; Fax: +33 2 23 23 40 51
Adresse électronique : [Link]@[Link]

CIFEM1-2012 -1-
 mécanismes de la diffusion moléculaire. Ce type de
transfert est défini par une prépondérance des collisions
1. INTRODUCTION entre molécules par rapport aux collisions entre les
molécules et les parois internes de la structure poreuse
Deux des enjeux majeurs des prochaines décennies sont [9].
l’accès à l’eau potable, notamment en site isolé, et Le système de thermofrigopompe peut également être
l’utilisation rationnelle de l’énergie. Dans un but social et utilisé pour produire du froid positif pour la conservation
industriel, il parait donc vital de développer de nouvelles de denrées alimentaires ou le conditionnement d’air, qui
installations de production d’eau potable à forte efficacité représente un besoin important dans les régions proches
énergétique. Plusieurs systèmes étudiés dans la de l’équateur. En restant dans le cas d’un cycle
littérature scientifique sont dans cette optique et monoétagé, l’écart de température entre les sources
s’approchent du système étudié dans la présente froide et chaude impose de choisir un fluide frigorigène
publication. Par exemple, Cherif et al. propose une ayant une température critique suffisamment élevée
solution de couplage d’un champ photovoltaique, d’une (environ 130°C) pour réaliser le réchauffage de la
éolienne avec une unité de dessalement par osmose solution à des températures entre 40 et 70°C. Dans cette
inverse [1]. Slesarenko présente un système de étude, le fluide choisi est l’isobutane (R600a).
dessalement d’eau de mer par évaporation puis Le bloc compresseur de la machine frigorifique peut être
condensation à l’aide d’une pompe à chaleur [2]. Notre entraîné par un moteur à aimants permanents, alimenté
équipe de recherche travaille sur des procédés de en courant continu 12 ou 24 V par des panneaux
chauffage et de rafraîchissement simultanés par photovoltaïques associés à un système de batteries [10].
thermofrigopompe (TFP) sur les aspects de modélisation Les moteurs à aimants permanents sont aujourd’hui
[3], [4] et expérimental [5], [6]. Une thermofrigopompe l’objet d’un développement notable du fait de leur forte
(machine frigorifique produisant de la chaleur et du froid efficacité énergétique, de leurs faibles nuisances sur le
utiles) à haute température peut permettre de produire la plan acoustique et de leur longévité grâce à des forces
chaleur nécessaire au dessalement d’eau de mer ou plus de frottement quasi-nulles entre le rotor et le stator.
hal-00729502, version 1 - 7 Sep 2012

généralement à la distillation d’eau saumâtre par un Certaines petites installations de réfrigération utilisent
procédé membranaire et simultanément rafraîchir un déjà ce genre de moteur à courant continu pour entrainer
espace. Un lieu d’application d’un tel système serait par le compresseur. L’association d’une thermofrigopompe à
exemple un dispensaire dans lequel se trouveraient une des panneaux solaires photovoltaïques confèrera une
salle de soins rafraîchie ou une armoire réfrigérée et une autonomie de fonctionnement à ce système. Une telle
unité de production d’eau potable. installation combinant ces technologies innovantes peut
trouver des applications dans la marine, l’aide médicale
Apport de
d’urgence ou l’amélioration des conditions de vie des
chaleur populations en site isolé.
Cet article présente en détail la modélisation de l’unité de
Τs1 Plaque froide Membrane Τe2 dessalement à l’aide d’EES (Engineering Equation
Canal de Solver). Le logiciel EES est un solveur d’équations qui
condensation
dispose de fonctions internes de calcul des propriétés
Τcd Τsfm Τscm thermophysiques des fluides et des matériaux. La
δf modélisation des transferts couplés de chaleur et de
Fluide Fluide
froid chaud
masse du système membranaire est réalisée d’après les
Vapeur
équations de la publication de Alklaibi et Lior [11]. Le
modèle est validé par comparaison avec les résultats de
la thèse de Banat [12]. Des modèles simples de
δa thermofrigopompe (réalisé avec EES), de besoins de
δp δm froid, de panneaux photovoltaïques et de batteries de
Τe1 Eau Τs2
stockage sont couplés à l’unité de dessalement à l’aide
m& e potable m& s du logiciel Trnsys. Le modèle de thermofrigopompe est
Eau de mer Jv A Saumure tiré des caractéristiques d’un compresseur ATEX
disponible dans le commerce. Notre étude aboutit à
Fig. 1. Principe de la distillation membranaire l’évaluation des consommations électriques, de la
production d’eau potable et de la performance globale de
Différentes méthodes sont utilisées pour le dessalement ce système par rapport à d’autres systèmes présentés
[7] : la distillation, l’humidification-déshumidification d’air, dans la littérature.
l’osmose inverse, l’électrodialyse... La distillation
membranaire (figure 1), choisie pour cette étude, est
apparue à la fin des années 1960 grâce aux progrès 2. PRESENTATION DU MODELE
réalisés dans la fabrication de membranes [8]. Elle fait
partie des dernières évolutions du dessalement par 2.1. MODELE GLOBAL
procédé thermique. Cette technique repose sur le
principe de migration de vapeur d’eau à travers une Le schéma global de l’installation est présenté sur la
membrane microporeuse hydrophobe par différence de figure 2. Le système de production d’électricité est
pression de vapeur entre une solution surchauffée et l’air modélisé sous Trnsys à l’aide du type 194 pour les
contenu dans un canal de condensation. Ce dernier est panneaux photovoltaïques et du type 47a pour le
refroidi par la solution à température initiale ou par un système de stockage d’électricité. Un modèle de
autre milieu à faible température. Pour des membranes réfrigérateur placé dans un local a été réalisé à l’aide du
ayant des diamètres de pores de l’ordre du micromètre, type 56 (modèle de bâtiment multizone). Une enceinte
le transfert de masse à travers la membrane est 3
réfrigérée d’1 m constitue une première zone thermique
considéré comme étant géré principalement par les

CIFEM1-2012 -2-
 3
contenue dans un local de 27 m constituant une Les évolutions du coefficient de performance et des
seconde zone. Les parois de l’enceinte réfrigérée sont puissances thermiques et électrique sont présentées
modélisées par 5 cm de polystyrène expansé de respectivement sur les figures 3 et 4. Les valeurs sont
-1 -1
conductivité 0,039 W.m .K . Les caractéristiques des calculées pour une température d’évaporation de 0°C ,
parois du local sont présentées dans le tableau 1. Un une surchauffe de 5 K et un sous-refroidissement de 2 K.
2
double vitrage d’1 m de surface et de coefficient de La température d’évaporation est supposée constante
-2 -1
déperdition Uw égal à 2,95 W.m .K est placé sur la car l’évaporateur se situe à l’intérieur de l’enceinte du
façade sud du local. réfrigérateur modélisée comme étant parfaitement
régulée à une température de 5°C. De ce fait, lorsq ue la
Tableau 1. Caractéristiques des parois du local température de condensation varie de 45 à 75°C, les
Epaisseur Conductivité performances diminuent. L’écart entre la température de
Matériaux
(cm) (W.m-1.K-1) condensation et la température de sortie du fluide chaud
Enduit intérieur plâtre 1,5 0,35 correspondant à Te2 (pincement de l’échangeur) est
Bloc béton 20 1,05 supposé constant et égal à 5 K.

Polystyrène expansé 7 0,039

Enduit extérieur mortier 1,5 1,15

Un lecteur de données climatiques est connecté aux

modèles de panneaux photovoltaïques et au modèle de
bâtiment. Le climat choisi est celui d’Ouagadougou
(Burkina-Faso). Le reste des composants est modélisé à
l’aide du solveur d’équation EES : le compresseur et son
moteur électrique, le condenseur fournissant la chaleur
nécessaire au dessalement, le détendeur et l’évaporateur
hal-00729502, version 1 - 7 Sep 2012

fournissant de l’énergie frigorifique au modèle de

réfrigérateur.

Unité de
Fig. 3. Evolution du COP en fonction de la
dessalement
température de condensation de la TFP
Chaleur
Panneaux
photovoltaïques Condenseur

Moteur à
courant
continu

Evaporateur
Batteries
Froid

Réfrigérateur
+ Bâtiment

Fig. 2. Schéma du système de dessalement par TFP Fig. 4. Evolution des puissances en fonction de la
température de condensation de la TFP
2.2. MODELE DE TFP
2.3. MODELE D’UNITE DE DESSALEMENT
En l’état actuel de ce projet de recherche, le moteur à
courant continu n’a pas été modélisé. Les performances Le modèle d’unité de dessalement est bâti d’après l’étude
d’un compresseur semi-hermétique fonctionnant au expérimentale des travaux de thèse de Banat [12] et les
R600a ont été obtenues à partir du logiciel de équations de l’étude simulée sont tirées de la publication
dimensionnement d’un fabricant de compresseurs ATEX. d’Alklaibi et Lior [11] qui s’appuient sur ces mêmes
Dans la gamme de compresseurs disponible, le choix travaux de thèse. Le modèle est monodimensionnel dans
s’est porté sur celui de plus faible puissance. Le temps le sens du flux de vapeur et fonctionne en régime
de fonctionnement de la TFP sur chaque pas de temps permanent.
permet de gérer le surdimensionnement du compresseur
par rapport aux besoins en froid. Il permet de calculer la 2.3.1. Transfert de vapeur
consommation électrique du compresseur sur chaque Le flux de vapeur à travers la membrane est défini par la
pas de temps en fonction des besoins en froid d’après les loi de Stefan (équation (3)).
équations (1) et (2). J v = K ∆p v (3)
B (1) Celui-ci dépend de la perméabilité de la membrane K et
t fct = f de la différence de pression de vapeur saturante de part
Pf
et d’autre de celle-ci ∆Pv. Les pressions de vapeur
C el = Pel t fct (2) saturante sont calculées grâce à la fonction interne

CIFEM1-2012 -3-
 d’EES faisant appel à la base de données des propriétés A(Tsfm − Tcd )
thermophysiques des fluides. La loi de Raoult est utilisée m& eCp (Ts1 − Te1 ) = J v A∆hcd + (14)
Rla
pour prendre en compte l’effet de la présence de sel
dans la solution en contact avec la surface chaude de la Puis, le modèle de TFP délivre une puissance calorifique
membrane (équation (4)). Celle-ci dépend de la fraction afin d’augmenter la température de cette solution jusqu’à
molaire en sel csel. la consigne définie par l’utilisateur. On simule ici le
pv ,scm = (1 − c sel ) pv (Tscm ) (4) fonctionnement d’un échangeur parfait (Equation (15)).
m& e Cp (Te 2 − Ts1 ) = Pc (15)
La perméabilité de la membrane est donnée par
l’équation (5) dans laquelle δm (l’épaisseur) ε (la porosité) Le second canal dispose de la paroi poreuse. La solution
perd de la masse et de la chaleur à travers la membrane
et χ (la tortuosité) sont des constantes du matériau
(équation (16)).
A(Tscm − Tsfm )
membranaire, Dv/a est le coefficient de diffusion de la
 m& e + m& s 
vapeur dans l’air (équation (6)), Pa,moy et Ta,moy sont  Cp (Te 2 − Ts 2 ) = J v A∆hev + (16)
respectivement les moyennes de la pression partielle  2  RmT
d’air et de la température dans la membrane et PT est la Le bilan de masse est donné par l’équation 17.
pression totale. m& e = m& s + J v A (17)
εDv / a M v pT
K= (5)
χδ m pa ,moy R(Tm ,moy + 273,15) 2.3.3. Validation du modèle
Le modèle d’unité de dessalement que nous avons
λv
Dv / a = (6) développé est construit sur la base du prototype
ρ vCpv expérimental de la thèse de Banat [12]. La courbe de
Ce modèle de diffusion, de type moléculaire, n’est valide validation du flux surfacique de vapeur en fonction de la
que pour des valeurs du nombre de Knudsen Kn température d’entrée d’eau chaude se trouve sur la figure
(équation (7)) inférieures à 0,01 [12]. Ce nombre est le 5. Le modèle est en très bonne concordance avec les
rapport du libre parcours moyen lpm (équation (8)) sur le résultats des travaux de Banat. Les valeurs minimales et
hal-00729502, version 1 - 7 Sep 2012

diamètre des pores dp. Dans notre étude, le libre maximales (Banat mini et Banat maxi) sont celles d’une
parcours moyen est inférieur à 90 nm. Ceci impose donc série de tests de reproductibilité réalisés pour plusieurs
un diamètre de pore supérieur à 9 µm. membranes de même facture avec de l’eau pure. Les
l pm valeurs expérimentales (Banat exp) correspondent aux
Kn = (7) essais avec de l’eau salée. La salinité de la solution en
dp entrée est de 35 g/kg. Celle-ci a peu d’influence sur le
Ru Tm ,moy flux de vapeur.
l pm = (8)
π 2d mol
2
pT

2.3.2. Transfert de chaleur

Le fluide chaud transfère une partie de sa chaleur au
fluide froid à travers la membrane, la lame d’air, le film
liquide et la plaque froide. L’autre partie de la chaleur
cédée est contenue dans la masse d’eau potable
produite, évacuée du système. La résistance thermique à
travers la membrane RmT (équation (9)) dépend de deux
flux de chaleur identifiables. Celui passant par la partie
solide de la structure poreuse (correspondant à la
résistance Rm de l’équation (10)) et la chaleur contenue
dans le flux de vapeur (Rv de l’équation (11)).
Rm Rv
RmT = (9) Fig. 5. Courbe de validation du modèle d’unité de
Rm + Rv dessalement
δm
Rm = (10)
ελa + (1 − ε )λm Les dimensions de la cellule de distillation membranaire
modélisée sont répertoriées dans le tableau 2.
1
Rv = (11)
J vCpv Tableau 2. Dimensions de la cellule
A travers la lame d’air, il est supposé que le flux convectif Longueur 16 cm
est totalement porté par la vapeur d’eau et le flux de Largeur 10 cm
chaleur par conduction se produit uniquement à travers la Epaisseur du canal chaud 1/16 pouce = 1,5875 mm
couche d’air. La résistance correspondante est calculée
grâce aux équations (11), (12) et (13). Epaisseur de la membrane 0,11 mm
Ra Rv Epaisseur du canal de condensation 10 mm
Rla = (12)
Ra + Rv Epaisseur de la plaque froide 1,5 mm
δa Epaisseur du canal froid 1/16 pouce = 1,5875 mm
Ra = (13)
λa
Suivant le schéma de la figure 1, la solution saline subit, Les caractéristiques des matériaux de la membrane et de
dans le premier canal, un réchauffement régi par la plaque froide sont données dans le tableau 3. Dans la
l’équation (14). simulation, la membrane est supposée être en

CIFEM1-2012 -4-
 polyfluorure de vinylidène (PVDF) et la plaque froide est en diminuant le débit d’alimentation de la cellule car
en cuivre. d’après l’équation (15), la puissance calorifique est
directement proportionnelle au débit d’eau salée. Le ratio
Tableau 3. Caractéristiques des matériaux de production d’eau potable est donc plus sensible à
Porosité de la membrane 0,75 l’évolution de la puissance calorifique et électrique de la
Tortuosité de la membrane 1,1 TFP qu’à celle du flux de vapeur.
Conductivité de la membrane 0,18 W.m-1.K-1
Conductivité de la plaque froide 400 W.m-1.K-1

2.3.3. Adaptation du modèle au cadre de l’étude

Afin d’augmenter les performances de l’unité de
dessalement modélisée, la longueur de la cellule est
portée à 2 m, une dimension adapté à une installation
dans un local technique. Le débit utilisé dans la thèse de
Banat est de 5,5 l/min. Dans le cadre de cette étude, un
débit de 3,2 l/min est également modélisé afin
d’augmenter le rapport de la masse d’eau potable
produite par unité d’énergie électrique consommée tout
en restant dans les limites du critère de validité de la
corrélation utilisée pour le calcul du nombre de Nusselt
en régime laminaire. Le critère de validité est donné par Fig. 6. Evolution du coefficient de performance en
l’équation (18) et la corrélation par l’équation (19). fonction de la température d’entrée du fluide chaud
dh
RePr > 10 (18)
L
hal-00729502, version 1 - 7 Sep 2012

0 ,14
d   µ moy 
13

Nu = 1,86 [Link]. h   (19)
 L   µ paroi 
Un facteur de taille est appliqué au système de
dessalement afin d’adapter la puissance nécessaire au
modèle d’unité de dessalement Pc-ref à la puissance
disponible au condenseur Pc de la TFP suivant l’équation
(20).
P (20)
ft = c
Pc − ref
Le facteur de taille et le temps de fonctionnement de la
TFP rentrent en compte dans le calcul de la masse d’eau
potable produite sur chaque pas de temps de simulation Fig. 7. Evolution du ratio de production d’eau potable
(équation (21)). en fonction de la température d’entrée du fluide
meau = J v Af t t fct (21) chaud

Enfin, une simulation annuelle menée avec le climat

3. RESULTATS d’Ouagadougou est présentée sur la figure 8. L’évolution
de la production d’eau potable suit les évolutions des
L’évolution du COP frigorifique et du flux de vapeur en températures extérieure et intérieure. Les besoins de
fonction de la température d’entrée du fluide chaud est froid du réfrigérateur augmentent avec la température
présentée sur la figure 6. Lorsque la température intérieure et le temps de fonctionnement de l’unité de
augmente, le COP diminue du fait de l’augmentation de dessalement augmente en conséquence. La production
l’écart de pressions de fonctionnement. Cette courbe est totale annuelle d’eau potable atteint 108 kg pour une
identique à celle de la figure 3. Il est néanmoins consommation électrique égale à 574 kWh. Le ratio de
intéressant de la confronter à l’évolution du flux de consommation revient à 5320 kWh/m . L’énergie
3

vapeur. Le flux de vapeur augmente avec la température électrique sert à produire de la chaleur pour le
du fluide chaud. Plus cette température est élevée, plus dessalement mais avant tout à produire du froid. On
les écarts de température et de pressions de vapeur pourrait croire que la chaleur produite est gratuite par
saturante entre les deux surfaces de membrane sont rapport à un réfrigérateur classique. Cependant, il faut
élevés. Le flux de vapeur est également sensible au débit garder en mémoire que la performance du système de
d’alimentation en eau de l’unité Qv. Lorsque le débit TFP pour le dessalement est inférieure à celle d’un
diminue, l’écart des températures moyennes sur la réfrigérateur du fait de la température de source chaude
hauteur de la cellule diminue et le flux de vapeur généré plus élevée. Par ailleurs, un autre point de vue consiste à
diminue également. considérer que le système de production d’électricité est
Les effets de la diminution du COP et de l’augmentation solaire et donc que l’énergie est gratuite. Pour autant,
du flux de vapeur avec la température d’entrée du fluide nous ne devons pas oublier de progresser sur le terrain
chaud se compensent pour donner un ratio de production de l’efficacité énergétique. Etant actuellement aux
d’eau potable quasi-constant (figure 7). Ce ratio est défini prémices du projet, aucun système de récupération de
comme étant le rapport de la masse d’eau produite par chaleur ou d’économie d’énergie n’a encore été mis en
unité d’énergie électrique consommée. Il est calculé sur place dans le modèle de simulation. Si bien que si le
la première journée de simulation. Sa valeur augmente compresseur était connecté au réseau électrique (et non

CIFEM1-2012 -5-
 à des panneaux solaires) avec un tarif d’énergie de 0,01 pistes d’étude paraissent intéressantes : le
€/kWh, le coût de revient de la production d’eau potable comportement du système avec d’autres types de
3
serait de 532 €/m . Au bilan, les ratios de consommation besoins en froid tels que le rafraîchissement d’espaces,
et de coût ne sont pas satisfaisants. L’objectif que nous la mise en place d’un échangeur récupérateur entre les
ambitionnons est de tendre vers des valeurs de ratios fluides chaud et froid ou d’une recirculation (avec
largement inférieures telles que celles publiées dans la mélange ou non) de la solution saline du fait de la faible
littérature. Alklaibi et Lior [8] présentent des ratios de sensibilité du système à la salinité [8]. L’étude d’autres
3
consommations de 1,25 kWh/m sur la base d’une cycles frigorifiques tels que le cycle transcritique du CO2
récupération de la chaleur d’un procédé industriel. En ou des cycles bi-étagés doit également être réalisée afin
termes de coût, Charcosset [7] cite une évaluation d’optimiser les performances de la TFP.
3
d’environ 11 à 14 €/m avec une production de chaleur
par énergie solaire en incluant le coût de l’équipement.
REFERENCES

[1] Cherif H, Belhadj J. Large-scale time evaluation for

energy estimation of stand-alone hybrid photovoltaic-wind
system feeding a reverse osmosis desalination unit.
Energy 2011;36:6058-6067.
[2] Slesarenko VV. Heat pumps as a source of heat
energy for desalination of seawater. Desalination
2001;139:405-410.
[3] Byrne P, Miriel J, Lénat Y. Design and simulation of a
heat pump for simultaneous heating and cooling using
HFC or CO2 as a working fluid. International Journal of
Refrigeration 2009;32:1711-1723.
hal-00729502, version 1 - 7 Sep 2012

[4] Byrne P, Miriel J, Lénat Y. Modelling and simulation of

a heat pump for simultaneous heating and cooling.
Fig. 8. Evolution de la masse d’eau produite et des Building Simulation: An International Journal 2012;5
températures intérieure et extérieure lors d’une [5] Byrne P, Miriel J, Lénat Y. Experimental study of a
simulation annuelle heat pump for simultaneous heating and cooling – Part 1:
Basic concepts and performance verification. Applied
Energy 2011;88:1841-1847.
4. CONCLUSION [6] Byrne P, Miriel J, Lénat Y. Experimental study of a
heat pump for simultaneous heating and cooling – Part 2:
Cette étude présente un système de dessalement couplé Dynamic behaviour and two-phase thermosiphon
à une thermofrigopompe alimentée par des panneaux defrosting technique. Applied Energy 2011;88:3072–
photovoltaïques. La production de froid par la TFP 3078.
conditionne la production d’eau potable. Ce système a [7] Charcosset C. A review of membrane processes and
été modélisé à l’aide d’un couplage entre les logiciels renewable energies for desalination. Desalination
Trnsys et EES. L’effort de développement s’est porté 2009;245:214–231.
principalement sur l’unité de dessalement par distillation [8] Alklaibi AM, Lior N. Membrane-distillation desalination:
membranaire. Une bonne correspondance entre les status and potential. Desalination 2004;171:111-131.
résultats de simulation et les résultats expérimentaux de [9] Loussif N, Orfi J, Sene M, Sow O, Maré T, Adj M.
la thèse de Banat a été obtenue. L’étude simulée montre Transferts couplés de chaleur et de masse dans une
que lorsque le débit d’alimentation de l’unité augmente, unité de dessalement par distillation membranaire. Actes
le flux de vapeur généré augmente mais le ratio de la du IXème Colloque Interunversitaire Franco-Québécois,
production d’eau potable par unité d’énergie électrique Lille, France, 2009.
consommée diminue. A l’avenir, les limites de l’étude [10] Axaopoulos PJ, Theodoridis MP. Design and
doivent être repoussées en termes de débit experimental performance of a PV Ice-maker without
d’alimentation et de caractéristiques de membrane par la battery. Solar Energy 2009;83:1360–1369.
réalisation d’un prototype afin d’améliorer les [11] Alklaibi AM, Lior N. Transport analysis of air-gap
performances du système jusqu’à atteindre des niveaux membrane distillation. Journal of Membrane Science
de ratios de consommation et de coût du même ordre de 2005;255:239-253.
grandeur que ceux des autres systèmes de dessalement [12] Banat FA. Membrane distillation for desalination and
(distillation membranaire couplée à l’énergie solaire, removal of volatile organic compounds from water.
osmose inverse, humidification-déshumidification…). En Doctoral Thesis, McGill University, Montreal, Canada;
l’état actuel de ce projet de recherche, le moteur à 1994.
courant continu n’a pas été modélisé. La mise au point
du prototype permettra de valider cette option. D’autres

CIFEM1-2012 -6-