FR 2932262 A1 20091211 FR 0803187 A 20080609 PROCEDE D'ESTIMATION DE LA PUISSANCE ABSORBEE PAR UN COMPRESSEUR AU DEMARIRAGE D'UNE INSTALLATION DE CLIMATISATION La présente invention concerne un procédé d'estimation de la puissance absorbée par un compresseur de fluide frigorigène d'un système de climatisation lors du démarrage de cette installation de climatisation. L'invention trouve une application particulièrement avantageuse dans le domaine de la climatisation à l'intérieur de l'habitacle des véhicules automobiles. Dans les véhicules à moteur classiques, le compresseur du circuit de climatisation est entraîné par le moteur et consomme donc une partie de la puissance fournie par ce dernier au détriment de la puissance disponible pour io le roulage du véhicule. Bien qu'elle ne soit pas très importante, la puissance absorbée par le compresseur lorsqu'il est en marche n'en réduit pas moins le rendement du moteur. Cette diminution de rendement a pour conséquence d'augmenter la consommation de carburant et la pollution générée par les gaz d'échappement du véhicule. 15 Pour optimiser le rendement du moteur, une solution consiste à estimer la puissance réellement absorbée par le compresseur. En effet, la connaissance de cette information permet à une unité centrale de gestion du véhicule d'adapter un certain nombre de paramètres aux besoins réels, notamment les paramètres d'injection du moteur, ainsi que ceux du groupe 20 moto-ventilateur (GMV) impliqué dans le refroidissement du fluide frigorigène à la sortie du compresseur. Dans ce contexte, la demande de brevet français n° 2 840 259 propose d'estimer la puissance WS,, absorbée par le compresseur à partir du débit massique tnrcp,. du fluide réfrigérant comprimé par le compresseur, ces deux 25 grandeurs étant reliées par une expression de la forme : Ws/, = a.mrcpr •W.s + b.N , où W;s est le travail de compression isentropique du compresseur, N la vitesse de rotation du compresseur, et a et b des paramètres caractéristiques du compresseur. En fait, dans la demande brevet précitée, le débit massique dans le compresseur n'est pas estimé en tant que tel mais approché par le débit massique rh,.eXp du fluide frigorigène dans l'organe de détente de l'installation de climatisation. La raison pour cela est que le débit massique m,.eXp dans l'organe de détente, outre des caractéristiques propres à l'organe de détente lui-même, est fonction de grandeurs thermodynamiques facilement mesurables qui sont la température Texpin du fluide frigorigène à l'entrée de l'organe de détente et les pressions Pexpiä el: Pexpo du fluide frigorigène respectivement à l'entrée et la sortie de l'organe de détente. La demande de io brevet français n° 2 840 259 fournit en détail les relations entre ces différents paramètres. Ce modèle, qui suppose l'égalité entre les deux débits massiques rh,.exp et rhrC,p. , est valable en régime stationnaire ou dans des régimes faiblement transitoires. Par contre, il n'est plus applicable lors de variations rapides du 15 débit du compresseur, en particulier au démarrage de l'installation de climatisation. Dans ce cas en effet, on observe d'une part, que la différence des pressions Pexpi,, et Pexpont au niveau de l'organe de détente est très faible, d'où un débit massique rh,.eXp quasi nul, voire nul, alors que, d'autre part, la variation de la pression dans le circuit haute pression, situé entre la sortie du 20 compresseur et l'entrée de l'organe de détente, augmente de manière brutale, ce qui indique un terme d'accumulation de masse non négligeable. Dans ces conditions, le débit massique mrcpr dia compresseur ne peut plus être assimilé au débit massique rrlrexp de l'organe de détente. On peut même écrire que riircp. est bien supérieur à 1h,.exp . L'estimation de la puissance absorbée par le 25 compresseur devient alors complètement erronée. Or, il y a avantage à pouvoir continuer à fournir une estimation de la puissance absorbée par le compresseur aussi précise que possible durant toute la phase de démarrage, même si elle est relativement courte. Cette avantage réside dans le fait qu'on évite tous à-coups sur le moteur consécutif 30 à la mise en marche et à l'arrêt intempestif du compresseur. Le confort est donc amélioré. Aussi, un but de l'invention est de proposer un procédé qui permettrait d'estimer avec une bonne précision la puissance absorbée par le compresseur au démarrage de l'installation de climatisation. Ce but est atteint, conformérnent à l'invention, grâce à un procédé s d'estimation de la puissance Pi/s, absorbée par un compresseur de fluide réfrigérant au démarrage d'une installation de climatisation, ladite installation comprenant un circuit dit haute pression situé entre la sortie du compresseur et l'entrée d'un organe de détente du fluide réfrigérant, remarquable en ce que ledit procédé comprend les étapes consistant, après démarrage de io l'installation de climatisation, à : (a) calculer la dérivée dHP/dt des variations de la haute pression HP du fluide réfrigérant dans le circuit haute pression, et en déduire une valeur 5 mr de la variation correspondante du débit massique du fluide frigorigène dans le circuit haute pression, 15 (b) estimer le débit massique mreXp du fluide réfrigérant dans ledit organe de détente, et en déduire le débit massique mrcpr du fluide réfrigérant dans le compresseur par la relation : mrcp,. = m,.eXp + 8 mr , (c) estimer la puissance absorbée yv,, par le compresseur à partir du débit massique mrcpr 20 Ainsi, comme on le verra en détail plus loin, le procédé conforme à l'invention fait usage d'un modèle dynamique faisant le bilan de la masse de fluide réfrigérant dans le circuit haute pression de l'installation de climatisation, de manière à obtenir une bonne estimation du débit massique mrcpr du compresseur en corrigeant le débit massique mrexp de l'organe de détente de 25 la quantité 8 mr évaluée à partir des variations de la pression HP dans le circuit haute pression durant la phase de démarrage. Selon un mode de réalisation particulier, le procédé conforme à l'invention comprend une étape préalable consistant à comparer, après démarrage de l'installation de climatisation, la valeur HP(t) de la pression du 30 fluide réfrigérant dans le circuit haute pression à la valeur HP(0) au démarrage de l'installation, et à effectuer les étapes (a) à (c) à partir d'un instant initial ts pour lequel la différence HP(tr)-HP(ts -80 est supérieure à un seuil donné AHPo . De cette manière, on s'assure, avant de mettre en oeuvre les étapes (a) à (c) du procédé, que l'augmentation constatée de la pression HP est bien représentative du démarrage de l'installation de climatisation et non pas due à des phénomènes fortuits. L'invention prévoit en outre que les étapes (a) à (c) sont effectuées pendant une durée donnée To après ledit instant initial ts. Cette durée To est limitée à la période de temps pendant laquelle le modèle dynamique de bilan de masse mentionné plus haut est valable, c'est-à-dire quelques secondes en pratique. io Enfin, l'invention propose d'utiliser pour le fluide frigorigène une équation d'état de la forme P = p,..RT , où p,. est la masse volumique du fluide, de sorte que la variation 8 m, du débit massique est reliée à la dérivée dHP/dt par la relation : S rim, = A.dHP/dt , où A = VHP /RT , VHP et T étant le volume et la température du circuit haute pression, et R une constante du 15 fluide frigorigène. Il suffit donc déterminer et de stocker dans une base de données les valeurs du paramètre A pour obtenir très simplement une estimation de la correction 8 m,. cherchée. Il faut également signaler que l'invention s'applique tout aussi bien aux installations de climatisation comprenant un condenseur et un organe de 20 détente à orifice de section fixe qu'aux installations plus spécialement dédiées aux fluides réfrigérants supercritiques, comprenant un refroidisseur de gaz, un échangeur de chaleur interne et un organe de détente à orifice de section variable. La description qui va suivre en regard des dessins annexés, donnés à 25 titres d'exemples non limitatifs, fera bien comprendre en quoi consiste l'invention et comment elle peut être réalisée. La figure 1 est un schéma d'une installation de climatisation apte à mettre en oeuvre le procédé conforme à l'invention. La figure 2 est un organigrarnme des différentes étapes du procédé 30 conforme à l'invention. La figure 3 est un graphe donnant en fonction du temps les variations du couple du compresseur mesuré (A), estimé (B) par le procédé de l'art antérieur et estimé (C) par le procédé conforme à l'invention. Sur la figure 1 est représentée une installation de climatisation comprenant un circuit dit haute pression situé entre la sortie d'un compresseur 11 de fluide frigorigène et l'entrée d'un organe '12 de détente. On notera HP la pression régnant dans le circuit haute pression. La sortie de l'organe 12 de détente et l'entrée du compresseur 11 définissent un circuit dit basse pression de l'installation. La pression dans le circuit basse pression sera notée BP. Bien que cela ne soit pas essentiel pour l'invention, l'installation de la figure 1 concerne plus particulièrement un fluide réfrigérant supercritique comme le dioxyde de carbone CO2. Dans ce cas, l'organe 12 de détente est io généralement constitué par un orifice de section variable. Le fluide réfrigérant à la pression BP reçu du circuit basse pression sous forme gazeuse est porté à la haute pression HP par le compresseur 11. Dans le cas du dioxyde de carbone (CO2), la basse pression BP est par exemple de l'ordre de 30 à 40 bars, la haute pression HP pouvant atteindre 15 des valeurs de l'ordre de 90-133 bars. Un refroidisseur 13 de gaz (connu sous le terme anglo-saxon de � gas cooler �) permet de refroidir, à pression sensiblement constante HP, le fluide réfrigérant comprimé à l'état gazeux par le compresseur 11. Dans le mode de réalisation de la figure 1, la haute pression HP est mesurée à l'entrée du 20 refroidisseur 13 par un capteur 14. Le fluide frigorigène est refroidi dans le refroidisseur 13 par échange thermique avec un flux d'air extérieur 23 mis en mouvement par l'avancement du véhicule ou, éventuellement, au moyen d'un groupe moto-ventilateur (GMV) 24. L'organe 12 de détente a pour effet de faire passer la pression du fluide 25 frigorigène provenant du refroidisseur 13 de la haute pression HP à la basse pression BP en l'amenant, au moins en partie, à l'état liquide. Cette détente s'accompagne d'une diminution de la température du fluide frigorigène. L'air 21 de climatisation destiné à être envoyé vers l'habitacle du véhicule est refroidi par échange thermique dans un évaporateur 15 avec le 30 fluide réfrigérant provenant de l'organe 12 de détente. Le flux d'air 21 de climatisation est produit par un pulseur 22 alimenté par l'air extérieur ou en partie venant de l'habitacle (recyclage partiel ou total). Dans l'évaporateur 15, le fluide frigorigène passe, à pression sensiblement constante BP, de l'état liquide à un état au moins partiellement gazeux. Un accumulateur 16 disposé en sortie de l'évaporateur 15 permet la séparation du liquide. On peut voir sur la figure 1 la présence dans l'installation de climatisation d'un échangeur thermique interne 17 permettant au fluide frigorigène à haute température circulant du refroidisseur 13 vers l'organe 12 de détente de céder de la chaleur au fluide à plus basse température circulant de l'évaporateur 15 vers le compresseur 11. L'échange de chaleur ainsi réalisé présente l'avantage d'améliorer le rendement de la boucle de climatisation effectuée par le fluide réfrigérant. Il faut cependant noter que io l'utilisation d'un échangeur thermique interne n'est pas indispensable à l'invention, même si elle présente des avantages notamment dans le cas de fluides réfrigérants supercritiques. Comme cela a été mentionné plus haut, l'estimation de la puissance WS,, absorbée par le compresseur 11 nécessite de connaître l'information de 15 débit massique iiircp. du fluide réfrigérant dans le compresseur 11. En régime permanent ou faiblement transitoire, ce débit massique mrcpr peut être approché avec une bonne précision par le débit massique mrexp du fluide dans l'organe 12 de détente. Cependant, au démarrage de l'installation de climatisation, la différence 20 entre les pressions Pm,. et Pespon[ en entrée el: en sortie de l'organe 12 de détente est très faible, voire nulle, de sorte que le débit massique mrexp est également très faible, voire nul, puisque proportionnel à .JPexpi,t ù Pexpord , alors que, par ailleurs et dans le même temps, le débit massique mrcpr du compresseur 11 dans le circuit haute pression est non négligeable. Il est donc 25 impératif d'apporter une correction 8 m,. au débit massique mrexp de manière à pouvoir estimer le débit massique li'rcp,. cherché par la relation : mrcpr ùmrexp + 8 m,. . Pour cela, on utilise un modèle dynamique de bilan de masse qui va maintenant être présenté. 30 Ce modèle repose sur l'idée qu'au démarrage de l'installation où le débit massique rmrcp, du compresseur 11 est non nul et le débit massique mrexp de l'organe 12 de détente quasi nul, il se produit dans le circuit haute pression une accumulation de masse engendrant un variation temporelle rapide de la haute pression HP. Cette accumulation de masse étant directement reliée à la différence 8 th,. entre les débits rim,.cp,. et mrexp , il existe une relation 8 m,. = f(dHP/dt) entre les variations temporelles du débit massique différentiel 8 mr et la dérivée par rapport au temps t de la haute pression HP. II suffit donc, après avoir établi une loi f d'accumulation de masse, de mesurer la pression HP au cours du temps et d'en calculer la dérivée pour obtenir une estimation de la correction 6 Fhr à appliquer au débit massique rh,.eXp et en déduire ainsi le débit massique m,.cp,. cherché pour le compresseur 11. io A titre d'exemple, on peut prendre pour le fluide réfrigérant l'équation d'état suivante : P= p,..RT où pr est la masse volumique du fluide et R une constante caractéristique du fluide qui dépend de la pression P et de la température T . Si on appelle VHP le volume du circuit haute pression, on peut écrire : 15 8 mr = 1*rcp,. mrexp = VHP.dpr /dt . En utilisant l'équation d'état ci-dessus, on obtient : 6 mr = (VHP/RT).dHP/dt qui représente la loi d'accumulation de masse dans le circuit haute pression. On posera A = VHP/RT . La haute pression HP étant mesurée par le capteur 14 à l'entrée du 20 refroidisseur 13 de gaz et la température T étant prise par exemple égale à la température ambiante, on en déduit les valeurs cle la dérivée dHP/dt et donc celles de 8 rh,., le volume VHP étant connu de même que la constante R . Le tableau ci-dessous donne un certain nombre de valeurs typiques des différents paramètres pour le dioxyde de carbone (CO2) comme fluide 25 réfrigérant. HP (bar) T (°C) p, (kg/m3) RT (Pa/kg/ m3) VHP (dm3) A (kg/Pa) 40 30 90 44602 0,35 7,8.10-9 50 30 134 40367 0.35 8,7.109 30 60 35 158 37859 0,35 9,2.10-9 70 40 198 35428 0,35 9,9.109 80 40 278 28743 0,35 1,2.10-9 90 45 339 26562 0,35 1,3.10-9 L'organigramme de la figure 2 montre comment le procédé conforme à l'invention est mis en oeuvre pratiquernent. A l'instant t=0 où l'installation de climatisation démarre, le capteur 14 5 de pression est mis en en fonctionnement et fournit des valeurs HP(t) de la haute pression au cours du temps selon un pas de mesure 8 t . Dans une étape préalable au procédé proprement dit, les valeurs HP(t) de la haute pression sont comparées à la valeur de démarrage HP(0) jusqu' à ce que, à un instant donné noté ts , la différence HP(t)ûHP(t s ûSt) io dépasse un seuil AHP0 donné. La valeur du seuil AHPo est fixée de manière qu'au-delà de cette valeur les variations de la haute pression HP puissent être considérées comme significatives, c'est-à-dire effectivement produites par un débit du compresseur 11. En pratique, pour une valeur du seuil AHPo de l'ordre de 2 bar si St est d'environ 0.1 seconde. 15 A partir de l'instant ts , la dérivée dHP/dt est calculée à chaque pas 8 t de mesure par la formule (HP(t + 8 t) -- HP(t))/8 t . On en déduit la valeur de la correction 8 rihr en appliquant la loi d'accumulation de masse établie plus haut. Enfin, le débit massique m,.cpr cherché est donné par mrcprùmrexp+8 mr, le débit massique m,.exp dans l'organe 12 de détente étant estimé par toute 20 technique connue, notamment celle décrite dans la demande de brevet français n° 2 840 259. Finalement, la puissance YVS,, absorbée par le compresseur 11 peut être estimée à chaque pas de mesure S t par une relation du type WS,,=a.rim,.cp,.W,S +b.N, où 4Ç est le travail de compression isentropique du 25 compresseur, N la vitesse de rotation du compresseur, et a et b des paramètres caractéristiques du compresseur. Le procédé se déroule ainsi pendant une durée To au-delà de laquelle le modèle dynamique de bilan de masse n'est plus applicable du fait que l'organe 12 de détente débite du fluide réfrigérant, mettant fin au régime 30 d'accumulation de masse dans le circuit haute pression. La durée To est de l'ordre de quelques secondes, 3 secondes par exemple. La demanderesse a testé sur un exemple réel le procédé d'estimation qui vient d'être décrit et a obtenu les résultats résumés sur le graphe de la figure 3 qui donne les variations au cours du temps, après l'instant t=0 de démarrage de l'installation de climatisation, du couple fourni par le compresseur 11, la puissance absorbée étant directement liée au couple développé. La courbe (A) correspond au couple mesuré, la courbe (B) au couple estimé selon un modèle classique de régime stationnaire et la courbe (C) au couple estimé selon le procédé conforme à l'invention prenant en compte le phénomène d'accumulation de masse au démarrage. On voit clairement sur ce graphe que le procédé de l'invention conduit à une réduction significative io de l'erreur introduite par le modèle utilisé dans l'art antérieur puisque la courbe (C) illustrant l'invention est plus proche de la courbe mesuré (A), c'est-à-dire l'état réel du compresseur, que la courbe (B). 30 REVENDICATIONS1. Procédé d'estimation de la puissance YVs,, absorbée par un compresseur (11) de fluide réfrigérant au démarrage d'une installation de climatisation, ladite installation comprenant un circuit dit haute pression situé entre la sortie s du compresseur (11) et l'entrée d'un organe (12) de détente du fluide réfrigérant, caractérisé en ce que ledit procédé comprend les étapes consistant, après démarrage de l'installation de climatisation, à : (a) calculer la dérivée dHP/dt des variations de la pression HP du fluide réfrigérant dans le circuit haute pression, et en cléduire une valeur 8 m,. de la io variation correspondante du débit massique du fluide réfrigérant dans le circuit haute pression, (b) estimer le débit massique m,.exp du fluide réfrigérant dans ledit organe (12) de détente, et en déduire le débit massique du fluide réfrigérant dans le compresseur (11) par la relation : rrt.p =mrexp+8 rn,., 15 (c) estimer la puissance absorbée Yi s,, par le compresseur (11) à partir du débit massique mrcpr 2. Procédé selon la revendication 1, comprenant une étape préalable consistant à comparer, après démarrage de l'installation de climatisation, la 20 valeur HP(t) de la pression du fluide réfrigérant dans le circuit haute pression à la valeur HP(t ù 6t) au démarrage de l'installation, et à effectuer les étapes (a) à (c) à partir d'un instant initial ts pour lequel la différence HP(t)-HP(t s ù 5t) est supérieure à un seuil donné AH1 ô . 25 3. Procédé selon la revendication 2, dans lequel les étapes (a) à (c) sont effectuées pendant une durée donnée To après ledit instant initial ts . 4. Procédé selon la revendication 3, dans lequel les étapes (a) à (c) sont effectuées pendant la durée To par pas temporel 8 t de mesure. 5. Procédé selon l'une quelconque des revendications 1 à 4, dans lequel la variation 8 riz, du débit massique est reliée à la dérivée dHP/dt par la 2932262 Il relation : 8 m,. = A. dHP/dt , où A = VHP /RT , VHP et T étant le volume et la température du circuit haute pression, et R une constante du fluide réfrigérant. 6. Procédé selon l'une quelconque des revendications 1 à 5, dans lequel la puissance Yi/s,, absorbée par le compresseur (11) est reliée au débit massique mrcpr par la relation : YVS, = a.m .~p..WS + b.N , où W,s est le travail de compression isentropique du compresseur, N la vitesse de rotation du compresseur, et a et b des paramètres caractéristiques du compresseur. io 7. Installation de climatisation mettant en oeuvre le procédé d'estimation de la puissance WS,, absorbée par un compresseur selon l'une quelconque des revendications précédentes.