Brevet FR2948148A1 — publié le 20110121

Le domaine d'application principal de l'invention est le traitement des gaz d'échappement des moteurs Diesel. En particulier la méthode concerne le contrôle de la température en aval d'un catalyseur d'oxydation. En particulier le contrôle de cette température sert à assurer la régénération active d'un filtre à particules. 15 État de la technique Tel qu'on le sait dans l'art antérieur, les véhicules Diesel soumis aux normes sur les émissions polluantes sont ou vont être équipés de filtres à particules pour satisfaire ces normes. Les suies s'accumulent sur les parois du filtre à particules. Pour assurer la bonne circulation des gaz dans le système d'échappement et donc le bon fonctionnement du 20 moteur, ces suies doivent être brûlées régulièrement lors d'une période de � régénération active �. Cette régénération active se fait en élevant la température du filtre à particules. Dans la plupart des configurations, cette élévation de température est réalisée à l'aide de l'oxydation d'hydrocarbures injectés dans le système d'échappement : un catalyseur d'oxydation est placé en amont d'un filtre à particules ; les hydrocarbures peuvent 25 notamment être fournis par le système d'injection principal du moteur ou par un injecteur spécifique localisé entre le collecteur d'échappement et l'entrée du catalyseur d'oxydation. Dans le système de post-traitement des moteurs Diesel comprenant un filtre à particules, il est très important de bien maîtriser la température du gaz entrant dans le filtre à particules sous peine de risquer de détériorer ce filtre. II est donc très important de maîtriser 30 la température du gaz sortant du catalyseur d'oxydation. Le problème technique est le suivant. Il s'agit de maintenir, pendant les périodes de régénération active, la température du catalyseur d'oxydation proche d'une température de 1 consigne dans toutes les conditions d'utilisation du véhicule. Lors de l'utilisation du véhicule, les conditions externes imposées au catalyseur d'oxydation sont rapidement et fortement variables. Pour cette raison, il est difficile de résoudre le problème mentionné ci-dessus. Des méthodes actuellement utilisées existent pour résoudre ce problème. La méthode la plus couramment utilisée consiste à calibrer un contrôleur classique (PI, PID, prédicteur de Smith, modèle interne, ...). Cette calibration est effectuée pour une combinaison importante de conditions opératoires stationnaires. Lors de l'utilisation du contrôle, les conditions opératoires sont détectées et les paramètres correspondants du régulateur sont appliqués. A cette méthode, un terme de précompensation est couramment ajouté. Ce terme est évalué notamment en fonction de la différence entre la température d'entrée et la température de consigne, comme décrit dans le document US 7,047,729 B2. Cependant, les méthodes actuelles présentent certaines limitations de performance car : elles ne tiennent pas compte finement des différences phénoménologiques entre l'effet des hydrocarbures sur le catalyseur d'oxydation et l'effet de la température en amont du catalyseur d'oxydation ; elles prennent en compte de façon limitée l'aspect dynamique des effets liés aux variations de débit de gaz ; elles prennent en compte de façon limitée les variations de débit de gaz pour l'adaptation dynamique des paramètres du contrôleur. De plus, les méthodes actuelles présentent des limitations en termes d'effort de calibration, car elles n'expriment pas les paramètres du contrôleur à partir de grandeurs issues de la modélisation physique du catalyseur d'oxydation. L'objet de l'invention est un procédé alternatif de traitement de gaz d'échappement d'un moteur à combustion interne, dans lequel on contrôle la régénération active d'un filtre à particule au moyen d'un contrôle dynamique de la température en sortie d'un catalyseur d'oxydation. Le procédé s'affranchit des limitations des méthodes antérieures au moyen d'une loi de commande comportant : - un terme prenant en compte des phénomènes transitoires ayant lieu dans le catalyseur d'oxydation lors d'une variation de débit de gaz, ou - un terme de précompensation prenant en compte les différences de dynamiques entre l'effet de la perturbation de la température en entrée du catalyseur d'oxydation et l'effet de la commande de l'injection d'hydrocarbures en amont du catalyseur d'oxydation, - ou un terme de rétroaction basé sur une modélisation physique du catalyseur d'oxydation permettant l'expression des paramètres du contrôleur en fonction des paramètres du modèle physique et prenant en compte de manière dynamique l'adaptation des paramètres du contrôleur lors des variations des conditions de fonctionnement, ou toute combinaison de ces trois termes. L'invention concerne également un moteur à combustion interne comprenant un système de contrôle adapté à la mise en oeuvre des procédés selon l'invention. La méthode selon l'invention Un objet de l'invention concerne un procédé de traitement de gaz d'échappement d'un moteur à combustion interne équipé d'un filtre à particules et d'un catalyseur d'oxydation placé en amont du filtre, dans lequel on contrôle la régénération active du filtre au moyen d'un contrôle dynamique d'une température T2 du gaz sortant du catalyseur d'oxydation. Le procédé comporte les étapes suivantes : i- on détermine une température de consigne Tc dudit gaz sortant dudit catalyseur pour assurer une régénération active dudit filtre ; ii- on détermine une demande de température nécessaire pour contrôler ladite température T2 ; iii- on détermine un gain stationnaire G, défini comme le rapport entre la valeur stationnaire d'une variation de T2 et une valeur stationnaire d'une variation de débit massique d'hydrocarbures ayant induit ladite variation de T2 ; iv- on corrige ledit gain G de façon à prendre en compte une variation de débit de gaz en entrée dudit catalyseur, le gain étant transitoirement diminué par rapport à G pour une augmentation de débit de gaz et transitoirement augmenté par rapport à G pour une diminution de débit de gaz ; v- on définit une loi de commande d'un débit d'hydrocarbures en entrée dudit catalyseur d'oxydation, par le rapport entre ladite demande de température et ledit gain corrigé ; vi- on modifie le débit d'hydrocarbures en appliquant la loi de commande pour que la température dudit gaz sortant dudit catalyseur respecte ladite température de consigne. Selon ce procédé, on peut corriger le gain G en calculant une dérivée du débit de gaz traversant le catalyseur, ou alors, en appliquant un retard au gain selon la formule suivante : Ginst(t) = G(t û D(t)) avec Ginst le gain corrigé, D le retard déterminé en fonction du débit de gaz Q et de ses variations, t le temps. Un objet de l'invention concerne un second procédé de traitement de gaz d'échappement d'un moteur à combustion interne équipé d'un filtre à particules et d'un catalyseur d'oxydation placé en amont du filtre, dans lequel on contrôle la régénération active du filtre au moyen d'un contrôle dynamique d'une température T2 du gaz sortant du catalyseur d'oxydation. Le procédé comporte les étapes suivantes : on détermine une demande de température TFF nécessaire pour compenser un effet d'une variation de température Ti dudit gaz en entrée dudit catalyseur, TFF étant une fonction d'une différence entre T1(t û 8(t)) et une valeur de référence Tbar, où T1(t -8(t� est une expression d'une mesure de T1(t) retardée d'un retard de valeur 8(t) ; ii- on détermine un gain stationnaire G, défini comme le rapport entre la valeur stationnaire d'une variation de T2 et une valeur stationnaire d'une variation de débit massique d'hydrocarbures ayant induit ladite variation de T2, tous les autres paramètres restants constants par ailleurs ; iii- on définit une loi de commande uFF d'un débit d'hydrocarbures en entrée dudit catalyseur d'oxydation, par le rapport entre ladite demande de température TFF et ledit gain ; iv- on modifie le débit d'hydrocarbures en appliquant la loi de commande pour que la température dudit gaz sortant dudit catalyseur respecte ladite demande de température TFF. Selon ce second procédé, le retard 8 peut être défini en fonction d'un retard 82 obtenue en résolvant l'équation implicite suivante : t L - Lu = fv(w)dw t-82(t) avec L : une longueur du catalyseur d'oxydation. Lu : une longueur de catalyseur d'oxydation permettant de rapprocher au mieux d'une part une réponse indicielle normalisée à une variation de Ti sur un catalyseur de longueur Lu, et d'autre part une réponse indicielle normalisée à une variation de débit d'hydrocarbures sur un catalyseur d'oxydation de longueur L. v une vitesse moyenne du gaz dans ledit catalyseur d'oxydation. Selon ce second procédé, on peut corriger le gain G de façon à prendre en compte une variation de débit de gaz en entrée dudit catalyseur, le gain étant transitoirement augmenté par rapport à G pour une augmentation de débit de gaz et transitoirement diminué par rapport à G pour une diminution de débit de gaz. Un objet de l'invention concerne un troisième procédé de traitement de gaz d'échappement d'un moteur à combustion interne équipé d'un filtre à particules et d'un catalyseur d'oxydation placé en amont du filtre, dans lequel on contrôle la régénération active du filtre au moyen d'un contrôle dynamique d'une température T2 du gaz sortant du catalyseur d'oxydation. Le procédé comporte les étapes suivantes : i- on détermine une température de consigne Tc dudit gaz sortant dudit catalyseur pour assurer une régénération active dudit filtre ; ii- on détermine une demande de température TFB nécessaire pour que la température T2 dudit gaz sortant dudit catalyseur respecte ladite température de consigne Tc, ladite demande TFB étant calculée à partir d'un contrôleur dont les paramètres sont déterminés au moyen d'un modèle physique du catalyseur d'oxydation, basé sur un bilan d'énergie d'une phase gaz et un bilan d'énergie d'une phase solide, en fonction des conditions opératoires du moteur ; iii- on détermine un gain stationnaire G, défini comme le rapport entre la valeur stationnaire d'une variation de T2 et une valeur stationnaire d'une variation de débit massique d'hydrocarbures ayant induit ladite variation de T2; iv- on définit une loi de commande uFB d'un débit d'hydrocarbures en entrée dudit catalyseur d'oxydation, par le rapport entre ladite demande de température TFF et ledit gain ; v- on modifie le débit d'hydrocarbures en appliquant la loi de commande pour que la température dudit gaz sortant dudit catalyseur respecte ladite température de consigne. Selon ce troisième procédé, on peut corriger le gain G de façon à prendre en compte une variation de débit de gaz en entrée dudit catalyseur, le gain étant transitoirement augmenté par rapport à G pour une augmentation de débit de gaz et transitoirement diminué par rapport à G pour une diminution de débit de gaz. Selon l'invention, la loi de commande appliquée pour modifier le débit d'hydrocarbures peut être définie par la somme des lois de commande uFF et uFB. Ces lois de commandes peuvent être définies à partir d'un gain corrigé de façon à prendre en compte une variation de débit de gaz en entrée dudit catalyseur, le gain étant transitoirement augmenté par rapport à G pour une augmentation de débit de gaz et transitoirement diminué par rapport à G pour une diminution de débit de gaz. Enfin, un autre objet de l'invention concerne un moteur à combustion interne comprenant un système de traitement de gaz d'échappement, un catalyseur d'oxydation (15) placé en amont d'un filtre à particules (16), et un système d'injection (80) agissant dans le moteur. Le moteur comporte un système de contrôle (50) pour piloter le système d'injection de façon à fournir un débit d'hydrocarbures en amont du catalyseur d'oxydation, en accord avec un des procédés précédemment décrits. D'autres caractéristiques et avantages de la méthode selon l'invention, apparaîtront à la lecture de la description ci-après d'exemples non limitatifs de réalisations, en se référant aux figures annexées et décrites ci-après. Présentation succincte des figures La figure 1 décrit un moteur à combustion interne selon l'invention, comprenant un système de traitement de gaz d'échappement. La figure 2 présente sommairement le système à commander et des variables utilisées. 25 La figure 3 est une présentation schématique du contrôleur utilisé par la méthode selon l'invention. La figure 4 illustre schématiquement un contrôleur avec modèle interne MU utilisé alternativement dans la méthode selon l'invention. La figure 5 illustre schématiquement un contrôleur de type prédicteur de Smith 30 représentant un cas particulier du contrôleur présenté sur la figure 4. La figure 6 illustre schématiquement un contrôleur avec modèle interne M utilisé alternativement dans la méthode selon l'invention.20 La figure 7 illustre les différences entre la réponse à une variation de température TI en amont du catalyseur d'oxydation et la réponse à une variation de débit d'hydrocarbures. La figure 8 montre une comparaison d'une réponse à des variations de température Tl sur une longueur fictive et à une réponse des variations de débit d'hydrocarbures. La figure 9 illustre la stratégie de synchronisation selon l'invention. La figure 10 montre l'évolution du rendement de conversion des hydrocarbures en fonction du débit de gaz. La figure 11 illustre l'évolution du débit d'hydrocarbures u à appliquer lors d'une variation de débit de gaz. Description détaillée de la méthode Les notations suivantes sont introduites pour décrire l'invention : TI température représentative de la température du gaz en amont du catalyseur d'oxydation (DOC). T2 température représentative de la température du gaz en aval du catalyseur d'oxydation (en amont du filtre à particule (DPF)). u débit d'hydrocarbures en amont du catalyseur d'oxydation. Il tient compte de tous les hydrocarbures présents. L longueur du catalyseur d'oxydation. Lu longueur de catalyseur d'oxydation permettant de rapprocher au mieux d'une part la réponse indicielle normalisée à une variation de TI sur un catalyseur de longueur Lu, et d'autre part la réponse indicielle normalisée à une variation de débit d'hydrocarbures sur un catalyseur d'oxydation de longueur L. On note que la notion de longueur différente entre L et Lu est équivalente à la notion de temps de réponse différent entre la réponse à une variation de TI et la réponse à une variation de débit d'hydrocarbures. On considère que Lu Selon une seconde alternative, la commande u comporte un terme de précompensation uFF. Ce terme uFF sert à précompenser les variations de température Tl, c'est-à-dire à rendre minimal l'effet des variations de la température T1 sur la température T2. C'est une fonction non linéaire des variations de Tl par rapport à une référence Tbar choisie arbitrairement. Selon une troisième alternative, la commande u comporte un terme de rétroaction uFB. Ce terme uFB sert à faire tendre la température T2 vers la température de consigne Tc par une méthode de rétroaction. Selon une quatrième alternative, la commande u se décompose en la somme du terme de précompensation uFF et du terme de rétroaction uFB. La figure 3 est une présentation schématique du contrôleur utilisé par le procédé de l'invention selon cette alternative. Les trois dernières alternatives peuvent être combinés avec la première alternative. Dans toutes les alternatives, le calcul de la commande u peut être adapté en fonction des conditions opératoires ENG (débit d'air, régime moteur, consigne de couple, pression du rail ...) ou des conditions ambiantes EXT (température extérieure,...). Définition du gain instantané Ginst Selon une première alternative, afin de convertir les demandes de variation de température dans l'unité de commande de l'actionneur (débit massique exprimé par exemple en kg/h), la loi de commande suivante est utilisée : T u= G avec : T représente la demande de température. u représente le débit d'hydrocarbures qu'il faut fournir pour satisfaire la demande de température T. G gain stationnaire, défini comme le rapport entre la valeur stationnaire d'une variation de T2 et la valeur stationnaire de la variation de débit massique d'hydrocarbures u ayant induit cette variation de T2, tous les autres paramètres restants constants par ailleurs. Le gain G peut être évalué expérimentalement ou à partir de l'enthalpie de réaction des hydrocarbures, par exemple comme présenté en annexe 1. A chaque valeur stationnaire du débit de gaz peut correspondre une valeur de ce gain. La valeur calculée classiquement avec le gain stationnaire ne tient pas compte des phénomènes transitoires ayant lieu dans le catalyseur d'oxydation lors des variations de débit de gaz, ce qui provoque des variations importantes de température T2, alors même que la température Tl et que la demande de température T sont constantes. Les méthodes actuelles présentent des limitations de performance car elles prennent en compte de façon limitée l'aspect dynamique des effets liés aux variations de débit de gaz. Selon une première alternative de l'invention, afin de convertir les demandes de variation de température dans l'unité de commande de l'actionneur (débit massique exprimé par exemple en kg/h), le calcul suivant est utilisé : T u= Ginst avec : T représente la demande de température. u : représente le débit d'hydrocarbures qu'il faut fournir pour satisfaire la demande de température T. Ginst : gain instantané à appliquer à l'instant courant. II est basé sur le gain stationnaire G auquel on applique une correction pour pallier les effets des variations de débit de gaz Q. Le principe général de cette correction est que, par rapport à la valeur de gain G calculée classiquement, suite à des augmentations de débit de gaz Q, le gain instantané Ginst est transitoirement augmenté, et en conséquence. Le débit d'hydrocarbures est transitoirement diminué par rapport à la valeur calculée classiquement. Par rapport à la valeur de gain G calculée classiquement, suite à des réductions de débit de gaz Q, le gain instantané Ginst est transitoirement réduit, et en conséquence, le débit d'hydrocarbures est transitoirement augmenté par rapport à la valeur calculée classiquement. Pour réaliser ce principe général, une autre méthode proposée est de retarder l'effet du gain stationnaire selon Ginst(t) = G(t ù D(t)) Ginst gain instantané à appliquer à l'instant courant. D retard à l'application du gain stationnaire. Il est déterminé en fonction du débit de gaz Q et ses variations. G gain stationnaire. G(tùD(t)) représente la valeur de G retardée d'un retard de valeur D(t) Ces deux exemples de réalisation du principe général sont schématisés dans la figure 11 pour le cas d'une augmentation de débit de gaz. La figure 11 illustre l'évolution du débit d'hydrocarbures u à appliquer lors d'une variation de débit de gaz pour une demande de température constante. La valeur ur est le débit d'hydrocarbures appliqué dans le cas classique se basant sur les variations de débit Q par l'intermédiaire du gain stationnaire G. Les valeurs ui et ud sont des exemples de réalisation de la commande u à appliquer selon l'invention, c'est-à-dire par l'intermédiaire du gain instantané Ginst. Dans le cas présenté d'une augmentation de débit de gaz Q, les débits d'hydrocarbures ui ou ud sont inférieurs au débit d'hydrocarbures ur pendant une période transitoire faisant suite à cette variation de débit de gaz Q. Définition du terme de précompensation uFF Le terme de précompensation uFF sert à compenser les effets des variations de température Tl. Le terme uFF représente un débit d'hydrocarbures (HC) nécessaire pour compenser les effets des variations de Tl. Le débit d'hydrocarbure uFF a pour but d'être tel que les variations de T2 engendrées par les variations de TI soient minimales. C'est une fonction non linéaire des variations de TI par rapport à une référence Tbar choisie arbitrairement. La stratégie de précompensation consiste à retarder la compensation d'un effet d'une variation de Tl, dans le but de synchroniser l'effet de TI avec l'effet engendré par une variation du débit d'hydrocarbures uFF. Le calcul du retard de synchronisation est réalisé à débit de gaz variable (Q). En effet, en se basant sur des modèles validés expérimentalement, on peut établir que les réponses à des variations de TI et à des variations de débit d'hydrocarbure sont différentes. La figure 7 illustre ce phénomène. La figure 7 illustre des différences entre les réponses à une variation de température de Tl et à une variation de débit d'hydrocarbures. La valeur TI SR représente la réponse indicielle à une variation de Tl. La valeur HC SR représente la réponse indicielle à une variation de débit d'hydrocarbures. La valeur NT est la variation normalisée de la température T2. La valeur t représente le temps. La stratégie de précompensation de la méthode fonctionne sur le principe que la réponse TI SR est plus lente que la réponse HC SR et consiste à retarder la commande de débit d'hydrocarbures lors de la détection d'une variation de Tl pour synchroniser les effets T1 SR avec ceux HC SR afin de compenser les effets TI SR par ceux de HC SR. La figure 8 montre une comparaison d'une réponse à des variations de température TI sur une longueur fictive et à une réponse à des variations de débit d'hydrocarbures. La valeur AT1 SR représente la réponse indicielle normalisée à une variation de TI pour un catalyseur d'oxydation ayant une longueur adaptée Lu. La valeur HC-SR représente la réponse indicielle normalisée à une variation débit d'hydrocarbures pour un catalyseur d'oxydation ayant une longueur L, L étant supérieure ou égale à Lu. La figure 9 illustre la stratégie de synchronisation selon l'invention. La valeur TI R représente l'effet de la variation de TI en sortie du catalyseur d'oxydation. La valeur N HC représente le débit normalisé des hydrocarbures oxydés dans le catalyseur d'oxydation. La valeur HC R représente l'effet de la variation de ces hydrocarbures N HC. La valeur T2 représente la température du gaz en sortie du catalyseur d'oxydation. Enfin, la valeur 8 représente le retard appliqué sur la commande de débit d'hydrocarbures lors d'une variation de la température Tl. L'utilisation de ce retard permet de rejeter l'effet Tl R de cette variation sur la température T2. Selon l'invention, afin de convertir les demandes de variation de température dans l'unité de commande de l'actionneur (débit massique exprimé par exemple en kg/h), le calcul suivant est utilisé : TFF f(T1ùThar ) uFF = = GFF GFF avec : Tbar une valeur de référence quelconque. Elle peut par exemple être égale à la valeur de TI au moment on l'on active le contrôleur. TFF : demande de variation de température pour compenser la variation de TI. TFF est fonction de la différence entre TI et Tbar GFF : gain instantané à appliquer à l'instant courant pour convertir la demande de variation de température TFF en débit d'hydrocarbures uFF. En particulier, on peut prendre GFF = Ginst. Selon un mode préféré de réalisation, on calcule le terme TFF comme une fonction de la différence entre T1(t ù 8(t)) et Tbar, où T1(t -8(t� est l'expression de la mesure de T1(t) retardée d'un retard de valeur 8(t) . En particulier, TFF se calcule par la formule : TFF = ùKd (T 1(t ù 8(t)) ù Tbar) En particulier le retard 8 se calcule sur la base d'un retard 82. La valeur 82 est obtenue en résolvant l'équation implicite suivante : L-Lu= fv(w)dw t-82(t) Le retard 8est obtenu par une fonction linéaire ou non linéaire de 82. En particulier, on peut exprimer 8 par : 8 = a.82 + b , où a et b sont des constantes à définir. On peut prendre 8= 82. D'autre part, le gain Kd est défini comme le rapport entre la valeur stationnaire d'une variation de T2 et la valeur stationnaire de la variation de Ti ayant induit cette variation de T2. On prend préférentiellement Kd=1. La valeur du gain instantané GFF est la valeur du gain à appliquer à l'instant courant. Cette elle peut être prise égale à la valeur du gain stationnaire G. Alternativement, elle peut être prise égale à la valeur du gain stationnaire G corrigée par une fonction de la dérivée filtrée du débit de gaz traversant le catalyseur, ou à la valeur du gain stationnaire G retardé, comme décrit selon le premier mode de réalisation. Un exemple de méthode pour calculer le gain est présenté en annexe 1. En pratique, on calcule en continu la différence entre Ti et Tbar. Si cette différence est non nulle, on souhaite modifier le débit d'hydrocarbures uFF pour compenser cette variation de température. L'instant de modification de ce débit est cadencé par le pas de calcul du contrôleur du moteur. A chaque pas de calcul, on calcule le retard 8(t) et on cherche quelles variations de TI ont eu lieu à l'instant t-8(t). Ce sont ces variations qu'on veut compenser. Pour cela, on modifie la demande de température TFF (et donc le débit uFF) selon ces variations de TI ayant eu lieu à l'instant t-8(t). Ainsi, lorsque l'on détecte une variation de Tl, elle n'est pas compensée immédiatement. On "attend" un temps 5 (qui n'est pas déterminé à l'instant où la variation est détectée) avant de faire varier la demande de variation de température TFF et donc le débit uFF. L'exemple ci-après illustre le principe du terme de précompensation uFF pour un débit de gaz Q constant: à t=0 s : T1=400°C TFF =O°C uFF=0 kg/h (Thar=T1=400°C) à t=1 s : T1=390°C TFF =0°C uFF=0 kg/h Le calcul du terme TFB est basé sur une modélisation physique du catalyseur d'oxydation. Ce calcul est réalisé par un contrôleur dont les paramètres sont déterminés au moyen d'un modèle physique du catalyseur d'oxydation, basé sur un bilan d'énergie d'une phase gaz et un bilan d'énergie d'une phase solide. Un exemple de modèle physique du catalyseur d'oxydation est présenté ci-après. Un tel modèle se base sur le bilan d'énergie d'une phase gaz et le bilan d'énergie d'une phase solide. Ces bilans d'énergie s'obtiennent en équilibrant les échanges d'énergie dans un volume de contrôle selon le principe suivant : le taux de variation de l'énergie stockée dans le volume de contrôle est égale à la différence entre le taux de variation de l'énergie entrante dans le volume de contrôle et le taux de variation de l'énergie sortante du volume de contrôle. Plusieurs formulations existent. Elles sont passées en revue dans le document suivant : C. Depcik, D. Assanis, One-dimensional automotive catalyst modeling, Progress In Energy Combustion and Science., vol. 31, 2005, pp 308-369. Le principe est le suivant : Le bilan énergétique en phase gaz s'écrit : A+B=C où - A est le terme de stockage d'énergie en phase gaz. Il est conservé dans notre exemple mais peut être très petit et donc parfois négligé. II s'écrit couramment A = Pg CPg e aT B est le terme de transport d'énergie en phase gaz lié au déplacement du gaz lui-même. Il s'écrit couramment B=PgCPgv ôz C est un terme d'échange d'énergie entre la phase gaz et la phase solide. Il s'écrit couramment C = hg Ga (Tm ù T) Le bilan énergétique en phase solide s'écrit : - Am=Cm+Dm+Em+Fm - Am est le terme de stockage d'énergie en phase solide. II s'écrit couramment Am = (1ù e) pm Cpm aTm at Cm est un terme d'échange d'énergie entre la phase gaz et la phase solide. Il s'écrit couramment comme l'opposé du terme C Cm = hg Ga (T ùTm) Dm est un terme d'échange d'énergie par conduction dans la phase solide. Il est négligé dans notre exemple mais peut être pris en compte si nécessaire. Il s'écrit couramment z Tm Dm=(1ùs)km a aZ2 Em est un terme représentant l'enthalpie produite ou absorbée par l'aspect exothermique ou endothermique des réactions se produisant dans le catalyseur d'oxydation. Il s'écrit couramment NM Em = GcaE Rj hj i=l Fm est un terme d'échange thermique avec le milieu extérieur. ll représente couramment les � pertes thermiques �. Il est négligé dans notre exemple mais pourrait être inclus s'il s'avérait nécessaire. Il s'écrit Fm = ùhoc Asoo (Tm ù Too) avec : - T température locale du gaz. Dans le document précité de Depcik, la notation Tg est adoptée à la place de la notation T. - Tm température locale du solide. pg masse volumique du gaz. Pm masse volumique du solide. Cpg capacité calorique du gaz. - Cpm capacité calorique du solide. E fraction de vide du pain catalytique. v vitesse représentative de la vitesse du gaz dans les canaux du pain catalytique. Dans le document précité de Depcik, la notation u est adoptée à la place de la notation v. z variable de la dimension axiale du pain catalytique. Dans le document précité de Depcik, la notation x est adoptée à la place de la notation z. - hg coefficient d'échange convectif d'enthalpie entre le gaz et le solide. Ga coefficient représentatif de la surface d'échange entre le gaz et le solide. km coefficient d'échange conductif du solide. - Gca coefficient représentatif de la surface d'échange entre le gaz et la surface catalytique. - NM nombre d'espèces chimiques prises en considération dans les réactions catalytiques. - Rj taux de réaction de l'espèce j. - hj enthalpie de réaction de l'espèce j. hoo coefficient d'échange convectif d'enthalpie entre le solide et le milieu ambiant. Asoo surface extérieure du solide en contact avec le milieu ambiant. Too température représentative de la température du milieu ambiant. Selon un exemple de réalisation, on utilise un modèle basé sur les équations de bilan précédentes, simplifiées en ne prenant en compte que les termes A, B, C, Am et Cm, puis normalisées de la manière suivante at = k2 (T ûTm) 0T aT at +v ù = ûkl (T ûTm) (sys) : 5Tm avec klù hg Ga s pg Cpg k2 ù hg Ga (1ù e) pm Cpm Ainsi, les paramètres de ce modèle physique simplifié sont : • kl, k2 deux constantes caractéristiques du catalyseur d'oxydation • Les paramètres géométriques du catalyseur permettant le calculer la vitesse v à partir du débit Q • La longueur L. • Le rendement r (qui peut être une fonction de v) • La variable d'adaptation Lu Les entrées de ce modèle physique simplifié sont • Le débit massique de gaz Q. • La température Tl en amont du catalyseur d'oxydation. • Le débit d'hydrocarbures de commande u (débit d'hydrocarbures présents en amont du catalyseur d'oxydation et corrigé du rendement). Le terme Em qui décrit la génération d'enthalpie par l'oxydation des espèces chimiques telles que les hydrocarbures n'est pas pris en compte directement dans le modèle physique simplifié. On note que la prise en compte des hydrocarbures se fait à travers le modèle simplifié mais en utilisant une longueur Lu. L'obtention de cette longueur Lu par la méthode décrite en annexe est similaire à la prise en compte du terme Em. Dans la méthode selon l'invention les paramètres kl et k2 sont calculés ou bien identifiés en utilisant la réponse indicielle expérimentale à une variation de Tl. Cette réponse est donnée de manière analytique par la formule (EQREP) donnée en annexe 1 qui correspond à la résolution du système d'équations (sys) pour une entrée échelon de Tl et pour un catalyseur de longueur L. Cette réponse est schématisée sur la figure 7 en tant que Tl SR. Le paramètre Lu est un paramètre de réglage qui permet d'assimiler la réponse indicielle à HC à une réponse indicielle à Tl sur une longueur fictive Lu. Lu peut notamment être une fonction de v. 18 Le paramètre Lu est utilisé pour le calcul de uFB dans l'adaptation des paramètres des contrôleurs. Une particularité de l'invention consiste à proposer une méthode pour adapter automatiquement les paramètres de ces contrôleurs en fonction des conditions opératoires du moteur. Plusieurs structures classiques de contrôle sont envisagées pour le calcul du TFB : • Le contrôleur de type proportionnel intégral (PI) • Le contrôleur de type intégral simple • Le contrôleur de type modèle interne • Le contrôleur de type prédicteur de Smith Ces méthodes sont détaillées en annexe 2. La figure 4 illustre schématiquement un contrôleur avec modèle interne MU utilisé alternativement dans la méthode selon l'invention. La figure 5 illustre schématiquement un contrôleur de type prédicteur de Smith représentant un cas particulier du contrôleur présenté sur la figure 4. La figure 6 illustre schématiquement un contrôleur avec modèle interne M utilisé alternativement dans la méthode selon l'invention. Prenons par exemple le contrôleur le plus courant : le contrôleur PI. Ce contrôleur 20 possède 2 paramètres. Dans une approche classique, le choix de ces paramètres se fait de la manière suivante. Le moteur est placé dans une condition de fonctionnement donnée. Les 2 paramètres sont réglés (optimisés) "à la main" (ou de manière automatique) en essayant plusieurs combinaisons de réglage. Une combinaison de réglage est arrêtée. Le moteur est 25 placé dans une autre condition fonctionnement. Une nouvelle combinaison de réglage est arrêtée. Et ainsi de suite pour chacune des conditions de fonctionnement du moteur. Cette approche soulève de nombreux problèmes. D'abord elle nécessite beaucoup de temps pour effectuer cette calibration. Ensuite, la cohérence des résultats entre les différents réglages n'est pas forcément assurée. Enfin, cette approche ne donne les valeurs des réglages que 30 lorsque le moteur est sur un point stationnaire : quelle valeur adopter lorsque le moteur passe d'une condition opératoire à une autre ? Selon l'invention, on y remédie en exprimant la valeur des paramètres du contrôleur à l'aide d'un nombre très limité de paramètres du modèle physique. A partir de ces paramètres, on calcule automatiquement par des formules (données en annexe 2) les valeurs zc et 5c correspondant à la constante de temps et au retard pur d'un système de premier ordre retardé. Connaissant par ailleurs le gain, on en déduit par les formules connues de l'homme du métier les paramètres Kp et Ti du contrôleur PI en fonction des conditions opératoires. L'invention concerne également un moteur à combustion interne comprenant un système de traitement de gaz d'échappement. Le moteur (10), en particulier un moteur Diesel, débite des gaz d'échappement dans une ligne d'échappement (20). Le système comprend un catalyseur d'oxydation (15) et un filtre à particules (16) placés dans la ligne d'échappement. Le catalyseur d'oxydation est placé en amont du filtre à particules. Des éléments tels qu'une turbine (30) ou une prise de recirculation des gaz d'échappement (40) peuvent être présents entre le moteur et le catalyseur d'oxydation. Le système inclut un système d'injection (80) agissant dans le moteur et piloté par un système de contrôle (50). Alternativement, le système de contrôle peut agir sur un système d'injection annexe (85) injectant directement dans la ligne d'échappement entre le moteur et le catalyseur d'oxydation. Le système de contrôle selon l'invention est adapté à la mise en oeuvre des procédés selon l'invention, en ce qu'il agit sur le système d'injection pour fournir un débit d'hydrocarbures souhaité en amont du catalyseur d'oxydation, selon les lois de commande de l'invention. Le système inclut préférentiellement un capteur de température (17) en amont du catalyseur d'oxydation qui fournit une valeur de température Tl représentative de la température des gaz en amont du catalyseur d'oxydation. Alternativement, la grandeur Ti peut être construite en utilisant l'information d'autres capteurs disponibles dans le moteur et la ligne d'échappement. Le système inclut préférentiellement un capteur de température (18) en aval du catalyseur d'oxydation qui fournit une valeur de température T2 représentative de la température des gaz en amont du filtre à particules. Alternativement, la grandeur T2 peut être construite en utilisant l'information d'autres capteurs disponibles. Cette grandeur T2 est la grandeur que vise à rapprocher le système de contrôle d'une température désirée Tc par la méthode de la présente invention. Le système de contrôle peut prendre en compte les conditions ambiantes (90) et les conditions opératoires (11) du moteur. Annexe 1 : Gains, valeurs k1, k2, Lu et r Introduction des qains G, Ginst, et Kd Le gain stationnaire G est défini comme le rapport entre la valeur stationnaire d'une variation de T2 et la valeur stationnaire de la variation correspondante de débit massique d'hydrocarbures u ayant induit cette variation de T2, tous les autres paramètres (notamment T1 et Q) restant constants par ailleurs. G AT2 r AQHC1 A exprime une variation par rapport à une valeur de référence quelconque. Le gain stationnaire G peut être calculé en prenant en compte notamment l'enthalpie de réaction des hydrocarbures oxydés dans le catalyseur d'oxydation, le débit de gaz entrant dans le catalyseur d'oxydation, la température Tl, la capacité calorifique de ce gaz, le rendement r, les conditions ambiantes EXT et les conditions opératoires du moteur ENG. En particulier, le gain stationnaire G peut être estimé par la formule suivante : Ah G= Q Cpg' où Ah est l'enthalpie de la réaction d'oxydation des hydrocarbures oxydés par unité de masse d'hydrocarbures oxydés et s'exprime préférentiellement en kJ/kg. La valeur du gain instantané Ginst est la valeur du gain à appliquer à l'instant courant. Cette valeur peut être prise égale à la valeur du gain stationnaire G. Alternativement, elle peut être prise égale à la valeur gain stationnaire G corrigée par une fonction de la dérivée filtrée du débit de gaz traversant le catalyseur. En particulier, Ginst est défini de la manière suivante : G 1 f(dQsat Q dt où Qsat = min(max(Q,lmnQ),lmxQ), lmnQ étant une valeur positive quelconque 25 préférentiellement égale à 0 et lmxQ étant une valeur positive quelconque strictement Ginst = 1ùKB supérieure à lmnQ ; KB étant une valeur quelconque ; f étant un filtre avec une fonction passe-bas ; f étant préférentiellement un filtre de type premier ordre dont la constante de temps vaut préférentiellement 1/k2.. Introduction des valeurs kl et k2 Les valeurs kl et k2 sont des valeurs caractéristiques du catalyseur d'oxydation. Elles peuvent être identifiées à l'aide de la formule EQREP en mesurant T2 lors d'une variation échelon de T1 L" T 2(t) =H t-- exp V) kl L1 x EQREP : i t-L / v J exp( k2 z) o où H est la fonction échelon définie par H(t) kl k2 L kl k2 L I12 z dz vz- v 1sit0 0 sit Les valeurs kl et k2 peuvent être identifiées à toute forme équivalente de cette formule. Alternativement, elles peuvent être calculées à l'aide des formules suivantes I1(x) = x2 Ek=o k! F(2 + k) ; F est la fonction gamma. 4 hg Ga Cpg Cpg hg Ga (1- s) ,mon Cpm Ga est le rapport géométrique surface sur volume représentatif de la surface d'échange entre le gaz et le solide. hg coefficient d'échange convectif d'enthalpie entre le gaz et solide 8 est la fraction de vide du pain catalytique ,og est la masse volumique moyenne du gaz pin est la masse volumique moyenne du solide Cpg est la capacité calorique moyenne du gaz Cps est la capacité calorique moyenne du solide Pour les définitions détaillées des valeurs, on peut se reporter à la publication suivante : C. Depcik, D. Assanis, One-dimensional automotive catalyst modeling, Progress ln Energy Combustion and Science., vol. 31, 2005, pp 308-369. Introduction de Lu En se basant sur des résultats expérimentaux, on peut établir que les réponses à des variations de T1 et à des variations de débit d'hydrocarbures sont différentes (FIG. 7). En adaptant la longueur du catalyseur d'oxydation dans l'expression de la réponse indicielle normalisée à une variation de T1 , on obtient une réponse très similaire à la réponse indicielle normalisée à une variation de débit d'hydrocarbures (FIG. 8), tous les autres paramètres restants égaux par ailleurs. On définit la longueur Lu comme la longueur de catalyseur d'oxydation permettant de rapprocher au mieux d'une part la réponse à une variation de T I sur un catalyseur de longueur Lu , et d'autre part la réponse indicielle à une variation de débit d'hydrocarbures sur un catalyseur d'oxydation de longueur L 24 kl = k2 = On note que la notion de longueur différente entre L et Lu est équivalente à la notion de temps de réponse différent entre la réponse à une variation de Ti et la réponse à une variation de débit d'hydrocarbures. Les variations de Ti sont considérées comme des perturbations. La stratégie de précompensation consiste à synchroniser les effets des variations de Ti sur T2 avec les effets des variations de débit d'hydrocarbures u sur T2 afin de rejeter les effets des variations de Ti (FIG. 9). Le calcul du retard de synchronisation est réalisé dans toutes les conditions opératoires. Introduction du rendement r De manière classique, le rendement r est défini comme le rapport entre le débit d'hydrocarbures consommés dans le catalyseur d'oxydation et le débit total d'hydrocarbures u en amont du catalyseur d'oxydation. u - QHC2 catalyseur d'oxydation ; la valeur QHC2 est le débit d'hydrocarbures présents en aval du catalyseur d'oxydation. II est communément constaté que, dans la gamme de température nécessaire pour effectuer une régénération, le rendement dépend essentiellement du débit de gaz traversant le catalyseur d'oxydation et est décroissant avec le débit Q , tel qu'illustré dans la FIG. 10. r = ; la valeur u est le débit d'hydrocarbures présents en amont du Annexe 2 : Calcul de TFB 1. Contrôleur de type proportionnel intégral (PI) Une première méthode consiste à calculer TFB à l'aide d'un contrôleur PI défini de la façon suivante : Entrée Tc-T2 Sortie V Fonction de transfert : Kp 1+ 1 Zz s s étant la variable de Laplace ; Kp et Ti étant les paramètres du contrôleur PI. On calcule alors TFB comme la valeur V limitée selon l'équation : TFB = max(min(V, lim max), lim min Les valeurs lim_max et lim min sont définies de manière à prendre en compte les contraintes physiques du système. Un algorithme d'anti-emballement est préférentiellement utilisé en supplément. En particulier, les paramètres Kp et ZZ du contrôleur PI peuvent être adaptés avec le débit de gaz, en se basant sur la constante de temps TC , le retard CSC et le gain Ginst d'un système de premier ordre retardé, à l'aide des réglages automatiques de Tavakoli-Fleming ou de Ziegler-Nichols ou de toutes autres formules de réglage automatique de contrôleurs PI basées sur un système de premier ordre retardé. L'homme du métier sait que tous ces réglages expriment Kp et ii en fonction des paramètres TC , & et Ginst En particulier la valeur TC est préférentiellement exprimée par aLu Tc exp( 2 kl Lu /v) k2 I1(2 kl Lu /v) ou alternativement par zc = 2 aLu~ kl Lu k2 v ou alternativement par = aLu exp(ù 2 k2 (5pist) k2 I1(2 k2 5pist) ou alternativement par rc = 2 -aLu-Fr./k2 Spist k2 En particulier la valeur Sc est préférentiellement exprimée par Lu 1 kl Lu (5c= +ù ùDC v k2 v ou alternativement par 6c = Spist ù zc ou alternativement par Lu Sc= +6pistùw v Dans les équations précédentes de & et Tc, la valeur t t est obtenue en résolvant l'équation implicite suivante t Lu = Jv(w)dw tùSpist(t) Dans les équations précédentes de Sc et DC , la valeur aLu peut être une valeur quelconque et en particulier une valeur proche de 0,45. Alternativement la valeur aLu peut être une valeur évaluée à l'aide de la formule suivante 27 28 kl Lu aLu =1ù exp k2 v 1 + J exp(ùk2 'r).\ o klk2Lu klk2Lu I12 dz vZ v v x 2. Contrôleur de type intégral simple Alternativement, une seconde méthode peut être utilisée pour calculer TFB . Le terme TFB se calcule à l'aide un contrôleur de type intégrateur simple défini ainsi : Entrée Tc ù T2 Sortie V - Fonction de transfert : 1 ZI s S étant la variable de Laplace ; ZI étant le paramètre du contrôleur. On calcule TFB comme la valeur V limitée selon l'équation : TFB = max(min(V, lim max), lim min Les valeurs limmax et lim min sont définies de manière à prendre en compte les contraintes physiques du système. Un algorithme d'anti-emballement est préférentiellement utilisé en supplément. La valeur ZI s'exprime comme une fonction de la valeur Spist définie précédemment. En particulier la valeur ZI peut être calculée selon la formule ZI =2 5pist. Alternativement, ZI peut être calculée comme une fonction des valeurs Zc et & définies précédemment. En particulier, la valeur TI peut être calculée selon la formule ZI=~c+3 zc. 3. Contrôleur de type modèle interne du premier ordre Alternativement, et en complément des différentes configurations présentées pour le calcul de uFF , le calcul de TFB peut se faire de manière classique avec un modèle interne, de préférence selon la configuration présentée dans la FIG. 4, et en calculant la sortie du modèle interne T2 int de la manière suivante rc(t) .zc(t) = ùxc(t) + Ginst(t) uFB(t) T2int(t) = xc(t ù sc(t)) avec les valeurs W , & et Ginst définies précédemment. On définit .x comme la dérivée de la variable x par rapport au temps t. Le terme TFB se calcule à l'aide de la sortie d'un système dynamique CT dont la sortie est V et dont l'entrée est la différence entre la consigne Tc et TA ; TA étant la différence filtrée entre la température T2 et la sortie T2 int d'un modèle interne d'entrée uFB . Le terme TFB correspond à la sortie limitée de CT selon l'équation TFB = max(min(V, lim max), limmin) les valeurs 11m max et llm_min étant définies de manière à prendre en compte les contraintes physiques du système. Le système dynamique CT est choisi de manière classique dans la synthèse d'un contrôleur à modèle interne de type premier ordre retardé. En particulier, il assure le rejet de perturbations additives constantes en sortie et le suivi de la consigne Tc. Le système dynamique CT est une fonction des conditions opératoires ENG ou des conditions ambiantes EXT . En particulier CT peut être choisi comme un scalaire de valeur 1. Le bloc C de la figure 4 est composé du bloc CT précédemment défini en série avec l'inverse du gain Ginst. 4. Contrôleur de type prédicteur de Smith De manière équivalente à la commande par modèle interne avec un modèle interne de type premier ordre retardé décrite précédemment, le calcul de TFB peut se faire de manière classique avec un contrôleur de type prédicteur de Smith, de préférence selon la configuration présentée dans la FIG. 5, et en calculant la sortie d'un modèle de premier ordre retardé utilisant ZC , 8C et Ginst pour paramètres. La sortie du modèle de premier ordre non retardé (bloc N) est donnée par T 2N int = xc avec rc(t) .xc(t) = ùxc(t) + Ginst(t) uFB(t) . La sortie du modèle de premier ordre retardé DT2Nint correspond à la valeur T2Nint retardée par l'opérateur retard DN d'un retard de valeur CYC(t) . Elle est donnée par DT2Nint(t) = T2Nint(t ù &(t)) . Un filtre de robustesse F est préférentiellement utilisé dans la structure du prédicteur de Smith. En particulier, ce filtre peut être un filtre de type premier ordre. Le terme TFB se calcule préférentiellement à l'aide un contrôleur PI d'entrée Tc ù T 2N intù F(T 2 ù DT 2NINT) et de sortie V ; la fonction de transfert du contrôleur PI étant Kp ; S étant la variable de Laplace ; Kp et Ti étant les paramètres du contrôleur PI ; Le paramètre Kp vaut de préférence 1 et le paramètre zi vaut de préférence W défini précédemment. On calcule TFB comme la valeur V limitée selon l'équation TFB = max(min(V, lim_ max), lim_ min) Les valeurs 11m_max et hm min sont définies de manière à prendre en compte 20 les contraintes physiques du système. Un algorithme d'anti-emballement est préférentiellement utilisé en supplément. Le bloc CS de la figure 5 correspond au terme TFB calculé comme décrit précédemment et transformé en débit d'hydrocarbures uFB à l'aide du gain Ginst. 25 5. Contrôleur de type Modèle interne du type Erf u En complément des différentes configurations présentées pour le calcul de uFF , le calcul de TFB peut se faire de manière classique avec un modèle interne, de préférence selon la configuration présentée dans la FIG. 4, et en calculant la sortie du modèle interne T2 int à l'instant t à l'aide de la valeur de sortie T2int à l'instant tN ; ti étant le i-ème échantillon de la valeur du temps t courant depuis un temps initial to ; i étant un entier naturel positif ; l'ordre des échantillons respectent l'écoulement du temps t de telle sorte que pour N > i > 0 on a tN > ti > t._1 > to ; tN étant l'échantillon de t le plus récent. La valeur de sortie T 21nt(tN) à l'instant tN se calcule pour N 1 suivant N T2 int(tN) = Ginst(to) uFB(to) + AnTn (tN ) n=1 An = Ginst(tN_n+1) uFB(tN_n+l) ù Ginst(tN_n) uFB(tN_n ) Tn (tN) = 1 1ùerf n H(tN ù tn ) 2\ 2 zn ~n = Luù tN J!v(s)ds kl tN-n zn = k2 2 (s )ds k12 2 _ 2 où erf est la fonction d'erreur définie par erf (x) = Je'°dw E n 0 particulier, on utilise préférentiellement T2int(t) = T2int(tN) avec t tN La sortie du modèle interne T2int peut aussi être calculée par toute approximation des formules ci-dessus. En particulier, on peut utiliser l'approximation qui consiste à tronquer la somme à un nombre limité D de termes à l'aide de la formule D T2 int(t) = Ginst(tN_D) uFB(tN_D) + E AnTn (tN ) n=1 tN-n 10 15 6. Contrôleur de type Modèle interne du premier ordre u + T1 Alternativement, et en complément des différentes configurations présentées pour le calcul de uFF , le calcul de TFB peut se faire de manière classique avec un modèle 5 interne de type deux entrées et une sortie, de préférence selon la configuration présentée dans la FIG. 6, et en calculant la sortie du modèle interne T2int de la manière suivante. T2int(t) = xd(t ù &1(t))+ xe(t ù &(t)), où les valeurs CSC et ZC sont décrites précédemment, la valeur xe est décrite par rc(t) ze(t) = ùxe(t) + Ginst(t) u(t), et la valeur xd est décrite par rd (t) .d (t) = -xd(t) + Kd (t) (T 1(t) - Thar) avec les valeurs 8d , rd , Kd et Spist définies de la manière suivante : La valeur rd est préférentiellement donnée par 'rd - exp( 2 kl L/v) k2 I1(2 kl L/v) ou alternativement par 2 kl L rd = ùaL 7r k2 v ou alternativement par aL exp(- 2 k2 6pistL) k2 I1(2 k2 SpistL) ou alternativement par 'rd = ? aL.. Jk2 SpistL k2 aL rd= 20 La valeur &d est préférentiellement donnée par &1=L ù+ 1 kl Lùrd v k2 v ou alternativement par &d =ôpistLùrd ou alternativement par (d = L +8pistLùrd v La valeur 8pistL est préférentiellement donnée par résolution de l'équation implicite suivante t L = Jv(w)dw tù8pistL(t) La valeur aL peut être une valeur quelconque et en particulier une valeur proche de 0,45. Alternativement la valeur aL peut être une valeur évaluée à l'aide de la formule suivante k1 Ll aL=1ùexp x 33 kl L k2 v 1+ J exp(k2 o k1 k2 L k1 k2 L I1 2 z dr vz v Le terme TFB se calcule à l'aide de la sortie d'un système dynamique CT défini précédemment. En particulier CT peut être choisi comme un scalaire de valeur 1. Le bloc C de la figure 6 est composé du bloc CT précédemment défini en série avec l'inverse du gain Ginst. 7. Contrôleur de type Modèle interne de type type Erf_u + Erf_T 10 15 Alternativement, et en complément des différentes configurations présentées pour le calcul de uFF , le calcul de TFB peut se faire de manière classique avec un modèle interne de type deux entrées et une sortie, de préférence selon la configuration présentée dans la FIG. 6, et en calculant la sortie du modèle interne T2int à l'aide des sorties xq et xp de deux modèles intermédiaires ; dans ces modèles, ti étant le i-ème échantillon de la valeur du temps t courant depuis un temps initial to ; i étant un entier naturel positif ; l'ordre des échantillons respectent l'écoulement du temps t de telle sorte que pour N > i > 0 on a tN > ti > t._1 > to ; tN étant l'échantillon de t le plus récent. La valeur de sortie xq à l'instant tN se calcule pour N 1 suivant N xq(tN) = Ginst(t0) u(t0) + Aqn Tn n=1 Aqn = Ginst(tN_n+1) u(tN_n+l) ù Ginst(tNùn) u(tN_n ) ( r \\ Tn = 1 1ù erf H(tN ù tn ) 2 2 Zn tN ;n = Lu ù Jv(s)ds k2 kl tN_n zn tNf k2 V2(S)ds J k1 tNùn et la valeur de la sortie xp à l'instant tN se calcule pour N 1 suivant N xp(tN) = Kd (to) (T 1(to) ù Tbar) + Apn Tn n=1 Apn = Kd (tN_n+l) (T l(tN_n+l) ù Tbar) ù Kd (tN_n)(T 1(tN_n) ù Tbar) / / \\ Tn = 1 1ù erf n H(tN ù tn ) 2 2 zn "n =Lù J---v(s)ds k2 tN-n kl t" k2 v2 kl(s)ds 2 x où erf est la fonction d'erreur définie par erf (x) _ je2dw En 0 particulier, on utilise préférentiellement T2int(t) = xq(tN) + xp(tN) avec t tN . La sortie du modèle interne T2int peut aussi être calculée par toute approximation des formules ci-dessus. En particulier, on peut utiliser l'approximation qui consiste à tronquer la somme à un nombre limité D de termes à l'aide des formules D 35 tN tN-n xq(tN) = Ginst(tN_D ) u(tN-D) + Aqn Tn (tN ) n=1 D xp(tN) = Kd (tN-D) (T 1(tN-D) ù Tbar) + E Apn Tn (tN ) n=1 Le terme TFB se calcule à l'aide de la sortie d'un système dynamique CT défini précédemment. En particulier CT peut être choisi comme un scalaire de valeur 1. Le bloc C de la figure 6 est composé du bloc CT précédemment défini en série avec l'inverse du gain Ginst.15 REVENDICATIONS1. Procédé de traitement de gaz d'échappement d'un moteur à combustion interne équipé d'un filtre à particules et d'un catalyseur d'oxydation placé en amont dudit filtre, dans lequel on contrôle une régénération active dudit filtre au moyen d'un contrôle dynamique d'une température T2 dudit gaz sortant dudit catalyseur d'oxydation, caractérisée en ce que: i- on détermine une température de consigne Tc dudit gaz sortant dudit catalyseur pour assurer une régénération active dudit filtre ; ii- on détermine une demande de température nécessaire pour contrôler ladite température T2 ; iii- on détermine un gain stationnaire G, défini comme le rapport entre la valeur stationnaire d'une variation de T2 et une valeur stationnaire d'une variation de débit massique d'hydrocarbures ayant induit ladite variation de T2 ; iv- on corrige ledit gain G de façon à prendre en compte une variation de débit de gaz en entrée dudit catalyseur, le gain étant transitoirement augmenté par rapport à G pour une augmentation de débit de gaz et transitoirement diminué par rapport à G pour une diminution de débit de gaz ; v- on définit une loi de commande d'un débit d'hydrocarbures en entrée dudit catalyseur d'oxydation, par le rapport entre ladite demande de température et ledit gain corrigé ; vi- on modifie le débit d'hydrocarbures en appliquant la loi de commande pour que la température dudit gaz sortant dudit catalyseur respecte ladite température de consigne. 2. Procédé selon la revendication 1, dans lequel on corrige ledit gain G en calculant une dérivée du débit de gaz traversant ledit catalyseur. 3. Procédé selon la revendication 1, dans lequel on corrige ledit gain G en appliquant un retard au gain selon la formule suivante : Ginst(t) = G(t ù D(t)) avec Ginst le gain corrigé, D le retard déterminé en fonction du débit de gaz Q et de ses variations, t le temps. 4. Procédé de traitement de gaz d'échappement d'un moteur à combustion interne équipé d'un filtre à particules et d'un catalyseur d'oxydation placé en amont dudit filtre, dans lequel on contrôle une régénération active dudit filtre au moyen d'un contrôle dynamique d'une température T2 dudit gaz sortant dudit catalyseur d'oxydation, caractérisée en ce que : i- on détermine une demande de température TFF nécessaire pour compenser un effet d'une variation de température Tl dudit gaz en entrée dudit catalyseur, TFF étant une fonction d'une différence entre T1(t û 8(t)) et une valeur de référence Tbar, où T1(t û 8(t)) est une expression d'une mesure de Tl (t) retardée d'un retard de valeur 8(t) ; ii- on détermine un gain stationnaire G, défini comme le rapport entre la valeur stationnaire d'une variation de T2 et une valeur stationnaire d'une variation de débit massique d'hydrocarbures ayant induit ladite variation de T2, tous les autres paramètres restants constants par ailleurs ; iii- on définit une loi de commande uFF d'un débit d'hydrocarbures en entrée dudit catalyseur d'oxydation, par le rapport entre ladite demande de température TFF et ledit gain ; iv- on modifie le débit d'hydrocarbures en appliquant la loi de commande pour que la température dudit gaz sortant dudit catalyseur respecte ladite demande de température TFF. 5. Procédé selon la revendication 4, dans lequel le retard 8 est défini en fonction d'un retard 82 obtenue en résolvant l'équation implicite suivante : L ù Lu = Jv(w)dw t-82(t) avec L une longueur du catalyseur d'oxydation. Lu : une longueur de catalyseur d'oxydation permettant de rapprocher au mieux d'une part une réponse indicielle normalisée à une variation de TI sur un catalyseur de longueur Lu, et d'autre part une réponse indicielle normalisée à une variation de débit d'hydrocarbures sur un catalyseur d'oxydation de longueur L.v : une vitesse moyenne du gaz dans ledit catalyseur d'oxydation. 6. Procédé selon l'une des revendications 4 et 5, dans lequel on corrige ledit gain G de façon à prendre en compte une variation de débit de gaz en entrée dudit catalyseur, le gain étant transitoirement augmenté par rapport à G pour une augmentation de débit de gaz et transitoirement diminué par rapport à G pour une diminution de débit de gaz. 7. Procédé de traitement de gaz d'échappement d'un moteur à combustion interne équipé d'un filtre à particules et d'un catalyseur d'oxydation placé en amont dudit filtre, dans lequel on contrôle une régénération active dudit filtre au moyen d'un contrôle dynamique d'une température T2 dudit gaz sortant dudit catalyseur d'oxydation, caractérisée en ce que: i- on détermine une température de consigne Tc dudit gaz sortant dudit catalyseur pour assurer une régénération active dudit filtre ; ii- on détermine une demande de température TFB nécessaire pour que la température T2 dudit gaz sortant dudit catalyseur respecte ladite température de consigne Tc, ladite demande TFB étant calculée à partir d'un contrôleur dont les paramètres sont déterminés au moyen d'un modèle physique du catalyseur d'oxydation, basé sur un bilan d'énergie d'une phase gaz et un bilan d'énergie d'une phase solide, en fonction des conditions opératoires du moteur ; iii- on détermine un gain stationnaire G, défini comme le rapport entre la valeur stationnaire d'une variation de T2 et une valeur stationnaire d'une variation de débit massique d'hydrocarbures ayant induit ladite variation de T2; iv- on définit une loi de commande uFB d'un débit d'hydrocarbures en entrée dudit catalyseur d'oxydation, par le rapport entre ladite demande de température TFF et ledit gain ; v- on modifie le débit d'hydrocarbures en appliquant la loi de commande pour que la température dudit gaz sortant dudit catalyseur respecte ladite température de consigne. 8. Procédé selon la revendication 7, dans lequel on corrige ledit gain G de façon à prendre en compte une variation de débit de gaz en entrée dudit catalyseur, le gain étant transitoirement augmenté par rapport à G pour une augmentation de débit de gaz et transitoirement diminué par rapport à G pour une diminution de débit de gaz. 9. Procédé selon les revendications 4 et 7, dans lequel ladite loi de commande appliquée pour modifier le débit d'hydrocarbures est définie par la somme desdites lois de commande uFF et uFB. 10. Procédé selon la revendication 9, dans lequel on corrige ledit gain G de façon à prendre en compte une variation de débit de gaz en entrée dudit catalyseur, le gain étant transitoirement augmenté par rapport à G pour une augmentation de débit de gaz et transitoirement diminué par rapport à G pour une diminution de débit de gaz. 11. Moteur à combustion interne comprenant un système de traitement de gaz d'échappement, un catalyseur d'oxydation (15) placé en amont d'un filtre à particules (16), et un système d'injection (80) agissant dans le moteur, caractérisé en ce qu'il comporte un système de contrôle (50) pour piloter ledit système d'injection de façon à fournir un débit d'hydrocarbures en amont du catalyseur d'oxydation, en accord avec un des procédés selon l'une des revendications 1 à 9.