Brevet FR3082002A1 — publié le 20191206

PROCEDE DE DETERMINATION D'UN MODELE DE LA COMPOSITION DES GAZ D'ECHAPPEMENT D'UN MOTEUR DIESEL La présente invention concerne la phase de développement de moteurs à combustion interne de véhicules automobiles et en particulier de moteur Diesel. Les normes réglementaires sur la pollution devenant de plus en plus sévères, ceci implique d'être toujours plus précis dans la détermination des paramètres optimaux permettant de minimiser les émissions polluantes. Pendant la phase de développement de moteurs à combustion interne de véhicules automobiles, on utilise des plans d'expérience. Pour chaque point de fonctionnement moteur du plan d'expérience, on modélise les émissions polluantes (CO, 002, NOx, HO) en fonction des paramètres d'entrées de calibration de ces points de fonctionnement, tels que l'avance à l'injection, le débit d'air, la pression rail, la pression de suralimentation, le débit de gaz récirculés, etc...). Cette modélisation peut se faire via des modèles quadratiques mathématiques de la forme : Y = ao + a1X1 + a2X2 + a3X3 + a11X + a22JC + a33X12_ + a12X1X2 + a13X1X2 + a23X2X3 Avec Xi les paramètres d'entrées de calibration du point de fonctionnement et Y un polluant. Ces modèles quadratiques sont classiquement obtenus à partir de transformations mathématiques génériques. Cependant, les transformations utilisées sont complexes et décorrélées des phénomènes physiques de la combustion. De plus, la manipulation par les metteurs aux points des coefficients des polynômes n'est pas intuitive. Par conséquent, le problème à la base de l'invention est de proposer, à partir de mesures de polluants sur les points de fonctionnement moteur d'un plan d'expérience, une modélisation simple mais qui prend en compte la dépendance physique entre les différents polluants mesurés, en particulier dans le cas d'une combustion dans un moteur Diesel. Pour résoudre ce problème, il est prévu selon l'invention un procédé de détermination d'un modèle de la composition des gaz d'échappement d'un moteur Diesel comprenant : 2 -une étape de mesures de la composition des gaz d'échappement sur un ensemble de points de fonctionnement du moteur, ces points étant définis par un ensemble de paramètres d'entrée de fonctionnement, et dans lequel : 5 -on détermine, pour chacun des points de fonctionnement, les coefficients d'équilibre d'une équation bilan de combustion prédéterminée à partir des mesures de composition des gaz d'échappement, -on détermine, pour chaque point de fonctionnement, des paramètres chimiques en fonction de ces coefficients d'équilibre et de l'équation bilan de combustion prédéterminée dans 10 laquelle, la réaction de combustion étant vue comme la somme d'une zone riche et d'une zone pauvre de combustion, les coefficients d'équilibre sont exprimés en fonction de ces paramètres chimiques, -on exprime la composition des gaz d'échappement en fonction de ces paramètres chimiques, 15 -on détermine, à partir des valeurs des paramètres chimiques obtenues pour les points de fonctionnement, un modèle quadratique de ces paramètres chimiques en fonction des paramètres de fonctionnement, -On détermine le modèle de composition des gaz d'échappement en remplaçant dans l'expression de la composition des gaz d'échappement les paramètres chimiques par leur 20 modèle quadratique. L'effet technique est d'avoir un modèle plus représentatif de la réalité en ce qu'il prend en compte les interdépendances entre les polluants. 25 Diverses caractéristiques supplémentaires peuvent être prévues, seules ou en combinaisons : Selon une réalisation, l'ensemble de point de fonctionnement moteur forme un plan d'expérience, ce qui permet de limiter au juste nécessaire le nombre d'essais physiques à 30 réaliser. Selon une réalisation, les paramètres chimiques de l'équation bilan sont : la proportion de la zone riche, aR, la richesse de la zone riche, R, la richesse de la zone pauvre L, la proportion de dioxyde de carbone dans la zone riche fi, la proportion de monoxyde de carbone dans la zone riche 8, la proportion d'hydrocarbures imbrûlés dans la zone pauvre y, la proportion de suie dans la zone riche C. 3 Selon une réalisation, la mesure de la composition des gaz d'échappement est réalisée par une baie d'analyse donnant la proportion du dioxyde de carbone, 002, du monoxyde de carbone, CO, des hydrocarbures imbrûlés, HO, de l'oxygène, 02, des oxydes d'azote, NOx, dans un volume de fumées sèche, VFS. 5 Selon une réalisation, l'équation bilan de combustion prédéterminée est : CHyOz + n(02+ xN2) a CO2+ b CO + d C'HvOw + eH2 + f 02+ g H20 + h N Ox + i N2 + t C8Fli -Avec des hydrocarbures imbrûlés modélisés sous la forme générale CuH2Ow, et des suies 10 modélisées sous la forme 05H1, -Avec n, a, b, d, e, f, g, h, i, t, les coefficients d'équilibre respectivement de l'air, du gaz carbonique, 002, du monoxyde de carbone CO, des hydrocarbures imbrûlés, de l'hydrogène H2, de l'oxygène en excès 02, de l'eau H20, des oxydes d'azotes NOx, de l'azote N2, des suies, et x la proportion d'azote dans l'air. 15 Selon une réalisation, les coefficients d'équilibre, n, a, b, d, e, f, g, h, i, t, sont déterminés par itération en : a) initialisant les coefficients d'équilibre avec des valeurs aléatoires, b) définissant la suite VFSk pour tout k>1, tel que VFSk+i = ak + bk + dk + ek + fk + hk + 20 ik, c) déterminant le rang k+1 des suites (a)k,(b)k,(d)k,(e)k,(f)k,(g)k,(h)k,(i)k,(t)k grâce aux relations : VFSk VFSk VFSk FSk VFSk ak - 100 %Vc02,bk - 100 %Vco, dk - k k 100 = 100 %V02, hk - -%VNox 100 et : akek hk 1 - ak - bk - Udk gk = Kx1- ik = Xnk - ;tk = 2 "k 8 Avec Kx la constant d'équilibre du gaz à l'eau, d) en itérant sur les étapes b) et c) jusqu'à ce que les variations entre deux itérations deviennent inférieures à un seuil déterminé. 25 30 D'autres particularités et avantages apparaîtront à la lecture de la description ci-après d'un mode particulier de réalisation, non limitatif de l'invention. 4 La présente invention a pour but de modéliser de manière interdépendante les émissions polluantes. Cette invention s'applique à la mise au point de moteur à allumage par compression, de type Diesel. 5 Selon cette invention, on réalise des mesures de polluants sur un ensemble de points de fonctionnement du moteur en régime et charge stabilisés, ces points de fonctionnement moteur formant un plan d'expérience. Chacun des points de fonctionnement est défini par des paramètres d'entrées de 10 fonctionnement Xi, tels que l'avance à l'injection, le débit d'air, la pression rail, la pression de suralimentation, le débit de gaz récirculés, etc...). En pratique, ces mesures de polluants peuvent être effectuées à l'aide d'une baie d'analyse cinq gaz. Une baie d'analyse cinq gaz fournit pour le gaz carbonique, CO2, le monoxyde de 15 carbone, CO, les hydrocarbures imbrûlés, HC, l'oxygène 02 et des oxydes d'azote NOx, leurs proportions volumiques en pourcentage de vapeur sèche, respectivement (3/0Vc02, cY0Vco, cY0VHc, %V02, °/oVN0N- dans un volume de vapeur (ou fumées) sèche(s) que nous noterons VFS (pour Volume de Fumées Sèches). 20 On définit ensuite une équation de bilan de combustion : CHyOz + n(02 + xN2) a CO2 + b CO + d CuH2Ow + eH2 + f 02 + g H20 + h NOx + i N2 + t C8111 (1) 25 Avec n, a, b, d, e, f, g, h, i, t, les coefficients permettant d'équilibrer l'équation de combustion pour respectivement l'air, le gaz carbonique CO2, le monoxyde de carbone CO, les hydrocarbures imbrûlés, HC, sont modélisés sous la forme générale CuH2Ow, l'hydrogène H2, l'oxygène en excès 02, l'eau , les oxydes d'azotes, NON, considérés ici sous leur forme la plus générale possible, l'azote N2 et les suies modélisée ici sous la forme C8I-11. 30 La variable x, liée à la composition de l'air est connue. Pour les variables u, y, w , liée à la nature des hydrocarbures imbrûlés, on fait une hypothèse qui permet de les écarter des inconnues à déterminer. Une hypothèse possible est de considérer que la formulation des hydrocarbures imbrûlés est la même que celle du carburant initial c'est-à-dire que Culi3O, 35 = CHyOz. D'autre hypothèses sont possibles, par exemple prendre les hydrocarbures imbrûlés sur base de méthane CH4. Pour la variable x, liés à la répartition NO/NO2, on fait également une hypothèse qui permet de les écarter des inconnues à déterminer. Une 5 hypothèse possible est de considérer que l'on a que du NO et pas de NO2, soit x = 1. D'autre hypothèses sont possibles, par exemple considérer classiquement une répartition de 80% de NO et de 20% de NO2. 5 On décrit ensuite la méthode pour obtenir les inconnues abdef ghi et t via l'analyse des résultats de la baie d'analyse des gaz. Comme déjà dit, cette baie exprime la proportion des polluants CO2, CO, HC, 02 et NOx dans un volume de vapeur (ou fumées) sèche(s) que nous noterons VFS (pour Volume de Fumées Sèches). Ce volume ne prend pas en compte : 10 - l'eau, car il s'agit d'un volume sec. - les particules C81-11, absentes de la phase gazeuse car à l'état solide. Le Volume de Fumées Sèches, VFS peut ainsi s'écrire : VFS=a+b+d+e+f+h+i (2) 15 Avec a = --Yovco2, = d = %Viic, f = h = %Vco, too t00%17Nox VFS VFS VFS VFS too too 100*u (3) La formulation particulière obtenue pour les HC provient de la spécificité de la mesure des 20 HC sur la baie qui sont convertis en équivalent CH4. Pour exemple : en envoyant un échantillon à 100ppm de C3H8, la baie retournera une valeur (3/0VHc de 300ppm. La quantité d doit donc tenir compte de la formulation des HC choisi. A ce stade, nous avons 10 inconnues (n, a, b, d, e, f, g, h, i, et t), pour 5 équations. 25 L'équilibre du gaz à l'eau permet l'obtention d'une sixième équation : CO2 + H2 e co + H20 (4) Cet équilibre est fortement lié à la température. Néanmoins, il est admis que la réaction est bloquée cinétiquement à partir de 1700K lorsque la chambre de combustion se refroidit. La 30 constante d'équilibre a alors pour valeur Kx = 3,5 : Kx = = 3,5 (5) Ce qui permet d'en déduire g: g = Kxceb' (6) Les quatre dernières équations pour fermer le système sont les équations de conservation 35 atomique. 6 Pour le carbone : 1=a+ b+ud+8t. La quantité t de suie s'en déduit : (7) 1-a-b-ud t= 8 Pour l'hydrogène : y = dy + 2e + 2g +t. On en déduit e en introduisant les quantités précédemment calculées : y-t-vd e= a 2+2-bKx Pour l'azote : xn = i + On en déduit la quantité i d'azote final : h = 10 Enfin la conservation atomique de l'oxygène donne : 2n + z = 2a +b+2f+g+xh+wd. D'où : n=a++f +1+.1h+vd- (10) 2 2 2 2 2 Ainsi avec 10 équations - 5 mesures de la baie, 4 équations de conservation atomique et 1 équation sur l'équilibre du gaz à l'eau - nous déterminons les 10 inconnues du système : 9 15 quantités portant sur les produits (a, b, d, e, f, g, h, i, et t) et 1 quantité portant sur l'air initialement disponible (n). La connaissance de la quantité d'air disponible nous permet de calculer notre propre richesse à partir de la mesure des 5 gaz (002, CO, HO, 02 et NON). Par définition de la 20 richesse, : Où mcarb et /flair Cncarb ,- ) Mau- reel ncarb M air (i) = rn = IP S ,,,, (-carbair I",Ai carb ) , e Mair stoech sont respectivement les masses de carburant et d'air, n -carb et nain le nombre de moles de carburant et d'air, Ips est le pouvoir comburivore du carburant s'exprimant : = (1 + z) M carb 4 2 M air 25 Ceci conduit à: =ncarb pair+ - Dans notre cas, nous consommons 1 mole de carburant pour n moles d'air : = (1 + n = (1 + y z 1 4 2) n ) 30 4 2 2 22 2 2 5 (8) (9) 7 Lors du post-traitement des résultats de la baie d'analyse, la difficulté est l'obtention des quantités a, b, d, f, et h, qui nécessite la connaissance de la valeur de VFS. Or VFS ne peut être calculé qu'après avoir déterminé les différentes quantités a, b, d, e, f, h, et i. 5 On utilise alors une approche itérative qui consiste à initialiser avec des valeurs aléatoires le système a, b, d, e, f, g, h, i, t. Le volume de fumées sèches, VFS, est ensuite calculé à partir de ces valeurs. Connaissant VFS, la résolution du système à 10 équations permet de retrouver a, b, d, e, f, g, h, i, t, et n. 10 Empiriquement, cette boucle itérative converge relativement rapidement (entre dix et vingt itérations). L'algorithme peut être décrit de la manière suivante. 15 1) On initialise le vecteur [cto, bo, do, eo, fo, go, ho, io, to] 2) On définit la suite (VFS)k pour tout k>1, tel que VFSk+1 = ak + bk + dk + ek + fk + hk + ik 3) On détermine le rang k+1 des suites (a)k,(b)k,(d)k,(e)k,(f)k,(g)k,(h)k,(i)k,(t)k grâce aux relations : 20 VFSk VFSk VFSk VFSk VFSk ak = %Vc02,bk - %Vco, dk = %17Hc, - - %V02, hk - %17Nox 100 100 100 100 100 Et: akek hk 1 - ak - bk - udk g k = Kx- - n k /1 y = Uk2 8 25 Et n à partir de : nk k itère sur les étapes 2) et 3) jusqu'à convergence. Lorsque les variations de chaque suite (a)k,(b)k,(d)k,(e)k,(f)k, (g)k,(h)k,(i)k,(t)k, entre deux pas, deviennent inférieures à eliFS = 10' (valeur arbitrairement petite), on stoppe la boucle itérative. 30 Dans une autre étape, nous construisons un physico-chimique de la réaction. Ce modèle présuppose que la combustion Diesel est fortement hétérogène et globalement pauvre. Nous faisons l'hypothèse que cette hétérogénéité peut être vue comme la somme de deux 35 zones homogènes : l'une riche, l'autre pauvre. 8 Concernant la zone pauvre, l'équation globale classique d'un mélange pauvre s'écrit : y z Y CHyOz + nL(02 + xN2) CO2 + - (1 + -4 - -2)) 02 -2H20 + xni,N2 (1+H) Avec : = où (h représente la richesse de cette zone pauvre. (PL 5 A cette équation nous allons rajouter deux polluants. Tout d'abord on prend en compte la formation des NOx qui a lieu en zone pauvre. A des fins de simplification, les NOx sont uniquement modélisés sous la forme NO. La réduction du mécanisme de Zeldovitch conduit à l'équation simplifiée: 1 1 - N2 - 02 NO 2 2 10 La production de NO correspond à une réaction parasite se produisant dans une partie réduite de la zone pauvre (la partie la plus chaude, proche de la stoechiométrie). Un facteur aLFIT (LHT pour « Lean High Temperature » en anglais) est introduit pour modéliser la proportion de production de NOx au sein de la zone pauvre. 15 D'autre part, on fait comme hypothèse que la majeure partie des HC sont produits en zone pauvre. Cette hypothèse est introduite pour deux raisons: Les HC sont produits proches des parois et des crevasses (extinction de flamme due à une température trop basse), dans des zones éloignées du jet, donc pauvre Les imbrûlés peuvent également être produits dans les zones extrêmement pauvres 20 en dessous de la limite d'inflammabilité basse. L'hypothèse implicite est que les imbrûlés présents dans les zones très riches (limite d'inflammabilité haute) sont minoritaires. Cette hypothèse est communément admise dans la littérature. Le modèle le plus probant est celui prenant en compte les hydrocarbures 25 imbrûlés, HC, en zone pauvre. On obtient (équation 12) : CHyOz + nL(02 + xN2) + (1 + 7z) + y (u +17; v+i) cruiT) 02 (1 - uy)C 02 + yCutivOw + 30 + (377 - y) H20 + aLHTNO + (xnL a121T) N2 Concernant la zone riche, l'équation usuelle de zone riche est donnée par : CHyOz + nR(02 + xN2) 9 fiC 02 + 8C0 + (2nR + z - 2fl - 8)H20 + (377 - 2nR - z + 2fl + ô ) H2 + XnRN2 (1-g-1) Avec nR = 4 2 , où oR représente la richesse de cette zone riche et avec fi + 8 = 1. OR 5 II convient également de prendre en compte la formation de suies modélisées sous la forme C81-11, ce qui donne (équation 13) : CHyOz + nR(02 + XN2) 13CO2 + 8C0 + (2nR + z 213 8)H20 + (Y c 7 - 2nR - z + 2fl + ô- 7) H2 + X71,RN2 + EC8111 Avec fi + 8 + 8E = 1 10 En multipliant chaque zone par sa proportion, aL pour la zone pauvre et aR pour la zone riche, on obtient une équation globale finale prenant en compte la coexistence d'une combustion en milieu riche et en milieu pauvre. On notera que aL = 1 - aR. 15 L'équation bilan est alors (équation 14) : CHyOz + ((1 - aR)ni, + aRnR)(02 + xN2) ((1 - aR)(1 - uy) + aRfl)C 02 + aR8C0 + (1 - aR)yCull,O, +[aR(77 - 2nR - z + 2fl + 8 - 76.)]112 + (1 - aR) (ni, - (1 + 3774 7z) + y (u +1.9 ID aulT) 02 + [(1 - aR)(7 - 71 9 y) + aR(2nR + z - 2fl - 8)1H20 + (1 - aR)aLHTNO + aRXnRN2 +[(1 - aR)(xnL aulT) 1 2 / + aRec8111 Les paramètres chimiques à modéliser sont donc : aR : Proportion de la zone riche, 20 (PR : Richesse de la zone riche, Richesse de la zone pauvre, 13 : Proportion de CO2 dans la zone riche, Ô: Proportion de CO dans la zone riche, y : Proportion de HO dans la zone pauvre, 25 E : Proportion de suie dans la zone riche. 10 La baie d'analyse nous fournit une indication sur cinq polluants. Nous voulons modéliser les sept paramètres chimiques précités. L'équation de conservation des atomes de carbone en zone riche (fl + 8 +8c = 1) nous fournit une équation supplémentaire. Il nous manque 5 néanmoins une équation pour fermer le système. On choisit dans ce premier modèle d'imposer la richesse en zone pauvre comme étant le ratio carburant/air en ne prenant pas en compte la formation de CO et HO (avec x=1 car NO x = NO), ceci nous fournit la septième équation nécessaire : 4 2) a+f+H y zh- a (15) 10 II ne reste plus que 6 paramètres chimiques à modéliser. On évalue ces paramètres chimiques en identifiant l'équation bilan initiale (équation 1) à l'équation 14, on obtient en identifiant sur 02 : 1 1 aR = 1 - (f + -2 h - d (u + v w)) nL-(1+11-1 4 2 4 2) (16) En identifiant sur NOx : 15 En identifiant sur 05H1 : En identifiant sur CO2 : 20 En identifiant sur CO : h aLHT = l_aR t E = - aR = 1 + a-l+ud fi' aR (17) (18) (19) (20) En identifiant sur les hydrocarbures imbrûlés: d Y = 1-an On notera que l'on vérifie bien que [3+15+8E = a+b+8t = 1 En identifiant sur H20 : g-(1-aR)(112-'22y) 1 1 nR = 2 Z 2aR D'où la richesse de la zone riche : 25 (21) (22) 11 y z (23) 1+7-7 - g - (1- a Re -Ey) 2 2 2CCR 2 r- 2 Ces différentes équations permettent d'identifier les coefficients aR, (th, p, 8,y et e sur l'ensemble des points de fonctionnement du moteur joués en plans d'expérience. Une fois leurs valeurs obtenues, ils peuvent être modélisés sous forme de polynômes comme dans 5 l'approche classique, autrement dit par des modèles quadratiques déjà évoqués. aR = (Xi);fi = f2(Xi);ô = f3(X1); etc ... Lorsqu'on récupère la modélisation des coefficients aR,0R,(pL, fi, 8,y et e, la transformation inverse doit être effectuée pour récupérer la concentration des différents polluants. On 10 utilise les équations suivantes. y-t-vd = y-aR(e-vy) A partir de l'équation 8: e = 2+2K1-ccR)+cco bX 2+2 aR8 A partir de l'équation 6 : g = K -ae = Kx(1-aR)+aRre X 31-aR(e-vY) x b Le volume de vapeur sèche est proportionnel à (équation 24) : Y DSV = (1- aR)(1 - uy) + aRfl + aR8 + (1 - aR)y + aR (7- 2nR - z +2fl + - 7E) D+ + (1 - aR)()(nL aLHT) aRXnR 2 Les produits de combustion sont obtenus en pourcentage de volume de vapeur sèche via les formules (équation 27) : 100 %vco2 = -[(1- otR)(1 - uy)+ aRfil 100 %17c° i;sv [aR8] 100 %viic =,517[(1-c(R)Y] 100 %17Nox - Dsv[(1 aR)auul 2 100 %Vo z v w 1 y 2 4 2 2 4 DSV aRS 2+2 IL-ccR)+ccRfl aR8 15 20 25 12 En utilisant dans ces équations l'expression des coefficients aR, oR, (pL,fl,(5,y et E par leur modèles quadratiques on peut ensuite retrouver la composition des gaz d'échappement à la source, en fonction des valeurs de paramètres de calibration Xi que l'on souhaite tester. 5 Avec cette invention, la modélisation des polluants étant liée, les compromis à réaliser prennent naturellement en compte les interdépendances et permettent de mieux répondre aux exigences réglementaires en vigueur. 10 L'invention a pour avantage de réduire le nombre d'essais nécessaire pour un plan d'expérience classique. En effet les modélisations obtenues sont non linéaires et ne correspondent pas aux modèles quadratiques initiaux. Les nouveaux modèles obtenus sont donc d'ordres supérieurs aux modèles quadratiques classiques et devraient à ce titre nécessiter des plans d'expérience plus volumineux, mais la non linéarité des modèles 15 proposés est issue de formules précises, et ce sont les coefficients aR, oR, (pL,fl,8,y REVENDICATIONS1. Procédé de détermination d'un modèle de la composition des gaz d'échappement d'un moteur Diesel comprenant : -une étape de mesures de la composition des gaz d'échappement (%17o02 ; %Vco; %Viic; %17Nc,,,; %V02) sur un ensemble de points de fonctionnement du moteur, ces points étant définis par un ensemble de paramètres d'entrée de fonctionnement (Xi), caractérisé en ce que : -on détermine, pour chacun des points de fonctionnement, les coefficients d'équilibre (n, a, b, d, e, f, g, h, i, t) d'une équation bilan de combustion (1) prédéterminée à partir des mesures de composition des gaz d'échappement (%17c02 ; %17c0; %Vilc; %Vivo,» %V02), -on détermine, pour chaque point de fonctionnement, des paramètres chimiques en fonction de ces coefficients d'équilibre (n, a, b, d, e, f, g, h, i, t) et de l'équation bilan (1) de combustion prédéterminée dans laquelle, la réaction de combustion étant vue comme la somme d'une zone riche et d'une zone pauvre de combustion, les coefficients d'équilibre sont exprimés en fonction de ces paramètres chimiques (aR, oR, (h, fi, ô, y, c.), -on exprime la composition des gaz d'échappement (%17c02 ; %17c0; %Viic; %V/vox; %V02) en fonction de ces paramètres chimiques (aR,oR, (th, fi, 8,y, c.), -on détermine, à partir des valeurs des paramètres chimiques (aR, 0 R, 0L, fi, 8,y, c) obtenues pour les points de fonctionnement, un modèle quadratique de ces paramètres chimiques (aR,oR, (th, fl , 8,y, c. ) en fonction des paramètres de fonctionnement (Xi), -on détermine le modèle de composition des gaz d'échappement en remplaçant dans l'expression de la composition des gaz d'échappement (%17o02 ; %Voo; %Viic; %VN,,x; %V02) les paramètres chimiques (aR, oR,(pL, p, 8,y, c.) par leur modèle quadratique. 2. Procédé selon la revendication 1, caractérisé en ce que l'ensemble de point de fonctionnement moteur forme un plan d'expérience. 3. Procédé selon la revendication 1 ou la revendication 2, caractérisé en ce que les paramètres chimiques de l'équation bilan sont : la proportion de la zone riche (aR), la richesse de la zone riche (OR), la richesse de la zone pauvre (oL), la proportion de dioxyde de carbone dans la zone riche (fi), la proportion de monoxyde de carbone dans la zone riche (8), la proportion d'hydrocarbures imbrûlés dans la zone pauvre (y), la proportion de suie dans la zone riche (E). 14 4. Procédé selon l'une des revendication 1 à 3, caractérisé en ce que la mesure de la composition des gaz d'échappement est réalisée par une baie d'analyse donnant la proportion du dioxyde de carbone (002), du monoxyde de carbone (CO), des hydrocarbures imbrûlés (HO), de l'oxygène (02), des oxydes d'azote (N0x) dans un volume 5 de fumées sèche (VFS). 5. Procédé selon la revendication 4, caractérisé en ce que l'équation bilan de combustion prédéterminée est : CHyOz + n(02+ xN2) a CO2+ b CO + d C'HvOw + eH2 +f 02+ g H20 + h NOx + i N2 10 + t C81-11 -Avec des hydrocarbures imbrûlés modélisés sous la forme générale CuH2Ow, et des suies modélisées sous la forme 05H1, -Avec n, a, b, d, e, f, g, h, i, t, les coefficients d'équilibre respectivement de l'air, du gaz carbonique (002), du monoxyde de carbone (CO), des hydrocarbures imbrûlés, de 15 l'hydrogène (H2), de l'oxygène en excès (02), de l'eau (H20), des oxydes d'azotes (N0x), de l'azote (N2), des suies, et x la proportion d'azote dans l'air. 6. Procédé selon la revendication 5, caractérisé en ce que les coefficients d'équilibre (n, a, b, d, e, f, g, h, i, t) sont déterminés par itération en : a) initialisant les coefficients d'équilibre avec des valeurs aléatoires, 20 b) définissant la suite (VFS)k pour tout k>1, tel que VFSk+i = ak + bk + dk + ek + fk + hk + ik, c) déterminant le rang k+1 des suites (a)k, (b)k, (d)k, (e)k, (f)k, (g)k, (h)k, (i)k, (0 k grâce aux relations : VFSk VFS k _7FSk 7FSk VFSk 100%7c02,bk - 100%7c0,'k -100%7Hc, 10o%V02,hk =100 25 et : akek hk 1 - ak - bk - Udk gk = Kx-1, ik = Xnk - -;tk = u 2 k 8 Avec Kx la constant d'équilibre du gaz à l'eau, d) en itérant sur les étapes b) et c) jusqu'à ce que les variations entre deux itérations 30 deviennent inférieures à un seuil déterminé. ak