Brevet FR3078745A1 — publié le 20190913

PROCEDE DE COMMANDE D'UN MOTEUR THERMIQUE [0001] La présente invention concerne de manière générale la gestion d'un moteur thermique monté sur un véhicule automobile thermique ou hybride. En particulier, l'invention concerne la gestion du moteur thermique en fonction d'une efficacité de traitement d'effluents présents dans les gaz 5 d'échappement. [0002] Il est connu dans l'art antérieur que des dispositifs de traitement des effluents de gaz d'échappement ont une efficacité qui dépend de la température d'un catalyseur placé dans une ligne d'échappement. Il est donc important de connaître la température du catalyseur (autrement appelé 10 pain catalytique) pour correctement estimer l'efficacité de traitement des gaz, et pour éventuellement décider de lancer un pilotage spécifique du moteur thermique pour augmenter la température des gaz d'échappement et donc du catalyseur. [0003] Afin de connaître la température du catalyseur, il est 15 possible soit d'en mesurer directement la température interne, mais cela impose de fortement complexifier la structure de la ligne d'échappement pour implanter une sonde de température au sein même du catalyseur. On peut alternativement estimer (c'est-à-dire calculer, simuler, ou déterminer par un modèle théorique) la température du catalyseur, comme divulgué dans le 20 document EP1108124 qui propose de calculer une température du catalyseur, mais ces méthodes sont peu précises, notamment car une seule valeur de température globale est estimée, ce qui peut poser problème notamment en cycle transitoire, lorsque le catalyseur présente des gradients internes de température. 25 [0004] Un but de la présente invention est de répondre aux inconvénients de l'art antérieur mentionnés ci-dessus et en particulier, tout d'abord, de proposer un procédé de commande de moteur thermique qui présente un pilotage précis d'un moteur thermique pour optimiser l'efficacité de traitement des effluents gazeux des gaz d'échappement. 2 [0005] Pour cela un premier aspect de l'invention concerne un procédé de commande d'un moteur thermique de véhicule équipé : - d'un moteur thermique, - d'au moins un catalyseur pour traiter des gaz d'échappement du moteur 5 thermique, - d'au moins une sonde aval de température de gaz d'échappement agencée en aval du catalyseur, le procédé de commande comprenant les étapes consistant à : - estimer, à un instant t du fonctionnement du moteur thermique, une 10 température de gaz d'échappement en une pluralité de points le long du catalyseur, en prenant en compte au moins : - une valeur de température initiale du catalyseur, - une valeur de température de gaz d'échappement en amont du catalyseur à l'instant t et/ou à un instant précédent t-dt, et 15 - un débit massique de gaz d'échappement, à l'instant précédent t-dt uniquement, - estimer à l'instant t une température de paroi du catalyseur en une pluralité de points le long du catalyseur, en prenant en compte au moins la température de gaz d'échappement le long du catalyseur, 20 - mesurer avec la sonde aval de température à l'instant t, une température de gaz d'échappement en aval du catalyseur, - corriger la température de paroi du catalyseur le long du catalyseur, en prenant en compte la mesure de température de gaz d'échappement en aval du catalyseur, si une erreur d'estimation de la température de gaz 25 d'échappement est détectée, - commander au moins un organe du moteur thermique en fonction de la température, estimée ou corrigée, de paroi du catalyseur le long du catalyseur. [0006] Le procédé selon la mise en oeuvre ci-dessus calcule une température de gaz d'échappement et une température de paroi au sein même 30 du catalyseur en une pluralité de points, ce qui permet de garantir que le pilotage ultérieur tient compte de l'état de température de l'ensemble du catalyseur. De plus, le procédé comporte une étape de mesure physique de la 3 température des gaz d'échappement en aval du catalyseur, pour recaler ou corriger les estimations effectuées, en cas de constatation d'une dérive. Cette étape de vérification permet de garantir que les températures estimées sont correctes, ce qui améliore la fiabilité du procédé de commande. [0007] Il est à noter que le procédé de commande prend en compte le débit massique de gaz d'échappement, c'est-à-dire la chaleur transmise par les gaz d'échappement, mais uniquement à l'instant précédent t-dt, ce qui ne fausse pas la prédiction. Le calcul est donc possible, et tient compte des échanges réels qui se passent au sein du catalyseur. [0008] Avantageusement, la température de gaz d'échappement et la température de paroi au sein du catalyseur en une pluralité de points sont estimées le long du catalyseur selon une seule dimension, en particulier selon la direction d'écoulement des gaz d'échappement. [0009] En d'autres termes, le procédé de commande ne calcule pas une cartographie des températures en deux ou trois dimensions. En effet, la demanderesse s'est aperçue que le calcul le long de la direction axiale du catalyseur procurait une précision et une représentativité suffisantes, sans pour autant nécessiter de grosses ressources de calcul. [0010] Avantageusement, la température de paroi du catalyseur en une pluralité de points le long du catalyseur à l'instant t est estimée en prenant également en compte une température ambiante au véhicule à l'instant t et/ou à l'instant précédent t-dt, et/ou un débit massique d'espèces gazeuses à catalyser à l'instant précédent t-dt, et optionnellement une valeur de richesse. [0011] Il est à noter que le procédé de commande prend en compte le débit d'espèces gazeuses à catalyser, c'est-à-dire que la chaleur dégagée par les réactions exothermiques de catalyse est prise en compte, mais uniquement à l'instant précédent t-dt, ce qui ne fausse toutefois pas la prédiction. Le calcul est donc possible, et tient compte des réactions réelles qui se passent au sein du catalyseur. 4 [0012] Autrement dit, la température de gaz d'échappement Tgau temps t et en une pluralité de points le long du catalyseur (axe x) est calculée avec une première fonction fi qui prend en compte les différentes variables de la manière suivante, pour un point d'abscisse x : 5 7:9(X, = 1.1(7:9(X, t - dt), Tg (X - dx, t - dt), Tg (X dx, t - dt), Ts(x, t - dt), Qech(t - dt)) La température de paroi du catalyseur Ts au temps t et en une pluralité de points le long du catalyseur (axe x) est calculée quant à elle avec une deuxième fonction f2 qui prend en compte les différentes variables de la 10 manière suivante, pour un point d'abscisse x : Ts(x, t) = f2(Ts(x, t - dt), Ts(x - dx, t - dt), Ts(x + dx, t - dt), Tg (X, t - dt),Tamh(t - dt), Qp 01(x, t - dt), Qech(t - dt)) Avec : x: position le long du catalyseur (mm), 15 t: instant (s) Tg: température de gaz d'échappement (K) Ts: température de paroi du catalyseur (K) Tamb température ambiante au système (K) Qech : débit massique en amont du catalyseur (Kg/s), obtenu par mesure ou par 20 estimation avec un modèle issu du contrôle moteur Qpoi : débit massique d'espèces gazeuses à catalyser (Kg/s) (par exemple, monoxyde de carbone (CO), oxydes d'azote (N0x), hydrocarbures imbrûlés (HC)...) Tg(0, t) = Tech(t) : température de gaz d'échappement en amont du catalyseur 25 (K), obtenue par mesure ou par estimation avec un modèle issu du contrôle moteur Tdx, 0) = Tim : température initiale du système (K). Ts(x, 0) = Tim : température initiale du système (K). 5 [0013] Dans le détail, la température de gaz d'échappement Tg en une pluralité de points le long du catalyseur est calculée de la manière suivante : Tg(X,t) = Tg(X,t - dt) + Qech(t dt) (Tg (x, t - dt) - Tg(X - dx, t - dt)) + 5 G2 (Ts(X,t dt) - Tg(X,t - dt)) + G3 (7:9(X dx, t dt) - 2 Tg(X,t dt) + 7:9(X + dx, t dt)) Avec G1, G2, et G3 qui sont des coefficients qui représentent les caractéristiques thermodynamiques du système : dt = pg E A, dx 10 dt G2 =pg E Cp hs,tnt Ageo 9 dt £L G3 = log E Cpg dx2 15 Où : pg : Masse volumique des gaz d'échappement (Kg/mm3) E : Ratio de surface ouverte dans le catalyseur Acs : Aire ouverte du catalyseur (mm2) CP : Capacité thermique massique des gaz d'échappement (J/K/Kg) 20 hs,int : Coefficient de convection gaz/paroi à l'intérieur du catalyseur (W/mm2/K) Ageo : Surface d'échange gaz/paroi à l'intérieur du catalyseur par unité de volume du catalyseur (mm2/mm3) 4 : Conductivité thermique des gaz d'échappement (W/mm/K) 25 6 Dans le détail, la température de paroi du catalyseur Ts en une pluralité de points le long du catalyseur est calculée de la manière suivante : Ts(x, t) = Ts(x, t - dt) + (Ts(x - dx, t - dt) - 2Ts(x, t - dt) + Ts(x + dx, t dt)) + S2 (Ts(X, t - dt) - Tg (X, t - dt)) + Ei [co,Fic,Nxl = Qpoi,i(x - dx, t - 5 dt) th(Ts(x, t dt), 'h(t - dt), Ram't) + S4 (Ts(X, t - dt)- - dt)) amb (t Avec : Si, S2, S3, et S4 qui sont des coefficients qui représentent les caractéristiques thermodynamiques du système : dt 2s = Ps Cs d x2 dt S2 - h tnt A Ps Cs (1 - E)s ' Ageo S3 i = Ps s E) 2 Vau dt Ki lli,max A out hout Ps Cs (1 E) Vcat Où ÀS Conductivité thermique du pain catalytique (W/mm/K) Ps Masse volumique du pain catalytique (Kg/mm3) Cs : Capacité calorifique massique du pain catalytique (J/K/Kg) 20 Ki : Coefficient de production de chaleur dû à la conversion de l'espèce gazeuse à catalyser i (J/Kg) 71 i,max Rendement de conversion maximal de l'espèce gazeuse à catalyser i Vcat : Volume du catalyseur (mm3) rii : Rendement de conversion de l'espèce gazeuse à catalyser i, ( Th E [0,1] ) 25 Q ech,max Débit massique à l'échappement maximum (Kg/s) S4 = dt 10 15 7 A't : Surface extérieure du pain catalytique (mm2) h't : Coefficient de convection entre le pain catalytique et l'air ambiant (W/mm2/K) Ramont valeur de richesse en amont du catalyseur (optionnel), obtenue par 5 mesure ou par estimation avec un modèle de calcul du contrôle moteur. [0014] Avantageusement, le procédé de commande comprend une étape consistant à : - estimer un débit massique local d'espèces gazeuses à catalyser en une 10 pluralité de points le long du catalyseur à l'instant t, en fonction de la température de paroi corrigée. Cette étape consistant à estimer un débit massique local d'espèces gazeuses à catalyser permet de calculer ensuite par intégration l'efficacité globale du traitement catalytique des gaz d'échappement, et optionnellement une valeur de richesse. 15 [0015] Autrement dit, le débit massique local d'espèces gazeuses à catalyser en une pluralité de points le long du catalyseur est calculé avec une troisième fonction f3 qui prend en compte les différentes variables de la manière suivante : Qpol,i(x,t dt) - f3(Qpol,i(x Clt,t - Clt),Ts(X,t dt), ( ech(t dt), Ramont) 20 Dans le détail : t - dt) = Qpoi,t(x - dt, t - dt) (1 - t(Ts(x, t - dt), 0 t ,ech\- - dt), Ramont) c±) Avec : Lcat : Longueur du catalyseur (pain catalytique) (mm) 25 [0016] Avantageusement, l'étape, consistant à estimer le débit massique local d'espèces gazeuses à catalyser en une pluralité de points le 8 long du catalyseur à l'instant t, comprend une étape consistant à estimer le débit d'espèces gazeuses à catalyser en amont du catalyseur, en fonction du débit massique de gaz d'échappement à l'instant précédent t-dt, et/ou d'un régime du moteur thermique à l'instant précédent t-dt et optionnellement une 5 valeur de richesse. [0017] Autrement dit, le débit d'espèces gazeuses à catalyser en amont du catalyseur est calculé avec une quatrième fonction [4 qui prend en compte les différentes variables de la manière suivante : Qpoi(0, t) = f4(Q 'h(t dt), wmot(t - dt), R amont) 10 Dans le détail : t) = ai,o(wmot(t - dt), R amont) + ai,i(wmot(t - dt), Ramont) Qech(t - dt) ai,2(wmot(t dt), Ramont)(Q ech(t dt))2 Avec : ai3O, au, ai,2 : Coefficients du modèle quadratique de la production de l'espèce 15 gazeuse à catalyser i par le moteur thermique w : vitesse de rotation du moteur thermique (rad s-1) [0018] Avantageusement, l'étape consistant à corriger la température de paroi du catalyseur le long du catalyseur comprend les étapes 20 consistant à : - estimer à l'instant tune température de gaz d'échappement en extrémité aval du catalyseur, en fonction : de la température de gaz d'échappement en extrémité aval du catalyseur à l'instant précédent t-dt, 25 de la température de paroi du catalyseur en extrémité aval du catalyseur à l'instant précédent t-dt, de la température ambiante au véhicule à l'instant précédent t-dt, - comparer la température de gaz d'échappement en extrémité aval du catalyseur estimée à l'instant t avec la mesure à l'instant t de la température 9 de gaz d'échappement en aval du catalyseur, - si une différence entre la température de gaz d'échappement en extrémité aval du catalyseur estimée à l'instant t et la mesure à l'instant t de la température de gaz d'échappement en aval du catalyseur excède un seuil 5 prédéterminé (par exemple 10 K) : - corriger la température de gaz d'échappement estimée à l'instant t en une pluralité de points le long du catalyseur, - ré-estimer à l'instant t la température de paroi du catalyseur en une pluralité de points le long du catalyseur, en prenant en compte au 10 moins la température de gaz d'échappement corrigée le long du catalyseur. [0019] Autrement dit, la température de gaz d'échappement en extrémité aval du catalyseur, est estimée avec une cinquième fonction f5 : T aval (t) = f5 (Tg (I, cap t - dt), Taval(t - dt), Tamb(t - dt)) 15 Et : T aval (t) = Taval(t - dt) + Al (T aval (t - dt) - Tg (L cat, t - dt)) + A2 (Tavai(t dt) - Tamb(t - dt)) Avec Al, A2 qui sont des coefficients qui représentent les caractéristiques thermodynamiques du système : 20 Al = dt haval,int Saval,int A2 = raval Vaval raval : Inertie thermique de l'élément entre la sortie du catalyseur et la mesure de température aval (J/K/mm3) raval Vaval dt haval,ext S aval,ext 25 havai : Coefficient d'échange thermique entre le gaz à la sortie du catalyseur et l'élément en aval (W/mm2/K) -10- Saval,int : Surface d'échange thermique entre le gaz à la sortie du catalyseur et l'élément en aval (mm2) haval,ext: Coefficient d'échange thermique entre l'élément en aval du catalyseur et l'air ambiant (W/mm2/K) 5 Saval,ext : Surface d'échange thermique entre l'élément en aval du catalyseur et l'air ambiant (mm2) [0020] Avantageusement, si la température de paroi du catalyseur estimée en l'un des points le long du catalyseur est inférieure à une 10 température seuil, alors au moins un des organes du moteur thermique, tel qu'un dispositif d'injection et/ou un dispositif d'allumage, est commandé pour augmenter ladite température de paroi du catalyseur. [0021] Avantageusement : - un intervalle de temps entre l'instant t et l'instant t-dt est inférieur ou égal à 5 15 ms, et avantageusement inférieur ou égal à 1 ms, et/ou - un intervalle géométrique entre deux points, auxquels la température de paroi du catalyseur est estimée, est inférieur ou égal à 5 mm, et avantageusement inférieur ou égal à 1 mm. [0022] Avantageusement, le véhicule est équipé d'une sonde 20 amont de température de gaz d'échappement agencée en amont du catalyseur, et : - la valeur de température initiale du catalyseur est une température mesurée par la sonde amont de température de gaz d'échappement et/ou la sonde aval de température de gaz d'échappement, avant mise 25 en fonction du moteur thermique, et/ou - la valeur de température de gaz d'échappement en amont du catalyseur à l'instant t et/ou à un instant précédent t-dt est une température mesurée par la sonde amont de température de gaz d'échappement. [0023] Il est alors possible de mettre en place un système d'équations pour un modèle thermique du catalyseur discrétisé en N noeuds selon la direction longitudinale x du catalyseur, comme suit : Ô(t) = g (0 (t - dt), u(t), t) ty(t) = [0 0 ... 0 0 0 ... 0 1] 0 5 Avec : OT = [Tgl, Tg2, ..., TgN, Ts1, n2, ... TSN, Tavai] : états du système uT = [Q 'h, T'h, Tamb] : entrées du système YT = Tavai : sortie du système 10 [0024] Pour obtenir une représentation d'état du système, il vient é(t) = A(0 (t - dt), u(t)) 0(t - dt) + B (0 (t - dt), u(t)) u (t) + W y (t) = C 0 (t) + V Avec 15 A : matrice jacobienne détaillant la variation de la fonction g par rapport aux états 0 sur un point de fonctionnement (0(t - dt), u(t)) donné : agi - agi- a01 aoN A(0 (t - dt), u(t)) = agN - a,gN a 0 1 a ° N E kt- dt),u(t) B : matrice jacobienne détaillant la variation de la fonction g par rapport aux entrées u sur un point de fonctionnement 0(t - dt), u(t) donné : agi agi - a ui a u3 a.gN a.gN a u 1 a U3 - E kt- dt),u(t) 20 B (0 (t - dt), u(t)) = -12- W : variance du modèle, qui peut être déterminée à l'aide d'un essai de validation du modèle ; C = [0 0 ... 0 0 0 ... 0 1] : vecteur liant la sortie aux états ; V : variance de la mesure, qui peut être obtenue dans les données 5 caractéristiques du (des) capteurs de mesure. Il est à noter que les valeurs des coefficients ci-dessus peuvent être déterminées par des essais sur banc. [0025] Une fois la représentation d'état obtenue, on peut estimer les états présents dans la représentation d'état en fusionnant Taval du modèle 10 et de la mesure. L'estimateur se compose de deux étapes, une étape de prédiction, et une étape d'estimation : Pour l'étape de prédiction : 6pred(t) = est(t - dt) + A Cd est(t - dt), u(t)) est(t - dt) dt + B est(t - dt), u(t)) u(t) dt 15 Ou bien : Cipred(t) = est(t - dt) + g C est(t - dt), U(t), t) dt Avec 6pred: états de l'estimateur obtenus durant la phase de prédiction C'est : états de l'estimateur obtenus durant la phase d'estimation 20 Pour l'étape d'estimation : est(t) = prece) K (t)(Y (t) - C 6precl(t)) Avec K : vecteur des gains déterminant le niveau de confiance accordé aux modèles 25 par rapport à la mesure ; -13- Y : sortie mesurée (ici Tay,/ mesurée par la première sonde de température (typiquement un thermocouple)). [0026] Avantageusement, les étapes consistant à : 5 - estimer la température de gaz d'échappement en une pluralité de points le long du catalyseur, - estimer la température de paroi du catalyseur en une pluralité de points le long du catalyseur, - mesurer la température de gaz d'échappement en aval du catalyseur, 10 - corriger la température de paroi du catalyseur le long du catalyseur, sont effectuées pendant une phase de conduite du véhicule, et idéalement par un calculateur embarqué du véhicule. [0027] Avantageusement, le catalyseur est un catalyseur trois voies, qui traite simultanément monoxyde de carbone (CO), oxydes d'azote 15 (N0x), et hydrocarbures imbrûlés (HC). Le modèle de calcul selon l'invention est bien adapté à ce type de catalyseur, car il prend en compte les réactions de catalyse qui peuvent impacter la température du catalyseur, en raison de leur caractère exothermique. [0028] Un second aspect de l'invention concerne un moteur 20 thermique de véhicule automobile, comprenant au moins un organe commandé par le procédé de commande selon le premier aspect. [0029] Un troisième aspect de l'invention concerne un véhicule automobile thermique ou hybride, comprenant au moins un moteur thermique selon le deuxième aspect. 25 [0030] D'autres caractéristiques et avantages de la présente invention apparaîtront plus clairement à la lecture de la description détaillée qui suit d'un mode de réalisation de l'invention donné à titre d'exemple nullement limitatif et illustré par les dessins annexés, dans lesquels : -14- - la figure 1 représente une vue simplifiée d'un moteur thermique avec une ligne d'échappement comprenant un catalyseur dont la température est estimée avec le procédé de commande selon l'invention ; - la figure 2 représente un schéma du catalyseur représenté figure 5 1 ; - la figure 3 représente un diagramme de fonctionnement d'un superviseur de commande du véhicule équipé du moteur thermique de la figure 1. [0031] La figure 1 représente un moteur thermique 10 équipé d'un 10 dispositif de catalyse 20. Typiquement, le moteur thermique 10 comprend un ensemble piston/cylindre 11, une commande d'admission 12, un débitmètre d'air 13, un dispositif d'injection 14, et un dispositif d'allumage 15 pour commander une combustion d'un mélange air/carburant. [0032] Le dispositif de catalyse comprend quant à lui un 15 catalyseur 21, une sonde aval de température 23, et en option une sonde amont de température 22. On peut également prévoir une sonde lambda pour estimer une richesse. Avantageusement, le catalyseur 21 est un catalyseur trois voies, qui provoque trois réactions simultanées : Une réduction des oxydes d'azote en diazote et en dioxyde de carbone : 20 2NO + 2C0 -> N2 + 2002 Une oxydation des monoxydes de carbone en dioxyde de carbone : 2C0 + 02 -> 2CO2 Une oxydation des hydrocarbures imbrulés (HC) en dioxyde de carbone et en eau : 25 4CxHy + (4x+y)02 -> 4xCO2 + 2yH2O [0033] Il est à noter que l'exemple ici est celui d'un moteur à allumage commandé (moteur essence), mais on peut mettre en oeuvre l'invention pour un moteur à allumage par compression, ou moteur Diesel. -15- [0034] La figure 2 représente un schéma du catalyseur 21, pour montrer la structure interne qui comprend typiquement un pain catalytique en structure nid d'abeille par exemple, et qui est traversé par un flux de gaz d'échappement, caractérisé par un débit massique de gaz d'échappement que 5 l'on peut noter Û -ech, et qui contient des espèces gazeuses à catalyser (comme vu ci-dessus), que l'on peut caractériser par un débit massique d'espèces gazeuses à catalyser que l'on peut noter QpoI. [0035] Le catalyseur 21 présente une longueur Lcat, et l'invention propose d'estimer (par calcul ou simulation) la température du catalyseur 21 10 sur sa longueur, c'est-à-dire le long d'un axe x représenté figure 2. [0036] En effet, il est important de connaître la température du catalyseur, car l'efficacité de catalyse est fortement dépendante de la température de paroi Ts du catalyseur. Par exemple, et de manière non limitative, on peut considérer que l'efficacité de catalyse n'est que de 50% à 15 une température de 600 K environ, supérieure à 80% pour une température de 700 K, et supérieure à 90% pour une température de 800 K. [0037] Il est donc important de savoir si l'ensemble du catalyseur est bien au dessus d'une certaine température seuil (par exemple, et de manière non limitative 750 K) pour garantir que le moteur thermique 10 peut 20 être commandé selon un mode de fonctionnement nominal. Par contre, si le catalyseur 21 n'est pas ou pas totalement monté en température, on peut envisager de mettre en place soit une stratégie de commande du moteur thermique pour faire monter l'ensemble du catalyseur 21 en température (en jouant sur la richesse ou sur l'avance à l'allumage, avec des commandes 25 spécifiques du dispositif d'allumage 15 ou du dispositif d'injection 14), ou alors moins utiliser le moteur thermique 10 (en privilégiant un moteur de traction électrique pour un véhicule hybride par exemple). [0038] A cet effet, l'invention propose de prédire la température des gaz d'échappement et de la paroi du catalyseur 21, tout le long de l'axe x 30 de la figure 2, pour ensuite proposer la commande adéquate du moteur thermique 10. -16- [0039] Comme vu ci-dessus, la température de gaz d'échappement Tg en une pluralité de points le long du catalyseur est calculée avec une première fonction fi qui prend en compte les différentes variables de la manière suivante : 5 7:9(X, = 1.1(7:9(X, t - dt), 7:0(X - dx, t - dt), 7:0(X dx, t - dt), t - dt), Qech(t - dt)) La température de paroi du catalyseur Ts en une pluralité de points le long du catalyseur est calculée quant à elle avec une deuxième fonction f2 qui prend en compte les différentes variables de la manière suivante : 10 t) = f2(Ts(x, t - dt), Ts(x - dx, t - dt), Ts(x + dx, t - dt), 7:0 (X, t - dt), Tamb (t - dt), Qp 01(x, t - dt), Qech(t - dt)) Avec : x: position le long du catalyseur (mm), t: instant (s) 15 Tg : température de gaz d'échappement (K) Ts : température de paroi du catalyseur (K) Tamb : température ambiante au système (K) Qech : débit massique en amont du catalyseur (Kg/s), obtenu par mesure ou par estimation avec un modèle issu du contrôle moteur 20 Qpoi : débit massique d'espèces gazeuses à catalyser (par exemple, monoxyde de carbone (CO), oxydes d'azote (N0x), hydrocarbures imbrûlés (HC)...) Tg(0, t) = Tech (t) : température de gaz d'échappement en amont du catalyseur (K), obtenue par mesure ou par estimation avec un modèle issu du contrôle moteur 25 Tg(x, 0) = Tim : température initiale du système (K). Ts(x, 0) = Tim : température initiale du système (K). [0040] Le détail des fonctions et calculs est exposé ci-dessus et ne sera pas revu ici. À l'aide d'un calcul itératif convergent, il est donc possible de connaître en tout point du catalyseur 21 la température de paroi et donc 30 l'efficacité globale de traitement des espèces gazeuses à catalyser. -17- [0041] La figure 3 représente un exemple de structure de traitement de l'information de calcul de la température du catalyseur 21. Un calculateur 110 estime la température de gaz et de paroi (E-T°) le long du catalyseur 21, en recevant comme données d'entrée le débit massique de gaz 5 d'échappement, la température de gaz en amont du catalyseur 21, le débit massique des espèces gazeuses à catalyser (ces données peuvent être issues de mesure, ou être issues de calculs effectués par des calculateurs moteur utilisant des modèles de prédiction). [0042] La sonde aval de mesure de température 23 envoie 10 également une mesure de température de gaz effectuée en aval du catalyseur 21, pour corriger si nécessaire les estimations de température de gaz et de paroi. [0043] En sortie, le calculateur 110 envoie les valeurs de température du catalyseur 21 estimées et/ou corrigées selon l'axe des x à un 15 calculateur 120, qui estime l'efficacité globale n (E-Poil) de traitement des espèces gazeuses à catalyser. [0044] En fonction de cette efficacité (ou directement en fonction de la température du catalyseur 21 estimée et/ou corrigée), un superviseur 130 (S) peut alors mettre en place une stratégie visant à répartir un ratio de 20 puissance entre un moteur de traction électrique (Pmei) et le moteur thermique (Pmth), ou encore instruire spécifiquement un calculateur moteur 140 (C-Mth) pour piloter le moteur thermique selon un mode nominal ou selon une stratégie spécifique pour générer des gaz d'échappement plus chauds, de sorte à faire monter en température le catalyseur 21 pour se placer dans un domaine de 25 catalyse efficace. [0045] On comprendra que diverses modifications et/ou améliorations évidentes pour l'homme du métier peuvent être apportées aux différents modes de réalisation de l'invention décrits dans la présente description sans sortir du cadre de l'invention REVENDICATIONS1. Procédé de commande d'un moteur thermique (10) de véhicule équipé : - d'un moteur thermique (10), - d'au moins un catalyseur (21) pour traiter des gaz d'échappement du moteur thermique (10), - d'au moins une sonde aval de température (23) de gaz d'échappement agencée en aval du catalyseur (21), 10 le procédé de commande comprenant les étapes consistant à : - estimer, à un instant (t) du fonctionnement du moteur thermique (10), une température de gaz d'échappement en une pluralité de points le long du catalyseur (21), en prenant en compte au moins : - une valeur de température initiale du catalyseur (21), 15 - une valeur de température de gaz d'échappement (Tg) en amont du catalyseur (21) à l'instant (t) et/ou à un instant précédent t-dt, et - un débit massique de gaz d'échappement (Qech), à l'instant précédent t-dt uniquement, - estimer à l'instant (t) une température de paroi (Ts) du catalyseur (21) en une 20 pluralité de points le long du catalyseur (21), en prenant en compte au moins la température de gaz d'échappement (Tg) le long du catalyseur (21), - mesurer avec la sonde aval de température (23) à l'instant t, une température de gaz d'échappement en aval du catalyseur (21), - corriger la température de paroi (Ts) du catalyseur (21) le long du catalyseur 25 (21), en prenant en compte la mesure de température de gaz d'échappement en aval du catalyseur (21), si une erreur d'estimation de la température de gaz d'échappement (Tg) est détectée, - commander au moins un organe (14, 15) du moteur thermique (10) en -19- fonction de la température, estimée ou corrigée, de paroi (Ts) du catalyseur (21) le long du catalyseur (21). 2. Procédé de commande selon la revendication précédente, dans lequel la température de paroi (Ts) du catalyseur (21) en une pluralité de points 5 le long du catalyseur (21) à l'instant t est estimée en prenant également en compte une température ambiante au véhicule à l'instant (t) et/ou à l'instant précédent t-dt, et/ou un débit massique d'espèces gazeuses à catalyser (QpoI) à l'instant précédent t-dt. 3. Procédé de commande selon l'une des revendications précédentes, comprenant une étape consistant à : - estimer un débit massique local d'espèces gazeuses à catalyser (QpoI) en une pluralité de points le long du catalyseur (21) à l'instant (t), en fonction de la température de paroi (Ts) corrigée. 4. Procédé de commande selon la revendication précédente, dans lequel l'étape, consistant à estimer le débit massique local d'espèces gazeuses à catalyser (Q0) en une pluralité de points le long du catalyseur (21) à l'instant (t), comprend une étape consistant à estimer le débit d'espèces gazeuses à catalyser (QpoI) en amont du catalyseur (21), en fonction du débit massique de gaz d'échappement à l'instant précédent t-dt, et/ou d'un régime du moteur thermique (10) à l'instant précédent t-dt. 5. Procédé de commande selon l'une des revendications précédentes, dans lequel l'étape consistant à corriger la température de paroi (Ts) du catalyseur (21) le long du catalyseur (21) comprend les étapes consistant à : - estimer à l'instant (t) une température de gaz d'échappement en extrémité aval (Taval) du catalyseur (21), en fonction : de la température de gaz d'échappement en extrémité aval (Tavai) du catalyseur (21) à l'instant précédent t-dt, de la température de paroi (Ts) du catalyseur (21) en extrémité aval du catalyseur (21) à l'instant précédent t-dt, -20- de la température ambiante au véhicule à l'instant précédent t-dt, - comparer la température de gaz d'échappement en extrémité aval (Tavai) du catalyseur (21) estimée à l'instant t avec la mesure à l'instant t de la température de gaz d'échappement en aval (Taval) du catalyseur (21), 5 - si une différence entre la température de gaz d'échappement en extrémité aval (Tavai) du catalyseur (21) estimée à l'instant t et la mesure à l'instant t de la température de gaz d'échappement en aval (Taval) du catalyseur (21) excède un seuil prédéterminé : - corriger la température de gaz d'échappement (Tg) estimée à l'instant 10 (t) en une pluralité de points le long du catalyseur (21), - ré-estimer à l'instant (t) la température de paroi (Ts) du catalyseur (21) en une pluralité de points le long du catalyseur (21), en prenant en compte au moins la température de gaz d'échappement (Tg) corrigée le long du catalyseur (21) . 15 6. Procédé de commande selon l'une des revendications précédentes, dans lequel, si la température de paroi (Ts) du catalyseur (21) estimée en l'un des points le long du catalyseur (21) est inférieure à une température seuil, alors au moins un des organes du moteur thermique (10), tel qu'un dispositif d'injection et/ou un dispositif d'allumage, est commandé 20 pour augmenter ladite température de paroi (Ts) du catalyseur (21). 7. Procédé de commande selon l'une des revendications précédentes, dans lequel : - un intervalle de temps entre l'instant t et l'instant t-dt est inférieur ou égal à 5 ms, et avantageusement inférieur ou égal à 1 ms, et/ou 25 - un intervalle géométrique entre deux points, auxquels la température de paroi (Ts) du catalyseur (21) est estimée, est inférieur ou égal à 5 mm, et avantageusement inférieur ou égal à 1 mm. 8. Procédé de commande selon l'une des revendications précédentes, dans lequel le véhicule est équipé d'une sonde amont de -21 - température (22) de gaz d'échappement agencée en amont du catalyseur (21), et dans lequel : - la valeur de température initiale du catalyseur (21) est une température mesurée par la sonde amont de température (22) de gaz d'échappement 5 et/ou la sonde aval de température (23) de gaz d'échappement, avant mise en fonction du moteur thermique (10), et/ou - la valeur de température de gaz d'échappement en amont du catalyseur (21) à l'instant t et/ou à un instant précédent t-dt est une température mesurée par la sonde amont de température (22) de gaz d'échappement. 10 9. Moteur thermique (10) de véhicule automobile, comprenant au moins un organe commandé par le procédé de commande selon l'une des revendications précédentes. 10. Véhicule automobile thermique ou hybride, comprenant au moins un moteur thermique (10) selon la revendication précédente. 15