FR 2916231 A1 20081121 FR 0755139 A 20070518 DEMANDE DE BREVET B07-0137FR MSA PJ7250 Société par Actions Simplifiée dite : RENAULT s.a.s. Procédé et dispositif d'adaptation d'un estimateur de température d'un système de post-traitement des gaz d'échappement Invention de : BARRILLON Pascal POILANE Emmanuel Procédé et dispositif d'adaptation d'un estimateur de température d'un système de post-traitement des gaz d'échappement L'invention concerne le domaine des systèmes de post traitement, et plus particulièrement, l'estimation des températures dans les systèmes de post traitement. Les moteurs à combustion interne produisent des gaz d'échappement qui contiennent des espèces polluantes, telles que des oxydes d'azote, des hydrocarbures imbrûlés, du monoxyde de carbone et des particules de suie. Les taux de rejet admis de ces différentes espèces polluantes sont strictement réglementés et sujets à révisions périodiques à la baisse. Pour abaisser les émissions de ces différentes espèces polluantes, on utilise depuis de nombreuses années des systèmes de post traitement insérés dans la ligne d'échappement. Les systèmes de post traitement comprennent un matériau catalytique apte à décomposer les différentes espèces polluantes en composés neutres pour la santé, et à stocker les suies pour combustions périodiques. Ces systèmes de post traitement fonctionnent de manière discontinue, alternant les phases de stockage et les phases de décomposition. Pour obtenir une bonne efficacité du système de post traitement, une certaine température est nécessaire. Cette température est atteinte par l'augmentation de la chaleur déposée par le flux de gaz d'échappement traversant le système de post traitement, par l'énergie libérée lors de la décomposition des espèces polluantes et par l'utilisation de l'inertie thermique du système de post traitement. Le système de post traitement peut également être pollué par les oxydes de soufre. La régénération par désulfatation requiert des températures supérieures aux températures de régénération des systèmes de post traitement et sont proches des limites d'endommagement du système de post traitement. Le contrôle de la température du système de post traitement est ainsi crucial, une régulation de la température en sortie n'étant pas suffisante pour contrôler les effets combinés de l'inertie thermique du système de post traitement, de la libération d'énergie thermique lors des réactions de régénération et de l'alternance rapide des modes de combustion nécessaires à la régénération. La demande de brevet FR 2 853 691 décrit une méthode de modélisation de la température d'un piège catalytique en fonction des paramètres de fonctionnement du moteur, de l'environnement du véhicule et du piège catalytique lui-même. Cependant cette demande de brevet présente plusieurs inconvénients. D'une part, elle est limitée aux systèmes de traitement des oxydes d'azote. D'autre part, les calculs à réaliser et les paramètres à calibrer sont nombreux et complexes, nécessitant une forte puissance de calcul et une grande quantité de mémoire. De plus, aucun moyen n'est prévu pour estimer et corriger les dérives et dispersions pouvant apparaître lors de l'alternance des modes de combustion. L'objet de l'invention est d'estimer les dérives des modes de combustion utilisés lors de la régénération du système de post traitement. Un autre objet de l'invention est de corriger des paramètres des modes de combustion afin de tenir compte des dérives déterminées. Selon un mode de mise en oeuvre, on définit un procédé de correction de l'estimation par modèle physique embarqué de l'état thermique d'un système de post traitement implanté dans la ligne d'échappement d'un moteur à combustion interne. Le procédé permet de corriger l'effet des dérives et dispersions des différents organes ayant un impact sur l'état thermique du système de post traitement, en déterminant une énergie potentielle corrigée dépendant des énergies potentielles mesurées et modélisées. On peut déterminer l'énergie potentielle mesurée en réalisant un bilan thermique à partir des mesures de la température en entrée du système de post traitement, de la température en sortie du système de post traitement et du mode de fonctionnement du moteur à combustion interne. On peut déterminer l'énergie potentielle modélisée d'après le mode de fonctionnement du moteur à combustion interne et d'après le modèle physique mémorisé. On peut déterminer une énergie potentielle corrigée en réalisant le produit de l'énergie potentielle modélisée et d'un facteur de correction issu de la somme pondérée du coefficient de correction du modèle physique de l'itération précédente et du rapport entre les énergies potentielles mesurées et modélisées. Selon un mode de réalisation, on définit un dispositif de correction de l'estimation par modèle physique embarqué de l'état thermique d'un système de post traitement implanté dans la ligne d'échappement d'un moteur à combustion interne pouvant fonctionner selon différents modes. Le dispositif de correction comprend un moyen de correction du modèle physique apte à modifier le résultat des calculs issus dudit modèle physique en fonction du signal reçu du moyen de calcul contenant le modèle physique et du signal émis par le moyen de détermination de l'énergie potentielle mesurée. Le dispositif peut comprendre un moyen de calcul contenant un modèle physique de l'énergie libérée lors de la régénération capable de déterminer l'énergie libérée lors de la régénération en fonction du mode de fonctionnement du moteur à combustion interne. Le dispositif peut également comprendre un capteur de température en amont du système de post traitement. Le dispositif peut comprendre en outre un capteur de température en aval du système de post traitement. Le moyen de détermination de l'énergie potentielle mesurée peut recevoir des signaux des capteurs de température en amont et en aval du système de post traitement afin de déterminer l'énergie libérée lors de la régénération en fonction de la température en amont et en aval du système de post traitement. Le système de post traitement peut comprendre un piège catalytique. Selon un autre aspect de l'invention, on définit un mode d'utilisation d'un dispositif pour diagnostiquer le vieillissement du système de post traitement et/ou la dérive des modes de fonctionnement du moteur. D'autres buts, caractéristiques et avantages de l'invention apparaîtront à la lecture de la description suivante, donnée uniquement à titre d'exemple non limitatif et faite en référence aux dessins annexés sur lesquels : -la figure 1 illustre le système d'adaptation d'un estimateur de températures et l'estimateur de température associé. - la figure 2 illustre la modélisation en réacteurs parfaitement agités du système de post traitement ainsi que le placement des capteurs de température. Sur la figure 2, les gaz d'échappement issus de la combustion, pénètrent dans le système de post traitement 5 par une conduite d'échappement 6 et sont évacués par une autre conduite 7. Un premier capteur de température 8 est disposé en amont du système de post traitement 5 et est relié au système de correction par la connexion 8a. Un deuxième capteur de température 9 est disposé en aval du système de post traitement 5 et est relié au système de correction par la connexion 9a. Lors du fonctionnement normal du moteur à combustion interne, le système de post traitement 5 piège différentes espèces polluantes. Périodiquement une régénération du système de post traitement est nécessaire afin de convertir les espèces polluantes en espèces neutres pour la santé, et pour régénérer les propriétés de stockage du système de post traitement 5. Lors de cette régénération, une alternance de deux modes de combustion est déclenchée afin d'élever la température à l'intérieur du système de post traitement 5 de façon à déclencher les réactions de régénération. Les réactions de régénération permettent de transformer les espèces polluantes piégées dans le système de post traitement 5 en espèces moins nocives qui peuvent alors être évacuées avec les gaz d'échappement. Suivant le type de réaction catalytique à déclencher, la température requise est plus ou moins élevée. Un endommagement du système de post traitement peut intervenir si la température dépasse certains seuils, par exemple une température supérieure à 750 C. Pour éviter cela, des moyens de calcul estiment la température à l'intérieur du système de post traitement 5 en fonction de la composition des gaz d'échappement déterminée pour le mode de fonctionnement courant. Le système de détermination de la température, décrit sur la figure 1, comprend un moyen 1 d'estimation des espèces polluantes, un moyen 2 de détermination de l'énergie libérée lors de la régénération, un moyen 3 de calcul de la répartition de l'énergie libérée entre les différents réacteurs et un moyen 4 de détermination de la température de sortie des différents réacteurs. Le moyen 1 d'estimation des espèces polluantes comprend un moyen de modélisation le recevant en entrée des paramètres caractérisant le fonctionnement du moteur à combustion interne, tels que le motif d'injection par la liaison la, le couple délivré par la liaison lb, le régime moteur par la liaison le et le mode de combustion par la liaison ld. Le mode de combustion représente un mode de fonctionnement particulier du moteur. On peut citer par exemple, un mode de régénération du filtre à particule ou un mode à mélange enrichi. Un mode de combustion peut comprendre une ou plusieurs phases requérant, par exemple, des quantités de carburant différentes, des nombres et des fréquences d'injections de carburant différentes. Ces paramètres forment un motif d'injection. La cartographie le des émissions polluantes peut prendre également en compte la quantité de carburant ou de réducteurs injectés dans la ligne d'échappement en aval du collecteur d'échappement. Le moyen de modélisation le émet en sortie par la liaison 1f la nature et le taux d'émission des différentes espèces produites lors de la combustion. Le moyen 2 de détermination de l'énergie libérée lors de la régénération comprend un moyen 21 de détermination des espèces limitantes, un moyen 24 de détermination de l'équivalent énergétique, un moyen 26 de calcul de l'énergie totale libérée, un moyen 29 de correction du modèle physique et un moyen de calcul 28. Le moyen 21 de détermination des espèces limitantes reçoit sur au moins une de ses entrées la nature et le taux d'émission des espèces polluantes par la liaison 1f provenant de la cartographie le des émissions polluantes. Lors de la régénération, les réactions chimiques mettent en jeu des quantités déterminées de réactifs. Un ou plusieurs de ces réactifs peuvent ne pas être présents en quantité suffisante pour que tous les réactifs soient consommés. On qualifie d'espèce limitante les réactifs présents en quantité insuffisante par rapport aux autres espèces, compte tenu de la réaction chimique. Le moyen 21 de détermination des espèces limitantes émet sur au moins une de ses sorties la nature de l'espèce limitante lors de la régénération par la liaison 23 et la quantité de cette espèce limitante qui va réagir par la liaison 22. Les liaisons 22 et 23 sont reliées au moyen 24. Le moyen 24 de détermination de l'équivalent énergétique émet sur au moins une de ses sorties, par la liaison 25 à destination du moyen 26, la quantité équivalente d'espèce limitante correspondant aux quantités des différentes espèces qui vont réagir compte tenu de la quantité et de la nature de l'espèce limitante. Le moyen 26 de calcul de l'énergie totale libérée détermine l'énergie totale libérée lors de la régénération AHmod et émet la valeur correspondante par au moins une de ses sorties à destination du multiplicateur 297 par la liaison 27. Le moyen de calcul 28 reçoit les mesures de température provenant du capteur 8 par la connexion 8a et du capteur 9 par la connexion 9a. Le moyen de calcul 28 détermine une valeur de l'énergie totale libérée AHmes lors des réactions de régénération à partir des températures mesurées et la transmet par la connexion 292 au diviseur 293. Le diviseur 293 reçoit la valeur AHmod par la dérivation 291 de la connexion 27. Le diviseur 293 détermine le coefficient a issu du rapport entre les valeurs estimées AHmod et mesurées AHmes de l'énergie totale libérée lors des réactions de régénération et émet la valeur correspondante par la connexion 294. Le moyen de calcul 295 détermine la correction mod(i) permettant de corriger les dérives et dispersions en fonction de la valeur de a. Le multiplicateur 297 reçoit sur une entrée la valeur de mod(i) par la connexion 296 et par au moins une autre entrée la valeur estimée AHmod de l'énergie totale libérée lors des réactions de régénération et émet en sortie une valeur corrigée AHcorr de l'énergie totale libérée lors des réactions de régénération par la connexion 298. Le moyen 2 de détermination de l'énergie totale libérée lors de la régénération est capable de déterminer l'énergie totale libérée lors de la régénération AHmod engendré par les réactions d'oxydation et de réduction des émissions polluantes lors de la régénération. Le moyen 2 de détermination peut prendre également en compte la quantité de carburant ou de réducteurs injectés dans la ligne d'échappement en aval du collecteur d'échappement. Le moyen 21 de détermination des espèces limitantes détermine la nature de l'espèce limitante. La décomposition des espèces polluante lors de la régénération fait appel à des réactions chimiques dont la stoechiométrie est connue. I1 est alors possible de déterminer les espèces sous-stoechiométriques. La connaissance des espèces sous-stoechiométriques permet également de déterminer la quantité de l'espèce limitante qui réagira. En fonction des données reçues, le moyen 24 de détermination de l'équivalent énergétique détermine le mode de fonctionnement courant et la quantité des différentes espèces qui vont réagir. Selon la teneur des différentes émissions polluantes, plusieurs modes de fonctionnement existent. Soit l'oxygène n'est pas l'espèce limitante, et le mode de fonctionnement est dit mode pauvre. Soit l'oxygène est l'espèce limitante, et le mode de fonctionnement est dit mode riche. Le moyen 24 de détermination estime alors une masse de l'espèce limitante dont la régénération libérera une quantité d'énergie comparable à la quantité d'énergie libérée lors de la régénération des différentes espèces qui vont réagir. La masse d'espèce limitante ainsi déterminée est appelée masse équivalente. Le moyen 24 émet la valeur de la masse équivalente vers le moyen 26 de calcul de l'énergie totale libérée. Le moyen 26 de calcul de l'énergie totale libérée comprend une zone mémoire comprenant les chaleurs de réaction des différents polluants pouvant être émis par le moteur à combustion interne. Selon la nature de l'espèce limitante reçue du moyen 24, le moyen 26 détermine l'énergie totale libérée AHmod. Ainsi, selon le mode de fonctionnement, par exemple en mode pauvre ou en mode riche, seule la valeur de la chaleur de réaction de l'espèce limitante change. Une seule constante est mémorisée, limitant l'espace mémoire nécessaire. En mode de fonctionnement pauvre, la puissance Wmod et l'énergie totale AHmod libérées lors de la régénération sont exprimées en fonction du débit massique et du pouvoir calorifique de l'espèce limitante, par exemple du monoxyde de carbone. Wmod = dmCOeq . PCI CO On obtient l'énergie en intégrant la puissance sur une durée, entre tO et tf. f f AHmod = J Wmod dt = J dmcoeq.PCI CO dt t0 t0 avec dmcoeq = débit massique de monoxyde de carbone équivalent au débit des différentes espèces polluantes PCIco = pouvoir calorifique du monoxyde de carbone lors de la régénération En mode de fonctionnement riche, la puissance Wmod et l'énergie totale AHmod libérée lors de la régénération sont exprimées en fonction du débit massique et de la chaleur de réaction de l'oxygène pour une composition des gaz donnée. Wmod = dm02eq . PCIO2 tf tf AHmod Wmod dt = J dm02eq •PCIO2 • dt t0 t0 avec dmO2 = débit massique de l'oxygène CRO2 = chaleur de réaction de l'oxygène pour une composition des gaz donnée L'énergie totale libérée AHmod permet de déterminer les températures en plusieurs points du système de post traitement. Cependant, des dérives et des dispersions des modes de combustion peuvent apparaître au cours de la vie du véhicule. Suivant leur ampleur, ces dérives et dispersions peuvent fausser l'estimation des températures dans le système de post traitement jusqu'au point où une partie des températures locales du système de post traitement dépasse la température d'endommagement, provoquant la ruine du système de post traitement. Pour compenser les dérives, on réalise un bilan d'énergie à partir de températures mesurées en amont et en aval du système de post traitement par des capteurs de températures 8 et 9, et on recale le modèle sur la situation observée au travers des capteurs de température. Bien que le taux d'échantillonnage des capteurs de température soit élevé devant la dynamique du système, le système fonctionne sur des quantités d'énergie et non des puissances instantanées de façon à réduire l'influence du bruit de mesure. Le moyen de calcul 28 estime l'énergie totale libérée AHmes lors de la régénération dans le système de post traitement 5, en réalisant le bilan de l'énergie apportée par les gaz d'échappement et l'énergie libérée lors de la régénération. Le bilan de conservation des puissances thermiques peut s'écrire de la façon suivante : Pr + Pe = Ps + APcata avec Ps=Qech(t)•CPgaz(t)•Ts(t) : la puissance véhiculée par les gaz sortant du système de post traitement 5 Pr=Wmod(t) : la puissance dégagée par les réactions de régénération Pe=Qech(t).Cpgaz(t).Te(t) : la puissance véhiculée par les gaz entrant dans le système de post traitement 5 APcata = 1 / dt • (meata (CPcata (t + dt)Tcata (t + dt) ù CPcata (t)Tcata (t)) + ... mgaz (CPgaz (t + dt)Tcata (t + dt) ù CPgaz (t)Tcata (t) la variation de puissance dans le système de post traitement 5 Qech = débit massique des gaz d'échappement Cpgaz= Capacité calorifique massique à pression constante des gaz d'échappement Tata = Température du système de post traitement Te= Température des gaz en entrée du réacteur courant RPAi mcata= Masse de système de post traitement dans le réacteur courant RPAi Cpcata= Capacité calorifique massique à pression constante du système de post traitement En négligeant la masse des gaz devant la masse du système de post traitement, on obtient : APcata = 1 / dt • meata (CPcata (t + dt)Tcata (t + dt) ù CPcata (t)Tcata (t)) Dans ce bilan d'énergie, les pertes sont négligées. En replaçant, on obtient : Wexo (t) _ Q (t)c (t)T (t) + w (t) ù Q (t)c (t)T (t) + m (CPcata (t + dt)Tcata (t + dt) ù CPcata (t)Tcata (t)) (1) ech Pgaz s pertes ech Pgaz e tata dt Pour les deux modes de fonctionnement, le moyen de calcul 28 calcule l'énergie potentielle d'après les mesures des capteurs de température 8 et 9. Pour cela, le moyen de calcul 28 intègre la relation (1) entre tO et tf, généralement la durée entre deux mesures de température. On obtient alors les expressions suivantes. OHmes = Sirex, ' dt t0 f f OHmes = mcata ' (Cpcata (tf ). Tcata (tf )ù Cpcata (t0) . Tcata (t0 ))+ f Qech ' C1 gaz . Ts . dt ù JQech Cpgaz . Te . dt (2) t0 t0 Le bilan d'énergie potentielle (2) permet d'obtenir une valeur d'énergie potentielle dite mesurée AHmes, c'est-à-dire estimée d'après les mesures de température effectuées sur le système. Afin de déterminer la présence d'une dérive ou d'une dispersion des modes de combustion, le moyen de calcul 293 compare l'énergie potentielle mesurée AHmes à l'énergie potentielle calculée AHmod par le modèle en déterminant un coefficient a reflétant l'erreur sur l'énergie libérée calculée due à une dérive des injecteurs. Le coefficient a peut être déterminé par exemple en réalisant le rapport des deux énergies. a = Mimes AHmod tf tf tf mcata (Cpcata (t f ) Tcata (t f ) ù Cpcata (t0 ) Tcata (t0 ))+ f Qech Cpgaz dt ù f Qech Cpgaz Te • dt + J Wpertes.dt t0 t0 t0 -"mod Si a n'est pas égal à un, dans une certaine gamme de tolérances, une dérive est présente. a peut alors être utilisé pour recaler le modèle. Pour cela, le moyen de calcul 295 définit une relation de récurrence définissant un coefficient Mod(i). Mod(i) =(1ù pond) . Mod (i -1) + pond • a (i) avec i = itération de calcul en cours pond = coefficient de relaxation Le coefficient Mod(i) obtenu est utilisé par le moyen de calcul 297 pour corriger l'estimation de la puissance libérée lors de la régénération Wexo. Wexo corn )ù cxo(i) Mod(i) En mode de fonctionnement pauvre, la puissance totale libérée lors de la régénération après correction sera Wexocorr (i) = Mod(i) • dmcoeq . PCICO En mode de fonctionnement riche, la puissance totale libérée lors de la régénération après correction sera : Wexo corn (i) = Mod(i) • dmO2 • PCIO2 En mode de fonctionnement pauvre, l'énergie libérée 25 s'écrit alors : tf tf 12 OHcorr = f to = f Mod(i) • dmcoeq •PCICO • dt to En mode de fonctionnement riche, l'énergie libérée s'écrit de même : tf tf 4Hcorr = f Wexo_corr . dt = f Mod(i) • dmO2.PCIO2 • dt to to Afin de modéliser le champ de température dans le système de post traitement, on modélise ledit système en une succession de réacteurs parfaitement agités (RPAi) comme on peut le voir sur la figure 2. Ces réacteurs sont disposés en série, la température de sortie d'un réacteur étant la température d'entrée du réacteur suivant. On modélise généralement un système de post traitement avec moins d'une dizaine de réacteurs parfaitement agités. I1 est à noter que ce nombre n'est dépendant que de la puissance de calcul qui peut être investie dans la modélisation. Le dispositif 3 de répartition de l'énergie totale libérée comprend un moyen 3a de répartition de l'énergie totale libérée dans le système de post traitement. Le moyen 3a de répartition de l'énergie totale libérée reçoit en entrée la valeur corrigée de l'énergie totale libérée AHcorr par la connexion 298 et émet la fraction gexo1 de la quantité d'énergie totale libérée AHcorr correspondant au premier réacteur parfaitement agité par la liaison 3b, la fraction gexo_2 correspondant au deuxième réacteur parfaitement agité par la liaison 3c et la fraction gexo_3 correspondant au troisième réacteur parfaitement agité par la liaison 3d. La répartition de la quantité totale d'énergie libérée est déterminée empiriquement. Elle correspond à une valeur gexo1 égale à 80% de AHcorr, une valeur gexo_2 égale à 20% de AHcorr et des valeurs de gexoi égales à 0% de AHcorr pour i supérieur ou égal à 3. La répartition est donnée ici à titre d'exemple, d'autres répartitions seraient possibles. Un nombre différent de réacteurs pourrait également être pris en compte. Le moyen 4 de détermination de la température de sortie des différents réacteurs comprend un moyen de modélisation pour chaque réacteur parfaitement agité. Dans l'exemple illustré sur la figure 1, la modélisation du réacteur RPA1 est assurée par un moyen de modélisation 4c, la modélisation du réacteur RPA2 par un moyen de modélisation 4e et la modélisation du réacteur RPA3 par un moyen de modélisation 4g. Le moyen de modélisation 4c du réacteur parfaitement agité RPA1 reçoit en entrée le débit massique des gaz d'échappement Qech par la liaison 4a, la température des gaz en entrée Te du système de post traitement par la liaison 4b et la fraction de l'énergie totale libérée gexo1 à dissiper dans le réacteur parfaitement agité RPA1 par la liaison 3b. Le débit massique des gaz d'échappement Qech et la température Te sont déterminés par des moyens d'estimation non représentés. Le moyen de modélisation 4c émet en sortie par la liaison 4d, la valeur de la température en sortie Tsl du réacteur parfaitement agité RPA1. La même description s'applique aux deux autres réacteurs. Le moyen de modélisation 4e du réacteur parfaitement agité RPA2 reçoit en entrée le débit des gaz d'échappement Qech par une dérivation de la liaison 4a, la température des gaz Tsl en sortie du réacteur parfaitement agité RPA1 par la liaison 4d et la fraction de l'énergie totale libérée gexo_2 à dissiper dans le réacteur parfaitement agité RPA2 par la liaison 3c. Le moyen 4e émet en sortie par la liaison 4f, la valeur de la température en sortie Ts_2 du réacteur parfaitement agité RPA2. Le moyen de modélisation 4g du réacteur parfaitement agité RPA3 reçoit en entrée le débit des gaz d'échappement Qech par une dérivation de la liaison 4a, la température des gaz Ts_2 en sortie du réacteur parfaitement agité RPA2 par la liaison 4f et la fraction de l'énergie totale libérée gexo_3 à dissiper dans le réacteur parfaitement agité RPA3 par la liaison 3d. Le moyen 4g émet en sortie par la liaison 4h, la valeur de la température en sortie Ts_3 du réacteur parfaitement agité RPA3. Pour calculer la température de sortie Ts_i d'un réacteur parfaitement agité courant RPAi, le moyen de modélisation réalise d'abord le bilan de l'énergie thermique dans le réacteur. Teata (t -1) . meata - Cpeata + Te(t) - Qech(t) - Cpgaz - dt + gexo_i(t) - dt Ts_i(t)= meata Cpeata + Qech(t) Cpgaz dt (3) Avec mcata= Masse de système de post traitement dans le réacteur courant RPAi gexo i= Fraction de l'énergie totale libérée à dissiper sur le réacteur courant RPAi Les différents termes de l'équation (3) peuvent être identifiés de la façon suivante. Le terme T ata(t -1)• m • Cpeata est l'énergie emmagasinée dans le système de post traitement lors de l'itération précédente. Le terme Te(t) • Qech(t) •Cpgaz • dt est la quantité de chaleur véhiculée par les gaz d'échappement. Le terme q,o i(t)•dt est la quantité de chaleur dégagée lors de la régénération dans le réacteur courant RPAi, gexoi étant la fraction de l'énergie totale libérée AHcorr à dissiper dans le réacteur parfaitement agité courant RPAi. Enfin, le terme m• Cpcata +Qech(t) • Cpgaz • dt est la capacité calorifique d'un réacteur courant RPAi, c'est-à-dire la somme du terme m • Cpcata représentant la capacité calorifique du système de post traitement 5 et du terme Qech(t) •Cpgaz • dt représentant la capacité calorifique des gaz d'échappement contenus dans le réacteur courant RPAi. La température de sortie Ts_i d'un réacteur courant RPAi est la température d'entrée Te_i+l d'un réacteur RPAi+l suivant le réacteur courant RPAi. Le premier réacteur de la série RPA1 utilise comme température d'entrée Te, la température des gaz en entrée du système de post traitement. Le dernier réacteur de la série permet d'obtenir la température de sortie Ts du système de post traitement 5. Alternativement, la température de sortie Ts du système de post traitement 5 peut être déterminée en prenant en compte les pertes thermiques du système de post traitement 5. On ajoute alors le terme de pertes suivant : Text ' Sext • Hext (Vveh)• dt avec Hext(Vveh) = un facteur de convection externe väeh = la vitesse du véhicule Text = la température extérieure Sext = la surface de contactavec l'extérieur Ts(t) _ m • Cpcata + Qech(t) • Cpgaz • dt + Sext • Hext wveh Cependant, la correction de décalage n'est pas appliquée aux deux modes de fonctionnement utilisés lors d'une régénération d'un élément catalytique. En effet, les dispersions et dérives constatées en sortie de l'élément catalytique peuvent avoir plusieurs sources plus ou moins bien identifiées. Cependant la correction a pour but de s'assurer que l'estimation de la température de l'élément catalytique par un modèle physique est conforme aux températures réellement présentes dans l'élément catalytique. Pour cela, on part du principe qu'une seule variable du modèle présente des dispersions et des dérives, les autres étant considérées comme conformes aux mesures. La variable prise en compte est alors corrigée de façon à ce que le résultat de la 15 modélisation soit conforme aux mesures effectuées, même si la variable corrigée ne nécessite pas de correction. De fait, elle devient une variable d'ajustement du modèle. Dans le cas présent, un des deux modes de combustion devient la variable d'ajustement du modèle. Un contrôle permanent de ce mode de combustion est assuré par mesure et 20 afin de permettre l'ajustement du modèle. Cela est facilité par le fait que les deux modes de combustion se succèdent à un rythme élevé, de l'ordre de quelques dizaines de secondes. Un effet de moyenne entre les deux modes apparaît avant la prise de mesure. L'ajustement du modèle est alors sensible aux effets combinés des deux modes de 25 combustion et reflète donc bien l'état thermique de l'élément catalytique. Le dispositif d'adaptation peut être appliqué à la compensation des dispersions de fabrication ou à la correction d'une dérive relativement lente du système. Mais le dispositif d'adaptation pourrait 30 également être utilisé pour le diagnostic de l'état de l'élément catalytique, un élément catalytique en fin de vie ayant une efficacité moindre et générant donc une quantité moindre d'énergie lors de la régénération. Dans ce cas, l'équation (3) permettant de déterminer la température du système devient : Tcata m • Cpcata + Te(t) • Qech(t) •Cpgaz • dt +T ext ' Sext ' Hext (vveh )' dt + gexo (t) dt 10 REVENDICATIONS 1. Procédé de correction de l'estimation par modèle physique de l'état thermique d'un système de post traitement (5) implanté dans la ligne d'échappement d'un moteur à combustion interne, caractérisé par le fait que l'on corrige l'effet des dérives et dispersions des différents organes ayant un impact sur l'état thermique du système de post traitement (5), en déterminant une énergie potentielle corrigée dépendant des énergies potentielles mesurées et modélisées. 2. Procédé selon la revendication 1, dans lequel on détermine l'énergie potentielle mesurée en réalisant un bilan thermique à partir des mesures de la température en entrée du système de post traitement (5), de la température en sortie du système de post traitement (5) et du mode de fonctionnement du moteur à combustion interne. 3. Procédé selon l'une des revendications précédentes, dans lequel on détermine l'énergie potentielle modélisée d'après le mode de fonctionnement du moteur à combustion interne et d'après le modèle physique mémorisé. 4. Procédé selon l'une des revendications précédentes, dans lequel on détermine une énergie potentielle corrigée en réalisant le produit de l'énergie potentielle modélisée et d'un facteur de correction issu de la somme pondérée du coefficient de correction du modèle physique de l'itération précédente et du rapport entre les énergies potentielles mesurées et modélisées. 5. Dispositif de correction de l'estimation par modèle physique de l'état thermique d'un système de post traitement implanté dans la ligne d'échappement d'un moteur à combustion interne pouvant fonctionner selon différents modes, caractérisé par le fait qu'il comprend un capteur (8) de température en amont du système de post traitement (5), un capteur (9) de température en aval du système de post traitement (5), un moyen (28) de détermination de l'énergie potentielle mesurée et un moyen (29) de correction du modèle physique apte à modifier le résultat des calculs issus dudit modèle 2916231 l$ physique en fonction du signal reçu du moyen (26) de calcul contenant le modèle physique, du moyen (28) de détermination de l'énergie potentielle mesurée et des capteurs de température (8,9). 6. Dispositif selon la revendication 5 comprenant un moyen 5 (26) de calcul contenant un modèle physique de l'énergie libérée lors de la régénération capable de déterminer l'énergie libérée lors de la régénération en fonction du mode de fonctionnement du moteur à combustion interne. 7. Dispositif selon l'une des revendications 5 ou 6, dans lequel 10 le moyen (28) de détermination de l'énergie potentielle mesurée reçoit des signaux des capteurs de température en amont (8) et en aval (9) du système de post traitement afin de déterminer l'énergie libérée lors de la régénération en fonction de la température en amont et en aval du système de post traitement (5). 15 8. Dispositif selon l'une des revendications 5 à 7, dans lequel le système de post traitement comprend un piège catalytique. 9. Utilisation d'un dispositif selon l'une des revendications 5 à 8 pour diagnostiquer le vieillissement du système de post traitement (5) et/ou la dérive des modes de fonctionnement du moteur.