modelisation dynamique d'une pompe a chaleur air/eau ... - iusti
modelisation dynamique d'une pompe a chaleur air/eau ... - iusti
modelisation dynamique d'une pompe a chaleur air/eau ... - iusti
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
12èmes Journées Internationales de Thermique<br />
------------------------------------------------------------------------------------------------------------------------------------------------<br />
MODELISATION DYNAMIQUE D’UNE POMPE A CHALEUR AIR/EAU<br />
FONCTIONNANT AU CO 2<br />
Imed GUITARI, Philippe HABERSCHILL, André LALLEMAND<br />
CETHIL–UMR 5008/INSA-UCBL-CNRS, 20 av Albert Einstein, 69621 VILLEURBANNE cedex, France<br />
philippe.haberschillnsa-lyon.fr<br />
RESUME<br />
Le modèle proposé est basé sur l’intégration<br />
temporelle des variables d’états qui caractérisent le<br />
fonctionnement de la <strong>pompe</strong> à <strong>chaleur</strong>. Il permet de<br />
visualiser l’évolution des températures et des pressions à<br />
la sortie de chaque composant de l’installation, du débit<br />
massique du CO 2 dans le détendeur et des températures de<br />
l’<strong>air</strong> et de l’<strong>eau</strong> à l’aval de l’évaporateur et de l’échangeur<br />
haute pression. Ces variables ont fait l’objet de mesures<br />
sur le banc d’essais de la machine modélisée, ce qui a<br />
permis la validation du modèle. Certains paramètres de<br />
fonctionnement de l’installation, inaccessibles expérimentalement<br />
telle que l’évolution des températures le<br />
long des échangeurs, peuvent l’être grâce au modèle.<br />
INTRODUCTION<br />
Le CO 2 est l’une des meilleures alternatives aux<br />
fluides frigorigènes classiques. Les premiers travaux<br />
effectués au CETHIL concernaient le développement d’un<br />
banc d’essais d’une <strong>pompe</strong> à <strong>chaleur</strong> <strong>air</strong>/<strong>eau</strong> utilisant ce<br />
fluide et la mise au point d’un modèle en régime statique<br />
pour caractériser le fonctionnement de ces installations<br />
[1]. Ces travaux sont poursuivis actuellement par le<br />
développement d’un modèle en régime variable décrit<br />
dans ce présent travail.<br />
ARCHITECTURE DU MODELE<br />
Mêmes si les installations fonctionnant au CO 2 ne<br />
présentent pas les mêmes modes de fonctionnement elles<br />
sont constituées des même composants (Figure 1) que les<br />
installations classiques : un compresseur semi hermétique,<br />
un échangeur haute pression du type tubes et calandre, un<br />
échangeur intermédi<strong>air</strong>e, une vanne de détente, un<br />
évaporateur à <strong>air</strong> et une bouteille de réserve.<br />
se<br />
ee<br />
2 Ehp<br />
3<br />
expérimentales. Les variables de sortie du modèle sont les<br />
paramètres de fonctionnement de chacun des éléments de<br />
l’installation, pour chaque volume de contrôle des<br />
échangeurs, une résolution par volumes finis ayant été<br />
adoptée. Ainsi, les variables représentées en gras (Figure<br />
2) sont des grandeurs vectorielles. Les résultats de calcul<br />
sont enregistrés dans des fichiers relatifs au<br />
fonctionnement de chacun des éléments de l’installation<br />
pour chaque pas de temps. Parmi ces paramètres de sortie<br />
on note :<br />
• la température du CO 2 dans chacun des volumes de<br />
contrôle des échangeurs T<br />
hp<br />
et T<br />
ev<br />
;<br />
• le coefficient de transfert au sein du CO 2 dans chaque<br />
volumes de contrôle des échangeurs, α<br />
hp<br />
et α<br />
ev<br />
;<br />
• la température des parois pour chaque volume<br />
contrôle des échangeurs T<br />
tev<br />
et T<br />
thp<br />
et la température de<br />
la calandre T<br />
c<br />
;<br />
• les pressions d’aspiration et de refoulement P 1<br />
et<br />
P<br />
2<br />
;<br />
• le débit massique dans le compresseur et dans le<br />
détendeur m&<br />
cp<br />
et m&<br />
d<br />
;<br />
• la température de sortie de l’<strong>eau</strong> T<br />
se<br />
;<br />
• la température de l’<strong>air</strong> à l’aval de l’évaporateur T<br />
sa<br />
;<br />
• la puissance de compression W & cp<br />
;<br />
• la température des parois du compresseurT pc<br />
.<br />
S<br />
T<br />
ee<br />
m&<br />
e<br />
T ea<br />
V & Na<br />
Modèle en<br />
régime<br />
<strong>dynamique</strong><br />
T ,T hp thp<br />
T<br />
T<br />
e<br />
, Tc<br />
, α<br />
hp<br />
ev<br />
, Ttev<br />
, α<br />
ev<br />
m &<br />
, &<br />
cp<br />
m d<br />
P1 , P 2<br />
T<br />
se<br />
, T sa<br />
W& , T c pc<br />
t<br />
7<br />
Bouteille<br />
de réserve<br />
1<br />
6<br />
ea<br />
Eint<br />
Evap<br />
sa<br />
4<br />
Vanne<br />
de détente<br />
5<br />
Figure 1 : Installation modélisée<br />
Les paramètres d’entrée du modèle principal sont ceux<br />
qui définissent les conditions de fonctionnement de<br />
l’installation. Ces variables d’entrée sont les suivantes :<br />
• la température d’entrée et le débit d’<strong>eau</strong> T ee et m& ;<br />
e<br />
Paramètres dimensionnels de l’installation<br />
Figure 2 : Variables d’entrée et de sortie du modèle<br />
Chaque élément de l’installation est étudié séparément<br />
afin de définir les fonctions reliant ses variables d’entrée<br />
et ses variables de sortie en fonction du temps. Ces<br />
fonctions sont issues de l’analyse physique du<br />
fonctionnement de chacun de ces éléments, des principes<br />
fondamentaux de la thermo<strong>dynamique</strong> et de corrélations<br />
des transferts thermiques convectifs issues de la<br />
littérature :<br />
• la corrélation de Gneilinski pour les transferts dans le<br />
domaine supercritique [2] ;<br />
• la corrélation de Chen pour la condensation du CO 2 ;<br />
• la corrélation de Bennet-Chen modifiée pour les<br />
transferts à l'ébullition ;<br />
• la corrélation de Dittus-Boelter pour l’<strong>eau</strong> ;<br />
• un coefficient constant du coté de l’<strong>air</strong>, déterminé à<br />
• la température et la vitesse de l’<strong>air</strong> T<br />
ea<br />
et v<br />
a<br />
;<br />
• l’ouverture de la vanne de détente S ;<br />
• la vitesse de rotation du compresseur N .<br />
Les variables d’entrée évoluent en fonction du temps<br />
et sont lues à partir d’un tabl<strong>eau</strong> de valeurs<br />
Tanger, Maroc du 15 au 17 Novembre 2005 81
12èmes Journées Internationales de Thermique<br />
------------------------------------------------------------------------------------------------------------------------------------------------<br />
partir des essais sur l’installation modélisée.<br />
P2<br />
Les modèles élément<strong>air</strong>es sont intégrés dans un où δ = est le taux de compression, et v1<br />
est le volume<br />
P1<br />
modèle principal les réunissant pour simuler le cycle<br />
décrit par le CO 2 dans l’installation (Figure 3). A chaque<br />
massique du CO 2 à l’aspiration du compresseur.<br />
pas de temps, les variables de sortie d’un élément de<br />
Le coefficient polytropique k dépend du volume<br />
l’installation constituent les variables d’entrée du modèle massique à l’aspiration v 1<br />
et du taux de compression<br />
suivant en respectant le sens de circulation du CO 2 dans selon l’équation (3).<br />
l’installation.<br />
n<br />
k = ( a v1 + b) δ<br />
(3)<br />
Début T ea<br />
où a , b et n sont des coefficients déduits du<br />
dT<br />
pc<br />
dt<br />
dm hp<br />
dt<br />
dP 2<br />
dt<br />
dh hp<br />
dt<br />
dT thp<br />
dt<br />
dT e<br />
dt<br />
dT c<br />
dt<br />
dT tev<br />
dt<br />
dh ev<br />
dt<br />
dP 1<br />
dt<br />
dm ev<br />
dt<br />
Calcul du pas de<br />
temps<br />
Integration des<br />
fonctions d’état<br />
T<br />
pc<br />
FIN<br />
N<br />
Compresseur<br />
T c<br />
T thp T e<br />
h P m& h 2<br />
m&<br />
hp 2<br />
e<br />
c<br />
Echangeur<br />
haute presssion<br />
P P h 3 h 7 m& m&<br />
2<br />
1<br />
c d<br />
Echangeur intermédiaure<br />
P<br />
2<br />
P<br />
1 h 4 S<br />
Détendeur<br />
P m h ev T tev h 5 m& m&<br />
1 ev<br />
m&<br />
d c a T ea<br />
Evaporateur<br />
h 6<br />
dt<br />
8<br />
P<br />
2<br />
P<br />
1<br />
h<br />
1<br />
Saturation<br />
1<br />
Figure 3 : Schéma synoptique du modèle global<br />
Afin de diminuer le temps du calcul, les propriétés<br />
thermo<strong>dynamique</strong>s du CO 2 sont calculées par<br />
interpolation liné<strong>air</strong>e à partir de tables générées par le<br />
logiciel REFPROP6 ® . Sur la Figure 3, les paramètres<br />
représentés dans deux cercles concentriques sont des<br />
variables d’états sujets d’intégration temporelle selon la<br />
relation (1), les paramètres dans des cercles sont des<br />
variables issues des calculs, les paramètres représentés<br />
dans des carrés sont des constantes liées aux dimensions<br />
et aux conditions de fonctionnement et les traits gras<br />
représentent des sous modèles ou des étapes de calcul.<br />
⎛ df ⎞<br />
f ( t + dt) = f ( t) + ⎜ ⎟ dt<br />
(1)<br />
⎝ dt ⎠<br />
f est l’une des fonctions d’état du modèle.<br />
MODELES ELEMENTAIRES<br />
Compresseur<br />
Le calcul de l’enthalpie de sortie du CO 2 h<br />
2<br />
avec le<br />
logiciel REFPROP6 ® nécessite de connaître un paramètre<br />
complément<strong>air</strong>e à la pression de refoulement P<br />
2<br />
. Le<br />
paramètre choisi est le volume massique au refoulement.<br />
En considérant la compression polytropique, il est calculé<br />
selon l’équation (2).<br />
v<br />
2<br />
− 1<br />
k<br />
= v δ (2)<br />
1<br />
fonctionnement du compresseur. Le débit massique du<br />
CO 2 m&<br />
cp<br />
est calculé selon l’équation (4) en considérant<br />
que la vitesse de rotation du compresseur est celle de<br />
synchronisme du moteur électrique ( N = 1500 tr/min).<br />
N Vc<br />
m& cp<br />
= η v<br />
(4)<br />
60 v 1<br />
Dans cette expression, le rendement volumique η<br />
v<br />
est<br />
calculé à partir du taux de volume mortτ et du taux de<br />
compression δ.<br />
1<br />
= +<br />
⎛ ⎞<br />
⎜ −<br />
k<br />
η<br />
v<br />
1 τ 1 δ ⎟<br />
(5)<br />
⎝ ⎠<br />
Cependant, pour être validé par les expériences, le<br />
taux de volume mort ne doit pas être constant. Il obéit à<br />
l’équation :<br />
m<br />
τ = cδ<br />
(6)<br />
c et m sont des coefficients déduits de l’expérience.<br />
La puissance électrique d’entraînement du compresseur<br />
est calculée selon l’équation (7).<br />
W & cp<br />
= m&<br />
cp<br />
( h − h ) + Q&<br />
+ pc<br />
Q&<br />
2 1<br />
amb<br />
(7)<br />
où Q & amb<br />
est la puissance thermique échangée entre le<br />
compresseur et le milieu ambiant. Cette puissance<br />
englobe les pertes dues au frottement mécanique, la<br />
puissance dissipée par la <strong>pompe</strong> de lubrification, les<br />
pertes fer et les échanges thermiques avec l’ambiance.<br />
Elle est donnée par la relation empirique suivante :<br />
Qamb<br />
= d + e( Tpc<br />
− Tamb<br />
)<br />
(8)<br />
où d et e sont des coefficients empiriques, T<br />
pc<br />
et T amb<br />
sont respectivement la température de la paroi du<br />
compresseur et la température ambiante. Q & pc<br />
est la<br />
puissance thermique accumulée dans la paroi du<br />
compresseur et est calculée selon la relation :<br />
⎛ dT ⎞<br />
Q& pc<br />
= m cp<br />
⎜<br />
⎟<br />
pc<br />
(9)<br />
⎝ dt ⎠<br />
m est la masse du compresseur, cp est sa capacité<br />
calorifique. La variation de la température des parois du<br />
compresseur en fonction du temps est régie par un bilan<br />
thermique selon la relation :<br />
dTpc<br />
⎛ ⎛ T2<br />
+ T1<br />
⎞<br />
⎞ −<br />
= ⎜ K ⎜ − T ⎟ − ( − ) ( ) 1<br />
pc<br />
e Tpc<br />
Tamb<br />
2<br />
⎟ m cp (10)<br />
dt ⎝ ⎝<br />
⎠<br />
⎠<br />
K et e sont les conductances thermiques respectivement<br />
entre le gaz carbonique et la paroi du compresseur, et<br />
entre la paroi et l'ambiance, T<br />
1<br />
la température<br />
d’aspiration, T<br />
2<br />
la température de refoulement.<br />
Échangeurs de <strong>chaleur</strong><br />
L’échangeur haute pression et l’évaporateur sont<br />
discrétisés en n hp<br />
et n<br />
ev<br />
éléments de volume (Figure 4)<br />
dans lesquels les propriétés physiques sont considérées<br />
Tanger, Maroc du 15 au 17 Novembre 2005 82
12èmes Journées Internationales de Thermique<br />
------------------------------------------------------------------------------------------------------------------------------------------------<br />
uniformes que ce soit pour le fluide, pour les parois des En considérant que le coefficient de transfert avec la<br />
tubes, pour l’<strong>eau</strong> ou pour la calandre. L’évaporateur est calandre et les tubes est la même, on a pour le bilan sur<br />
un échangeur à courants croisés, sa discrétisation est plus l’<strong>eau</strong> :<br />
simple. L’échangeur intermédi<strong>air</strong>e est de petite taille on dTei<br />
me<br />
cpe<br />
= m&<br />
ecpe<br />
( Te i+ 1<br />
− Te i<br />
) + αe i<br />
( Ste( Tti<br />
− Te i<br />
) + Sc( Tci<br />
− Te i<br />
)<br />
néglige son inertie et il est considéré comme un seul dt<br />
élément de volume.<br />
La calandre est suffisamment isolée thermiquement,<br />
pour calculer sa variation de température on négligeant les<br />
Tube<br />
Calandre<br />
m& V i<br />
m& m&<br />
i<br />
i−1<br />
V i+1 i+ 1<br />
i+1<br />
i<br />
V i-1<br />
∆X i-1<br />
m&<br />
e<br />
CO 2<br />
Figure 4 : Discrétisation de l’échangeur haute pression<br />
Pour chaque élément de volume les équations de<br />
conservation de la masse et de l’énergie sont appliquées.<br />
En négligeant la variation de l’enthalpie massique liée à la<br />
dissipation et par exemple pour l’échangeur haute<br />
pression, l’équation de conservation de l’énergie pour<br />
chaque élément de volume i de l’écoulement de CO 2 est :<br />
dh<br />
dP<br />
−1<br />
(<br />
−1<br />
−1<br />
( )<br />
⎟ ⎞<br />
= m& i<br />
hi<br />
− m&<br />
i<br />
hi<br />
+ αi<br />
St<br />
T<br />
i<br />
−Tt<br />
i<br />
+ Vi<br />
mi<br />
(11)<br />
dt<br />
dt ⎠<br />
Vi<br />
et<br />
St<br />
sont le volume et la surface interne du volume de<br />
contrôle. La pression est supposée uniforme dans<br />
l’échangeur. Ainsi, la variation de la pression de<br />
refoulement en fonction du temps est calculée en<br />
considérant que la masse volumique du CO 2 pour chaque<br />
volume de contrôle de l’échangeur haute pression ρ<br />
i<br />
est<br />
une fonction d’état dépendant de l’enthalpie massique et<br />
de la pression :<br />
dρ<br />
⎛ ⎞ dh<br />
i hpi<br />
⎜<br />
∂ρi<br />
⎟ ⎛ ∂ρi<br />
⎞ dP2<br />
=<br />
+ ⎜ ⎟<br />
dt<br />
h<br />
(12)<br />
⎝<br />
∂<br />
hpi ⎠<br />
dt ⎝ ∂P<br />
⎠h<br />
dt<br />
P<br />
Pour chaque volume de contrôle, les variations de la<br />
masse volumique par rapport à l’enthalpie massique, à<br />
pression constante, et par rapport à la pression, à enthalpie<br />
massique constante, sont calculées à partir des tables du<br />
fluide frigorigène pour des transformations isobare et<br />
isenthalpique. La variation de l’enthalpie massique par<br />
rapport au temps est calculée à partir de l’équation (11).<br />
En multipliant l’équation (12) par le volume considéré, en<br />
sommant ces termes pour tous les éléments de volume de<br />
l’échangeur haute pression et en appliquant le principe de<br />
conservation de la masse, on obtient :<br />
⎛ n<br />
⎞<br />
⎜<br />
⎛ ⎞<br />
n<br />
⎜<br />
∂ρ<br />
dh<br />
i hp i<br />
⎟<br />
dP ⎛ ∂ρi<br />
⎞<br />
m& − =<br />
+ ⎜ ⎟<br />
⎟<br />
cp<br />
m&<br />
d<br />
V<br />
⎜∑<br />
∑ ⎟<br />
(13)<br />
1<br />
⎝ ⎝<br />
∂hhp i ⎠<br />
dt dt i ⎝ ∂P<br />
⎠<br />
P<br />
h<br />
⎠<br />
A partir de cette équation on calcule la seule inconnue<br />
qui est la variation de la pression en fonction du temps.<br />
Les tubes de l’échangeur haute pression et de<br />
l'évaporateur sont en contact avec le CO 2 d’un côté et<br />
avec de l’<strong>eau</strong> ou de l’<strong>air</strong> de l’autre. La variation de la<br />
température de ces tubes dépend des flux échangés avec<br />
ces deux fluides :<br />
dTi<br />
mt<br />
cpt<br />
= ( α<br />
i<br />
Sti<br />
( Ti<br />
− Tti<br />
) + αei<br />
Ste( Tei<br />
− Tt i<br />
)) (14)<br />
dt<br />
avec : S<br />
ti<br />
et Ste<br />
les surfaces interne et externe du tube, m<br />
t<br />
la masse du tube, pour chaque volume de contrôle. α<br />
i<br />
et<br />
α<br />
ei<br />
les coefficients de transfert du côté du CO 2 et du côté<br />
des fluides calo et frigoporteurs.<br />
<strong>eau</strong><br />
transferts thermiques avec l’extérieur :<br />
dTci<br />
−1 −1<br />
= α<br />
ei<br />
Sc( Tc i<br />
− Te i<br />
) mc<br />
cpc<br />
(15)<br />
dt<br />
avec m<br />
c<br />
et cp<br />
c<br />
la masse d’un volume élément<strong>air</strong>e et la<br />
capacité calorifique de la paroi, S<br />
c<br />
la surface d’échange<br />
<strong>eau</strong>-calandre correspondante.<br />
Bouteille anti-coup de liquide<br />
La variation de la basse pression par rapport au temps<br />
est calculée au niv<strong>eau</strong> de la bouteille anti-coup de liquide<br />
placée à la sortie de l’évaporateur. L’équation (16) est<br />
issue de la dérivée par rapport au temps de l’équation<br />
d’état des gaz parfaits sur le volume de la bouteille de<br />
réserve.<br />
dP ( m&<br />
d<br />
v − m&<br />
cp<br />
v )<br />
1<br />
6<br />
7<br />
=<br />
dt ⎛ 1 1 dT ⎞<br />
(16)<br />
⎜ ⎟<br />
− Vbt<br />
P T dP<br />
⎝ 1 7 ⎠<br />
dT à la pression P 1<br />
est obtenue à partir des tables à la<br />
dP<br />
saturation. V<br />
bt<br />
est le volume de la bouteille de réserve.<br />
L’inertie thermique de ce composant est négligée.<br />
Détendeur<br />
Le modèle de détendeur permet de calculer le débit<br />
massique m&<br />
d<br />
qui traverse le détendeur au cours d’une<br />
détente isenthalpique du fluide frigorigène d’une<br />
température d’entrée T 4<br />
à la température d’évaporation<br />
T<br />
5<br />
et d’une pression P 2<br />
à l’entrée à une pression<br />
d’évaporation P<br />
1<br />
à sa sortie. Pour ce composant, les<br />
inerties thermiques sont négligées [6].<br />
m<br />
⎛ ⎛ ⎞ ⎞<br />
⎜ ⎜<br />
p<br />
= S P − A − B ⎟ P ⎟<br />
⎜<br />
⎝ ⎝ ⎠ ⎠<br />
1<br />
&<br />
d<br />
2 1<br />
ρ<br />
⎜ ⎟<br />
4 (17)<br />
P ⎟<br />
cr<br />
RESULTATS DE SIMULATION<br />
Les résultats représentés sur les figures (5, 6, 7 et 8)<br />
donnent l’évolution des paramètres de fonctionnement<br />
pour trois démarrages de l’installation à des ouvertures<br />
différentes de la vanne de détente. La température de l’<strong>air</strong><br />
est de 10 °C, sa vitesse est de 3 m/s, la température de<br />
l’<strong>eau</strong> est de 20 ; 30 et 40 °C et son débit massique est de<br />
0,15 kg/s. Les résulats représentés par des points sont<br />
ceux issus des esssais et ceux qui sont représentés par des<br />
traits continus sont issus des calculs.<br />
Tanger, Maroc du 15 au 17 Novembre 2005 83
12èmes Journées Internationales de Thermique<br />
------------------------------------------------------------------------------------------------------------------------------------------------<br />
120<br />
la température de l’<strong>eau</strong> dans l’échangeur haute pression<br />
pour les 25 nœuds de calcul.<br />
P (bar)<br />
100<br />
80<br />
60<br />
34<br />
32<br />
1 2 3 4 5 6 7 8 9<br />
10 11 12 13 14 15 16 17 18<br />
19 20 21 22 23 24 25<br />
40<br />
20<br />
0<br />
P2 : 10, 20 °C P1 : 10, 20 °C P2 : 10, 30 °C<br />
p1 : 10, 30 °C P2 : 10, 40 °C P1 : 10, 40 °C<br />
0 20 40 60 80 100<br />
temps (s)<br />
Figure 5 : Variation des pressions expérimentales et<br />
calculées au cours du démarrage<br />
débit (kg/s)<br />
0.06<br />
0.05<br />
0.04<br />
0.03<br />
0.02<br />
92 bar, 40 °C<br />
0.01<br />
120 bar, 30 °C<br />
76 bar, 20 °C<br />
0<br />
0.00 20.00 40.00 60.00 80.00 100.00<br />
temps (s)<br />
Figure 6 : Variation du débit massique dans le détendeur<br />
au cours du démarrage<br />
110<br />
T (°C)<br />
30<br />
28<br />
26<br />
24<br />
22<br />
20<br />
0 20 40 60 80 100<br />
Temps (s)<br />
Figure 9 : Variation de la température de l’<strong>eau</strong> dans<br />
l’échangeur haute pression<br />
CONCLUSIONS<br />
Les résultats de ce modèle sont satisfaisants en régime<br />
stabilisé. Si le régime permanent s’établit plus rapidement<br />
selon les calculs par rapport aux essais c’est parce que les<br />
inerties de la bouteille de réserve de l’échangeur<br />
intermédi<strong>air</strong>e et des tubes qui relient différents éléments<br />
de l’installation ne sont pas prises en compte. Un<br />
développement futur de ce travail consisterait à prendre<br />
en compte ces paramètres.<br />
REFERENCES BIBLIOGRAPHIQUES<br />
[1] GUITARI I., HABERSCHILL P., LALLEMAND A.<br />
« Étude expérimentale et modélisation d’une <strong>pompe</strong> à<br />
<strong>chaleur</strong> fonctionnant avec du CO 2 ». Congrès français de<br />
Thermique, SFT 2003, Grenoble, 857-862, 3-6, (2003).<br />
T (°C)<br />
90<br />
70<br />
50<br />
T2<br />
T3<br />
T4<br />
[2] OLSON D. A. « Heat transfer of supercritical carbon<br />
dioxide flowing in a cooled horizontal tube ». IIF-IIR,<br />
Commission B1, B2 with E1, E2, Purdue University,<br />
USA, 261-269, (2000).<br />
30<br />
10<br />
0 20 40 60 80 100<br />
temps (s)<br />
Figure 7 : Variation des températures dans la partie<br />
haute pression au cours du démarrage<br />
Pression (bar)<br />
140<br />
120<br />
100<br />
80<br />
60<br />
40<br />
20<br />
-10°C<br />
0°C<br />
10°C<br />
20°C<br />
30°C<br />
40°C<br />
50°C<br />
Tee=20 °C<br />
Tee=30 °C<br />
Tee=40 °C<br />
-50 0 50 100 150 200 250 300 350 400<br />
Enthalpie (kJ/ kg)<br />
70°C<br />
90°C<br />
110°C<br />
130°C<br />
150°C<br />
Figure 8 : Cycles de fonctionnement simulé à t=100 s<br />
À partir des ces résultats on note que les variables<br />
expérimentales et calculées évoluent vers des valeurs très<br />
proches et que les résultats issus des calculs évoluent plus<br />
rapidement que ceux issus des essais. La figure 8<br />
représente les cycles de fonctionnement à t =100 s issus<br />
des calculs pour un fonctionnement stabilisé pour les trois<br />
cas de la Figure 5. La Figure 9 représente la variation de<br />
[3] JANG J., HRNJAK P. S. « Condensation of carbon<br />
dioxide at low temperatures ». IIR-IIF, Commission B1,<br />
B2, E1 and E2, Guangzhou, China, 64-74, (2002).<br />
[4] HWANG Y., KIM B. H., RADERMACHER R. «<br />
Boiling heat transfer correlation for carbon dioxide ». IIF-<br />
IIR, Commission B1, with E1 and E2, College Park,<br />
USA, 81-88, (1997).<br />
[5] BOURKE P. J., PULLING D. J., GILL L. E.,<br />
DENTON W. H. « Forced convective heat transfer to<br />
turbulent CO 2 in supercritical region ». International<br />
Journal of Heat and Mass Transfer, Vol. 13, 1339-1348,<br />
(1970).<br />
[6] SKAUGEN G., SVENSSON M. C. « Dynamic<br />
modelling and simulation of a trans critical CO 2 heat<br />
pump unit ». IIF-IIR, Commission B2, with B1, E1 and<br />
E2, Oslo, Norway, 230-239, (1998).<br />
Tanger, Maroc du 15 au 17 Novembre 2005 84