Pompe à chaleur résidentielle avec stockage
Pompe à chaleur résidentielle avec stockage
Doctorat ParisTech
THÈSE
pour obtenir le grade de docteur délivré par
Samer MAARAOUI
le 13 mars 2013
Jury
M. Dominique MARCHIO, Professeur, MINES ParisTech Président
T
M. Jocelyn BONJOUR, Professeur, INSA Lyon Rapporteur
M. Joseph VIRGONE, Professeur, INSA Lyon Rapporteur
M. Denis CLODIC, Directeur de recherche émérite, MINESParisTech Directeur de thèse
H
È
S
E
MINESParisTech
Centre d’Efficacité énergétique des Systèmes (CES)
5, Rue Léon Blum, 91120 Palaiseau
Thoughts Become Things…
J’adresse également mes remerciements à tous les membres et invité du jury, M. Dominique
Marchio, M. Joseph Virgone, M. Jocelyn Bonjour et M. Pascal Dalicieux de m’avoir fait
l’honneur d’être dans mon jury de thèse.
Je remercie Anne-Marie Bonnet, d’être toujours là pour les aides aimables et pour les
relectures minutieuses des rapports.
Je remercie Elias Bou Lawz, pour son guidage technique et pour le rapport d’amitié que j’ai
pu entretenir avec lui.
Je remercie Maroun Nemer, responsable du CEP, qui m’a permis d’accomplir ce travail.
Je remercie Gérard Helie pour son accueil chaleureux souriant chaque matin.
Je remercie Youssef Larabi et Rachid Dib pour leur assistance dans le montage du banc
d’essais et du prototype.
Je remercie également tous les membres du CEP pour les rapports amicaux que j’ai pu
entretenir avec eux et pour tous les bons moments que nous avons passés ensemble.
Un grand merci est adressé à ma sœur Rouba, son époux Micho, ma cousine Grace, mon
cousin Georges et sa femme Randa qui m’ont soutenu dès mon arrivée en France.
Un merci du fond de mon cœur à ma chère Liliane qui a vécue avec moi les moments de
stress et les moments de joie de cette période.
i
Sommaire
V.8 Résumé propriétés des familles MCP ...................................................................................... 44
V.9 MCP commercialisés ....................................................................................................................... 45
VI. Modélisation du phénomène de changement de phase ............................................................ 47
VI.1 Etat de l’art ......................................................................................................................................... 47
VI.2 Thermodynamique du changement de phase...................................................................... 49
VI.3 Cinétique du changement de phase.......................................................................................... 49
VI.4 Modélisation par position du front........................................................................................... 50
VI.5 Modélisation par méthode numérique (volumes finis) ................................................... 51
VII. Calorimétrie différentielle à balayage ....................................................................................... 53
VII.1 Calorimétrie des MCP ..................................................................................................................... 53
VII.2 Calorimétrie mesurée des candidats ....................................................................................... 54
VII.3 Calcul des calorimétries vraies................................................................................................... 56
VII.4 Comparaison des calorimétries ................................................................................................. 57
VIII. Conclusions........................................................................................................................................... 57
Chapitre III. Conception d’un Echangeur Stockeur - Déstockeur .................59
I. Introduction ................................................................................................................................................ 59
II. Etat de l’art sur les échangeurs stockeurs ...................................................................................... 59
III. Modèle échangeur complet ................................................................................................................... 61
III.1 Langage Modélica ............................................................................................................................ 61
III.2 Modèles élémentaires principaux ............................................................................................. 61
IV. Structure de l’échangeur stockeur ..................................................................................................... 64
IV.1 Géométrie « tube dans tube » ..................................................................................................... 64
IV.2 Optimisation de la structure tube dans tube ........................................................................ 76
IV.3 Echangeur de première génération.......................................................................................... 76
IV.4 Tubes entre plaques ....................................................................................................................... 80
IV.5 Structure d’un prototype optimisé ........................................................................................... 86
IV.6 Structure du prototype .................................................................................................................. 88
V. Conclusions .................................................................................................................................................. 92
Chapitre IV. Validation Expérimentale de la PAC à Stockage ........................95
I. Introduction ................................................................................................................................................ 95
II. Schémas du banc ....................................................................................................................................... 95
III. Acquisition et commande ...................................................................................................................... 96
IV. Sélection des éléments principaux..................................................................................................... 99
IV.1 Le compresseur ................................................................................................................................ 99
IV.2 Les détendeurs .................................................................................................................................. 99
V. Première génération de l’échangeur stockeur ........................................................................... 100
V.1 Construction de la première génération ............................................................................. 100
ii
Sommaire
V.2 La décharge thermique ............................................................................................................... 103
VI. Le deuxième prototype ........................................................................................................................ 106
VI.1 Construction du deuxième prototype .................................................................................. 106
VI.2 Validation expérimentale du prototype .............................................................................. 109
VII. Conclusions........................................................................................................................................ 113
Chapitre V. Pompe A Chaleur avec Stockage ...............................................115
I. Introduction ............................................................................................................................................. 115
II. Pompe à chaleur dans un bâtiment ................................................................................................ 115
II.1 Régime cyclé des pompes à chaleur ...................................................................................... 116
II.2 Compresseur à vitesse variable .............................................................................................. 117
II.3 Givrage/dégivrage ........................................................................................................................ 119
II.4 Variation des COP de PAC en fonction de la température extérieure ..................... 120
III. PAC à condenseur stockeur ............................................................................................................... 121
III.1 Effacement lors du pic de demande d’électricité ............................................................. 121
III.2 Système de stockage gagnant .................................................................................................. 122
III.3 Avantages de la PAC avec stockage ....................................................................................... 125
IV. Conclusions ............................................................................................................................................... 131
Conclusions et Perspectives ................................................................................133
Références bibliographiques .................................................................................................................... 137
iii
Nomenclature
Nomenclature
1
Nomenclature
Lettres grecques
titre liquide
µ viscosité Pa.s-1
Δ différence
λ conductivité thermique W.m-1.K-1
ρ masse volumique kg.m-3
φ nombre de phases
Indices
ail ailette
bat bâtiment
comp compresseur
cond condenseur
cont contact
crit critique
dil dilution
diss dissolution
E eutectique
ech échangeur
ef énergie finale
ent entrée
ep énergie primaire
eq équivalent
ev évaporateur
ext externe
f fusion
f frottement
h hydraulique
i changement de phase
int interne
L ou l liquide
l0 tout liquide
max maximal
min minimal
P péritectique
p paroi
r réaction
r référence
S/s solide
sort sortie
th thermique
therm thermodynamique
tot total
v vapeur
v0 tout vapeur
2
Nomenclature
Abréviations
BBC bâtiment à basse consommation énergétique
Bbio besoins bioclimatiques
CEP centre énergétique et procédés
CL chaleur latente
COP coefficient de performance
CS chaleur sensible
DSC differential scanning calorimetry
Ep epaisseur
M Microtek
MCP matériau à changement de phase
PAC pompe à chaleur
RT règlementation thermique
RT Rubitherm
SCS stockage par chaleur sensible
SHON surface hors œuvre nette
Nombres adimensionnels
Bo nombre d’ébullition
Co nombre de convection
Fr nombre de Froude
j facteur de Colburn
Pr nombre de Prandtl
Nu nombre de Nusselt
Re nombre de Reynolds
We nombre de Webber
3
4
Introduction générale
Introduction générale
La production et la demande en énergie électrique sont soumises à des aléas (indisponibilité
des moyens de production, climat, énergies renouvelables intermittentes) et l’électricité ne se
stockant pas facilement, il est difficile d’assurer en toutes circonstances l’adéquation de la
production et de la demande, surtout si on veut limiter le recours aux centrales à flamme.
Aujourd’hui, en France, le principal aléa pesant sur la consommation d’électricité est la
température extérieure avec une sensibilité atteignant 2 100 MW.°C-1 en hiver contre moins
de 1 500 MW.°C-1 il y a quelques années. Le principal responsable de cette thermo-
sensibilité étant le chauffage par effet Joule, il est urgent de lui trouver des solutions de
substitution et mieux encore des solutions pilotables, en particulier par des Pompes à
Chaleur (PAC) qui doivent diviser au moins par 3 la consommation d’électricité associée au
chauffage du fait de leur coefficient de performance saisonnier élevé. De plus, le constat
d’une certaine érosion du potentiel d’effacement consenti volontairement par les
consommateurs amène à chercher de nouveaux gisements d’effacement de consommation,
dans le cadre d’une gestion intelligente des systèmes.
Un premier objectif de la thèse est de doter une pompe à chaleur air/air destinée à une
maison individuelle basse consommation d’un dispositif « condenseur stockeur/déstockeur »
permettant l’effacement total ou partiel de la PAC pendant au moins deux heures. La
conception de cette pompe à chaleur s’effectue dans le cadre des maisons à Basse
Consommation (BBC), voire à très basse consommation (BBC+).
La toute première étape consiste à établir les besoins de chauffage horaires d’une maison
individuelle BBC ou BBC+ pendant la saison de chauffe. Ces calculs sont présentés dans le
chapitre 1 en abordant plusieurs scénarios et variantes qui définissent la puissance de
chauffage d’un tel logement. La quantité de chaleur à stocker a été évaluée afin d’assurer
l’effacement de la PAC entre 18h00 et 20h00 durant toute la saison de chauffe. La puissance
thermique maximale à fournir lors de l’effacement a été également évaluée pour les maisons
BBC, mais aussi en réhabilitation du parc existant.
Le chapitre 2 traite les phénomènes de changement de phase pour les différents types de
matériaux. Les critères d’un MCP (Matériau à Changement de Phase) idéal ont été identifiés
permettant de définir un cahier des charges pour choisir un tel matériau. Les MCP ont été
classifiés selon leurs natures, avec une description pour chaque famille et une liste la plus
exhaustive possible des matériaux de chaque famille. Un modèle dynamique du phénomène
de changement de phase a été développé dans ce chapitre. Il est basé sur la forme
enthalpique de l’équation de l’énergie. Une DSC (Differential Scaning Calorimetry) a été
réalisée sur des MCP commercialisés, candidats pour l’utilisation dans la PAC avec
stockage. Cette DSC a été corrigée en utilisant une méthode inverse et le modèle développé
de changement de phase.
Une fois le comportement des MCP analysé, le chapitre 3 présente la conception de
structures possibles du condenseur stockeur/déstockeur. La modélisation ayant permis
d’identifier les conditions adaptées pour les transferts thermiques avec le MCP, l’air et le
fluide frigorigène, les structures ont été optimisées afin d’obtenir une taille raisonnable. Il est
à noter que la conductivité thermique des MCP généralement faible (de l’ordre de
0,2 W.m -2.k-1 pour les paraffines), amène à concevoir des agencements innovants entre le
5
Introduction générale
MCP et ses surfaces d’échange à densité de flux très différentes. Dans ce chapitre deux
échangeurs stockeurs optimisés ont été dimensionnés.
Le chapitre 4 décrit le banc d’essais qui a été construit afin de tester la PAC comportant
l’échangeur stockeur. Une première structure optimisée a été réalisée et testée sur le banc.
Cet échangeur de première génération a donné des résultats intéressants. En se basant sur
l’expérience acquise de la première génération, la deuxième génération a été réalisée. Cet
échangeur a été intégré dans une PAC prototype qui a été également testé sur le banc
d’essais et a donné des résultats favorables qui respectent le cahier des charges d’une PAC
efficace.
Une fois la PAC conçue, son comportement est étudié expérimentalement dans le chapitre
4. L’effet du cycle stockage/déstockage a été évalué en termes d’efficacité et d’émissions
CO2. Dans le chapitre 5, un modèle simulant le fonctionnement de trois PAC air/eau avec
compresseur à vitesse variable, air/air avec compresseur à vitesse variable et air/air avec
condenseur stockeur a été développé. Ce modèle, basé sur des résultats d’essais, permet
de prendre en compte les irréversibilités réelles associées au cycle de givrage et dégivrage
et aux fonctionnements des PAC en cycles courts associés aux faibles charges. Ce modèle
a permis ainsi de comparer la PAC à stockage telle que conçue à une PAC air/eau et une
PAC air/air à vitesse variable. Le modèle prévoit une amélioration entre 20 % et 30 % du
COP saisonnier par rapport aux deux autres technologies et ce du fait de l’intégration du
stockage et déstockage dans la régulation globale de cette PAC.
6
Chapitre I. Le stockage d’énergie intégré dans une PAC
I. Introduction
L’offre et la demande en énergie électrique étant soumise à des aléas (indisponibilité des
moyens de production, climat) et l’électricité ne se stockant pas facilement, il est difficile
d’assurer en toutes circonstances leur adéquation. Aujourd’hui, en France, le principal aléa
pesant sur la consommation d’électricité est la température extérieure avec une sensibilité
atteignant 2 300 MW.°C-1 en hiver contre moins de 1 500 MW.°C-1 il y a quelques années. Le
principal responsable de cette thermo-sensibilité étant le chauffage par effet Joule, il est
urgent de lui trouver des solutions de substitution et mieux encore des solutions pilotables.
En effet, le constat d’une certaine érosion du potentiel d’effacement consenti volontairement
par les consommateurs amène à chercher de nouveaux gisements d’effacement de
consommation, dans le cadre d’une gestion intelligente des systèmes.
Dans ce contexte, la Direction Optimisation Amont-Aval Trading d’EDF a défini plusieurs
formats d’effacement susceptibles de rendre un service intéressant du point de vue de
l’optimiseur (capter l’essentiel des pics de consommation et couvrir l’aléa climatique) :
soir 2 h de 18h00 à 20h00
matin 8 h de 6h00 à 14h00
jour 16 h de 6h00 à 22h00
permanent pour réduire le différentiel été-hiver
Cette thèse vise à mettre au point une solution couplant une pompe à chaleur air/air et un
stockage thermique par des matériaux à changement de phase (solide /liquide et vice versa)
pour répondre à la demande d’effacement de la pointe entre 18 h et 20 h (Figure I-1). Le fait
de savoir si cette solution technologique peut être utilisée aussi bien en rénovation, par
remplacement du chauffage effet Joule du parc existant, que pour les bâtiments à basse
consommation de demain dépend, bien sûr, des questions de coûts des équipements de
chauffage qui sont en concurrence sur ces marchés. Concevoir des pompes à chaleur
7
Chapitre I. Le stockage d’énergie intégré dans une PAC
effaçables en y associant des stockages thermiques locaux introduit un nouveau levier
d’optimisation aval de l’équilibre offre et demande d’électricité. Bien à l’opposé du stockage
saisonnier, la gestion délocalisée de la pointe journalière, dont la durée est de l’ordre de une
à deux heures, présente des enjeux multiples, notamment en termes d’impact
environnemental, de dimensionnement et de sécurité des réseaux.
Il est important de noter aussi le fait que cette thèse cherche, en plus de l’effacement des
pics d’appel de puissance, à améliorer l’efficacité des systèmes de Pompe A Chaleur (PAC).
Ceci est en effet possible, et compatible avec l’effacement, en faisant fonctionner la PAC
pendant les périodes les plus favorables de la journée, c’est-à-dire la fin de l’après-midi
(Figure I-2), juste avant le pic d’appel. De plus, la capacité de stockage de l’énergie
thermique devrait améliorer le COefficient de Performance (COP) saisonnier en faisant
fonctionner la PAC de manière plus continue pour les faibles charges thermiques et pendant
les heures où la température d’air extérieur est la plus élevée. C’est l’objet du chapitre 4
d’analyser ce que le stockage thermique peut amener pour améliorer l’efficacité énergétique
du système PAC, alors qu’a priori l’introduction d’un Matériaux à Changement de Phase
(MCP) entre l’air et le fluide frigorigène introduit non seulement, un écart de température
supplémentaire entre le fluide frigorigène et l’air, mais aussi un seuil en température qui n’est
pas forcément optimum pour l’obtention d’un COP saisonnier au moins égal à celui de la
même PAC sur l’ensemble de la saison de chauffe.
8
Chapitre I. Le stockage d’énergie intégré dans une PAC
l’ensemble des usages électriques dits spécifiques (cuisson, réfrigération, lavage, séchage,
télévision, communication, informatique). Pour toutes ces raisons, l’isolation des parois,
l’étanchéité à l’air, les vitrages font l’objet de développements tels que la consommation de
chauffage dans la région parisienne peut être réduite à moins de 1 kWhep.m-2.an-1. Ce niveau
change totalement la conception des pompes à chaleur et amène à les intégrer dans le
système aéraulique.
De là, à la fois pour des raisons de coûts (de matériels et d'installation) mais aussi d’inertie
thermique, les systèmes centralisés avec des radiateurs sur boucle d’eau chaude font de
moins en moins sens. La rénovation du parc existant aura à trouver des solutions
compatibles avec la vectorisation existante de la chaleur. Pour ce travail de thèse, nous
prendrons comme référence un système de chauffage par pompe à chaleur air /air et ce
pour une maison individuelle. La question qu’il faudra évaluer par modélisation et essais sera
de savoir si cette PAC sera centralisée avec séparation des flux d’air dans chaque pièce ou
sera un système multi-split avec des condenseurs dans chaque pièce.
L’arrivée programmée d’un nouveau dispositif de comptage qui permet notamment de
programmer, relever, et vérifier le fonctionnement du compteur à distance, offre de nouvelles
perspectives en termes de télé-services et de pilotage des usages. Effacer, au moins
partiellement, certains équipements constitue un moyen aval de régulation de l’équilibre
offre-demande qui allègerait les contraintes du système électrique. Il s’agit donc d’étudier les
possibilités théoriques de foisonnement et de gestion en décalé de ces effacements courts et
diffus (sur un très grand nombre de sites) pour obtenir, par agrégation, un effet de lissage de
la consommation d’électricité, et donc de limiter les pointes supportées par le réseau,
notamment en période de grand froid.
Sous réserve de disposer d’un système de stockage thermique compact, le chauffage est un
usage potentiellement « dispatchable ». L’allure du signal France d’une journée d’hiver
(Figure I-1) montre tout l’intérêt d’une gestion délocalisée de la pointe.
Le Tableau I-2 [RTE] permet également de confirmer les consommations exceptionnelles
enregistrées pendant l’hiver 2009-2010, qui se situent effectivement toutes pendant la pointe
du soir, qui est bien la cible de cette étude.
Tableau I-2 – Puissances instantanées relevées en janvier 2009 et 2010
Date Heure Puissance instantanée
Mercredi 7 janvier 2009 19h00 92 400 MW
Mardi 12 janvier 2010 18h58 91 537 MW
Mardi 6 janvier 2009 19h01 91 500 MW
Jeudi 8 janvier 2009 19h03 91 402 MW
Vendredi 8 janvier 2010 18h54 91 383 MW
9
Chapitre I. Le stockage d’énergie intégré dans une PAC
La problématique consiste à stocker de la chaleur lors du creux observé en fin d’après-midi,
i.e. au moment où la température extérieure, et donc le coefficient de performance des PAC
aérothermiques, sont généralement les plus élevés, pour la restituer lors des pics de
consommation. Une durée d’effacement de quelques heures est déjà pertinente pour que le
système électrique en tire un bénéfice significatif.
Sur l’exemple de la Figure I-2 montrant l’évolution de la température extérieure une journée
d’hiver en région parisienne, un écart de 7 à 8 degrés est observé entre les heures les plus
chaudes et la tranche de pointe (autour de 19h00). Le COP des pompes à chaleur sur l’air
variant d’environ 0,1 point par degré de température d’air, un gain potentiel est estimé sur la
performance de la machine, pouvant atteindre 20 à 25 % en privilégiant son fonctionnement
durant ces heures les plus chaudes. De plus, lorsque les charges thermiques sont faibles
(typiquement pour des températures extérieures variant entre 9 et 14 °C), le stockage
thermique permet de faire fonctionner la PAC de manière continue en évitant ainsi les courts
cycles qui constituent une pénalité énergétique majeure. Cet effacement des courts cycles
est évalué au chapitre 5.
Aujourd’hui, en prenant pour exemple l’application Pompe à Chaleur (PAC) air/air, en
utilisant un matériau à changement de phase optimisé pour une température de transition
correspondant à l’application, soit de 30 à 45 °C, un effacement de 4 kWh (2 heures de
chauffage pour une pièce de séjour pas très bien isolée) nécessite environs 60 litres de
matériau (en prenant pour capacité de stockage 60 kWh/m3).
Figure I-3 – Le réservoir MCP conçu par F. Agyenim et al. Figure I-4 – Le réservoir MCP conçu par
[AGY10] U. Stritih [STR04]
10
Chapitre I. Le stockage d’énergie intégré dans une PAC
U. Stirih a conçu deux réservoirs contenant de la paraffine, pour stocker de la chaleur latente
utile dans le domaine des bâtiments (30 °C). Le premier système [STR03] consiste à
incorporer de la paraffine dans les parois externes du bâtiment. Durant la journée, la
paraffine, fondant à 30 °C, absorbe la chaleur solaire transmise par radiation. La restitution
de cette chaleur se déroule durant la soirée par convection et par conduction. Le second
réservoir présenté par U. Stritih [STR04] (Figure I-4) ayant une géométrie rectangulaire,
échange de la chaleur avec de l’eau circulant sur une de ses facettes. Il existe dans la
littérature un grand nombre d’études sur les designs et les optimisations des réservoirs à
MCP comme [LIU06][MAR10][SHU09][BRO96][CON96][HAM06].
Figure I-5 – Echangeur air - MCP conçu par P. Figure I-6 – Echangeur air - MCP conçu par M.
Dolado et al. [DOL11] Zukowski [ZUK07a]
L’intégration des MCP dans les gaines d’air est aussi utilisée pour le stockage par chaleur
latente. P. Dolado et al. [DOL11] (Figure I-5) par exemple ont modélisé et expérimenté un
échangeur air - MCP. Le MCP est « macro-encapsulé » entre des plaques en aluminium et
l’air circule parallèlement à ces plaques. Une autre géométrie réalisée par M. Zukowski
[ZUK07a] (Figure I-6), consiste à placer le MCP en amont des gaines de distribution d’air
dans des microcapsules en contact avec des matrices métalliques. Il existe également
plusieurs études de stockage thermique dans les gaines comme [LAZ09a] [LAZ09b]
[ZUK07b]. Le problème majeur de l’intégration des MCP dans les gaines est le coefficient
d’échange médiocre entre l’air et le MCP dans les deux régimes, stockage et déstockage.
L’air doit être soufflé avec une plus grande vitesse pour améliorer ce coefficient. Ces faits
conduisent à une plus grande consommation électrique pour la ventilation, un niveau de bruit
plus élevé et de mauvais rendements exergétique et énergétique.
11
Chapitre I. Le stockage d’énergie intégré dans une PAC
permettant l’effacement total ou partiel de la PAC pendant au moins deux heures. Les
travaux comportent :
une modélisation dynamique du champ de température représentant le matériau à
changement de phase dans différents types de structure d’échange
la modélisation ayant permis d’identifier les conditions adaptées pour les transferts
thermiques avec le MCP, l’air et le fluide frigorigène ; la conception et le
dimensionnement d’échangeurs stockeurs/déstockeurs de chaleur peuvent alors être
entrepris
les essais sur un banc ad hoc pour valider la conception des échangeurs prototypes
enfin l’intégration dans le système PAC, le développement des logiques de pilotage
et les tests en conditions réelles.
Il est à noter que la conductivité thermique des MCP, généralement faible (de l’ordre
de 0,2 W.m -2.K-1 pour les paraffines), amène à concevoir des agencements innovants entre
le MCP et ses surfaces d’échange à densité de flux très différentes.
Les condenseurs prototypes seront testés sur un banc d’essais qui permettra de mesurer les
températures d’intérêt, les écarts de pressions, le débit d’air et les grandeurs permettant de
caractériser la PAC comportant un stockage.
La toute première étape consiste à établir les besoins de chauffage horaires d’une maison
individuelle BBC pendant la saison de chauffe.
12
Chapitre I. Le stockage d’énergie intégré dans une PAC
quelques exigences de moyens, limitées au strict nécessaire, pour refléter la volonté
affirmée de faire pénétrer significativement une pratique (test d’étanchéité à l’air…).
Les exigences de résultats imposées par la RT 2012 sont de trois types.
o L’efficacité énergétique du bâti
L’exigence d’efficacité énergétique minimale du bâti est définie par le coefficient « Bbiomax »
(besoins bioclimatiques du bâti). Cette exigence impose une limitation simultanée du besoin
en énergie pour les composants liés à la conception du bâti (chauffage, refroidissement et
éclairage), imposant ainsi son optimisation indépendamment des systèmes énergétiques mis
en œuvre.
o La consommation énergétique du bâtiment
L’exigence de consommation conventionnelle maximale d’énergie primaire se traduit par le
coefficient « Cepmax », portant sur les consommations de chauffage, de refroidissement,
d’éclairage, de production d’eau chaude sanitaire et d’auxiliaires (pompes et ventilateurs).
Conformément à l’article 4 de la loi Grenelle 1, la valeur du Cepmax s’élève à
50 kWhEP/(m².an) d’énergie primaire en moyenne, modulée selon la localisation
géographique, l’altitude, le type d’usage du bâtiment et la surface moyenne des logements.
Cette exigence impose, en plus de l’optimisation du bâti exprimée par le Bbio, le recours à
des équipements énergétiques performants, à haut rendement.
o Le confort d’été dans les bâtiments non climatisés
A l’instar de la RT 2005, la RT 2012 définit des catégories de bâtiments dans lesquels il est
possible d’assurer un bon niveau de confort en été sans avoir à recourir à un système actif
de refroidissement. Pour ces bâtiments, la réglementation impose que la température la plus
chaude atteinte dans les locaux, au cours d’une séquence de cinq jours très chauds d’été,
n’excède pas un seuil.
II.2.1 Enveloppe
Dans le cadre du projet Stock-aiR2 des simulations ont été réalisées sous Dymola. Bien que
modélisée en mono-zone (nœud d’air d’un volume de 252,15 m3), l’inertie liée aux parois
internes est prise en compte. Par souci de simplicité, le garage attenant et les combles non
aménagés ne sont pas pris en compte. La géométrie finale du modèle est donc
13
Chapitre I. Le stockage d’énergie intégré dans une PAC
parallélépipédique 12,3 m x 8,2 m avec allongement du bâtiment selon les directions Est-
Ouest, les façades étant orientées selon les quatre points cardinaux.
Le Tableau I-3 présente les différentes caractéristiques de l’enveloppe de la maison Mozart
ainsi que les recommandations Effinergie (label BBC [EFF]) pour comparaison. Ce tableau
permet aussi le calcul de la valeur analytique du coefficient de transfert thermique global du
bâtiment Ubat :
(Eq. I-1)
Plafond Laine de verre+ 980 0.191 5.2 6.5 à 10 100.9 1.0 19.3
0.25 100 0.05 5.00
Solives en bois
Plaques de plâtre 0.01 850 800 0.35 0.03
intérieur hconv 10.00 0.10
double vitrage
Vitrages 1.43 0.7 15 1.0 21.45
argon peu émissif
La valeur analytique du Ubat a été vérifiée par une simulation en régime statique. Il est à noter
que cette formule ainsi que la modélisation proposée ne tiennent pas compte des ponts
thermiques linéiques et ponctuels. Ils jouent un rôle minime pour les bâtiments BBC. Dans la
littérature, les valeurs de Ubat recommandées pour être Bâtiment Basse Consommation
(BBC), selon les exigences de 2012, sont comprises entre 0,2 W.m-2.K-1 et 0,4 W.m-2.K-1
[TIT08] [CHL08] [THI08].
La composition des parois internes est présentée Tableau I-4. Ces parois sont de deux types,
les refends et les cloisons légères.
Tableau I-4 – Caractéristiques des parois intérieures
Ep ρ Cp λ S
Type paroi Matériau -3 -1 -1 -1 -1
(m) (kg.m ) (J.kg .K ) (W.m .K ) (m²)
Plâtre 0,01 900 837 0,35
Refends Blocs creux 0,2 1600 920 0,8 19
Plâtre 0,01 900 837 0,35
Plâtre 0,01 900 837 0,35
Cloisons
Air 0,03 1,218 1000 0,19 67,13
légères
Plâtre 0,01 900 837 0,35
14
Chapitre I. Le stockage d’énergie intégré dans une PAC
II.2.2 Ventilation
Le système de ventilation retenu pour cette étude est un système de type hygroréglable pour
lequel le débit de ventilation dépend de l’humidité intérieure. Or, dans cette étude, les
comportements du bâti et des systèmes sont uniquement modélisés en thermique pure (pas
de scénario de production de vapeur ou tampon hygroscopique…). Des études plus
approfondies [KOF09] [AGC] ont montré qu’un système de ventilation hygroréglable sous
conditions standards d’humidité intérieure pouvait être modélisé par un débit de ventilation
moyen d’environ 0,35 vol.h-1.
II.3 Sénarios
II.3.1 Météo
La météo choisie pour cette étude est celle de Trappes. Elle est la météo la plus
représentative du point de vue du nombre d’habitations. Elle correspond à la zone H1a
(Figure I-8) qui est une des zones les plus froides du climat français.
15
Chapitre I. Le stockage d’énergie intégré dans une PAC
apports thermiques sont fixés à 70 % des apports maximum (soit 189 W th pour trois
personnes). Le reste du temps, les apports sont maximums (3 * 90 W th = 270 W th).
o Usages spécifiques
Les usages spécifiques sont définis par un coefficient de 5,7 W.m-2habitable et par une
corrélation en fonction de la présence et de l’activité des occupants. Lors des périodes de
vacances, ils sont fixés à 0 %. Lors des absences en semaine et pendant les périodes de
sommeil (de 22h00 à 6h00) les apports sont de 20 %. Le reste du temps, les apports sont
fixés à 100 %.
o Eclairage
Les apports thermiques liés à l’éclairage sont aussi corrélés aux activités et à la présence
des occupants. Pendant les périodes de vacances, d’absence en semaine et de sommeil, les
apports sont nuls. Le reste du temps ils sont estimés à 2 W.m-2habitable.
(Eq. I-2)
(Eq. I-3)
En considérant une température de consigne fixe plus réaliste de 21 °C (au lieu de 19 °C) et
les mêmes scénarios d’apports internes et de réduits que précédemment, le besoin en
énergie finale de chauffage sera de :
(Eq. I-4)
16
Chapitre I. Le stockage d’énergie intégré dans une PAC
Les calculs RT avec température intérieure de 19 °C sous-estiment les besoins de chauffage
par rapport à la réalité du terrain (température moyenne plutôt autour de 21 °C). De même, il
est plus sûr pour le dimensionnement du système de privilégier la solution avec les
scénarios de réduits évoqués précédemment plutôt que ceux réglementaires du paragraphe.
Par ailleurs, les scénarios d’apports internes sont toujours pris en compte. La simulation
retenue correspond donc au besoin en énergie finale de chauffage Besoin3.
35 25.0%
30
20.0%
25
Nombre de jour [ j ]
Occurence [%]
15.0%
20
23.6% 22.9%
15 20.7%
10.0%
10
11.4% 12.1%
5.0%
5
4.3% 4.3%
0.7%
0 0.0%
0
00
0
50
50
00
50
00
50
00
10
0-
-1
-2
-2
-3
-3
-4
0-
00
00
00
00
00
00
50
10
15
20
25
30
35
Afin d’effacer le système de chauffage sans perte de confort, i.e. maintien de la température
intérieure à la consigne, le stockage doit être en mesure de fournir la puissance thermique
nécessaire. Pour ce faire, une analyse en puissance doit être réalisée.
17
Chapitre I. Le stockage d’énergie intégré dans une PAC
II.5.2 Puissance de déstockage
Afin de quantifier la puissance de déstockage nécessaire pour maintenir la température
intérieure à sa consigne, une analyse sur le jour le plus froid doit être menée. Il est défini
comme le jour nécessitant la plus grande énergie de chauffage entre 18h00 et 20h00. En
théorie, la puissance de chauffage maximale le jour le plus froid n’est pas nécessairement la
puissance de chauffage maximale sur la période de 18h00 à 20h00, cependant elle en
constitue une bonne indication. Ici, c’est le même jour.
La Figure I-10 représente l’évolution des températures intérieure et extérieure le jour le plus
froid. La courbe rouge représente la puissance thermique nécessaire au système de
chauffage pour maintenir la température intérieure à sa valeur de consigne. Il est à noter que
le 362ème jour de l’année est un jour de la semaine de vacances d’hiver, par conséquent la
température de consigne est de 16 °C (courbe verte).
La Figure I-10 permet aussi de mettre en évidence la puissance de déstockage nécessaire
qui est la valeur maximale de la puissance de chauffage sur la période de 18h00 à 20h00. La
puissance de déstockage permettant de maintenir la température de consigne doit être, au
minimum, d’environ 1 910 W.
20 2500
15
2000
Effacement 1500
Text
5
Tint
Puissance chauffage 1000
0
00
00
00
00
00
00
00
00
00
10 0
00
11 0
12 0
13 0
14 0
15 0
16 0
17 0
18 0
19 0
20 0
21 0
22 0
23 0
0
0
:0
:0
:0
:0
:0
:0
:0
:0
:0
:0
:0
:0
:0
:0
0:
1:
2:
3:
4:
5:
6:
7:
8:
9:
0:
500
-5
-10 0
Temps [h]
Figure I-10 – Evolution des températures et puissance le jour le plus froid (362ème jour)
La quantité d’énergie nécessaire pour effacer le système de chauffage le jour le plus froid est
représentée par l’aire orange. Cette quantité d’énergie est l’intégrale de la puissance de
chauffage entre 18h00 et 20h00. Cette quantité d’énergie est égale à 3 822 Wh le jour le
plus froid.
Le Tableau I-5 illustre la quantité d’énergie nécessaire pour effacer le système de chauffage
entre 18h00 et 20h00 lors des vingt jours les plus froids de l’année (définis comme les jours
où les besoins de chauffage entre 18h00 et 20h00 sont les plus importants).
Les résultats figurant Tableau I-5 ont été obtenus avec les mêmes hypothèses et scénarios
que ceux utilisés pour les simulations.
Afin de réduire le dimensionnement de la puissance du système, une autre possibilité
consiste à prendre le jour le plus froid ayant un besoin d’effacement inférieur à 3 000 Wh.
D’après le Tableau I-5, c’est le 16ème jour de l’année. La Figure I-11 représente l’évolution
18
Chapitre I. Le stockage d’énergie intégré dans une PAC
des températures intérieure et extérieure le jour le plus froid. La courbe rouge représente la
puissance thermique nécessaire au système de chauffage pour maintenir la température
intérieure à sa valeur de consigne. Le 16ème jour de l’année ne fait pas partie de la période
de vacances et par conséquent la température de consigne est de 21 °C. Il résulte que la
puissance de déstockage doit être d’environ 1 640 W. L’énergie de stockage est de
2 959 kWh. Il est à noter que la forte variation de la puissance de chauffage est due
majoritairement aux scénarios d’apports internes.
Tableau I-5 – Besoins énergétiques d’effacement sur les 20 jours les plus froids
Besoin Effacement Text 18h00 Text 19h00 Text 20h00
Jour (W.h) (°C) (°C) (°C)
362 3822 -0,5 -0,4 -0,2
365 3799 1,9 1,6 1,2
360 3741 -0,7 -0,9 -1,1
359 3691 -3,6 -3,6 -3,6
361 3611 -0,2 -0,3 -0,5
15 3582 -3,4 -3,1 -2,9
363 3097 3,4 3,4 3,4
16 2959 0,6 0,6 0,5
12 2952 -5,2 -5,6 -6,0
14 2912 -1,3 -1,4 -1,5
45 2796 -3,2 -3,7 -4,1
357 2541 -0,5 -0,8 -1,0
35 2513 -0,5 -0,3 0,0
9 2493 1,4 1,2 1,0
355 2442 -0,7 -0,6 -0,6
46 2416 -0,7 -0,5 -0,4
37 2407 1,9 1,7 1,5
40 2400 -0,3 -0,4 -0,4
347 2375 1,0 1,2 1,6
10 2329 2,1 2,0 1,8
25 2500
20
2000
Puissance thermique [W]
15
Température [°C]
1500
10
Effacement
5 Text
1000
Tint
0 Puissance chauffage
00
00
00
00
00
00
00
00
00
10 00
00
11 0
12 0
13 0
14 0
15 0
16 0
17 0
18 0
19 0
20 0
21 0
22 0
23 0
0
:0
:0
:0
:0
:0
:0
:0
:0
:0
:0
:0
:0
:0
:0
500
0:
1:
2:
3:
4:
5:
6:
7:
8:
9:
0:
-5
-10 0
Temps [h]
Figure I-11 – Evolution des températures et puissance le 16ème jour
19
Chapitre I. Le stockage d’énergie intégré dans une PAC
II.5.3 Stratégies de fonctionnement
Le dimensionnement de l’effacement, tel qu’il est présenté ici, laisse apparaître plusieurs
stratégies tant en terme de puissance de déstockage que d’énergie d’effacement :
stratégie 1 : effacement de 2 h garanti et pas d’inconfort,
stratégie 2 : effacement de 2 h garanti pour 95 % des jours et aucun inconfort,
stratégie 3 : effacement de 2 h garanti pour 95 % des jours et pas d’inconfort ces
mêmes jours.
La stratégie 1 repose sur un dimensionnement du système « au maximum », ainsi aucun
inconfort n’est autorisé, ni une durée d’effacement inférieure à 2 h. Pour ce faire, l’énergie de
stockage doit être d’environ 3 820 Wh et la puissance de déstockage doit être d’environ
1 910 W.
Pour la stratégie 2, il est souhaité de s’effacer pendant 2 h (au minimum) 95 % des jours
sans perte de confort. Les autres jours, il n’y aura pas d’inconfort mais la période
d’effacement sera nécessairement inférieure à 2 h. L’énergie stockée doit être d’environ
3 000 Wh et la puissance de décharge d’environ 1 910 W.
Dans la stratégie 3, un inconfort et une durée d’effacement inférieure à 2 h sont autorisés
dans les 5 % des jours les plus froids. Le reste du temps il n’y a pas d’inconfort et la durée
d’effacement de 2 h est garantie. L’énergie stockée doit être d’environ 3 000 Wh et la
puissance de décharge d’environ 1 640 W.
D’autres stratégies de dimensionnement sont envisageables suivant les paramètres les plus
critiques pour le système de stockage MCP (quantité d’énergie à stocker et/ou puissance de
décharge).
(Eq. I-5)
20
Chapitre I. Le stockage d’énergie intégré dans une PAC
III. Conclusions
Le chauffage domestique a une grande influence sur la demande d’électricité surtout les
heures de pointe durant l’hiver. D’où l’idée d’intégrer un système de stockage thermique
journalier dans les pompes à chaleur, qui si cette technologie est utilisée massivement pour
le chauffage domestique permettra alors de réduire de pic de demande électrique entre 18
et 20h et cela en décalant le temps de production de celui de consommation.
Le système PAC avec stockage intégré au condenseur n’est pas forcément moins efficace
qu’une PAC sans stockage. Ceci est dû au fait que les heures creuses de demande
d’électricité, qui sont bien sûr les heures les plus favorables pour un fonctionnement en
mode stockage de la PAC correspondent aux heures pendant la journée aux heures où la
température extérieure est la plus élevée et donc ou le COP de la PAC est potentiellement
le plus élevé. Au chapitre 5, des simulations réalistes ont montré l’impact du
stockage/déstockage sur l’efficacité énergétique de la PAC pour différents scénarios de
régulation et différentes journées types.
La réglementation thermique française dite RT 2012 impose une consommation en énergie
primaire de 50 kWhEP.m-2.an-1 pour couvrir les cinq besoins (chauffage, Eau chaude
Sanitaire, ventilation, climatisation, éclairage) des logements neufs. Une maison Mozart de
100 m² habitable conçue pour respecter la RT 2012, nécessitera entre 24 et 36 kWh.m-2.an-1
en énergie finale pour le chauffage, selon les scénarios et les températures de consigne
considérés. Le jour le plus froid a été identifié comme le jour ayant les besoins les plus
élevés entre 18h et 20h (les heures d’effacement).
L’effacement entre 18h et 20h peut être réalisé avec une PAC stockant 800 Wh apte à
déstocker une puissance de 400 W en continu. Bien sûr, cet effacement pourra durer plus
longtemps pour les besoins thermiques plus faibles.
21
Chapitre I. Le stockage d’énergie intégré dans une PAC
22
Chapitre II. Matériaux à Changement de Phase – Paraffines, Sels ou Eutectiques
I. Introduction
Un matériau à changement de phase - ou MCP - est un matériau apte à changer d'état
physique dans une plage de température restreinte. Dans cet intervalle de température, les
changements de phase vont s’effectuer dans le sens de la fusion ou de la solidification
suivant l’apport ou, au contraire, l’emport d’énergie.
Tout matériau, solide ou liquide (ou gazeux) possède une capacité à stocker ou céder de
l'énergie sous forme de chaleur. On distingue deux types de transfert de chaleur (ou transfert
thermique) :
Le transfert thermique par Chaleur Sensible (CS) - Dans ce cas, le matériau peut céder
ou stoker de l'énergie thermique par variation de sa propre température, sans pour
autant changer d'état. La grandeur utilisée pour quantifier la CS échangée par un
matériau est la chaleur massique, notée Cp et exprimée en J.kg-1.K-1 [LAN83].
(Eq. II-1)
Le transfert thermique par Chaleur Latente (CL) - Le matériau peut stocker ou céder de
l'énergie par changement d'état, tout en conservant une température constante, celle du
changement d'état. La grandeur utilisée pour quantifier la CL échangée par un matériau
est la chaleur latente de changement de phase notée hLS pour un changement de phase
liquide/solide, et hLV pour un changement de phase liquide/vapeur. Celle-ci est exprimée
en J.kg-1 [ABH83].
(Eq. II-2)
Dans un transfert thermique où le MCP est endothermique, une première phase est marquée
par un échange par chaleur sensible (SCS), la température du MCP s'élève par absorption
de la chaleur. Contrairement aux SCS classiques, lorsque les MCP atteignent la température
à laquelle ils changent de phase (leur température de fusion), ils absorbent de grandes
quantités de chaleur à cette même température presque constante. Le MCP continue à
absorber la chaleur sans une augmentation significative de sa température jusqu'à ce que
tout le matériau se transforme en liquide. Lorsque la température ambiante autour d'un
matériau liquide chute, le MCP se solidifie et libère sa chaleur latente stockée.
23
Chapitre II. Matériaux à Changement de Phase – Paraffines, Sels ou Eutectiques
(Eq. II-3)
Le stockage thermique par chaleur latente peut être obtenu par des transformations solide-
solide, solide-liquide, solide-gaz et liquide-gaz.
En transition solide - solide, la chaleur est stockée lorsque le matériau change de structure
cristalline. Les changements de phase solide-solide sont généralement très lents et
présentent une chaleur de transformation faible. Par contre, les MCP solide-solide réduisent
les contraintes dues aux changements de volume et offrent une plus grande souplesse de
conception [CLA81]. Les matériaux les plus prometteurs sont les solutions solides
organiques de Penta-érythritol, Penta-glycerine, Li2SO4 et KHF2 [PAP99]. Un mur Trombe
rempli de ces matériaux présente de meilleures performances d’inertie qu'un mur Trombe
béton ordinaire.
Les changements de phase solide-gaz et liquide-gaz génèrent une chaleur latente élevée.
Par contre, cette transition est accompagnée d’un grand changement de volume, ce qui rend
leur utilisation pour le stockage thermique complexe et peu pratique.
Les transformations solide-liquide, comparativement aux liquide-gaz, ont une chaleur latente
plus faible. Cependant, ces transformations impliquent seulement un petit changement en
volume (de l'ordre de 15 %). Les transitions solide-liquide sont prouvées être
économiquement intéressantes pour l’utilisation dans les systèmes de stockage d'énergie
thermique. Par contre, les MCP mêmes ne peuvent pas être utilisés en tant que milieu de
transfert de chaleur. L'échangeur de chaleur à utiliser doit être conçu spécialement pour
l’application, en vue de la faible diffusivité thermique des MCP en général. Ceci constitue le
défi majeur des systèmes utilisant ces matériaux. Le changement de volume du MCP lors de
la fusion nécessite également une conception spéciale du conteneur. Il devrait être en
mesure d'absorber ces variations de volume et devrait également être compatible avec le
MCP utilisé. Tout système de stockage par chaleur latente doit donc posséder au moins les
trois propriétés suivantes : un MCP adapté à son point de fusion dans la plage de
température souhaitée, une surface d'échange de chaleur appropriée et un récipient
compatible avec le MCP.
De nombreux MCP solide-liquide sont disponibles dans les plages de température d’intérêt.
Pour les températures de 20 °C à 50 °C associées au confort thermique humain, certains
MCP sont très efficaces. Ils stockent de 5 à 14 fois plus de chaleur par unité de volume de
stockage que les matériaux classiques à chaleur sensible tels que l'eau, la maçonnerie, ou
les roches.
24
Chapitre II. Matériaux à Changement de Phase – Paraffines, Sels ou Eutectiques
(Eq. II-4)
Cette règle est illustrée pour les corps purs ou pour les solutions. L’hypothèse prise est que
la pression est constante et égale à la pression atmosphérique. D’après la règle des phases
ci-dessus, toute modification de la pression devrait entraîner une modification de la
température ou des concentrations d’équilibre. Il faut remarquer que les propriétés physiques
des phases condensées varient relativement peu avec la pression
25
Chapitre II. Matériaux à Changement de Phase – Paraffines, Sels ou Eutectiques
III.2.1 Solutions salines
Les solutions salines sont composées d’un solvant et d’un soluté. Le solvant le plus utilisé
est l’eau où la solution sera dite solution aqueuse. Ces solutions admettent deux degrés de
liberté, la pression sera considérée imposée égale à la pression atmosphérique. La
température de fusion et la composition physique (fraction des phases) seront données alors
par un diagramme de phases spécifiques à chaque solution [DUM]. La figure II-1 représente
le diagramme de phases isobare d’un mélange binaire où les phases solides sont
immiscibles. C’est le cas des mélanges eau-sel.
T TF B
T0 Solution
GB
TFA
TF(x)
TF(x0) GA
TE
E
0 x0 xE x 1 x
A B
Figure II-1 – Diagramme de phases schématique d’un mélange binaire avec phases solides
immiscibles
(Eq. II-5)
Un diagramme de phase d’une solution binaire est constitué de deux liquidus et deux
solidus. Les liquidus GA et GB présentent, pour une fraction massique en constituant B, la
température d’équilibre avec les composants A et B respectivement. Les solidus, qui
représentent les températures d’équilibre en fonction des concentrations en phase solide,
ont ici une forme un peu particulière puisque les phases solides sont pures. Elles sont donc
constituées des deux droites verticales à x = 0 et x = 1, limitées vers le haut par les
températures de fusion des corps purs TFA et TFB.
Un point très particulier d’un diagramme de phase est celui de l’intersection des liquidus : le
point E. En effet, à la température TE, il y a équilibre entre trois phases : le solide A pur, le
solide B pur et une solution à la fraction molaire xE. Dans ce cas, appelé mélange eutectique,
trois phases existent alors ϕ = 3 entraînant w = 1. Si la pression est fixée, toutes les autres
variables intensives sont fixées par la nature conjointe des constituants A et B. Ainsi sont
fixées la température TE et la concentration de la solution xE. Cela veut dire que, tant qu’il y a
présence des trois phases, la température de fusion reste rigoureusement fixe : c’est la
fusion eutectique. Au cours de cette fusion, les quantités de A et de B nécessaires se
transforment pour que la fraction molaire du liquide soit exactement xE. C’est un autre moyen
d’obtenir un MCP fondant à une température fixe. Les MCP eutectiques aqueux sont utilisés
pour le stockage du froid à des niveaux de température inférieurs à 0 °C. Comme les
températures de fusion des sels dissouts sont très élevées, la concentration eutectique est
généralement basse. Il faut donc assez peu de sel (qui coûte le plus cher) pour disposer d’un
MCP. De plus, la chaleur latente de fusion eutectique sera très proche de celle de fusion de
26
Chapitre II. Matériaux à Changement de Phase – Paraffines, Sels ou Eutectiques
la glace (qui est très intéressante par ailleurs). Pour avoir différents MCP à différentes T E, il
suffira de changer le sel utilisé.
Le processus de fusion d’une solution à une concentration x0 différente de celle de
l’eutectique se déroule d’une manière plus compliquée. En considérant que x0 < xE, au
départ deux phases solides immiscibles A et B sont présentes, dont les proportions sont
respectivement (1 – x0) et x0. Cet état restera tant que la température du mélange est
inférieure à TE. Une fois cette dernière température atteinte, une solution de concentration xE
commence à se former en présence des quantités des solides A et B. La température du
mélange ne dépassera pas TE tant qu’il y a encore du solide B (x0 < xE). Une fois tout le
solide B fondu, la hausse de température est accompagnée par une augmentation de la
masse liquide (diminution de la concentration de la solution en B) et une diminution de la
masse du solide A. Lorsque le mélange atteint la température d’équilibre (donnée par le
liquidus GA pour une concentration x0) il sera complètement à l’état liquide de concentration
massique x0.
Le thermogramme de la figure II-2 est obtenu en calorimétrie au cours de la fusion d’une
solution cristallisée.
Puissance
T F (x0)
TE
27
Chapitre II. Matériaux à Changement de Phase – Paraffines, Sels ou Eutectiques
(Eq. II-6)
En réalité, ΔHdil (TE) est négligeable (de l’ordre de 1 %), la relation II-7 sera réduite à :
(Eq. II-7)
Il faut, par contre, tenir compte de la variation de la chaleur latente de fusion de la glace avec
la température parce que celle-ci est relativement importante du fait de la différence très
importante entre les capacités thermiques massiques de l’eau liquide et de la glace. Cette
différence est beaucoup moins marquée pour le sel. En effet
(Eq. II-8)
Il est bien entendu que CpL eau ne peut que représenter la capacité thermique massique de
l’eau surfondue (ou une extrapolation de ces valeurs pour les basses températures). En
supposant ces grandeurs sensiblement indépendantes de la température et en prenant leur
valeur à 0 °C il vient :
(Eq. II-9)
Comme CpL eau – CpS eau > 0, la relation II-9 confirme que la chaleur latente à l’eutectique est
inférieure à la chaleur latente de fusion de la glace. Cette valeur est plus basse que celle de
l’eau pure mais est quand même du même ordre de grandeur, ce qui justifie que l’utilisation
de solutions aqueuses eutectiques est toujours intéressante sur le plan énergétique. De plus,
elles ont un coût relativement faible.
(Eq. II-10)
Chaque particule de sel est liée à n molécules d’eau. Il faut souligner que la réaction II-10 est
fortement exothermique.
28
Chapitre II. Matériaux à Changement de Phase – Paraffines, Sels ou Eutectiques
T (°C)
60
0 H2O
Solution
40
10 H2O
20
Na2SO4, 10 H2O
0
E 0,1 0,2 0,3 0,4
x
-20
Le diagramme de phases de tels types de sels est de même type que celui de la figure II-3
qui présente le cas du sulfate de sodium, dit sel de Glauber, qui a été le premier MCP
proposé pour le stockage de chaleur [LAN83].
Ce graphe présente trois liquidus :
A. La ligne horizontale au point P qui représente l’équilibre :
(Eq. II-11)
29
Chapitre II. Matériaux à Changement de Phase – Paraffines, Sels ou Eutectiques
T2(x0)
GL dT(x0)
T1(x0)
GS
0 x0 1 x
Figure II-4 – Diagramme de phases avec solution solide
(Eq. II-12)
T A
Refroidissement
du liquide
Solidification
TF
Refroidissement
du solide Stable
Surfusion Métastable
30
Chapitre II. Matériaux à Changement de Phase – Paraffines, Sels ou Eutectiques
Le liquide est dit surfondu ou métastable, état qui peut se maintenir longtemps (des heures
ou des jours) et à des températures assez basses. La rupture de la métastabilité ou le
déclenchement de la cristallisation se fait de deux manières, soit provoquée, soit spontanée.
1. Provoquée :
par incorporation d’additifs tensio-actifs
par facilitation de la nucléation au moyen de cristaux solides stables autour de la
température de cristallisation : ceux-ci servent de points d'accroche pour la
cristallisation du matériau
par maintient d’une zone froide à l'intérieur du matériau
par mise en vibration du liquide métastable.
2. Spontanée : lorsque le refroidissement sera suffisant pour déclencher la cristallisation.
La rupture de surfusion provoquée ne sera pas forcément facile à mettre en œuvre et
entraîne une opération supplémentaire à gérer.
Plusieurs lois de cristallisation spontanée sont proposées :
Caractère erratique ou stochastique : plusieurs échantillons liquides apparemment
identiques et produit dans un même lot, ne cristallisent pas à la même température
au cours d’un même refroidissement.
Analyse statistique : il est possible de définir l’intervalle de température où se produit
la majorité des ruptures de surfusion ainsi qu’une température moyenne de
cristallisation, T*, qui sera représentative de l’importance de la surfusion.
Probabilité de rupture de métastabilité : il est aussi possible de définir une probabilité
de cristallisation de l’échantillon en fonction de la température par unité de temps.
(Eq. II-13)
31
Chapitre II. Matériaux à Changement de Phase – Paraffines, Sels ou Eutectiques
32
Chapitre II. Matériaux à Changement de Phase – Paraffines, Sels ou Eutectiques
33
Chapitre II. Matériaux à Changement de Phase – Paraffines, Sels ou Eutectiques
Tableau II-2 – Liste de quelques MCP organiques [FAR04] [ZAL03] [SHA99b] [HAS98] [HAL71]
Nom Tfusion (°C) hLS (kJ.kg-1) Nom Tfusion (°C) hLS (kJ.kg-1)
N-Tetradeca 5,5 226 N-Penta Cosane 53,7 164
Acide Formique 7,8 247 Myristic Acid 54,0 199
N-Pentadeca 10,0 205 Oxolate 54,3 178
Acide Acétique 16,7 273 Tristearin 54,5 191
N-Hexadeca 16,7 237 O-Xylene 55,0 121
Caprilone 40,0 260 β-Chloroacetic 56,0 147
Docasyle 40,0 201 N-Hexacosane 56,3 255
N-Henicosan 40,5 161 Nitro Naphthalene 56,7 103
Phénol 41,0 120 α-Chloroacetic 61,2 130
N-Laurique 43,0 183 N-Octacosane 61,4 134
P-Joluidine 43,3 167 Palmitic Acid 61,8 164
Cynamide 44,0 209 Bees Wax 61,8 177
N-Docosane 44,5 157 Glyolic Acid 63,0 109
N-Tricosane 47,6 130 P-Bromophenol 63,5 86
Hydrocinna 48,0 118 Azobenzene 67,1 121
Cetyl 49,3 141 Acrylic Acid 68,0 115
O-Nitroanili 50,0 93 Dintro Toluene 70,0 111
Camphene 50,0 239 Acide Phénylacetique 76,7 102
Diphenyl 52,9 107 Thiosinamine 77,0 140
34
Chapitre II. Matériaux à Changement de Phase – Paraffines, Sels ou Eutectiques
Tableau II-3 – Liste de quelques MCP inorganiques [LAN83] [ABH83] [FAR04] [ZAL03]
[SHA99b] [HAS98] [HAL71]
Nom Tfusion (°C) hLS (kJ.kg-1) Nom Tfusion (°C) hLS (kJ.kg-1)
H2O 0,0 333 BI3 31,8 10
POCl3 1,0 85 SO3(β) 32,3 151
D2O 3,7 318 TiBr4 38,2 23
SbCl5 4,0 33 H4P2O6 55,0 213
H2SO4 10,4 100 SO3 (γ) 62,1 331
IC l (β) 13,9 56 SbCl3 73,4 25
MOF6 17,0 50 NaNO3 307 172 - 199
SO3(α) 17,0 108 KNO3 333 - 380 116 - 266
IC l (α) 17,2 69 KOH 380 149
P4O6 23,7 64 MgCl2 714 - 800 452 - 492
H3PO4 26,0 147 NaCl 800 - 802 466 - 492
Cs 28,3 15 Na2CO3 854 275
Ga 30,0 80 KF 857 452
AsBr3 30,0 38 K2CO3 897 235
Cependant aucun matériau ne pourra avoir toutes les propriétés requises pour un système
de stockage. Pour cela il faut adapter une conception adéquate du système pour essayer de
compenser les propriétés non idéales du MCP choisi. Par exemple, des ailettes métalliques
peuvent être utilisées pour augmenter la conductivité thermique du MCP, le sous-
refroidissement peut être largement réduit ou même supprimé par l'introduction d'un agent
de nucléation dans le matériau de stockage, et une fusion non congruente peut être inhibée
par l'utilisation d’une épaisseur convenable du MCP. Les comportements thermique et
chimique sont différents pour deux MCP ayant une même température de fusion, une même
chaleur latente et étant de deux sous-groupes différents. D’où la conception d’un système
est aussi liée au sous-groupe du MCP.
V.1 Paraffines
Les paraffines normales de type CnH2n +2 sont une famille d'hydrocarbures saturés avec des
propriétés très similaires. Les paraffines entre C5 et C15 sont des liquides, et le reste sont des
solides cireux. La paraffine cire est le MCP le plus utilisé pour les systèmes de stockage
d’énergie thermique ([LAN83], [HAL71]). Les cires se composent principalement
d'hydrocarbures à chaîne droite qui ont des températures de fusion allant de 23 à 67 °C
([ABH83]). La paraffine commerciale est obtenue par distillation du pétrole, elle n'est pas une
substance pure, mais une combinaison de différents hydrocarbures. En général, plus la
longueur moyenne de la chaîne hydrocarbonée est grande, plus la température de fusion est
élevée ([HIR94]). Les propriétés de plusieurs paraffines sont présentées dans le tableau II.4
([DAW09], [TAN09]) Les paraffines sont facilement disponibles auprès de nombreux
fabricants et sont généralement plus chères que les hydrates de sel ([LAN83], [HAL71]).
35
Chapitre II. Matériaux à Changement de Phase – Paraffines, Sels ou Eutectiques
pas de tendance à la surfusion et n’ont donc pas besoin d’agents de nucléation ([LAN83],
[BUD94]). Elles sont sûres et non-réactives ([HAL71], [HAS98]). Elles sont compatibles avec
tous les conteneurs en métal et facilement incorporées dans des systèmes de stockage de
chaleur ([LAN83]). Des précautions doivent être prises lors de l'utilisation des récipients en
plastique comme les paraffines ont tendance à infiltrer et adoucir certains plastiques
([LAN83]).
36
Chapitre II. Matériaux à Changement de Phase – Paraffines, Sels ou Eutectiques
Contrairement aux sels hydratés, les paraffines commerciales n'ont généralement pas de
points de fusion bien définis, puisqu’elles sont composées d’un mélange de matières. Les
paraffines sont inflammables, mais cela peut être limité par le confinement dans un récipient
approprié ([HIR94], [HAL71], [HAL98]). Ce risque doit être soigneusement analysé si de très
grandes quantités de paraffines sont intégrées dans des parois.
V.2 Non-paraffines
Il s'agit de la catégorie la plus nombreuse de matériaux candidats pour le stockage de
chaleur latente. Lane ([LAN78], [LAN83], [LAN89]), Abhat [ABH83] et Buddhi et al. ([BUD94])
ont mené des enquêtes approfondies sur les matières organiques et ont identifié un certain
nombre d'esters, d'acides gras, des alcools et des glycols appropriés pour le stockage de
l'énergie. Ces matières organiques sont classées en sous-groupes : acides gras et autres
matières organiques non-paraffines. Les composés organiques non-paraffiniques sont des
MCP ayant des propriétés très variées. Chacun de ces matériaux a ses propres propriétés,
contrairement aux paraffines, qui ont des propriétés très similaires. Ces matériaux sont
inflammables et ne doivent pas être exposés à des températures trop élevées, des flammes
ou des agents oxydants. Les non-paraffines sont présentées Tableau II-5.
Tableau II-5 – Propriétés de quelques non- paraffines [LAN83] [ABH83] [FAR04] [ZAL03]
[SHA99b] [HAS98] [HAL71]
Nom Tfusion Masse volumique hLS (kJ.kg-1)
(°C) (kJ.kg-3)
Acide Formique 7,8 1 226,715C 247
Acide Acétique 16,7 1 05020C 187
Glycérine 17,9 1 26020C 198,7
Lithium Chloride 21 - 188
Ethanolate
Polyéthylène Glycol 600 20-25 1 10020C 146
D – Acide Lactique 26 1 24915C 184
1-3 Méthyle Pentacosane 29 - 197
Camphenilone 39 - 205
Docasyl Bromide 40 - 201
Caprylone 40 - 259
Heptadecanone 41 - 201
1-Cyclohexylooctadecane 41 - 218
4-Heptadacanone 41 - 197
Cyanamide 44 1 08020C 209
Méthyl Eicosanate 45 85179C 230
3-Heptadecanone 48 - 218
2-Heptadecanone 48 - 218
Camphene 50 84254C 238
9-Heptadecanone 51 - 213
Méthyl Behenate 52 - 234
Acide Pentadecanoique 52,5 - 178
Acide Hypophosphorique 55 - 213
Acide Chloroacetique 56 1 58020C 130
Trimyristin 33-57 86220C 201-213
Acide Heptaudecanoique 60,6 - 189
Cire d’abeilles 61,8 950 177
37
Chapitre II. Matériaux à Changement de Phase – Paraffines, Sels ou Eutectiques
Tableau II-5 – Propriétés de quelques non- paraffines [LAN83] [ABH83] [FAR04] [ZAL03]
[SHA99b] [HAS98] [HAL71] (suite)
Nom Tfusion Masse volumique hLS (kJ.kg-1)
(°C) (kJ.kg-3)
Acide Glycolique 63 - 109
Cire Oxazoline TS 970 74 - -
Acide Arachique 76,5 - 227
Bromcamphor 77 1 44981C 174
Durene 79,3 83820C 156
Acétamide 81 1 159 241
Méthyl Brombrenzoate 81 - 126
Alpha Napthol 96 1 09598C 163
Acide Glautarique 97,5 1 429 156
p-Xylene Dichloride 100 - 138,7
Méthyl Fumarate 102 1 045 242
Catechol 104,3 1 37015C 207
Quinone 115 1 31820C 171
Acetanilide 115 1 2104C 142
Succinic Anhydride 119 1 104 204
Acide Benzoique 121,7 1 26615C 142,8
Stibene 124 1 16415C 167
Benzamide 127,2 1 341 169,4
Phenacetin 137 - 136,7
Alpha Glucose 141 1 544 174
Acétyle – p- Toluidene 146 - 180
Phenylhdrazone 155 - 134,8
Benzaldehy
Acide Salicylique 159 1 44320C 199
Benzanilide 161 - 162
O-Mannitol 166 1 48920C 294
Hydroquinone 172,4 1 35820C 258
p-Acide Aminobenzoique 187 - 153
38
Chapitre II. Matériaux à Changement de Phase – Paraffines, Sels ou Eutectiques
utilisateurs et comme les utilisations sont très larges et beaucoup ne requièrent pas un haut
niveau de pureté, les fabricants d’acide stéarique n’offrent pas de garantie pour la précision
des températures de changement de phase
C'est la raison pour laquelle les chaleurs latentes de fusion de l'acide stéarique diminuent
généralement avec le nombre de cycles thermiques.
Tableau II-6 – Propriétés de quelques acides gras [LAN83] [ABH83] [FAR04] [ZAL03]
[SHA99b] [HAS98] [HAL71]
Nom Tfusion Masse volumique Conductivité hLS (kJ.kg-1)
(°C) (kJ.kg-3) thermique
(W.m-1.kg-1)
Propyl Palmiate 10 - - 186
Isopropyl Palmiate 11 - - 100
Oleic Acid 13,5 – 16,3 86360C - -
Isopropyl Stearate 14 - 19 - - 140 - 142
Caprylic Acid 16 90130C 0,14939C 148
Butyl Stearate 19 - - 140
Dimethyl Sabacate 21 - - 120 - 135
25C 1C 39C 61C
Myristic Acid + 24 888 , 1018 0,164 , 0,154 147,7
Capric
Vinyl Stearate 27 - 29 - - 122
Methyl Palmitate 29 - - 205
Capric Acid 32 87845C 0,15338C 152,7
Erucic Acid 33 85370C - -
Lauric Acid 42 - 44 86260C, 100724C - 178
Elaidic Acid 47 85179C - 218
Pelargoinc Acid 48 - - -
Myristic Acid 49 - 51 86155C - 205
Palmitic Acid 64 85065C 0,16268C 185,4
Stearic Acid 69 84870C 0,17270C 202,5
Valporic Acid 120 - - -
39
Chapitre II. Matériaux à Changement de Phase – Paraffines, Sels ou Eutectiques
Tableau II-7 – Propriétés de quelques sels hydratés [LAN83] [ABH83] [FAR04] [ZAL03]
[SHA99b] [HAS98] [HAL71] [HIR94]
Nom Tfusion Masse Conductivité hLS Fusion
(°C) volumique thermique (kJ.kg-1)
(kJ.kg-3) (W.m-1.kg-1)
LiClO3.3H2O 8 - - 253 congruente
NH4Cl.Na2SO4.10H2O 11 - - 163 -
K2HO4.6H2O 14 - - 108 congruente
NaCl.Na2SO4.10H2O 18 - - 286 -
KF.4H2O 18 - - 330 congruente
K2HO4.4H2O 18.5 1 44720C, - 231 -
1 45518C
Mn(NO3)2.6H2O 25 1 73820C [12] - 148 -
LiBO2.8H2O 25,7 - - 289 -
FeBr3.6H2O 27 - - 105 -
CaCl2.6H2O 29 - 30 1 56232C, 0,561 , 1,00823C
61C
170 - 192 incongruente
1 80224C
LiNO3.3H2O 30 - - 189 - 296 congruente
Na2SO4.10H2O 32 1 48524C 0,544 251 - 254 incongruente
Na2CO3.10H2O 33 - 36 1 442 - 247 incongruente
KFe(SO4)2.12H2O 33 - - 173 incongruente
CaBr2.6H2O 34 1 95635C, - 115 - 138 -
2 19424C
LiBr.2H2O 34 - - 124 incongruente
Na2HPO4.12H2O 35 1 522 - 256 - 281 incongruente
Zn(NO3)2.6H2O 36 1 82836C, 0,464 , 0,46961C
40C
134 - 147 congruente
1 93724C
Mn(NO3)2.4H2O 37 - - 115 -
FeCl3.6H2O 37 - - 223 congruente
CaCl2.4H2O 39 - - 158 incongruente
CoSO4.7H2O 40.7 - - 170 -
CuSO4.7H2O 40.7 - - 171 -
KF.2H2O 42 - - 162 - 266 congruente
MgI2.8H2O 42 - - 133 -
CaI2.6H2O 42 - - 162 -
Ca(NO3)2.4H2O 43 - 47 - - 106 - 140 congruente
Zn(NO3)2.4H2O 45 - - 110 -
K3PO4.7H2O 45 - - 145 -
Fe(NO3)3.9H2O 47 - - 155 - 190 incongruente
Mg(NO3)3.4H2O 47 - - 142 -
Na2SiO3.5H2O 48 - - 168 -
Na2HPO4.7H2O 48 - - 135 - 170 incongruente
Na2S2O3.5H2O 48 1 600 - 209 -
K2HPO4.3H2O 48 - - 99 -
MgSO4.7H2O 48,4 - - 202 -
Ca(NO3)2.3H2O 51 - - 104 -
Na(NO3)2.6H2O 53 - - 158 -
40
Chapitre II. Matériaux à Changement de Phase – Paraffines, Sels ou Eutectiques
Tableau II-7 – Propriétés de quelques sels hydratés [LAN83] [ABH83] [FAR04] [ZAL03]
[SHA99b] [HAS98] [HAL71] [HIR94] (suite)
Nom Tfusion Masse Conductivité hLS Fusion
(°C) volumique thermique (kJ.kg-1)
(kJ.kg-3) (W.m-1.kg-1)
Zn(NO3)2.2H2O 55 - - 68 congruente
FeCl3.2H2O 56 - - 90 -
CO(NO3)2.6H2O 57 - - 115 -
Ni(NO3)2.6H2O 57 - - 168 -
MnCl2.4H2O 58 - - 151 -
CH3COONa.3H2O 58 - - 270 - 290 incongruente
LiC2H3O2.2H2O 58 - - 251-377 -
MgCl2.4H2O 58,0 - - 178 -
NaOH.H2O 58 - - 272 -
Na(CH3COO).3H2O 58 - - - -
Cd(NO3)2.4H2O 59 - - 98 -
Cd(NO3)2.1H2O 59,5 - - 107 -
Fe(NO3)2.6H2O 60 - - 125 -
NaAl(SO4)2.12H2O 61 - - 181 incongruente
FeSO4.7H2O 64 - - 200 -
Na3PO4.12H2O 65 - - 168 -
Na2B4O7.10H2O 68 - - - -
Na3PO4.12H20 69 - - - -
LiCH3COO.2H2O 70 - - 150 - 251 congruente
Na2P2O7.10H2O 70 - - 186 - 230 incongruente
Al(NO3)2.9H2O 72 - - 155 - 176 incongruente
Ba(OH)2.8H2O 78 1 93784°C, 0,653 , 0,67898C
86C
265 - 280 congruente
2 07024°C
Al2(SO4)3.18H2O 88 - - 218 incongruente
Sr(OH)2.8H2O 89 - - 370 incongruente
Mg(NO3)2.6H2O 89 - 90 1 55094°C, 0,49095C, 162 - 167 congruente
1 63625°C 0,502110C
KAl (SO4)2.12H2O 91 - - 184 -
(NH4)Al(SO4).6H2O 95 - - 269 -
Na2S.51/2H2O 97,5 - - - -
LiCl.H2O 99 - - 212 incongruente
CaBr2.4H20 110 - - - -
Al2(SO4)2.16H2O 112 - - - -
MgCl2. 6H2O 115 - 117 14 50120°C, 0,570120C, 165 - 169 -
1 44278°C 0,598140C
NaC2H3O2.3H2O 137 1 450 - 172 -
41
Chapitre II. Matériaux à Changement de Phase – Paraffines, Sels ou Eutectiques
changement de volume inférieur aux autres MCP. Cela rend plus facile la conception du
conteneur du MCP.
V.5 Eutectiques
Un MCP eutectique est composé d’au moins deux composants ; chacun de ces matériaux
fond et se solidifie congrûment, et lorsqu’ils sont mélangés ils forment des cristaux
homogènes lors de leur cristallisation conjointe ([LAN89]). Généralement les eutectiques
changent de phase sans ségrégation grâce à l’homogénéité des cristaux formés, ce qui
laisse peu d'opportunité à la séparation des composants. Lors de la fusion, les deux
composants se liquéfient en même temps, encore une fois avec une séparation peu
probable. Une liste des eutectiques développés est donnée tableau II-8.
42
Chapitre II. Matériaux à Changement de Phase – Paraffines, Sels ou Eutectiques
Tableau II-8 – Propriétés de quelques eutectiques organiques et inorganiques [LAN83]
[ABH83] [GAR85] [FAR04] [ZAL03] [SHA99b] [HAS98] [HAL71] (suite)
Nom Composition (% m) Tfusion(°C) hLS (kJ.kg-1)
Napthalene + Benzoic Acid 67.1 + 32,9 67 234
AlCl3 + NaCl + ZrCl2 79 + 17 + 4 68 95
AlCl3 + NaCl + KCl 66 + 20 + 14 70 209
NH2CONH2 + NH4Br 66.6 + 33,4 76 151
LiNO3 + NH4 NO3 + NaNO3 25 + 65 + 10 80,5 113
AlCl3 + NaCl + KCl 60 + 26 + 14 93 213
AlCl3 + NaCl 66 + 34 93 201
NaNO2 + NaNO3+KNO3 40 + 7 + 53 142 -
Tableau II-9 – Propriétés de quelques MCP solide - solide [LAN83] [ABH83] [GAR85] [ZAL03]
[SHA99b] [HAS98] [HAL71]
Nom Tfusion(°C) hLS (kJ.kg-1)
38,2 %NPG/61,8 %PE 26 - 32 18 - 75
38,2 %NPG/61,8 %TAM 22 - 35 27 - 33
76,4 %NPG/33,6 %TAM 28 - 38 75 - 80
76,4 %NPG/33,6 %PE 31 - 37 35 - 46
91 %NPG/9 %PE 31 - 36 68
91 %NPG/9 %TAM 30 - 39 143 - 150
Neopentyl Glycol(NPG) 43 130
Diamnopentacrythritol 68 184
2-Amino - 2 - Methyl - 1, 3 - Propanediol 78 264
2 - Methyl - 2 - Nitro - 1, 3 - Propanediol 79 201
Trumethylolethane 81 192
Pentaglycerin 81 192
2-Hydroxymethyl-2-Methyl-1, 3 - 81 192
Propanediol
Monoaminopentaerythritol 86 192
Polyéthylène réticulé 110 - 115 125 - 146
Tris(Hydroxymethyl)Acide Acétique 124 205
2-Amino-2-Hydroxymethyl-1, 3 - 131 285
Propanediol
Trihydroxy Methyl Aminométhane (TAM) 135 295
Polyéthylène réticulé HDPE 125 - 146 167 - 201
2,2-Bis(Hydroxymethyl) Acide propionique 152 289
38.2%NPG/61.8%PE 170 147
Penterythritol(PE) 185 303
43
Chapitre II. Matériaux à Changement de Phase – Paraffines, Sels ou Eutectiques
V.6.1 Avantages du polyéthylène réticulé
Bien qu’ils soient plus coûteux que les autres MCP, les polyéthylènes n'ont pas besoin
d'encapsulation. Par conséquent, ils offrent d’une part un coût acceptable pour le système de
stockage thermique complet, et d’autre part une mise en œuvre facilitée [LAN83]. Les
polyéthylènes réticulés sont généralement atoxiques et chimiquement inertes.
V.7 Polyalcools
Les polyalcools absorbent de l’énergie thermique lors du changement de structure interne
d'un état hétérogène (basse température) à une forme cubique à faces centrées (haute
température) ([WAN00]). Ils présentent de nombreux avantages par rapport aux MCP solide-
liquide, comme le faible changement de volume, pas de fuite et aucune dégradation de la
matière ([WAN00]). Ces avantages sont cependant associés à de forts inconvénients : faible
chaleur latente, haute température de changement de phase et un coût élevé. Wang et al.
([WAN00]) ont fait un mélange de polyalcools pour ajuster la température de changement de
phase. Certains de leurs résultats sont présentés dans le Tableau II-9 comme le néopentyl
glycol (NPG), le Penaerythritol (PE), et le Trihydroxy Methyl Aminométhane (TAM).
Les propriétés thermodynamiques privilégient les MCP ayant les chaleurs latentes les plus
élevées et étant les plus conductifs. En outre, un MCP avec un grand sous-refroidissement
est fortement déconseillé à moins d’intégrer un déclencheur de solidification dans le
système, ce qui le rend plus complexe. D’autre part, la fusion non-congruente est également
un problème qui conduit à éviter ce type de MCP non fiable. Comme il a été précisé
précédemment, le MCP doit être compatible avec le conteneur, d’où un matériau ayant un
effet corrosif ne peut pas être utilisé dans une structure métallique. Deux solutions peuvent
remédier plus ou moins à ce défaut :
44
Chapitre II. Matériaux à Changement de Phase – Paraffines, Sels ou Eutectiques
1. Encapsulation du matériau, c’est-à-dire ensacher les matériaux dans des sacs en
polymère, ce qui rend leur utilisation moins pratique en cas de macro-encapsulation
(>10 mm), et moins efficaces en cas de mini et micro-encapsulation (<10 mm).
2. Revêtement des surfaces métalliques par une couche isolante, ce qui est
techniquement complexe et imparfait.
Finalement le coût et la disponibilité peuvent aussi jouer un rôle pour le choix du MCP.
45
Chapitre II. Matériaux à Changement de Phase – Paraffines, Sels ou Eutectiques
Tableau II-11 - Propriétés de quelques MCP commercialisés [CRI] [TEA] [RUB] [CLI] [EPS] [MIC]
(suite)
TFusion Fournisseur Nom hLS (kJ.kg-1) Masse Conductivité Famille
(°C) volumique thermique
(kJ.kg-3) (W.m-1.kg-1)
27 Rubitherm RT 27 179 760-880 0,2 Organique
27 Mitsubishi STL 27 213 1090 - Sel Hydraté
28 Climator C28 162 1450 0,5-0,7 Sel Hydraté
28 EPS Ltd. A28 230 789 0,21 Organique
29 TEAP 29T 175 1500 1 Sel Inorganique
30 EPS Ltd. S30 201 1304 0,58 Sel Hydraté
31 Rubitherm RT 31 130 760-880 0,2 Organique
32 Climator C32 162 1450 0,5-0,7 Sel Hydraté
32 EPS Ltd. A32 215 845 0,21 Organique
32 EPS Ltd. S32 186 1460 0,51 Sel Hydraté
32 TEAP 32S 200 1450 0,6 Sel Inorganique
34 EPS Ltd. S34 115 2100 0,52 Sel Hydraté
34 TEAP 34S 220 1450 0,6 Sel Inorganique
36 Rubitherm RT 36 130 760-880 0.2 Organique
36 TEAP 36S 230 1450 0,6 Sel Inorganique
37 Microtek MPCM 37D 195 900 - Organique
39 EPS Ltd. A39 190 900 0,22 Organique
40 EPS Ltd. PlusICE X40 131 1046 0,209 Solide-Solide
41 Rubitherm PX 42 114 694 0,1 Organique
42 EPS Ltd. A42 195 905 0,21 Organique
42 Rubitherm RT 41 174 760-880 0,2 Organique
43 Microtek MPCM 43D 105 900 - Organique
44 EPS Ltd. S44 105 1584 0,43 Sel Hydraté
45 TEAP 45S 230 1450 0,6 Sel Inorganique
46 EPS Ltd. S46 190 1587 0,45 Sel Hydraté
47 Mitsubishi STL 47 221 1340 - Sel Hydraté
48 Climator C48 227 1360 0,5-0,7 Sel Hydraté
49 Rubitherm PX 50 110 694 0,1 Organique
50 EPS Ltd. S50 104 1601 0,43 Sel Hydraté
50 Rubitherm RT 49 168 760-880 0,2 Organique
52 Rubitherm RT 52 173 760-880 0,2 Organique
52 Mitsubishi STL 52 201 1300 - Sel Hydraté
53 EPS Ltd. A53 210 910 0,22 Organique
55 EPS Ltd. A55 195 905 0,22 Organique
55 Rubitherm RT 55 172 770-880 0,2 Organique
55 Mitsubishi STL 55 242 1290 - Sel Hydraté
58 Climator C58 288 1460 0,5-0,7 Sel Hydraté
58 EPS Ltd. A58 210 910 0,22 Organique
58 EPS Ltd. S58 167 1505 0,69 Sel Hydraté
58 Rubitherm RT 58 178 770-880 0,2 Organique
58 Rubitherm PX 58 117 694 0,1 Organique
60 EPS Ltd. A60 197 910 0,22 Organique
60 Rubitherm RT 60 144 770-880 0,2 Organique
62 EPS Ltd. A62 202 910 0,22 Organique
62 Rubitherm RT 61 146 780-880 0,2 Organique
65 Rubitherm RT 65 152 780-880 0,2 Organique
70 Climator C70 281 1700 0,6-0,7 Sel Hydraté
70 EPS Ltd. A70 220 890 0,23 Organique
72 EPS Ltd. S72 140 1666 0,58 Sel Hydraté
80 EPS Ltd. PlusICE X80 192 1193 0,335 Solide-Solide
46
Chapitre II. Matériaux à Changement de Phase – Paraffines, Sels ou Eutectiques
Tableau II-11 - Propriétés de quelques MCP commercialisés [CRI] [TEA] [RUB] [CLI] [EPS] [MIC]
(suite)
TFusion Fournisseur Nom hLS (kJ.kg-1) Masse Conductivité Famille
(°C) volumique thermique
(kJ.kg-3) (W.m-1.kg-1)
82 EPS Ltd. A82 170 850 0,22 Organique
82 Rubitherm RT 82 176 770-880 0,2 Organique
82 Rubitherm PX 82 116 694 0,1 Organique
83 EPS Ltd. S83 152 1600 0,62 Sel Hydraté
86 Mitsubishi PCM 86 246 - - Organique
89 EPS Ltd. S89 163 1550 0,67 Sel Hydraté
95 EPS Ltd. A95 202 900 0,22 Organique
100 Rubitherm RT 100 124 770-880 0,2 Organique
110 EPS Ltd. H100 243 2145 - Sel Inorganique
117 EPS Ltd. S117 169 1450 0,7 Sel Hydraté
118 Mitsubishi PCM 118 339,8 1480 - Organique
Des réservoirs de stockage de froid en forme de plaques planes ont été étudiés
numériquement et expérimentalement dans Ismail et al. [ISM99]. Des arrangements
tubulaires, notamment des tubes d'avion, tubes à ailettes axiales et radiales, tubes à ailettes
internes et externes ont également été étudiés comme dans Ismail [ISM98b] et al Ismail et
al. [ISM00]. Des réservoirs de stockage à froid du type glace sur serpentin ont été étudiés
par Jeckel et al. [JEC93], Ismail et al. [ISM70] et Benta [BEN01].
D'autres géométries dans lesquelles le MCP est encapsulé dans des conteneurs cylindriques
ou sphériques ont également été étudiées. Ismail et Henriquez [ISM02] présentent un
modèle numérique pour simuler un système de stockage composé de capsules sphériques
remplis de MCP (eau) placés à l'intérieur d'un réservoir cylindrique. Un modèle transitoire
simplifié unidimensionnel a été utilisé en divisant le réservoir en un nombre de couches dont
l'épaisseur axiale est toujours égale ou supérieure au diamètre d’une capsule. Tieyi et al.
[TIE96] ont étudié un système de stockage de chaleur latente, dans lequel le réservoir de
stockage de froid est rempli d'enceintes rectangulaires contenant du MCP. Ils ont présenté
un modèle numérique simplifié pour étudier les effets des facteurs opérationnels tels que la
puissance frigorifique, la surfusion du MCP et ont comparé leurs résultats avec les résultats
expérimentaux disponibles.
47
Chapitre II. Matériaux à Changement de Phase – Paraffines, Sels ou Eutectiques
Kurklu et al. [KUR96] ont présenté un modèle mathématique pour prédire le rendement
thermique d'un réservoir de stockage contenant de polypropylène comme MCP dans un tube
de 1 m de long et 38 mm de diamètre échangeant de la chaleur avec de l’air. Le modèle est
basé sur l’équation de conservation de l’énergie. Les résultats des simulations ont été
cohérents avec les résultats expérimentaux. Velraj et al. [VEL97] ont réalisé une étude
numérique et expérimentale sur la solidification d’un MCP à l'intérieur de tubes verticaux
ailetés intérieurement et ont proposé un modèle bidimensionnel basé sur la formulation
enthalpique de l’équation de l’énergie. Le modèle a été utilisé pour étudier l'influence du
diamètre du tube interne et du nombre d'ailettes sur l’intensité du flux thermique superficiel
échangé en dynamique et sur la vitesse de solidification.
Tao [TAO67] présente une méthode pour l'analyse de la solidification dans des géométries
cylindriques et sphériques. Cho et Sunderland [CHO70] utilisent une transformation sous la
forme d'un produit de la température et du rayon pour obtenir un terme de correction à
appliquer à la géométrie plane pour représenter le cas d'une géométrie sphérique. Shih et
Chou [SHI71] ont proposé une méthode itérative d'approximations successives pour
représenter la solidification à l'intérieur d’une sphère. Leurs résultats sont comparables aux
résultats numériques disponibles, à l’exception des cas présentant des grandes valeurs du
nombre Stefan et des petites valeurs du nombre de Biot. Pedroso et Domoto [PED73] ont
appliqué la technique de perturbation pour étudier la solidification d’un MCP à l'intérieur
d’une sphère dont la température de surface est uniforme. Riley et al. [RIL74] ont étudié
analytiquement le processus de solidification à l'intérieur des cylindres et des sphères de
pour des petits nombres de Stefan. Moore et Bayazitoglu [MOO82] ont étudié le processus
de fusion du MCP à l’intérieur des enveloppes sphériques et ont été en mesure de prédire la
position de l'interface et les profils de température pour différentes valeurs du nombre de
Stefan et de Fourier. Hill et Kucera [HIL83] ont développé une méthode semi-analytique pour
étudier la solidification dans un conteneur sphérique, y compris les effets du rayonnement de
la surface du récipient. Ils ont comparé leurs prédictions de la position de l'interface avec des
résultats disponibles et ont obtenu un bon accord. Milanez [MIL84a] et Milanez et Ismail
[MIL84b] ont présenté un modèle pour la solidification à l'intérieur des sphères en utilisant
l’approche du front mobile pour déterminer la position du front mouvant et les profils de
températures. Les résultats des simulations ont été validés par des expériences. D'autres
études numériques et expérimentales sur la solidification et la fusion en géométrie sphérique
ont été réalisés par Prud’homme et al. [PRU89], Cho et Choi [CHO00] et Chen et al.
[CHE99].
Une approche analytique a été réalisée par Mosaffa et al. [MOS12] afin de suivre le profil de
température lors du changement de phase dans une géométrie rectangulaire 2D. Dans cette
même étude, un modèle numérique a aussi été développé et a donné des résultats
conformes au modèle analytique. El Qarnia et al. [ELQ12] et Ho et Chu [HO95] ont étudié un
système de protection thermique par des couches MCP/air rectangulaires verticales. Les
effets de convection dans la partie liquide ont été pris en compte. La durée de protection
assurée est exprimée en fonction du nombre de Rayleigh modifié, du nombre de Stefan
modifié, de la surfusion, de l’épaisseur de la couche MCP/air et du rapport des épaisseurs
MCP/air.
Une synthèse sur les travaux de modélisation du changement de phase est présentée dans
les références [YAO89] et [REG08] concernant les aspects théoriques et pratiques.
48
Chapitre II. Matériaux à Changement de Phase – Paraffines, Sels ou Eutectiques
Dans le cas de la fusion (Figure II-7 a), une première partie la chaleur absorbée de
l’extérieur par le matériau fait chauffer le liquide fondu, une deuxième partie (généralement la
plus grande) est responsable de faire fondre une masse du solide et la troisième partie est
responsable du chauffage du solide restant. Lors de ce processus, une gaine de liquide se
forme entre la surface active et le matériau encore solide. Les transferts s’y produisent par
conduction, puis dans de nombreux cas par convection naturelle et, éventuellement, par
rayonnement entre la surface et le front solide. La convection naturelle améliore le transfert
de chaleur, mais sa formulation mathématique en est complexe, car elle nécessite le
traitement du couplage entre les équations de Navier-Stokes dans le liquide et de Fourier
dans le solide.
Dans le cas de la solidification (Figure II-7 b), la chaleur cédée à l’extérieur provient de trois
phénomènes : refroidissement du solide, refroidissement du liquide non encore solidifié et
essentiellement de la solidification d’une masse de liquide. L’énergie est extraite à travers
une croûte solide qui s’épaissit progressivement et dont la résistance thermique augmente
avec l’épaisseur. Le liquide se refroidit progressivement pour atteindre la température de
solidification. Ces problèmes sont généralement traités en conduction pure, car le rôle de la
convection naturelle y est en général peu important. Cette courte description permet de
comprendre que ce problème de transferts couplés est fortement non linéaire. De ce fait, les
solutions analytiques sont rares et, bien qu’elles permettent de fournir certains ordres de
grandeur, elles sont souvent peu utilisables dans les configurations réelles. En particulier, la
résolution des cas multidimensionnels, où les conditions aux limites sont celles rencontrées
dans les applications pratiques, n’a été possible que par le développement des méthodes
numériques.
49
Chapitre II. Matériaux à Changement de Phase – Paraffines, Sels ou Eutectiques
(Eq. II-14)
(Eq. II-15)
Les conditions aux limites au niveau des surfaces frontières du domaine peuvent être soit
des températures imposées, soit des flux imposés, ou être décrites par un coefficient
d’échange à la paroi. Une analyse des lois d’échelle permet d’introduire le paramètre
pertinent des problèmes de changement de phase, le nombre de Stefan.
En considérant le cas de la fusion en géométrie plane dont la phase initiale est un solide à la
température de fusion TF avec une température de paroi TP imposée et en prenant
l’hypothèse que le transitoire thermique est négligeable à l’échelle du temps caractéristique
du processus de fusion, l’équation II-15 devient :
(Eq. II-16)
50
Chapitre II. Matériaux à Changement de Phase – Paraffines, Sels ou Eutectiques
La résolution de cette équation différentielle ordinaire, en faisant apparaître la diffusivité du
liquide
(Eq. II-17)
(Eq. II-18)
(Eq. II-19)
appelé nombre de Stefan est défini comme le rapport des enthalpies sensible et de
changement d’état mis en jeu.
Le Tableau II-12 donne quelques valeurs usuelles du nombre de Stefan pour des matériaux
courants.
Tableau II-12 – Valeurs du nombre de Stephan pour différents MCP [CLA81]
MCP Température de fusion (°C) Nombre de Stéphan
Glace 0 ~0,0125 ΔT
n-Octadécane 28 ~0,0091 ΔT
Etain 232 ~0,0038 ΔT
Soude 297 à 321 ~0,0055 ΔT
Aluminium 660 ~0,0027 ΔT
Chlorure de Sodium 801 ~0,0023 ΔT
51
Chapitre II. Matériaux à Changement de Phase – Paraffines, Sels ou Eutectiques
propriétés physiques du matériau comme la masse volumique, la conductivité thermique et
la chaleur massique dépendent de son état physique. Pour cela, il faut introduire une
variable décrivant l’état physique de chaque volume élémentaire : le titre de liquide, qui est le
rapport de la masse du liquide sur la masse totale.
(Eq. II-20)
(Eq. II-21)
Les propriétés physiques dépendront de cette variable et sont égales aux moyennes
pondérées de celles en états solide et liquide [NAK95].
(Eq. II-22)
Une simplification de cette complexité consiste à exprimer l’équation de l’énergie sous forme
« enthalpique » (Equation II-23), du fait que la variable enthalpie reste monotone (strictement
croissante ou décroissante) durant le changement de phase. L’équation de la chaleur dans
les phases solide et liquide est :
(Eq. II-23)
(Eq. II-24)
(Eq. II-25)
52
Chapitre II. Matériaux à Changement de Phase – Paraffines, Sels ou Eutectiques
(Eq. II-26)
Les figures II-8 et II-9 montrent deux maillages des éléments 2D d’épaisseur e,
respectivement rectangulaire et cylindrique, ainsi que quelques échanges thermiques
internes et des conditions limites possibles.
T∞
Mi Hi Ti ki
Mi Hi Ti ki
q
T∞
Figure II-8 – Maillage d’un élément 2D Figure II-9 – Maillage d’un élément 2D
rectangulaire cylindrique
(Eq. II-27)
Les mailles des extrémités peuvent subir un flux thermique externe imposé et/ou avoir une
température de paroi imposée. L’effet de la convection dans la partie liquide est négligé,
étant donné que le MCP est mis en place sous une faible épaisseur dans un échangeur où le
liquide sera bloqué. Le transfert radiatif est également négligé à cause des écarts de
températures relativement faibles. Ces deux hypothèses seront validées dans le chapitre
suivant.
Ces équations sont appliquées sur les volumes élémentaires quelles que soient leurs
géométries. Plusieurs modèles de conduction lors du changement de phase ont été
développés sur Dymola (langage Modelica) pour différentes géométries. Ce langage permet
d’assembler les modèles élémentaires à d’autres modèles plus complexes pour obtenir la
géométrie réelle. Des exemples seront montrés dans le chapitre suivant.
53
Chapitre II. Matériaux à Changement de Phase – Paraffines, Sels ou Eutectiques
mesure calorimétrique donne le flux mesuré sur la paroi inférieure du conteneur du matériau
en fonction de la température de cette même paroi (Figure II-10). Or les MCP présentent une
diffusivité thermique très basse, la température mesurée est différente de la température
réelle du matériau qui est en plus non homogène. Pour cela le modèle de conduction
développé sera utilisé pour déterminer une calorimétrie plus précise utilisant une méthode
inverse.
Un modèle de la cellule de mesure comportant le MCP est conçu sur Dymola. Les
paramètres physiques du MCP seront balayés afin de retrouver les valeurs qui permettent de
retracer la même DSC donnée par la machine. La conductivité thermique est connue
puisqu’elle a été mesurée en régime permanent monophasique. La masse volumique,
mesurée par picométrie, est également connue. Les paramètres à trouver sont : la
température de déclenchement du changement de phase, l’hystérésis à la cristallisation, le
glissement de température, les capacités calorifiques solide et liquide et la chaleur latente.
Comme décrit précédemment, le retard à la cristallisation est un effet aléatoire, pourtant
cette analyse permettra de donner un ordre de grandeur de l’hystérésis. Afin de réduire les
degrés de liberté des paramètres à trouver, un bilan global d’énergie échangée est appliqué
entre les températures initiale et finale. Cette équation permet de lier les capacités
calorifiques monophasiques et la chaleur latente du matériau. L’organigramme de la Figure
II-11 montre le fonctionnement en boucle inverse, retrouvant les bons paramètres du MCP.
Modèle cellule + échantillon
Qp = f(Tp)
Initiation Propriétés MCP
54
Chapitre II. Matériaux à Changement de Phase – Paraffines, Sels ou Eutectiques
La DSC de ces matériaux a été réalisée par l’équipe MMC de EDF dans le cadre du projet
Stock-aiR2. Le cycle thermique effectué sur les échantillons pour les mesures d’enthalpie
sera (Température de paroi imposée):
Départ à 25 °C
Montée rapide de 25 °C à 70 °C à 1 °C.mn-1
Stabilisation à 70 °C durant 600 s
Stabilisation à 70 °C durant 1 800 s
Descente de 70 °C à 0 °C à 0 ,15 °C.mn-1 -Cycle refroidissement
Stabilisation à 0 °C durant 1 800 s
Stabilisation à 0 °C durant 900 s
Stabilisation à 0 °C durant 1 800 s
Montée de 0 °C à 70 °C à 0,15 °C.mn-1 - Cycle chauffage
Stabilisation à 70 °C durant 1 800 s
Retour rapide à 25 °C à 1 °C.mn-1
L’évolution de la température est présentée Figure II-12.
80
70
60
50
T ( C)
40
30
20
10
0
0 5 10 15 20
t (h)
Figure II-12 – Cycle thermique effectué sur les échantillons pour la DSC
40
30
20
10
0
0 100 200 300 400
H(kJ.kg-1)
Figure II-13 – Courbes mesurées de solidification dans le plan (T,H)
55
Chapitre II. Matériaux à Changement de Phase – Paraffines, Sels ou Eutectiques
50 50
T ( C)
T ( C)
40 40
30 30
20 20
10 10
0 0
0 100 200 300 400 0 100 200 300 400
H(kJ.kg-1) H(kJ.kg-1)
Figure II-14 – Comparaison entre DSC mesurée Figure II-15 – Comparaison entre DSC mesurée
et calculée lors de la fusion de la RT 36 et calculée lors de la solidification de la RT 36
Une première méthode simplifiée, qui permet d’obtenir une DSC précise, est de calculer les
chaleurs massiques en monophasique dans un premier temps. De préférence, la chaleur
massique liquide est calculée au début de la phase de refroidissement de la DSC et la
chaleur massique solide au début de la phase de chauffage de la DSC. En suite, le point de
début de fusion est identifié, la chaleur latente est calculée en faisant le bilan énergétique
entre l’état initial et l’état final. En traçant les deux droites représentant le chauffage liquide et
solide, le glissement de température sera évalué. L’hystérésis est évaluée par la courbe de
refroidissement, où la courbe montre le début de la solidification. Il est préférable de l’évaluer
en tant que valeur variable dans une plage obtenue par plusieurs essais.
Une deuxième méthode permettant d’obtenir une DSC précise consiste à chauffer ou à
refroidir, à une vitesse assez lente telle que la variation de la température appliquée sur le
plateau puisse atteindre tout le MCP. Cette vitesse dépendra de la diffusivité du matériau.
Cette méthode risque de donner une hystérésis plus grande puisqu’elle dépend de la vitesse
de refroidissement.
56
Chapitre II. Matériaux à Changement de Phase – Paraffines, Sels ou Eutectiques
40
30
20
10
0
0 100 200 300 400
H(kJ.kg-1)
Figure II-16 – Courbes calculées de solidification dans le plan (T,H)
Toutes ces courbes montrent un glissement de température plus ou moins grand lors de la
solidification. Ceci est dû à plusieurs effets, dont les deux plus importants sont :
les impuretés existantes dans le MCP ; la paraffine n’est pas purifiée à 100 %
les paraffines commerciales sont généralement un mélange de paraffines donnant la
température de fusion désirée.
D’autre part ces quatre paraffines présentent une surfusion assez faible qui ne dépasse pas
1 K pour les échantillons testés (de l’ordre du cm3).
VIII. Conclusions
Les systèmes de stockage par chaleur latente offrent une densité d'énergie thermique de
stockage élevée, par unité de masse et par unité de volume, à une température constante
correspondant à la température de transition de phase du MCP.
57
Chapitre II. Matériaux à Changement de Phase – Paraffines, Sels ou Eutectiques
Le stockage thermique par chaleur latente peut être obtenu par des transformations solide-
solide, solide-liquide, solide-gaz et liquide-gaz. Le changement de phase adapté à
l’effacement d’une pompe à chaleur est celui de la transition solide-liquide grâce à la quantité
de chaleur délivrée et avec changement de volume acceptable lors de cette transformation.
Les MCP peuvent être classés en trois grandes familles selon leurs natures chimiques :
Organiques, inorganiques et présentant un point d’eutexie. Les propriétés de deux MCP
d’une même famille sont semblables.
La calorimétrie différentielle à balayage avec ajustement des propriétés par méthode inverse
permet d’évaluer la chaleur latente, les températures de fusion et de solidification et les
glissements de température lors de la transformation avec plus grande précision. Elle permet
aussi d’évaluer l’hystérésis du matériau.
Dans cette thèse une méthode numérique par volume finis a été développée, basée sur la
forme enthalpique de l’équation de l’énergie. Le modèle de changement de phase développé
est connecté à d’autres modèles pour former le système global afin de concevoir l’échangeur
stockeur et le système PAC.
58
Chapitre III. Conception d’un Echangeur Stockeur - Déstockeur
I. Introduction
Le condenseur est l’élément clé de la pompe à chaleur intégrant un stockage thermique. Cet
échangeur est le siège de l’échange entre les trois « matériaux actifs » du système : le fluide
frigorigène, l’air ambiant et le matériau de stockage. Dès lors, la faisabilité et la rentabilité de
ce système sont liées directement à la performance de cet échangeur stockeur. Il sera conçu
principalement pour être capable de stocker la quantité d’énergie thermique nécessaire
pour l’effacement, la délivrer avec la puissance requise à un niveau de température
assurant le confort thermique et un écart de températures acceptable entre stockage et
déstockage. Il devra respecter également d’autres critères comme un faible débit de
soufflage d’air, la plus grande compacité possible et un fonctionnement sans hystérésis.
Selon les calculs des besoins de la maison « Mozart Passive », un stock de 800 Wh
d’énergie thermique est suffisant pour assurer le confort entre 18 h et 20 h durant toute
l’année. La stratégie choisie est celle du dimensionnement maximal, sans inconfort et une
durée d’effacement d’au moins deux heures. Le système doit être capable de décharger une
puissance thermique maximale de 400 W. Le système de chauffage qui va être utilisé est
une PAC monosplit à une unité intérieure (Figure III-1).
MCP
59
Chapitre III. Conception d’un Echangeur Stockeur - Déstockeur
que les ailettes dans le MCP pourraient améliorer largement le transfert de chaleur du
processus de fusion de l'unité de stockage thermique.
Plusieurs brevets présentent des dispositifs d’échangeurs tubes plats ailetés, souvent utilisés
pour la climatisation des véhicules, permettant un stockage par chaleur latente. La
Figure III-2 montre le schéma du brevet [EP1221389A2] qui propose d’insérer des réservoirs
de MCP entre le tube où circule le fluide caloporteur et les ailettes externes qui sont brasées
sur les réservoirs. Le brevet [FR2878613A1] décrit une méthode d’intégration plus simple au
niveau de l’assemblage de l’échangeur. Les réservoirs sont adjacents aux tubes des
caloporteurs, comme le montre la Figure III-3, les ailettes externes sont brasées en partie sur
les tubes et en partie sur les réservoirs. Dans le brevet [EP1424531B1], les tubes sont
composés de faisceaux de circulation du caloporteur adjacents à d’autres faisceaux
comportant du MCP. Le dispositif du brevet [DE1020004035818] présente des petits
réservoirs rectangulaires de MCP brasés d’un côté et de l’autre de chaque tube d’un
échangeur à tubes plats ailetés.
Sections de passage
du caloporteur
Ailette
Réservoirs MCP
L’inconvénient de ces dispositifs est le ratio volumique associé à chacun des trois fluides :
caloporteur, fluide externe et fluide de stockage qui est inadapté. Les caloporteur et l’air sont
en écoulement, le volume associé est le volume de passage. Par contre, le volume associé
au MCP permet un stockage thermique très limité. En d’autres termes la densité énergétique
entre les trois composants de l’échangeur est mal équilibrée.
Le brevet [WO2006059005] propose des réservoirs parallélépipédiques de MCP ayant une
petite épaisseur. Ces derniers peuvent être insérés dans des tubes de circulation d’un
caloporteur dans un échangeur ou bien entre les plaques d’un échangeur à plaques. De
même les brevets [DE102006007543], [DE3910356] et [DE10240246] présentent des
réservoirs de stockage par chaleur latente où le MCP occupe le circuit du premier fluide d’un
échangeur à plaques et le caloporteur circule dans le deuxième. Le réservoir du brevet
[DE19953113] est identique sauf que le caloporteur circule dans des tubes perpendiculaires
aux plaques.
Le brevet [DE10256665] porte sur un système de climatisation doté d’un réservoir de
stockage par chaleur latente. Le froid est généré par un circuit frigorifique et stocké dans un
réservoir de MCP où un échangeur intermédiaire peut absorber ce froid pour le transmettre à
un autre échangeur et refroidir de l’air injecté dans l’habitacle d’un véhicule. Les brevets
60
Chapitre III. Conception d’un Echangeur Stockeur - Déstockeur
[DE10233415] et [DE10108152] présentent également des systèmes de climatisation
comportant un réservoir de stockage thermique.
Le déstockage par l’intermédiaire d’une boucle supplémentaire de caloporteur entraîne des
écarts de températures supplémentaires entre source et « puits » d’où des pertes
exergétiques, ce qui affecte le rendement du système.
Une étude menée par M. Medrano et al. [MED09] compare les performances de différents
échangeurs commercialisés intégrant des MCP. La meilleure performance est fournie par
celui où le dopage du côté MCP est le plus efficace.
ρ λ (Eq. III-1)
61
Chapitre III. Conception d’un Echangeur Stockeur - Déstockeur
ρ λ (Eq. III-2)
62
Chapitre III. Conception d’un Echangeur Stockeur - Déstockeur
Pour assembler les modèles de conduction, ces derniers doivent comporter des connecteurs
qui communiquent la température et le flux thermique échangé aux extrémités. D’autres
géométries peuvent être modélisées aussi selon la structure de l’échangeur. Ces modèles
peuvent aussi être utilisés pour les métaux qui sont des MCP avec une très haute
température de fusion.
(Eq. III-3)
63
Chapitre III. Conception d’un Echangeur Stockeur - Déstockeur
III.2.3 Assemblage des modèles
Les deux modèles élémentaires principaux sont utilisés pour composer le modèle principal.
D’autres modèles secondaires sont aussi nécessaires, selon la structure, pour compléter le
modèle. La Figure III-6 montre un exemple d’assemblage représentant un simple échange
entre un volume de MCP contenu dans un tube métallique et de l’air circulant sur la paroi
externe de ce tube.
Puits air
Source air
MCP
Air
64
Chapitre III. Conception d’un Echangeur Stockeur - Déstockeur
Ce condenseur stockeur va être étudié en prenant en compte trois modes de
fonctionnement :
Constitution du stock d’énergie latente : le système de compression fonctionne, la
température de condensation du fluide frigorigène est supérieure à la température de
fusion du MCP qui passe donc en phase liquide. Dans ce cas, le ventilateur du
condenseur est arrêté, la chaleur n’est pas diffusée dans la pièce, toute l’énergie de
condensation du fluide frigorigène sert à la fusion du MCP. Ce mode est appelé mode
de « charge ». En cas de besoin thermique faible, l’échangeur peut fonctionner en
mode de charge combiné, où le stockage se fait avec un chauffage simultané de l’air.
Restitution de l’énergie stockée : le système de compression est arrêté, le fluide
frigorigène ne circule donc pas, le MCP est entièrement fondu, le ventilateur du
condenseur fonctionne, l’air entrant est à une température inférieure à la température
de solidification, le MCP restitue sa chaleur à l’air en se solidifiant, appelé ici mode
« décharge ».
Fonctionnement en condensation du fluide frigorigène avec MCP fondu : c’est un
mode où le changement de phase du MCP n’est pas utilisé, la masse de MCP rend
plus constante la température de surface de l’échangeur. C’est un mode secondaire
intéressant à analyser du fait de la constance de la température de soufflage et donc
du confort ressenti.
65
Chapitre III. Conception d’un Echangeur Stockeur - Déstockeur
Les paramètres de ce sous-modèle sont :
L Masse volumique du fluide en état liquide
V Masse volumique du fluide en état vapeur
CpL Capacité calorifique du fluide en état liquide
CpV Capacité calorifique du fluide en état vapeur
λL Conductivité thermique du fluide en état liquide
λV Conductivité thermique du fluide en état vapeur
Tev Température d’évaporation (considérée égale à celle de condensation)
hev Chaleur latente massique du fluide
Qimp Puissance source donnée par le fluide (lors du stockage)
hin Coefficient d’échange convectif avec le tube interne
hex Coefficient d’échange convectif avec le tube externe
Rex Rayon externe de la partie
Rin Rayon interne de la partie
θ Angle complémentaire de celui du secteur
Nr Nombre de volumes suivant la direction radiale
Ntt Nombre de volumes suivant la direction tangentielle
66
Chapitre III. Conception d’un Echangeur Stockeur - Déstockeur
IV.1.2 Déstockage tube dans tube
La phase de déstockage (décharge) est donc simulée en premier lieu, car elle est la phase la
plus délicate du processus d’échange. Pour obtenir les résultats de déstockage, il suffit
d’imposer des conditions de température ou de flux à la surface externe du grand tube. On
rappelle que le système de compression est à l’arrêt, il n’y a donc pas de débit circulant de
fluide frigorigène. Dans le condenseur, le fluide frigorigène est considéré comme étant à
l’état gazeux, la phase liquide étant stockée dans un réservoir situé à l’aval du condenseur ;
des simulations sont faites aussi en considérant qu’il reste du liquide.
Une étude de sensibilité sera faite sur les paramètres suivants :
a) l’état du fluide frigorigène
b) le maillage du MCP
c) la conductivité thermique du MCP
d) le coefficient d’échange entre le tube externe et l’air
e) les diamètres des deux tubes
Les hypothèses considérées sont :
l’effet de dilatation du MCP est très faible
l’effet de la convection dans le MCP est négligé
la masse volumique du fluide frigorigène est de 1 300 kg.m-3 en liquide et 6 kg.m-3 en
vapeur, celle du MCP est de 760 kg.m-3
la conductivité thermique du fluide frigorigène est de 90 mW.m-1.K-1 en liquide et de
10 mW.m-1.K-1 en vapeur, celles du MCP sont égales à 0,2 W.m-1.K-1 en liquide et en
solide, sauf indication contraire
les capacités calorifiques du fluide frigorigène et de l’air sont égales à
1 kJ.kg-1, celle du MCP est considérée égale à 2 kJ.kg-1
la chaleur latente du MCP est considérée égale à 130 kJ.kg-1
la température de solidification du MCP est considérée de 36 °C
le coefficient d’échange de convection entre le fluide frigorigène stagnant et les tubes
est considéré égal à 20 W.m-2.K-1 en état gazeux et 1 000 W.m-2.k-1 en état liquide
le coefficient d’échange équivalent entre l’air et le tube externe rapporté à la surface
de la paroi, est considéré par défaut égal à 20 W.m-2.k-1 et la surface ailetée est
considérée 5 fois plus grande que celle du tube seul. L’épaisseur des tubes est
toujours considérée égale à 0,35 mm
la résistance de contact au brasage est négligeable
l’angle que forme la surface de brasage est de 75°
la température d’entrée d’air ambiant est de 21 °C
le débit d’air soufflé par unité de longueur du tube est de 1 g.s-1.m-1
le MCP est initialement à une température uniforme de (Tf +1)
Le diamètre externe du tube interne sera noté Dint, et celui du tube externe sera noté Dext. La
première simulation est celle du petit tube où stagne le fluide frigorigène. Ce petit tube est
brasé sur la génératrice du grand tube aileté qui contient le MCP. Le fluide frigorigène est en
état gazeux, ses mouvements éventuels sont dus à des effets de convection naturelle
puisqu’il n’y a pas de débit circulant. Le diamètre du tube interne (Dint) est de 4 mm, celui du
tube externe (Dext) est de 8 mm. Le MCP est discrétisé en cinq zones suivant la direction
radiale et cinq zones suivant la direction angulaire, soit au total 25 éléments. Ce cas est
considéré comme le cas de référence pour la décharge thermique.
67
Chapitre III. Conception d’un Echangeur Stockeur - Déstockeur
40
T1
38
T2
36
34 T3
32 T4
T ( C)
30 T5
28 Tint
26 Text
24
22
20
0 100 200 300 400
t (s)
Figure III-10 – Températures du MCP et des tubes externe et interne dans le cas de référence
40 20
38 18
Text
36 16
34 14
32 Tsort 12
T ( C)
P (W.m-1)
30 10
28 8
Tent
26 6
24 4
22 2
20 0
0 100 200 300 400 0 100 200 300 400
t (s) t (s)
Figure III-11 – Températures du tube externe,
Figure III-12 – Flux thermique échangé avec l’air
d’entrée et de sortie de l’air dans le cas de
dans le cas de référence
référence
La Figure III-10 montre les températures du MCP dans la direction radiale (T1 au centre, T5
à la paroi), celle des parois des tubes interne et externe. Les emplacements des
températures mesurées sont montrés Figure III-13.
68
Chapitre III. Conception d’un Echangeur Stockeur - Déstockeur
La Figure III-11 montre la température du tube externe ainsi que celles d’entrée et de sortie
de l’air après l’échange avec le tube. La Figure III-12 montre la puissance absorbée par l’air
par mètre linéaire de tube. Le MCP étant initialement en état liquide, se refroidit et entre
rapidement en état diphasique liquide-solide pour céder la chaleur latente stockée. Au fur et
à mesure que le front solide se déplace de l’extrémité vers le centre du petit tube, la
résistance thermique entre l’air et le MCP diphasique augmente. Ceci se traduit par la baisse
de la température de sortie de l’air, ainsi que par la baisse du flux thermique allant du tube
vers l’air (Figure III-12). Le point critique du MCP, c’est-à-dire le dernier volume de MCP à se
solidifier, est celui au centre du tube interne. Il se solidifie après 120 secondes. Après
400 secondes, l’échangeur entier atteint la température d’entrée de l’air, 21 °C. L’écart de
températures entre le point critique et la température de sortie d’air est dû principalement
aux deux résistances thermiques : celle du MCP qui augmente avec l’épaisseur de la couche
du MCP solide et celle associée à la convection entre l’air et le tube externe.
69
Chapitre III. Conception d’un Echangeur Stockeur - Déstockeur
20 20
Ptot Ptot
18 18
Pbras Pbras
16 16
Pgaz Pliq
14 14
12 12
P (W.m-1)
P (W.m-1)
10 10
8 8
6 6
4 4
2 2
0 0
0 100 200 300 400 0 100 200 300 400
t (s) t (s)
Figure III-14 – Flux thermiques passant à travers Figure III-15 – Flux thermiques passant à travers
le métal et à travers le fluide en état gazeux le métal et à travers le fluide en état liquide
Dans ces conditions d’absence de débit de frigorigène, le flux passant du fluide frigorigène
directement vers le MCP est négligeable, ce qui est confirmé par les simulations sur le
modèle de référence en supprimant le sous-modèle représentant le fluide. Les Figure III-16
et Figure III-17 représentant les mêmes courbes que celles des Figure III-11 et Figure III-12
du cas de référence, sont presque identiques. Dans les simulations ultérieures du mode
déstockage, le sous-modèle du fluide frigorigène sera supprimé.
40 20
T1
38 18
T2
36 16
T3
34 14
T4
32 12
T5
T ( C)
P (W.m-1)
30 Tint 10
28 Text 8
26 Tsort 6
24 4
22 2
20 0
0 100 200 300 400 0 100 200 300 400
t (s) t (s)
Figure III-16 – Températures du MCP, des tubes
Figure III-17 – Flux thermique échangé avec l’air
interne et externe et de sortie de l’air en
en supprimant le modèle du fluide frigorigène
supprimant le modèle du fluide frigorigène
stagnant
stagnant
70
Chapitre III. Conception d’un Echangeur Stockeur - Déstockeur
40 20
T1 Nr5
38 18
T2 Nr10
36 16
T3
34 14
T4
32 12
T5
T ( C)
P (W.m-1)
30 10
Tint
28 8
Text
26 6
Tsort
24 4
22 2
20 0
0 100 200 300 400 0 100 200 300 400
t (s) t (s)
Figure III-18 – Températures du MCP, des tubes Figure III-19 – Flux thermique échangé avec l’air
interne et externe et de sortie de l’air dans le cas dans les cas de maillage de référence et du
du maillage raffiné maillage raffiné
P (W.m-1)
Tsort_1,0
30 10
28 T1_2,0 8
26 Tsort_2,0 6
24 4
22 2
20 0
0 100 200 300 400 0 100 200 300 400
t (s) t (s)
Figure III-20 – Températures du MCP au centre du Figure III-21 – Flux thermiques échangés avec
tube, et celles de sorties de l’air pour différentes l’air pour différentes valeurs de conductivité
valeurs de conductivité thermique du MCP thermique du MCP
71
Chapitre III. Conception d’un Echangeur Stockeur - Déstockeur
IV.1.6 L’effet du coefficient d’échange externe
Le transfert convectif entre le tube externe et l’air est le deuxième facteur limitant du transfert
thermique global. La sensibilité de cet effet est étudiée en faisant varier le coefficient
d’échange équivalent entre le tube externe et l’air par rapport au cas de référence. La
Figure III-22 montre l’écart de températures entre le tube externe et la sortie d’air lors de la
phase de décharge des cas où le coefficient d’échange externe est de 10 W.m-2.K-1, 20 W.m-
2 -1
.K et 30 W.m-2.K-1. Cet écart devient acceptable (pour cette géométrie) au-delà de 20 W.m-
1 -1
.K (moins de 2 K). L’échange global lors de la solidification est d’autant plus intense que le
coefficient de convection entre le tube externe et l’air est élevé (Figure III-23).
6 20
H20 18 H20
5
16
H10 H10
14
4
H30 12 H30
ΔT ( C)
P (W.m-1)
3 10
8
2
6
4
1
2
0 0
0 100 200 300 400 0 100 200 300 400
t (s) t (s)
Figure III-22 – Ecarts de température entre le Figure III-23 – Flux thermique échangé avec l’air
tube externe et la sortie de l’air pour différentes pour différentes valeurs de coefficient d’échange
valeurs de coefficient d’échange externe externe
72
Chapitre III. Conception d’un Echangeur Stockeur - Déstockeur
40
38
36 T1_D4
34 Tsort_D4
32
T1_D5
T ( C)
30
28 Tsort_D5
26 T1_D7
24
Tsort_D7
22
20
0 100 200 300 400
t (s)
Figure III-24 – Températures du MCP au centre du tube et celles de sorties de l’air
pour différentes valeurs de Dint
40
38
36 T1_D8
34 Tsort_D8
32
T1_D10
T ( C)
30
28 Tsort_D10
26 T1_D12
24
Tsort_D12
22
20
0 100 200 300 400
t (s)
Figure III-25 – Températures du MCP au centre du tube, et celles de sorties de l’air pour différentes
valeurs de Dext
73
Chapitre III. Conception d’un Echangeur Stockeur - Déstockeur
latente de 200 kJ.kg-1, l’énergie qui peut être stockée dans ce tube sera de 5,65 kJ. Pour
stocker 800 Wh, l’échangeur doit être composé d’environ 512 m linéaires de tube, ce qui
n’est pas réaliste. D’autre part, la section de passage du fluide frigorigène est assez grande
pour les débits nécessaires. Ceci risque de dégrader le coefficient d’échange entre le fluide
et les parois des tubes. Ces faits, mènent à aborder la deuxième option d’intégrer le MCP
entre les deux tubes et faire circuler le fluide frigorigène dans le petit tube.
Le modèle a été refait en inversant les places du fluide frigorigène et du MCP. Le sous-
modèle du fluide frigorigène sera toujours éliminé pour la phase de décharge. La phase de
charge sera également simulée. Afin de rendre l’échangeur plus compact, une première
configuration sera de prendre un diamètre de tube externe de 8 mm et celui du tube interne
de 4 mm.
40 16
38 Tcrit
14
36
Tsort 12
34
32 10
T ( C)
P (W.m-1)
30 8
28 6
26
4
24
22 2
20 0
0 200 400 600 800 0 200 400 600 800
t (s) t (s)
Figure III-26 – Température du point critique et
Figure III-27 – Flux thermique absorbé par l’air
celle de sortie d’air pour le cas où le MCP est
pour le cas où le MCP est dans l’inter-tubes
dans l’inter-tubes
Les Figure III-26 et Figure III-27 montrent les résultats des simulations du mode décharge.
La Figure III-26 montre la température du point critique du MCP, et celle de sortie d’air. Le
temps de déstockage est de 390 s, et la température de sortie d’air est de 31 °C à la fin de la
phase de décharge. La Figure III-27 montre la puissance absorbée par l’air, elle décroît de
14 à 9 W par mètre linéaire de tube durant la décharge. Ces valeurs sont proches de celles
du cas de référence sauf que le volume de MCP dans cette disposition est trois fois plus
grand que celui du cas de référence. En comparant un mètre linéaire d’un tube avec cette
disposition avec trois mètres de tube de référence, l’énergie de stockage est identique. Les
mêmes courbes sont représentées pour le cas d’un débit de soufflage trois fois plus grand et
un coefficient d’échange 2,5 fois plus grand. La Figure III-28 montre un temps de décharge
de 175 s et une température de sortie d’air de 28 °C à la fin de la décharge. La Figure III-29
montre une diminution de la puissance de 35 à 22 W durant la décharge. Cette puissance
est à peu près 1,5 fois plus petite que celle de trois mètres de tubes de référence. Cette
disposition permet de gagner en compacité avec une perte de puissance. Cette perte doit
être récupérée en dopant le MCP, ce qui peut être modélisé en augmentant sa conductivité
thermique.
74
Chapitre III. Conception d’un Echangeur Stockeur - Déstockeur
40 45
38 Tcrit 40
36 35
Tsort
34
30
32
25
T ( C)
P (W.m-1)
30
20
28
15
26
24 10
22 5
20 0
0 200 400 600 800 0 200 400 600 800
t (s) t (s)
Figure III-28 – Température du point critique et Figure III-29 – Flux thermique absorbé par l’air
celle de sortie d’air pour le cas où le MCP est pour le cas où le MCP est dans l’inter-tubes à
dans l’inter-tubes à débit d’air trois fois plus grand débit d’air trois fois plus grand
46
44
42
Tcrit
40
Tcond
T ( C)
38
36
34
32
30
0 50 100 150 200
t (s)
Figure III-30 – Température critique du MCP et celle de condensation du fluide frigorigène en mode
de charge
Le fait marquant des simulations des deux phases est que le point critique est le même dans
les deux cas. La température est donc homogène dans toute la structure métallique quelle
que soit la direction du flux thermique global absorbé par l’air ou cédé par le fluide
frigorigène. Ceci peut aussi valider le fait que l’effet dominant du transfert thermique global
est celui de la conduction dans le MCP.
75
Chapitre III. Conception d’un Echangeur Stockeur - Déstockeur
IV.1.1 Dimensionnement de l’échangeur tube dans tube
En considérant un tube de Dint de 4 mm et de Dext de 8 mm, et le MCP inséré entre les deux
tubes, le volume de MCP que peut contenir un mètre d’un tel tube est de 40 ml. Si le MCP a
une masse volumique de 1 000 kg.m-3 et une chaleur latente de 200 kJ.kg-1, l’énergie qui
peut être stockée dans ce tube sera de 8 kJ. Pour stocker 800 Wh, l’échangeur doit être
composé d’environ 360 m linéaires de tube. Afin d’obtenir une taille réaliste, il faut par
exemple concevoir un tube externe de diamètre de 30 mm, où la longueur de tubes
nécessaire pour stocker les 800 Wh sera d’environ 20 m linéaires. Cette conception a été
réalisée sur une première maquette.
76
Chapitre III. Conception d’un Echangeur Stockeur - Déstockeur
Figure III-31 – Dimensions en mm du tube dans tube avec des ailettes internes et externes
(Eq. III-4)
(Eq. III-5)
(Eq. III-6)
77
Chapitre III. Conception d’un Echangeur Stockeur - Déstockeur
IV.3.1.1 Déstockage à débit constant
45 1200
Tcrit
1000
40
Ttube
800
35
Tsort
T ( C)
P (W)
600
30
400
25
200
20 0
0 1000 2000 3000 4000 5000 6000 0 1000 2000 3000 4000 5000 6000
t (s) t (s)
Figure III-32 – Température du MCP au point
ème
critique, celle des tubes du 3 rang et celle de Figure III-33 – Puissance de déstockage pour un
3 -1
sortie d’air pour un débit de soufflage de débit de soufflage de 150 m .h avec la RT31
3 -1
150 m .h avec la RT31
45 1200
Tcrit
1000
40
Ttube
800
35
Tsort
T ( C)
P (W)
600
30
400
25
200
20 0
0 1000 2000 3000 4000 5000 6000 0 2000 4000 6000
t (s) t (s)
Figure III-34 – Température du MCP au point
ème
critique, celles des tubes du 3 rang et celle de Figure III-35 – Puissance de déstockage pour un
3 -1
sortie d’air pour un débit de soufflage de débit de soufflage de 150 m .h avec la RT36
3 -1
150 m .h avec la RT36
78
Chapitre III. Conception d’un Echangeur Stockeur - Déstockeur
La vitesse de l’air sur les ailettes est limitée à 1 m.s-1, ce qui limite le débit de soufflage à
750 m3.h-1. Le deuxième régime de déstockage est celui du déstockage à puissance
maximale, où l’air est soufflé avec un débit de 750 m3.h-1. Les mêmes types de courbes sont
représentés Figure III-36, Figure III-37, Figure III-38 et Figure III-39 pour les paraffines RT31
et RT36. Pour la paraffine RT31 le déstockage se déroule en 2 100 secondes. Les courbes
ont les mêmes allures que celles du cas précédent. La puissance de déstockage chute de
3 000 W à 300 W. Pour la paraffine RT36, la puissance chute de 5 000 W à 500 W en 1 200
secondes.
45 3000
Tcrit
2500
40
Ttube
2000
35
Tsort
T ( C)
P (W)
1500
30
1000
25
500
20 0
0 1000 2000 3000 4000 5000 6000 0 2000 4000 6000
t (s) t (s)
Figure III-36 – Température du MCP au point
ème
critique, celles des tubes du 3 rang et celle de Figure III-37 – Puissance de déstockage pour un
3 -1
sortie d’air pour un débit de soufflage de débit de soufflage de 750 m .h avec la RT31
3 -1
750 m .h avec la RT31
45 5000
Tcrit 4500
40 4000
Ttube
3500
35 3000
Tsort
T ( C)
P (W)
2500
30 2000
1500
25 1000
500
20 0
0 1000 2000 3000 4000 5000 6000 0 2000 4000 6000
t (s) t (s)
Figure III-38 – Température du MCP au point
ème
critique, celles des tubes du 3 rang et celle de Figure III-39 – Puissance de déstockage pour un
3 -1
sortie d’air pour un débit de soufflage de débit de soufflage de 750 m .h avec la RT36
3 -1
750 m .h avec la RT36
Le choix de « dopage » est valorisé par des simulations sur l’échangeur sans les ailettes
internes. La Figure III-40 montre la température du point critique qui est au centre du
troisième tube, la température de ce même tube et la température de sortie d’air lors du
déstockage pour un échangeur sans dopage interne et en utilisant la paraffine RT 36. Le
temps de déstockage est d’environ 4 000 s, contre 1 200 s pour le cas avec dopage. L’air
sort de l’échangeur entre 33 °C et 21,7 °C. La Figure III-41 montre la puissance de
déstockage pour ce cas. La puissance de déstockage chute de 3 000 W à 180 W. Cette
simulation montre la nécessité du dopage pour assurer le besoin lors du déstockage avec
une compacité acceptable.
79
Chapitre III. Conception d’un Echangeur Stockeur - Déstockeur
45 3000
Tcrit
Ttube 2500
40
Tsort 2000
35
T ( C)
P (W.m-1)
1500
30
1000
25
500
20 0
0 1000 2000 3000 4000 5000 6000 0 2000 4000 6000
t (s) t (s)
Figure III-40 – Température du MCP au point
ème Figure III-41 – Puissance de déstockage pour un
critique, celles des tubes du 3 rang et celle de 3 -1
débit de soufflage de 750 m .h avec la RT36
sortie d’air pour un débit de soufflage de
3 -1 sans les ailettes internes
750 m .h avec la RT36 sans les ailettes internes
900
800
700
600
500
Débit (m3.h-1)
400
300
200
100
0
0 1000 2000 3000 4000 5000 6000 7000 8000
t (s)
Figure III-42 – Débit d’air nécessaire pour un déstockage régulé pour assurer un besoin de 500 W
Il est donc essentiel de disposer d’un débit d’air variable pour adapter la puissance
thermique en fonction des besoins cibles.
80
Chapitre III. Conception d’un Echangeur Stockeur - Déstockeur
plaques ailetées (Figure III-44). Le MCP est intégré dans l’inter-volume entre les tubes et les
plaques. Un modèle dynamique de cet échangeur a été réalisé également sous Dymola.
Figure III-43 – Schéma représentant une partie Figure III-44 – Schéma représentant une partie
d’un échangeur tubes entre plaques (sans d’un échangeur tubes entre plaques (avec
ailettes) ailettes)
Figure III-45 – Coupe horizontale des tubes avec Figure III-46 – Modèle simplifié de la coupe
les plaques horizontale des tubes avec les plaques
Le modèle complet est l’agrégation des motifs élémentaires. La Figure III-47 montre
l’assemblage des motifs élémentaires en attribuant une couleur à chaque matière : bleue
pour le MCP, verte pour le fluide frigorigène en état gazeux, bronze pour le cuivre et grise
pour l’aluminium. La Figure III-48 montre la section élémentaire de l’échangeur pour les
simulations, limitée par deux tubes et deux plaques. Cette section élémentaire est elle aussi
un motif reproductible par symétrie permettant de représenter des échangeurs de cette
structure et de taille variable.
81
Chapitre III. Conception d’un Echangeur Stockeur - Déstockeur
IV.4.2 Déstockage tubes entre plaques
Il est clair que le point critique du MCP (celui qui se solidifie le dernier) est dans la zone du
point A (Figure III-48). L’écart de températures entre la température de solidification (ou
solidification) et celle de sortie de l’air est notée ΔTfa.
Une étude de sensibilité sera faite sur les paramètres suivants :
a) l’espacement entre deux tubes
b) la conductivité thermique du MCP
c) le coefficient d’échange entre les plaques et l’air
d) le débit de soufflage d’air
Les hypothèses considérées sont :
l’effet de dilatation du MCP est très faible
l’effet de la convection dans le MCP est négligé
le fluide frigorigène est considéré en état gazeux et stagnant
la masse volumique du fluide frigorigène est de 6 kg.m-3 en vapeur, celle du MCP est
de 760 kg.m-3
la conductivité thermique du fluide frigorigène est de 10 mW.m-1.K-1 en vapeur, celles
du MCP sont égales à 0,2 W.m-1.K-1 en liquide et en solide, sauf indication contraire
les capacités calorifiques du fluide frigorigène, de l’air sont égales à 1 kJ.kg-1, celle du
MCP est considérée égale à 2 kJ.kg-1
la chaleur latente du MCP est égale à 130 kJ.kg-1
la température de solidification du MCP est de 36 °C
le coefficient d’échange de convection entre le fluide frigorigène et les tubes est
considéré par défaut égal 20 W.m-2.K-1
le coefficient d’échange équivalent entre l’air et les plaques externes, rapporté à la
surface de la paroi, est considéré par défaut égal à 10 W.m-2.k-1 et la surface ailetée
est 10 fois plus grande que celle de la plaque nue
l’épaisseur des tubes est toujours considérée égale à 0,35 mm
l’épaisseur des plaques est considérée égale à 0,5 mm
la résistance de contact au brasage est négligeable
l’angle que forme la surface de brasage est de 75°
la température d’entrée d’air ambiant est de 21 °C
le débit d’air soufflé par unité de longueur est considéré par défaut de 1 g.s-1.m-1
le MCP est initialement à une température uniforme de (Tf +1)
Dans un premier cas, considéré comme la référence, le diamètre interne du tube est égal à
8 mm et l’espacement entre les centres des tubes est de 20 mm. La Figure III-49 montre la
température du MCP au point critique, celle d’entrée et de sortie d’air et celle des plaques.
Le MCP étant initialement en état liquide, se refroidit et entre rapidement en état diphasique
pour céder la chaleur latente stockée. Au fur et à mesure que la partie solide se développe,
la résistance thermique entre l’air et le MCP diphasique augmente. Ceci se traduit par
l’augmentation du ΔTfa, ainsi que par la baisse du flux thermique absorbé par l’air
(Figure III-50). L’écart de températures entre les plaques et l’air en sortie de l’échangeur est
négligeable. Le MCP est totalement solidifié après environ 1 100 secondes. Après
2 500 secondes, l’échangeur entier atteint la température d’entrée de l’air, 21 °C.
82
Chapitre III. Conception d’un Echangeur Stockeur - Déstockeur
40 20
38 Tparoi 18
36
Tent 16
Tcrit
34 14
Tsort
32 12
T ( C)
P (W.m-1)
30 10
28 8
26 6
24 4
22 2
20 0
0 500 1000 1500 2000 2500 3000 0 500 1000 1500 2000 2500 3000
t (s) t (s)
Figure III-49 – Température du MCP au point
Figure III-50 – Flux thermique absorbé par l’air
critique, celle des plaques et celles d’entrée et de
dans le cas de référence
sortie d’air dans le cas de référence
40 20
38 Tparoi 18
36
Tent 16
Tcrit
34 14
Tsort
32 12
T ( C)
P (W.m-1)
30 10
28 8
26 6
24 4
22 2
20 0
0 500 1000 1500 2000 2500 3000 0 500 1000 1500 2000 2500 3000
t (s) t (s)
Figure III-51 – Température du MCP au point
critique, celles des plaques et celles d’entrée et Figure III-52 – Flux thermique absorbé par l’air où
de sortie d’air où la distance inter-centres est de la distance inter-centres est de 15 mm
15 mm
Pour une distance de 35 mm entre les centres des tubes, la quantité d’énergie déstockée est
encore plus grande donc la durée de déstockage est plus longue (Figure III-53). La
puissance de décharge pour ce cas est supérieure à celles des deux autres cas
(Figure III-54), toujours pour la raison des surfaces d’échange entre le point critique et l’air.
83
Chapitre III. Conception d’un Echangeur Stockeur - Déstockeur
40 20
38 Tparoi 18
36
Tent 16
Tcrit
34 14
Tsort
32 12
T ( C)
P (W.m-1)
30 10
28 8
26 6
24 4
22 2
20 0
0 500 1000 1500 2000 2500 3000 0 500 1000 1500 2000 2500 3000
t (s) t (s)
Figure III-53 – Température du MCP au point
critique, celles des plaques et celles d’entrée et Figure III-54 – Flux thermique absorbé par l’air où
de sortie d’air où la distance inter-centres est de la distance inter-centres est de 35 mm
35 mm
84
Chapitre III. Conception d’un Echangeur Stockeur - Déstockeur
40
38 Cond 0,2
36
Cond 0,5
34
32 Cond 1,0
T ( C)
30
Cond 2,0
28
26
24
22
20
0 500 1000 1500 2000 2500 3000
t (s)
Figure III-56 – Température de sortie de l’air pour des différentes valeurs de conductivité thermique
du MCP
40 20
38 Tparoi 18
36
Tent 16
Tcrit
34 14
Tsort
32 12
T ( C)
P (W.m-1)
30 10
28 8
26 6
24 4
22 2
20 0
0 500 1000 1500 2000 2500 3000 0 500 1000 1500 2000 2500 3000
t (s) t (s)
Figure III-57 – Température du MCP au point
Figure III-58 – Flux thermique absorbé par l’air
critique, celles des plaques et celles d’entrée et
avec un coefficient d’échange externe
de sortie d’air avec un coefficient d’échange -2 -1
-2 -1 de 50 W.m .K
externe de 50 W.m .K
85
Chapitre III. Conception d’un Echangeur Stockeur - Déstockeur
30 40
Débit 1,0 38 Débit 1,0
25
Débit 1,5 36 Débit 1,5
34
20 Débit 2,0 Débit 2,0
32
P (W.m-1)
T ( C)
15 30
28
10
26
24
5
22
0 20
0 500 1000 1500 2000 2500 3000 0 500 1000 1500 2000 2500 3000
t (s) t (s)
Figure III-59 – Flux thermique absorbé par l’air Figure III-60 – Température de sortie d’air pour
pour différentes valeurs de débit d’air soufflé différentes valeurs de débit d’air soufflé
55
53
51
49 Tcrit
47
Tcond
T ( C)
45
43
41
39
37
35
0 200 400 600
t (s)
Figure III-61 – Température critique du MCP et celle de condensation du fluide frigorigène en mode
de charge du cas de référence de tubes entre plaques
86
Chapitre III. Conception d’un Echangeur Stockeur - Déstockeur
Le deuxième critère est celui du niveau de bruit : la vitesse de l’air dans l’échangeur ne doit
pas dépasser 1 m.s-1 pour limiter le niveau de bruit. Pour minimiser le volume du passage
libre de l’air dans l’échangeur, il faut conserver une vitesse d’air constante sur tout
l’échangeur. Ceci amène à choisir des passages à section constante, en considérant que la
variation de la masse volumique de l’air est négligeable entre l’entrée et la sortie.
Le troisième critère est celui de la température de sortie d’air : l’air doit sortir de l’échangeur
avec une température minimale de 26 °C, ce qui amène à un débit de soufflage maximal de
250 m3.h-1. En mettant ce débit en rapport avec la vitesse de 1 m.s-1, la section de passage
d’air sera de 0,07 m². En considérant une longueur de passage d’air égale à L pair, le volume
associé de ce compartiment sera de (70 x Lpair) litres.
Le quatrième critère est celui des tubes du fluide frigorigène. Ces tubes seront en contact
avec la paroi des conteneurs qui échange avec l’air.
Le cinquième critère est celui de la disponibilité des conteneurs et des ailettes métalliques.
Figure III-62 – Mouse d’aluminium avec différentes formes, de gauche à droite : cellules fermées peu
poreuses [MBG], semi fermées grandes pores [ALU], ouvertes structurées [CTIF], ouvertes non
structurées [FRA]
87
Chapitre III. Conception d’un Echangeur Stockeur - Déstockeur
l’isolation des ondes électromagnétiques spécialement dans les zones urbaines ou pour les
salles informatiques.
Dans le domaine de l’architecture, la structure poreuse et les nombreux grains donnent un
effet esthétique. La faible masse permet également un maniement facile des meubles et
autres objets de décoration similaires, même s’ils sont de grande taille.
Dans le domaine thermique, cette mousse est utilisée comme isolant thermique à haute
température et comme conducteur thermique dans un milieu peu conductif. En effet, la
mousse d’aluminium admet une conductivité thermique qui peut varier, selon ses propriétés,
entre 3 et 15 W.K-1.m-1pour une porosité entre 90 et 95 % [MES05] [MEI09]. Cette valeur est
considérée basse comparativement à celles des métaux, pourtant cette mousse admet une
haute température de fusion (~750 °C), pour cela elle peut servir à l’isolation thermique à
haute température. Ces valeurs de conductivité sont considérées élevées par rapport à
celles des MCP (de l’ordre de 0,2 W.K-1.m-1 pour les paraffines), et la mousse peut servir à
favoriser les échanges thermiques dans ces matériaux.
88
Chapitre III. Conception d’un Echangeur Stockeur - Déstockeur
La Figure III-64 montre une vue de dessus et en perspective d’un conteneur. Les ailettes
internes (côté MCP) sont de la mousse d’aluminium insérée dans les conteneurs. Les
conteneurs de MCP sont des tubes rectangulaires en laiton de section 50 x 100 mm² de
400 mm de longueur et 2 mm d’épaisseur de parois. L’échangeur est composé de
10 conteneurs. La mousse est fabriquée par la société canadienne Alusion [ALU]. Les
mousses fabriquées pour les utilisations mécaniques et thermiques sont de 10 à 100 fois
plus coûteuses. La mousse choisie présente une porosité d’environ 90 % avec une taille
moyenne de cellule de l’ordre de 10 mm. La paraffine utilisée comme indiqué ci-dessus, est
la Microtek MPCM 37 qui dispose des propriétés les plus adaptées pour cette application
(niveau de température et chaleur latente). L’échangeur contient 12 kg de cette paraffine. Le
circuit frigorifique est composé de tubes en cuivre standard de section ¼ pouce. Deux tubes
sont brasés sur toute la longueur des conteneurs. L’air circule parallèlement aux côtés de
100 mm des conteneurs. Sur ces côtés sont collées des ailettes en aluminium. La surface
d’échange externe totale est de l’ordre de 13 m². Le circuit frigorifique de l’échangeur est
composé de 20 tubes verticaux (deux par conteneur), connectés en 5 circuits parallèles avec
quatre tubes par circuit. Le choix du nombre de circuits a été fait de sorte à limiter les pertes
de pression du côté fluide en évaporation à 10 kPa pour un débit de 5 g.s-1 à une
température d’évaporation de 45 °C. De nombreuses études ont été réalisées sur les pertes
de pression en diphasique. Celles de Friedel (Eq. III-11) et Zhang et al. (Eq. III-12) [BIG01]
sont parmi les plus utilisées.
(Eq. III-7)
89
Chapitre III. Conception d’un Echangeur Stockeur - Déstockeur
(Eq. III-8)
Circuit frigorifique
MCP +
mousse
d’aluminium
Ailettes air
Air
90
Chapitre III. Conception d’un Echangeur Stockeur - Déstockeur
(Eq. III-9)
Ce coefficient est considéré constant sur toute la longueur de passage d’air, et égal à sa
valeur en x = 50 mm.
L’épaisseur de la colle entre les ailettes et les conteneurs est considérée de 0,5 mm, et sa
conductivité de 1 W.m-1.K-1.
(Eq. III-10)
Les simulations montrent que cet échangeur peut transférer 67 W pendant 12 heures. La
Figure III-65 montre les températures : du point critique du MCP, d’entrée et sortie d’air, de la
paroi du conteneur et de la paroi interne et externe de la colle. Après 12 heures, l’air sort de
l’échangeur avec une température de 26,8 °C. Le point critique se solidifie après 14h16. Les
courbes des parois et de la colle sont presque confondues, ce qui montre le faible effet de la
résistance de contact due à la colle. La Figure III-66 montre la variation du débit d’air au
cours du temps pour maintenir une puissance constante délivrée de 67 W. Pour l’échangeur
réel, le débit sera environ six fois ce débit calculé. Le débit de soufflage après 12 heures est
de 34 m3.h-1, d’où, dans le cas réel, ce débit devra être d’environ 200 m3.h-1après 2 heures
de déstockage.
40
Tsort 140
38
36 Tcrit 120
34
Tent 100
32
T ( C)
Q (m3.h-1)
Tparoi 80
30
28 Tco_ext 60
26 Tco_int 40
24
20
22
20 0
0 20000 40000 60000 80000 0 10000 20000 30000 40000 50000
t (s) t (s)
Figure III-66 - Débit de soufflage nécessaire pour
Figure III-65 - Températures lors du déstockage
assurer les 67 W lors du déstockage du prototype
du prototype sans prendre en compte l’effet de la
sans prendre en compte l’effet de la mousse
mousse d’aluminium
d’aluminium
91
Chapitre III. Conception d’un Echangeur Stockeur - Déstockeur
40 700
38 Tsort
600
36 Tcrit
34 500
Tent
32
400
T ( C)
Q (m3.h-1)
Tparoi
30
Tco_ext 300
28
26 Tco_int 200
24
100
22
20 0
0 5000 10000 0 2000 4000 6000 8000 10000
t (s) t (s)
Figure III-68 - Débit de soufflage nécessaire pour
Figure III-67 - Températures lors du déstockage
assurer les 400 W lors du déstockage du
du prototype en considérant une conductivité
prototype en considérant une conductivité
équivalente
équivalente
V. Conclusions
92
Chapitre III. Conception d’un Echangeur Stockeur - Déstockeur
Deux échangeurs stockeurs ont donc été conçus et présentés dans ce chapitre. Les
performances de ces échangeurs seront validées expérimentalement et les résultats
présentés dans le chapitre suivant.
93
Chapitre III. Conception d’un Echangeur Stockeur - Déstockeur
94
Chapitre IV. Validation Expérimentale de la PAC à Stockage
I. Introduction
Pour tester ce système de PAC intégrant un stockage par chaleur latente, un banc d’essais a
été conçu. La PAC avec son condenseur stockeur/déstockeur sera étudiée dans des
conditions réalistes en mesurant les grandeurs pertinentes pendant son fonctionnement. Le
banc permet de valider les résultats des simulations, de mieux comprendre les phénomènes
physiques et les problèmes pratiques de mise en œuvre.
Le banc est conçu de façon à tester plusieurs échangeurs stockeurs et chacun permet de
tester plusieurs MCP. De plus, l’échangeur sera testé en deux étapes : la première étape
consiste à faire fondre le MCP par de l’air chaud, seule l’étape de déstockage par l’air est
alors étudiée. Ceci a pour but de simplifier les problèmes d’essais et de se concentrer sur
l’étude de l’échangeur seul. La deuxième étape consiste à étudier le système PAC complet
en phases de stockage, de déstockage et en fonctionnement normal.
95
Chapitre IV. Validation Expérimentale de la PAC à Stockage
96
Chapitre IV. Validation Expérimentale de la PAC à Stockage
En ce qui concerne la boucle d’air, des sondes PT100 sont intégrées en amont et en aval de
chaque composant. La mesure de température au niveau de l’échangeur stockeur nécessite
une plus grande précision. Pour cela, cinq sondes PT100 ont été disposées sur les surfaces
de chaque côté de cet élément. Deux hygromètres à capteur de température intégré ont été
placés également à l’entrée et à la sortie du condenseur. Le débit volumique de l’air circulant
est mesuré à l’aide d’un anémomètre à fil chaud. Afin d’assurer un régime établi au niveau
de la mesure du débit, une gaine de longueur égale à cinq fois le diamètre est installée en
amont de l’anémomètre et la longueur en aval est égale à trois fois le diamètre. La régulation
de la température d’entrée est faite par contrôle du débit d’eau passant par le refroidisseur et
par la puissance fournie par le réchauffeur électrique. Le débit d’air circulant est contrôlable
grâce au ventilateur à vitesse variable et au registre motorisé. La puissance de l’échange
entre l’air et l’échangeur est calculée par l’équation :
(Eq. IV-1)
Les points 1 et 2 sont définis à partir des températures et des pressions à l’entrée et à la
sortie du compresseur respectivement. Le point 3 est défini à partir de la température du
fluide (supposée identique) à la sortie des deux condenseurs. En considérant une détente
adiabatique, le point 4 sera défini par la basse pression et l’enthalpie au point 3. La
puissance thermique absorbée par le condenseur stockeur (Pcond1) est calculée par
l’équation :
(Eq. IV-2)
97
Chapitre IV. Validation Expérimentale de la PAC à Stockage
Les deux détendeurs électroniques sont commandés par deux contrôleurs PID afin de
conserver une surchauffe et une répartition de débits entre les deux circuits. En effet, la
proportion des puissances absorbées par les condenseurs varie selon le régime : stockage,
chauffage avec stock et chauffage sans stock. De plus, la puissance absorbée par le
condenseur stockeur lors de stockage décroit du début à la fin de cette phase. Les
températures de condensation et d’évaporation sont régulées par les températures et les
débits d’eau glycolée entrant respectivement au condenseur secondaire et à l’évaporateur.
La température de l’eau glycolée peut être régulée à partir du groupe frigoifique. Le débit est
contrôlé par des vannes motorisées. Dans chacun des circuits d’eau sont placés : un débit
mètre électromagnétique (volumique), une sonde de température à l’entrée de l’échangeur et
une à la sortie. Les puissances absorbées par le condenseur secondaire et l’évaporateur
sont calculées à partir du bilan énergétique du côté de l’eau glycolée par l’équation suivante :
(Eq. IV-3)
(Eq. IV-4)
L’acquisition est réalisée par des modules Field Point. La Figure IV-4 montre le schéma de
l’armoire d’acquisition et de commande du banc.
FP 1601 ch2
ch3
ch4
T Sortie air condenseur 1 d
T Sortie air condenseur 1 e
T Entrée réchauffeur air
FP 1601 ch2
ch3
ch4
ch5 T Sortie réchauffeur air ch5
ch6 T Entrée air condenseur 1 a ch6
ch7 T Entrée air condenseur 1 b ch7
98
Chapitre IV. Validation Expérimentale de la PAC à Stockage
calcul. Refprop 8.0 est utilisé pour calculer les propriétés thermodynamiques du fluide
frigorigène.
99
Chapitre IV. Validation Expérimentale de la PAC à Stockage
cela, la détente est faite en deux étages : par le détendeur électronique, puis par une vanne
manuelle.
100
Chapitre IV. Validation Expérimentale de la PAC à Stockage
101
Chapitre IV. Validation Expérimentale de la PAC à Stockage
102
Chapitre IV. Validation Expérimentale de la PAC à Stockage
Figure IV-12 – L’échangeur complet avec le circuit frigorifique et les trémies d’air
0,9
0,8
0,7
163
Erécupérée (kWh)
0,6
255
0,5
0,4 370
0,3 460
0,2 575
0,1
790
0
0 1000 2000 3000 4000 5000 6000 7000 8000
t (s)
Afin de comparer les résultats expérimentaux aux simulations, un modèle a été réalisé sur
Dymola qui prend les valeurs de température et de débit d’air mesurées à l’entrée de
l’échangeur pour les injecter dans le modèle, puis comparer la température de sortie
calculée par le modèle à celle mesurée (Figure IV-16).
Tmodèle
T D
Air Echangeur
Tmesurée
Figure IV-16 – Modèle permettant de comparer les résultats expérimentaux à ceux du modèle dans
les mêmes conditions
Les courbes des Figure IV-17 et Figure IV-18 représentent les énergies restituées
expérimentales et calculées respectivement pour les cas de déstockage à 163 m3.h-1 et
255 m3.h-1. L’erreur relative maximale est de 2,5 %. Le modèle de changement de phase
utilisé dans cette première comparaison ne tient pas compte du glissement de température.
104
Chapitre IV. Validation Expérimentale de la PAC à Stockage
0,9 0,9
0,8 0,8
Erécupérée (kWh) 0,7 0,7
Erécupérée (kWh)
0,6 0,6
0,5 0,5
0,4 0,4
0,3 0,3
0,2 Réel 0,2 Réel
0,1 0,1
Mod Mod
0 0
0 2000 4000 6000 8000 0 1000 2000 3000 4000 5000 6000
t (s) t (s)
Figure IV-17 – Energies restituée expérimentale Figure IV-18 – Energies restituée expérimentale
3 -1 3 -1
et calculée pour déstockage à 163 m .h et calculée pour déstockage à 255 m .h
45 45
Tmod Tmod
40 Tent 40
Tent
35
Tcrit 35 Tcrit
T ( C)
T ( C)
Treelle Treelle
30 30
25 25
20 20
0 1000 2000 3000 4000 0 1000 2000 3000 4000
t (s) t (s)
Figure IV-19 – Température du MCP au point Figure IV-20 – Températures du MCP au point
critique et celles de l’air entrant et sortant, critique et celles de l’air entrant et sortant,
calculées et mesurées, de l’échangeur (modèle calculées et mesurées, de l’échangeur (modèle
sans glissement) avec glissement)
La Figure IV-21 montre la puissance et le débit d’air soufflé pour un déstockage régulé pour
assurer un besoin de 400 W. Il montre que le système a pu fonctionner pendant
6 800 secondes de déstockage à cette puissance (au lieu des 7 200 s escomptées), l’objectif
a donc été rempli à 95 %. Le fait de n’avoir pu utiliser que les deux tiers des tubes chargés
en MCP, a changé le ratio entre l’énergie latente et l’énergie sensible. Dès lors, la
température du stock chute plus rapidement comparativement à une masse de MCP
supérieure avec donc une plus grande capacité en chaleur latente. De là le constat d’une
température de sortie d’air plus basse et un débit de soufflage plus fort. D’autre part, le débit
maximal fourni par le ventilateur dans ce circuit est de 800 m3.h-1 d’où le temps de
déstockage a été réduit de 7 200 secondes à 6 800 secondes.
105
Chapitre IV. Validation Expérimentale de la PAC à Stockage
1200 900
P Débit
1000 750
Débit (m3/h)
800 600
P (W)
600 450
400 300
200 150
0 0
0 1000 2000 3000 4000 5000 6000 7000 8000
t (s)
Figure IV-21 – Déstockage à puissance constante régulée à 400 W
106
Chapitre IV. Validation Expérimentale de la PAC à Stockage
L’échangeur et la paraffine sont chauffés à environ 50 °C. La paraffine est versée
sous forme liquide dans les conteneurs. Un bilan de masse est fait avant et après
cette étape pour vérifier la masse du MCP versée (Figure IV-28).
L’échangeur est fermé de tous les côtés et prêt pour les tests (Figure IV-29).
107
Chapitre IV. Validation Expérimentale de la PAC à Stockage
108
Chapitre IV. Validation Expérimentale de la PAC à Stockage
(Eq. IV-5)
La Figure IV-30 montre les températures d’entrée et de sortie de l’air de l’échangeur ainsi
que TA et TB en régime établi avec une puissance échangée de 400 W. En se basant sur ces
résultats, la résistance de contact due à la colle sera de l’ordre de 0,002 K.W -1. Cette valeur
est quatre fois plus grande que la valeur estimée, et cela est dû à plusieurs facteurs comme :
l’imprécision au niveau de l’épaisseur de la couche de colle et de la valeur de sa
conductivité, la présence de bulles d’air emprisonnées à l’intérieur de la colle et la marge de
l’erreur des mesures au niveau des PT100 (de l’ordre de 0,2 K pour ces sondes). La
Figure IV-31 montre le coefficient d’échange externe expérimental et calculé en fonction du
débit de soufflage. La corrélation utilisée surestime légèrement ce coefficient.
40 6
38 TA TB Tsort Tent
5
36
hext (W.K-1.m-2)
34
4
32
T ( C)
30 3
28
2
26
24 hext_exp
1
22 hext_calc
20 0
0 500 1000 1500 2000 2500 100 120 140 160 180 200
t (s) Débit (m3.h-1)
Figure IV-31 - Comparaison entre hext
Figure IV-30 - Températures TA, TB, d’entrée et
expérimental et mesuré en fonction du débit de
de sortie d’air en régime établi
soufflage
109
Chapitre IV. Validation Expérimentale de la PAC à Stockage
pendant le déstockage. Après 7 200 secondes, l’air sort de l’échangeur à 25,2 °C, et le débit
de soufflage est de 310 m-3.h-1 (Figure IV-33). En répétant les simulations du modèle sans la
mousse d’aluminium, le facteur de réduction du temps de déstockage obtenu est de l’ordre
de 5, et la conductance du MCP à la fin de la décharge est de l’ordre 40 W.K-1.
40 900
38 TS
800
36 TE 700
34
600
Q (m-3.h-1)
32
500
T ( C)
30
400
28
300
26
24 200
22 100
20 0
0 2000 4000 6000 8000 0 2000 4000 6000 8000
t (s) t (s)
Figure IV-32 - Températures d’entrée et de sortie Figure IV-33 - Débit de soufflage d’air lors du
d’air lors du déstockage à 400 W déstockage à 400 W
110
Chapitre IV. Validation Expérimentale de la PAC à Stockage
condensation peut être déduite, c’est celle correspondant à la pression de saturation
à l’entrée du condenseur secondaire. La différence entre les deux températures de
condensation est due à la perte de pression créée par le débitmètre Coriolis.
- Puis le compresseur est arrêté, et l’air est soufflé à 21 °C afin d’absorber toute la
chaleur pendant les dernières 1 600 secondes de cette phase. L’échangeur est
entièrement au niveau de référence de 21 °C.
- La phase 2 commence quand le compresseur redémarre en absorbant environ
1 100 W et en délivrant une puissance variable entre 1 000 W et 1 200 W dans le
condenseur stockeur et entre 1 800 W et 2 800 W dans le condenseur secondaire.
Une partie de Pco1 est absorbée directement par l’air pour chauffer (400 W) et l’autre
partie est stockée. La température de condensation dans le condenseur stockeur
monte de 45 °C à 50 °C. Ces niveaux de températures de condensation sont assez
élevés du fait du surdimensionnement du compresseur. La puissance absorbée par
le condenseur stockeur est presque trois fois plus grande que sa puissance nominale
(400 W). La température d’évaporation varie de -5 °C à +5 °C. Cette phase dure à
peu près deux heures, et se termine en déclenchant la troisième phase.
- Durant la troisième phase, le détendeur du circuit stockeur est fermé et la puissance
de stockage moyenne en valeur absolue est égale à celle absorbée par l’air. La
puissance du côté fluide fluctue autour de la valeur zéro, du fait de la migration du
fluide dans le condenseur stockeur car il n’y a pas une vanne qui l’isole
complètement de la PAC. Cette phase est décrite dans le régime « déstockage à
400 W » du côté de l’échangeur stockeur. La puissance absorbée par le compresseur
varie entre 500 et 800 W, et celle du condenseur secondaire varie entre 1 000 et
2 500 W. La température de condensation varie entre 35 et 25 °C et celle
d’évaporation entre -5 et -15 °C.
- La quatrième phase commence quand le détendeur primaire est ré-ouvert et la
puissance Pco1 remonte, en assurant toujours le besoin de chauffage de 400 W,
mais sans stocker de l’énergie. La température de condensation varie entre 31 °C et
35 °C. Ce régime sera utilisé dans le fonctionnement de la PAC en mode normal,
c’est-à-dire. avant de stocker et après le déstockage. Durant cette phase, la
puissance au niveau du condenseur secondaire est de l’ordre de 2 000 W et la
puissance électrique absorbée par le compresseur est de 900 W.
- La cinquième phase commence quand le besoin de chauffage est de 900 W, et la
puissance Pco1 augmente pour assurer ce besoin. En augmentant la température de
condensation, le stock recommence à se former. La puissance de stockage monte de
200 à 1 000 W avec l’augmentation de la température de condensation de 35 °C à
45 °C. Cette phase décrit le comportement de la PAC durant le régime de chauffage
avec stockage. A la fin de cette phase l’échangeur stockeur est chargé.
111
Chapitre IV. Validation Expérimentale de la PAC à Stockage
Phase 1 2 3 4 5
2000
1500
Pair
1000
500
P (W)
0 Pstock
-500
-1000
Pco1
-1500
-2000
0 5000 10000 15000 20000 25000 30000
4000
Pev
3000
2000
Pco2
P (W)
1000
Pcomp
0
-1000
Pco1
-2000
0 5000 10000 15000 20000 25000 30000
65
55
Tcond1
45
35
T ( C)
25 Tcond2
15
5
Tev
-5
-15
0 5000 10000 15000 20000 25000 30000
t (s)
Figure IV-34 – Régime de fonctionnement
La Figure IV-35 présente les pertes de pression dans le condenseur, mesurées et calculées
par les corrélations de Fridel et de Zhang et al.. Les pertes de pression sont de l’ordre de
7 kPa pour une température de condensation de 45 °C et une puissance thermique de
1 100 W. Les deux corrélations surestiment les pertes de pression, mais celle de Friedel
reste la plus adaptée.
112
Chapitre IV. Validation Expérimentale de la PAC à Stockage
Phase 1 2 3 4 5
500
450
DP_MES
400
350
300
DP (mbar)
250 DP_FRI
200
150
100
DPZHA
50
0
0 5000 10000 15000 20000 25000 30000
t (s)
Figure IV-35 – Comparaison pour la perte de pression dans le condenseur
La Figure IV-36 montre le COP réel mesuré et le COP Carnot (entre température de
condensation et d’évaporation) de la PAC durant les huit heures de fonctionnement étudié.
Le COP thermodynamique est calculé par l’équation III-10, la consommation du ventilateur
n’est pas prise en compte. Ce graphe montre un rapport variant entre 0,45 et 0,50 entre le
COP réel et le COP Carnot.
Phase 1 2 3 4 5
10
8
COPcar
7
6
COP
3
COPmes
2
0
0 5000 10000 15000 20000 25000 30000
t (s)
Figure IV-36 – Comparaison entre COP thermodynamique de la PAC et COP Carnot
VII. Conclusions
Un banc d’essais a été construit afin de tester l’échangeur stockeur et la PAC complète.
Un échangeur de première génération, qui est une amélioration de la géométrie tube dans
tube, a été testé. Les résultats des simulations sont conformes aux résultats expérimentaux.
L’échangeur de « première génération » est capable de stocker la quantité d’énergie et de
décharger la puissance nécessaire. Par contre, il fonctionne à des débits de soufflages très
élevés et des températures de sorties relativement basses ; en plus il est volumineux.
113
Chapitre IV. Validation Expérimentale de la PAC à Stockage
Le prototype a été testé et validé sur le banc d’essais pour tous les régimes de
fonctionnement : charge, décharge et chauffage.
La PAC intégrant un stockage par chaleur latente est conçue. La prochaine étape est
d’étudier son comportement lors d’une saison complète et trouver la bonne stratégie de
contrôle en bénéficiant de l’avantage de l’inertie pour conclure sur son efficacité par rapport
à une PAC normale.
114
Chapitre V. Pompe A Chaleur avec Stockage
I. Introduction
En France, dans les années 90, la volonté de réduire les émissions de CO2 et de limiter
l'effet de serre a conduit à la recherche d’équipements de chauffage plus sobres
énergétiquement, ce qui a conduit au développement de pompes à chaleur résidentielles qui
répondent à ces attentes.
Les deux incitations fiscales successives de 2004 et 2005 ont encore accéléré le nombre de
pompes à chaleur (PAC) vendues en France pour le chauffage et/ou la production d’eau
chaude sanitaire. Cependant la RT 2012 et le coefficient de conversion d’énergie finale en
énergie primaire de 2,58 modifient la compétition entre les chaudières à gaz et les pompes à
chaleur.
Il existe trois principaux types de pompes à chaleur en fonction des technologies mises en
œuvre à l’extérieur et à l’intérieur des bâtiments :
eau/eau, influencé principalement par les différents niveaux de température des lois
d’eau
air/eau, influencé par les conditions climatiques et les lois d’eau
air/air, influencé par les conditions climatiques et la régulation de la température d’air
intérieur.
Pour notre étude, il est nécessaire de comprendre le fonctionnement réel des pompes à
chaleur air-air pour comprendre ce qu’une pompe à chaleur air/air dotée d’un système de
stockage au condenseur peut amener comme améliorations énergétiques ou, au contraire,
comme pénalité énergétique.
115
Chapitre V. Pompe A Chaleur avec Stockage
1200
1000
Puissance (W)
800
PuissPAC
600
Besoin
400
200
0
-10 -5 0 5 10 15
Text ( C)
Figure V-1 – Courbes théoriques de la puissance calorifique d’une PAC et du besoin d’un
bâtiment
La maison étudiée est une maison « Mozart passive » de type BBC. Le besoin est fortement
influencé par les apports internes et les flux solaires et pas uniquement donc par les
déperditions. La Figure V-2 montre les vrais besoins de la maison modélisée dans le
chapitre 1 durant toute l’année et celle de la puissance calorifique théorique d’une PAC à
vitesse constante en fonction de la température externe. Cette figure montre que les besoins
sont aléatoires et que même si la tendance globale de la baisse des charges avec la
température extérieure est bien vérifiée, la variation des besoins n’est pas uniquement
corrélée aux déperditions. En fait, il faut prendre en compte la surface du plan (Text, Besoins)
pour dimensionner le système de chauffage et sa régulation.
116
Chapitre V. Pompe A Chaleur avec Stockage
température de sortie du condenseur fixée à 35 °C montre une fluctuation entre 33 °C et
42 °C en régime cyclé. Les pertes énergétiques sont dues à :
une température moyenne de condensation supérieure au point de consigne du fait des
cycles courts pour les fortes surpuissances de la PAC vis-à-vis des besoins
une partie de la puissance est utilisée pour remettre les parois des échangeurs en
température et pour redistribuer le fluide frigorigène dans le circuit
lors des arrêts et en absence d’électrovanne, la pression tend à s’équilibrer dans
l’ensemble du circuit par migration du fluide frigorigène du condenseur vers
l’évaporateur. Le système doit compenser ce transfert de fluide à chaque démarrage
pour atteindre un régime stabilisé
la remise en pression stable des échangeurs s’effectue avec des variations initiales
fortes : pression de condensation plus élevée au démarrage et ce pendant plusieurs
minutes, pression d’évaporation plus basse ; ces conditions transitoires entraînent une
baisse du COP qui est d’autant plus marquée que le nombre de démarrages est grand
un appel important de courant à chaque démarrage équivalent à 5 à 7 fois le courant
nominal de fonctionnement. Ne durant que quelques millisecondes, cette puissance est
négligeable par rapport à la consommation totale, mais elle pénalise directement
l’installation électrique.
Ces effets dépendent directement de la fréquence de cyclage. Pour limiter les cycles de
fonctionnement trop courts (1 à 2 minutes de fonctionnement), la régulation de ces PAC
autorise généralement un dépassement de la consigne supérieure de la température de
surchauffe pour les températures d’air extérieur élevées, là où la PAC est surpuissante. On
doit donc noter que tous les cycles de fonctionnement dont la durée est inférieure à
5 minutes entraînent une baisse significative du COP réel due à l’ensemble des
irréversibilités transitoires présentées ci-dessus.
Une PAC air/air utilisant un compresseur à vitesse fixe, ayant une inertie thermique assez
faible, profite de l’inertie du bâtiment pour fonctionner en régime cyclé. Elle s’autorise une
fluctuation de la température du bâtiment. Cette fluctuation est mieux amortie tant que
l’inertie du bâtiment est plus importante. En effet, la surchauffe du bâtiment mène à des
déperditions plus importantes et donc à des consommations plus élevées. Le sous-
refroidissent est limité afin d’éviter l’inconfort dans le bâtiment. Pour cela les compresseurs à
vitesse fixe ne sont plus utilisés dans les PAC air/air. Ces dernières sont équipées de
compresseurs à vitesse variable.
117
Chapitre V. Pompe A Chaleur avec Stockage
Figure V-3 – Variation du COP en fonction de la puissance pour les machines à vitesse fixe
et à vitesse variable [PER04]
Ces résultats sont cohérents avec des travaux de SP (SP Technical Research Institute of
Sweden) [SP] où une campagne d’essais a été menée sur une multitude de pompes à
chaleur à vitesse variable, de marques différentes, pour étudier la variation du COP pour
différents taux de charge. Les résultats montrent une amélioration du COP variant entre 20
et 50 % pour un taux de charge qui descend jusqu’à 50 %.
Les performances des compresseurs à vitesse variable, pour les conditionneurs d’air utilisant
des compresseurs Mitsubishi, Hitachi et Sanyo de type piston rotatif, ont été analysées sur la
base de données expérimentales [SHA04] (Figure V-4). Cette analyse a permis de
développer un modèle pour calculer la performance pour des températures d’évaporation et
de condensation ainsi que des vitesses de rotation variables. L’évolution du COP en fonction
de la fréquence montre bien une chute significative du COP pour les faibles vitesses de
rotation. Cette chute est de 30 % pour un régime à fonctionnement minimal et de 23 % pour
le fonctionnement à puissance maximale.
En se basant sur les performances annoncés Figure V-4, la plage de besoins que peut
engendrer un compresseur à vitesse variable en fonctionnement continu est présentée
Figure V-5. La puissance electrique nominale de ce compresseur est de 110 W et un
rendement total de 0,5 pour cette puissance. Ce rendement baisse de 30 % lors du
fonctionnement à la puissance minimale (55 W électrique) et de 23 % à la puissance
maximale (220 W électrique). Cette dégradation est la seule pénalisation sur le COP lors du
118
Chapitre V. Pompe A Chaleur avec Stockage
fonctionnement entre les deux courbes Puiss Min, et Puiss Max (à part les cycles
givrage/dégivrage). Durant la saison de chauffage étudiée, cette PAC fonctionnera 32 % du
temps seulement en régime cyclé (besoin inférieur à la puissance calorifique minimale).
Figure V-5 – Plage de besoins que peut engendrer une PAC à compresseur à vitesse
variable
II.3 Givrage/dégivrage
Lorsque la température de surface de l’évaporateur est à la fois inférieure à celle de rosée
de l’air externe et inférieure à celle de cristallisation de l’eau (0 °C), une couche de givre se
forme au niveau de l'évaporateur des PAC à source air. Ce givre constitue une couche
isolante sur les parois de l’évaporateur et réduit les sections de passage d’air, ce qui réduit la
puissance des échanges et donc affecte négativement la performance de la PAC. Le
phénomène de givrage intervient là où les besoins de chauffage sont les plus importants (Text
basse), ce qui accentue son importance dans le bilan énergétique. La prise en compte de ce
phénomène dans la modélisation du comportement d'une PAC à air est donc essentielle.
Selon plusieurs études, le phénomène de givrage et les cycles de dégivrage associés
peuvent représenter jusqu'à 20 % de la consommation saisonnière d'une PAC [MIL82]
[MIL84] [KON95].
La croissance du givre dépend de la vitesse d’air ainsi que de la teneur de l’air en humidité.
Afin de limiter l’effet de la couche de givre, la régulation de la PAC procède à un dégivrage,
généralement par inversion du cycle. Une vanne à quatre voies sert à inverser le cycle de la
PAC, et une puissance calorifique sera injectée dans l’évaporateur afin de faire fondre le
givre. Cette procédure affecte aussi le COP, mais elle est nécessaire pour limiter l’effet du
givre qui devient de plus en plus pénalisant.
Un modèle, simplifiant les pertes par les cycles de givrage et dégivrage, est présenté par
[KIN09] (Figure V-6). Les cycles de givrage/dégivrage sont effectués généralement pour des
températures extérieures inférieures à 7 °C. La Figure V-7 montre la puissance calorifique
fournie par une PAC à vitesse fixe en fonction de la température extérieure en tenant compte
de l’effet de givrage. Cet effet oblige à sur-dimensionner la PAC afin de couvrir les besoins.
119
Chapitre V. Pompe A Chaleur avec Stockage
1200
1000
Puissance (W)
800
600
400
Théorique
200
Givrage
0
-7 -4 -1 2 5 8 11 14
Text ( C)
Figure V-6 – Pertes sur le COP dues aux cycles Figure V-7 – Puissance calorifique d’une PAC avec
givrage/dégivrage et sans effet de givrage
120
Chapitre V. Pompe A Chaleur avec Stockage
Zone 4 : décroissance continue du COP jusqu’à la valeur 1, où la résistance thermique
assure tout le besoin.
4 30
COP
3,5 20
COP Carnot
3 10
COP Carnot
COP mesuré
2,5 1 0
2
2 -10
3
1,5 -20
4
1 -30
-7 -6 -5 -4 -3 -2 -1 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14
Température extérieure (°C)
Figure V-8 – COP mesuré d’une PAC air/eau en fonction de la température extérieure
[DGUP]
De même, la Figure V-9 montre l’évolution du COP d’une PAC air/air à vitesse variable en
fonction de la température extérieure. Cette PAC est dimensionnée pour assurer les besoins
aux températures extérieures positives. Le fait marquant de cette courbe est la décroissance
du COP pour les températures extérieures élevées. Ceci est à cause de la faible inertie du
système et du faible rendement du compresseur à vitesse variable lors du fonctionnement à
puissance réduite.
Figure V-9 – COP mesuré d’une PAC air/air en fonction de la température extérieure
[DGALN]
121
Chapitre V. Pompe A Chaleur avec Stockage
100000 100
90000 90
80000 80
0 0
00:00
01:15
03:45
06:15
07:30
08:45
10:00
12:30
15:00
16:15
18:45
21:15
22:30
23:45
16:15
17:30
18:45
20:00
21:15
22:30
23:45
02:30
05:00
11:15
13:45
17:30
20:00
00:00
01:15
02:30
03:45
05:00
06:15
07:30
08:45
10:00
11:15
12:30
13:45
15:00
Figure V-10 – Production d’électricité en France Figure V-11 – Emissions de CO2 par kWh
durant une journée d’hiver [RTE] électrique produit durant un jour d’hiver
En outre, le stockage d’énergie durant les heures creuses et la restitution lors des heures de
pointe évite la consommation d’électricité provenant de la production par fioul. Ceci entraîne
bien sûr une réduction des émissions de CO2 de la production d’électricité. La Figure V-10
montre la production d’électricité par type de combustibles ou de renouvelables durant une
journée d’hiver. Les émissions de CO2 montrées Figure V-11 sont calculées sur la base des
données RTE France. Ces graphes montrent que la sur-demande entre 18h00 et 20h00 est
assurée par une production par fioul. La production de fioul, peut donc être limitée, voire
évitée, en aplatissant la courbe de consommation de 16h00 à 20h00. Selon RTE, 1 kWh
électrique coute 800 g de CO2 s’il est produit par fioul, 360 g par gaz et 0 g par hydraulique.
D’où, une PAC capable de stocker 0,8 kWh et de s’effacer entre 18h00 et 20h00 pourra
réduire les émissions jusqu’à 320 g de CO2 par jour. En considérant une valeur moyenne de
160 g sur chaque jour de la saison froide, cet effacement évitera environ 35 kg de CO2 émis
par an.
122
Chapitre V. Pompe A Chaleur avec Stockage
9%
8%
7%
(Eq. V-1)
Afin de calculer la perte totale, la température extérieure journalière a été moyennée sur
toute la saison de chauffage et présentée Figure V-14. En considérant par exemple un
stockage entre 16h00 et 18h00, suivi d’un déstockage de 18h00 à 20h00, la différence de
température moyenne entre ces deux phases sera de 1,5 °C. D’où la perte totale sera entre
1 % et 3 %. Plus le besoin est élevé, plus le temps de charge est long, plus le temps de
décharge est court. Ce cas est le plus pénalisant pour la PAC avec stockage, vu qu’elle
fonctionne plus longtemps pour charger, et d’autre part la différence de température est
minimale. En prenant le cas de besoin minimal, la PAC stockera l’énergie sur une durée
minimale, et déchargera durant une durée maximale. La température moyenne durant la
décharge sera plus faible, ce qui est gagnant pour la PAC à stockage.
123
Chapitre V. Pompe A Chaleur avec Stockage
6%
DT1 DT2
4% DT3 DT4
2%
-2%
-4%
-6%
-8%
-10 -5 0 5 10 15
Text ( C)
Figure V-13 – Pertes sur le COP pour plusieurs valeurs d’écart de températures lors de la
production et de la restitution du stockage
12
11
10
Text ( C)
9
8
7
6
5
4
0 2 4 6 8 10 12 14 16 18 20 22 24
Heure
Figure V-14 – Température moyenne horaire durant la saison de chauffage
124
Chapitre V. Pompe A Chaleur avec Stockage
stockage fonctionne aussi en régime cyclé, mais la fréquence de cyclages sera bien moins
importante que celle d’un PAC simple.
125
Chapitre V. Pompe A Chaleur avec Stockage
80
70
CHAUF REF
60
50
T( C)
40
30
20
10
0
0 100 200 300 400
H(kJ.kg-1)
Figure V-15 – Evolution de la température de la M37, en fonction de l’enthalpie lors du
chauffage et du refroidissement, obtenue par DSC ajustée
En revenant à l’analyse de l’effet du cyclage des PAC simples à faibles charges thermiques,
on rappelle qu’elles sont principalement dues à deux effets :
- le premier est dû à la montée de la température de condensation afin de surchauffer
la pièce et de limiter la fréquence des cycles
- le deuxième est dû aux phénomènes dynamiques lors de l’arrêt et du démarrage
avant de stabiliser les températures des parois de la PAC et les pressions du fluide.
Le premier effet dépend des limites sur la température de condensation imposées par la
régulation. Ces limites peuvent être variables dans une PAC simple afin de limiter le nombre
de cycles pour répondre aux besoins thermiques les plus faibles. Le deuxième effet dépend
principalement de la fréquence du cyclage. Le COP est d’autant plus pénalisé que le nombre
de cycles est élevé. L’hypothèse prise sera de considérer que les cinq premières minutes de
chaque démarrage concentrent les pertes dynamiques.
126
Chapitre V. Pompe A Chaleur avec Stockage
mêmes paramètres des PAC testés dans [DGUP] et [DGALN], le modèle donne une courbe
de variation de COP similaire à celle obtenue en expérimentation (Figure V-16 et
Figure V-17). Pour les températures extérieures élevées, le contrôle de la PAC air/eau
autorise une grande surchauffe de l’eau pour limiter la fréquence du cyclage. Ceci justifie
l’allure décroissante de la courbe réelle qui est différente de celle donnée par le modèle.
La Figure V-17 montre que cette PAC air/air n’est pas adaptée au besoin sollicité. Le modèle
calcule le COP thermodynamique ne prenant pas en compte la consommation des
auxiliaires. Il est clair qu’une pompe d’une PAC air/eau consomme plus d’énergie qu’un
ventilateur d’une PAC air/air.
4,0 5,5
5,0
3,5
4,5
3,0 4,0
COP
COP
3,5
2,5
3,0
2,0 Modèle 2,5 Modèle
2,0 Mesuré
1,5
Mesuré 1,5
1,0 1,0
-7 -4 -1 2 5 8 11 14 -7 -4 -1 2 5 8 11 14
Text ( C) Text ( C)
Figure V-16 – COP mesuré et modélisé d’une Figure V-17 – COP mesuré et modélisé d’une
PAC air/eau en fonction de la température PAC air/air en fonction de la température
extérieure extérieure
On remarque sur les deux figures, Figure V-16 et Figure V-17, et que les hypothèses
choisies pour la modélisation de la PAC air-eau d’une part et de la PAC air/air à
compresseur à vitesse variable donnent des COP supérieurs aux COP mesurés. Ces
hypothèses plus favorables correspondent à des machines mieux conçues.
Les hypothèses considérées sont :
le rendement global nominal des compresseur égal à 0,5 ; il suit la courbe de la
Figure V-4 pour le compresseur à vitesse variable
l’écart de température entre la température extérieure et la température d’évaporation
est de 5 K
la PAC avec stockage utilise la totalité de l’inertie du condenseur lors du cyclage
la température de condensation pour les systèmes à air est de 35 °C et pour le
système à eau de 40 °C
cette température augmente jusqu’à 5 K lors du cyclage
la puissance électrique du compresseur de la PAC à eau est de 220 W, celle de la
PAC à stockage de 200 W
la puissance nominale du compresseur à vitesse variable est de 110 W et varie entre
55 W et 220 W
le condenseur stockeur admet une inertie de 800 Wh, le circuit d’eau de la PAC
air/eau admet une inertie de 160 Wh, et le condenseur de la PAC air/air admet une
inertie de 20 Wh.
En considérant les mêmes sollicitations et en optimisant les puissances des compresseurs,
la Figure V-18 montre le COP thermodynamique des trois PAC à comparer : PAC air/eau
vitesse fixe, PAC air/air vitesse variable et PAC stockage (vitesse fixe).
127
Chapitre V. Pompe A Chaleur avec Stockage
6,0
5,5
5,0
4,5
COP
4,0
3,5
3,0
2,5 Air/Eau
2,0 Air/Air
1,5 Stock
1,0
-7 -4 -1 2 5 8 11 14
Text ( C)
Figure V-18 – COP thermodynamiques obtenus par le modèle pour les trois PAC
en fonction de la température extérieure
128
Chapitre V. Pompe A Chaleur avec Stockage
Air/Eau
750 2,5 2,5 2,6 2,7 2,3 2,4 2,4 2,5 2,6 2,7 2,8 2,9 3,0 3,2 4,0 4,0 4,1 4,2 4,3 4,4 4,5 4,6
700 2,2 2,2 2,3 2,3 2,3 2,4 2,5 2,5 2,6 2,7 2,8 2,9 3,0 3,2 3,9 4,0 4,1 4,1 4,2 4,3 4,4 4,5
650 2,2 2,3 2,3 2,3 2,4 2,4 2,5 2,5 2,6 2,7 2,8 2,9 3,0 3,1 3,8 3,9 4,0 4,1 4,2 4,3 4,3 4,4
600 2,3 2,3 2,3 2,3 2,4 2,4 2,5 2,6 2,6 2,7 2,8 2,9 3,0 3,1 3,8 3,9 3,9 4,0 4,1 4,2 4,3 4,4
Besoin (W)
550 2,3 2,3 2,3 2,3 2,4 2,4 2,5 2,5 2,6 2,7 2,8 2,9 3,0 3,1 3,7 3,8 3,9 3,9 4,0 4,1 4,2 4,3
500 2,3 2,3 2,3 2,3 2,4 2,4 2,5 2,5 2,6 2,7 2,8 2,9 3,0 3,1 3,7 3,7 3,8 3,9 4,0 4,0 4,1 4,2
450 2,3 2,3 2,3 2,3 2,4 2,4 2,5 2,5 2,6 2,7 2,8 2,8 2,9 3,0 3,6 3,7 3,7 3,8 3,9 3,9 4,0 4,1
400 2,3 2,3 2,3 2,3 2,4 2,4 2,5 2,5 2,6 2,6 2,7 2,8 2,9 3,0 3,5 3,6 3,6 3,7 3,8 3,9 3,9 4,0
350 2,2 2,3 2,3 2,3 2,3 2,4 2,4 2,5 2,6 2,6 2,7 2,8 2,9 3,0 3,4 3,5 3,6 3,6 3,7 3,8 3,9 3,9
300 2,2 2,2 2,3 2,3 2,3 2,4 2,4 2,5 2,5 2,6 2,7 2,7 2,8 2,9 3,4 3,4 3,5 3,5 3,6 3,7 3,8 3,8
250 2,2 2,2 2,2 2,3 2,3 2,3 2,4 2,4 2,5 2,6 2,6 2,7 2,8 2,9 3,3 3,3 3,4 3,5 3,5 3,6 3,7 3,7
200 2,2 2,2 2,2 2,2 2,3 2,3 2,3 2,4 2,5 2,5 2,6 2,7 2,7 2,8 3,2 3,2 3,3 3,4 3,4 3,5 3,6 3,6
150 2,1 2,2 2,2 2,2 2,2 2,3 2,3 2,4 2,4 2,5 2,5 2,6 2,7 2,8 3,1 3,2 3,2 3,3 3,3 3,4 3,5 3,5
100 2,1 2,1 2,1 2,2 2,2 2,2 2,3 2,3 2,4 2,4 2,5 2,6 2,6 2,7 3,0 3,1 3,1 3,2 3,2 3,3 3,4 3,4
50 2,1 2,1 2,1 2,1 2,2 2,2 2,2 2,3 2,3 2,4 2,4 2,5 2,6 2,7 2,9 3,0 3,0 3,1 3,1 3,2 3,3 3,3
1 2,0 2,0 2,1 2,1 2,1 2,1 2,2 2,2 2,3 2,3 2,4 2,5 2,5 2,6 2,8 2,9 2,9 3,0 3,0 3,1 3,2 3,2
-7 -6 -5 -4 -3 -2 -1 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14
Text (°C)
Figure V-19 – COP thermodynamique obtenu par le modèle pour la PAC air/eau
Le COP de la PAC air/air met en évidence une zone de fonctionnement optimale, qui est
celle de fonctionnement en continu entre les puissances minimale et maximale du
compresseur. En dehors de cette zone, le COP est fortement pénalisé par le cyclage
(besoins faibles) et par l’appoint de chauffage (besoins élevés). Une forte sensibilité à l’effet
de givrage est aussi clairement identifiable.
Air/Air
750 1,4 1,5 1,5 1,5 1,5 1,6 1,6 1,7 1,8 1,8 2,0 2,1 2,3 2,6 4,7 4,8 5,0 5,1 5,3 5,5 5,7 5,9
700 1,5 1,5 1,5 1,6 1,6 1,6 1,7 1,8 1,9 2,0 2,1 2,3 3,3 3,5 4,7 4,8 5,0 5,1 5,3 5,5 5,7 5,9
650 1,5 1,6 1,6 1,6 1,7 1,7 1,8 1,9 2,0 2,1 2,3 3,2 3,3 3,5 4,7 4,8 5,0 5,1 5,3 5,5 5,7 5,9
600 1,6 1,6 1,7 1,7 1,8 1,8 1,9 2,0 2,2 2,9 3,0 3,2 3,3 3,5 4,7 4,8 5,0 5,1 5,3 5,5 5,7 5,9
Besoin (W)
550 1,7 1,7 1,8 1,8 1,9 2,0 2,6 2,7 2,8 2,9 3,0 3,2 3,3 3,5 4,7 4,8 5,0 5,1 5,3 5,5 5,7 5,9
500 2,4 2,4 2,4 2,5 2,5 2,6 2,6 2,7 2,8 2,9 3,0 3,2 3,3 3,5 4,7 4,8 5,0 5,1 5,3 5,5 5,7 5,9
450 2,4 2,4 2,4 2,5 2,5 2,6 2,6 2,7 2,8 2,9 3,0 3,2 3,3 3,5 4,7 4,8 5,0 5,1 5,3 5,5 5,7 5,9
400 2,4 2,4 2,4 2,5 2,5 2,6 2,6 2,7 2,8 2,9 3,0 3,2 3,3 3,5 4,7 4,8 5,0 5,1 5,3 5,5 5,7 5,9
350 2,4 2,4 2,4 2,5 2,5 2,6 2,6 2,7 2,8 2,9 3,0 3,2 3,3 3,5 4,7 4,8 5,0 5,1 5,3 5,5 5,7 5,9
300 2,4 2,4 2,4 2,5 2,5 2,6 2,6 2,7 2,8 2,9 3,0 3,2 3,3 3,5 4,7 4,8 5,0 5,1 5,3 5,5 5,7 5,9
250 2,4 2,4 2,4 2,5 2,5 2,6 2,6 2,7 2,8 2,9 3,0 3,2 3,3 3,5 4,7 4,8 5,0 5,1 3,7 3,9 4,0 4,1
200 2,4 2,4 2,4 2,5 2,5 2,6 2,6 2,7 2,8 2,9 3,0 3,2 3,3 3,5 3,3 3,3 3,4 3,5 3,6 3,7 3,8 3,9
150 2,4 2,4 2,4 2,5 2,5 2,6 2,6 2,7 2,8 2,9 3,0 2,2 2,3 2,4 3,1 3,1 3,2 3,3 3,3 3,4 3,5 3,6
100 1,7 1,7 1,7 1,7 1,8 1,8 1,8 1,9 1,9 2,0 2,0 2,1 2,2 2,3 2,7 2,8 2,8 2,9 3,0 3,0 3,1 3,1
50 1,5 1,5 1,5 1,6 1,6 1,6 1,6 1,7 1,7 1,7 1,8 1,8 1,9 2,0 2,3 2,4 2,4 2,5 2,5 2,6 2,6 2,7
1 1,2 1,2 1,2 1,2 1,3 1,3 1,3 1,3 1,4 1,4 1,4 1,5 1,5 1,6 1,9 1,9 1,9 2,0 2,0 2,1 2,1 2,1
-7 -6 -5 -4 -3 -2 -1 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14
Text (°C)
Figure V-20 – COP thermodynamique obtenu par le modèle pour la PAC air/air
Le COP de la PAC à stockage montre une légère sensibilité au besoin et à l’effet du givrage.
129
Chapitre V. Pompe A Chaleur avec Stockage
Stockage
750 2,6 2,7 2,8 2,9 2,6 2,6 2,7 2,8 2,9 3,0 3,2 3,3 3,5 3,7 4,6 4,8 4,9 5,1 5,2 5,4 5,6 5,8
700 2,4 2,4 2,5 2,5 2,6 2,7 2,7 2,8 2,9 3,1 3,2 3,3 3,5 3,7 4,6 4,8 4,9 5,1 5,2 5,4 5,6 5,8
650 2,4 2,5 2,5 2,6 2,6 2,7 2,8 2,9 3,0 3,1 3,2 3,4 3,5 3,7 4,6 4,8 4,9 5,0 5,2 5,4 5,6 5,8
600 2,5 2,5 2,5 2,6 2,6 2,7 2,8 2,9 3,0 3,1 3,2 3,4 3,6 3,7 4,6 4,7 4,9 5,0 5,2 5,4 5,5 5,7
Besoin (W)
550 2,5 2,5 2,6 2,6 2,7 2,7 2,8 2,9 3,0 3,1 3,3 3,4 3,6 3,8 4,6 4,7 4,9 5,0 5,2 5,3 5,5 5,7
500 2,5 2,6 2,6 2,6 2,7 2,8 2,8 2,9 3,0 3,2 3,3 3,4 3,6 3,8 4,6 4,7 4,8 5,0 5,1 5,3 5,5 5,7
450 2,6 2,6 2,6 2,7 2,7 2,8 2,9 3,0 3,1 3,2 3,3 3,4 3,6 3,8 4,5 4,7 4,8 5,0 5,1 5,3 5,5 5,6
400 2,6 2,6 2,6 2,7 2,7 2,8 2,9 3,0 3,1 3,2 3,3 3,5 3,6 3,8 4,5 4,6 4,8 4,9 5,1 5,2 5,4 5,6
350 2,6 2,6 2,7 2,7 2,8 2,8 2,9 3,0 3,1 3,2 3,3 3,5 3,6 3,8 4,5 4,6 4,8 4,9 5,1 5,2 5,4 5,6
300 2,6 2,6 2,7 2,7 2,8 2,8 2,9 3,0 3,1 3,2 3,4 3,5 3,7 3,8 4,5 4,6 4,7 4,9 5,0 5,2 5,4 5,5
250 2,6 2,7 2,7 2,7 2,8 2,9 2,9 3,0 3,1 3,2 3,4 3,5 3,7 3,8 4,4 4,6 4,7 4,8 5,0 5,2 5,3 5,5
200 2,7 2,7 2,7 2,8 2,8 2,9 3,0 3,0 3,1 3,3 3,4 3,5 3,7 3,8 4,4 4,5 4,7 4,8 5,0 5,1 5,3 5,5
150 2,7 2,7 2,7 2,8 2,8 2,9 3,0 3,1 3,2 3,3 3,4 3,5 3,7 3,9 4,4 4,5 4,6 4,8 4,9 5,1 5,3 5,4
100 2,7 2,7 2,7 2,8 2,8 2,9 3,0 3,1 3,2 3,3 3,4 3,5 3,7 3,9 4,4 4,5 4,6 4,7 4,9 5,0 5,2 5,4
50 2,7 2,7 2,8 2,8 2,9 2,9 3,0 3,1 3,2 3,3 3,4 3,5 3,7 3,9 4,3 4,4 4,6 4,7 4,9 5,0 5,2 5,4
1 2,7 2,7 2,8 2,8 2,9 2,9 3,0 3,1 3,2 3,3 3,4 3,5 3,7 3,9 4,3 4,4 4,5 4,7 4,8 5,0 5,1 5,3
-7 -6 -5 -4 -3 -2 -1 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14
Text (°C)
Figure V-21 – COP thermodynamique obtenu par le modèle pour la PAC avec stockage
Les besoins obtenus par les simulations de la maison « Mozart passive » sont utilisés pour
calculer le COP saisonnier. La Figure V-22 montre les 30 200 occurrences de besoin en
fonction de la température extérieure avec un pas de temps de 10 minutes durant toute la
saison de chauffe.
Figure V-22 – Variation des besoins de la maison en fonction de Text durant la saison de
chauffe
En affectant les valeurs du COP aux données du besoin et de température extérieure durant
la saison de chauffe, les COP thermodynamiques saisonniers de la maison sont présentés
Tableau V-1.
Ce tableau montre que la PAC à stockage intégré au condenseur permet un gain potentiel
de 32 % et de 21 % sur le COP thermodynamique comparativement à une PAC air/eau et
une PAC air/air respectivement. En outre une PAC air/air surchauffe l’air ambiant et les
130
Chapitre V. Pompe A Chaleur avec Stockage
parois afin de limiter le cyclage. Ceci conduit à l’augmentation des déperditions, donc
l’augmentation des besoins, qui ne sont pas prises en compte dans le modèle.
La PAC avec stockage permet un effacement durant les heures de pic de demande
d’électricité, et donc émet moins de CO2 qu’une PAC simple. Elle offre un confort thermique
continu grâce à la grande inertie de son condenseur. Elle peut utiliser un compresseur à
vitesse fixe facile à contrôler et plus fiable. Elle présente un COP saisonnier potentiellement
plus élevé que les autres technologies.
Le Tableau V-2 présente les différentes valeurs du COP thermodynamique saisonnier en
variant les hypothèses suivantes : le rendement du compresseur, l’inertie du condenseur, la
portion de l’inertie utilisée et l’écart entre la température extérieure et la température
d’évaporation.
Tableau V-2 – Etude de sensibilité du COP
Cas Rendement Inertie du Portion Ecart COP
compresseur condenseur utilisée évaporation saisonnier
(-) (Wh) (%) (K) (-)
1 0,5 800 100 5 3,79
2 0,4 800 100 5 3,05
3 0,7 800 100 5 5,20
4 0,8 800 100 5 5,88
5 0,5 400 100 5 3,62
6 0,5 1600 100 5 3,89
7 0,5 800 50 5 3,61
8 0,5 800 100 3 4,02
9 0,5 800 100 7 3,59
Le cas de référence est celui pour un compresseur de rendement de 0,5, un condenseur
d’inertie de 800 Wh qui est utilisée entièrement durant le cyclage et un écart de températures
de 5 K sur l’évaporateur. Le COP de cette PAC peut atteindre 5,88 pour un compresseur de
rendement global de 0,8. D’autre part, le COP baisse de 0,17 en divisant la quantité de stock
par deux, et augmente de 0,9 en doublant l’inertie. Ce COP baisse aussi dans le cas où la
PAC n’utilise que la moitié du stockage thermique de son condenseur, même si la limite de
la température de condensation est plus basse, l’effet de cyclage reste dominant. L’efficacité
de l’évaporateur est représentée par l’écart entre la température extérieure et la température
d’évaporation. La variation 1 K de cet écart entraîne une variation du COP de l’ordre de 0,1.
Ce système peut aussi être amélioré en l’intégrant dans un système de ventilation double
flux. Ce système permet de récupérer la chaleur de l’air rejeté, ce qui réduit les besoins. Or,
le COP de cette PAC n’est quasiment pas affecté par le besoin, le gain de la récupération ne
sera pas pénalisé. D’autre part, ayant une grande inertie, cette PAC peut profiter du système
de dégivrage par soufflage d’air chaud sans créer un inconfort dans le bâtiment. En tenant
compte de ces deux effets (15 % de puissance récupérée), le COP thermodynamique
saisonnier peut atteindre 4,7 pour les mêmes paramètres que le cas de référence.
IV. Conclusions
Le problème majeur des systèmes de PAC pour le secteur résidentiel est que la puissance
calorifique générée est plus faible lorsque la température extérieure est plus basse, à cause
de l’augmentation de l’écart entre les températures de condensation et d’évaporation. Les
PAC sont dès lors surdimensionnées pour l’ensemble des occurrences de température et de
besoins hors le point nominal.
131
Chapitre V. Pompe A Chaleur avec Stockage
Les besoins de chauffage d’un logement BBC sont fortement aléatoires du fait qu’ils sont
fortement affectés par les apports internes et les apports solaires, ce qui rend l’adaptation
entre puissance thermique générée et besoins plus complexe.
Plusieurs solutions sont disponibles sur le marché, comme le fonctionnement en régime
cyclé et les compresseurs à vitesse variable. Ces solutions ont un effet négatif sur le COP.
Ces effets ont été évalués en fonction de la température extérieure et des besoins.
L’effet des cycles givrage/dégivrage sur le COP a été évalué en fonction de la température
extérieure.
Le cycle stockage/déstockage thermique dans la PAC à stockage peut conduire à un gain
énergétique en bénéficiant de la possibilité de stockage lors des heures les plus favorables
de la journée. Ce cycle est gagnant à partir de 2 °C de différence sur la température
extérieure durant le stockage et celle durant le déstockage. En outre ce cycle permet de
réduire le temps de fonctionnement de cette PAC, ce qui limite les pertes dues aux cycles
givrage/dégivrage et surtout aux marches/arrêts du régime cyclé à faible charge thermique.
La grande inertie du condenseur stockeur permet de réduire le temps de fonctionnement de
cette PAC et de limiter ainsi les pertes dues aux cycles givrage/dégivrage et marche/arrêt du
régime cyclé.
Un modèle permettant de simuler le COP réel d’une PAC a été développé en prenant en
compte tous les effets pratiques lors de son fonctionnement. Ce modèle a été validé par
comparaison à des travaux expérimentaux disponibles dans la littérature. Ce modèle a
prouvé la réduction de l’effet de variation des besoins et de l’effet de givrage sur le COP de
la PAC à stockage.
La PAC à stockage thermique au condenseur présente une amélioration du COP
thermodynamique saisonnier entre 20 et 30 % par rapport au COP d’une PAC air/air à
vitesse variable et d’une PAC air/eau. Le COP thermodynamique saisonnier de la PAC à
stockage à Trappes peut atteindre 5,88 en adaptant la bonne quantité de stockage et un
compresseur de rendement égal à 0,8.
132
Conclusions et Perspectives
Conclusions et Perspectives
Le chauffage domestique a une grande influence sur la demande d’électricité surtout aux
heures de pointe durant l’hiver. D’où l’idée d’intégrer un système de stockage thermique
journalier dans les pompes à chaleur qui, si cette technologie est utilisée massivement pour
le chauffage domestique, permettra alors de réduire le pic de demande électrique entre 18 et
20h00 et cela en décalant le temps de production de celui de consommation.
Le système PAC avec stockage intégré au condenseur n’est pas forcément moins efficace
qu’une PAC sans stockage. Ceci est dû au fait que les heures creuses de demande
d’électricité, qui sont bien sûr les heures les plus favorables pour un fonctionnement en
mode stockage de la PAC, correspondent aux heures pendant la journée aux heures où la
température extérieure est la plus élevée et donc où le COP de la PAC est potentiellement le
plus élevé. Au chapitre 5, des simulations réalistes ont montré l’impact du
stockage/déstockage sur l’efficacité énergétique de la PAC pour différents scénarios de
régulation et différentes journées types.
La réglementation thermique française, dite RT 2012, impose une consommation en énergie
primaire de 50 kWhEP.m-2.an-1 pour couvrir les cinq besoins (chauffage, Eau Chaude
Sanitaire, ventilation, climatisation et éclairage) des logements neufs. Une maison Mozart de
100 m² habitables, conçue pour respecter la RT 2012, nécessitera entre 24 et
36 kWh.m-2.an-1 en énergie finale pour le chauffage, selon les scénarios et les températures
de consigne considérés. Le jour le plus froid a été identifié comme le jour ayant les besoins
les plus élevés entre 18h00 et 20h00 (les heures d’effacement).
L’effacement entre 18h00 et 20h00 peut être réalisé avec une PAC stockant 800 Wh apte à
déstocker une puissance de 400 W en continu. Bien sûr, cet effacement pourra durer plus
longtemps pour les besoins thermiques plus faibles.
Le stockage thermique par chaleur latente peut être obtenu par des transformations
solide-solide, solide-liquide, solide-gaz et liquide-gaz. Le changement de phase adapté à
l’effacement d’une pompe à chaleur est celui de la transition solide-liquide grâce à la quantité
de chaleur délivrée et avec changement de volume acceptable lors de cette transformation.
Les critères de choix d’un MCP sont des critères chimiques (corrosion, dilatation/rétractation,
dénaturation), thermodynamiques et thermiques (capacité calorifique, chaleur latente,
conductivité thermique, diffusivité thermique), de cinétique de changement de phase (vitesse
de changements de phase, hystéresis) et économiques.
Les MCP peuvent être classés en trois grandes familles selon leurs natures chimiques :
organiques, inorganiques et présentant un point d’eutexie. Les propriétés de deux MCP
d’une même famille sont semblables.
La calorimétrie différentielle à balayage avec ajustement des propriétés par méthode inverse
permet d’évaluer la chaleur latente, les températures de fusion et de solidification et les
133
Conclusions et Perspectives
glissements de température lors de la transformation avec plus grande précision. Elle permet
aussi d’évaluer l’hystérésis du matériau.
Dans cette thèse, une méthode numérique par volumes finis a été développée, basée sur la
forme enthalpique de l’équation de l’énergie. Le modèle de changement de phase développé
est connecté à d’autres modèles pour former le système global afin de concevoir l’échangeur
stockeur et le système PAC.
Le condenseur stockeur est l’élément clé de la pompe à chaleur intégrant un stockage
thermique. Cet échangeur est le siège de l’échange entre les trois « matériaux actifs » du
système : le fluide frigorigène, l’air ambiant et le matériau de stockage.
Cet échangeur doit être capable de stocker 800 Wh et de décharger une puissance
thermique maximale de 400 W pour assurer l’effacement de la PAC entre 18h00 et 20h00.
La modélisation des échangeurs a été réalisée sur Dymola. Plusieurs structures ont été
modélisées et testées.
La géométrie des tubes insérés dans des tubes donne des résultats satisfaisants au niveau
de la répartition des flux thermiques. Par contre, elle présente un problème de compacité,
mais elle peut être prometteuse en apportant des améliorations de conception.
La géométrie tubes entre plaques a aussi donné des résultats satisfaisants au niveau de la
répartition des flux thermiques avec des améliorations limitées au niveau de la compacité par
rapport à la géométrie tube dans tube.
Le facteur qui a le plus grand impact sur l’efficacité d’un échangeur avec stockage par
chaleur latente, est la conductivité thermique du MCP. Une bonne structure est obtenue
lorsque les surfaces d’échange sont dimensionnées de sorte à équilibrer les trois
conductances surfaciques.
L’échangeur stockeur étant le siège de transfert thermique entre trois matériaux actifs (air,
MCP et fluide frigorigène), sa structure doit assurer l’équilibre des trois conductances afin de
limiter les écarts de températures. L’efficacité et même la rentabilité de cette PAC dépendent
de ces écarts de températures.
Un banc d’essais a été construit afin de tester l’échangeur stockeur et la PAC complète.
Un échangeur de première génération, qui est une amélioration de la géométrie tube dans
tube, a été conçu, modélisé, construit et testé. Les résultats des simulations sont conformes
aux résultats expérimentaux.
L’échangeur « première génération » est capable de stocker la quantité d’énergie et
décharger la puissance nécessaire. Par contre, il fonctionne à des débits de soufflage très
élevés et des températures de sorties relativement basses. De plus, il est volumineux.
Un échangeur de deuxième génération a été conçu et construit. Il est considéré comme
premier prototype puisqu’il vérifie complètement le cahier des charges d’un bon échangeur
stockeur.
Le prototype a été testé et validé sur le banc d’essais sur tous les régimes de
fonctionnement : charge, décharge et chauffage.
134
Conclusions et Perspectives
La PAC intégrant un stockage par chaleur latente a été conçue. La prochaine étape sera
d’étudier son comportement in situ lors d’une saison complète et de valider la stratégie de
contrôle tirant partie de l’inertie pour conclure sur son efficacité énergétique par rapport à
une PAC sans condenseur stockeur.
Le problème majeur des systèmes de PAC résidentielles est que la puissance calorifique
générée baisse lorsque la température extérieure elle-même baisse, et ce comportement
thermodynamiquement fatal est contradictoire avec l’accroissement du besoin de chauffage
avec cette baisse de la température extérieure. Ceci entraîne de fait un surdimensionnement
des PAC lorsque les besoins thermiques sont plus faibles que le point nominal de
conception.
De plus, les besoins de chauffage dans une maison BBC sont fortement variables et
aléatoires du fait qu’ils sont affectés par les comportements des occupants, les apports
internes et les apports solaires, ce qui rend l’adaptation entre puissance thermique générée
et besoins plus complexe.
Les régulations des PAC sont délicates à mettre au point et la perte d’efficacité énergétique
à faible charge thermique est notable, même avec des PAC à vitesse variable. Ces solutions
ont un effet négatif sur le COP. Ces effets ont été évalués en fonction de la température
extérieure et des besoins.
L’effet des cycles givrage/dégivrage sur le COP a été évalué en fonction de la température
extérieure.
Le cycle stockage/déstockage thermique dans la PAC à stockage peut conduire à un gain
énergétique en bénéficiant de la possibilité de stockage lors des heures les plus favorables
de la journée (températures d’air extérieur plus élevées). Le fonctionnement en mode
stockage est énergétiquement gagnant dès que la différence de températures entre le mode
stockage et déstockage est supérieure à 2 K. En outre, ce cycle permet de réduire le temps
de fonctionnement de cette PAC qui peut limiter les pertes dues aux cycles
givrage/dégivrage et aux cycles marche/arrêt du régime cyclé, du fait que la grande inertie
du condenseur stockeur permet d’accroître le temps de fonctionnement en régime
permanent stable.
Un modèle permettant de simuler le COP réel de PAC air/eau, air/air et air/air avec stockage
au condenseur a été développé en prenant en compte les irréversibilités associées au
niveau de température extérieure. Ce modèle a été validé par comparaison à des travaux
expérimentaux disponibles dans la littérature scientifique.
Ce modèle a prouvé la réduction de l’effet de variation des besoins et de l’effet de givrage
sur le COP de la PAC à stockage.
La PAC à stockage présente une amélioration du COP thermodynamique saisonnier entre
20 et 30 % respectivement par rapport au COP d’une PAC air/air à vitesse variable et une
PAC air/eau.
Le COP thermodynamique saisonnier de la PAC à stockage à Trappes peut atteindre 5,2 en
adaptant la bonne quantité de stockage et un compresseur de rendement égal à 0.7, ce qui
est un rendement atteignable.
Vu les avantages du fonctionnement en régime permanent continu d’une PAC, un système
de stockage améliore le COP saisonnier. Cette PAC devra comporter un échangeur capable
de stocker une chaleur suffisante pour assurer le besoin pendant des durées de 2 à 4
heures, voire plus. Plusieurs stratégies seront à étudier. La première stratégie serait de
produire le stock durant la nuit où la demande d’électricité est faible et les tarifs de
l’électricité aussi. Cette stratégie est pénalisée par le niveau de température extérieure lors
135
Conclusions et Perspectives
du fonctionnement de la PAC en mode stockage, car la température moyenne de nuit est
inférieure à la moyenne diurne. La deuxième stratégie consiste à stocker l’énergie durant le
pic de température extérieure qui se situe entre 12h00 et 15h00. D’où la chaleur est produite
avec le plus haut COP possible. La troisième stratégie sera de sous-dimensionner le
compresseur et de le faire fonctionner plus longtemps à faible puissance.
Un condenseur stockeur est aussi un réservoir de chaleur où des chaleurs peuvent être
récupérées d’autres sources, comme la chaleur solaire. Un quatrième fluide pourra alors être
mis en œuvre avec les trois compartiments principaux de l’échangeur stockeur. Ce
quatrième fluide peut être un caloporteur venant de capteurs solaires.
Un système de stockage à plus grande échelle peut être exploité aussi, le stockage
saisonnier. Il consiste à produire la chaleur durant l’été, et la stocker dans un réservoir à
MCP. Cette chaleur stockée sera utilisée pour chauffer le bâtiment durant l’hiver. La PAC ne
fonctionnera que durant l’été avec un COP très élevé. De plus, la chaleur absorbée par
l’évaporateur durant la charge pourra être utilisée pour climatiser le bâtiment durant l’été.
L’ensemble de ces solutions et leurs combinaisons constituent le futur des stockages
d’énergie thermiques froide ou chaude qui vont modifier nos schémas d’utilisation des
énergies thermiques et électrique.
136
Références bibliographiques
Références bibliographiques
M.A. ABDOLLY, D. RAPP, Theoretical and experimental numerical studies of
[ABD82] stratified thermocline storage hot water, Energy Convers. Manage. 23 (1982)
275.
Abhat, A., 1978, "Investigation of physical and chemical properties of phase
[ABH78] change materials for space heating/cooling applications", Proc. of the
International solar energy congress, New Delhi, India.
Abhat, A., 1981, " Development of a modular heat exchanger with an
[ABH81] integrated latent heat storage", Report No. BMFT FBT 81-050, Germany
Ministry of Science and Technology Bonn.
[ABH83] Abhat, A., 1983, "Low temperature latent heat thermal energy storage: heat
storage materials", Solar Energy, 30, 4, 313-331.
[AGC] AGC – Etude sur l’impact du verre dans le bilan environnemental d’une
maison individuelle – www.yourglass.com
F. Agyenim, N. Hewitt - « The development of a finned phase change material
(PCM) storage system to take advantage of off-peak electricity tariff for
[AGY10]
improvement in cost of heat pump operation » - Energy and Buildings, 2010,
42, 1552-1560
[ALU] www.alusion.com
E.S. Benta, Estudo experimental e numérico da solidificacao da agua ao redor
[BEN01] de tubos curvos, Doctoral Thesis, Universidade Estadual de Campinas,
Campinas, Sao Paulo, Brazil, 2001.
G. Bigot, G. – Etude et conception de systèmes air/air inversables utilisant des
[BIG01] mélanges à glissement de température – Thèse de doctorat / CEP - Mines
ParisTech – 2001.
[CHO70] S.H. Cho, J.E. Sunderland, Phase change of spherical bodies, International
Journal of Heat and Mass Transfer 13 (1970) 1231–1233.
137
Références bibliographiques
138
Références bibliographiques
[EPS] www.pcmproducts.net
Farid, M.M., Khudhair, A.M., Razack, S.A.K and Al-Hakkaj, S., 2004, "A review
[FAR04] on phase change energy storage:materials and applications", Energy
Conversion and Management, 45, 1597 - 1615.
J. FERNANDEZ-SEARA , F. UHIA, J. SIERES - Experimental analysis of a
[FER07] domestic electric hot water storage tank. Part I: Static mode of operation -
Applied Thermal Engineering 27 (2007) 129–136
[FRA] www.ifam-dd.fraunhofer.de
Francesco G., Marco P., Matthias H., ―A numerical model to evaluate the
[FRA12] thermal behaviour of PCM glazing system configurations‖ , Energy and
Buildings, Volume 54 ( 2012) P. 141-153
[GAR85] Garg, H.P., Mullick, S.C. and Bhargava, A.K., 1985, "Solar thermal energy
storage", D. Reidel Publishing Co.
[HAL71] Hale, D.V., Hoover, M.J., and O'Neill, M.J., 1971, "Phase change materials
Hand Book", Marshal Space Flight Center, Alabama.
Z. HAMMOU, M. LACROIX - A hybrid thermal energy storage system for
[HAM06] managing simultaneously solar and electric energy - Energy Conversion and
Management 47 (2006) 273–288
Hasnain, S., 1998, ―Review on sustainable thermal energy storage
[HAS98] technologies, Part I: heat storage materials and techniques,‖ Energy
Conservation and Management, 39, 1127-1138.
H. Henderson, D. Parker, Y. Huang - Improving DOE-2’s RESYS routine : user
[HEN00] defined functions to provide more accurate part load energy use and humidity
predictions.
J.R. Henriquez, Estudo numérico e experimental sobre vidros térmicos, M.Sc.
[HEN96] Thesis, Universidade Estadual de Campinas, Campinas, Sao Paulo, Brasil,
1996.
[HIL83] J.M. Hill, A. Kucera, Freezing a saturated liquid inside a sphere, International
Journal of Heat and Mass Transfer 26 (1983) 1631–1637.
Hiran, S., Suwondo, A., Mansoori, G.,1994, ―Characterization of alkanes and
[HIR94] paraffin waxes for application as phase change energy storage medium‖,
Energy Sources, 16, 117-128.
C. J. HO, C. H. CHU, Numerical simulation of heat penetration through a
[HO95] vertical rectangular phase change material/air composite cell, International
Journal of Heat Mass Transfer. Vol. 39, No. 9, pp. 1785-1795, 1996
M. HUANG , P. EAMES , B. NORTON - Thermal regulation of building-
[HUA04] integrated photovoltaics using phase change materials - International J. of
Heat and Mass Transfer 47, pp. 2715-2733
[INC] F. Incropera, D. Dewitt – Fundamental of heat and mass transfer 4th Edition –
John Wiley & Sons, 1996.
K.A.R. Ismail, J.R. Henrıquez, L.F.M. Moura, M.M. Ganzarolli, Ice formation
[ISM00] around isothermal radial finned tubes, Energy Conversion and Management
41 (2000) 585–605.
K.A.R. Ismail, J.R. Henriquez, Numerical and experimental study of spherical
[ISM02] capsules packed bed latent heat storage system, Applied Thermal
Engineering 22 (2002) 1705–1716
139
Références bibliographiques
[ISM97] K.A.R. Ismail, J.N.C. Castro, PCM thermal insulation in buildings, International
Journal of Energy Research 21 (1997) 1281–1296.
K.A.R. Ismail, J.R. Henriquez, U-values, optical and thermal coefficients of
[ISM98a] composite glass systems, Solar Energy Materials and Solar Cell 52 (1998)
155–182.
[JEC93] T.B. Jeckel, J.W. Mitchell, S.A. Klein, Modelling of ice-storage tanks, ASHRAE
Transaction 99 (1) (1993) 1016– 1023.
[KAU91] Kauranen, P., Peippo, K. and Lund, P.D, 1991, ―An organic PCM storage
system with adjustable melting temperature,‖Solar Energy, 46, 275-278.
[KAY00] Sari, A. and Kaygusuz, K., 2000, ―Energy and exergy calculations of latent
heat storage systems,‖ Energy Sources, 22, 117-126.
[KAY99] Kaygusuz, K., 1999, ―The viability of thermal energy storage‖, Energy Sources,
21, 745-755.
E. Kinab – Optimisation des performances non nominales des pompes à
[KIN09] chaleur réversibles pour le secteur tertiaire - Thèse Ecole des Mines de Paris,
Décembre 2009.
[LAN78] Lane, G.A., 1978, " Macro–encapsulation of PCM", Report No. oro/5117-8,
Dow Chemical Company, Midland, Michigan, 152.
[LAN83] Lane G., 1983, "Latent heat materials", Volume 1, CRC Press, Boca Raton,
Florida.
[LAN89] Lane, G.A., 1989, " Hand book of thermal design, Chapter 1:Phase change
thermal storage materials’, Edited by C. Guyer, McGraw Hill Book Co.
A. LAZARO , P. DOLADO , J. MARIN , B. ZALBA - PCM-air heat exchangers
[LAZ09a] for free-cooling applications in buildings: Empirical model and application to
design - Energy Conversion and Management 50 (2009) 444–449
A. LAZARO , P. DOLADO , J. MARIN , B. ZALBA - PCM–air heat exchangers
[LAZ09b] for free-cooling applications in buildings: Experimental results of two real-scale
prototypes - Energy Conversion and Management 50 (2009) 444–449
Z.LIU , Z. WANG, C. MA - An experimental study on the heat transfer
[LIU06] characteristics of a heat pipe heat exchanger with latent heat storage - Energy
Conversion and Management 47 (2006) 967–991
140
Références bibliographiques
[MIL84] W. Miller – Frosting experiments for a heat pump having a one-row spine-fine
outdoor coil – ASHRAE Transaction, 4 (1984), p. 1009-1025
L.F. Milanez, Desenvolvimento teorico e verificacao experimental de modelos
[MIL84a] para a solidificacao em geometria esferica, Doctoral Thesis, Universidade
Estadual de Campinas, Campinas, Sao Paulo, Brasil, 1984.
L.F. Milanez, K.A.R. Ismail, Solidification in spheres, theoretical and
[MIL84b] experimental investigation, in: 3rd International Conference on Multi-Phase
and Heat Transfer, Miami, 1984.
[NAK95] A.E. Nakhi, Adaptive construction modelling within whole building dynamic
simulation, Ph.D. thesis, University of Strathclyde, Glasgow 1995.
A. Nuntaphan , T. Kiatsiriroat , C.C. Wang – Heat transfer and friction
[NAN04] characteristics of crimped spiral finned heat exchangers with dehumidification
– Applied Thermal Engineering 25 (2005) 327–340
141
Références bibliographiques
[PAP99] PAPON (P.), LEBLOND (J.) et MEIJER (P.). – Physique des Transitions de
Phase - Dunod (1999)
[PED73] R.I. Pedroso, G.A. Domoto, Perturbation solutions for spherical solidification of
saturated liquids, Journal of Heat Transfer 95 (1973) 42–46.
T. Perrotin – Etude côté d’air d’échangeurs à haute efficacité énergétique
[PER04] utilisés dans les systèmes de conditionnement d’air résidentiels - Thèse Ecole
des Mines de Paris, Juin 2004
[PIL76] Pillai, K.K. and Brinkwarth, B.J.,1976, "The storage of low grade thermal
energy using phase change materials,"Applied Energy, 2, 205-216.
M. Prud’homme, T.H. Nguyen, D.L. Nguyen, A heat transfer analysis for
[PRU89] solidification of slabs, cylinders, and spheres, Journal of Heat Transfer 111
(1989) 699–705.
Regin, A. F., Solanki, S. C., and Saini, J. S., 2008, ―Heat Transfer
[REG08] Characteristics of Thermal Energy Storage System Using PCM Capsules: A
Review,‖ Renewable Sustainable Energy Rev., 12, pp. 2438–2458.
D.S. Riley, S.T. Smith, G. Poots, The inward solidification of spheres and
[RIL74] circular cylinders, International Journal of Heat and Mass Transfer 17 (1974)
1507–1516.
[RT12] www.rt-batiment.fr
[RTE] www.rte-france.com
[RUB] www.rubitherm.com
Sharma, A., Sharma, S.D. and Buddhi, D., 2002, " Accelerated thermal cycle
[SHA02] test of acetamide, stearic acid and paraffin wax for solar thermal latent heat
storage applications", Energy Conversion and Management, 43, 1923 - 1930.
, S. Shao, W. Shi, X. Li, H. Chen, Performance representation of variable
[SHA04] speed compressor for inverter air-conditioners based on experimental data
International journal of refrigeration 27 (2004) 805-815.
[SHA99a] Sharma, S.D., Buddhi, D. and Sawhney, R.L., 1999, ―Accelerated thermal
cycle test of latent heat-storage materials,‖Solar Energy, 66, 483-490.
Sharma, S.D., 1999, "Study of thermal energy storage in phase change
[SHA99b] materials for low temperature solar applications", Ph.D. Dissertation, Devi
Ahilya University, Indore, India.
[SHI71] Y.P. Shih, T.C. Chou, Analytical solution for freezing a saturated liquid inside
or outside spheres, Chemical Engineering Science 26 (1971) 1787–1793.
A. SHUKLA , D. BUDDHI, R.L. SAWHNEY - Solar water heaters with phase
[SHU09] change material thermal energy storage medium: A review - Renewable and
Sustainable Energy Reviews 13 (2009) 2119–2125
[SP] www.sp.se
U. Stritih – « An experimental study of enhanced heat transfer in rectangular
[SRT04] PCM thermal storage » - International Journal of Heat and Mass Transfer,
2004, 47, 2841-2847
U. Stritih - « Heat transfer enhancement in latent heat thermal storage system
[STR03]
for buildings » - Energy and Buildings, 2003, 35, 1097-1104
142
Références bibliographiques
[TH12] www.bulletin-officiel.developpement-durable.gouv.fr
S. THIERS – Bilans énergétiques et environnementaux de bâtiments à
[THI08] énergie positive – Thèse de doctorat à l’Ecole Nationale Supérieure des Mines
de Paris – 2008
G. Tieyi, Z. Yinping, G. Xinshi, Theoretical analysis of heat transfer in a cool
[TIE96] storage tank filled with rectangular enclosures containing phase change
material, in: W. Buxuan (Ed.), Heat Transfer Science and Technology, Higher
Education Press, Beijing, China, 1996, pp. 734–739.
P. TITTELEIN – Environnements de simulation adaptés à l’étude du
[TIT08] comportement énergétique des bâtiments basse consommation – Thèse de
doctorat à l’Université de Savoie – 2008
R. Velraj, R.V. Seeniraj, B. Hafner, C. Faber, K. Schwarzer, Experimental
[VEL97] analysis and numerical modeling of inward solidification on a finned vertical
tube for a latent heat storage unit, Solar Energy 60 (1997) 281–290.
Veraj, R., Seeniraj, B., Hafner B., Faber, C. and Schwarzer, K.,1998, ―Heat
[VER98] transfer enhancement in a latent heat storage system,‖ Solar Energy,65, 171-
180.
Wang, X., Lu, E., Lin, W., Liu, T., Shi, Z., Tang, R. and Wang, C., 2000. ―Heat
[WAN00] storage performance of the binary systems neopentyl glycol/pentaerythritol
and neopentyl glycol/trihydroxy menthyl-aminomethane as solid phase change
materials,‖ Energy Conservation and Management, 41, 129-134.
[WIL89a] M.W. Wildin, C.R. Truman, Performance of stratified vertical cylindrical thermal
storage tanks—Part I (scale model tank), ASHRAE Trans. 95 (1) (1989) 1086.
[YAO89] Yao, L. S., and Prusa, J., 1989, ―Melting and Freezing,‖ Adv. Heat Transfer,
19, pp. 1–95.
Zalba, B., Marin, J.M., Cabeza, L.F. and Mehling, H., 2003, "Review on
[ZAL03] thermal energy storage with phase change:materials, heat transfer analysis
and applications, Applied Thermal Engineering, 23, 251 - 283.
M. Zukowski – « Experimental study of short term thermal energy storage unit
[ZUK07a] based on enclosed phase change material in polyethylene film bag » - Energy
Conversion and Management, 2007, 48, 166-173
M. ZUKOWSKI - Mathematical modeling and numerical simulation of a short
[ZUK07b] term thermal energy storage system using phase change material for heating
applications - Energy Conversion and Management 48 (2007) 155–165
143
ETUDE ET CONCEPTION D’UNE POMPE A CHALEUR RESIDENTIELLE
INTEGRANT UN STOCKAGE PAR CHALEUR LATENTE
Mots clés : pompe à chaleur, stockage thermique, MCP, chauffage résidentiel, Modelica