FR 2852558 A1 20040924 FR 0303362 A 20030319 Installation de climatisation de véhicule pour fluide supercritique L'invention concerne les circuits de climatisation des véhicules à moteur. 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 du moteur. La puissance absorbée par le compresseur, quand il est en marche, diminue le rendement du moteur, et par suite, augmente la consommation de carburant et la pollution générée par les gaz d'échappement du véhicule. Cet inconvénient est en particulier gênant pour les compresseurs mécanique à contrôle externe, dont l'utilisation se généralise. Par ailleurs, dans des réalisations existantes, le calculateur d'injection du véhicule ne dispose pas de la valeur instantanée de la puissance réellement absorbée par le compresseur et choisit donc, pour le fonctionnement du compresseur, des paramètres d'injection par défaut correspondant à la valeur maximale de la puissance absorbée, valeur qui est rarement atteinte en pratique. Par conséquent, une solution pour optimiser le rendement du moteur consiste à estimer la valeur instantanée de cette puissance réellement absorbée par le compresseur. La connaissance de cette information peut alors permettre d'adapter les paramètres d'injection du moteur aux besoins réels. La demande de brevet français n 01 16568 propose une installation de climatisation propre à estimer le débit d'un fluide frigorigène sous-critique dans un circuit de climatisation pour calculer la puissance instantanée absorbée par le compresseur à partir de cette estimation. L'installation décrite dans la demande de brevet français n 0116568 comprend des organes de mesures permettant d'établir des valeurs relatives à la température du flux d'air extérieur à l'entrée du condenseur et à la pression du fluide frigorigène en sortie du compresseur. L'installation comprend en outre une carte électronique qui met en oeuvre une équation pour calculer une estimation du débit de fluide frigorigène dans le circuit de climatisation àpartir des valeurs de la température du flux d'air extérieur soufflé à l'entrée du condenseur, et de la haute pression. L'estimation du débit de fluide frigorigène proposée par cette installation s'adapte mal aux fluides supercritiques et ne permet pas d'estimer leur débit avec une précision satisfaisante. L'utilisation des fluides frigorigènes supercritiques, tel que le C02, s'est développée dans les circuits de climatisation des véhicules pour limiter les effets néfastes des fluides frigorigènes sur l'environnement. Les fluides frigorigènes supercritiques possèdent en effet un potentiel de réchauffement de la planète nettement inférieur à celui des fluides frigorigènes sous-critiques tels que les composés fluorés. Un circuit de climatisation utilisant un fluide supercritique comporte un compresseur, un refroidisseur de gaz ("gas cooler"), un échangeur thermique interne, un organe de détente et un évaporateur parcourus, dans cet ordre, par le fluide frigorigène. Dans un tel circuit, le refroidissement du fluide après compression n'entraîne pas de changement de phase. Le fluide ne passe à l'état liquide qu'au cours de la détente. Cette propriété des fluides supercritiques ne permet pas d'utiliser l'installation de la demande de brevet n 0116568 pour estimer le débit de fluide supercritique et la puissance consommée par le compresseur. C'est donc un but de l'invention de proposer une installation de climatisation permettant d'estimer le débit d'un fluide frigorigène supercritique, avec une précision satisfaisante. C'est encore un but de l'invention de fournir une telle estimation du débit de fluide frigorigène supercritique à partir de grandeurs représentatives du fonctionnement du refroidisseur de gaz. A cet effet, l'invention propose une installation de climatisation pour véhicule à moteur, munie d'un circuit de fluide frigorigène supercritique comprenant un compresseur, un refroidisseur de gaz de longueur totale donnée, un organe de détente et un évaporateur, le refroidisseur de gaz recevant un flux d'air passant par un groupe moto-ventilateur. L'installation comprend en outre un dispositif électronique de contrôle destiné à interagir avec le circuit de fluide frigorigène, pour estimer une grandeur relative au fluide frigorigène à partir de températures relatives au refroidisseur de gaz. Avantageusement, le dispositif électronique de contrôle comprend une fonction de calcul utilisant deux écarts de température, dont l'un au moins est fondé sur la température du fluide frigorigène en un point intermédiaire choisi du refroidisseur de gaz. Selon une caractéristique de l'invention, le point intermédiaire choisi est situé à une distance de l'entrée du refroidisseur de gaz, comprise entre 5 % et 35% de la longueur totale du refroidisseur de gaz. Selon une autre caractéristique de l'invention, les deux écarts de température sont en outre fondés sur : - la température du fluide frigorigène à l'entrée du refroidisseur de gaz, - la température du flux d'air reçu par le refroidisseur de gaz. Dans un mode de réalisation particulier, la grandeur relative au fluide frigorigène est le débit massique du fluide frigorigène. Selon ce mode de réalisation particulier, les deux écarts de températures utilisés par la fonction de calcul comprennent : - l'écart entre la température du fluide frigorigène à l'entrée du refroidisseur de gaz et la température du flux d'air reçu par le refroidisseur de gaz; - l'écart entre la température du fluide frigorigène au point intermédiaire choisi du refroidisseur de gaz et la température du flux d'air reçu par le refroidisseur de gaz. L'invention propose également un procédé d'estimation du débit massique d'un fluide frigorigène supercritique circulant dans un circuit de climatisation, ledit circuit comprenant un compresseur, un refroidisseur de gaz de longueur totale donnée, un organe de détente et un évaporateur, le refroidisseur de gaz recevant un flux d'air passant par un groupe moto-ventilateur. Avantageusement, le procédé comprend les étapes suivantes, pour une valeur initiale donnée du débit massique du fluide frigorigène : a) calculer - une estimation de l'écart entre la température du fluide frigorigène à l'entrée du refroidisseur de gaz et la température du flux d'air reçu par le refroidisseur de gaz; - une estimation de l'écart entre la température du fluide frigorigène au point intermédiaire choisi du refroidisseur de gaz et latempérature du flux d'air reçu par le refroidisseur de gaz b) calculer une grandeur de sortie, à partir d'une grandeur d'entrée et des valeurs fournies par l'étape a); c) déterminer si l'écart entre la grandeur de sortie calculée à l'étape b) et la grandeur d'entrée est inférieure à une valeur de seuil prédéfinie. En particulier, la grandeur d'entrée comprend la valeur initiale du débit de fluide frigorigène, à la première itération de l'étape a). Selon un aspect de l'invention, les étapes a) à c) sont réitérées avec une grandeur d'entrée égale à la grandeur de sortie calculée à l'étape c) précédente, si l'étape c) n'est pas vérifiée. Selon un autre aspect de l'invention, la grandeur de sortie comprend la valeur du débit massique de fluide frigorigène, si l'étape c) est vérifiée. L'invention couvre également un programme-produit, qui peut être défini comme comprenant les fonctions mises en oeuvre pour estimer le débit de fluide frigorigène. D'autres caractéristiques et avantages de l'invention apparaîtront à l'examen de la description détaillée ci-après, et des dessins annexés sur lesquels: - la figure la est un schéma d'un circuit de climatisation de véhicule automobile à moteur, fonctionnant avec un fluide supercritique, - la figure lb est un schéma d'une installation de climatisation de véhicule automobile à moteur, selon l'invention, - la figure 2 est un diagramme thermodynamique dans lequel est représentée l'évolution de la pression du fluide C02 en fonction de son enthalpie, durant un cycle, - la figure 3 est un schéma du refroidisseur de gaz, - la figure 4 représente l'évolution de la température le long du refroidisseur de gaz, - la figure 5 est un organigramme illustrant les différentes étapes mises en oeuvre par l'installation selon une première forme de réalisation selon l'invention, - les figures 6 à 9 représentent l'évolution de l'efficacité normale en fonction du débit de fluide C02, pour différentes valeurs du débit d'air reçu par le refroidisseur de gaz, selon l'invention, et - les figures 10 et 11représentent l'évolution de l'efficacité normale en fonction du débit de fluide C02, pour différentes valeurs du débit de l'air entrant dans le refroidisseur de gaz. Les annexes A, B et C comportent des équations mathématiques principales utilisées pour mettre en oeuvre l'installation. Les dessins contiennent, pour l'essentiel, des éléments de caractère certain. Ils pourront donc non seulement servir à mieux faire comprendre la description, mais aussi contribuer à la définition de l'invention, le cas échéant. En référence à la figure la, l'installation de climatisation selon l'invention comporte un circuit de climatisation parcouru par un fluide frigorigène supercritique. Par la suite, la description sera faite en référence au fluide frigorigène supercritique C02, à titre d'exemple non limitatif. La description du circuit de climatisation ci-après sera également faite en référence au diagramme de la figure 2 représentant l'évolution de la pression du fluide C02 en fonction de son enthalpie. Un circuit fonctionnant en mode supercritique comporte: - un compresseur 14 propre à recevoir le fluide à l'état gazeux à une pression PO et à le comprimer jusqu'à la pression P1, - un refroidisseur de gaz 11("gaz cooler") propre à refroidir le gaz comprimé par le compresseur, à pression P1sensiblement constante, - un organe de détente 12 propre à abaisser la pression du fluide issu du refroidisseur de gaz 11de la pression P1jusqu'à la pression PO, en l'amenant au moins en partie à l'état liquide, et - un évaporateur 13propre à faire passer le fluide à l'état liquide provenant de l'organe de détente à l'état gazeux, à pression sensiblement constante PO, pour produire un flux d'air climatisé 21 qui est envoyé vers l'habitacle du véhicule. Le circuit peut comporter en outre un échangeur thermique interne 23, permettant au fluide circulant du refroidisseur de gaz vers le détendeur de céder de la chaleur au fluide circulant de l'évaporateur vers le compresseur. Le circuit peut en outre comporter un accumulateur 17 placé entre la sortie de l'évaporateur et l'entrée du compresseur pour éviter les coups de liquide. Le refroidisseur de gaz 11reçoit un flux d'air extérieur 16 pour évacuer la chaleur prélevée dans l'habitacle, qui dans certaines conditions de fonctionnement est mis en mouvement par un groupe moto-ventilateur 15. L'évaporateur 13reçoit un flux d'air d'un pulseur alimenté par un flux d'air extérieur pour produire un flux d'air climatisé 21. L'organe de détente 12 peut être par exemple une vanne électrique ou un orifice calibré. Dans un circuit de climatisation fonctionnant en mode supercritique, le fluide frigorigène supercritique est refroidi par l'air avec lequel il échange de la chaleur, dans le refroidisseur de gaz 11, mais l'échange n'est pas effectué de manière isotherme, comme cela se produit pour un fluide sous-critique circulant dans un condenseur. Plus précisément, le fluide frigorigène supercritique est comprimé en phase gazeuse et amené à une pression élevée P1par le compresseur 14. Le refroidisseur de gaz 11refroidit ensuite le fluide frigorigène grâce au flux d'air entrant 16. A la différence des circuits de climatisation fonctionnant avec un fluide sous-critique, le refroidissement du fluide après compression n'entraîne pas de changement de phase. Le fluide ne passe à l'état liquide qu'au cours de la détente. L'échangeur thermique interne 23 permet de refroidir très fortement, voire de liquéfier le fluide sortant du refroidisseur de gaz. On se réfère maintenant à la figure lb qui représente une installation, selon un premier mode de réalisation de l'invention, mise en place dans un véhicule automobile, éventuellement en mouvement selon une vitesse d'avancement Va. Le véhicule automobile est animé par un moteur 43, qui peut être commandé par un calculateur d'injection 42. Le calculateur d'injection 42 reçoit des informations de divers capteurs qu'il interprète pour ajuster les paramètres d'injection. Le calculateur d'injection 42 peut en outre fournir des informations sur les conditions intérieures ou extérieures du véhicule (informations fournies par un capteur solaire, nombre d'occupants, etc. ). En particulier, il peut fournir des informations sur des valeurs instantanées relatives au fonctionnement du véhicule, et notamment sur la vitesse de rotation du compresseur N. Le véhicule est également équipé de l'appareil de climatisation 10 décrit ci-dessus, représentée schématiquement sur la figure lb. De plus, l'installation est munie d'un calculateur de climatisation 40, comprenant un régulateur d'habitacle 41 et un régulateur de boucle de climatisation 402. Le régulateur d'habitacle 41 est destiné à fixer la consigne de température de l'air extérieur soufflé à l'entrée de l'évaporateur 13. Le calculateur d'injection du moteur peut agir sur l'appareil de climatisation grâce à un régulateur de climatisation 402. Cette liaison peut interdire la mise en marche de l'appareil de la climatisation lorsque le moteur est fortement sollicité. La Demanderesse propose d'optimiser le fonctionnement d'une telle installation de climatisation fonctionnant en mode supercritique, en fournissant, avec une précision satisfaisante, une estimation d'une grandeur relative au fluide frigorigène, en particulier du débit massique de fluide frigorigène, à partir de températures relatives au refroidisseur de gaz. La demande de brevet n 01 16568 évoquée précédemment s'applique aux circuits de climatisation fonctionnant en mode sous-critique et propose une estimation du débit massique de fluide frigorigène à partir de : - l'écart =T -Ta, entre la température de condensation Tk et la température de l'air reçu par le condenseur Tai, - de la vitesse d'avancement du véhicule Va, et - de la tension du groupe moto-ventilateur Ugmv. Une solution analogue pour un circuit fonctionnant en mode supercritique reviendrait à estimer le débit massique d'un fluide frigorigène supercritique en fonction du pincement A=TgC0-Ta], qui représente l'écart entre la température Tgo en sortie du refroidisseur de gaz et la température de l'air reçu par ce dernier Tal. Or cette solution fournit des résultats peu précis, en raison des propriétés spécifiques des circuits fonctionnant en mode supercritique. En effet, ce manque de précision résulte du fait que l'échange de chaleur au niveau du refroidisseur de gaz n'est pas isotherme et en particulier au faible écart de température entre la sortie du refroidisseur de gaz et l'air de refroidissement. En outre, lorsque le véhicule est à l'arrêt, un retour d'air chaud peut se produire dans le refroidisseur de gaz, ce qui nuit encore à l'estimation. Ce faible écart de température entre la sortie du refroidisseur de gaz et l'air de refroidissement est illustré sur la figure 3 en relation avec la figure 4. La figure 3 est un schéma d'un refroidisseur de gaz sur lequel sont représentées: -un point x, à l'entrée du refroidisseur de gaz, au niveau duquel le fluide frigorigène est à la température Tgcl, - un point intermédiaire x, choisi entre l'entrée x, et la sortie Xo du refroidisseur de gaz, dont le positionnement sera détaillé ci-après. Au niveau de ce point, le fluide frigorigène est à la température Tg,x, et - un point Xo à la sortie du refroidisseur de gaz, au niveau duquel le fluide frigorigène est à la température Tgco. Le refroidisseur de gaz 11 a une longueur totale L entre son entrée x, et sa sortie xo. Cette longueur L correspond à la longueur développée du refroidisseur de gaz 11, c'est-à-dire à la somme des longueurs des différentes branches qui le composent. Par exemple, sur la figure 3, la longueur totale L est égale à la somme des longueurs des branches horizontales LI, L2, L3, L4, L5 et des branches de raccord verticales L,2, Lu, L34, et L45' Le flux d'air reçu par le refroidisseur de gaz 16 est également représenté sur la figure 3 avec une température Ta,. La figure 4 représente un exemple d'évolution de la température du fluide frigorigène C02 dans le refroidisseur de gaz 11en fonction de la position où est mesurée cette température sur le refroidisseur de gaz. Cette figure montre que l'écart entre la température Tgco du fluide frigorigène en sortie Xo du refroidisseur de gaz 11et la température Tal de l'air entrant 16 est faible. Par suite, une estimation du débit massique d'un fluide frigorigène à partir de cet écart A=TgC0-Tai, de manière analogue au brevet n 01 16568 précité, n'est pas suffisamment précise. L'invention propose une estimation du débit massique de fluide frigorigène fondée sur un bilan énergétique optimisé sur le refroidisseur de gaz. Le bilan énergétique au niveau du refroidisseur de gaz peut être optimisé lorsque les deux conditions suivantes sont réalisées : - l'écart (Tg,-TgX), lié au bilan sur le fluide C02, est maximal, et - l'écart (Tgcx-Tai), lié au bilan sur l'air, est maximal. Lorsque x=x,, l'écart (Tgcx-Ta,) est maximal, alors que l'écart (Tgcl-Tgcx) est minimal. Au contraire, lorsque x=x0, l'écart (TgCl-TgcJ est maximal alors que l'écart (Tgcx - Ta) est minimal. Il existe une valeur de x qui permet de vérifier les deux conditions à la fois. La Demanderesse a trouvé que ces deux conditions sont simultanément réalisées lorsque le point x est un point intermédiaire choisi du refroidisseur de gaz. Plus précisément, le bilan énergétique est optimal lorsque ce point est situé à une distance D de l'entrée x,, comprise entre 5% et 35% de la longueur totale L du refroidisseur de gaz. En particulier, la distance D=15% de la longueur L fournit un bilan énergétique optimal. L'installation de climatisation selon l'invention exploite cette optimisation du bilan énergétique pour estimer une grandeur relative au fluide frigorigène, en particulier le débit massique de fluide frigorigène mC02' Cette estimation est obtenue à partir de deux écarts de températures relatives au refroidisseur de gaz, dont l'un au moins est fondé sur la température du fluide frigorigène T gcx au point intermédiaire choisi x du refroidisseur de gaz. Les températures relatives au refroidisseur de gaz comprennent en outre : - la température Tgcl du fluide frigorigène, à l'entrée du refroidisseur de gaz, et - la température Tat de l'air 16 reçu par le refroidisseur de gaz. L'installation de climatisation comprend un dispositif électronique de contrôle, par exemple une carte électronique 401 destinée à interagir avec le circuit de climatisation 10, via les liaisons 30/31, et le calculateur d'injection 42, via les liaisons 32/33, pour estimer la grandeur relative au fluide frigorigène. La description ci-après sera faite en référence à l'estimation du débit massique de fluide frigorigène mC02, comme exemple non limitatif de grandeur relative au fluide frigorigène. La carte électronique 401 peut-être considérée comme partie intégrante du calculateur de climatisation 40 du véhicule. La carte électronique 401 comprend une fonction de calcul qui utilise deux écarts de températures relatives au fluide frigorigène, dont l'un au moins est fondé sur latempérature du fluide frigorigène Tgcx au point intermédiaire choisi x du refroidisseur de gaz. La carte électronique 401 peut récupérer des informations 30 provenant de capteurs mis en place sur le circuit de climatisation 10. Elle peut recevoir également des informations du calculateur d'injection du moteur 42 par la liaison 33, en particulier la vitesse d'avancement du véhicule Va, la tension du groupe moto-ventilateur Ugmv, la vitesse de rotation du compresseur N et la température de l'air extérieur Text. Dans une première forme de réalisation, la fonction de calcul de la carte électronique utilise les deux écarts de températures suivants : - l'écart entre la température du fluide frigorigène à l'entrée du refroidisseur de gaz Tgcl et la température du flux d'air reçu par le refroidisseur de gaz Ta,; -l'écart entre la température du fluide frigorigène Tgcx au point intermédiaire choisi du refroidisseur de gaz x et la température du flux d'air reçu par le refroidisseur de gaz Tal. La figure 5 représente les différentes étapes mises en oeuvre par la carte électronique pour estimer le débit de fluide frigorigène mC02' selon la première forme de réalisation de l'invention. Selon cette première forme de réalisation, le dispositif électronique de contrôle 401 est agencé pour appeler itérativement la fonction de calcul, à partir d'une valeur initiale donnée du débit massique du fluide frigorigène mo2 ,,t,Jusqu'à satisfaire un critère de fin d'itération, comprenant une condition de seuil sur l'écart entre la grandeur d'entrée mco2 in et la grandeur de sortie de la fonction de calcul mC02~out> Ainsi, à l'étape 200, l'installation estime ou mesure les valeurs de la température du fluide frigorigène à l'entrée du refroidisseur de gaz T,,,, de la température Ta, de l'air reçu par le refroidisseur de gaz 11 et de la température du fluide frigorigène T gcx au point intermédiaire choisi du refroidisseur de gaz. Pour cela, l'installation comprend des organes de mesures propres à fournir ces grandeurs. Plus précisément, la température du fluide frigorigène T,,,, à l'entrée du refroidisseur de gaz 11peut être directement mesurée par une sonde 130 placée à l'entrée du refroidisseur de gaz. En variante, la carte électronique peut comporter une fonction d'estimation de température de fluide capable de déterminer la température du fluide frigorigène Tgcl à l'entrée du refroidisseur de gaz lia partir de la pression PO et de la température T su en entrée du compresseur et de la pression en sortie du compresseur Pl. La température Tsu peut être mesurée par la sonde 180 de la figure lb placée en entrée du compresseur. En variante, la température Tsu peut également être estimée. Cette estimation sera détaillée lors de la description de la deuxième forme de réalisation de l'invention. La température Tal de l'air reçu par le refroidisseur de gaz 11peut être mesurée par une sonde de température 120, en référence à la figure 1b, placée entre le groupe motoventilateur 15 et le refroidisseur de gaz 11. En variante, la carte électronique peut comporter une fonction d'estimation de température d'air capable de déterminer la température Ta, de l'air reçu par le refroidisseur de gaz à partir de la vitesse d'avancement Va du véhicule et de la température de l'air extérieur Text en entrée du groupe moto-ventilateur. La vitesse d'avancement Va du véhicule et de la température de l'air extérieur Text en entrée du groupe moto-ventilateur peuvent être fournies par le calculateur d'injection du moteur 42 et par le régulateur d'habitacle 41. Selon cette variante, l'installation peut comprendre une sonde de température qui fournit une mesure de la température T ext de l'air extérieur au véhicule. La température du fluide frigorigène T gcx peut être directement mesurée par une sonde 110 placée au point x choisi du refroidisseur de gaz. A l'étape 202, la carte électronique 401 choisit une valeur initiale du débit de fluide frigorigène mC02~'mt' par exemple mCO2 init= 100kg/h. Le choix de cette valeur peut être fondé sur des estimations ou des mesures antérieures du débit de fluide frigorigène. La carte électronique appelle ensuite une première itération de la fonction de calcul avec une grandeur d'entrée mC02~ln égale à la valeur initiale mCO2 imt. La fonction de calcul effectue les étapes 204 à 214 décrites ci-après. Les étapes 204 à 208 fournissent une estimation de la puissance Pk(x) échangée entre le fluide C02 circulant dans la portion du refroidisseur de gaz, comprise entre x, et x, et l'air reçu à partir de: - la grandeur d'entrée mC02~ln; - l'écart entre la température du fluide frigorigène à l'entrée du refroidisseur de gaz Tgcl et la température du flux d'air reçu par le refroidisseur de gaz Tal; - l'écart entre la température du fluide frigorigène Tgcx au point intermédiaire choisi du refroidisseur de gaz x et la température du flux d'air reçu par le refroidisseur de gaz Tal. Plus précisément, à l'étape 204, la fonction de calcul calcule le coefficient d'échange thermique global entre l'air et le fluide frigorigène supercritique ksg à partir du coefficient d'échange thermique côté air ksa et du coefficient d'échange thermique côté C02, ksC02' selon l'équation A3. 1 de l'annexe A. Le coefficient d'échange thermique côté air ksa est une fonction du débit d'air refroidissant la portion du refroidisseur de gaz comprise entre l'entrée x, et le point intermédiaire x, et donc de la tension du groupe moto-ventilateur UGMV et de la vitesse d'avancement du véhicule Va, selon la relation A3.21 de l'annexe A. Un exemple d'équation A3.22 est donnée dans l'annexe A, où a, b, c, d et e sont des constantes. Le coefficient d'échange thermique côté C02, kSC02, est une fonction du débit de C02, selon la relation A3.31 de l'annexe A. Un exemple d'équation A3.32 est donné dans l'annexe A, où a' et b' sont des constantes qui dépendent de la technologie du refroidisseur de gaz et de ses caractéristiques. La grandeur d'entrée de la fonction de calcul mC02~ln peut donc être utilisée pour calculer le coefficient ksC02' A l'étape 206, la fonction de calcul détermine la moyenne logarithmique #TLM de la différence de température globale, entre l'entrée du refroidisseur de gaz x, et le point x, à partir : - de l'écart entre la température du fluide frigorigène à l'entrée du refroidisseur de gaz Tgcl et la température du flux d'air reçu par le refroidisseur de gaz Tal; - de l'écart entre la température du fluide frigorigène Tgcx au point intermédiaire choisi du refroidisseur de gaz x et la température du flux d'air reçu par le refroidisseur de gaz Tal. Ces écarts sont estimés à partir des températures Tgcl, Tgex et Ta, obtenues à l'étape initiale 200. Le calcul de la moyenne logarithmique #TLM en fonction des écarts de température précités est établi conformément à la relation A5 de l'annexe A. A l'étape 208, la fonction de calcul détermine la puissance Pk(x) échangée entre le fluide C02 circulant dans la portion du refroidisseur de gaz, comprise entre x, et x, et l'air reçu. La puissance Pk(x) est calculée à partir de la valeur de moyenne logarithmique #TLM obtenue à l'étape 206, du coefficient d'échange thermique global ksg obtenu à l'étape 204, conformément à l'équation A2 de l'annexe A. La fonction de calcul utilise également le débit d'air mai qui dépend de la tension du groupe moto-ventilateur Ug,r, et de la vitesse Va d'avancement du véhicule, ainsi que la capacité calorifique de l'air Cpa. À l'étape 210, la fonction de calcul détermine la valeur de la différence d'enthalpie Ahk(x), entre l'entrée x, du refroidisseur de gaz et le point x. Comme exprimé par la relation A4 de l'annexe A, cette valeur est calculée à partir de la température à l'entrée du refroidisseur de gaz T,,,, de la température au point x du refroidisseur de gaz Tacx et de la haute pression PI en sortie du compresseur. La haute pression P 1 peut par exemple être mesurée par le capteur 160 de la figure 1 b. Ce capteur 160 qui mesure la valeur instantanée de la haute pression P 1 peut-être placé à tout endroit approprié entre la sortie du compresseur et l'entrée du refroidisseur de gaz 11ou entre la sortie du refroidisseur de gaz et l'entrée de l'organe de détente 12. A partir de la valeur de la puissance Pk(x) obtenue à l'étape 208 et de la valeur de la différence d'enthalpie Ahk(x) obtenue à l'étape 210, la fonction de calcul de la carte électronique 401 fournit la grandeur de sortie, MC02-..t, selon l'équation Al de l'annexe A. A l'étape 214, la carte électronique 401détermine si le critère de fin d'itération est satisfait. Le critère de fin d'itération comprend une condition de seuil sur l'écart entre la grandeur de sortie MC02-,Ut et la grandeur d'entrée MC02-,. de la fonction de calcul. Plus précisément, la carte électronique détermine si la différence entre l'écart entre la grandeur de sortie mC02~out et la grandeur d'entrée mC02~1n est inférieure à une valeur de seuil prédéfinie e. Cette valeur de seuil est très petite, par exemple, sensiblement égale à 5%. Si c'est le cas, le critère de fin d'itération est satisfait et la grandeur de sortie mCO2~out de la fonction de calcul comprend la valeur du débit de fluide mCO2 (étape 218). Sinon, la carte électronique appelle itérativement la fonction de calcul (étapes 204 à 214) jusqu'à ce que le critère de fin d'itération soit satisfait. A chaque nouvelle itération de la fonction de calcul, la carte électronique affecte, à l'étape 216, la valeur de la grandeur de sortie de l'itération précédente à la grandeur d'entrée mCO2~m de la nouvelle itération. La valeur du débit mCO2 obtenue à l'étape 218 peut alors être utilisée par une fonction d'estimation de puissance de la carte électronique pour estimer la puissance réelle consommée par le compresseur, comme on le verra plus loin. Dans une deuxième forme de réalisation de l'invention, on estime le débit massique de fluide frigorigène mCO2, à partir d'un calcul préalable de l'efficacité du refroidisseur de gaz entre l'entrée x, et le point x. L'efficacité # du refroidisseur de gaz sur la portion comprise entre x, et x est exprimée en fonction des deux écarts de températures suivants: - l'écart entre la température Tgcl à l'entrée du refroidisseur de gaz et la température Tgcx du fluide C02 au point x, et - l'écart entre la température Tg, à l'entrée du refroidisseur de gaz et la température de l'air reçu par le refroidisseur de gaz Tal. L'équation B2. 1 qui permet d'obtenir # en fonction de ces paramètres est donnée dans l'annexe B. On définit l'efficacité normale #N par l'équation B2. 2 de l'annexe B, dans laquelle F est un facteur de correction qui prend en compte la haute pression P1. L'efficacité # est une fonction décroissante du débit de fluide qui traverse le refroidisseur de gaz, selon l'équation de l'annexe B 1. l, où a, et s sont des constantes qui dépendent des caractéristiques du refroidisseur de gaz et du débit d'air reçu par ce dernier. Le débit massique de C02 est donc lié à l'efficacité normale selon la relation B 1. 2 de l'annexe B. Les figures 6 à 10 représentent des exemples d'évolution de l'efficacité normale % en fonction du débit de C02, dans lesquels le point x est choisi à une distance D= 15% de L de l'entrée x,. La figure 11met en évidence le manque de précision de l'estimation du débit de fluide obtenu avec un positionnement x=xo. Sur chacune de ces figures, la courbe représentée en traits pleins représente l'évolution de exp( 11N) en fonction de mC02 et la courbe constituée de petits losanges représente l'évolution de % en fonction de mCO2. Selon la deuxième forme de réalisation de l'invention, la fonction de calcul de la carte électronique effectue dans une première étape une estimation/mesure de la température Tgcl d'entrée du refroidisseur de gaz, de la température Tgcx du fluide C02 au point x choisi et de la température Tal de l'air reçu par le refroidisseur de gaz, comme dans l'étape 200 de la première forme de réalisation. Dans une deuxième étape, la fonction de calcul détermine l'efficacité normale selon les équations Blet B3 de l'annexe B. Dans une dernière étape, la fonction de calcul résout l'équation de l'annexe B1pour déterminer le débit massique de C02, mCO2, à partir de la valeur obtenue de l'efficacité normale. L'estimation du débit massique de fluide frigorigène, obtenue selon les deux formes de réalisation de l'invention, peut être utilisée pour calculer, selon les relations de l'annexe C, la puissance mécanique absorbée. Pour cela, la carte électronique comprend une fonction d'estimation de puissance capable d'estimer la puissance absorbée par le compresseur à partir du travail de compression isentropique Wis, et de la vitesse de rotation du compres- seur N. Les constantes k4, k5 et k6 sont liées à des paramètres de fonctionnement du circuit de climatisation. La fonction d'estimation de puissance met en oeuvre l'équation CIde l'annexe C pour calculer le travail de compression isentropique Wis à partir : - du taux de compression Pr représentant le rapport entre la haute pression P 1 et la basse pression PO, et - de la température d'aspiration Tsu du compresseur. La pression P1 peut être mesurée comme décritprécédemment pour l'étape 210 de la figure 7. La température d'aspiration Tsu peut être directement mesurée par la sonde 180 de la figure lb, placée à l'entrée du compresseur. En variante, la température d'aspiration Tsu peut être estimée à partir de l'efficacité #Hxl de l'échangeur interne 23, de la température d'évaporation Tev et de la température en sortie du refroidisseur de gaz Tgco, selon l'équation de l'annexe C3. Selon cette variante, pour estimer la température Tsu, l'installation peut comporter une sonde placée à l'entrée de l'évaporateur 13pour mesurer la température d'évaporation Tev, désignée par la référence 170 de la figure 1b. Cette sonde peut aussi être une sonde à thermistance, de structure classique, placée dans les ailettes de l'évaporateur. En variante, la température Tev peut être estimée à partir de la température d'air soufflé en sortie de l'évaporateur fournie par le régulateur d'habitacle 41. L'installation peut en outre comporter un capteur de température 150 pour mesurer Tgco. Ce capteur est placé en sortie du refroidisseur externe 11, comme indiqué sur la figure 1b. Dans une autre variante encore, la température d'aspiration Tsu peut être estimée en fonction de la température de refoulement Td, du taux de compression Pr et d'une constante a selon la relation C4 de l'annexe A4. La valeur de la basse pression PO, utilisée pour le calcul du travail isentropique Wis, peut être mesurée directement par un capteur 140 placé entre l'évaporateur et le compresseur. Cette mesure est transmise à la carte électronique 401, par la liaison 30. La valeur de la basse pression PO peut également être estimée à partir de la valeur de la température d'évaporation Tev, selon la loi de saturation des fluides. A partir de la valeur du travail isentropique Wis, de la valeur du débit de fluide supercritique mCO2 estimée et de la vitesse de rotation N, la fonction d'estimation de puissance détermine une estimation de la puissance absorbée par le compresseur Pa, selon la relation C2 de l'annexe C. La vitesse de rotation du compresseur N est fournie à la carte électronique par le calculateur d'injection du moteur 42, par la liaison 33 en référence à la figure lb. Le calculateur utilise alors la valeur estimée de la puissance réelle consommée par le compresseur pour ajuster les paramètres d'injection. Par suite la consommation de carburant est réduite. La présente invention vise également le code logiciel qu'elle fait intervenir, tout particulièrement lorsqu'il est mis à disposition sur tout support lisible sur un ordinateur. L'expression "support lisible par ordinateur" couvre un support de stockage, par exemple magnétique ou optique, aussi bien qu'un moyen de transmission, tel qu'un signal numérique ou analogique. ANNEXE A Al. Débit massique du fluide frigorigène supercritique moco2= Pk(x)/hk(x) A2. Puissance thermique échangée entre x, et x Pk(x) = kS.OTLM .1/(1 + ksl(2ma,.Cpa) ) A3. Coefficient d'échange thermique global entre l'air et le C02 A3.1- Coefficient d'échange thermique global ksg = 1 (1/kas, + A + 1 /kso2) A3.2- Coefficient d'échange côté air A3. 21-Forme générale kSa = fi (U gmv' Va) A3.22-Exemple ksa = (a + b(UgmvY + d.Ve A3. 3- Coefficient d'échange côté fluide A3. 31-Forme générale kSC02 = f2(mCO2) A3. 32-Exemple k SC02 - -, a. ( IDC02 )b' A4.différence d'enthalpie du fluide frigorigène entre x, et x Ah,(x)=h(Pl,T)-h(Pl,T,J AS.moyenne logarithmique de la différence de température globale entre x, et x ATLM = [(TgCl - Ta) - (Tgcx - TJ] / log[(Tg, - TJ /(tacx - TJ] ANNEXE B Bl.Relation entre l'efficacité normale #N et le débit massique mCO2- Bl.1- efficacité 11 en fonction du débit mCO2 11 = a1.mc02b B1.2-efficacité 11N en fonction du débit MC02 exp('IN) = kl-exp(-k2-MC02) B2.Efficacité du fluide frigorigène supercritique # B2.1-Efficacité du fluide frigorigène # = (Tgci - Tgcx)/ (Tgc. 'Tal) B2. 2-Efficacité normale du fluide frigorigène #N = F.# B3.Exemples de relations entre l'efficacité normale et le débit massique mCO2 avec y= exp (#N) B3. 1-débit d'air= 2200 kg/h et x=15%.L y = 0,9296.exp(-0,003-MC02) B3.2-débit d'air= 1500 kg/h et x=15%.L y = 0,9119.exp(-0,0036.mco2) B3. 3-débit d'air= 1200 kg/h et x=15%.L y = 0,898.exp(-0,004.mco2) B3. 4-débit d'air= 760 kg/h et x=15%.L y = 0,9296.exp(-0,003'MC02) B3.5-débit d'air= 760 kg/h et x=100%.L y = 1,0247.exp(-0,0013.mcoz) ANNEXE C Cl. Estimation du travail de compression fourni par le compresseur Ws = k6.Tsu.(Pr-l) C2. Puissance mécanique absorbée par la climatisation Pa = k4.mC02'Ws + ks.N C3. Température d'aspiration Tu dans un circuit à fluide supercritique Tsu = 1HX (Tgco - Tev) + Tel C4. Température d'aspiration Tsu dans un circuit à fluide supercritique Td = T su* Pr", où a est une constante Revendications 1. Installation de climatisation pour véhicule à moteur, munie d'un circuit (10) de fluide frigorigène supercritique comprenant un compresseur ( 14), un refroidisseur de gaz ( 11 ) de longueur totale donnée (L), un organe de détente (12) et un évaporateur (13), le refroidisseur de gaz recevant un flux d'air (16) passant par un groupe moto-ventilateur, l'installation comprenant en outre un dispositif électronique de contrôle (401) destiné à interagir avec le circuit de fluide frigorigène (10) pour estimer une grandeur relative au fluide frigorigène (mCO2) à partir de températures relatives au refroidisseur de gaz, caractérisée en ce que ledit dispositif électronique de contrôle (401) comprend une fonction de calcul utilisant deux écarts de températures relatives au refroidisseur de gaz, dont l'un au moins est fondé sur la température du fluide frigorigène (Tgcx) en un point intermédiaire choisi du refroidisseur de gaz (x). 2. Installation de climatisation selon la revendication 1, caractérisée en ce que le point intermédiaire choisi est situé à une distance de l'entrée (x,) du refroidisseur de gaz, comprise entre 5 % et 35% de la longueur totale du refroidisseur de gaz (L). 3. Installation de climatisation selon l'une des revendications 1 et 2, caractérisée en ce que les deux écarts de température sont en outre fondés sur : - la température du fluide frigorigène à l'entrée du refroidisseur de gaz (Tgcl), - la température du flux d'air reçu par le refroidisseur de gaz (Tal). 4. Installation de climatisation selon l'une des revendications 1 à 3, caractérisée en ce que la grandeur relative au fluide frigorigène est le débit massique du fluide frigorigène (mCO2). 5. Installation de climatisation selon la revendication 4, caractérisée en ce que les deux écarts de températures utilisés par la fonction de calcul comprennent : - l'écart entre la température du fluide frigorigène à l'entrée du refroidisseur de gaz (Tgcl) et la température du flux d'air reçu par le refroidisseur de gaz (Ta,); -l'écart entre la température du fluide frigorigène (Tgcx) au point intermédiaire choisi du refroidisseur de gaz (x) et la température du flux d'air reçu par le refroidisseur de gaz (Tal). 6. Installation de climatisation selon la revendication 5, caractérisée en ce que le dispositif électronique de contrôle (401) est agencé pour appeler itérativement la fonction de calcul, à partir d'une valeur initiale donnée du débit massique du fluide frigorigène (MC02 init), jusqu'à satisfaire un critère de fin d'itération, comprenant une condition de seuil sur l'écart entre la grandeur d'entrée (mCO2~m) et la grandeur de sortie de la fonction de calcul (MC02-out), 7. Installation de climatisation selon la revendication 6, caractérisée en ce que la grandeur d'entrée de la première itération de la fonction de calcul comprend ladite valeur initiale (MC02-init)* 8. Installation de climatisation selon l'une des revendications 6 et 7, caractérisée en ce que la grandeur d'entrée (mCO2 in) d'une itération donnée de la fonction de calcul comprend la grandeur de sortie (mC02 out) de l'itération précédente de la fonction de calcul. 9.Installation de climatisation selon l'une des revendications 6 à 8, caractérisée en ce que le critère de fin d'itération est satisfait et en ce que la grandeur de sortie comprend la valeur du débit massique du fluide frigorigène (mCO2). 10. Installation de climatisation selon l'une des revendications 6 à 9, caractérisée en ce que la fonction de calcul (401) est capable de calculer la grandeur de sortie d'une itération donnée, à partir: - de la puissance (Pk(x)) échangée entre le fluide frigorigène et le flux d'air reçu, sur une portion du refroidisseur de gaz comprise entre l'entrée du refroidisseur de gaz (xi) et le point intermédiaire choisi (x), et - de la différence d'enthalpie du fluide frigorigène entre l'entrée du refroidisseur de gaz et ledit point intermédiaire choisi (#hk(x)). 11 .Installation de climatisation selon l'une des revendications 6 à 10, caractérisée en ce que la condition de seuil comprend le fait que l'écart entre la grandeur d'entrée (mCO2~m) et la grandeur de sortie (mC02~out) de la fonction de calcul est inférieur à une valeur de seuil (e). 12.Installation de climatisation selon la revendication 11, caractérisée en ce que la valeur de seuil (e) est sensiblement égale à 5%. 13.Installation de climatisation selon l'une des revendications 10 à 12, caractérisée en ce que la fonction de calcul est apte à calculer la puissance échangée entre le fluide frigorigène et le flux d'air reçu sur ladite portion, à partir de la moyenne logarithmique de la différence de température globale (ATLM) sur ladite portion et du coefficient d'échange thermique global (ksg) entre l'air et le fluide frigorigène. 14. Installation de climatisation selon la revendication 13, caractérisée en ce que la fonction de calcul est apte à calculer le coefficient d'échange thermique (ksg) à partir de la grandeur d'entrée (Il1c02~m), de la vitesse d'avancement (Va) du véhicule et de la tension du groupe moto-ventilateur (UGMV)' 15. Installation de climatisation selon l'une des revendications 12 à 14, caractérisée en ce que le dispositif électronique de contrôle (401) est apte à calculer la moyenne logarithmique de la différence de température globale (ATLM) sur ladite portion à partir des deux écarts de température. 16. Installation de climatisation selon l'une des revendications 10 à 15, caractérisée en ce que la fonction de calcul est apte à calculer la différence d'enthalpie du fluide frigorigène, entre l'entrée du refroidisseur de gaz et le point intermédiaire choisi (x), à partir : - de la pression du fluide en sortie du compresseur (P1), - de la température du fluide frigorigène à l'entrée du refroidisseur de gaz (Tgcl), et - de la température du fluide frigorigène au point intermédiaire choisi du refroidisseur de gaz (Tgcx). 17. Installation de climatisation selon la revendication 4, caractérisée en ce que les deux écarts de températures utilisés par la fonction de calcul comprennent : - l'écart entre la température du fluide frigorigène à l'entrée du refroidisseur de gaz (Tgcl) et la température du fluide frigorigène (Tgcx) au point intermédiaire choisi du refroidisseur de gaz (x), - l'écart entre la température du fluide frigorigène à l'entrée du refroidisseur de gaz (Tgcl) et la température du flux d'air reçu par le refroidisseur de gaz (Tal). 18. Installation de climatisation selon l'une des revendications précédentes, caractérisée en ce qu'elle comporte une sonde de température placée au point intermédiaire choisi du refroidisseur de gaz (x) et propre à fournir une mesure de la température du fluide frigorigène (Tgcx) au point intermédiaire choisi. 19. Installation de climatisation selon l'une des revendications 3 à 18, caractérisée en ce qu'elle comporte une sonde de température (130), placée à l'entrée du refroidisseur de gaz et propre à fournir une mesure de la température du fluide frigorigène à l'entrée du refroidisseur de gaz (Tgcl). 20. Installation de climatisation selon l'une des revendications 3 à 18, caractérisée en ce que le dispositif de contrôle comprend une fonction d'estimation de température de fluide capable de calculer une estimation de la température du fluide frigorigène à l'entrée du refroidisseur de gaz (Tg,) à partir : - de la pression du fluide en entrée du compresseur (PO), - de la température du fluide à l'entrée du compresseur (Tsu)' et - de la pression du fluide en sortie du compresseur (P1). 21.Installation de climatisation selon l'une des revendications 3 à 20, caractérisée en qu'elle comporte une sonde (120), placée entre le groupe moto-ventilateur (15) et le refroidisseur de gaz (11), propre à fournir une mesure de la température du flux d'air reçu par le refroidisseur de gaz (Tal). 22. Installation de climatisation selon l'une des revendications 3 à 20, caractérisée en ce que le dispositif électronique de contrôle (401) comprend une fonction d'estimation de température d'air capable de calculer la température du flux d'air reçu par le refroidisseur de gaz (Ta,) à partir : - de la température de l'air extérieur (Text) et -de la vitesse d'avancement (Va) du véhicule. 23. Installation de climatisation selon la revendication 22, caractérisée en ce en qu'elle comporte une sonde propre à fournir une mesure de la température du flux d'air extérieur (Text). 24. Installation de climatisation selon l'une des revendications 4 à 23, caractérisée en ce que le dispositif électronique de contrôle comporte une fonction d'estimation de puissance capable d'estimer la puissance absorbée par le compresseur à partir: - du taux de compression (Pr), et - de la température du fluide frigorigène en 'entrée du compresseur (Tsu). 25. Procédé d'estimation du débit massique d'un fluide frigorigène supercritique circulant dans un circuit de climatisation (10), ledit circuit comprenant un compresseur (14), un refroidisseur de gaz (11) de longueur totale donnée (L), un organe de détente (12) et un évaporateur (13), le refroidisseur de gaz recevant un flux d'air (16) passant par un groupe moto-ventilateur, caractérisé en ce qu'il comprend les étapes suivantes, pour une valeur initiale donnée du débit massique du fluide frigorigène (mC02~lmt) : a) calculer - une estimation de l'écart entre la température du fluide frigorigène à l'entrée du refroidisseur de gaz (Tgct) et la température du flux d'air reçu par le refroidisseur de gaz (Ta,); - une estimation de l'écart entre la température du fluide frigorigène (Tgcx) au point intermédiaire choisi du refroidisseur de gaz (x) et la température du flux d'air reçu par le refroidisseur de gaz (Ta,) b) calculer une grandeur de sortie(MC02-.,,,) , à partir d'une grandeur d'entrée (mC02m) et des valeurs fournies par l'étape a); c) déterminer si l'écart entre la grandeur de sortie calculée à l'étape b) et la grandeur d'entrée est inférieur à une valeur de seuil prédéfinie (#). 26. Procédé selon la revendication 25, caractérisé en ce que la grandeur d'entrée (mCO2 m) comprend la valeur initiale du débit de fluide frigorigène (mC02~'mt)' à la première itération de l'étape a). 27. Procédé selon l'une des revendications 25 et 26, caractérisé en ce que les étapes a) à c) sont réitérées avec une grandeur d'entrée (mC02~ln) égale à la grandeur de sortie (mC02~out) calculée à l'étape c) précédente, si l'étape c) n'est pas vérifiée. 28. Procédé selon l'une des revendications 26 à 27, caractérisé en ce que la grandeur de sortie comprend la valeur du débit massique de fluide frigorigène (mCO2), si l'étape c) est vérifiée.