La présente invention concerne le domaine de la détermination de la masse de gaz enfermée dans une chambre de combustion d'un moteur à combustion interne.

La masse de gaz enfermée dans une chambre de combustion est un paramètre utile notamment pour le contrôle des moteurs à combustion interne, pour les tests des moteurs à combustion interne, afin d'analyser notamment la combustion, d'améliorer les stratégies de contrôle de la combustion, d'améliorer le comportement cycle à cycle et entre cylindres d'un moteur à combustion interne. La connaissance de ce paramètre est important en particulier pour les moteurs à allumage commandé et pour les moteurs Diesel équipés d'une recirculation des gaz d'échappement (EGR) et/ou d'une distribution variable.

Aucun capteur ne permet de mesurer directement la masse de gaz enfermée dans une chambre de combustion. C'est pourquoi, plusieurs méthodes ont été développées pour estimer ce paramètre. Toutefois, aucune méthode ne permet de disposer d'une information fiable.

En particulier, certaines méthodes sont basées sur des mesures de la pression d'admission, qui peut ne pas représenter de manière fiable les phénomènes physiques mis en jeu au sein de la chambre de combustion. D'autres méthodes nécessitent des techniques de traitement de signal qui peuvent être complexes, longues et qui peuvent également générer des imprécisions sur la masse de gaz enfermée estimée.

La demande de brevet WO 2007060349 décrit une méthode basée sur l'utilisation d'un capteur de pression dans le cylindre et d'un capteur de température en aval de la soupape d'échappement. Ainsi, cette méthode nécessite une instrumentation spécifique du moteur à combustion interne. En outre, cette méthode reconstruit la température du gaz dans le cylindre à partir de la mesure du capteur de température en aval de la soupape d'échappement. Cette reconstruction génère des approximations, qui influent sur la précision de la masse de gaz enfermée estimée.

La demande de brevet WO 2015082731 décrit une méthode basée sur l'estimation de la fréquence de résonance de la pression dans la chambre de combustion au moyen d'une transformée de Fourier, qui est un traitement du signal issu du capteur de pression. Cette estimation de la fréquence de résonance manque de précision pour une estimation précise de la masse de gaz enfermée. En outre, cette méthode utilise la résolution des équations de Bessel, ayant pour hypothèse que la chambre de combustion est un cylindre parfait. Cependant, les chambres de combustion des moteurs à combustion interne ne sont pas purement cylindriques. L'estimation obtenue par cette méthode se retrouve par conséquent erronée en raison de cette hypothèse. De plus, cette méthode nécessite des oscillations d'amplitudes importantes. La thèse de Pau Bares "In-cylinder pressure resonance analysis for trapped mass estimation in automotive engines" (30 juin 2017 (2017-06-30), Valencia) divulgue une calculation de la masse d'air et du combustible dans la chambre de combustion.

La présente invention a pour but de déterminer de manière précise la masse de gaz enfermée dans une chambre de combustion d'un moteur à combustion interne, avec une instrumentation simple et classique, et sans traitement du signal complexe. Dans ce but, le procédé de détermination de la masse de gaz enfermée dans une chambre de combustion selon l'invention se base sur la seule mesure de la pression au sein de la chambre de combustion. Le procédé met en œuvre un modèle de la chambre de combustion, dont un des paramètres, la fréquence d'oscillations de pression, est déterminée avec la mesure de la pression de gaz dans la chambre de combustion. Ainsi, le procédé selon l'invention est précis, fiable, et nécessite une instrumentation simple.

L'invention selon la revendication 1 concerne un procédé de détermination de la masse de gaz enfermée dans une chambre de combustion d'un moteur à combustion interne, dans lequel on met en œuvre les étapes suivantes :

1. a) On mesure la pression P dudit gaz dans ladite chambre de combustion ;
2. b) On détermine les oscillations de ladite pression de gaz P à partir de ladite mesure de pression P dudit gaz ;
3. c) On détermine ladite masse de gaz enfermée m au moyen de la formule $m=\frac{\mathit{\gamma PV}}{{\lambda }^2{f}^2}$ et par détermination de la fréquence d'oscillations f de pression avec lesdites oscillations mesurées, V étant le volume de ladite chambre de combustion, λ la longueur d'onde dans ladite chambre de combustion, P ladite pression mesurée dudit gaz dans ladite chambre de combustion, et γ le rapport thermique spécifique.

Selon un mode de réalisation, ledit volume V de ladite chambre de combustion est obtenu au moyen d'une cartographie en fonction de l'angle vilebrequin dudit moteur à combustion interne.

Selon l'invention, ladite longueur d'onde λ est obtenue au moyen d'une cartographie en fonction de l'angle vilebrequin dudit moteur à combustion interne.

Selon l'invention, ladite cartographie de ladite longueur d'onde λ est construite préalablement numériquement.

Selon l'invention, ladite cartographie de ladite longueur d'onde λ est construite par une approche par éléments finis en trois dimensions.

Selon un aspect de l'invention, on détermine ledit rapport thermique spécifique γ en prenant en compte la composition de ladite masse enfermée, et au moyen de la capacité thermique des composants de ladite masse enfermée.

Conformément à une caractéristique, on détermine ladite masse de gaz enfermée en temps réel.

Selon l'invention, on détermine la fréquence d'oscillations de pression f en mettant en œuvre les étapes suivantes de calibration :

1. i) On considère au moins deux masses de gaz enfermées prédéfinies m₁, ..., m_n ;
2. ii) On détermine un rapport thermique spécifique γ pour chaque masse de gaz enfermée prédéfinie m₁, ..., m_n ;
3. iii) On détermine au moyen de ladite formule $m=\frac{\mathit{\gamma PV}}{{\lambda }^2{f}^2}$ une fréquence d'oscillations de pression f pour chaque masse de gaz enfermée prédéfinie m₁, ..., m_n en fonction de ladite pression de gaz P mesurée, dudit volume V de la chambre de combustion, de ladite longueur d'onde λ et dudit rapport thermique spécifique γ déterminé ; et
4. iv) On compare chaque fréquence d'oscillations de pression f déterminée avec lesdites oscillations de pression mesurées, et on détermine ladite masse de gaz enfermée en tant que la masse prédéfinie pour laquelle la fréquence d'oscillations de pression f minimise ladite comparaison.

De manière avantageuse, on détermine chaque rapport thermique spécifique γ de chaque masse de gaz enfermée prédéfinie m₁, ..., m_n en mettant en œuvre les étapes suivantes :

1. (1) On détermine la composition de ladite masse enfermée ;
2. (2) On détermine une relation qui relie ledit rapport thermique spécifique γ et la capacité du gaz c_p dans ladite chambre de combustion : $\gamma =\frac{C_p}{C_p-R}$ avec R la constante du gaz, et la capacité du gaz c_p étant fonction de la température dudit gaz et de ladite composition de ladite masse enfermée ;
3. (3) On détermine la température dudit gaz au moyen de la loi des gaz parfaits appliquée à chaque masse enfermée prédéfinie m₁, ..., m_n ; et
4. (4) On en déduit ledit rapport thermique spécifique γ pour chaque masse enfermée prédéfinie m₁, ..., m_n.

Alternativement, on détermine la fréquence d'oscillations de pression f au moyen d'une moyenne sur une fenêtre temporelle desdites oscillations déterminées.

Selon une mise en œuvre, on détermine une longueur d'onde λ pour chaque mode d'oscillations de ladite pression dudit gaz.

D'autres caractéristiques et avantages du procédé 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.

• La figure 1 illustre un cylindre d'un moteur à combustion interne.
• La figure 2 illustre les étapes du procédé selon un mode de réalisation de l'invention.
• La figure 3 illustre les étapes de la calibration de la fréquence d'oscillations de pression selon un mode de réalisation de l'invention.
• La figure 4 illustre un graphique de la masse relative de gaz estimée par le procédé selon l'invention pour un exemple.

La présente invention concerne un procédé de détermination de la masse de gaz enfermée dans une chambre de combustion d'un moteur à combustion interne.

De manière classique, un moteur à combustion interne comprend au moins un cylindre, un piston coulissant dans ce cylindre en un mouvement rectiligne alternatif, des moyens d'admission d'un comburant (gaz), des moyens d'échappement de gaz brûlés, une chambre de combustion, et des moyens d'injection pour injecter un combustible dans la chambre de combustion. La chambre de combustion est agencée dans la partie haute de cylindre, et son volume varie par le déplacement du piston dans le cylindre. La figure 1 illustre, schématiquement et de manière non limitative, un cylindre 1 d'un moteur à combustion interne selon un mode de réalisation de l'invention. Sur cette figure, les moyens d'admission, d'échappement, d'injection, et les éventuels moyens d'allumage ne sont pas représentés. Au sein du cylindre 1 se déplace un piston 6 selon un mouvement rectiligne alternatif. La chambre de combustion 2 est la zone où se produit la combustion, elle est limitée par la partie supérieure du piston 6, la paroi latérale 4 du cylindre, et le toit du cylindre 5. Cette zone correspond à la zone blanche de la figure 1. En outre, pour le procédé selon l'invention, le cylindre comporte un capteur de pression 3, qui mesure la pression du gaz dans la chambre de combustion 2.

Le moteur à combustion interne peut être de tout type, à allumage commandé ou à auto-inflammation, avec ou sans recirculation des gaz d'échappement, avec ou sans suralimentation. Toutefois, le procédé selon l'invention est particulièrement adapté aux moteurs à allumage commandé, et aux moteurs à auto-inflammation équipés d'une recirculation des gaz d'échappement (EGR) et/ou d'une distribution variable.

Par les termes gaz ou comburant, il est compris de l'air à pression ambiante ou de l'air suralimenté ou encore un mélange d'air (suralimenté ou non) avec des gaz brûlés.

Le procédé de détermination de la masse de gaz enfermée selon l'invention utilise un unique capteur : un capteur de pression dans la chambre de combustion. L'utilisation d'un unique capteur permet de limiter l'instrumentation du moteur. De plus, la mesure de la pression dans la chambre de combustion permet d'avoir directement des informations relatives au comportement du gaz dans la chambre de combustion, contrairement à un capteur placé à l'admission ou à l'échappement qui nécessiterait de reconstruire la pression dans la chambre de combustion. Un tel capteur présente également l'avantage d'équiper classiquement un banc moteur, et les moteurs à combustion interne des véhicules récents, en particulier pour les nouveaux moteurs à auto-inflammation, ce qui facilite la mise en œuvre du procédé, et ce qui limite son coût.

Le procédé selon l'invention comporte les étapes suivantes :

• Mesure de la pression du gaz dans la chambre de combustion ;
• Détermination des oscillations de la pression de gaz mesurée ; et
• Détermination de la masse de gaz enfermée par détermination de la fréquence d'oscillations avec les mesures.

Ces étapes seront détaillés dans la suite de la description.

Selon un mode de réalisation de l'invention, ces étapes peuvent être mises en œuvre en temps réel, afin de déterminer en temps réel la masse de gaz enfermée, facilitant l'exploitation de ce paramètre.

De plus, ces étapes peuvent être mises en œuvre en ligne, sur un moteur à combustion interne d'un véhicule ou sur un banc moteur.

Les étapes de détermination des oscillations de la pression de gaz et de la détermination de la masse de gaz enfermée peuvent être mises en œuvre par des moyens informatiques, notamment un calculateur équipant le moteur à combustion interne.

La figure 2 illustre, schématiquement et de manière non limitative, les étapes du procédé selon l'invention. La première étape consiste à mesurer la pression du gaz MES P dans la chambre de combustion, au moyen d'un capteur de pression. La deuxième étape consiste à déterminer les oscillations de pression de gaz OSC P à partir de la mesure de pression de gaz réalisée précédemment. Ensuite, on détermine la masse de gaz enfermée m par détermination DET f de la fréquence d'oscillations f de pression avec les oscillations de pression déterminées à l'étape précédente.

Lors de cette étape, on mesure la pression du gaz dans la chambre de combustion, au moyen d'un capteur de gaz placé dans la chambre de combustion.

Lors de cette étape, on détermine les oscillations de la pression de gaz, au cours du temps, à partir de la pression du gaz. Ainsi, on peut connaître précisément les phénomènes physiques mis en œuvre dans la chambre de combustion.

Selon un mode de réalisation, cette étape peut être mise en œuvre au moyen d'un filtre, de préférence au moyen d'un filtre passe bande autour de la fréquence d'intérêt. Selon un exemple non limitatif, on peut mettre en œuvre un filtre passe bande entre 5 et 9 kHz. De manière non limitative, on peut utiliser un filtre Butterworth.

Lors de cette étape, on détermine la masse de gaz enfermée m au moyen de la formule $m=\frac{\mathit{\gamma PV}}{{\lambda }^2{f}^2}$ (appelée également modèle de la chambre de combustion) et par détermination de la fréquence d'oscillations f de pression avec les oscillations mesurées déterminées à l'étape précédente. Dans cette formule, V est le volume de la chambre de combustion (variable dans le temps), λ est la longueur d'onde dans la chambre de combustion, P la pression mesurée du gaz dans la chambre de combustion, et γ le rapport thermique spécifique.

La formule $m=\frac{\mathit{\gamma PV}}{{\lambda }^2{f}^2}$ est obtenue à partir de la loi des gaz parfaits : PV = mRT, de l'équation de la fréquence d'oscillations $f=\frac{c}{\lambda }$ et de l'équation de la vitesse du son $c=\sqrt{\mathit{\gamma RT}}$, c étant la vitesse instantanée du son dans la chambre de combustion. Ainsi, la formule est représentative des phénomènes physiques du gaz dans la chambre de combustion, ce qui permet l'obtention précise de la masse de gaz enfermée.

Selon un mode de réalisation de l'invention, le volume V de la chambre de combustion peut être obtenu au moyen d'une cartographie ou une table du volume de la chambre de combustion en fonction de l'angle vilebrequin du moteur à combustion interne. Cette cartographie traduit les variations du volume lié au mouvement du piston. L'utilisation d'une telle cartographie permet de rendre la détermination du volume de la chambre de combustion rapide, et utilisable en temps réel.

La longueur d'onde λ dépend de la position du piston (elle est donc variable dans le temps) et de la forme de la chambre de combustion. Conformément à une mise en œuvre de l'invention, la longueur d'onde λ peut être obtenue au moyen d'une cartographie de la longueur d'onde λ en fonction de l'angle vilebrequin du moteur à combustion interne. L'utilisation d'une telle cartographie permet de rendre la détermination de la longueur d'onde λ rapide, et utilisable en temps réel.

Pour cette mise en œuvre, la cartographie de la longueur d'onde λ peut être construite préalablement et numériquement. La détermination numérique de la longueur d'onde λ, permet de déterminer une longueur d'onde λ pour toute forme de la chambre de combustion.

Ainsi, la détermination de la masse enfermée ne comporte pas d'approximation liée à la forme de la chambre de combustion.

Selon un aspect préféré de cette mise en œuvre, la cartographie de la longueur d'onde λ peut être construite par une approche par éléments finis en trois dimensions (appliquée à la chambre de combustion) pour résoudre l'équation d'onde en trois dimensions à chaque angle vilebrequin, et pour ainsi déterminer le rapport entre la fréquence d'oscillations de pression et la vitesse imposée du son. Les entrées de l'approche par éléments finis peuvent être :

• La géométrie de la chambre de combustion à chaque angle vilebrequin,
• L'équation d'onde à résoudre : $\mathrm{\Delta }P=\frac{1}{c^2}\frac{{\delta }^{2}P}{{\mathit{\delta t}}^2}$ avec la vitesse du son c fixe,
• Les conditions aux limites : ∇P. n = 0, n étant le vecteur normal à la paroi du cylindre,
• Les conditions initiales : source, c'est-à-dire un gradient de pression localisé utilisé pour initialiser le calcul.

La résolution de l'équation d'onde permet de surveiller les oscillations de pression à différents endroits de la chambre de combustion ; par la suite, les fréquences d'oscillations peuvent être dérivées via une simple analyse de densité spectrale et les longueurs d'onde associées à chaque mode d'oscillation peuvent être calculées, pour chaque angle vilebrequin.

Alternativement, la longueur d'onde λ peut être déterminée par d'autres méthodes, par exemple analytiquement.

De préférence, on peut déterminer une longueur d'onde λ pour chaque mode d'oscillations de la pression de gaz, de manière à obtenir une information plus précise.

Le rapport thermique spécifique y dépend du temps et varie d'un point de fonctionnement à l'autre du moteur à combustion interne.

Selon un mode de réalisation de l'invention, on peut déterminer le rapport thermique spécifique γ en prenant en compte la composition de la masse enfermée dans la chambre de combustion, et au moyen de la capacité thermique des composants de la masse enfermée dans la chambre de combustion.

Pour ce mode de réalisation, on peut prendre en compte la fraction volumique des composants du carburant injecté dans la chambre de combustion (par exemple la fraction volumique d'éthanol dans une essence), le rapport quantité d'air par quantité de carburant dans la chambre de combustion par rapport au rapport stœchiométrique, et la proportion de gaz imbrûlés dans la chambre de combustion. De plus, ce mode de réalisation peut mettre en œuvre les étapes suivantes :

• On détermine les coefficients de la formule chimique du carburant, par exemple les coefficients x, y, z pour un carburant de la forme CxHyOz (lorsque le carburant n'est pas un mélange, et lorsque le carburant est un mélange on peut pondérer les coefficients des composants de base du carburant, en fonction de la fraction volumique des composants du carburant) ;
• On détermine la capacité thermique du mélange non brûlé dans la chambre de combustion, en prenant en compte les coefficients de la formule chimique du carburant (par exemple x, y, z), et le rapport quantité d'air par quantité de carburant dans la chambre de combustion par rapport au rapport stœchiométrique. On peut se servir de l'équation bilan de la combustion (pour déterminer la concentration des composants de l'air et du carburant) et d'une équation de capacité Cp thermique de chaque composant de l'équation bilan, cette équation pouvant être de la forme Cp(T) = a + b × T + c × T² + d × T³ + e × T⁴ avec T la température dans la chambre de combustion, et a, b, c, d et e des coefficients obtenus par des tables (ex : Table de JANAF) pour chaque composant du mélange non brûlé,
• On détermine la capacité thermique du mélange brûlé dans la chambre de combustion, en prenant en compte les coefficients de la formule chimique du carburant (par exemple x, y, z), la méthodologie pouvant être identique à celle décrite pour la capacité thermique du mélange non brûlé. On peut se servir de l'équation bilan de la combustion (pour déterminer la concentration des composants du mélange brûlé) et des équations de capacité Cp thermique de chaque composant.
• On détermine la capacité thermique du mélange carburé par association de la capacité thermique du mélange brûlé et de la capacité thermique du mélange non brûlé, selon un mode de réalisation de cette étape, on peut, par exemple, considérer que le mélange carburé comprend 50% de mélange brûlé et 50% de mélange non brûlé (ou toute autre répartition),
• On détermine le rapport de capacité thermique γ par l'équation $\gamma =\frac{c_p}{C_p-R}$ avec c_p la capacité thermique du mélange carburé déterminée à l'étape précédente, et R la constante des gaz parfaits.

Pour la mise en œuvre de ce mode de réalisation, la température T du mélange peut être déduite de la loi des gaz parfaits PV = mRT en sélectionnant au moins une masse de gaz enfermée hypothétique, P, V et R étant connus par ailleurs.

Selon un deuxième mode de réalisation, le rapport de capacité thermique peut être obtenu par une formule du type $\gamma =\frac{-\mathit{PdV}}{\mathit{VdP}}$ sur la base d'une hypothèse d'un comportement adiabatique.

Conformément à une mise en œuvre de l'invention, on peut déterminer la fréquence d'oscillations de pression en mettant en œuvre une calibration au moyen des étapes suivantes:

• On considère au moins deux masses de gaz enfermées prédéfinies m₁, ..., m_n, il s'agit d'au moins deux masses de gaz enfermées possibles, n est un entier supérieur ou égal à 2, selon un aspect de l'invention, les masses de gaz enfermées prédéfinies m₁, ..., m_n sont contenues dans une liste prédéfinie préalablement au déroulé du procédé,
• On détermine un rapport thermique spécifique γ pour chaque masse de gaz enfermée prédéfinie m₁, ..., m_n , par exemple on peut déduire de chaque masse de gaz enfermée prédéfinie m₁, ..., m_n une température T dans la chambre de combustion, et en appliquant les étapes décrites ci-dessus, on peut en déduire la capacité thermique du mélange c_p, qui est ensuite utilisée pour obtenir le rapport thermique spécifique γ en fonction de l'angle vilebrequin,
• On détermine au moyen de la formule $m=\frac{\mathit{\gamma PV}}{{\lambda }^2{f}^2}$ une fréquence d'oscillations de pression f pour chaque masse de gaz enfermée prédéfinie m₁, ..., m_n en fonction de la pression mesurée P et des paramètres déterminés, le volume V de la chambre de combustion et la longueur d'onde λ, et le rapport thermique spécifique γ déterminé à l'étape précédente, en d'autres termes on inverse la formule $m=\frac{\mathit{\gamma PV}}{{\lambda }^2{f}^2}$ pour obtenir une équation $f=\sqrt{\frac{\mathit{\gamma PV}}{{\mathit{m\lambda }}^2}}$ dont tous les paramètres sont déterminés (P, V, y, λ) ou prédéfini (m),
• On compare chaque fréquence d'oscillations de pression f obtenue à l'étape précédente avec les oscillations de pression mesurée (déterminées à l'étape 2), et on détermine la masse de gaz enfermée en tant qu'une des masses de gaz enfermées prédéfinies m₁, ..., m_n , pour laquelle la fréquence d'oscillations de pression f minimise la comparaison, en d'autres termes, la masse de gaz enfermée déterminée par le procédé est une des masses de gaz enfermées prédéfinies m₁, ..., m_n : celle pour laquelle la différence entre la fréquence d'oscillation de pression f et les oscillations de pression mesurée est minimale. La comparaison peut consister en une comparaison des oscillations. De préférence, seule la fréquence des oscillations de pression mesurée est utilisée pour la calibration, l'amplitude de celles-ci n'est pas utile pour la comparaison. De préférence, cette étape peut consister en une génération des oscillations de pression avec les différentes fréquences déterminées puis on compare les oscillations obtenues avec celles déterminées à l'étape 2, et on choisit le signal pour lequel l'oscillation du signal est comparable aux oscillations de la pression mesurée.

La figure 3 illustre, schématiquement et de manière non limitative, cette mise en œuvre de cette étape de l'invention. Dans un premier temps, on prédéfinit au moins deux masses de gaz enfermées m₁, ..., m_n. On en déduit au moins deux rapports thermiques γ₁, ..., γ_n qui correspondent aux masses de gaz enfermées prédéfinies m₁, ..., m_n. Ensuite, on calcule au moins deux fréquences d'oscillations de pression f₁, ..., f_n qui correspondent aux masses de gaz enfermées prédéfinies m₁, ..., m_n., au moyen de l'équation $f=\sqrt{\frac{\mathit{\gamma PV}}{{\mathit{m\lambda }}^2}}$ qui découle de la formule $m=\frac{\mathit{\gamma PV}}{{\lambda }^2{f}^2}$. Ces valeurs sont comparées (COMP) avec les oscillations de pression mesurée OSC, pour choisir la fréquence d'oscillations f parmi les fréquences d'oscillations de pression f₁, ..., f_n qui se rapproche le plus de la fréquence des oscillations de pression mesurée OSC. La masse de gaz enfermée déterminée par le procédé selon l'invention est alors la masse de gaz enfermée parmi les masses de gaz enfermées m₁, ..., m_n qui correspond à la fréquence d'oscillations de pression f choisie.

En variante, la fréquence d'oscillation f peut être obtenue par toute autre méthode de calibration, par exemple par une méthode de minimisation d'une fonction objectif.

Alternativement, la fréquence d'oscillation peut être déterminée au moyen d'une moyenne sur une fenêtre temporelle desdites oscillations déterminées lors de l'étape 2).

Comme il va de soi, l'invention ne se limite pas aux seules formes de réalisation des étapes du procédé décrites ci-dessus à titre d'exemple, elle embrasse au contraire toutes les variantes de réalisation.

Les caractéristiques et avantages du procédé selon l'invention apparaîtront plus clairement à la lecture de l'exemple d'application ci-après.

L'exemple concerne la comparaison de la masse de gaz enfermée déterminée selon le procédé selon l'invention, avec la masse de gaz enfermée déterminée par simulation LES. Pour cet exemple, on considère un moteur à combustion interne suralimenté ayant une vitesse de rotation à 5500 tr/min, et on considère 23 cycles critiques qui ont été identifiés comme ayant des niveaux d'oscillations de pressions variables. Pour cet exemple, la chambre de combustion du moteur à combustion interne est équipée de deux capteurs de pression à des positions différentes au sein de la chambre de combustion.

La figure 4 est un graphique illustrant la masse m de gaz enfermée dans la chambre de combustion en fonction du numéro N du cycle considéré. La masse m utilisée dans ce graphique est une valeur relative qui correspond au ratio de la masse de gaz enfermée dans la chambre de combustion obtenue par le procédé selon l'invention par rapport à la masse de gaz enfermée dans la chambre de combustion obtenue par simulation. Ainsi, une valeur proche de 1 indique une bonne corrélation entre les valeurs simulées et les valeurs obtenues par le procédé selon l'invention. A chaque N correspond deux masses de gaz enfermées, chacune étant obtenue par un des deux capteurs de pression utilisé. On remarque que les valeurs (représentés par les points) et leur moyenne représentée par le segment pointillé sont très proches de la valeur un. Ainsi, le procédé selon l'invention permet de déterminer de manière précise la masse de gaz enfermée dans la chambre de combustion. De plus, les points et les moyennes des deux capteurs de pression sont proches. Par conséquent, le procédé selon l'invention est bien indépendant de la position du capteur dans la chambre de combustion.