GB2481364A - A method for estimating the combustion of a quantity of fuel injected in a combustion chamber of a compression ignition internal combustion engine - Google Patents

Abstract

Disclosed is a method for estimating the combustion of a quantity of fuel injected in a combustion chamber of a compression ignition internal combustion engine during an engine cycle, wherein said quantity of fuel is injected through at least an injection pulse i, and wherein the method comprises the step of cal­culating the heat release rate dQi(t ) due to said at least an injection pulse i according to the following equation: dQ i (t) = k i •[Q FUELi (t-τ i ) -Q(t)] wherein ki is an index representing the speed of the combustion process during said at least an injection pulse i, Q(t) is a total heat released up to the time t, QFUELi(t-τi) is the energy of the fuel injected at the time (t-τi) by said at least an injection pulse i, τi is an index representing the delay between the time at which a certain fuel quantity of said at least an injection pulse i is injected and the time at which it actually burns.

Description

METHOD FOR ESTI!ING A CCtv1BUSTIC PROCESS IN A C4PRESSION I(JITION INTERNIL C4BUSTION ENGINE

TECHNICAL FIELD

The present invention relates to a method for estimating the combus- tion of a quantity of fuel injected in a combustion chamber of a corn- pression ignition internal combustion engine of a motor vehicle, typ-ically a Diesel engine, during an engine cycle.

BAOU

In current Diesel engines, the fuel is injected in the combustion chamber by means of a plurality of injection pulses per engine cycle according to a rr&ilti-injection pattern, typically by means of at least a pilot injection pulse and a following main injection pulse.

The combustion phasing is controlled acting on the start of injection (SOl) of said injection pulses. Current injectors have reached a very good quality and permit a high repeatability of the needle opening actuation across production dispersion and mileage. Nonetheless, the start of injection is not univocally related to the combustion phas- ing, indeed according to engine operating conditions and environmen-tal conditions, even with the same start of injection (SQl) and the same injection flow, the combustion could develop with very different speed and characteristics. This variability is further increased by the fact that modern Diesel engines are using increasingly higher rate of exhaust gas recirculation, which further increase the time delay between the start of injection and the start of combustion.

To monitor the combustion phasing, cylinder pressure sensors have been added to the engine control architecture. These sensors allow to measure the develoent of the combustion process and to calculate some more relevant indexes based on the heat released by the combus-tion. The drawback of this approach is in the very high cost of this sensors and in the fact that they give information on the combustion phasing after the combustion has happened. Therefore they cannot be used for feed-forward control of the injection and combustion but on-ly for closed loop corrections.

In order to overcome this drawback, it is advisable to estimate one or more relevant indexes of the combustion phasing by means of exist-ing information (no need to add cylinder pressure sensors), prior than combustion happens.

Several models exist in literature to predict the combustion event and the heat release, by way of example from Rolf Egnell "A Simple Approach to Study the Relation Between Fuel Rate, Heat Release Rate and NOx-Formation in Diesel Engines", SAE Technical Paper 1999-01- 3548, published on Octobe 1999, pages 4-6; or from Fabrizio Ponti, Enrico Corti, Gabriele Serra and Matteo De Cesare "Common Rail Multi-Jet Diesel Engine Combustion Model Development for Control Purposes" SAE Technical Paper 2007-01-0383, published on April 2007, pages 2-4; or from Franz G. Chinela and Gerhard C. Orthaber "Rate of Heat Release Prediction for Direct Injection Diesel Engines Based on Purely Mixing Controlled Combustion", SAE Technical Paper 1999-01-0186, published on March 1999, pages 2-5; or from Claes Ericson, Bjorn Westerberg, Magnus Andersson and Rolf Egnell "Modelling Diesel Engine Combustion and NOx Formation for Model Based Control and Simulation of Engine and Exhaust Aftertreatment Systems", SAE Technical Paper 2006-01- 0687, published on April 2006, pages 4-5.

The major limitations of these existing models are either in a too high complexity, leading to a not feasible implementation in real time control system applications, or in a too rough simplification, leading to an inaccuracy of the results outside very strict applica-tion boundaries.

In view of the above, it is an object of the present invention to provide a model based combustion estimation strategy capable to esti- mate the combustion dynamic with enough accuracy and acceptable com- putational requirements, such as to be effectively usable in a con-trol system of an engine.

Another object of the present invention is to achieve the above men- tioned purpose with a simple, rational and rather inexpensive solu-tion.

DISCLOSURE

These and/or other objects are attained by the characteristics of the embodiments of the invention as reported in independent claims. The dependent claims recite preferred and/or especially advantageous fea-tures of the embodiments of the invention.

Pn embodiment of the invention provides a method for estimating the combustion of a quantity of fuel injected in a combustion chamber of a compression ignition internal combustion engine during an engine cycle, wherein said quantity of fuel is injected through at least an injection pulse i, and wherein the method comprises the step of cal- culating the heat release rate dQ1(t), due to said at least an injec-tion pulse i, according to the following equation: dQ1(t) =k1 1QI/EL,I(t-rI)-Q(t)] wherein k1 is an index representing the speed of the combustion process during said at least an injection pulse i, Q(t) is a total heat released up to the time t, Q,1(t-r1) is the energy of the fuel injected at the time t-r1 by said at least an injection pulse i, r is an index representing the delay between the time at which a certain fuel quantity of said at least an injection pulse i is injected and the time at which it actually bums.

Taking into account the ignition delay between the time at which a certain fuel quantity is injected and the time at which it actually burns, this new strategy allows a reliable estimation of the heat re- lease rate with the aid of a simple equation, which requires a li- mited computational and memory efforts, such that it can be advanta-geously used for engine control purposes.

According to an embodiment of the invention, the energy of the in-jected fuel at the time t-r is calculated according to the following equation: QFUEL,I (t -= H fin, (t -soil wherein H0 is the lower heating value of the fuel, SOL. is the start of injection and rn is the fuel injection rate of said at least an injection pulse i.

This embodiment has the advantage of providing a simple way to calcu- late the energy of the injected fuel, which further reduces the can-putational and memory efforts.

The index r can be determined implementing various strategies that are currently available in literature, by way of example from H. Kim and N. Sung "Ccznbustion and Emission Modeling for a Direct Injection Diesel Engine", SAE Technical paper 2004-01-0104, published on 2004, page 5.

Nonetheless, a preferred embodiment of the invention provides for de- termining the index r1 on the base of at least a parameter representa- tive of the energy released by the combustion and a parameter repre-sentative of the air density in the combustion chamber, this latter being preferably considered at the time in which the piston is at the Top Death Center (TDC).

This embodiment of the invention provides a feasible strategy for de-termining the index r, because it has been found that an effective correlation exists between the heat released by the combustion, the air density in the combustion chamber and the ignition delay of an injected fuel quantity.

According to an aspect of this embodiment, the index r is thus deter-mined through an empirically determined relationship correlating the index r2 to at least the above mentioned combustion energy and air density related parameters.

This relationship can be advantageously determined through an experi- mental activity and then implemented as a subroutine, thereby allow-ing the determination of the index r with a minimum of computational effort.

Since different kinds of injection pulse, such as for example pilot injection pulse or main injection pulse, may be characterized by a different ignition delays, also the empirical relationship between r, combustion energy and air density could be different.

As a consequence, the above mentioned empirical relationship is pre-ferably determined for each injection pulse i independently, with the aid of a dedicated experimental activity focused on said specific in-jection pulse i, thereby allowing a more accurate estimation of the heat release rate caused by the fuel injection.

In order to simplify the experimental activity and reduce the compu-tational and memory efforts, the index r can determined through an empirically determined relationship correlating the index r1 to at least a single factor, which is a function of both the above men-tioned caribustion energy and air density related parameters.

According to an embodiment of the invention, the index k is deter-mined on the base of at least a parameter representing the engine speed and a parameter representing the quantity of fuel to be in-jected by said at least an injection pulse i.

This embodiment of the invention provides a reliable strategy for de-termining the index k1, because it has been experimentally found that an effective correlation exists between combustion speed, engine speed and fuel injected quantity.

According to an aspect of this embodiment, the index k1 is thus deter-mined through an empirically determined relationship correlating the index k1 to at least the above mentioned engine speed and injected fuel quantity related parameters.

This relationship can be advantageously determined through an experi- mental activity and then implemented in a subroutine, thereby allow-ing the determination of the index k1 with a minimum of computational effort.

Also in this case, the empirical relationship between k, engine speed and fuel injected quantity could be different for different kinds of injection pulse, so that this empirical relationship is preferably determined for each injection pulse i independently, by means of a dedicated experimental activity specially focused on said specific injection pulse i.

According to an embodiment of the invention, the method comprises the further step of time integrating the calculated heat release rate dQ1(t), so as to calculate a total heat Q released by said at least an injection pulse i.

In this way, the method advantageously allows a reliable estimation not only of the heat release rate, but also of the heat that is to-tally released during a combustion process.

According to another embodiment of the invention, when the fuel is injected through a plurality of injection pulses, the method provides for calculating the heat release rate dQ1(t) of each injection pulse separately.

Thanks to this solution, it has been found that the estimation of the heat release rate during the combustion process results more accu-rate.

According to an aspect of this embodiment, the method comprises the step of adding the heat release rates dQ1 (t) calculated for each in-jection pulse, so as to calculate a cumulative heat release rate dQ(t).

In this way, the cumulative heat release rate results a reliable in- dex of the combustion process, which takes into account the contribu-tions of all injection pulses.

According to another aspect of this embodiment, the method comprises the further step of time integrating the calculated cumulative heat release rate dQ(t), so as to calculate a cumulative heat Q totally released by said quantity of injected fuel.

In this way, the cumulative heat results still another reliable index of the combustion process, which represents the heat that is totally released during the combustion process with the contributions of all injection pulses.

The method according to the invention can be carried out with the help of a computer program comprising a program-code for carrying out all the steps of the method described above, and in the form of a computer program product cariprising the computer program.

The computer program product can be embodied as an internal cornbus-tion engine comprising an engine control unit (ECU), a data carrier associated to the ECU, and the computer program stored in the data carrier, so that, when the ECU executes the computer program, all the steps of the method described above are carried out.

The method can be also embodied as an electromagnetic signal, said signal being modulated to carry a sequence of data bits which represent a computer program to carry out all steps of the method.

BRIEF DESCRIPTION OF THE DRAWINGS

The present invention will now be described, by way of example, with reference to the accanpanying drawing.

Figure 1 is a flowchart representing a combustion estimating method according to an embodiment of the invention.

DETAILED DESCRIPTION

An embodiment of the present invention is hereinafter disclosed with reference to a conventional Diesel engine of a motor vehicle (not shown), which generally comprises a plurality of cylinders that indi-vidually defines a ccxnbustion chamber.

The combustion chamber is closed at the bottom by a reciprocating piston, whose movement is caused by the combustion process of a quan-tity of fuel that is injected into the combustion chamber during each engine cycle.

The present embodiment of the invention provides a method for esti-mating the dynamic of said combustion process.

During each engine cycle, the fuel is conventionally injected into the combustion chamber by means of a plurality of injection pulses that follows one another, according to a predetermined multi-injection pattern.

The fuel quantity that actually burns into the combustion chamber is usually injected by means of at least a pilot injection pulse and a following main injection pulse.

However, the present embodiment of the invention does not exclude that the total burning fuel quantity may be injected through a single injection pulse, or through more than two injection pulses of the multi-injection pattern.

Accordingly, the generic injection pulse of burning fuel is hereinaf-ter indicated with i, wherein i can vary from 1 to n, wherein 1 represents the first injection pulse of burning fuel and n represents the last injection pulse of burning fuel, and wherein n can be equal or greater than one.

The combustion estimating method provides a first routine, represented by the block 10 of figure 1, for separately calculating the heat release rate dQ1(t) caused by each injection pulse i at the generic time t.

Each heat release rate dQ (t) is calculated according to the following equation: dQ1(t) = k. 1QIEl,I(t-t)-Q(1)i wherein k1 is an index representing the speed of the combustion process calibrated for the 1th injection pulse, Q(t) is the total heat released up to the time t, Q,(t-r1) is the energy of the fuel in-jected by the th injection pulse at the time t-r1, r is an index representing the delay between the time at which a certain fuel quan-tity is injected by the 1th injection pulse and the time at which it actually burns.

The energy Q,(t-r) can be calculated according to the following equation: QIJIL,1 (t -t,) = H fin, (1 -r, )dt so/i wherein I-la is the lower heating value of the fuel, SOI is the start of injection of the 1th injection pulse, and m1 is the fuel injection rate of the 1th injection pulse.

In the present embodiment, each index r is determined on the base of a single factor, also referred as Waiting Factor WE', which is a func-tion of two key parameters that affect the ignition delay, namely a parameter representing the energy released by the combustion during an engine cycle and the air density p in the combustion chamber when the piston is at the TDC.

The air density p can be estimated according to the following equa-tion: p TD( RTrn IVTDC wherein P is the pressure of the air in the intake manifold of the engine, R is the specific gas constant, T is the temperature of the air in the intake manifold, yam is the combustion chamber volume when the piston is at the Bottom Death Center, and V is the combustion chamber volume when the piston is at the TDC.

The combustion energy related parameter can be either an Indicated Mean Effective Pressure (IMEP) or equivalently a fuel injected quan-tity qjnj,j to be injected by the th injection pulse.

In view of the above, WE" can be calculated according to the following equation: WF=C.1pT+D.1IMEP 1\SPTDCO) IMEP0 or equivalently a ( \h WF=C*1 PTDC +D.l P) q,10 wherein a, b, C and D are empirical constant dependent by the specif-ic engine; p-r, II'P0 and q,jo are respectively the air density at the TDC, the indicated mean effective pressure and the fuel injected quantity in a reference condition 0.

By way of example, for a turbocharged CIDI (Compression Ignition Di- rect Injection) engine having a displacement of 2 liters and a Com-pression Ratio of 15.5, the following parameter can be assumed: a=2, b=2, C=O.5, D=0.3, p=35kg/m3, Ifr.'P0=25bar and q1,10=80mrrrn3.

Returning now to the combustion estimating method, each index r is determined through a dedicated subroutine that correlates the index r to the waiting factor WE'.

For the sake of simplicity, these dedicated subroutines are globally represented by a single block 11 in figure 1.

The relationship between WF and r1 implemented in each subroutine is empirically determined through an experimental activity focused on the correspondent th injection pulse.

This experimental activity generally provides for operating a test Diesel engine, having the same characteristic of that to which the combustion estimating method is directed, in a plurality of different conditions, mainly at a plurality of different values of F (i.e. at a plurality of different values of the combustion energy related pa-rameter and of the air density); for measuring the heat release rate dQ1 due to the th injection pulse in each condition with the aid of a proper sensor (i.e. a cylinder pressure sensor); for calculating the index r1 in each condition; and for interpolating the calculated in-dexes r, in order to determine the empirical relationship between WF and r, which is finally implemented in the subroutine.

As a matter of fact, this experimental activity must be performed for each kind of injection pulse i separately, by way of example for both the pilot injection pulse and the main injection pulse, in order to obtain a plurality of subroutine, each of which implements a respec-tive relationship between WF and the index r, which is specifically calibrated for the correspondent injection pulse i.

By way of example, an experimental activity of this kind, performed for the main injection pulse of the above mentioned turbocharged CIDI (Compression Ignition Direct Injection) engine having a displacement of 2 liters and a Compression Ratio of 15.5, has returned the follow-ing empirical relationship: = >KJ WF if WF<0.6 and t*ftZU1fl[,U]=12O if WF»=0.6 wherein is the index referred to the main injection pulse, Ka=1801.3, I<i=16443, K=68524, 1<3=153979, K4=194421, 1<5=129332 and K5=35118.

In the present embodiment, each index k1 is determined on the base of the engine speed co and the quantity of fuel qinj,i to be injected by the th injection pulse.

In greater detail, each index k1 is determined through a dedicated subroutine, which is calibrated for the correspondent 1th injection pulse, and which correlates the index k to the engine speed co and the quantity of fuel qinj,i.

For the sake of simplicity, these dedicated subroutines are globally represented by a single block 12 in figure 1.

Each dedicated subroutine implements a relationship between co, qj,,j and k1, which is empirically determined through an experimental activ-ity focused on the correspondent 1th injection pulse.

This experimental activity generally comprises the steps of: operat-ing a test Diesel engine, having the same characteristic of that to which the combustion estimating method is directed, in a plurality of different conditions, mainly at a plurality of different values of the quantity of fuel qj,j and of the engine speed cdi; measuring the heat release rate dQ due to the i injection pulse in each condition with the aid of a proper sensor (i.e. a cylinder pressure sensor); calculating the index k in each condition; and interpolating the cal-culated indexes k, in order to determine the empirical relationship between u, qjnj,i and k, which is finally implemented in the dedicated subroutine.

As a matter of fact, this experimental activity must be performed for each kind of injection pulse i separately, by way of example for both the pilot injection pulse and the main injection pulse, in order to obtain a plurality of subroutines, each of which implements a rela-tionship between w, qjnj,j and k1 that is specifically calibrated on the correspondent injection pulse 1.

By way of example, an experimental activity, performed for the main injection pulse of the above mentioned turbocharged CIDI (Ccmpression Ignition Direct Injection) engine having a displacement of 2 liters and a Compression Ratio of 15.5, has returned the following empirical relationship: km,n =O.4819*ai+848.45 if q,1<20 mm3; = O.2992w+532.17 if rrm3; kmain (O.48l9*w+848.45)_(O.2992*w+532.l7).(q -20)+04819w+84845 if 20 rrm3<q, <40 rrrn3.

After the individual heat release rate dQ(t) have been calculated, the combustion estimating method provides a second routine, represented by the block 13 of figure 1, for calculating a cumulative heat release rate dQ(t) caused by all the injection pulses at the ge-neric time t.

The cumulative heat release rate dQ(t) is calculated according to the following equation: dQ(t)=dQ,(t) wherein dQ(t) is the heat release rate caused by each injection pulse i at the time t.

After the cumulative heat release rate dQ(t) has been calculated, the combustion estimating method provides a third routine, represented by the block 14 in figure 1, for calculating a cumulative heat Q global-ly released by the injected burning fuel.

The cumulative heat Q is calculated by time integrating the cumula-tive heat release rate dQ(t) according to the following equation: Q = jdQ(t)dt or equivalently according to the following equation: Q=Qi wherein Q is the total heat released by the th injection pulse and is calculated by time integrating dQ1(t) as follow: Q, JdQ,(r)dt.

While at least one exemplary embodiment has been presented in the foregoing summary and detailed description, it should be appreciated that a vast number of variations exist. It should also be appreciated that the exemplary embodiment or exemplary embodiments are only exam- ples, and are not intended to limit the scope, applicability, or con- figuration in any way. Rather, the forgoing summary and detailed de-scription will provide those skilled in the art with a convenient road map for implementing at least one exemplary embodiment, it being understood that various changes may be made in the function and ar-rangement of elements described in an exemplary embodirrnt without departing from the scope as set forth in the appended claims and in their legal equivalents.

REFERENCE NUMBERS

First routine ii ii Subroutine 12 Ki subroutine 13 Second routine 14 Third routine ti Ignition delay index Engine speed related parameter qinj,i Injected fuel quantity ki Combustion speed index dQi(t) Heat release rate dQ(t) Cumulative heat release rate Q Cumulative heat released WF Waiting factor CLD4S

Claims (15)

1. Method for estimating the combustion of a quantity of fuel in-jected in a combustion chamber of a compression ignition internal combustion engine during an engine cycle, wherein said quantity of fuel is injected through at least an injection pulse i, and wherein the method comprises the step of calculating the heat release rate dQ1(t) due to said at least an injection pulse i according to the fol-lowing equation: dQ1(t) Ic [QEL,I(t --Q(t)j wherein k is an index representing the speed of the combustion process during said at least an injection pulse i, Q(t) is a total heat released up to the time t, Q,1(t-r) is the energy of the fuel injected at the time -r by said at least an injection pulse i, r is an index representing the delay between the time at which a certain fuel quantity of said at least an injection pulse i is injected and the time at which it actually bums.

2. Method according to claim 1, wherein the energy of the injected fuel at the time t-r1 is calculated according to the following equa-tion: QI;UEL,i (t -= H Jrn1 (t -r1)dt soli wherein H is the lower heating value of the fuel, SOI. is the start of injection and rn is the fuel injection rate of said at least an injection pulse i.

3. Method according to claim 1, wherein the index r is determined on the base of at least a parameter (IMEP, qjnj,i) representing the energy released by the combustion and a parameter (p) representing the air density in the combustion chamber.

4. Method according to claim 3, wherein the index r Is determined through an empirically determined relationship correlating the index r to at least the combustion energy related parameter (IMEP, qj,j) and air density related parameter (p).

5. Method according to claim 3, wherein the index r is determined through an empirically determined relationship correlating the index r to at least a single factor (WE), which is a function of both the combustion energy related parameter (llP, qinj,i) and air density re-lated parameter (pc).

6. Method according to claim 1, wherein the index k is determined on the base of at least of a parameter representing the engine speed (w) and a parameter representing the quantity of fuel (qj,j) to be injected by said at least an injection pulse i.

7. Method according to claim 5, wherein the index k is determined through an empirically determined relationship correlating the index k to at least the engine speed related parameter (w) and injected fuel quantity related parameter (qthj,j).

8. Method according to claim 1, comprising the step of time inte-grating the calculated heat release rate dQ(t), so as to calculate a total heat Q released by said at least an injection pulse i.

9. Method according to claim 1, wherein said quantity of fuel is in- jected through a plurality of injection pulses, and wherein the me- thod provides for calculating the heat release rate dQ1 (t) of each in-jection pulse separately.

10. Method according to claim 9, comprising the step of adding the heat release rates dQ1(t) calculated for each injection pulse, so as to calculate a cumulative heat release rate dQ(t).

11. Method according to claim 10, comprising the step of time inte-grating the calculated cumulative heat release rate dQ(t), so as to calculate a cumulative heat Q totally released by said quantity of injected fuel.

12. Computer program comprising a computer-code for carrying out a method according to any of the preceding claims.

13. Computer program product on which the computer program according to claim 12 is stored.

14. Compression ignition internal combustion engine caiiprising an ECU, a data carrier (14) associated to the ECU, and a computer pro-gram according to claim 12 stored in the data carrier.

15. An electromagnetic signal modulated as a carrier for a sequence of data bits representing the computer program according to claim 12.Amendment to the claims have been filed as follows 1. Method of controlling the combustion of a quantity of fuel in-jected in a combustion chamber of a compression ignition internal combustion engine during an engine cycle, wherein said quantity of fuel is injected through at least an injection pulse i, and wherein the method comprises the following steps: a) calculating, prior to a combustion of injected fuel, the heat release rate dQ1(t) due to said at least an injection pulse i according to the following equation: dQ,(t) = k. [QFUEL, (t -z-,) -Q(t)j wherein k1 is an index representing the speed of the combustion process during said at least an injection pulse i, Q(t) is a to-tel heat released up to the time t, Q,(t-r1) is the energy of the fuel injected at the time t-r by said at least an injection pulse i, r is an index representing the delay between the time at which a certain fuel quantity of said at least an injection :.: pulse i is injected and the time at which it actually burns, b) controlling the start of injection with said calculated heat re-lease rate.2. Method according to claim 1, wherein the energy of the injected fuel at the time t-r is calculated according to the following equa-tion: QI,lL.I(t -r,) = H, fin1 (t -r1)dt SO!, wherein H is the lower heating value of the fuel, SOIL is the start of injection and m is the fuel injection rate of said at least an injection pulse i.3. Method according to claim 1, wherein the index r1 is determined on the base of at least a parameter (IIVIEP, qinj,i) representing the energy released by the combustion and a parameter (p) representing the air density in the combustion chamber.4. Method according to claim 3, wherein the index r1 is determined through an empirically determined relationship correlating the index r to at least the combustion energy related parameter (IME'P, qjnj,j) and air density related parameter (p).5. Method according to claim 3, wherein the index r. is determined through an empirically determined relationship correlating the index S.....* r to at least a single factor (WE), which is a function of both the ***.*.* combustion energy related parameter (IMEP, qj,j) and air density re-n..' 15 lated parameter (p).6. Method according to claim 1, wherein the index Ic1 is determined on the base of at least of a parameter representing the engine speed (&) and a parameter representing the quantity of fuel (q,j) to be injected by said at least an injection pulse 1.7. Method according to claim 5, wherein the index k is determined through an empirically determined relationship correlating the index k1 to at least the engine speed related parameter (w) and injected fuel quantity related parameter (qj,j).8. Method according to claim 1, comprising the step of time inte-grating the calculated heat release rate dQ(t), so as to calculate a total heat Q1 released by said at least an injection pulse i.9. Method according to claim 1, wherein said quantity of fuel is in- jected through a plurality of injection pulses, and wherein the me- thod provides for calculating the heat release rate dQ1(t) of each in-jection pulse separately.10. Method according to claim 9, comprising the step of adding the heat release rates dQ(t) calculated for each injection pulse, so as to calculate a cumulative heat release rate dQ(t).U. Method according to claim 10, comprising the step of time inte-grating the calculated cumulative heat release rate dQ(t), so as to calculate a cumulative heat Q totally released by said quantity of injected fuel. ** 12. Computer program comprising a computer-code for carrying out a S.....* method according to any of the preceding claims.n... 15 13. Computer program product on which the computer program according to claim 12 is stored.14. Compression ignition internal combustion engine comprising an ECU, a data carrier (14) associated to the ECU, and a computer pro-gram according to claim 12 stored in the data carrier.15. An electromagnetic signal modulated as a carrier for a sequence of data bits representing the computer program according to claim 12.*::r: INTELLECTUAL 25 . ... PROPERTY OFFICE Application No: GB 1003856.0 Examiner: Mr Alastair Kelly Claims searched: 1-15 Date of search: 19 October 2011 Patents Act 1977: Search Report under Section 17 Documents considered to be relevant: Category Relevant Identity of document and passage or figure of particular relevance to claims X,E 1-15 GB2478647A [GM] See whole document A -W02010/100965A1 [HONDA] See abstract and figures A -U52006/102147A1 [NORIMOTO] See abstract and figures A -U52009/266336A1 [MORIMOTO] See abstract and figures A -JP2007309309 A [HONDA] See figures and note EPODOC & WPI abstracts attached A -JP11148410A [ISUZU] See figures and note EPODOC & WPI abstracts attached A -JP11125141A [ISUZU] See figures and note EPODOC & WPI abstracts attached Categories: X Document indicating lack of novelty or inventive A Document indicating technological background and/or state step of the art.Y Document indicating lack of inventive step if P Document published on or after the declared priority date but combined with one or more other documents of before the filing date of this invention.same category.& Member of the same patent family E Patent document published on or after, but with priority date earlier than, the filing date of this application.Field of Search:Search of GB, EP, WO & US patent documents classified in the following areas of the UKCX: Worldwide search of patent documents classified in the following areas of the IPC FO2D The following online and other databases have been used in the preparation of this search report EPODOC, WPI Intellectual Property Office is an operating name of the Patent Office www.ipo.gov.uk *.:r: INTELLECTUAL 26 . ... PROPERTY OFFICE International Classification: Subclass Subgroup Valid From None Intellectual Property Office is an operating name of the Patent Office www.ipo.gov.uk