FR 2928417 A1 20090911 FR 0851477 A 20080306 f PROCEDE DE DETERMINATION DE LA QUANTITE DE CARBURANT A INJECTER AU DEMARRAGE D'UN MOTEUR A INJECTION INDIRECTE [0001 ] La présente invention concerne un procédé de détermination de la quantité de 5 carburant à injecter au démarrage d'un moteur à injection indirecte. [0002 ] La combustion par les moteurs montés dans les véhicules automobiles produit des émissions polluantes. Ces polluants sont susceptibles d'être libérés dans l'atmosphère et de dégrader l'environnement. Des normes anti-pollution de plus en plus sévères sont donc mises en place. notamment pour les moteurs à allumage commandé et injection indirecte. Afin de 10 respecter les normes, des systèmes de post-traitement des gaz d'échappement se sont développés sur les véhicules. Les systèmes de post-traitement catalytiques sont basés principalement sur l'emploi de métaux précieux et permettent de réduire de l'émission de polluants hors du véhicule. [0003 ] Cependant_ le bon fonctionnement des catalyseurs n'est assuré que pour des 15 températures supérieures à une température d'amorçage. De ce fait.. lors du démarrage du moteur à froid, le traitement des polluants n'est pas optimum. Ainsi, 80% des émissions polluantes du véhicule sont émises pendant les 40 à 80 premières secondes après le démarrage du moteur. [0004] Dans ces conditions, il apparaît intéressant de diminuer la quantité de polluants émis 20 au démarrage du véhicule. Un des moyens est d'utiliser le moins de carburant possible pour le démarrage du moteur. [0005 1 Mais, la détermination de la quantité minimale de carburant à injecter pour assurer le démarrage du moteur est délicate. Cette évaluation dépend en effet beaucoup des caractéristiques du carburant. En particulier, une des difficultés provient de la variation 25 importante de la composition chimique du carburant d'un pays à l'autre et selon la saison considérée. [0006 ] Actuellement, pour assurer un démarrage du moteur pour tous les carburants, la quantité de carburant injectée est supérieure à la quantité juste nécessaire qui devrait idéalement être injectée. Une fois que le démarrage a été effectué, des sondes situées dans la ligne 30 d'échappement du véhicule permettent de détecter que trop cle carburant est injecté dans le moteur. Ce défaut est corrigé en diminuant progressivement la quantité de carburant injectée jusqu'à la disparition de la détection du défaut par les sondes. [0007 ] Cette technique qualifiée de � a posteriori implique une sur-consommation de carburant au démarrage par rapport à une situation optimisée. Cette sur-consommation engendre une émission d'une quantité plus importante de polluants par rapport à la situation optimale. [000s 1 Il existe aussi des démarches de détermination dite � a priori � de la quantité de carburant nécessaire au démarrage. Ces méthodes sont basées sur des paramètres de volatilité du carburant comme les points de la courbe de distillation ou toute combinaison linéaire de ces paramètres. Ces derniers procédés sont communément appelés � indice d'agrément de conduite � ou � driveabilit;index � en anglais. [0009 ] Toutefois, les modèles existants utilisant une information sur la volatilité du carburant ne prennent pas bien en considération les interactions entre espèces. Les modèles existants ne sont pas capables de déterminer correctement la quantité de carburant nécessaire au démarrage. Ceci est d'autant plus vrai pour des carburants contenant des composés polaires comme les alcools. Compte tenu de la variabilité des carburants de part le monde. cela engendre des problèmes de robustesse de la prestation de démarrage et de mise en action du moteur : non-démarrage. mauvais démarrage, calage, agrément de conduite dégradé après démarrage. [0010] Il existe donc un besoin pour un procédé qui assure le démarrage d'un moteur à 20 allumage commandé injection indirecte en utilisant la quantité de carburant la plus juste possible. [0002 ] La présente invention a ainsi pour premier objet un procédé de détermination de la quantité de carburant à injecter au démarrage d'un moteur à injection indirecte, un mélange d'une phase gaz d'air et de carburant et d'une phase liquide de carburant étant injecté dans le 25 moteur, ce procédé comprenant le calcul de la composition chimique de la phase gaz et de la phase liquide à l'interface entre la phase gaz et la phase liquide. par utilisation d'un modèle et, à l'aide de la composition chimique de la phase gaz, le calcul du taux de transport de la phase gaz, par utilisation d'un modèle. [0003 ] Dans une variante, associée à un moteur comprenant une chambre de combustion et 30 un conduit d'admission clans le moteur du mélange. la phase liquide étant déposée sur les parois du conduit d'admission. [0004 ] Dans une variante le modèle de calcul de la composition chimique de la phase gaz et de la phase liquide prend en considération la composition chimique du carburant préalablement déterminée. Le nombre d'espèces du carburant pour la détermination de la composition chimique du carburant est réduit. de préférence à quinze. La composition chimique du carburant est de préférence déterminée par un organe d'analyse de laboratoire ou par un organe d'analyse embarqué à bord d'un véhicule. de préférence un organe d'analyse embarqué du type capteur opto-électronique utilisant la spectroscopie proche infrarouge. [0005 ] Dans une variante, le modèle de calcul de la composition chimique de la phase gaz et de la phase liquide est basé sur l'hypothèse des gaz parfaits et la loi de Raoult. [00061 Dans une variante. le modèle de calcul de la composition chimique de la phase gaz et de la phase liquide prend en considération le coefficient de fugacité de chacune des espèces du carburant, le coefficient de fugacité étant calculé ou obtenu dans des tables. [0007 ] Dans une variante, le procédé comprend en outre le calcul de la richesse du mélange obtenue par le rapport entre la valeur de l'intégrale du taux de transport de la phase gaz et la valeur du volume d'air injecté. De façon préférée, le modèle de calcul du taux de transport prend en considération la masse de carburant à injecter, le procédé comprenant en outre la détermination de la quantité de carburant à injecter au démarrage dans le moteur pour parvenir à une richesse prédéterminée du mélange par utilisation de tables ou par itération de la masse de carburant à injecter. [0008 ] Le procédé selon l'invention est plus spécialement utilisé dans un organe de contrôle d'un moteur ou dans un outil d'aide à la mise au point d'un moteur. [0001 ] D'autres caractéristiques et avantages de l'invention apparaîtront à la lecture de la description qui suit des modes de réalisation de l'invention. donnés à titre d'exemple et en référence à la figure annexée qui montre figure 1, un exemple d'organigramme du procédé. [0002 ] Un procédé de détermination de la quantité de carburant à injecter au démarrage est proposé. Ce procédé s'applique à un moteur à injection indirecte dans lequel un mélange de phase gaz d'air et de carburant et d'une phase liquide de carburant est injecté dans le moteur. Le procédé comprend le calcul de la composition chimique de la phase gaz et de la phase liquide à l'interface entre la phase gaz et la phase liquide, par utilisation d'un modèle. De plus, à l'aide de la composition chimique de la phase gaz, le procède calcule le taux de transport de la phase gaz, par utilisation d'un modèle. [00031 Le procédé rend ainsi possible. à partir de la composition chimique détaillée d'un carburant. la détermination de la composition chimique de la phase gaz qui est injectée au niveau de la chambre à combustion du moteur à allumage commandé injection indirecte. En conséquence; la quantité de carburant nécessaire au démarrage d'un moteur à allumage commandé injection indirecte peut être déterminée. [00041 La quantité de carburant injectée au démarrage peut donc être optimisée. De ce fait, le procédé permet de réduire les émissions polluantes de la voiture ce qui diminue la charge en métaux précieux du post-traitement catalytique. Le prix de revient de fabrication de la voiture est ainsi réduit. [0005 ) De plus, la qualité et la robustesse de la prestation démarrage et mise en action est améliorée, notamment parce que cette prestation ne dépend pas du pays ou de la saison considérés. [00061 Le procédé peut comprendre un certain nombre d'étapes comme le montre la figure 1. La figure 1 est un exemple d'organigramme du procédé de détermination de la quantité de carburant à injecter au démarrage. [0007 ] Préalablement à l'étape 6 de calcul de la composition chimique de la phase gazeuse et de la phase liquide, le procédé peut comprendre une étape 2 de détermination de la composition chimique détaillée du carburant. La composition chimique détaillée donne accès aux différentes espèces du carburant ainsi qu'à leur proportion dans le carburant. [0008 1 La composition chimique détaillée peut provenir d'une analyse réalisée par exemple au niveau du réservoir à essence du véhicule par l'une ou toute combinaison des techniques suivantes : chromatographie en phase gazeuse, spectrométrie de masse. résonance magnétique nucléaire (RMN), spectroscopie ultraviolet ou spectroscopie infrarouge. [0009 L'analyse du carburant peut être effectuée par un organe d'analyse de laboratoire. La chromatographie et la spectrométrie de masse peuvent notamment être utilisées dans le laboratoire. De tels instruments permettent en effet d'obtenir une composition chimique du carburant très précise pour la mise au point des moteurs. [00101 La composition chimique du carburant peut aussi être déterminée par un organe d'analyse embarqué à bord d'un véhicule. Cela permet d'obtenir la composition de l'essence dans le véhicule. En particulier, un capteur opto-électronique qui exploite la spectroscopie proche infrarouge (en anglais � near infrared spectroscopy � ou � spectroscopy NIR �) peut être utilisé. Un tel capteur a l'avantage d'être peu encombrant. [0011 1 La composition chimique détaillée ainsi obtenue sert (le donnée d'entrée pour le calcul de la composition dans la phase gaz qui est décrit dans la suite. [0012 1 Dans certains cas, la composition chimique détaillée du carburant peut impliquer un temps de calcul trop important pour effectuer le calcul de la composition dans la phase gaz. Afin de réduire le temps de calcul, il peut ainsi être intéressant de procéder à une étape 4 d'utilisation d'un nombre réduit d'espèces du carburant trouvées dans l'étape 2. [00131 A titre d'exemple, pour représenter une essence commerciale Sans Plomb 95 européenne actuelle, une composition chimique réduite à 15 espèces peut être utilisée : 12 hydrocarbures et 3 composés oxygénés. Le tableau suivant donne les 15 espèces classées par température d'ébullition (T°éb) croissante : Espèce chimique T°éb (K) n-Butane 273 Isopentane 301 n-Pentane 310 Méthy ]butène 312 MTBE (Méthyltertiobutyléther) 328 2-Méthylpentane 334 Hexène 337 ETBE (Ethyltertiobutyléther) 346 Ethanol 351 Cyclohexane 354 n-Heptane 37_ 2,2,4-Triméthylpentane 372 Toluène 384 m-Xylène 412 1.2.4- Triméthylbenzène 443 [00141 Par ailleurs, la composition de l'air peut aussi être simplifiée et être réduite à deux espèces : di-oxygène (02), di-azote (N2). [00151 La composition chimique réduite ou détaillée ainsi obtenue à l'étape 2 ou 4 est notamment utilisée pour une étape 6 du procédé. L'étape 6 est une étape dans laquelle la 5 composition chimique de la phase gazeuse du carburant est calculée à tout instant. [0016 ] Le calcul est réalisé à l'aide d'un modèle. Le modèle est un modèle thermodynamique et plus précisément un modèle d'équilibre thermodynamique chimique. Le modèle thermodynamique permet à la fois de déterminer les compositions de la phase liquide et de la phase gazeuse du carburant à l'équilibre thermodynamique. 10 [00171 Le modèle thermodynamique repose sur l'utilisation d'une condition d'équilibre thermodynamique et utilise des données thermodynamique (température. pression...) ainsi que la composition chimique du mélange entre l'air et le carburant réalisé au niveau du conduit d'admission du moteur. 100181 En thermodynamique. deux phases (par exemple liquide. gaz ou solide) sont en 15 équilibre lorsque l'enthalpie libre notée G du système composé des deux phases atteint son minimum, soit lorsque la différentielle de l'enthalpie libre notée Ki est nulle. [0019 ] Or, à pression et température fixée, la différentielle de l'enthalpie libre s'écrit dans notre cas particulier d'un mélange binaire gaz-liquide : dCi = L [t; dn (ln (E1) 20 où n désigne le nombre de moles, l'indice i l'espèce considérée, l'indice gaz la phase gazeuse, l'indice liq désigne la phase liquide et u le potentiel chimique. [0020 ] De plus, lorsque le système est fermé et qu'aucune réaction chimique n'a lieu, la conservation de la matière implique que: dn + dn,= 0 (E2) 25 [0021 ] Les conditions de l'équation E2 peuvent être considérées comme remplies dans le cas du carburant contenu dans le réservoir à essence du véhicule. [0022 II peut être déduit des deux équations El et E_ que la condition d'équilibre thermodynamique du modèle thermodynamique devient : (u, :z,_ ù )dn,,,,_ = o(E3) [0023 ] De ce fait, le modèle thermodynamique peut se traduire par l'égalité des potentiels chimiques de chaque espèce i dans chacune des phases liquides ou gazeuses. [0024] Or. pour toute phase (gazeuse ou liquide), le potentiel chimique est défini par la relation Ea : ,u,(T,P)=,LL°(T,P °)+ RT In f , f ° où T est la température, P; est la pression de l'espèce considérée. R la constante des gaz parfaits, l'exposant �°� désigne les conditions standard soit une pression de 1 bar et une température de 0°C et enfin f; est la fugacité de l'espèce i dans une phase donnée. [0025] Ainsi, la condition d'équilibre thermodynamique du modèle thermodynamique se 10 traduit par l'égalité des fugacités de chaque espèce i dans chacune des phases. La condition d'équilibre entre les phases est donc pour chaque espèce i : ,.r,,, (E5) [0026 ] La fugacité s'exprimant par définition comme : f...t = q),, .P.X,..V (E6) 15 avec 9; a le coefficient de fugacité de l'espèce i dans la phase x (gaz pour gazeux, liq pour liquide). X, la fraction molaire de l'espèce i dans la phase x (gaz pour gazeux, liq pour liquide) et P la pression dans le conduit d'admission du moteur. [0027 1 Le modèle thermodynamique se traduit donc par l'égalité E7 pour chaque espèce i: rP- , où R la constante des gaz parfaits. [0031 ] Les équations d'état cubiques figurent parmi les quelques types d'équations d'état existantes. Fortement utilisées dans l'industrie, les équations cubiques sont assez simples et peuvent s'exprimer sous forme de polynômes du troisième degré en Z. Les équations cubiques ont l'avantage de rendre facultative la résolution des volumes qui est une méthode itérative nécessitant un temps de calcul non négligeable. Les deux équations les plus utilisées dans l'ingénierie sont l'équation de Peng et Robinson (notée équation PR clans la suite) et l'équation de Soave, Redlich et Kwong (notée équation SRK dans la suite) . [0032 ] Une équation d'état cubique peut s'exprimer en fonction du facteur de 10 compressibilité. sous la forme de l'équation E9 : Z' ù(I+B*ùkB*)I_-' +(.1*+k'B `2 ùkB*ùkB*')Zù(;1*13*+k'B+k'B*')=0 (E9) où k, k', A* et B* sont des paramètres détaillés dans les paragraphes qui suivent. [0033 ] Les constantes k et k' permettent de distinguer les différentes équations d'état. Par 15 exemple, l'équation PR correspond à des valeurs de k = 2 et k' = -1 tandis que l'équation SRK est caractérisée par les valeurs k = 1 et k' = O. [0034 ] Les paramètres A* et B* sont définis par les relations L,o : aP hP = et 13*=--- (Lin) R-'T2 RT où a et h sont respectivement le paramètre énergétique et le covolume. Plusieurs expressions de 20 a et b existent dans la littérature. [0035 ] Le calcul des paramètres a et b implique de connaître des propriétés physiques supplémentaires des espèces d'hydrocarbures et en particulier leur température critique notée leur pression critique notée P, et le facteur acentrique des corps pur noté w. [0036 ] Par exemple, il est possible d'utiliser la corrélation de Lee-Kesler pour calculer le 25 facteur acentrique w de chaque espèce de l'essence. En définissant le paramètre 0 par la relation O=Téh/T, où Tot, est la température d'ébullition de l'espèce et T, la température critique, la relation de Lee-Kesler s'écrit Ln P,=-(f,+L) où les paramètres fi ou f2 sont donnés par les relations Eä : 1; = 5,92714 6,09' _1,288621n O + 0.1693478' 30 f 15,6875 _ =152518 -13,47211n8+0,435778" (El1) [0037] Par exemple, sachant que pour le butane n-C411].:) , la pression critique vaut P,= 38.017 bar. la température critique T = 425,2 K et la température d'ébullition Teh = 272,7 K. le calcul du facteur acentrique par la méthode de Lee-Kesler donne e = 0.197. [0038 Les coefficients de fugacité de chaque composé calculés par la relation El: : Y Ltt(rp } = ~ (Z -1) ù Ln(Z h r) + 4 --g ' Ln bI3X~kJù4k' dans chacune des phases sont 2Z+ B'' k+~~k2 -4k' 2Z_ ùB* k+ti!k2 -4k' 7 avec (S; = ù~cr' . {l) où k1 désigne les coefficients d'interaction binaire dans cette /ai expression : a, et b, le paramètre énergétique et le covolume de l'espèce i, a,,, et b,,, du mélange [0039] Avec la connaissance des coefficients de fugacité. il devient alors possible de résoudre l'équation E7 et de déterminer la fraction de l'espèce i dans chacune des phases. [0040 ] La méthode par le calcul des coefficients de fugacité de chacune des espèces du carburant est relativement précise puisque cette méthode permet de prendre en compte le phénomène de dissolution des gaz dans la phase liquide. [0041 ] En contrepartie. le calcul des coefficients de fugacité nécessite beaucoup de temps de calculs et la connaissance de beaucoup de coefficients différents pour chaque espèce. Ainsi, pour chaque espèce i, le calcul du coefficient de fugacité implique la connaissance de la température critique, de la pression critique et si possible du coefficient binaire d'interaction avec chacune des autres espèces gazeuses. [0042] De ce fait, la méthode de calcul des coefficients de fugacité peut s'avérer peu commode à implémenter directement sur un ordinateur embarqué à bord d'un véhicule. Pour pallier à cet inconvénient, d'autres méthodes pour résoudre l'équation E7 peuvent être envisagées. [0043] En particulier, il peut être avantageux de calculer préalablement les coefficients de fugacité par la méthode de calcul précédente afin de les tabuler. Il suffit alors de résoudre E7 avec les coefficients tabulés. Le temps de calcul est ainsi réduit. Il est ainsi possible d'utiliser un calculateur moins puissant pour implémenter la méthode. La méthode par tabulation des coefficients de fugacité est donc plus facilement mise en place dans un ordinateur embarqué que la méthode par calcul. [00441 Un autre exemple de méthode plus simple à mettre en place que la méthode de calcul des coefficients de fugacité est une méthode où deux hypothèses sont préalablement réalisées. Les deux hypothèses sont l'hypothèse des gaz parfaits et l'hypothèse de mélange idéal dans la phase liquide. Ces hypothèses reviennent à considérer les coefficients de fugacité comme unitaires , ce qui permet de ne pas réaliser le calcul des coefficient de fugacité pour chaque espèce du carburant. [0045 L'utilisation de la loi de Raoult permet alors d'obtenir des fractions molaires en phase gazeuse à l'état de saturation. par la relation E X =P (Np) J, où P,,,uest la pression de vapeur saturante et P est la pression du système. [0046] La pression de vapeur saturante peut être calculée par la relation d'Antoine : où A,. B; et C sont les coefficients d'Antoine et Ts,,,,icc est la température de la surface de l'interface entre la phase gaz et la phase liquide_ Les coefficients A;, B, et C, sont spécifiques de l'hydrocarbure i considéré et s'obtiennent par des corrélations. [00171 La méthode est valable pour des mélanges d'hydrocarbures. Cette méthode basée sur l'hypothèse des gaz parfaits et la loi de Raoult est plus simple que la méthode de calcul des coefficients de fugacité. De ce fait, la méthode est particulièrement adaptée à une utilisation dans des ordinateurs embarqués à bord de voiture. [00481 ytais, la méthode n'est pas applicable dans tous les cas. Ln particulier, la méthode n'est plus correcte dès lors que les différentes molécules du mélange ont de fortes interactions entre elles. Un mélange comprenant des espèces polaires comme les alcools est un exemple de mélange auquel la méthode basée sur l'hypothèse des gaz parfaits et la loi de Raoult ne peut pas s'appliquer. [0049 Quelle que soit la méthode choisie pour mener le calcul de l'étape 6, l'étape 6 permet d'obtenir la composition chimique de la phase liquide à l'interface de la phase gaz et la phase liquide à l'étape 8 ainsi que la composition chimique de la phase gaz à l'interface de la phase gaz et la phase liquide à l'étape 10. L'information de la composition chimique n'est pas un paramètre physico-chimique trop réducteur comme la pression de vapeur. La connaissance de la composition chimique dans les deux phases permet en effet de connaître les caractéristiques essentielles du système. [0050 A titre d'exemple, le tableau suivant présente les compositions chimiques des phases gazeuse et liquide obtenues aux étapes 8 et 10 du procédé pour une essence commerciale Sans Plomb 95 européenne actuelle, à une pression de 1 atm et une température de 20°C. I.,e calcul de l'étape 6 est effectué par la méthode de calcul des coefficients de fugacité. Espèce chimique Phase liquide Phase gazeuse (fraction molaire) (fraction molaire) n-Butane 0. 078 0.160 Isopentane 0.136 0.108 n-Pentane 0.040 0.023 Méthvlbutène 0.095 0.049 MTBE (méthyltertiobutv léther) 0.054 0.015 2-Méthylpentane 0.060 0.014 Hexène 0.043 0.009 ETBE (éthyltertiobutvléther) 0.000 0.000 Ethanol 0.000 0.000 Cyclohexane 0.025 0.003 n-Heptane 0.063 0.003 2,2,4-Triméthylpentane 0.073 0.004 Toluène 0.101 0.004 in-Xylène 0.109 0.001 1.2,4-Triméthylbenzène 0.119 0.000 D'azote ---~ 0.001 0.472 Dioxygène 0.000 0.133 [0051 1 Le tableau met en évidence que certaines espèces ont commencé à s'évaporer comme le n-Butane ou le n-Pentane alors que d'autre espèces comme le m-Xylène ou le 1.2.4-Triméthylbenzè.ne ont une fraction molaire en phase gazeuse nulle ou quasi-nulle. [0052 ] La composition chimique des phases gazeuse et liquide à l'interface entre la phase gaz et la phase liquide correspond à la composition chimique au niveau du conduit d'admission du moteur. Il ne s'agit néanmoins pas de la composition chimique du mélange entre le carburant et l'air qui est brûlé au niveau de la bougie de la chambre de combustion. En effet. le passage de la phase gazeuse de l'interface entre la phase gaz et la phase liquide vers le coeur du conduit d'admission et vers le coeur de la chambre de combustion modifie la composition chimique du mélange parce que les espèces n'ont pas la même vitesse de déplacement dans la phase gazeuse. Il en résulte que la quantité d'espèces présentes à l'interface entre la phase gaz et la phase liquide n'est pas la même dans la chambre de combustion. [00531 De ce fait, le procédé comprend en outre une étape 12 dans laquelle le transport en phase gaz est étudié. Cette analyse consiste à calculer le taux de transport de la phase gaz à l'aide de la composition chimique de la phase gaz à l'interface entre la phase gaz et la phase liquide. [00541 Le calcul est réalisé par utilisation d'un modèle. Le modèle est un modèle de transport qui est valable pour les moteurs à injection indirecte c'est-à-dire des moteurs ou le mélange carburant / air se fait avant l'admission dans la chambre de combustion. Cela signifie notamment que le mélange carburant `air est réalisé avant une soupape d'admission d'un des cylindres du moteur. Le modèle de transport s'appuie sur plusieurs hypothèses successivement exposées dans les paragraphes qui suivent. [00551 On suppose d'abord qu'au démarrage, le carburant se dépose sous la forme d'un film liquide sur les parois froides du conduit d'admission du moteur. On admet alors que toute l'évaporation provient de l'évaporation du film liquide. I 'évaporation dans le cylindre en fin de compression ainsi que l'évaporation liée à l'arrachage des gouttes du film liquide pendant l'ouverture de la soupape d'admission sont donc négligées. [0056 1 De plus, il est admis que l'évaporation a lieu pendant la phase d'admission (durant un aller du piston soit 180°V). c'est-à-dire le temps pendant lequel la soupape est ouverte. L'évaporation avant l'ouverture de la soupape, plus précisément entre l'injection et l'ouverture de la soupape, n'est pas prise en compte. [0057] Du fait de la très grande différence de capacité calorifique entre le métal du conduit d'admission et le carburant liquide, la température du film liquide est supposée être la température du conduit d'admission lors du premier cycle du moteur. [0058 1 En outre, on admet, conformément à une hypothèse classique dans la littérature, que la hauteur du film liquide notée hf;i,,, est proportionnelle à la masse injectée notée mmj de carburant à une certaine puissance notée /6 : h,;;,,, °` (m,,,Y [0059 ] La même relation peut être écrite entre la masse injectée et la surface du film liquide notée °c )~ù(3 `~Jil,n n1r;g [0060 1 Il est courant dans la littérature d'attribuer à ce coefficient fi une valeur proche de 0.5. Cependant, il ne s'agit que d'un exemple puisque le coefficientf dépend de diverses caractéristiques du moteur et notamment de la forme du conduit d'admission. de l'injecteur... Ce paramètre peut ainsi servir à caler le modèle pour chaque moteur. [0061 ] On fait aussi les hypothèses que la hauteur et la surface du film liquide ne dépendent pas de la température et qu'il existe une surface mouillée maximum pour chaque 5 conduit d'admission. - [0062] Enfin, on suppose que la diffusion des espèces est infiniment rapide dans le film liquide. En d'autres termes. les fractions molaires des espèces sont homogènes dans le film liquide. [0063 ] Le modèle de transport modélise ensuite le transport de l'espèce i dans la 10 phase gazeuse par une relation issue de la théorie de la mécanique des fluides. En particulier. il est possible d'utiliser la relation Ei issue de la littérature : =a*(0,552*Re"_,) �' )'-/1 *D... *ln (Ei) z z où m,,,,,, , a. Re. DIA et Y;0 sont des paramètres détaillés dans les paragraphes suivants. [00641 Le paramètre m:,. est la dérivée temporelle de la masse de l'espèce i et représente 15 le taux de transport de l'espèce i dans la phase gazeuse, depuis l'interface entre la phase gaz et la phase liquide vers le coeur du mélange gazeux. [0065j Le coefficient a est un coefficient à ajuster pour caler le modèle sur le moteur considéré. [0066] La constante Re est le nombre de Reynolds. La constante Re se calcule avec la 20 vitesse moyenne dans le conduit d'admission. [0067 ] v est le coefficient de diffusion binaire entre le carburant et l'air. Il est possible de considérer un coefficient global Di A pour le carburant ou un coefficient Dl.,:v par espèce. [00681 Enfin, Y, (,A/ est la fraction massique de l'espèce i en phase gazeuse à l'interface entre la phase gaz et la phase liquide. La fraction massique de l'espèce i est liée à la fraction 25 molaire calculée précédemment par application du modèle thermodynamique de l'étape 6. [0069 1 De ce fait. compte tenu des hypothèses précédemment décrites. l'étape 12 permet de déterminer à partir de la masse de carburant injectée le taux de transport m;,.,,� de chaque espèce i dans la phase gazeuse, depuis l'interface entre la phase gaz et la phase liquide vers le coeur du mélange gazeux. [0070 ] Une étape 14 supplémentaire du procédé peut en outre permettre de déduire la richesse du mélange air/carburant obtenu à l'aide d'un calcul. [0071 ] La richesse du mélange carburé dans le cylindre est obtenue par le rapport entre la valeur de l'intégrale de la somme des taux de transport clans la phase gaz sur 180°V et la valeur du volume d'air introduit. [0072 ] Le mélange carburé dans le cylindre étant considéré comme homogène. la richesse obtenue à l'étape 14 correspond à la richesse au point d'allumage. Le procédé permet ainsi de déterminer la richesse au point d'allumage à partir de conditions de pression, de température donnés et d'une masse de carburant injecté donné. [0073 ] Le procédé peut en outre comprendre une étape 16 de détermination de la quantité de carburant à injecter. [0074 ] Pour réaliser la détermination. la richesse obtenue à l'étape 14 est comparée à une valeur prédéterminée de richesse à atteindre. La valeur de la richesse à atteindre pour assurer un bon démarrage du véhicule dépend notamment du moteur considéré et de la puissance de la bougie. Par exemple. la littérature s'accorde usuellement autour d'une richesse de 0,75 à 0.8 pour permettre la première combustion. [0075 ] La quantité de carburant à injecter au démarrage peut être obtenue par itération successive sur la masse de carburant à injecter. Une masse initiale d'injection est arbitrairement choisie. Si cette masse initiale conduit à une richesse supérieure à la valeur prédéterminée de richesse à atteindre, cela signifie que la masse à injecter est trop importante. Le calcul de l'étape 12 est recommencé avec une masse à injecter plus petite. Au contraire, si cette niasse initiale correspond à une richesse inférieure à la valeur prédéterminée de richesse à atteindre. le calcul de l'étape 12 est repris avec la masse à injecter augmentée. 1 e processus est réitéré jusqu'à qu'une richesse satisfaisante soit obtenue. [0076 ] La quantité de carburant à injecter peut aussi être obtenue par utilisation de tables donnant la quantité de carburant liquide à utiliser au démarrage. Une manière commode d'implémenter cette méthode informatiquement peut être d'utiliser une base de données clans un ordinateur de bord du véhicule. [0077 ] La quantité de carburant minimale à injecter est ainsi connue à l'étape 18 du procédé. Ceci permet de réduire les émissions polluantes, notamment les émissions de CO2. De ce fait, la charge en métaux précieux du catalyseur placé en post-traitement de la ligne d'échappement peut être diminué. Le coût du post-traitement catalytique est donc réduit. Le prix de revient de fabrication de la voiture est ainsi diminué. [007$ ] La qualité et la robustesse de la prestation démarrage mise en action est aussi améliorée. De plus, cette prestation ne dépend plus du carburant. que ce soit d'un pa's à l'autre ou d'une saison à l'autre. [0079] A bord d'un véhicule, le procédé peut être utilisé dans un organe de contrôle du moteur. L'information � composition détaillée du carburant � est obtenue par un capteur embarqué, par exemple un capteur utilisant la spectroscopie proche infrarouge. immergé clans le réservoir carburant du véhicule. [0050 ] La donnée de sortie de l'outil de modélisation (quantité de carburant liquide à injecter au démarrage) est une donnée d'entrée des lois de commande du contrôle moteur. [0081 ] La qualité de la prestation de démarrage est améliorée. De plus, les défaillances et les coûts de garantie sont diminués. Tous ces éléments contribuent à augmenter la satisfaction de l'utilisateur. [0052 ] Le procédé peut aussi être utilisé dans un outil d'aide à la mise au point d'un moteur, notamment dans un laboratoire. Le temps de mise au point des modèles est ainsi diminué. Le nombre d'essais nécessaires à la mise au point du démarrage peut être diminué. De plus. cela permet de calibrer la masse initiale de carburant utilisée à l'étape 16 du procédé. [005 5 ] Le procédé de détermination de la quantité de carburant à injecter au démarrage d'un moteur à injection indirecte selon l'invention permettent donc de déterminer la quantité de carburant nécessaire au démarrage d'un moteur à allumage commandé injection indirecte. REVENDICATIONS I. Procédé de détermination de la quantité de carburant à injecter au démarrage d'un moteur à injection indirecte, un mélange d'une phase gaz d'air et de carburant et d'une phase 5 liquide de carburant étant injecté dans le moteur le procédé comprenant : - le calcul de la composition chimique de la phase gaz et de la phase liquide à l'interface entre la phase gaz et la phase liquide. par utilisation d'un modèle, - à l'aide de la composition chimique de la phase gaz, le calcul du taux de transport de la phase 10 gaz, par utilisation d'un modèle. 2. Procédé selon la revendication I. le moteur comprenant une chambre de combustion et un conduit d'admission dans le moteur du mélange, la phase liquide étant déposée sur les parois du conduit d'admission. 3. Procédé selon la revendication 1 ou 2, dans lequel le modèle de calcul de la composition 15 chimique de la phase gaz et de la phase liquide prend en considération la composition chimique du carburant préalablement déterminée. 4. Procédé selon la revendication 3, dans lequel le nombre d'espèces du carburant pour la détermination de la composition chimique du carburant est réduit. de préférence à quinze. 5. Procédé selon la revendication 3 ou 4, dans lequel la composition chimique du carburant 20 est déterminée par un organe d'analyse de laboratoire ou par un organe d'analyse embarqué à bord d'un véhicule. 6. Procédé selon la revendication 5, dans lequel l'organe d'analyse embarqué est un capteur opto-électronique utilisant la spectroscopie proche infrarouge. 7. Procédé selon l'une des revendications 1 à 6, dans lequel le modèle de calcul de la 25 composition chimique de la phase gaz et de la phase liquide est basé sur l'hypothèse des gaz parfaits et la loi de Raoult. g. Procédé selon l'une des revendications 1 à 6. dans lequel le modèle de calcul de la composition chimique de la phase gaz et de la phase liquide prend en considération lecoefficient de fugacité de chacune des espèces du carburant, le coefficient de fugacité étant calculé ou obtenu dans des tables. 9. Procédé selon l'une des revendications 1 à 4. comprenant en outre le calcul de la richesse du mélange obtenue par le rapport entre la valeur de l'intégrale du taux de transport de 5 la phase gaz et la valeur du volume d'air injecté. 10. Procédé selon la revendication 9, le modèle de calcul du taux de transport prend en considération la masse de carburant à injecter, le procédé comprenant en outre la détermination de la quantité de carburant à injecter au démarrage dans le moteur pour parvenir à une richesse prédéterminée du mélange par utilisation de tables ou par itération de la masse de carburant à 10 injecter. 11. Utilisation du procédé selon l'une des revendications 1 à 10 dans un organe de contrôle d'un moteur ou dans un outil d'aide à la mise au point d'un moteur.