12.07.2015 Views

Etude de l'influence du fluide caloporteur sur les dimensions des ...

Etude de l'influence du fluide caloporteur sur les dimensions des ...

Etude de l'influence du fluide caloporteur sur les dimensions des ...

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

∇P= − UkAprès intégration <strong>de</strong> cette équation on obtient :µlepξp,évap=k ρ π n Dlµ laLévap(6)(9)2.2. Le déroulement <strong>du</strong> calculDans ce calcul, outre la température <strong>de</strong> fonctionnement <strong>de</strong> la BFDPT et la puissancedissipée à l’évaporateur, <strong>de</strong> nombreux paramètres géométriques peuvent être choisis :perméabilité, porosité et épaisseur <strong>du</strong> poreux, nombre et diamètre <strong>de</strong>s artères dansl'évaporateur, longueurs <strong>de</strong> l'évaporateur et <strong>du</strong> con<strong>de</strong>nseur, fraction volumique <strong>de</strong> vapeur dansle con<strong>de</strong>nseur, différence <strong>de</strong> hauteur entre l'évaporateur et le con<strong>de</strong>nseur, longueurs <strong>de</strong> lacon<strong>du</strong>ite <strong>de</strong> liqui<strong>de</strong> et <strong>de</strong> la con<strong>du</strong>ite <strong>de</strong> vapeur et enfin coefficient <strong>de</strong> sécurité pour le calcul<strong>du</strong> volume. Les diamètres <strong>de</strong> la con<strong>du</strong>ite <strong>de</strong> liqui<strong>de</strong> ou <strong>de</strong> vapeur sont <strong>les</strong> paramètres <strong>du</strong>solveur. En effet c'est la taille <strong>de</strong>s con<strong>du</strong>ites et <strong>du</strong> con<strong>de</strong>nseur qui détermine le volume <strong>de</strong>flui<strong>de</strong> à l'intérieur <strong>de</strong> la boucle. Le diamètre <strong>du</strong> con<strong>de</strong>nseur est par défaut égal au diamètre <strong>de</strong>la ligne <strong>de</strong> vapeur. Le diamètre optimal pour la configuration étudiée est donc obtenu à l’issue<strong>du</strong> calcul. Le tableau 1 donne <strong>les</strong> caractéristiques géométriques <strong>de</strong> la BFDPT <strong>de</strong> cette étu<strong>de</strong> :Perméabilité k Porosité ε Rayon <strong>de</strong> pore Epaisseur <strong>du</strong> Nombre Diamètrer pores poreux e p d’artères d’artèrem -2 m m m5,0.10 -13 0,7 2,43.10 -6 8,0.10 -3 3 8,0.10 -3Longueur <strong>de</strong> la Longueur <strong>de</strong> Longueur <strong>du</strong> Différence <strong>de</strong> hauteur entreligne <strong>de</strong> l’évaporateur con<strong>de</strong>nseur évaporateur et con<strong>de</strong>nseurvapeurLongueur <strong>de</strong> laligne <strong>de</strong> liqui<strong>de</strong>m m m m m0,7 0,7 0,3 0,6 0,5Tableau 1 : Géométrie <strong>du</strong> problèmeLe volume <strong>de</strong> la boucle (réservoir non inclus) et la masse <strong>de</strong> flui<strong>de</strong> dans la BFDPT sontalors calculés pour six flui<strong>de</strong>s différents : ammoniac, éthanol, méthanol, acétone, hexane etpentane. Les propriétés thermodynamiques <strong>de</strong>s flui<strong>de</strong>s étudiés sont également disponib<strong>les</strong>.L’eau n’est pas adaptée à une utilisation automobile puisqu’elle gèle à 0°C et c’est pourquoielle n’a pas été retenue. Par contre ce choix <strong>de</strong> flui<strong>de</strong>s ne tient pas compte <strong>de</strong>s normes envigueur dans le milieu automobile en matière <strong>de</strong> toxicité. Un calcul préliminaire <strong>du</strong> volume <strong>de</strong>flui<strong>de</strong> est effectué pour plusieurs valeurs <strong>du</strong> diamètre <strong>de</strong> la con<strong>du</strong>ite <strong>de</strong> liqui<strong>de</strong> et pour chaquetempérature <strong>de</strong> l’intervalle [-40°C ; 100°C], plage <strong>de</strong> température classique <strong>du</strong> milieuautomobile. Le diamètre choisi en entrée <strong>du</strong> calcul principal est le plus petit diamètre <strong>de</strong>con<strong>du</strong>ite <strong>de</strong> liqui<strong>de</strong> qui permet <strong>de</strong> vérifier l’équation (3) <strong>sur</strong> l’ensemble <strong>de</strong> cette plage <strong>de</strong>température. Le calcul principal effectué à l’ai<strong>de</strong> <strong>du</strong> solveur Excel fournit ensuite la valeuroptimisée <strong>du</strong> diamètre <strong>de</strong> la con<strong>du</strong>ite <strong>de</strong> vapeur. Tous <strong>les</strong> paramètres nécessaires au calcul <strong>du</strong>volume <strong>de</strong>s con<strong>du</strong>ites et <strong>du</strong> con<strong>de</strong>nseur sont alors disponib<strong>les</strong>.3. Résultats et analyse3.1. Puissances à transférer <strong>de</strong> 400 et 1800 WLa figure 2 illustre la masse <strong>de</strong> flui<strong>de</strong> nécessaire en fonction <strong>de</strong> la température pour <strong>les</strong> sixflui<strong>de</strong>s étudiés et une puissance injectée <strong>de</strong> 1800 W. Cette masse découle directement <strong>du</strong>


volume <strong>de</strong> flui<strong>de</strong> calculé à partir <strong>de</strong>s diamètres optimisés décrits au paragraphe 2.2. Le flui<strong>de</strong>le plus intéressant dans notre contexte automobile sera celui qui pour une même puissancetransférée présentera le plus faible encombrement et la plus petite masse. Dans la suite <strong>de</strong> cetarticle nous exprimerons cette idée par <strong>les</strong> termes « flui<strong>de</strong> performant ». On voit dans lafigure 2 que la masse d'ammoniac nécessaire au bon fonctionnement <strong>de</strong> la boucle diphasiquereste quasi-constante avec la température au regard <strong>de</strong>s autres flui<strong>de</strong>s. Les autres flui<strong>de</strong>sprésentent une augmentation significative <strong>du</strong> volume (et donc <strong>de</strong> la masse) nécessaire quandla température diminue. Cette tendance est particulièrement prononcée pour <strong>les</strong> alcools :éthanol et méthanol. A basse température (-20°C par exemple) <strong>les</strong> flui<strong>de</strong>s sont ici touscompris entre 50 g (ammoniac) et 250 g (éthanol). Cependant le volume et la masse <strong>de</strong>sflui<strong>de</strong>s ont augmenté avec la puissance dans <strong>de</strong>s proportions différentes pour chaque flui<strong>de</strong>.La figure 3 représente γ, la variation relative <strong>de</strong> volume <strong>de</strong> flui<strong>de</strong> en fonction <strong>de</strong> la puissance<strong>de</strong> dimensionnement. Elle se définit ainsi :γ =V ( 400)Cette gran<strong>de</strong>ur montre que le dimensionnement d’une boucle à ammoniac est beaucoupmoins sensible à la puissance <strong>de</strong>mandée et à la température considérée par rapport aux autresflui<strong>de</strong>s. Dans la plage <strong>de</strong> basse température (-30°C/0°C) la variation <strong>de</strong> volume est <strong>de</strong> 20%pour l'ammoniac et est comprise entre 70% et 140% pour <strong>les</strong> autres flui<strong>de</strong>s.3.2. Comportement à basse température3.2.1. Nombre <strong>de</strong> mérite <strong>de</strong>s flui<strong>de</strong>s étudiés( 1800) ( 400)V −VDunbar et Ca<strong>de</strong>ll [4] ont défini un nombre <strong>de</strong> mérite (12) à partir <strong>de</strong> l’équation d’équilibre(3) pour un flui<strong>de</strong> <strong>de</strong> travail utilisé dans une BFDPT qui n’est pas soumise à la gravitéterrestre :1,75ρvσ hlvNm=0,25(12)µLe calcul <strong>de</strong> ce nombre est effectué grâce à l’hypothèse stipulant que <strong>les</strong> pertes <strong>de</strong> chargeliées au liqui<strong>de</strong> sont négligeab<strong>les</strong> par rapport à cel<strong>les</strong> <strong>du</strong>es à la vapeur dans une boucle. Lafigure 4 représente N m pour <strong>les</strong> flui<strong>de</strong>s <strong>de</strong> cette étu<strong>de</strong>. La supériorité <strong>de</strong> l’ammoniac en termes<strong>de</strong> performance par rapport aux autres flui<strong>de</strong>s, déjà observée <strong>sur</strong> <strong>les</strong> figures 2 et 3, est iciconfirmée. Cette figure montre aussi une perte <strong>de</strong> performances <strong>de</strong>s alcools à bassetempérature plus appuyée que celle <strong>de</strong> l’acétone. On retrouve ici <strong>de</strong>s résultats cohérents avec<strong>les</strong> volumes <strong>de</strong> flui<strong>de</strong> calculés précé<strong>de</strong>mment. La figure 5 montre l’évolution <strong>de</strong> la masse <strong>de</strong>flui<strong>de</strong> nécessaire à 1800 W en fonction <strong>du</strong> nombre <strong>de</strong> mérite. Cette représentation permet <strong>de</strong>bien situer <strong>les</strong> flui<strong>de</strong>s qui offrent <strong>les</strong> meilleures performances. En effet il faut s’intéresser à lalocalisation <strong>de</strong> l’ensemble <strong>de</strong>s points concernant un flui<strong>de</strong> donné. La position <strong>de</strong> chaquecourbe ainsi obtenue est alors représentative <strong>de</strong> la performance <strong>du</strong> flui<strong>de</strong> <strong>sur</strong> tout l’intervalle<strong>de</strong> température étudié. Les points concernant l’ammoniac sont concentrés dans la zone <strong>de</strong>faible masse et <strong>de</strong> nombre <strong>de</strong> mérite élevé qui constitue la zone <strong>de</strong> performance optimale. Enrevanche la courbe <strong>de</strong> l’éthanol est très mal située par rapport aux autres (faible nombre <strong>de</strong>mérite et masse élevée). La courbe <strong>du</strong> pentane apparaît en décalage dans une zone <strong>de</strong> plusfaible masse que <strong>les</strong> autres flui<strong>de</strong>s. Ceci montre que <strong>les</strong> propriétés thermophysiques <strong>du</strong> liqui<strong>de</strong>ont aussi une influence <strong>sur</strong> nos résultats. En effet, le pentane se démarque <strong>de</strong>s autres flui<strong>de</strong>spar ses propriétés avantageuses (cf. §3.2.2) sachant que le calcul a été effectué à diamètre <strong>de</strong>con<strong>du</strong>ite <strong>de</strong> liqui<strong>de</strong> fixé <strong>sur</strong> tout l’intervalle <strong>de</strong> température. Une caractérisation <strong>de</strong>s flui<strong>de</strong>suniquement basée <strong>sur</strong> le nombre <strong>de</strong> mérite doit donc être complétée.v(11)


3.2.2. Influence <strong>de</strong>s propriétés <strong>sur</strong> <strong>les</strong> résultatsLa légère diminution <strong>du</strong> nombre <strong>de</strong> mérite <strong>de</strong> l’ammoniac aux fortes températures (T >40°C, figure 4) peut s’expliquer par l'évolution <strong>de</strong> sa tension <strong>de</strong> <strong>sur</strong>face et <strong>de</strong> sa chaleurlatente. En effet, el<strong>les</strong> décroissent environ <strong>de</strong>ux fois plus vite avec la température par rapportaux autres flui<strong>de</strong>s. Par contre, la masse nécessaire évolue peu au regard <strong>de</strong>s autres. En effet, laforte augmentation <strong>de</strong> la masse volumique <strong>de</strong> sa vapeur relativement aux autres flui<strong>de</strong>spermet <strong>de</strong> minimiser <strong>les</strong> pertes <strong>de</strong> charge générées et ainsi expliquer le fait que la masse <strong>de</strong>flui<strong>de</strong> à intro<strong>du</strong>ire dans la boucle n’augmente pas. Il reste donc le flui<strong>de</strong> dont l'utilisation estla plus avantageuse (même au <strong>de</strong>ssus <strong>de</strong> 20°C). Les pertes <strong>de</strong> charge d’un flui<strong>de</strong> dans laboucle dépen<strong>de</strong>nt <strong>de</strong> sa viscosité cinématique <strong>de</strong> la vapeur (figure 6). Plus cette viscosité estgran<strong>de</strong>, plus <strong>les</strong> pertes <strong>de</strong> charge dans la vapeur sont importantes, ce qui entraîne uneaugmentation <strong>du</strong> diamètre calculé (et donc <strong>du</strong> volume). En effet, à -40°C par exemple, on peutremarquer que la viscosité cinématique <strong>de</strong> la vapeur <strong>du</strong> pentane est cinq fois plus gran<strong>de</strong> quecelle <strong>de</strong> l'ammoniac et celle <strong>de</strong> l'hexane vingt fois plus gran<strong>de</strong>. Ceci nous apporte un élément<strong>de</strong> réponse pour expliquer la différence <strong>de</strong> performances observée entre ces <strong>de</strong>ux flui<strong>de</strong>s etl'ammoniac à basse température. Quant aux alcools, elle est beaucoup plus élevée commedans le cas <strong>du</strong> liqui<strong>de</strong>, ce qui <strong>les</strong> pénalise en termes d’efficacité. Les variations <strong>de</strong> la chaleurlatente <strong>de</strong> vaporisation et <strong>de</strong> la tension <strong>de</strong> <strong>sur</strong>face en fonction <strong>de</strong> la température ont égalementété observées. Les alcools ont globalement une meilleure chaleur latente que l’acétone,l’hexane et le pentane. A basse température l’ammoniac présente <strong>les</strong> meilleures chaleurlatente et tension <strong>de</strong> <strong>sur</strong>face. Cel<strong>les</strong>-ci diminuent ensuite considérablement et sa tension <strong>de</strong><strong>sur</strong>face est même la plus faible au-<strong>de</strong>là <strong>de</strong> 80°C, ce qui détériore son nombre <strong>de</strong> mérite alorsque sa masse n’augmente pas.Le nombre <strong>de</strong> mérite étant considéré pour un environnement 0-g, il est aussi nécessaire <strong>de</strong>s’intéresser à l’évolution <strong>de</strong> la masse volumique <strong>du</strong> liqui<strong>de</strong>. En effet, même si cette massevolumique varie peu avec la température au regard <strong>de</strong>s autres propriétés, la perte <strong>de</strong> chargegravitaire vient diminuer le saut <strong>de</strong> pression capillaire disponible pour équilibrer <strong>les</strong> pertes <strong>de</strong>charge dans la boucle (3). La masse volumique <strong>du</strong> liqui<strong>de</strong> (figure 7) est plus élevée pour <strong>les</strong>alcools. Or on a vu que cette propriété agit directement <strong>sur</strong> le saut <strong>de</strong> pression gravitaire quivient s’ajouter aux pertes <strong>de</strong> charges <strong>du</strong> système. L’augmentation <strong>de</strong> la masse volumique <strong>du</strong>liqui<strong>de</strong> à basse température est donc un facteur limitant pour <strong>les</strong> performances <strong>de</strong>s BFDPTutilisées dans <strong>de</strong>s applications terrestres. Une fois encore, l’ammoniac est le plus avantagé <strong>de</strong>tous <strong>les</strong> flui<strong>de</strong>s étudiés ici et sa masse volumique <strong>de</strong> liqui<strong>de</strong> constitue un atout supplémentairepour ses performances.4. ConclusionL’efficacité d’un flui<strong>de</strong> pour <strong>de</strong>s applications <strong>de</strong> type bouc<strong>les</strong> diphasiques est déterminéedans cette étu<strong>de</strong> par l’intermédiaire <strong>du</strong> volume à injecter dans la BFDPT pour que celle-cifonctionne à une puissance appliquée donnée. Une étu<strong>de</strong> complémentaire <strong>de</strong>s nombres <strong>de</strong>mérite <strong>de</strong>s flui<strong>de</strong>s au sens <strong>de</strong> Dunbar et Ca<strong>de</strong>ll [4] permet <strong>de</strong> mettre en évi<strong>de</strong>nce l’action <strong>de</strong>spropriétés thermophysiques suivantes <strong>sur</strong> <strong>les</strong> performances observées : tension <strong>de</strong> <strong>sur</strong>face,chaleur latente, viscosité cinématique et masse volumique <strong>de</strong> la vapeur. Les résultats obtenussont <strong>de</strong>s ordres <strong>de</strong> gran<strong>de</strong>ur : la géométrie <strong>de</strong> la boucle et l'expression <strong>de</strong>s pertes <strong>de</strong> chargesont simplifiées dans ce calcul. Ces ordres <strong>de</strong> gran<strong>de</strong>urs obtenus pour <strong>les</strong> volumes et masses<strong>de</strong> flui<strong>de</strong> dans la boucle sont très raisonnab<strong>les</strong> et ne remettent pas en cause l'utilisation d'uneboucle diphasique dans le cadre d’une application automobile.


Remerciements - Cette étu<strong>de</strong> entre dans le cadre <strong>du</strong> projet européen HI-CEPS : « Highly IntegratedCombustion Electric Propulsion System » pour le co-développement entre différents partenairesin<strong>du</strong>striels et universitaires d’un véhicule automobile hybri<strong>de</strong>. Cf. http://www.hi-ceps.eu/Références[1] Yu. F. Maydanik, Loop heat pipes, Applied Thermal Engineering, 25 (2005), 635-657.[2] S. Launay, V. Sartre, J. Bonjour, Parametric analysis of loop heat pipe operation: a literaturereview, International Journal of Thermal Science, 46 (2006), 621-636.[3] S. Launay, V. Sartre, J. Bonjour, Effect of fluid thermophysical properties on loop heat pipeoperation, 14 th International Heat Pipe Conference (Florianopolis, Brésil, 22-27 avril 2007).[4] N. Dunbar, P. Ca<strong>de</strong>ll, Working fluids and figure of merit for CPL/LHP applications, TheAerospace Corporation (CPL-98 Workshop, 1998) 1.3-1–1.3-6.[5] S. Bories, M. Prat, Transferts <strong>de</strong> chaleur dans <strong>les</strong> milieux poreux, Techniques <strong>de</strong> l’ingénieur, traitégénie énergétique, B 8 250.masse (g)500450400350300250200150100500-60 -40 -20 0 20 40 60 80 100 120température <strong>de</strong> fonctionnement (°C)ammoniacethanolméthanolacétonehexanepentaneFigure 2 : Masse <strong>de</strong> flui<strong>de</strong> à 1800Wg1,51,41,31,21,11,00,90,80,70,60,50,40,30,20,10,0-60 -40 -20 0 20 40 60 80 100 120température <strong>de</strong> fonctionnement (°C)ammoniacethanolméthanolacétonehexanepentaneFigure 3 : Variation relative <strong>de</strong> volumeNm1,0E+121,0E+111,0E+101,0E+09ammoniacethanolméthanol1,0E+08acétonehexanepentane1,0E+07-60 -40 -20 0 20 40 60 80 100 120température <strong>de</strong> fonctionnement (°C)Figure 4 : Nombre <strong>de</strong> méritemasse (g)5004504003503002502001501005001,00E+07 1,00E+08 1,00E+09 1,00E+10 1,00E+11 1,00E+12Nombre <strong>de</strong> mérite NmammoniacethanolméthanolacétonehexanepentaneFigure 5 : Masse <strong>de</strong> flui<strong>de</strong> et nombre <strong>de</strong> mériteviscosité (m2/s)1,E-021,E-031,E-041,E-051,E-06ammoniacethanolméthanolacétonehexanepentanemasse volumique (kg/m3)900850800750700650600550ammoniacethanolméthanolacétonehexanepentane5001,E-07-60 -40 -20 0 20 40 60 80 100 120température <strong>de</strong> fonctionnement (°C)Figure 6 : Viscosité cinématique vapeur450-60 -40 -20 0 20 40 60 80 100 120température <strong>de</strong> fonctionnement (°C)Figure 7 : Masse volumique liqui<strong>de</strong>

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!