See discussions, stats, and author profiles for this publication at: [Link]

net/publication/336916118

Stockage d’énergie solaire thermique en utilisant les Nano-MCP

Conference Paper · October 2019

CITATIONS READS

0 320

3 authors:

Touati Boumediene Kerroumi Nordine

Université Tahri Mohammed Béchar Université Tahri Mohammed Béchar
34 PUBLICATIONS   78 CITATIONS    10 PUBLICATIONS   6 CITATIONS   

SEE PROFILE SEE PROFILE

Joseph Virgone
Claude Bernard University Lyon 1
208 PUBLICATIONS   2,932 CITATIONS   

SEE PROFILE

Some of the authors of this publication are also working on these related projects:

séchage solaire des plantes aromatiques et médicinales View project

Etude permettant d’optimiser au mieux le choix du triplet apports solaires/isolation/inertie en climat chaud View project

All content following this page was uploaded by Touati Boumediene on 17 November 2019.

The user has requested enhancement of the downloaded file.

 XIVème Colloque International Franco-Québécois en énergie
16-20 juin 2019, École de technologie supérieure, Université du Québec
Baie St-Paul, Québec, Canada

CIFQ2019 / ART-99

Stockage d’énergie solaire thermique en utilisant les Nano-

MCP

B. TOUATI1, N. kerroumi1, J. VORGONE2

1
laboratoire ENERGARID, équipe Gisement Solaire et ses Applications (GSA) Université de Tahri Mohammed, Bechar, Algérie

2
Centre de Thermique de Lyon (CETHIL, UMR 5008) Université Claude Bernard Lyon 1, Cedex, France

RÉSUMÉ

Cet article présente l’étude numérique de l’utilisation des nanoparticules intégrées dans les matériaux à changement de
phase pour stocker l'énergie thermique solaire. Les nanoparticules sont intégrées à une cire de paraffine pour augmenter
les performances thermiques de stockage du MCP. Le code ANSYS FLUENT 15 est utilisé pour la simulation
numérique basé sur la méthode des enthalpies pour résoudre les problèmes de changement de phase. L'effet des ailettes
est également étudié et trois User define qui ont été développées par le code C sont injectées dans le code ANSYS
fluent pour décrire les propriétés thermo-physiques des NMCP.

Mots Clés:NMCP, ANSYS FLUENT 15, ailettes, User define, stockage thermique, énergie solaire.

NOMENCLATURE

Symboles ρ: la densité (kg/m3)

H : l’enthalpie (J) μ : la viscosité dynamique (kg/m s)
T : température (K) λ : la conductivité thermique. (W/m2)
t : temps (S)  : concentration (%)
g : gravité (m s−2) Indices / Exposants :
C : la chaleur spécifique ([Link]−1.K−1) MCP: matériaux de changement de phase
L : chaleur latente NMCP nanomatériaux de changement de phase
Lettres grecques Sol : solide
 : la fraction liquide (%) Liqui : liquide, np : nanoparticules

Cette situation présente un déséquilibre entre l’offre

INTRODUCTION et la demande en énergie. Stocker l'énergie solaire peut
Le sud algérien est caractérisé par un climat aride, être une excellente solution pour combler cet écart. Vue
qui a une grande intensité solaire. Ce climat a une leurs capacités de stockage et leurs chaleurs latente
grande différence de température entre la période importantespar rapport aux matériaux sensibles tels que
diurne et nocturne. Pour les applications solaires, plutôt l’eau, les pierres, etc., Lesmatériaux de changement de
pose un réel problème, notamment dans le chauffage phase (MCPs)présentent une solution adéquate pour
des bâtiments et pour des applications d’eau chaude. stocker l’énergie solaire thermique.

-1- Copyright © 2019 CIFQ

 Les MCPS existent sous forme des MCPs utilisant.K. Biswas [6] présente une étude transitoire de
organique, les inorganiques, qui se trouvent avec une la fusion du Nano-MCP. Il at constaté que le temps de
large gamme dans le marché pour chaque température fusion total diminuait avec l’augmentation de la fraction
de fusion et capacité de stockage thermique. De plus volumétrique des nanoparticules à l’intérieur d’un
pour augmenter l’énergie solaire stockée, des plusieurs rectangle.
papiers de recherche présentent l’utilisation des ailettes L'effet des ailettes dans le système de stockage
pour augmenter la surface d’échange thermique avec thermique fait l'objet de nombreuses études. Jundika C.
des matériaux caractérisés par leur conductivité Kurnia et al. [7] évaluent numériquement diverses
thermique élevée. Finalement, une autre solution pour configurations des dispositifs de stockage d'énergie
augmenter l’énergie solaire emmagasinée durant la thermique par MCP : un tube en U, un tube en U avec
journée est l’utilisation des nanoparticules en vue du des ailettes en ligne, un tube en U avec des ailettesen
mélange avec le MCP pour améliorer quinconce .Les résultats indiquent que cette dernière
l’échangethermique (la conduction et la convection au accélère l’échange thermique entre le fluide caloporteur
niveau de l’MCP). Dans ce papier on présente une et le MCP en comparent avec le tube en U avec des
étude numérique de l’utilisation des NMCP dans une ailettes en lignes. De même Abduljalil A. Al-Abidi et al.
cuve de stockage pour but stocker de l’énergie solaire [8] présentent une étude numérique en utilisent le code
avec l’ajout des ailettes et une concentration de 8% des fluent 6.3.2 d’un tube triplex avec et sans ailettes
nanoparticules. internes et externes. Ils ont constaté que l’utilisation des
ailettes améliore le transfert de chaleur entre le fluide
MCP ET NMCP POUR STOCKER DE caloporteur et le MCP. Dans le même sens B. Touati et
L'ÉNERGIE al. [9] prouvent l’accélération du transfert de chaleur
Comme indiqué dans l’introduction, le stockage de dans le MCP avec la présence des ailettes ceci conduit à
l’énergie solaire devient une solution pour réduire réduire le temps total de décharge d’énergie à partir
l’écartentre la demande et l’offre d’énergie pendant la d’un MCP. Aussi Maryam Gharebaghiet I. Sezai [10]
[Link] nombreuxchercheurs ont étudié le stockage de ont arrivé a amélioré le taux de stockage d'énergie dans
l’énergie à l’aide du MCP. Soroush Ebadi et al. [1] ont une unité de stockage d'énergie thermique remplie par
étudié numériquement et expérimentalement la fusion un MCP, où ils en insérant un réseau d'ailettes dans le
de Nano-PCM dans un cylindre pour stocker l’énergie dispositif de stockage.
thermique. Ils ont dispersé les nanoparticules d’oxyde
de cuivre dans le MCP et ont examiné les effets des PRÉSENTATION DE LA GEOMETRIE
fractions volumiques des nanoparticules. Ils ont constaté
que les fractions en volume croissantes des Dans cette section, le schéma de l'échangeur de
nanoparticules augmentent l’énergie stockée. chaleur a tube triplex en présence d’ailettes est
Jasim M. Mahdi et Emmanuel C. Nsofor [2] ont présenté.Il s’agit d’un ballon de stockage remplit d’un
étudié la solidification du MCP avec des nanoparticules matériau de changement de phase avec des
d’alumine plutôt l'oxyde d'aluminium (AL2O3) dans un nanoparticules. Cette géométrie peut être utilisée dans
tube triplex pour un système de stockage d’énergie différents applications comme : un ballon de stockage
thermique. Leurs résultats montrent que les pour un chauffe-eau solaire, un échangeur de chaleur
concentrations volumétriques des nanoparticules pour un système de de chauffage du bâtiment, aussi un
comprises entre 3 et 8% permettent d’économiser entre échangeur de chaleur pour un système de dessiccation
8 et 20% du temps total de solidification. La figure 1 représente le domaine d’étude concerné.
Aussi Mohammad Parsazadeh et Xili Duan. [3] Les rayons intérieur et extérieur
présentent un papier dans lequel ils simulent une coque mesurentrespectivement 75 mm et 20mm. La longueur
et un tube contenant un Nano-MCP. Ils utilisent un nano des ailettes est de 50 pour une largeur de 1mm. Enfin
fluide CuO - eau dans un système de stockage d'énergie l’espace entre l’intérieur et l’extérieur contient la cire de
thermique. Leurs résultats présentent l’effet de la paraffine en tant que MCP avec 8% d'Al2O3 sous forme
concentration volumique desnanoparticules. De plus de nanoparticules
plusieurs cavités contenant du nano -MCP sont étudiés Le fluide caloporteur circulant dans le tube interne
analytiquement et numériquement par Nabeel S. est l’eau venant d’un chauffe-eau [Link] utilisant
Dhaidan et al. [4] .Ils utilisent différentes géométries, et cette géométrie simplifiée, cette étude se limite à la
comparent les résultats entre les deux études. moitié du domaine en raison de la symétrie.
Une étude expérimentale et numérique sur
l'utilisation de Nano-MCP dans un récipient carré est
présentée par Nabeel S. Dhaidan et al. [5]. Ils utilisent le
n-octadécane en tant que MCP avec une nanoparticule
de CuO avec un flux constant. Ainsi,la méthode des
éléments finis est utilisée pour validationnumérique

-2- Copyright © 2019 CIFQ

 T
hnmcp = h +  cp
ref
nmcp nmcp dT (5)
Tref

ref
Où hnmcp est l'enthalpie de référence, cpnmcp
chaleur spécifique à pression constante.

Enfin la fraction liquide est définie comme:

H npcm =  L (6)

Où L est la chaleur latente de fusion, β est la fraction

liquide.
La fraction liquide intervient au moment de changement
de phase entre l'état solide et l'état liquide lorsque la
température est Tliqui> T> Tsol, de sorte qu'il peut être
écrit comme suit:

Figure 1 : domaine étudié 


ÉQUATIONS GOUVERNANTES 0 pourT Tsol

 = 1 pourT Tliqui (7)
ANSYS FLUENT 15 utilise la méthode de 
l'enthalpie pour résoudre les équations qui régissent le  T − Tsol pourTsol T Tliqui
problème de la solidification – fusion [11]. Les  Tliqui − Tsol
équations suivantes sont :
Le NMCP considéré dans cette étude est la cire de
-L’équations de conservation de la masse paraffine où la densité, ρnpcm, est donnée par:

 nmcp l
+ (  nmcp u ) = 0 (1)  nmcp = (8)
t 0.001(T − Tl ) + 1
-l'équation de conservation de la quantité de Où la constant 0.001 représente le coefficient
mouvement: d’expansion thermique.
nmcp
+ ( nmcp u  u ) =  + ( nmcp g ) (2) La conductivité thermique du MCP est définit par
t
- l’équation de conservationde l’énergie: npcm
(s)
siT p Tsol
 (l ) (9)
 npcm siT f Tliqui
(  nmcp H nmcp ) + .(  nmcp uH nmcp ) = .(T ) (3)
t
La viscosité dynamique est donnée par la formule :
Où ρnpcm, μnpcm et Hnmcp sont la densité, la viscosité
dynamique et l'enthalpie du NMCP, g est la gravité et λ npcm = 0.001exp(−4.25 + 1790 / T )
la conductivité thermique.
La densité, la capacité thermique spécifique et
-L'enthalpie des NMCP Hnpcm peut être exprimée latente et la chaleur du Nano-PCM sont données par les
par: équations (10) à (12)

H nmcp = hnmcp + H nmcp (4)

 nmcp =  np + (1 −  )  mcp (10)

-3- Copyright © 2019 CIFQ

  (  c p ) np + (1 −  )(  c p ) mcp
cpnmcp = (11)
 npcm
MAILLIAGE ET SIMULATION
(1 −  )(  L) mcp
Lnpcm = (12) Le maillage du domaine d’étude est généré par un
 npcm CFD (computationalfluiddynamics).Pour le NPCM,
44940 éléments utilisée (nœuds), et pour les ailettes,
Où  est la concentration volumétrique de 1080 éléments.
nanoparticules. La température initiale du système est de 319 K et
[Link] et al.[12]mettent en équation la celle de l’eau est de l'eau est 328 K. La paroi du tube
viscosité dynamique ( équation (13)) et la conductivité extérieur est considérée comme parfaitement isolée.
thermique effective (équation (14)) Pour la simulation sous ANSYSFLUENT 15,
l’algorithme SIMPLE est utilisé pour le couplage
pression-vitesse. Le schéma d’optimisation de la
nmcp = 0.983e(12.958 ) pression (PRESTO) est utilisé pour l’équation de
(13) correction de pression, les valeurs des facteurs de sous-
relaxation sont 0.3, 0.2, 1 et 0.9 pour la pression, la
np + 2mcp − 2(mcp − np ) vitesse, l’énergie et fraction liquide [Link]
npcm =  pas de temps est d’ordre 0,1 s (choisi).
np + 2mcp + 2(mcp − np ) mcp Les trois USER DEFINE sont écrits sur le code C
(UDF) pour décrire les propriétés thermo-physiques des
BT NMCP.
+5*104  k  mcp C p
 npd np f (T , ) Afin de valider la géométrie et le maillage choisi,
les travaux de Abduljalil A. Al-Abidi et al. [8] ont été
(14) Où B est la constante de Boltzmann égal à 1,381 × repris.
10-23 J / K La figure 2présente la comparaison entre les
résultats obtenus dans cette étude et les résultats donnés
La deuxième partie de l’équation (14) explique le par Abduljalil A. Al-Abidi et al. [8]. Les courbes étant
mouvement brownien, qui entraîne la dépendance en proches, la géométrie et le maillage choisi sont justifiés
température de la conductivité thermique du NMCP.

f (T ,  ) = (2.8217*10( −2)  + 3.917*10( −3)T /Tref )

(15)
+(−0.669*10( −2)  − 3.91123*10( −3) )

Le tableau 1 présente les propriétés thermo-

physique du MCP,du cuivre et de l’eau

Tableau 1 : propriétés thermo-physiques du mcp,cuivre

et de l'eau.
Propriété Cire de Paraffin Al2O3 Eau

Densité ρ (kg/m3) l=736 s=751 3600 998.2

Capacité spécifique (J/kg K) 2890 765 4182

Chaleur latente de fusion, L 173400 - - Figure 2 : Comparaison entre la fraction liquide

(J/kg)
présentée dans [8] et celle déterminée dans cette étude
Température de fusion, (K) 321 - -

Température de solidification 319 RÉSULTATS

Condictivité thermique, (W/m 0.21 si T<Tsol , 36 0.6
K) 0.12 si T >Tliq La figure 3représente la température des
Viscosité dynamique, μ(Kg/m S=0.294 , l=0.286 - 0.1003 NMCP après 3 et 6 heures de simulation, en mode
s) [Link] transfert de chaleur par conduction et
convection naturelle existe et l'effet des ailettes est
observé.

-4- Copyright © 2019 CIFQ

 La figure 5 montre les lignes de courant dans le NMCP
après 4 heures de simulation.

(a)

Figure 5 : Lignes de courant dans le NMCP après 4

heures
(b)
La figure 6 présente la fraction liquide dans le
Figure 3 : Température des NMCP après (a) 3 heures NMCP en fonction du temps, aussiappelé processus de
(b) 6 heures chargement desNMCP. L’augmentation du liquide dans
le NMCP suit la courbe d’apprentissage suivante : après
La figure 4 montre la fraction liquide dans le 2 heures,45% de liquide est présent dans le tube , 60%
tube après 3 et 6 heures de simulation. Après 3 heures, après 3 heures et après 6 heures 88% .
60% de liquide est présent dans le NMCP, contre 88 %
après 6 heures. Les résultats obtenus indiquent que le
fond du tube contient plus de liquide .Ceci est expliqué
par l'effet de la convection naturelle et la différence de
densité.

(a)
Figure 4: Fraction liquide dans le NMCP après (a)3
heure(b)6heures

Figure 6 : Fraction liquide dans le NMCP en fonction

du temps

(b) CONCLUSION

La figure 4 : Fraction liquide dans le tube après (a) Dans cette étude, la simulation numérique d’un
3heures et (b) 6 heure. système de stockage d’énergie solaire thermique par un
tube contenant des NMCP est réalisée sur ANSYS
FLUENT 15. Des ailettes et une concentration en
particules de 8%. Les ailettes jointes à l’unité
augmentent la surface d’échange d’énergie entre l’eau et
les NPCM et par la suite réduisent le temps de fusion
total du NPCM. De plus, l’utilisation d’une
concentration volumique des nanoparticules

-5- Copyright © 2019 CIFQ

 augmentel’échanges de chaleur par convection et [8] A A Al-Abidi, Sohif Mat, K Sopian, MY
conductiondans les NPCM. Sulaiman and AbdulrahmanTh Mohammad, “
La convection a un effet remarquable dans la Numerical study of PCM solidification in a triplex
géométrie car il y a plus de fraction liquide dans le fond tube heat exchanger with internal and external fins
du tubequ’en haut. L’utilisation des NPCM présente une ”, International Journal of Heat and Mass
solution pour minimiser l’écart dans la demande Transfer, vol. 61, pp. 684–695, 2013.
d’énergie entre la période diurne et nocturne. [9] B. Touati, N. Kerroumi and J. Virgone,” Solar
Cette étude pourra être compléter prochainement thermal energy discharging from a multiple phase
par l’étude de l’effet du nombre et la taille des ailettes change materials storage tank,” Applied Solar
L’effet de la concentration des nanoparticules pourra Energy, Vol.53, pp 185–189, 2017.
être déterminé à travers d’autre matériaux, comme le [10] Maryam Gharebaghiand I Sezai, “Enhancement of
cuivre et l’or. Heat Transfer in Latent Heat Storage Modules
with Internal Fins”, Numerical Heat Transfer,
vol.A(53), pp. 749–765, 2008.
[11] M. Auriemma and A. Iazzetta, “Numerical
RÉFÉRENCES Analysis of Melting of Paraffin Wax with Al2O3,
ZnO and CuO Nanoparticles in Rectangular
Enclosure,” Indian Journal of Science and
[1] SoroushEbadi, Syeda Humaira Tasnim, Amir A. Technology, Vol. 9, Issue 3, January 2016.
Aliabadi and Shohel Mahmud, “Melting of nano- [12] Jasim M. Mahdi, Sina Lohrasbi, Davood D .Ganji
PCM inside a cylindrical thermal energy storage and Emmanuel C. Nsofor ,”Accelerated melting of
system: Numerical study with experimental PCM in energy storage systems via novel
verification ,”Energy Conversion and Management configuration of fins in the triplex-tube heat
, vol. 166, pp.241–259, 2018. exchanger,” International Journal of Heat and
[2] Jasim M. Mahdi and Emmanuel C. Nsofor, Mass Transfer ,Vol. 124, pp. 663-676, 2018
“Solidification of a PCM with nanoparticles in
triplex-tube thermal energy storage system ,”
Applied Thermal Engineering , vol. 108 ,pp. 596–
604, 2016
[3] Mohammad Parsazadehand Xili Duan, “Numerical
and statistical study on melting of nanoparticle
enhanced phase change material in a shell-and-
tube thermal energy storage system,”Applied
Thermal Engineering, vol. 111, pp. 950–960,
2017.
[4] Nabeel S. Dhaidan, “Nanostructures assisted
melting of phase change materials in various
cavities,” Applied Thermal Engineering, vol. 111,
pp.193–212, 2017.
[5] Nabeel S. Dhaidan, J.M. Khodadadi, Tahseen A.
Al-Hattab and Saad M. Al-Mashat, “Experimental
and numerical investigation of melting of phase
change material/nanoparticle suspensions in a
square container subjected to a constant heat flux,”
International Journal of Heat and Mass Transfer
,vol. 66, ppt. 672–683, 2013.
[6] [Link], “Nano-based phase change materials
for building energy efficiency,” Start-Up Creation
the Smart Eco-Efficient Built Environment, pp-
183-211, 2016.
[7] Jundika C. Kurnia, Agus P. Sasmito, Sachin V.
Jangam and Arun S. Mujumdar, “Improved design
for heat transfer performance of a novel phase
change material (PCM) thermal energy storage
(TES) ” Applied Thermal Engineering, vol.50, pp.
896-907, 2013.

-6- Copyright © 2019 CIFQ

View publication stats