METHOD OF CONTROLLING THE COMBUSTION OF A SPARK-IGNITION ENGINE USING COMBUSTION TIMING CONTROL

Abstract

A method of controlling the combustion of a spark-ignition engine having application to gasoline engines is disclosed. An engine control system controls actuators so that the values of physical parameters linked with the combustion of a mixture of gas and fuel in a combustion chamber are equal to their setpoint values, to optimize the combustion. A setpoint value is determined for an ignition crank angle of the fuel mixture which is then corrected before the physical parameters reach their setpoint values. A correction to be applied to this ignition angle setpoint value is calculated so that the crank angle CAy is equal to its setpoint value. Finally, the engine control system controls the ignition of the mixture in the combustion chamber when the crank angle is equal to the corrected setpoint value to optimize combustion.

Description

BACKGROUND OF THE INVENTION

1. Field of the Invention

The present invention relates to the field of engine control and more particularly to the combustion control of spark-ignition engines.

2. Description of the Prior Art

Operation of a (gasoline type) spark-ignition engine is based on the combustion of a mixture of air, burnt gas and fuel. The engine cycle can be broken down into four phases (FIG. 1):

The intake phase (ADM): the intake valve allows the mixture of air and of burnt gas into chamber CHB. The air is taken from the outside environment of the engine. The burnt gas is taken from exhaust manifold ECH and sent back to the intake manifold (exhaust gas recirculation EGR) and/or sucked back by the exhaust valve (internal exhaust gas recirculation iEGR). The fuel is injected during the intake phase. The Variable Valve Timing (VVT) device allows a time lag to be applied to the intake (VVTa) and exhaust (VVTe) valve lift profiles. This has a direct impact on the gas composition and on the turbulence in the combustion chamber;

The compression phase: after the intake valve closes (IVC: Intake Valve Closing), piston PIS compresses the gas;

The combustion phase: spark plug BOUG produces a spark that initiates the combustion of the mixture of air, burnt gas and fuel which ignites while releasing the chemical energy available in the fuel, thus creating an overpressure that pushes the piston backwards;

The expansion phase: once the piston has gone back down again, the exhaust valve opens and the gas mixture is then discharged through the exhaust manifold.

The goal of engine control is to supply the driver with the torque required while minimizing the noise and pollutant emissions. Control of the amounts of the different gases and of the fuel therefore has to be adjusted as finely as possible.

To carry out combustion control of a spark-ignition engine, there are known methods allowing determination of the combustion medium by means of detectors mounted in the engine. There are methods based on direct measurements in the combustion chamber, such as those provided by:

Cylinder Pressure Detectors: Paljoo Yvon et al., “Closed-loop Control of Spark Advance and Air-fuel Ratio in SI Engines Using Cylinder Pressure”, Society of Automotive Engineering World Congress, 2000-01-0933,

Ionization Current Detectors in the Combustion Chamber: Lars Eriksen et al., “Closed loop Ignition Control by Ionization Current Interpretation”, SAE 1997 Transactions, Journal of Engines, Vol. 106, Section 3, pp. 1216-1223, 1997.

However, using such detectors in standard vehicles is difficult due to their considerable cost. Furthermore, these detectors are generally subject to relatively fast drifts.

There are also methods wherein the amounts and timing are optimized on each static working point (speed and torque) so as to bring out an ideal strategy at each point. An engine test bench calibration is therefore performed in order to obtain the optimum values for the main three data sets:

Air Loop:

the mass of air M_air and of burnt gas M_gb required in the combustion chamber;

the pressure and the temperature of these gases in the combustion chamber; and

the position of the variable valve lift devices (VVT) and notably the intake valve closure angle denoted by θ_ivc.

These thermodynamic and physical variables X_air=(M_air, M_bg, P, T, θ_IVC) are represented by X_air.

Fuel Loop:

The mass of fuel M_f injected into the combustion chamber (injection directly into the chamber or indirectly into the intake pipe), X_fuel=(M_f).

Ignition Loop:

The crank angle θ_all at which the spark appears (via the plug), denoted by X_all=(θ_all).

However, these strategies are insufficient in transient phases. In fact, during transition phases from one working point to another (change in the vehicle speed or in the road profile), the engine control supervises the various actuators present in the engine to guarantee the desired torque while minimizing the noise, the pollutant emissions and the consumption. This is thus translated into the change from the values of the parameters of the initial point to the values of the parameters of the final point:

$\begin{array}{cc}\{\begin{array}{c}{X}_{\mathrm{air}}^{\mathrm{initial}}\to X_{\mathrm{air}}^{\mathrm{final}}\\ X_{\mathrm{fuel}}^{\mathrm{initial}}\to X_{\mathrm{fuel}}^{\mathrm{final}}\\ X_{\mathrm{all}}^{\mathrm{initial}}\to X_{\mathrm{all}}^{\mathrm{final}}\end{array}& \begin{array}{c}\begin{array}{c}\left(a\right)\\ \left(b\right)\end{array}\\ \left(c\right)\end{array}\end{array}$

Now, there are two time scales in the engine. The faster one (50 Hz) corresponds to the entire combustion phenomenon (1 engine cycle). On this scale, the injection (X_fuel) and the ignition (X_all) strategy can be changed to control the combustion. The slower one (1 Hz) corresponds to the gas dynamics in the engine manifolds (intake, exhaust, burnt gas recirculation) and the inertia of the actuators (turbocompressor TC). The strategy of this air loop (X_air) cannot be changed faster.

With current methods, the controlled variables (X_air, X_fuel, X_all) do therefore not reach at the same time their setpoint values because of the difference in dynamics. The objectives regarding torque production, namely, consumption, pollutants, and noise are thus met in the static phases (the dynamic loops are stabilized at their reference values). On the other hand, if precautions are not taken in the transient phases, part of the parameters reach nearly instantaneously their final setpoint value whereas the other part is still at the initial setpoint value. This results in the engine then producing more pollutant emissions or noise and can even cause stopping in some cases.

Furthermore, without cylinder pressure detectors, the known methods do not allow controlling combustion timing during the transient phases. Now, as illustrated by FIGS. 2 and 3, this is not sufficient to ensure proper operation of the engine under transient conditions.

SUMMARY OF THE INVENTION

The invention relates to a method providing control of the combustion of a spark-ignition engine, notably under transient conditions, while overcoming prior art problems. The method achieves this, on the one hand, by controlling three dynamic loops separately and, on the other hand, by correcting the reference value of the ignition angle via control of angle CA₅₀.

The invention thus relates to a method of controlling the combustion of a spark-ignition engine, comprising: determining setpoint values for physical parameters linked with the combustion of a mixture of gas and of fuel in a combustion chamber and a setpoint value (θ_all)_ref for an ignition crank angle for the mixture, the setpoint values being determined to optimize combustion, and an engine control system that controls actuators so that values of the physical parameters are equal to the setpoint values.

The method comprises the following stages:

correcting setpoint value (θ_all)_ref before the physical parameters reach their setpoint values, by calculating a correction dθ_all to be applied to the setpoint value (θ_all)_ref so that a crank angle CA_y at which y percent of the fuel is consumed during combustion is equal to a setpoint value of the angle for an optimized combustion;

the engine control system controls the ignition of the mixture in the combustion chamber when the crank angle is equal to the corrected setpoint value (θ_all)_ref in order to keep combustion optimal.

According to the invention, correction dθ_all can be determined by accounting for differences dp between real values p of the physical parameters and the setpoint values p^ref of the physical parameters. It is therefore possible to use a combustion model defined by a differential equation allowing modelling an evolution over time of a consumed fuel mass, and by linearizing the combustion model to p around setpoint values p^ref, then by calculating a first-order solution for the correction to be made so that correction dθ_all is proportional to differences dp.

According to an embodiment, correction dθ_all can thus be determined by applying the following stages:

determining the real values of the physical parameters;

calculating the differences dp between the real values and the setpoint values;

determining the setpoint value (CA_y)_ref of crank angle CAy by of a numerical integration of the combustion model by assigning to each parameter of the model the setpoint value;

calculating linearization matrix Λ of the combustion model by linearizing the combustion model to p around setpoint values p^ref; and

calculating correction dθ_all by the formula:

dθ _all=((CA_y)_ref−(θ_all)_ref)·Λ·dp

According to a preferred embodiment, crank angle CAy is the crank angle at which fifty percent of the fuel is consumed during combustion.

According to the invention, the physical parameters can be selected from among at least the following parameters upon valve closing: pressure in the combustion chamber (P_IVC), temperature in the combustion chamber (T_IVC), ratio (X_IVC) between a burnt gas mass and a total gas mass in the combustion chamber, air mass (M_IVC) in the cylinder and closure angle (θ_ivc) of an intake valve.

Finally, it is also possible to adapt a mass of fuel injected into the combustion chamber before the physical parameters reach their setpoint values by controlling the richness of the combustion.

BRIEF DESCRIPTION OF THE DRAWINGS

Other features and advantages of the method according to the invention will be clear from reading the description hereafter of embodiments given by way of non limitative example, with reference to the accompanying figures wherein:

FIG. 1 shows the various phases of a combustion cycle of a spark-ignition engine;

FIG. 2 illustrates a combustion chronology as a function of the crank angle according to three combustion control situations: optimum control (performed in stabilized phase), current control in transient phase without CA₅₀ control and desired control in transient phase with CA₅₀ control;

FIG. 3 illustrates the three energy release curves Q as a function of crank angle θ for the three situations described in FIG. 2; and

FIG. 4 illustrates a calculation scheme for correction dθ_all of the ignition angle.

DETAILED DESCRIPTION OF THE INVENTION

The method according to the invention allows controlling the combustion progress of a spark-ignition engine, in a static phase as well as in a transient phase. It comprises separate and independent control of the air loop (slow loop) and of the fuel and ignition loops (fast loops), through adaptation of the fast loop dynamics to be coherent with the air loop. The method thus allows adaptation of X_fuel and X_all to keep the characteristics of the combustion required (through the driver's torque request). The impact on emissions and noise is thus limited while ensuring the required torque to the driver.

According to this method, control of the combustion of a spark-ignition engine is carried out in five stages:

1—Determining Setpoint Values for Various Physical Parameters

During transition phases from one working point to another (change in the vehicle speed or in the road profile), the engine control supervises the various actuators present in the engine to guarantee the desired torque while minimizing the noise, the pollutant emissions and the consumption. This is thus translated into the change from the values of parameters X_air, X_fuel and X_all of an initial point to the values of the parameters of a final point:

$\begin{array}{cc}\{\begin{array}{c}{X}_{\mathrm{air}}^{\mathrm{initial}}\to X_{\mathrm{air}}^{\mathrm{final}}\\ X_{\mathrm{fuel}}^{\mathrm{initial}}\to X_{\mathrm{fuel}}^{\mathrm{final}}\\ X_{\mathrm{all}}^{\mathrm{initial}}\to X_{\mathrm{all}}^{\mathrm{final}}\end{array}& \begin{array}{c}\begin{array}{c}\left(a\right)\\ \left(b\right)\end{array}\\ \left(c\right)\end{array}\end{array}$

The final values are defined to optimize combustion, that is, to burn a maximum amount of fuel in order to minimize emissions and consumption while minimizing the noise. These final values optimizing the combustion are referred to as setpoint values. The engine control enforces these setpoint values.

The important physical parameters regulated by the air loop are the pressure, the temperature, the chemical composition of the gases in the chamber and the intake valve closing angle. Ideally, these parameters reach their setpoint value instantaneously. In reality, the slowness of the air loop results in an error on these parameters X_air between their setpoint value and their real value, throughout the transition phase. Consequently, the thermodynamic parameters (mass, pressure, temperature and burnt gas rate) of the gas feed sucked in the cylinder are different from their setpoint value. The fuel and ignition loop control is adapted to the errors on the following parameters:

P: The pressure in the combustion chamber. It depends on crank angle θ,

T: The temperature in the combustion chamber. It depends on crank angle θ,

X: The ratio between the burnt gas mass and the total gas mass in the combustion chamber (parameter between 0 and 1). It depends on crank angle θ,

M_air: The mass of air trapped in the cylinder.

The value of these parameters is determined upon valve closing (IVC):

P_IVC: The pressure in the combustion chamber upon valve closing;

T_IVC: The temperature in the combustion chamber upon valve closing;

X_IVC: The ratio between the burnt gas mass and the total gas mass in the combustion chamber upon valve closing;

M_IVC: The mass of air in the cylinder upon valve closing;

θ_IVC: The closing angle of the intake valve; it directly influences the turbulence in the combustion chamber.

The values of these five parameters are continuously determined. Therefore, it is assumed that composition (X_IVC) and pressure (P_IVC) in the cylinder upon valve closing are the same as those in the intake manifold where measurements are available (through detectors or estimators). T_IVC is estimated by means of the ideal gas law

$T_{\mathrm{IVC}}=\frac{P_{\mathrm{IVC}}V_{\mathrm{IVC}}}{{\mathrm{RM}}_{\mathrm{IVC}}}$

where R is the ideal gas constant (R=287) and M_IVC is the mass sucked by the cylinder that is measured by a flowmeter.

For these five physical parameters linked with the intake of gaseous oxidizer in the combustion chamber of the engine, the setpoint values are respectively denoted by: P_ref, T_ref, X_ref, M_ref and (θ_ivc)_ref.

These setpoint values are obtained from a setpoint map established on an engine test bench. The setpoint values of these parameters are given by the optimum point mapped on the test bench (values that these parameters must reach). These setpoint values are determined to optimize the combustion. These parameters are related by the ideal gas relation (PV=MRT) but, for simplicity reasons, this relation is not directly made explicit. This does not affect the method provided in any way.

According to the invention, the parameter which is controlled is the mixture ignition angle: θ_all. Its reference value (given by the optimum point mapped on the test bench) is denoted by (θ_all)_ref. The parameter to be kept constant is crank angle CAy, which is the angle at which y percent of the fuel is consumed during combustion. It is attempted to maintain this angle at a setpoint value (CA_y)_ref of this angle for an optimum combustion. According to a preferred embodiment, the half combustion angle CA₅₀ is used. It is the crank angle at which 50% of the fuel has been consumed during the optimized combustion (combustion obtained with the setpoint values).

2—Air Loop Control (Slow Loop)

Once setpoint values P_ref, T_ref, X_ref, M_ref and (θ_ivc)_ref are determined, an engine control system controls actuators so that the values of the physical parameters P_IVC, T_IVC, X_IVC, M_IVC and θ_ivc equal to their setpoint values P_ref, T_ref, X_ref, M_ref and (θ_ivc)_ref.

3—Fuel Loop Adaptation (Fast Loop)

Adapting the control of the fuel mass injected into the air loop dynamics is conventionally achieved by controlling the combustion richness: in fact, removing from the exhaust gas from gasoline engines can be accomplished by a three-way catalyst. It allows efficient treatment of the CO, HC and NOx produced by the combustion, provided that the exhaust gas is globally neither oxidizing nor reducing. The combustion richness (λ) is defined as the excess air mass M_air in relation to the fuel mass M_f brought to the same ratio in the case of stoichiometric combustion (this stoichiometric ratio is denoted by PCO). Thus,

$\lambda =\frac{M_f}{M_{\mathrm{air}}}\frac{1}{\mathrm{PCO}}.$

The removing of pollution is thus efficient when the richness is close to 1. The control strategy for the fuel mass that is injected is thus reduced to the estimation of the air mass sucked into the cylinder from the air loop parameters. Estimation of the air mass M_air then allows applying the command

$M_f=\frac{1}{\mathrm{PCO}}M_{\mathrm{air}}.$

4—Ignition Loop Adaptation

a. Calculating the Correction of the Ignition Angle Setpoint Value (θ_all)_ref

The conventional control strategy for ignition angle θ_all is a prepositioning depending on the engine speed and on the estimation of the air mass sucked in the cylinder (via mapping). Unlike the fuel mass control, this strategy is not optimal. In fact, if the fuel mass injected provides a torque potential, it is the ignition timing that guarantees good exploitation of this potential.

If ignition occurs too early, combustion also takes place early, causing a pressure increase during the compression phase, which goes against torque production;

if ignition occurs too late, combustion also takes place late and it is the entire expansion phase that will occur at a lower pressure. The torque produced is then lower.

The CA₅₀ (crank angle at which 50% of the fuel is burned) is the crank angle that allows accounting for this combustion timing. It is conventionally admitted that each engine has a fixed reference crank angle (CA₅₀)_ref depending on the engine's technical data. The ignition strategy is then optimal if the CA₅₀ is regulated to its reference value (CA₅₀)_ref.

It is therefore necessary to take into account all the thermodynamic and physical parameters that influence the combustion in order to best control the ignition angle. Conventional ignition angle prepositioning, depending on the engine speed and the air mass sucked, is thus the start of an adaptation of the ignition loop to the slow parameters of the air loop, but it is not complete. In fact, other parameters of the air loop influence the combustion:

the pressure in the combustion chamber;

the temperature in the combustion chamber;

the burnt gas rate in the combustion chamber; and the turbulence in the combustion chamber (via θ_ivc).

If the air loop control was perfect, the five parameters P_IVC, T_IVC, X_IVC, M_IVC and θ_ivc would reach their reference values P_ref, T_ref, X_ref, M_ref and (θ_ivc)_ref instantaneously. In reality, in the transient phase, parameters P_IVC, T_IVC, X_IVC, M_IVC and θ_ivc are different from their reference value. The content of the cylinder upon valve closing is therefore different from the reference content for which the ignition strategy was mapped.

The errors on these parameters upon valve closing therefore have to be taken into account to modify the ignition angle so as to keep a combustion that is as close as possible to the reference combustion (optimized combustion). Therefore, the following exists:

p the vector (P_IVC; T_IVC; M_IVC; X_IVC; θ_ivc);

p^ref the vector (P_ref; T_ref; M_ref; X_ref; (θ_ivc)_ref); and

dp the vector p−p^ref=(dP; dT; dM; dX; dθ_ivc), with:

• • dP=P_IVC−P_ref, dT=T_IVC−T_ref, dM=M_IVC−M_ref, dX=X_IVC−X_ref and dθ=θ_ivc−(θ_ivc)_ref.

A new corrected ignition angle (θ_all)_ref+dθ_all is therefore sought so that angle CA₅₀ is at its reference value (dCA₅₀=CA₅₀−(CA₅₀)_ref=0). Therefore, dθ_all is sought such that (see FIGS. 2 and 3 for the three situations):

if there is no error, that is, if all the parameters have reached their reference value (dp=0), the situation of the reference working pointportion is that dθ_all=0 (situation {circle around (1)}),

if the parameters have not reached their reference value (dp≠0), the combustion velocity is not identical to that of the reference combustion. Thus, a phase shift of the combustion and angle CA₅₀ has not reached its reference value (situation {circle around (2)});

in order to counterbalance the errors dp≠0, an angular correction dθ_all≠0 on the ignition angle is introduced to have the same phasing CA₅₀ (situation {circle around (3)}).

FIG. 2 illustrates a combustion chronology according to three situations. For each situation, the horizontal axis represents crank angle θ. These axes comprise: the setpoint value (θ_all)_ref of the ignition angle, the ignition angle θ_all and the corrective term dθ_all. FIG. 3 illustrates the three energy release curves Q as a function of crank angle θ for the three situations described above (FIG. 2).

In order to determine correction dθ_all, modelling of the combustion system is performed. According to a particular embodiment, it is possible to use a combustion model defined by a differential equation allowing modelling of the evolution over time of the mass of fuel consumed by the combustion. Such a combustion model can thus be written in the compact form as follows (appendix 1 illustrates an example of such a combustion model):

$\begin{array}{cc}\{\begin{array}{c}\frac{x}{\theta }=f\left(x,y,p,\theta \right)\\ \frac{y}{\theta }=g\left(x,y,p,\theta \right)\\ x\left({\theta }_{\mathrm{all}}\right)=0\\ y\left({\theta }_{\mathrm{all}}\right)=h\left(p\right)\end{array}\theta \in \left[{\theta }_{\mathrm{all}},{\mathrm{CA}}_{50}\right]& \left(1\right)\end{array}$

with:

• • x: mass fraction of burnt fuel (x has dimension 1×1)
  • y: other variables whose dynamics are necessary for the combustion model (pressure, temperature . . . ). y has dimension 1×n with nε*.
  • p: parameters of the air loop to be compensated during transient phases
    • p has dimension 1×n with nε*
  • θ_all: mixture ignition angle
  • f, g and h are entirely known functions (see Appendix 1 for example).

It can be noted that the ignition angle control method according to the invention is applicable to any combustion model in differential equation form.

Estimating the Correction of the Ignition Angle Setpoint Value (θ_all)_ref

Correction calculation is carried out by linearizing the combustion model to p around reference values p^ref by introducing differences dp.

Considering the complex form of the combustion model, it is difficult to find an analytical expression for correction dθ_all. Consequently, a first-order solution for the correction to be made is sought. The correction thus is proportional to the air loop error dp.

The correction is obtained as follows:

A sensitivity analysis of differential equation (1) with respect to the two variables p and θ_all is carried out. The first-order variation of the states of the differential equation (dx, dy) consecutive to any variations of p and θ_all around their reference values is obtained,

The fact that, in the presence of error dp and of correction dθ_all, the CA₅₀ should not be affected (dCA₅₀=0) is translated into a condition on the sensitivity of state x: dx(θ_all)=0. Thus, an equation is obtained in which dp and dθ_all appear. Through inversion of this equation, an expression for dθ_all as a function of dp and of all the variables present in the model is obtained.

The correction obtained is written as follows:

$d{\theta }_{\mathrm{all}}=\left({\left({\mathrm{CA}}_{50}\right)}_{\mathrm{ref}}-{\left({\theta }_{\mathrm{all}}\right)}_{\mathrm{ref}}\right)·\Lambda ·\mathrm{dp}$ $\mathrm{with}\text{:}$ $\Lambda =\frac{{\int }_{{\left({\theta }_{\mathrm{all}}\right)}_{\mathrm{ref}}}^{{\left({\mathrm{CA}}_{50}\right)}_{\mathrm{ref}}}\Gamma \left({\theta }^T\right)\left(\begin{array}{c}{B}^{11}\left(\theta \right)\\ B^{21}\left(\theta \right)\end{array}\right)\theta +{\Gamma \left({\left({\theta }_{\mathrm{all}}\right)}_{\mathrm{ref}}\right)}^T\left(\begin{array}{c}0\\ \frac{g}{p}\left(p^{\mathrm{ref}}\right)\end{array}\right)}{{\int }_{{\left({\theta }_{\mathrm{all}}\right)}_{\mathrm{ref}}}^{{\left({\mathrm{CA}}_{50}\right)}_{\mathrm{ref}}}{\Gamma \left(\theta \right)}^T\left(\begin{array}{c}{B}^{12}\left(\theta \right)\\ B^{22}\left(\theta \right)\end{array}\right)\theta }$ $B^{11}\left(\theta \right)=\frac{\partial f}{\partial p}\left(x^{\mathrm{ref}},y^{\mathrm{ref}},p^{\mathrm{ref}},\theta \right)$ $B^{12}\left(\theta \right)=f\left(x^{\mathrm{ref}},y^{\mathrm{ref}},p^{\mathrm{ref}},\theta \right)+\left(\theta -{\left({\mathrm{CA}}_{50}\right)}_{\mathrm{ref}}\right)\frac{\partial f}{\partial p}\left(x^{\mathrm{ref}},y^{\mathrm{ref}},p^{\mathrm{ref}},\theta \right)$ $B^{21}\left(\theta \right)=\frac{\partial g}{\partial p}\left(x^{\mathrm{ref}},y^{\mathrm{ref}},p^{\mathrm{ref}},\theta \right)$ $B^{22}\left(\theta \right)=g\left(x^{\mathrm{ref}},y^{\mathrm{ref}},p^{\mathrm{ref}},\theta \right)+\left(\theta -{\left({\mathrm{CA}}_{50}\right)}_{\mathrm{ref}}\right)\frac{\partial g}{\partial \theta }\left(x^{\mathrm{ref}},y^{\mathrm{ref}},p^{\mathrm{ref}},\theta \right)$ $\Gamma \left(\theta \right)=\exp \left({\int }_{\theta }^{{\left({\mathrm{CA}}_{50}\right)}_{\mathrm{ref}}}\left(\begin{array}{cc}\frac{\partial f}{\partial x}\left(x^{\mathrm{ref}},y^{\mathrm{ref}},p^{\mathrm{ref}},\theta \right)& \frac{\partial g}{\partial x}\left(x^{\mathrm{ref}},y^{\mathrm{ref}},p^{\mathrm{ref}},\theta \right)\\ \frac{\partial f}{\partial y}\left(x^{\mathrm{ref}},y^{\mathrm{ref}},p^{\mathrm{ref}},\theta \right)& \frac{\partial g}{\partial y}\left(x^{\mathrm{ref}},y^{\mathrm{ref}},p^{\mathrm{ref}},\theta \right)\end{array}\right)\theta \right)\left(\begin{array}{c}1\\ 0\end{array}\right)$

b. Ignition Loop Adaptation (Fast Loop)

The engine control system activates the fuel ignition system in the combustion chamber when the crank angle is equal to the corrected setpoint value (θ_all)_ref+dθ_all in order to keep combustion optimal.

One interest of the method is to directly relate the air loop errors to the correction to be applied to the ignition command via matrix Λ. The latter is entirely calculable: it only depends on the combustion model, on reference values P_ref, T_ref, X_ref and M_ref and on a certain number of known constants.

By applying the previous correction to the ignition angle, it can be ensured (first order) that angle CA₅₀ is at its reference value. Little by little, the air loop leads errors dP, dT, dM and dX move towards zero, and the correction therefore disappears in the stabilized static phases. The control strategy is diagrammatically shown in FIG. 4. This figure illustrates a calculation scheme for correction dθ_all of the ignition angle. After estimating or measuring (EST−ACT), the real values of parameters p, determining (DET−CONS) setpoint values p^ref of these parameters, and (θ_all)_ref, (CAL−MAT) is calculated which is the linearization matrix of the combustion model. Then the following coefficient is calculated: ((CA₅₀)_ref−(θ_all)_ref)·Λ.

The value of (CA₅₀)_ref (value of the half combustion angle of the reference combustion) is required. Thus, the differential system (1) (combustion model) with reference values p^ref and (θ_all)_ref is used as the initial conditions. The following system if obtained:

$\{\begin{array}{c}\frac{x^{\mathrm{ref}}}{\theta }=f\left(x^{\mathrm{ref}},y^{\mathrm{ref}},p^{\mathrm{ref}},\theta \right)\\ \frac{y^{\mathrm{ref}}}{\theta }=g\left(x^{\mathrm{ref}},y^{\mathrm{ref}},p^{\mathrm{ref}},\theta \right)\\ x\left({\left({\theta }_{\mathrm{all}}\right)}^{\mathrm{ref}}\right)=0\\ y\left({\left({\theta }_{\mathrm{all}}\right)}^{\mathrm{ref}}\right)=h\left(p^{\mathrm{ref}}\right)\end{array}\theta \in \left[{\left({\theta }_{\mathrm{all}}\right)}^{\mathrm{ref}},{\left({\mathrm{CA}}_{50}\right)}^{\mathrm{ref}}\right]$

By numerical integration of this system, the value of angle θ is determined when x_ref=0.5. This angle corresponds to the value of (CA₅₀)_ref.

Finally, (CAL−COR) is calculated which is the correction dθ_all:

dθ _all=((CA₅₀)_ref−(θ_all)_ref)·Λ·dp

The method according to the invention thus allows ensured combustion control of a spark-ignition engine by controlling the three dynamic loops separately and by correcting the reference value of the ignition angle. This correction is determined in such a way that angle CA₅₀ is at its reference value, and the same energy release as for the reference combustion (optimized) is consequently obtained.

APPENDIX 1

Combustion Model

An example of a combustion model defined by a differential equation is presented in the following document:

F.-A. Lafossas, et al., “Application of a New 1d Combustion Model to Gasoline Transient Engine Operation,” in Proc. SAE World Congress, no. 2005-01-2107, 2005.

This model represents the volume of the cylinder in two zones (the burnt zone and the unburnt zone) separated by the flame front (modelled as an infinitely thin layer). Throughout combustion, the flame propagates from the burnt zone to the unburnt zone. The equations of the model are as follows:

$\{\begin{array}{c}\frac{\mathrm{dx}}{d\theta }=\frac{1}{N_e}\left(\frac{C_1·p_2}{{V\left(t\right)}^{\mathrm{\kappa \gamma }+1}}·y^{\kappa +1/\gamma }+\frac{C_2p_1^{-1/\gamma }\sqrt{k\left(\theta ,p_3\right)}}{V\left(t\right)}y^{1/\gamma }\tanh \left(\frac{1}{r_0}{\left(\frac{3V_{\mathrm{fl}}}{4\pi }\right)}^{1/3}-1\right)\right)S_{\mathrm{geo}}\left(V_{\mathrm{fl}},t\right)\\ V_{\mathrm{fl}}=\max \left(f_{\mathrm{vol}},V\left(\theta \right)\left(1-\left(1-x\right)p_1^{1/\gamma }y^{-1/\gamma }\right)\right)\\ S_{\mathrm{geo}}\left(V_{\mathrm{fl}},\theta \right)=\{\begin{array}{cc}\sqrt[3]{36\pi }·V_{\mathrm{fl}}^{2/3}& \mathrm{if}V_{\mathrm{fl}}<\frac{\pi }{6}{\left(\frac{V\left(\theta \right)}{A}\right)}^3\\ 2\sqrt{2\pi V_{\mathrm{fl}}\frac{V\left(\theta \right)}{A}+\frac{4}{3}{{\pi }^2\left(\frac{V\left(\theta \right)}{A}\right)}^4}& \mathrm{if}V_{\mathrm{fl}}>\frac{\pi }{6}{\left(\frac{V\left(\theta \right)}{A}\right)}^3\end{array}\\ k\left(\theta ,p_3\right)=C_3p_4\left(-{\int }_{p_3}^{\theta }{\zeta }^1\left(z,p_3\right)e^{\frac{C_4}{N_e}\left(z-p_3\right)}\mathrm{dz}+C_5\right)e^{-\frac{C_4}{N_e}\left(\theta -p_3\right)}\\ \zeta \left(\theta ,p_3\right)={V\left(\theta \right)}^2{\left(\frac{\theta -{\theta }_{\mathrm{TDC}}}{p_3-{\theta }_{\mathrm{TDC}}}\right)}^2{\delta }_{\theta <{\theta }_{\mathrm{TDC}}}\\ \frac{\mathrm{dy}}{d\theta }=Q_{\mathrm{LHV}}M_f\left(\gamma -1\right){V\left(\theta \right)}^{\gamma -1}\frac{\mathrm{dx}}{d\theta }\\ x\left({\theta }_{\mathrm{all}}\right)=0\\ y\left({\theta }_{\mathrm{all}}\right)=p_1\end{array}$
with:
Symb.	Quantity	Unit
Ne	Engine speed	rpm
θ	Crank angle	[deg]
θivc	Crank angle at ivc	[deg]
θTDC	Crank angle of piston upper position	[deg]
θall	Crank angle of mixture ignition	[deg]
V(θ)	Cylinder volume	m3
Vivc	Cylinder volume at ivc	m3
P(θ)	Pressure in cylinder	Pa
Pivc	Pressure in cylinder at ivc	Pa
T(θ)	Temperature in mixture	K
Tivc	Temperature in mixture at ivc	K
Tu(θ)	Temperature of unburnt zone	K
AFRS	air/fuel ratio at stoichiometry	—
minj	mass of fuel injected	kg
mair	mass of air admitted	kg
mbg	burnt gas rate in combustion chamber	kg
mf	mass of fuel burnt during combustion (0 to Mf)	kg
ρu	density in the unburnt zone	kg/m3
(ρu)ivc	density in the unburnt zone at ivc	kg/m3
Yu	mass fraction of fuel in the unburnt zone	—
U	laminar flame velocity	m/s
$\Xi$	wrinkling due to turbulence	—
γ	adiabatic index	—
PMI	mean indicated pressure	bar
QLHV	mass energy available in the fuel	J/kg
A	piston surface area	m2
Sfl	flame surface area	m2
Sgeo	geometric flame surface area (without wrinkling)	m2
fvol	minimum flame volume (flame volume initiation)	m3
x	burnt fuel mass fraction	—
y	P.Vγ (variable representing the pressure in the chamber) Pa m3γ

The air loop parameters to be compensated are all grouped together in the three parameters as follows:

$p_1=P_{\mathrm{ivc}}V_{\mathrm{ivc}}^{\gamma }$ $p_2={\left(\frac{\mathrm{Tivc}}{P_{\mathrm{ivc}}^{\frac{\gamma }{\gamma -1}}}\right)}^{\alpha }\frac{1-2.1\frac{m_{\mathrm{bg}}}{m_{\mathrm{bg}}+m_{\mathrm{air}}}}{P_{\mathrm{ivc}}^{1/\gamma }V_{\mathrm{ivc}}}$ $p_3={\theta }_{\mathrm{ivc}}$

Finally, the parameters of the model are: C₁, C₂, C₃, C₄, α. A numerical value example (international system of units SI) for these parameters is given in the table below:

	Parameter
	C1	C2	C3	C4	α
	Value	2.92 e−5	2.11	5.34 e7	1.67 e−2	2.12

The model can thus eventually be written in the condensed form as follows:

$\{\begin{array}{c}\frac{x}{\theta }=f\left(x,y,p,\theta \right)\\ \frac{y}{\theta }=g\left(x,y,p,\theta \right)\\ x\left({\theta }_{\mathrm{all}}\right)=0\\ y\left({\theta }_{\mathrm{all}}\right)=C^Tp\end{array}\theta \in \left[{\theta }_{\mathrm{all}},{\mathrm{CA}}_{50}\right]$

with p=(p₁, p₂, p₃)^T, C=(1, 0, 0)^T and f and g defined by the equation of the combustion model at the beginning of the appendix.

In these equations, term CA₅₀ can be readily substituted for any angle CA_y.

Claims

1-7. (canceled)

8) A method of controlling combustion of a spark-ignition engine, comprising: determining setpoint values for physical parameters linked with the combustion of a mixture of gas and of fuel in a combustion chamber and a setpoint value for an ignition crank angle for the mixture, the setpoint values being determined to optimize combustion;using an engine control system to control actuators so that the values of the physical parameters are equal to the setpoint values by correcting the setpoint value of the ignition crank angle before the physical parameters reach the setpoint values, by calculating a correction dθall of the setpoint value of the ignition crank angle, so that a crank angle CAy at which y percent of the fuel is consumed during combustion is equal to a setpoint value of the angle for an optimized combustion by using a combustion model including a differential equation representing an evolution over time of a consumed fuel mass; andcontrolling with the engine control system the ignition of the mixture in the combustion chamber when the crank angle is equal to the corrected ignition crank angle setpoint value to optimize combustion.

9) A method as claimed in claim 8, wherein the correction dθall is determined by accounting for differences dp between real values p of the physical parameters and the setpoint values pref of the physical parameters.

10) A method as claimed in claim 9, wherein the correction dθall is determined by linearizing the combustion model to p around setpoint values pref and then calculating a first-order solution for the correction dθall to be made so that the correction dθall is proportional to the differences dp.

11) A method as claimed in claim 10, wherein the correction dθall is determined by a method comprising: determining the real values of the physical parameters;calculating the differences dp between the real values and the setpoint values;determining the setpoint value of crank angle CAy by a numerical integration of the combustion model by assigning to each parameter of the model a setpoint value thereof;calculating a linearization matrix Λ of the combustion model by linearizing the combustion model to p around setpoint values pref;calculating the correction dθall by the formula: dθall=((CAy)ref−(θall)ref)·Λ·dp where (θall)ref is the setpoint value of the ignition crank angle of the mixture and (CAy)ref is the setpoint value of crank angle CAy.

12) A method as claimed in claim 8, wherein the crank angle CAy is the crank angle at which fifty percent of the fuel is consumed during combustion.

13) A method as claimed in claim 9, wherein the crank angle CAy is the crank angle at which fifty percent of the fuel is consumed during combustion.

14) A method as claimed in claim 10, wherein the crank angle CAy is the crank angle at which fifty percent of the fuel is consumed during combustion.

15) A method as claimed in claim 11, wherein the crank angle CAy is the crank angle at which fifty percent of the fuel is consumed during combustion.

16) A method as claimed in claim 8, wherein the physical parameters are selected from among the following parameters upon valve closing: a pressure in the combustion chamber (PIVC), a temperature in the combustion chamber (TIVC), a ratio (XIVC) between a burnt gas mass and a total gas mass in the combustion chamber and an air mass (MIVC) in the cylinder and closure angle (θivc) of an intake valve.

17) A method as claimed in claim 9, wherein the physical parameters are selected from among the following parameters upon valve closing: a pressure in the combustion chamber (PIVC), a temperature in the combustion chamber (TIVC), a ratio (XIVC) between a burnt gas mass and a total gas mass in the combustion chamber and an air mass (MIVC) in the cylinder and closure angle (θivc) of an intake valve.

18) A method as claimed in claim 10, wherein the physical parameters are selected from among the following parameters upon valve closing: a pressure in the combustion chamber (PIVC), a temperature in the combustion chamber (TIVC), a ratio (XIVC) between a burnt gas mass and a total gas mass in the combustion chamber and an air mass (MIVC) in the cylinder and closure angle (θivc) of an intake valve.

19) A method as claimed in claim 11, wherein the physical parameters are selected from among the following parameters upon valve closing: a pressure in the combustion chamber (PIVC), a temperature in the combustion chamber (TIVC), a ratio (XIVC) between a burnt gas mass and a total gas mass in the combustion chamber and an air mass (MIVC) in the cylinder and closure angle (θivc) of an intake valve.

20) A method as claimed in claim 12, wherein the physical parameters are selected from among the following parameters upon valve closing: a pressure in the combustion chamber (PIVC), a temperature in the combustion chamber (TIVC), a ratio (XIVC) between a burnt gas mass and a total gas mass in the combustion chamber and an air mass (MIVC) in the cylinder and closure angle (θivc) of an intake valve.

21) A method as claimed in claim 13, wherein the physical parameters are selected from among the following parameters upon valve closing: a pressure in the combustion chamber (PIVC), a temperature in the combustion chamber (TIVC), a ratio (XIVC) between a burnt gas mass and a total gas mass in the combustion chamber and an air mass (MIVC) in the cylinder and closure angle (θivc) of an intake valve.

22) A method as claimed in claim 14, wherein the physical parameters are selected from among the following parameters upon valve closing: a pressure in the combustion chamber (PIVC), a temperature in the combustion chamber (TIVC), a ratio (XIVC) between a burnt gas mass and a total gas mass in the combustion chamber and an air mass (MIVC) in the cylinder and closure angle (θivc) of an intake valve.

23) A method as claimed in claim 15, wherein the physical parameters are selected from among the following parameters upon valve closing: a pressure in the combustion chamber (PIVC), a temperature in the combustion chamber (TIVC), a ratio (XIVC) between a burnt gas mass and a total gas mass in the combustion chamber and an air mass (MIVC) in the cylinder and closure angle (θivc) of an intake valve.

24) A method as claimed in claim 8, wherein a mass of fuel injected into the combustion chamber before the physical parameters reach their setpoint values is controlled by controlling combustion richness.

25) A method as claimed in claim 9, wherein a mass of fuel injected into the combustion chamber before the physical parameters reach their setpoint values is controlled by controlling combustion richness.

26) A method as claimed in claim 10, wherein a mass of fuel injected into the combustion chamber before the physical parameters reach their setpoint values is controlled by controlling combustion richness.

27) A method as claimed in claim 11, wherein a mass of fuel injected into the combustion chamber before the physical parameters reach their setpoint values is controlled by controlling combustion richness.

28) A method as claimed in claim 12, wherein a mass of fuel injected into the combustion chamber before the physical parameters reach their setpoint values is controlled by controlling combustion richness.

29) A method as claimed in claim 13, wherein a mass of fuel injected into the combustion chamber before the physical parameters reach their setpoint values is controlled by controlling combustion richness.

30) A method as claimed in claim 14, wherein a mass of fuel injected into the combustion chamber before the physical parameters reach their setpoint values is controlled by controlling combustion richness.

31) A method as claimed in claim 15, wherein a mass of fuel injected into the combustion chamber before the physical parameters reach their setpoint values is controlled by controlling combustion richness.

32) A method as claimed in claim 16, wherein a mass of fuel injected into the combustion chamber before the physical parameters reach their setpoint values is controlled by controlling combustion richness.

33) A method as claimed in claim 17, wherein a mass of fuel injected into the combustion chamber before the physical parameters reach their setpoint values is controlled by controlling combustion richness.

34) A method as claimed in claim 18, wherein a mass of fuel injected into the combustion chamber before the physical parameters reach their setpoint values is controlled by controlling combustion richness.

35) A method as claimed in claim 19, wherein a mass of fuel injected into the combustion chamber before the physical parameters reach their setpoint values is controlled by controlling combustion richness.

36) A method as claimed in claim 20, wherein a mass of fuel injected into the combustion chamber before the physical parameters reach their setpoint values is controlled by controlling combustion richness.

37) A method as claimed in claim 21, wherein a mass of fuel injected into the combustion chamber before the physical parameters reach their setpoint values is controlled by controlling combustion richness.

38) A method as claimed in claim 22, wherein a mass of fuel injected into the combustion chamber before the physical parameters reach their setpoint values is controlled by controlling combustion richness.

39) A method as claimed in claim 23, wherein a mass of fuel injected into the combustion chamber before the physical parameters reach their setpoint values is controlled by controlling combustion richness.