FR 2941267 A1 20100723 FR 0950396 A 20090122 La présente invention concerne un procédé de détermination de la pression en amont d'une turbine d'un turbocompresseur de suralimentation d'un moteur thermique. Dans le domaine de la mesure de pression il est généralement connu d'utiliser un capteur, par exemple du type piézoélectrique mesurant une variation de pression. Cependant l'implantation d'un tel capteur est coûteuse. La présente invention se propose de remplacer un capteur de pression par un estimateur. L'invention a pour objet un procédé de détermination, pour un turbocompresseur de suralimentation d'un moteur thermique comprenant une turbine entraînée par les gaz d'échappement issus dudit moteur thermique et mécaniquement solidaire en rotation d'un compresseur afin de comprimer l'air d'admission injecté dans le moteur thermique, de la pression en amont de la turbine en fonction du débit d'air d'admission au travers du compresseur, de la pression en amont du compresseur, de la température en amont du compresseur, de la pression en aval du compresseur, de la température en amont de la turbine et de la pression en aval de la turbine. D'autres caractéristiques, détails et avantages de l'invention ressortiront plus clairement de la description détaillée donnée ci-après à titre indicatif en relation avec des dessins sur lesquels : - la figure 1 illustre un moteur thermique avec un turbocompresseur de suralimentation, - la figure 2 illustre un moteur thermique équipé d'un dispositif de suralimentation comprenant deux 30 turbocompresseurs, - la figure 3 présente un schéma figurant les variables d'entrée et de sortie du procédé, - la figure 4 présente un bloc-diagramme d'un premier mode de réalisation du procédé selon l'invention, 35 - la figure 5 présente un bloc-diagramme d'un second mode de réalisation du procédé selon l'invention, - les figures 6-10 présentent les cartographies respectives des fonctions f1, f2, f3, f4 et f5, - les figures 11-14 présentent les définitions numériques respectives des fonctions f1-f4, - la figure 15 illustre la qualité du résultat produit par le procédé. Afin de faciliter la lecture de la description, des blocs-diagrammes et particulièrement des formules, il est utilisé les notations suivantes : variables: N : régime ou vitesse de rotation (du turbocompresseur), R rapport de pression (rapport de compression du compresseur, rapport de détente de la turbine), Q : débit, P : pression H : puissance, T : température, . rendement, Cp constante thermodynamique, capacité thermique à pression constante, Cv : constante thermodynamique capacité thermique à volume 20 constant, y : constante thermodynamique, coefficient égal à Cp/Cv, J : moment d'inertie (du turbocompresseur). indices . c : compresseur, 25 t : turbine, cor : grandeur corrigée, ref : grandeur de référence, u : amont, d : aval, 30 n : indice du temps, pas de calcul courant, n-1 : pas de calcul courant précédent. La figure 1 illustre le contexte de l'invention. Un moteur thermique 4 reçoit classiquement par des tubulures 35 d'admission 6 de l'air 5. Le moteur 4 produit des gaz d'échappement 7 qui s'échappent par des tubulures d'échappement 8. Un turbocompresseur 1 de suralimentation permet d'augmenter la quantité d'air 5 admise par le moteur thermique 4. Pour cela le turbocompresseur 1 comprend une turbine 2 et un compresseur 3. La turbine 2 est fluidiquement connectée sur les tubulures d'échappement 8 afin d'être entraînée par les gaz d'échappement 7 issus du moteur thermique 4. La turbine 2 est mécaniquement solidaire du compresseur 3 qu'elle entraîne en rotation. Le compresseur 3 est fluidiquement connecté sur les tubulures d'admission 6, afin que le compresseur 3 comprime l'air d'admission 5 avant son entrée dans le moteur thermique 4. Il est possible d'isoler la turbine 2 au moyen d'une vanne de by-pass 11. Il est possible d'isoler le compresseur au moyen d'une vanne de by-pass 10. Le repère 9 figure un capteur de débit d'air d'admission 5. Le schéma de la figure 3 illustre ce même environnement et présente les variables du système. Le turbocompresseur 1 est relié au moteur 4. La turbine 2 est disposée sur l'échappement 8. Le compresseur 3 est disposé sur l'admission 6. Le problème tel que posé suppose que l'on souhaite estimer la pression Put en amont de la turbine 2, figurée encadrée sur la figure 3. On suppose connues les grandeurs suivantes : débit Q, (non représenté) d'air d'admission au travers du compresseur 3, la pression Pu, en amont du compresseur 3, la température Tu, en amont du compresseur 3, la pression Pdc en aval du compresseur 3, la température Tut en amont de la turbine 2 et la pression Pdt en aval de la turbine 2. La connaissance de cette pression Put en amont de la turbine 2 est primordiale afin de piloter finement ledit turbocompresseur 1 afin d'éviter sa détérioration et d'améliorer le brio du véhicule lors des transitoires. Cependant, il n'est pas souhaité de recourir à un capteur de pression. L'objet de l'invention est ainsi une méthode d'estimation de cette pression en fonction des six autres grandeurs connues par ailleurs. La figure 2 illustre un cas particulier d'utilisation. Ici un second turbocompresseur 15 est ajouté en série. La suralimentation est alors réalisée par un double turbocompresseur étagé. Le second turbocompresseur 15 réalise une première compression de l'air d'admission 5. Il est encore dénommé turbocompresseur basse pression. Le premier turbocompresseur 1 réalise ensuite une seconde compression de l'air d'admission 5 issu du compresseur du turbocompresseur 15 basse pression. Le premier turbocompresseur 1 est encore dénommé turbocompresseur 1 haute pression. Une vanne de by-pass 12 permet d'isoler la turbine basse pression. L'invention s'applique particulièrement au cas du turbocompresseur 1 haute pression. La méthode est particulièrement adaptée à un turbocompresseur à géométrie fixe. Dans cette configuration particulière, les six grandeurs d'entrée du procédé selon l'invention sont avantageusement déterminées au moyen de capteurs pour le débit d'air Q, d'admission au travers du compresseur 3, la pression Pdc en aval du compresseur 3 et la température Tut en amont de la turbine 2, tandis que la pression Puc en amont du compresseur 3, la température Tuc en amont du compresseur 3 et la pression Pdt en aval de la turbine 2 sont déterminées par un estimateur déterminant les grandeurs du turbocompresseur 15 basse pression. Comme visible sur la figure 2, la pression Pdt en aval de la turbine 2 haute pression est égale à la pression en amont 25 de la turbine basse pression. Il peut être nécessaire de refroidir l'air d'admission 5. Le choix a été fait de n'utiliser qu'un unique échangeur thermique 13, le cas échéant, disposé en aval du compresseur 3. Ainsi, l'absence d'échangeur dans la tubulure d'admission 30 6 entre le compresseur basse pression et le compresseur 3 haute pression, permet de connaître la température Tuc en amont du compresseur 3 haute pression, puisqu'elle est égale à la température à l'aval du compresseur basse pression. Le principe du procédé selon l'invention est illustré 35 selon deux modes de réalisation par les blocs diagrammes des figures 4 et 5. Le procédé de détermination de la pression Put en amont de la turbine 2 peut arbitrairement être découpé en six étapes suivantes . 1) calcul du régime corrigé Ncor du turbocompresseur 1, en fonction du rapport de compression Rc du compresseur 3 et du débit corrigé Qccor d'air d'admission au travers du compresseur 3, 2) calcul du régime N du turbocompresseur 1 en fonction du régime corrigé Ncor du turbocompresseur 1 et de la température Tuc en amont du compresseur 3, 3) calcul de la puissance Hc du compresseur 3 en fonction du débit Qc d'air d'admission au travers du compresseur 3, du rendement rlc du compresseur 3, de la température Tuc en amont du compresseur 3 et du rapport de compression Rc du compresseur 3, 4) calcul de la puissance Ht de la turbine 2 en fonction 15 du régime N du turbocompresseur 1 et de la puissance Hc du compresseur 3, 5) calcul du rapport de détente Rt de la turbine 2, 6) calcul de la pression Put en amont de la turbine 2 en fonction de la pression Pdt en aval de la turbine 2 et du 20 rapport de détente Rt de la turbine 2. A noter que les étapes 1-4 et 6 sont identiques dans les deux modes de réalisations. Seule l'étape 5 les différencie. A l'étape 1) est calculé le régime corrigé Ncor du turbocompresseur 1, en fonction du rapport de compression Rc 25 du compresseur 3 et du débit corrigé Qccor d'air d'admission au travers du compresseur 3 au moyen d'une fonction f1. Cette fonction f1 du rapport de compression Rc du compresseur 3 et du débit corrigé Qccor d'air d'admission au travers du compresseur 3 est calculée au bloc f1. Cette fonction f1 est 30 définie par une cartographie à deux dimensions. Une cartographie est un moyen connu de définir une fonction f. Ladite fonction f est définie de manière graphique par une courbe (cartographie à une dimension) ou une surface (cartographie à deux dimensions). De manière 35 connue et classique le résultat z de la fonction f(x) = z (une dimension) ou f(x,y) = z (deux dimensions) est déterminé graphiquement à partir de la donnée de la courbe ou de la surface. Cette même fonction f peut encore de manière équivalente être définie par un tableau (mono ou bidimensionnel) de nombres. Ainsi la fonction f1 est, par exemple, définie par la surface de la figure 6 ou de manière équivalente par un tableau bidimensionnel de nombres. Ainsi la fonction f1 est parfaitement définie par le tableau de la figure 11 où x se lit sur la première colonne, y sur la première ligne et le résultat z au croisement de la ligne x et de la colonne y. De manière connue le résultat est déterminé par interpolation lorsque les valeurs x ou y ne sont pas directement présentes dans le tableau. Les différentes cartographies de fonction f1-f5, sont ainsi déterminées pour un compresseur 3 et une turbine 2 donnés à titre illustratif et représentés respectivement aux figures 6-10. En cas d'application à une turbine 2 ou à un compresseur 3 différents de ceux considérés ici, l'homme du métier sait déterminer les cartographies des fonctions f1-f5, directement ou en adaptant (mise à l'échelle, changement d'unité...) des cartographies de fonctionnement fournies par les constructeurs de ces machines tournantes 2, 3. Le rapport de compression Re du compresseur 3 est par définition égal au rapport de la pression amont Pu, du compresseur 3 à la pression aval Pde du compresseur 3 et est calculé au bloc 20. Le débit corrigé Qccor d'air d'admission du compresseur 3 est calculé selon la formule : Q e_cor P Zlc où Pc ref or est le débit corrigé d'air d'admission 5 au travers du compresseur 3, Tu, est la température en amont du compresseur 3, 30 Pu, est la pression en amont du compresseur 3, Tcref est une température de référence du compresseur 3, Pcref est une pression de référence du compresseur 3. Cette formule est mise en oeuvre au bloc 21. Les température Te ref et pression Pcref de référence sont 35 définies afin de permettre un calcul simplifié des différentes fonctions cartographiées fl-f5 en se ramenant toujours à des conditions de références permettant l'utilisation d'une unique cartographie pour chaque fonction fl-f5. Les températures et pressions de références sont dans les exemples illustratifs fournis, égales à : Tc ref = 298°K, Tt ref = 873°K, Pc ref = Pt ref = 1 atm. A l'étape 2) est calculé le régime N du turbocompresseur 1 selon la formule N = NeOr Tue où c ref N est le régime du turbocompresseur 1, Ncor est le régime corrigé du turbocompresseur 1, Tuc est la température en amont du compresseur 3, Tcref est la température de référence du compresseur 3, décrite précédemment. Cette formule est mise en oeuvre au bloc 22. A l'étape 3) est calculé la puissance Hc du compresseur 3 ( y,-1 \ , où i est la puissance du compresseur 3, est le débit d'air d'admission au travers du compresseur 3, rlc est le rendement du compresseur 3, Tuc est la température en amont du compresseur 3, Rc est le rapport de compression du compresseur 3, Cpc est une première constante thermodynamique de l'air d'admission, yc est une seconde constante thermodynamique de l'air d'admission. Cette formule est mise en oeuvre au bloc 23. Le rendement rlc du compresseur 3, qui est une entrée de ladite étape 3), est calculé en fonction du régime corrigé Ncor du turbocompresseur 1 et du débit corrigé Qccor d'air d'admission au travers du compresseur 3, au moyen d'une fonction f2 du régime corrigé Ncor du turbocompresseur 1 et du débit corrigé Qccor d'air d'admission au travers du compresseur 3, cette fonction est réalisée au bloc f2. Ladite fonction f2 est définie par une cartographie à deux selon la formule : He =QeCpe - Tue ReY` -1 lie dimensions. La figure 7 illustre la cartographie de la fonction f2. La fonction f2 est encore définie par le tableau de la figure 12. Dans la formule précédente, la première constante thermodynamique Cp, de l'air d'admission 5 est la capacité thermique de l'air d'admission 5 à pression constante et est égale à 1005 J/kg/K, et la seconde constante thermodynamique Yc de l'air d'admission 5 est le coefficient rapport Cp,/Cv, des capacités thermiques de l'air d'admission 5 respectivement à pression constante et à volume constant et est égale à 1,4. A l'étape 4) est ensuite calculée la puissance Ht de la turbine 2 selon la formule : Hz =JNdNùHe, où dt Ht est la puissance de la turbine 2, H, est la puissance du compresseur 3, N est le régime du turbocompresseur 1, d est l'opérateur de dérivation par rapport à la L'étape 5) a pour but de calculer le rapport de détente 25 Rt de la turbine 2. Il est ici proposé deux modes de réalisation de cette étape 5) conduisant respectivement au schéma bloc de la figure 4 et de la figure 5. Selon un premier mode de réalisation illustré au schéma 30 bloc de la figure 4 le rapport de détente Rt de la turbine 2 est calculé en fonction du débit corrige Qtcor de gaz d'échappement 7 au travers de la turbine 2 au moyen d'une fonction f4 du débit corrigé Qt cor de gaz d'échappement 7 au travers de la turbine 2, réalisée au bloc f4. Cette fonction 35 f4 est définie par une cartographie à une dimension. La figure 9 illustre la cartographie de la fonction f4. La dt variable temps et, J est une constante égale au moment d'inertie du 20 turbocompresseur 1. Cette formule, issue de la relation fondamentale de la dynamique est mise en oeuvre au bloc 24. fonction f4 est encore définie par le tableau de la figure 14. Ce débit corrigé Qt cor de gaz d'échappement 7 au travers de la turbine 2 est calculé selon la formule Lt_COY =Q t T ~Ut , où rut(n 1) Qtcor est le débit corrigé de gaz d'échappement 7 au travers de la turbine 2, Qt est le débit de gaz d'échappement 7 au travers de la turbine 2, Tut est la température en amont de la turbine 2, Put est la pression en amont de la turbine 2, l'indice n-1 indiquant ici qu'elle est déterminée à l'intervalle de temps n-1 précédent l'intervalle de temps courant n. Cette formule est mise en oeuvre au bloc 26. Le débit Qt de gaz d'échappement 7 au travers de la turbine 2 est calculé selon la formule 't' t' I t l ut Rt(nû1) i Qt est le débit de gaz d'échappement 7 au travers de la turbine 2, Ht est la puissance de la turbine 2, rit est le rendement de la turbine 2, Tut est la température en amont de la turbine 2, Rt est le rapport de détente de la turbine 2, l'indice n-1 indiquant ici qu'il est déterminé à l'intervalle de 25 temps n-1 précédent, Cpt est une première constante thermodynamique du gaz d'échappement 7, yt est une seconde constante thermodynamique du gaz d'échappement 7. 30 Le bloc 28 est un bloc délai 1/z permettant de mémoriser la valeur Put(n-1) de la grandeur Put de l'intervalle de temps précédent n-1. Le bloc 29 est un bloc multiplicatif permettant de calculer Rt (n-1) en multipliant Put (n-1) par Pdt. Ht a= 1 1 Selon un second mode de réalisation illustré au schéma bloc de la figure 5 le rapport de détente Rt de la turbine 2 est calculé en fonction de la puissance Ht de la turbine 2, du débit Qt de gaz d'échappement 7 au travers de la turbine 2, du rendement nt de la turbine 2, de la température Tut en amont de la turbine 2, selon la formule H Rt= 1û , où Qt (n -1)Cp,fl,T t Rt est le rapport de détente de la turbine 2, Ht est la puissance de la turbine 2, Qt est le débit de gaz d'échappement 7 au travers de la turbine 2, l'indice n-1 indiquant ici qu'il est déterminé à l'intervalle de temps n-1 précédent, nt est le rendement de la turbine 2, Tut est la température en amont de la turbine 2, Cpt est une première constante thermodynamique du gaz d'échappement 7, Yt est une seconde constante thermodynamique du gaz d'échappement 7. Cette formule est mise en oeuvre au bloc 30. Le débit Qt de gaz d'échappement 7 au travers de la turbine 2 est calculé en fonction du débit corrigé Qt cor de gaz d'échappement 7 au travers de la turbine 2 selon la Pt(n -1) formule : Q,(11ù1)=Q', cor , où ut Qt est le débit de gaz d'échappement 7 au travers de la turbine 2, l'indice n-1 indiquant ici qu'il est déterminé à l'intervalle de temps n-1 précédent, Qt~or est le débit corrigé de gaz d'échappement 7 au travers de la turbine 2, Put est la pression en amont de la turbine 2, l'indice n-1 indiquant ici qu'elle est déterminée à l'intervalle de temps n-1 précédent, et Tut est la température en amont de la turbine 2. Cette formule est mise en oeuvre au bloc 31. Le débit corrigé Qt~or de gaz d'échappement 7 au travers de la turbine 2 est calculé en fonction du rapport de détente Rt de la turbine 2 au moyen d'une fonction f5 du rapport de détente Rt de la turbine 2. Cette fonction est réalisée au bloc f5. Ladite fonction f5 est définie par une cartographie à une dimension. La figure 10 illustre la cartographie de la fonction f5. La fonction f5 est la fonction inverse de la fonction f4. La fonction f5 est encore définie par le tableau de la figure 14. Dans les formules précédentes des blocs 25 et 31, la première constante thermodynamique Cpt du gaz d'échappement 7 est la capacité thermique du gaz d'échappement 7 à pression constante et est égale à 1136 J/kg/K, et la seconde constante thermodynamique yt du gaz d'échappement 7 est le coefficient rapport Cpt/Cvt des capacités thermiques du gaz d'échappement 7 respectivement à pression constante et à volume constant et est égale à 1,34. Les deux variantes de l'étape 5) selon les deux modes de réalisation nécessitent une détermination du rendement nt de la turbine 2. Ce dernier est calculé en fonction du régime corrigé Ncor du turbocompresseur 1 et du rapport de détente Rt(n-1) de la turbine 2 déterminé à l'intervalle de temps n-1 précédent, au moyen d'une fonction f3 du régime corrigé Ncor du turbocompresseur 1 et du rapport de détente Rt de la turbine 2, réalisée au bloc f3. Ladite fonction f3 est définie par une cartographie à deux dimensions. La figure 8 illustre la cartographie de la fonction f3. La fonction f3 est encore définie par le tableau de la figure 13. L'étape finale 6) calcule le résultat, à savoir la 30 pression Put en amont de la turbine 2, selon la formule . Pt =PdtRt, issue de la définition de Rt, où Put est la pression en amont de la turbine 2, Pdt est la pression en aval de la turbine 2 et Rt est le rapport de détente de la turbine 2, 35 précédemment déterminé à l'étape 5). Cette formule est mise en oeuvre au bloc multiplicatif 27. L'invention concerne encore un estimateur réalisé au moyen d'un dispositif logique, mécanique, électronique, hydraulique ou encore d'un contrôleur et son programme logiciel, apte à mettre en œuvre le procédé selon l'un des modes de réalisation précédemment décrits. La figure 15 présente de manière comparative les résultats obtenus par le procédé ou l'estimateur selon l'invention. Pour un même événement (transitoire à 2000tr/min) sont représenté surs un même système d'axes, la pression Put en amont de la turbine 2 en fonction du temps. La courbe 16 figure le résultat obtenu avec le premier mode de réalisation. La courbe 17 figure le résultat obtenu avec le second mode de réalisation. Le résultat est très satisfaisant lorsqu'on le compare avec une courbe de référence 18. REVENDICATIONS1. Procédé de détermination, pour un turbocompresseur (1) de suralimentation d'un moteur thermique (4) comprenant une turbine (2) entraînée par les gaz d'échappement (7) issus dudit moteur thermique (4) et mécaniquement solidaire en rotation d'un compresseur (3) afin de comprimer l'air d'admission (5) injecté dans le moteur thermique (4), de la pression (Put) en amont de la turbine (2) en fonction du débit (Qc) d'air d'admission (5) au travers du compresseur (3), de la pression (Puc) en amont du compresseur (3), de la température (Tuc) en amont du compresseur (3), de la pression (Pdc) en aval du compresseur (3), de la température (Tut) en amont de la turbine (2) et de la pression (Pdt) en aval de la turbine (2), caractérisé en ce qu'il comprend les étapes suivantes - calcul du régime corrigé (Ncor) du turbocompresseur (1), en fonction du rapport de compression (Rc) du compresseur (3) et du débit corrigé (Qccor) d'air d'admission (5) au travers du compresseur (3), - calcul du régime (N) du turbocompresseur (1) en fonction du régime corrigé (Ncor) du turbocompresseur (1) et de la température (Tuc) en amont du compresseur (3), - calcul de la puissance (Hc) du compresseur (3) en fonction du débit (Qc) d'air d'admission (5) au travers du compresseur (3), du rendement (roc) du compresseur (3), de la température (Tuc) en amont du compresseur (3) et du rapport de compression (Rc) du compresseur (3), - calcul de la puissance (Ht) de la turbine (2) en fonction du régime (N) du turbocompresseur (1) et de la puissance (Hc) du compresseur (3), - calcul du rapport de détente (Rt) de la turbine (2), - calcul de la pression (Put) en amont de la turbine (2) en fonction de la pression (Pdt) en aval de la turbine (2) et du rapport de détente (Rt) de la turbine (2). 2. Procédé selon la revendication 1, où le débit corrigécor) d'air d'admission (5) du compresseur (3) est calculé selon la formule : f P où uc Q e_cor Pcref Qccor est le débit corrigé d'air d'admission (5) au travers du compresseur (3), Tuc est la température en amont du compresseur (3), Puc est la pression en amont du compresseur (3), Tcref est une température de référence du compresseur (3), Pcref est une pression de référence du compresseur (3). 3. Procédé selon la revendication 1 ou 2, où le régime corrigé (Ncor) du turbocompresseur (1) est calculé en fonction du rapport de compression (Rc) du compresseur (3) et du débit corrige ((Dccor) d'air d'admission (5) au travers du compresseur (3), au moyen d'une fonction (fi) du rapport de compression (PRc) du compresseur (3) et du débit corrigé (Qccor) d'air d'admission (5) au travers du compresseur (3), ladite fonction (fi) étant définie par une cartographie à deux dimensions. 4. Procédé selon l'une quelconque des revendications 1 à 3, où le régime (N) du turbocompresseur (1) est calculé selon la formule : N =N cor Tue , OÙ T ref N est le régime du turbocompresseur (1), Ncor est le régime corrigé du turbocompresseur (1), 25 Tuc est la température en amont du compresseur (3), Tcref est une température de référence du compresseur (3). 30 formule He =QeCpe 1 Tue Re'' -1 , où llc i He est la puissance du compresseur (3), Qc est le débit d'air d'admission (5) au travers du 5. Procédé selon l'une quelconque des revendications 1 à 4, où la puissance (Hc) du compresseur (3) est calculée selon la ( y~ù1 \compresseur (3), ric est le rendement du compresseur (3), Tu, est la température en amont du compresseur (3), R, est le rapport de compression du compresseur (3), Cp, est une première constante thermodynamique de l'air d'admission (5), Yc est une seconde constante thermodynamique de l'air d'admission (5). 6. Procédé selon la revendication 5, où le rendement ('k) du compresseur (3) est calculé en fonction du régime corrigé (Ncor) du turbocompresseur (1) et du débit corrigé (Qccor) d'air d'admission (5) au travers du compresseur (3), au moyen d'une fonction (f2) du régime corrigé (Ncor) du turbocompresseur (1) et du débit corrigé (Qccor) d'air d'admission (5) au travers du compresseur (3), ladite fonction (f2) étant définie par une cartographie à deux dimensions. 7. Procédé selon la revendication 5 ou 6, où la première constante thermodynamique (Cpc) de l'air d'admission (5) est égale à 1005 J/kg/K, et où la seconde constante thermodynamique (Yc) de l'air d'admission (5) est égale à 1,4. 8. Procédé selon l'une quelconque des revendications 1 à 7, où la puissance (Ht) de la turbine (2) est calculée selon la formule : Ht =JN dt ùHe, où Ht est la puissance de la turbine (2), 30 Hc est la puissance du compresseur (3), N est le régime du turbocompresseur (1), est l'opérateur de dérivation par rapport à la variable temps et, J est une constante égale au moment d'inertie du 35 turbocompresseur (1). d dt 9. Procédé selon l'une quelconque des revendications 1 à 8, où le rapport de détente (Rt) de la turbine (2) est calculé en fonction du débit corrigé (Qt cor) de gaz d'échappement (7) au travers de la turbine (2) au moyen d'une fonction (f4) du débit corrigé (Qtcor) de gaz d'échappement (7) au travers de la turbine (2), ladite fonction (f4) étant définie par une cartographie à une dimension. 10. Procédé selon la revendication 9, où le débit corrigé 10 (Qt cor) de gaz d'échappement (7) au travers de la turbine (2) est calculé selon la formule : Qt cor =Q t T t , où Put( \n ù 1) Qt cor est le débit corrigé de gaz d'échappement (7) au travers de la turbine (2), Qt est le débit de gaz d'échappement (7) au travers de la 15 turbine (2), Tut est la température en amont de la turbine (2), Put est la pression en amont de la turbine (2), l'indice (n-1) indiquant ici qu'elle est déterminée à l'intervalle de temps (n-1) précédent. 20 Rt(nù1)i i Qt est le débit de gaz d'échappement (7) au travers de la 25 turbine (2), Ht est la puissance de la turbine (2), nt est le rendement de la turbine (2), Tut est la température en amont de la turbine (2), Rt est le rapport de détente de la turbine (2), l'indice (n-30 1) indiquant ici qu'il est déterminé à l'intervalle de temps (n-1) précédent, Cpt est une première constante thermodynamique du gaz d'échappement (7), yt est une seconde constante thermodynamique du gaz 11. Procédé selon la revendication 10, où le débit (Qt) de gaz d'échappement (7) au travers de la turbine (2) est 11, calculé selon la formule : Qt= t où 1 CPt11tTw 1d'échappement (7). 12. Procédé selon l'une quelconque des revendications 1 à 8, où le rapport de détente (Rt) de la turbine (2) est calculé en fonction de la puissance (Ht) de la turbine (2), du débit (Qt) de gaz d'échappement (7) au travers de la turbine (2), du rendement (rit) de la turbine (2), de la température (Tut) en amont de la turbine (2), selon la formule ùYr ~ytù1 Rt= 1 Qt (n -1)Cptr1 tT t Rt est le rapport de détente de la turbine (2), Ht est la puissance de la turbine (2), Qt est le débit de gaz d'échappement (7) au travers de la turbine (2), l'indice (n-1) indiquant ici qu'il est déterminé à l'intervalle de temps (n-1) précédent, rit est le rendement de la turbine (2), Tut est la température en amont de la turbine (2), Cpt est une première constante thermodynamique du gaz d'échappement (7), yt est une seconde constante thermodynamique du gaz 20 d'échappement (7). 13. Procédé selon la revendication 12, où le débit (Qt) de gaz d'échappement (7) au travers de la turbine (2) est calculé en fonction du débit corrigé (Qt cor) de gaz 25 d'échappement (7) au travers de la turbine (2) selon la Pt(n -l) formule : Qt (n -1) = Qt cor , ou Tz Qt est le débit de gaz d'échappement (7) au travers de la turbine (2), l'indice (n-1) indiquant ici qu'il est déterminé à l'intervalle de temps (n-1) précédent, 30 Qt cor est le débit corrigé de gaz d'échappement (7) au travers de la turbine (2), Put est la pression en amont de la turbine (2), l'indice (n-1) indiquant ici qu'elle est déterminée à l'intervalle de temps (n-1) précédent, et 35 Tut est la température en amont de la turbine (2). 17 Hz , où 14. Procédé selon la revendication 13, où le débit corrigé (Qt cor) de gaz d'échappement (7) au travers de la turbine (2) est calculé en fonction du rapport de détente (Rt) de la turbine (2) au moyen d'une fonction (f5) du rapport de détente (Rt) de la turbine (2), ladite fonction (f5) étant définie par une cartographie à une dimension. 15. Procédé selon l'une quelconque des revendications 11 à 14, où la première constante thermodynamique (Cpt) du gaz d'échappement (7) est égale à 1136 J/kg/K, et où la seconde constante thermodynamique (yc) du gaz d'échappement (7) est égale à 1,34. 16. Procédé selon l'une quelconque des revendications 1 à 15, où le rendement (rit) de la turbine (2), est calculé en fonction du régime corrigé (Ncor) du turbocompresseur (1) et du rapport de détente (Rt(n-1)) de la turbine (2) déterminé à l'intervalle de temps (n-1) précédent, au moyen d'une fonction (f3) du régime corrigé (Ncor) du turbocompresseur (1) et du rapport de détente (Rt) de la turbine (2), ladite fonction (f3) étant définie par une cartographie à deux dimensions. 17. Procédé selon l'une quelconque des revendications 1 à 16, où la pression (Put) en amont de la turbine (2) est calculée selon la formule : Pz = Pd,R, , où Put est la pression en amont de la turbine (2), Pdt est la pression en aval de la turbine (2) et 30 Rt est le rapport de détente de la turbine (2). 18. Procédé selon l'une quelconque des revendications précédentes où le débit (Qc) d'air d'admission (5) au travers du compresseur (3), la pression (Pdc) en aval du compresseur 35 (3) et la température (Tut) en amont de la turbine (2) sont mesurés par des capteurs, et où la pression (Puc) en amont du compresseur (3), la température (Tuc) en amont du compresseur (3) et la pression (Pdt) en aval de la turbine (2) sontdéterminées par un estimateur. 19. Dispositif comprenant des moyens aptes à mettre en œuvre le procédé selon l'une quelconque des revendications 5 précédentes.