Science et Environnement 34 (2018) 74-83 74

ISSN 2409 -6245, [Link]

Dessalement de l’eau de mer par voie thermique solaire à

Djibouti
H. Osman Alia, A. Khamlichib, H. Ali Barkadc

a. Département de Physique, Faculté des Sciences de Tetouan, Tetouan 93030, Morocco,

IUT-I, Université de Djibouti, Djbouti, Djbouti, arhiba@[Link]
b. Département de Physique, Faculté des Sciences de Tetouan, Tetouan 93030, Morocco,
khamlichi7@[Link]
c. IUT-I, Université de Djibouti, Djibouti, Djibouti, hassan_ali_barkad@[Link]

Correspondance, courriel: arhiba@[Link]

Résumé :

Le manque d’eau douce représente un grand défi pour la République Djibouti. Face à ce
problème, le dessalement de l’eau de mer apparaît comme une solution très opportune. La
technique de dessalement par voie thermique solaire en provoquant l’évaporation de l’eau
saumâtre avant de la condenser pour récupérer l’eau douce est une solution qui semble être
très prometteuse. Dans le cadre de ce travail, l’attention est focalisée sur le distillateur actif
qui utilise un apport additionnel d’énergie thermique au moyen d’un capteur solaire
indépendant. Notre objectif est d’aboutir à la conception optimale de ce système pour
produire 100 m3 d’eau douce par jour à base de la simulation numérique. Nous avons pour
cela écrit les équations différentielles qui gouvernent les différents échanges thermiques se
produisant dans l’installation. Nous avons ensuite intégré ce système numériquement à l’aide
d’un schéma approprié. Une étude paramétrique a été menée en vue de déterminer la
conception optimale.
Mots clés: Dessalement, Distillateur, Thermique solaire, Evaporation, Optimisation

1 Introduction

Dans de nombreuses régions arides du monde à l’instar de Djibouti et du Maroc, les sources
d’eau douce sont inexistantes ou deviennent insuffisantes au regard de la croissance
démographique ou de la production industrielle. La République de Djibouti dispose d’un
important potentiel en énergies renouvelables. Vu sa position géographique, le climat est
tropical aride, enregistrant une faible pluviométrie et connaissant une forte évaporation en
raison de la chaleur : Djibouti manque d’eau douce. Avec l’important gisement solaire dont
Djibouti profite, il apparait opportun de s’intéresser à la technique thermique par évaporation
qui représente une solution prometteuse afin de subvenir à une partie des besoins en eau. D’où
l’intérêt du dessalement comme procédé durable de production de l’eau potable. Les enjeux
modernes relatifs à l’activité de dessalement de l’eau de mer tournent essentiellement autour
de la réduction du coût énergétique qui reste encore élevé, et particulièrement en cas de faible
besoin: consommation ne dépassant pas 100m3 / jour [1]. Pour satisfaire la demande en eau

H Osman Ali
 Science et Environnement 34 (2018) 74-83 75

potable à faible échelle, les distillateurs solaires passifs ont fait leur preuve depuis la première
réalisation pratique dans ce domaine qui remonte à la fin du dix-neuvième siècle : las Salinas
au Chili [2]. Plusieurs pays de tiers-monde utilisent les distillateurs solaires pour la production
d’eau potable à partir d’eau salée. Malik et al. [3] ont montré que la distillation solaire est la
méthode la plus économique pour la production d’eau douce en petites quantités. Dans ce
travail, nous étudions un distillateur solaire de type actif. Nous considérons que l’eau de mer
reçoit une quantité de chaleur supplémentaire de la part du fluide caloporteur qui circule en
bas d’une plaque en cuivre disposée au fond de la cuve. C’est dans le contexte du dessalement
à faible échelle que s’inscrit ce travail où le dessalement de l’eau saumâtre par voie thermique
est considéré via l’utilisation de l’énergie solaire comme unique source de chaleur. Les
équations du système ont été obtenues en effectuant le bilan énergétique pour les différents
éléments composant le distillateur [4]. Le modèle obtenu est résolu numériquement à l’aide
du logiciel Matlab. Le modèle complet ainsi construit permet de rendre compte de manière
explicite des nombreux paramètres en jeu. Ceux-ci comprennent les paramètres de
conception (longueur du distillateur, inclinaison de le la vitre, rapport d’aspect, longueur du
capteur solaire), les paramètres de fonctionnement (débit de l’eau de mer, débit du fluide
caloporteur, épaisseur de l’eau de mer dans le bassin), les paramètres environnementaux
(humidité relative, vitesse du vent, température de l’eau de mer, état du ciel) et les paramètres
qui conditionnent la puissance du rayonnement solaire capté. L’intégration de ces équations a
permis de conduire des études paramétriques extensives dans le but de déterminer l’influence
des différents facteurs qui interviennent dans le problème.

2 Description et modélisation du système

Le système de dessalement comprend deux circuits: le circuit primaire qui apporte la chaleur
et le circuit secondaire dans lequel circule l’eau saumâtre à évaporer. Le fluide caloporteur
chauffe à travers une paroi métallique en cuivre l’eau saumâtre à distiller. L’enveloppe en
verre sert à la fois d’apport en rayonnement solaire qui peut se rajouter à la quantité de
chaleur fournie par le circuit primaire. Cette enveloppe qui est au contact de l'air froide
permet de condenser la vapeur d'eau qui peut ruisseler sur la face intérieure et peut ainsi être
récupérée. La figure 1 présente le schéma de principe de cette installation.

H Osman Ali
 Science et Environnement 34 (2018) 74-83 76

Figure 1: Configuration d’un distillateur solaire actif [4]

Le modèle mathématique du système intégré peut être développé à base de bilan énergétique [4].

- Vitre extérieure

v Vv cp,v dTve  v
  Tvi  Tve   hc ve _ amb  Tve  Tamb   hrve _ ciel  Tve  Tciel  (1)
2A v dt ev

- Vitre intérieure

v Vv cp,v dTvi  h 
2A v dt
  m
 v I(t)  hcem _ vi  hrem _ vi  Tem  Tvi   ev v  v  Tvi  Tve 
Av ev
(2)

- Eau de mer

dTem
ρ em Vem c p,em =αem A em I(t) + hcc_em A em  Tc -Tem  + m  em cp,em Tem,en - m
 s cp,s Ts,sor
dt (3)
 ev h L - hcem_vi A em  Tem -Tvi  - hrem_vi A em  Tem -Tvi 
-m

- Plaque de cuivre

c Vc cp,c dTc
 c I(t)  hcfc _ c  Tfc  Tc   hcc _ em  Tc  Tem  (4)
Ac dt

- Fluide caloporteur

fc Vfc cp,fc dTfc 2fc Dfc cp,fc

Afc dt

Afc
Tfc,en  Tfc   hcfc _ c Tfc  Tc   Uperte Tfc  Tamb  (5)

H Osman Ali
 Science et Environnement 34 (2018) 74-83 77

où hc ve_amb est le coefficient d’échange par convection avec l’air ambiant, hrve_ciel le coefficient
d’échange par rayonnement avec le ciel, hcc_ em le coefficient de convection entre la plaque de
cuivre et l’eau de mer, hcem_vi le coefficient de convection entre l’eau de mer et la vitre
interne, hrem_ vi le coefficient d’échange par rayonnement entre l’eau de mer et la vitre interne,
hcfcc est le coefficient d’échange par convection entre le fluide caloporteur et le cuivre et
U perte le coefficient de convection de perte. Les coefficients de transferts de
chaleur peuvent être évalués par des corrélations précises. Ils sont donnés dans [4].

La mase d’eau distillée est donnée par :

tc
 m
 ev (t)dt (6)
0

où t c est la durée pour calculer la masse d’eau distillée.

3 Résultats et Discussion

La figure 2 montre le modèle numérique du système et précise les paramètres d’entrée et de

sortie. Tous les paramètres de conceptions ainsi que ceux du fluide caloporteur sont inclus
dans le modèle.

Vitesse Humidité Température Rayonnement

du vent relative ambiante Solaire

Température d’entrée de
Températures
l’eau de mer « Tem,en »
des sorties
Débit volumique de l’eau
de mer « Dem »
DISTILLATEUR SOLAIRE ACTIF
Masse d’eau
Epaisseur de l’eau de mer douce

Température Débit volumique

d’entrée du fluide du fluide
caloporteur « Tfc » caloporteur « Dfc »

Figure 2: Schéma du modèle numérique

 Nous simulons les évolutions des températures de chaque élément du système et nous
calculons la production de l’eau douce du distillateur.

H Osman Ali
 Science et Environnement 34 (2018) 74-83 78

360 10

9
350

debit massique de vapeur d eau en(kg/h)

8
340
7
Temperatures en (K)

330 6

5
320

Tvi 4
310 Tem
3
Tc
Tfc
300 2

1
290
0 5000 10000 15000
Temps en (s) 0
0 5000 10000 15000
Temps en (s)
Figure 3: Variation temporelle des différentes
températures des éléments du distillateur Figure 4 : Variation temporelle du débit
solaire actif produit d’eau distillée

La figure 3 montre l’évolution des différentes températures des éléments du distillateur solaire
actif de type four solaire. La température de fluide caloporteur est la température motrice de
notre système de dessalement. Dans le régime transitoire, la température du fluide caloporteur
atteint sa valeur minimale à l’instant 363s après le démarrage. Durant ce temps, la
température de la plaque de cuivre croit brusquement du fait de ses caractéristiques,
notamment l’épaisseur très faible, la couleur noire de la peinture et la conductivité thermique
très élevée. Par contre les températures de l’eau de mer et de la partie interne de la vitre
augmentent progressivement. Après 363s toutes les températures des différents composants
du système tendent à se stabiliser. La remontée de température du fluide caloporteur
s’explique par la circulation permanente du fluide caloporteur. Le système finit par atteindre
un état d’équilibre thermique au bout une durée de temps environ 5000s (1h23min).

La figure 4 décrit l’évolution du débit de l’eau distillée en fonction du temps. On voit bien
que la variation temporelle du débit d’eau distillée prend la même allure que celles des
températures de l’eau de mer et de la vitre interne: le débit de l’eau distillée augmente
progressivement et atteint le régime stationnaire après une durée de temps d’environ 5000s.

 Plusieurs simulations numériques ont été effectuées afin d’évaluer l’influence des
différents facteurs sur masse d’eau distillée  .

H Osman Ali
 Science et Environnement 34 (2018) 74-83 79

autres: < 1%

dem: 46%
Tfc: 51%

Dfc: 2%

Figure 5 : Effets de facteurs sur la masse d’eau distillée

Dfc [1.4;2.8] m3 /h , Tfc [353.15,363.15] K et d em [0.02;0.06] (m).

1.8

1.6

1.4
Masse d eau distillée (kg)

1.2

0.8

0.6

0.4

0.2

0
1.4 1.6 1.8 2 2.2 2.4 2.6 2.8
débit fluide caloporteur(m3/h)

Figure 6 : Variation de la production d’eau distillée en fonction du débit de fluide caloporteur

La figure 5 montre bien évidemment que la température du fluide caloporteur est le facteur
qui influence le plus le processus actif de production de l’eau douce par le distillateur solaire.
Cette température est suivie par le facteur d’épaisseur d’eau de mer dans la cuve
d’évaporation qui est aussi un facteur pesant dans un système du distillateur solaire.
La figure 6 montre que la masse d’eau distillée augmente légèrement avec l’augmentation du
débit de fluide caloporteur dans l’intervalle d’étude considéré.

H Osman Ali
 Science et Environnement 34 (2018) 74-83 80

300

masse d eau distllée(kg)

250

200

150

100

50
7
6 0.06
5 0.05
0.04
4 0.03
longueur du distillateur(m) 3 0.02
profondeur d eau de mer(m)

Figure 7: Evolution de la production de masse d’eau distillée en fonction de la profondeur de

l’eau de mer dans la cuve d’évaporation et la longueur du distillateur

210

200
Masse d eau distillée(kg)

190

180

170

160

150
0.9
0.8 6
0.7 5
4
0.6 3
Humidité relative 0.5 2
Vitesse de vent (m/s)

Figure 8: Evolution de la masse d’eau distillée en fonction de la vitesse du vent et de

l’humidité relative

La figure 7 montre l’évolution de production de masse d’eau distillée en fonction de

profondeur de l’eau de mer dans la cuve d’évaporation et la longueur du distillateur. On
remarque que dans l’intervalle d’étude, la productivité du distillateur augmente
considérablement avec la l’augmentation de la longueur de cuve et décroit lorsque l’épaisseur
de l’eau de mer augmente.

H Osman Ali
 Science et Environnement 34 (2018) 74-83 81

La figure 8 montre l’évolution de la production du système en fonction de vitesse de vent et

l’humidité relative. On remarque que la vitesse du vent influence considérablement la
production d’eau distillée. Par contre, on constante que l’humidité relative influence très peu
la production du distillateur solaire.
Le résultat obtenu permet de retrouver ce qui est déjà ressenti à partir de la physique du
problème, à savoir que l’augmentation de la vitesse du vent fait que la surface de
condensation devient plus frais par augmentation de la convection forcée, ce qui améliore la
condensation.

3 Conclusion

Nous avons étudié dans ce travail, un distillateur solaire actif en vue de conduire son
optimisation de manière paramétrique. Les équations physiques ayant servi à développer la
modélisation expriment les bilans thermiques pour les différents éléments composant
l’installation. Une étude paramétrique a été effectuée en fonction des paramètres de
conceptions. Nous pouvons synthétiser l’ensemble des résultats des études paramétriques
effectuées en concluant que lors de la conception d’un distillateur solaire actif, il est
préférable de prévoir son implantation dans une zone ventée, de prendre un bassin le plus
long possible et de veiller à limiter la profondeur de l’eau de mer dans ce bassin. Par ailleurs,
il convient mieux de travailler avec une température suffisamment élevée du fluide
caloporteur même si ceci exige de limiter au moyen d’un contrôle adéquat le débit du fluide
caloporteur.

Nomenclature:

Ac Aire de la surface de la plaque de cuivre ( m 2 )

A fc Aire de la surface de la cuve de fluide caloporteur ( m 2 )
Av Aire de la surface de l’enveloppe en verre ( m 2 )
A em Aire de la surface de la cuve de l’eau de mer ( m 2 )
c p,c Capacité thermique massique de cuivre ( [Link] 1.K 1 )
c p,em Capacité thermique massique de l’eau de mer ( [Link] 1.K 1 )
c p,fc Capacité thermique massique du fluide caloporteur ( [Link] 1.K 1 )
c p,s Capacité thermique massique de la saumure ( [Link] 1.K 1 )
c p,v Capacité thermique massique du verre ( [Link] 1.K 1 )
3 1
Dfc Débit volumique du fluide caloporteur ( m .s )

H Osman Ali
 Science et Environnement 34 (2018) 74-83 82

d em Profondeur de l’eau de mer ( m )

ev Epaisseur du verre ( m )
e is Epaisseur de l’isolant ( m )
2
I(t) Intensité solaire ( W.m )
hL Chaleur latente de vaporisation de l’eau ( [Link] 1 )
L Longueur du bassin ( m )
 em
m Débit massique de l’eau de mer entrant ( kg.s 1 )
 ev
m Débit massique de l’eau évaporée ( kg.s 1 )
s
m Débit massique de saumure sortant ( kg.s 1 )
Tamb Température ambiante ( K )
Ts,sor Température de saumure à la sortie ( K )
Tc Température de cuivre ( K )
Tfc Température du fluide caloporteur ( K )
Tfc,en Température d’entrée de fluide caloporteur ( K )
Tve Température du verre externe ( K )
Tvi Température du verre interne ( K )
Tem Température de l’eau de mer ( K )
Tem,en Température d’entrée de l’eau de mer ( K )
Tciel Température de ciel ( K )
t Temps ( s )
Vc 3
Volume du cuivre ( m )
3
Vfc Volume du fluide caloporteur ( m )
3
Vv Volume du verre ( m )

Vem Volume de l’eau de mer ( m3 )

 c Fraction de l’énergie solaire absorbée par une plaque de cuivre, sans dimension
v Fraction de l’énergie solaire absorbée par le verre, sans dimension
em Fraction de l’énergie solaire absorbée par l’eau de mer, sans dimension
 Masse d’eau distillée ( kg )
v -1 -1
Conductivité thermique du verre ( W.m .K )
c Masse volumique du cuivre ( kg.m 3 )
em Masse volumique de l’eau de mer ( kg.m 3 )
fc Masse volumique du fluide caloporteur ( kg.m 3 )
v Masse volumique du verre ( kg.m 3 )

H Osman Ali
 Science et Environnement 34 (2018) 74-83 83

Références

[1] L.G. Rodriguez, (2002), Seawater desalination driven by renewable energies: a review,
Desalination 143,103-113.
[2] E. Delyannis, (2003), Historic background of desalination and renewable energies, Solar
Energy, 75(5), 357-366.
[3] M.A.S. Malik, GN. Tiwari, A. Kumar, M.S. Sodha, (1982), Solar Distillation, Pergamon
Press Ltd, UK.
[4] H. Osman Ali, A. Khamlichi, (2014), Modelling an active solar still for sea water
desalination process optimization, Desalination 354, pp.1-8.

H Osman Ali