# Patent

## Title
Method to determine the mass of air trapped in each cylinder of an internal combustion engine

## Abstract

A method to determine the mass of air trapped in each cylinder of an internal combustion engine, which comprises determining, based on a model using measured and/or estimated physical quantities, a value for a first group of reference quantities; determining, based on the model, the actual inner volume of each cylinder as a function of the speed of rotation of the internal combustion engine and of the closing delay angle of the intake valve; and calculating the mass of air trapped in each cylinder as a function of the first group of reference quantities and of the actual inner volume of each cylinder.

## Background
CROSS-REFERENCE TO RELATED APPLICATIONS
This Patent application claims priority from Italian Patent Application No. 102018000010164 filed on Nov. 8, 2018, the entire disclosure of which is incorporated herein by reference.
BACKGROUND OF THE INVENTION
1. Field of the Invention
The invention relates to a method to determine the mass of air trapped in each cylinder of an internal combustion engine.
2. Description of the Related Art
As it is known, an internal combustion engine supercharged using a turbocharger supercharging system comprises a number of injectors injecting fuel into respective cylinders, each connected to an intake manifold by at least one respective intake valve and to an exhaust manifold by at least one respective exhaust valve.
The intake manifold receives a gas mixture comprising both exhaust gases and fresh air, i.e. air coming from the outside through an intake duct, which is provided with an air filter for the fresh air flow and is regulated by a throttle valve. Along the intake duct, preferably downstream of the air filter, there is also provided an air flow meter.
The air flow meter is a sensor connected to an electronic control unit and designed to detect the flow rate of fresh air taken in by the internal combustion engine. The flow rate of fresh air taken in by the internal combustion engine is an extremely important parameter for the engine control, in particular to determine the quantity of fuel to be injected into the cylinders so as to obtain a given air/fuel ratio in an exhaust duct downstream of the exhaust manifold. However, the air flow meter typically is a very expensive and fairly delicate component as oil vapours and dust can dirty it, thus altering the reading of the value of the flow rate of fresh air taken in by the internal combustion engine.
SUMMARY OF THE INVENTION
The object of the invention is to provide a method to determine the mass of air trapped in each cylinder of an internal combustion engine, said method being easy and economic to be implemented.
According to the invention, there is provided a method to determine the mass of air trapped in each cylinder of an internal combustion engine as claimed in the appended claims. Other objects, features and advantages of the present invention will be readily appreciated as the same becomes better understood after reading the subsequent description taken in connection with the accompanying drawings.



BRIEF DESCRIPTION OF THE DRAWINGS
The invention will now be described with reference to the accompanying drawings, showing a non-limiting embodiment thereof, wherein:
FIG. 1 schematically shows a preferred embodiment of an internal combustion engine provided with an electronic control unit implementing a method according to the invention;
FIG. 2 shows, in detail, a cylinder of the engine of FIG. 1;
FIGS. 3 and 4 schematically show the overlap phase of an intake valve and of an exhaust valve of the engine of FIG. 1; and
FIG. 5 shows the development of the function β used in the method according to the invention.



DETAILED DESCRIPTION OF THE INVENTION
In FIGS. 1 and 2, number 1 indicates, as a whole, an internal combustion engine, preferably supercharged using a turbocharger supercharging system.
The internal combustion engine 1 comprises a number of injectors 2, which directly inject fuel into four cylinders 3 (preferably four cylinders arranged in line), each connected to an intake manifold 4 by at least one respective intake valve 5 (shown in FIG. 3) and to an exhaust manifold 6 by at least one respective exhaust valve 7 (shown in FIG. 2). For each cylinder 3 there is provided a corresponding injector 2; according to the embodiment shown in FIG. 2, the injection is an indirect injection and, therefore, each injector 2 is arranged upstream of the cylinder 3 in an intake duct 8 connecting the intake manifold 4 to the cylinder 3. According to an alternative embodiment which is not shown herein, the injection is a direct injection and, therefore, each injector 2 is partially arranged inside the cylinder 3 through the crown end of the cylinder 3.
According to FIG. 1, each cylinder 3 houses a respective piston 9, which is mechanically connected, by a connecting rod, to a drive shaft 10, so as to transmit to the drive shaft 10 itself, in a known manner, the force generated by the combustion inside the cylinder 3.
The intake manifold 4 receives a gas mixture comprising both exhaust gases (as described more in detail below) and fresh air, i.e. air coming from the outside through the intake duct 8, which is preferably provided with an air filter for the fresh air flow and is regulated by a throttle valve 12, which preferably is an electronically controlled valve and is movable between a closing position and a maximum opening position. Furthermore, no air flow meter is provided along the intake duct 8.
The intake valves 5 and/or the exhaust valves 7 are controlled with a VVT (variable valve timing) device, which hydraulically acts upon the shaft operating the intake valves 5 and or the exhaust valves 7, respectively, changing the inclination thereof relative to a drive shaft.
In particular, the position of each exhaust valve 7 is directly controlled by a cam shaft 13, which receives the motion of the drive shaft 10; similarly, the position of each intake valve 5 is directly controlled by a cam shaft 14, which receives the motion of the drive shaft 10.
Along the intake pipe 8 there is preferably arranged an intercooler, which fulfils the function of cooling the air taken in and is preferably built-in in the intake manifold 4. The exhaust manifold 6 is connected to an exhaust duct 18, which feeds the exhaust gases produced by the combustion to an exhaust system, which releases the gases produced by the combustion into the atmosphere and normally comprises at least one catalytic converter (if necessary, provided with a diesel particulate filter) and at least one silencer arranged downstream of the catalytic converter.
The supercharging system of the internal combustion engine 1 comprises a turbocharger provided with a turbine, which is arranged along the exhaust duct 18 so as to rotate at a high speed due to the action of the exhaust gases expelled from the cylinders 3, and a compressor, which is arranged along the intake duct 8 and is mechanically connected to the turbine so as to be caused to rotate by the turbine itself in order to increase the pressure of the air present in the feeding duct 8.
The description above explicitly refers to an internal combustion engine 1 supercharged by using a turbocharger. Alternatively, the control method described above can find advantageous application in any supercharged internal combustion engine, for example an engine supercharged by using a dynamic or volumetric compressor.
According to one variant, along the exhaust duct 18 there is provided a bypass duct, which is connected in parallel to the turbine so as to have its ends connected upstream and downstream of the turbine itself.
The internal combustion engine 1 advantageously comprises, furthermore, a high-pressure exhaust gas recirculation circuit EGRHP, which comprises, in turn, a bypass duct connected in parallel to the assembly consisting of the four cylinders 3, the intake manifold 4 and the exhaust manifold 6. Along the bypass duct there is provided an EGR valve, which is designed to adjust the flow rate of the exhaust gases flowing through the bypass duct and is controlled by an electric motor. Along the bypass duct, downstream of the EGR valve, there is provided a heat exchanger, which fulfils the function of cooling the gases flowing out of the exhaust manifold.
Alternatively, along the intake duct 8 there is provided a bypass duct, which is connected in parallel to the compressor so as to have its ends connected upstream and downstream of the compressor itself; along the bypass duct there is provided a valve Poff, which is designed to adjust the flow rate of the air flowing through the bypass duct and is controlled by an electric actuator.
The internal combustion engine 1 is controlled by an electronic control unit 30, which controls the operation of all the components of the internal combustion engine 1. In particular, the electronic control unit 30 is connected to sensors which measure the temperature To and the pressure Po along the intake duct 8 upstream of the compressor and to sensors which measure the temperature and the pressure along the intake duct 8 upstream of the throttle valve 12 as well as to a sensor 31 which measures the temperature and the pressure of the gas mixture present in the intake manifold 4. Furthermore, the electronic control unit 30 is connected to a sensor which measures the angular position (and, hence, the rotation speed) of the drive shaft 10 and to a sensor (typically a UHEGO or UEGO linear oxygen sensor—which is known and not described in detail) which measures the air/fuel ratio of the exhaust gases upstream of the catalytic converter and, finally, to a sensor which measures the stroke of the intake and/or exhaust valves.
According to one variant, the internal combustion engine 1 finally comprises a low-pressure exhaust gas recirculation circuit EGRLP, which comprises, in turn, a bypass duct originating from the exhaust duct 18, preferably downstream of the catalytic converter, and leading into the intake duct 8, upstream of the compressor; the bypass duct is connected in parallel to the turbocharger. Along the bypass duct there is provided an EGR valve, which is designed to adjust the flow rate of the exhaust gases flowing through the bypass duct. Along the bypass duct, upstream of the EGR valve, there is also provided a heat exchanger, which fulfils the function of cooling the gases flowing out of the exhaust manifold 6 and into the compressor.
In the electronic control unit 30 there is stored a calculation model, which is used to determine, among other things, the mass m of air trapped in each cylinder 3 (for each cycle) and the mass MTOT of air taken in by the internal combustion engine 1.
The model includes a plurality of input parameters, among which there are: the number of revolutions (rpm), the pressure value in the intake manifold 4 and other side conditions (such as, for example, the temperature inside the intake manifold 4 and the temperature of the coolant fluid used in the supercharged internal combustion engine 1).
Since the VVT device varies the timing of the intake valves 5 and the timing of their crossing with the exhaust valves 7 (i.e. the phase during which the intake valve 5 and the exhaust valve 7 are simultaneously open), the model requires to know the following input parameters for each cylinder 3 as well; some parameters are illustrated schematically in FIG. 5 (with respect to the top dead center TDC and the bottom dead center BDC) where:
IVCref represents the reference closing angle of the intake valve 5;
IVOref represents the reference opening angle of the intake valve 5;
EVCref represents the reference closing angle of the exhaust valve 7;
EVOref represents the reference opening angle of the exhaust valve 7;
IVC represents the closing advance angle of the intake valve 5;
IVO represents the opening advance angle of the intake valve 5;
EVC represents the closing advance angle of the exhaust valve 7; and
EVO represents the opening advance angle of the exhaust valve 7.
Through the input parameters listed above, the following quantities are defined:

VVTI=IVC−IVCref=IVO−IVOref  [1]

VVTE=EVO−EVOref=EVC−EVCref  [2]

VVTI represents the angular extent of the opening or closing difference relative to the reference values concerning the intake valve 5; and
VVTE represents the angular extent of the opening or closing difference relative to the reference values concerning the exhaust valve 7.
In order to determine the mass m of air trapped in each cylinder 3 for each cycle, the model uses the ideal gas law (known from the literature), according to which

m=(P*V)/(R*T)  [3]

where:

P represents the mean of the pressure for the engine cycle inside the intake manifold 4;
T represents the temperature of the mixture of fresh air and/or exhaust gases inside the intake manifold 4;
R represents the constant of the mixture of fresh air and/or exhaust gases; and
V represents the inner volume of the cylinder 3, when the respective intake valve 5 and the respective exhaust valve 7 are closed).
The ideal gas law [3] was experimentally adjusted for the model by incorporating the constant R of the mixture of fresh air and/or exhaust gases, so that the mass m of air trapped in each cylinder 3 for each cycle is expressed as follows:

m=P*V*f1(T,P)*f2(TH2O,P)  [4]

wherein TH2O is the temperature of the internal combustion engine 1 (preferably expressed through the temperature of the coolant liquid of the internal combustion engine 1).

Parameters P, V, T, on the other hand, have the meaning described above for formula [3].
Finally, the ideal gas law [4] was further experimentally adjusted for the filling model so that the mass m of air trapped in each cylinder 3 for each cycle takes into account the gases produced by the combustion in the previous work cycle and present inside the cylinder 3 (because they did not flow out of the cylinder 3 or because they were re-sucked into the cylinder 3):

m=(P*V−OFF)*f1(T,P)*f2(TH2O,P)  [5]

wherein OFF is the variable (mass) taking into account the gases produced by the combustion in the previous work cycle and present inside the cylinder 3 (because they did not flow out of the cylinder 3 or because they were re-sucked into the cylinder 3).

Parameters P, V, T, again, have the meaning described above for formula [3].
In reference conditions, in order to calibrate the model, the temperature TH2O of the internal combustion engine 1, namely the temperature of the coolant liquid of the internal combustion engine 1, is assumed to be equal to 90° C. and the temperature T is assumed to be equal to 40° C.
Functions f1 and f2 mentioned above are defined in an experimental phase through (2d) maps as a function, respectively, of the pressure P inside the intake manifold 4 and of the temperature T inside the intake manifold 4 for function f1 and of the pressure P inside the intake manifold 4 and of the temperature TH2O of the internal combustion engine 1 for function f2. It is evident that, in reference conditions (for example, the reference temperature inside the intake manifold 4 is equal to 25° C.), functions f1 and f2 have a unitary value.
The inner volume V of the cylinder is variable 3 (from a geometrical point of view) as a function of the closing advance angle IVC of the respective intake valve 5. Indeed, the actual inner volume V of the cylinder 3 results from the sum of the dead volume VCC of the combustion chamber of the cylinder 3 (i.e. the volume that is not scavenged by the respective piston 9) and of the volume Vc scavenged by the respective piston 9 until the closing of the respective intake valve 5 (i.e. of the angle of rotation of the crank relative to the top dead centre PMS).
Hereinafter you can find the kinematic law (known from the literature and not described in detail) used to calculate the inner volume V of the cylinder 3 in the area of the crank angle indicated with α:



 
  
   
    
                    
     
      
       
        V
        ⁡
        (
        α
        )
       
       =
       
        
         V
         
          CC
             
          +
         
        
        ⁢
           
        
         
          V
          C
         
         (
         α
         )
        
       
      
      ⁢
      

      
       
        V
        ⁡
        (
        α
        )
       
       =
       
        
         V
         CC
        
        +
        
         
          S
          *
         
         ⁢
         
          r
          *
         
         ⁢
         
          
          
           [
           
            
             
              (
              
               1
               +
               
                1
                λ
               
              
              )
             
             *
             
              
               1
               -
               
                δ
                
                 
                  (
                  
                   1
                   +
                   λ
                  
                  )
                 
                 2
                
               
              
             
            
            -
            
             cos
             ⁢
                
             α
            
            -
            
             
              1
              λ
             
             *
             
              
               1
               -
               
                
                 (
                 
                  
                   
                    λ
                    *
                   
                   ⁢
                   sen
                   ⁢
                   α
                  
                  -
                  δ
                 
                 )
                
                2
               
              
             
            
           
           ]
          
         
        
       
      
     
    
   
   
    
     [
     6
     ]
    
   
  
 



where:

V represents the inner volume of the cylinder 3;
VCC represents the dead volume of the combustion chamber of the cylinder 3;
α represents the angle of rotation of the crank relative to the top dead centre PMS;
r represents the crank radius;
S represents the surface area of the piston 9;
L represents the length of the connecting rod;
d represents the offset between the axis of the cylinder 3 and the rotation axis of the drive shaft 10;
λ represents the r/L ratio; and
δ represents the d/L ratio.
According to one variant, generally speaking, the inner volume V of the cylinder 3 is variable as a function of a geometrical factor represented by the closing advance angle IVC of the respective intake valve 5, by a dynamic factor represented by the speed n of rotation of the internal combustion engine 1 (or number of revolutions rpm) and by the pressure P measured for the engine cycle inside the intake manifold 4.
In particular, the law [6] to determine the inner volume V of the cylinder 3 was experimentally adjusted for the model by introducing the two functions fv and fp and is expressed as follows:

V=fV(IVC,n)*fP(P,n)  [7]

Parameters P, n, IVC have the meaning already discussed above.
Furthermore, it should be taken into account that, at the beginning of the intake stroke of any engine cycle, inside the cylinder 3 there are also the residual gases of the combustion of the previous engine cycle.
From a geometrical point of view, the volume occupied by the residual gases of the combustion of the previous engine cycle can be expressed through the sum of the dead volume VCC of the combustion chamber of the cylinder 3 and of a volume VC scavenged by the respective piston 9 inside the cylinder 3.
The volume VC scavenged by the piston 9 inside the cylinder 3 is variable as a function of the parameter TVC, which is better described below.
In particular, according to a first variant, the volume VC scavenged by the piston 9 inside the cylinder 3 corresponds to the volume scavenged by the piston 9 until the instant in which the respective exhaust valve 7 closes, in case the respective intake valve 5 opens following the closing of the respective exhaust valve 7.
According to a second variant, the volume VC scavenged by the piston 9 inside the cylinder 3 corresponds to the volume scavenged by the piston 9 until the instant in which the respective intake valve 5 opens, in case the respective exhaust valve 7 closes following the opening of the respective intake valve 5.
According to a third variant, the volume VC scavenged by the piston inside the cylinder 3 corresponds to the volume scavenged by the piston 9 up to the top dead centre PMS, in case the opening instant of the respective intake valve 5 is prior to said top dead centre PMS. It is evident that, in this case, the volume VC scavenged by the respective piston inside the cylinder 3 in zero and the inner volume V of the cylinder 3 corresponds to the dead volume VCC of the combustion chamber of the cylinder 3.
In other words, the parameter TVC can alternatively correspond to the closing advance angle EVC of the exhaust valve 7 or to the greatest value between zero and the smallest value between the closing advance angle EVC of the exhaust valve 7 and the opening advance angle IVO of the intake valve 5.
Since the VVT system changes the timing of the intake valves 5 and of their overlap with the exhaust valves 7, the model also allows for a determination of the mass flow rate flowing during the overlap phase between each intake valve 5 and the respective exhaust valve 7. In the description below, the term overlap defines the phase (time interval) in which each intake valve 5 and the respective exhaust valve 7 are simultaneously open.
According to what is schematically shown in FIG. 4, the following geometrical quantities are defined (relative to the top dead centre TDC and to the bottom dead centre BDC):
OVL represents the duration of the overlap phase comprised between the closing advance angle EVC of the exhaust valve 7 and the opening advance angle IVO of the intake valve 5;
G represents the centre of gravity of the overlap phase between each intake valve 5 and the respective exhaust valve 7; and
g represents the difference between the top dead centre PMS and the centre of gravity G.
Hereinafter the law (known from the literature and not described in detail) used to calculate the mass flow rate through a section of a duct (or through an orifice) can be determined. In this case, the law is used to calculate the mass MOVL flowing from the exhaust to the intake through the intake valve 5 and the exhaust valve 7:



 
  
   
    
     
      M
      OVL
     
     =
     
      
       C
       D
      
      *
      A
      *
      
       
        P
        0
       
       
        
         R
         
          T
          0
         
        
       
      
      *
      
       B
       ⁡
       (
       
        P
        
         P
         0
        
       
       )
      
     
    
   
   
    
     [
     8
     ]
    
   
  
 



where:

A represents the area of the passage section;
CD represents the discharge coefficient;
P represents the pressure downstream of the passage section;
P0 represents the pressure at the inlet of the passage section;
T0 represents the temperature at the inlet of the passage section;
R represents the constant of the fluid flowing in the passage section; and
B represents the flow compressibility function expressed by the following equation [8′]:



 
  
   
    
     B
     =
     
      
       
        2
        ⁢
        K
       
       
        K
        -
        1
       
      
      *
      
       
        
         
          (
          
           P
           
            P
            0
           
          
          )
         
         
          2
          K
         
        
        -
        
         
          (
          
           P
           
            P
            0
           
          
          )
         
         
          
           K
           +
           1
          
          K
         
        
       
      
     
    
   
   
    
     [
     
      8
      ′
     
     ]
    
   
  
 



wherein K represents the ratio between the specific heat Cp at constant pressure and the specific heat Cv at constant volume.

The law [8] is experimentally adjusted for the model by integrating it between the instant t1 in which the overlap phase begins and the instant t2 in which the overlap phase ends according to the equation [9] below:



 
  
   
    
     
      M
      OVL
     
     =
     
      
       
        
         P
         0
        
        
         R
         /
         
          T
          0
         
        
       
      
      *
      
       B
       ⁡
       (
       
        P
        
         P
         0
        
       
       )
      
      *
      
       
        ∫
        
         t
         1
        
        
         t
         2
        
       
       
        
         
          A
          IS
         
         (
         t
         )
        
        ⁢
        dt
       
      
     
    
   
   
    
     [
     9
     ]
    
   
  
 



If the variable dt is replaced with dθ/ω (wherein θ represents the engine angle and ω represents the speed of rotation of the internal combustion engine 1), the following equation [10] is obtained:



 
  
   
    
     
      M
      OVL
     
     =
     
      
       
        
         P
         0
        
        
         R
         /
         
          T
          0
         
        
       
      
      *
      
       B
       ⁡
       (
       
        P
        
         P
         0
        
       
       )
      
      *
      
       ∫
       
        
         
          A
          IS
         
         (
         θ
         )
        
        *
        
         1
         ω
        
        ⁢
        d
        ⁢
        θ
       
      
     
    
   
   
    
     [
     10
     ]
    
   
  
 



Finally, assuming that the speed ω of rotation of the internal combustion engine 1 is constant during the overlap phase, equation [10] can be simplified in the following equation [11]:



 
  
   
    
     
      M
      OVL
     
     =
     
      
       
        
         P
         0
        
        
         R
         /
         
          T
          0
         
        
       
      
      *
      
       B
       ⁡
       (
       
        P
        
         P
         0
        
       
       )
      
      *
      
       ∫
       
        
         
          A
          IS
         
         (
         θ
         )
        
        ⁢
        d
        ⁢
        θ
       
      
     
    
   
   
    
     [
     11
     ]
    
   
  
 



In the preceding equations AIS represents the isentropic area.
Inside the electronic control unit 30, equation [11] is further experimentally adjusted for the model so as to obtain the mass MOVL as follows:



 
  
   
    
     
      M
      OVL
     
     =
     
      
       S
       id
      
      *
      
       β
       ⁡
       (
       
        
         P
         
          P
          0
         
        
        ,
        n
       
       )
      
      *
      
       
        P
        0
       
       
        P
        
         0
         ⁢
         
          _
          ⁢
          REF
         
        
       
      
      *
      
       
        
         T
         
          0
          ⁢
          
           _
           ⁢
           REF
          
         
        
        
         T
         0
        
       
      
      *
      
       1
       n
      
     
    
   
   
    
     [
     12
     ]
    
   
  
 



where:

Sid represents the ideal section;
n represents the speed of the internal combustion engine (1);
P0_REF represents the reference pressure upstream of the passage section;
T0_REF represents the reference temperature upstream of the passage section;
T0 represents the temperature upstream of the passage section;
P0, P represents the pressure upstream and downstream, respectively, of the passage section; and
B represents the compression ratio.
The ideal section Sid of the passage is obtained from the product of two functions, wherein the first function A is experimentally determined through the (2d) map variable as a function of the speed n of the internal combustion engine 1 and of the parameter OVL, whereas the second function G is experimentally determined through a (2d) map variable as a function of the speed n of the internal combustion engine 1 and of the parameter g.
The combustion chamber of the cylinder 3 is considered to be a passage section (preferably upstream and downstream of the respective valves 5, 7). In case the intake pressure is greater than the exhaust pressure, the “upstream” pressure and temperature to be taken into account are the pressure and the temperature upstream of the intake valve 5 (and, hence, measured by the sensor present in the intake manifold 4); whereas the “downstream” pressure and temperature to be taken into account are the pressure and the temperature downstream of the exhaust valves 7 and, hence, the pressure and the temperature of the exhaust gases (typically obtained from a model or, if possible, measured using a dedicated sensor).
If the exhaust pressure is greater than the intake pressure, the reverse logic applies; namely, the “downstream” pressure and temperature to be taken into account are the pressure and the temperature upstream of the intake valve 5 (and, hence, measured by the sensor present in the intake manifold 4); whereas the “upstream” pressure and temperature to be taken into account are the pressure and the temperature downstream of the exhaust valves 7 and, hence, the pressure and the temperature of the exhaust gases (typically obtained from a model or, if possible, measured using a dedicated sensor).
In both cases, we are dealing with mean values over the engine cycle, namely over the 720° of rotation of the drive shaft 10.
In case the pressure in the exhaust manifold 6 is greater than the pressure in the intake manifold 4, a portion of the exhaust gases produced by the combustion flows from the combustion chamber towards the intake manifold 4; during the following combustion cycle, the exhaust gas portion will be then reintroduced into the combustion chamber through the intake valve 5. This operating mode is indicated as “inner EGR” and formula [12] is adjusted by replacing the downstream pressure P0 with the exhaust pressure PEXH and by replacing the downstream temperature T0 with the exhaust temperature TEXH. Therefore, in this case, the mass MOVL is expressed as follows:



 
  
   
    
     
      M
      OVL
     
     =
     
      
       S
       id
      
      *
      
       β
       ⁡
       (
       
        
         P
         
          P
          EXH
         
        
        ,
        n
       
       )
      
      *
      
       
        P
        EXH
       
       
        P
        
         0
         ⁢
         
          _
          ⁢
          REF
         
        
       
      
      *
      
       
        
         T
         
          0
          ⁢
          
           _
           ⁢
           REF
          
         
        
        
         T
         EXH
        
       
      
      *
      
       1
       n
      
     
    
   
   
    
     [
     13
     ]
    
   
  
 



The mass MEGRI of “inner EGR” can be expressed as follows:

MEGRI=MOVL+PEXH*VCC/(R*TEXH)  [14]

Quantities MOVL, PEXH, VCC, R and TEXH have the meaning already discussed above.
In case the pressure in the intake manifold 4 is greater than the pressure in the exhaust manifold 6, a portion indicated with MSCAV of fresh air inside the intake manifold 4 during the overlap phase is directly directed towards the exhaust manifold 6 through the respective exhaust valve 7, also dragging towards the exhaust manifold 6 a residual flow rate MEXH_SCAV of exhaust gases present inside the combustion chamber. This phenomenon, on the other hand, is indicated as “scavenging” and formula [12] is adjusted by replacing the downstream pressure P0 with the pressure P of the incoming air (flowing into the intake manifold 4), by replacing the upstream pressure P with the exhaust pressure PEXH and by replacing the downstream pressure T0 with the temperature TAIR of the incoming air (flowing into the intake manifold 4). Therefore, in this case, the mass MOVL is expressed as follows:



 
  
   
    
     
      M
      OVL
     
     =
     
      
       S
       id
      
      *
      
       β
       ⁡
       (
       
        
         
          P
          EXH
         
         P
        
        ,
        n
       
       )
      
      *
      
       P
       
        P
        
         0
         ⁢
         
          _
          ⁢
          REF
         
        
       
      
      *
      
       
        
         T
         
          0
          ⁢
          
           _
           ⁢
           REF
          
         
        
        
         T
         AIR
        
       
      
      *
      
       1
       n
      
     
    
   
   
    
     [
     15
     ]
    
   
  
 



The residual flow rate MEXH_SCAV of exhaust gases present inside the combustion chamber and dragged towards the exhaust manifold 6 can be expressed as follows:

MEXH_SCAV=fSCAV(MOVL,n)*PEXH*VCC(R*TEXH)  [16]

Quantities MOVL, n, PEXH, VCC, R and TEXH have the meaning already discussed above. The function fSCAV is experimentally determined through a (2d) map variable as a function of the speed n of the internal combustion engine 1 and of the mass MOVL.
The portion MSCAV of fresh air inside the intake manifold 4 directly directed towards the exhaust manifold 6 through the respective exhaust valve 7 during the overlap phase can hence be expressed as follows:

MSCAV=MOVL−MEXH_SCAV  [17]

In other words, the portion MSCAV of fresh air inside the intake manifold 4 directly directed towards the exhaust manifold 6 is equal to the mass MOVL minus the residual flow rate MEXH_SCAV of exhaust gases present inside the combustion chamber and dragged towards the exhaust manifold 6.
The model is finally suited to determine the variable OFF, which takes into account the gases produced by the combustion in the previous work cycle and present inside the cylinder 3 (because they did not flow out of the cylinder 3 or because they were re-sucked into the cylinder 3). The calculation of the variable OFF changes as a function of the work conditions, in particular as a function of the ratio between the pressure in the intake manifold 4 and the pressure in the exhaust manifold 6.
In case the pressure in the exhaust manifold 6 is greater than the pressure in the intake manifold 4 (“inner EGR” operating mode), the variable OFF corresponds to the total mass MEGRI of “inner EGR” expressed through formula [14].
On the other hand, in case the pressure in the intake manifold 4 is greater than the pressure in the exhaust manifold 6 (“washing” operating mode), the variable OFF is expressed through the following formula [16]:

OFF=PEXH*VCC/(R*TEXH)−MEXH_SCAV  [18]

In case the pressure in the intake manifold 4 is greater than the pressure in the exhaust manifold 6, indeed, the gases produced by the combustion in the previous work cycle and present inside the cylinder 3 (because they did not flow out of the cylinder 3) are at least partially directly directed towards the exhaust manifold 6 during the overlap phase through the respective exhaust valve 7. The value assumed by the variable OFF is substantially positive or equal to zero in case the entire flow rate of the gases produced by the combustion in the previous work cycle and present inside the cylinder 3 is directly directed towards the exhaust manifold 6 during the overlap phase; the electronic control unit 30 is configured to saturate the variable OFF to the zero value.
According to a further variant, in case, due to dynamic and cooling effects of the combustion chamber of the cylinder 3, the variable OFF assumes a negative value, the electronic control unit 30 is configured to saturate the variable OFF to a negative value.
According to a further variant, the ideal gas law [5] can be further generalized in the way expressed by formulas [19] and [20] below in order to estimate the mass m of air trapped in the cylinder 3:

m=(P*V−OFF)*Kt*K1(VVTI,VVTE)*K2(VVTE,n)  [19]

m=(P*V(IVC,n)*K(P,n)−OFF)*Kt*K1*K2  [20]

where:

Kt represents the product of the previously discussed functions f1(T, P) and f2(TH2O, P);
OFF represents the variable (mass) taking into account the gases produced by the combustion in the previous work cycle and present inside the cylinder 3 (because they did not flow out of the cylinder 3 or because they were re-sucked into the cylinder 3);
K1(VVTI, VVTE) is a multiplying coefficient taking into account the angular extent VVTI of the difference relative to the reference values of the intake valve 5 and the angular extent VVTE of the difference relative to the reference values of the exhaust valve 7; and
K2(VVTE, n) is a multiplying coefficient taking into account the angular extent VVTE of the difference relative to the reference values of the exhaust valve 7 and the speed n of rotation of the internal combustion engine 1 (or number of revolutions rpm).
The law [19] used to obtain the mass m of air trapped in the cylinder 3 is used as model to calculate the quantity of fuel to be injected into the cylinder 3 in order to obtain an objective value of the air/fuel ratio of the exhaust gases. In other words, once the mass m of air trapped in each cylinder 3 for each cycle has been determined through the model, the electronic control unit 30 determines the quantity of fuel to be injected into the cylinder 3 allowing the objective value of the air/fuel ratio of the exhaust gases to be reached.
According to one embodiment, in the electronic control unit 30 there is also stored a calculation chain which, from the request for torque made by the user by acting upon the accelerator pedal, is capable of providing the mass mobj of combustion air needed by each cylinder 3 to fulfil the torque request. The calculation chain requires the user to act upon the accelerator pedal, thus determining, through maps stored in the electronic control unit 30 and knowing the speed n of rotation of the internal combustion engine 1 (or number of revolutions), the torque Cr requested to the drive shaft 10; the torque Cr requested to the drive shaft 10 is then preferably added to the pumping torques and to the torques of the auxiliary elements so as to obtain the total torque Ct requested to the drive shaft 10; then the torque Ct* requested for each cylinder 3 is calculated. Once the torque Ct* requested for each cylinder 3 has been determined, the calculation chain determines the mass mobj of combustion air needed by each cylinder 3 to obtain said torque value Ct*.
Once the mass mobj of combustion air needed by each cylinder 3 to obtain said torque value Ct* has been obtained, the electronic control unit 30 is designed to use law [19] or [20] of the model in a reverse manner relative to what discussed above. In other words, for a given value of the mass mobj of combustion air needed by each cylinder 3 (which, in this case, corresponds to the mass m of air trapped in each cylinder 3 for each cycle in formula [19] or [20]), law [19] or [20] is used to calculate the objective pressure value POBJ inside the intake manifold 4. In particular, by replacing the mass m of air trapped in each cylinder 3 for each cycle with the mass mobj of combustion air needed by each cylinder 3 and by replacing the mean P of the pressure for the engine cycle inside the intake manifold 4 with the objective pressure value POBJ inside the intake manifold 4 in formula [20], the following law [21] is obtained:

POBJ=[mobj/(Kt*K1*K2)+OFF]/(V(IVC,n)*K(P,n))  [21]

The throttle valve 12 is controlled by the electronic control unit 30 so as to obtain, inside the intake manifold 4, the objective pressure value POBJ determined through law [21].
The model stored inside the electronic control unit 30 uses measured and/or estimated physical quantities (such as, for example, the temperature and pressure values) and measured and/or objective physical quantities (such as, for example, the VVT timing of the intake valves 5 and of their overlap with the exhaust valves 7).
In case the internal combustion engine 1 comprises the low-pressure exhaust gas recirculation circuit EGRLP, the total mass MEGR_TOT recirculated through the low-pressure circuit EGRLP is calculated through formula [8], which was discussed in the description above.
On the other hand, the mass MEGR recirculated through the low-pressure circuit EGRLP for each cylinder 3 is calculated through the following formula:

MEGR=MEGR_TOT/(n*120*NCYL)  [22]

where:

n represents the speed of rotation of the internal combustion engine 1 (or number of revolutions rpm);
NCYL represents the number of cylinders 3; and
MEGR_TO represents the total mass recirculated through the low-pressure circuit EGRLP calculated by the electronic control unit 30 with a model or, alternatively, measured using a dedicated sensor.
MEGR represents the mass recirculated through the low-pressure circuit EGRLP for each cylinder 3.
Hence, laws [19] and [20] can be further generalized as follows in order to also take into account the mass MEGR recirculated through the low-pressure circuit EGRLP:

m=(P*V−OFF)*Kt*K1*K2−MEGR  [23]

m=(P*V(IVC,n)*K(P,n)−OFF)*Kt*K1*K2−MEGR  [24]

where:

Kt represents the product of the previously discussed functions f1(T, P) and f2(TH2O, P);
OFF represents the variable (mass) taking into account the gases produced by the combustion in the previous work cycle and present inside the cylinder 3 (because they did not flow out of the cylinder 3 or because they were re-sucked into the cylinder 3);
MEGR represents the mass recirculated through the EGR circuit for each cylinder 3; and
K1 K2 are the empirical multiplying coefficients taking into account the angular extent VVTI of the difference relative to the reference values of the intake valve 5, the angular extent VVTE of the difference relative to the reference values of the exhaust valve 7 and the speed n of rotation of the internal combustion engine 1 (or number of revolutions rpm).
The description above, which deals with the calculation of the mass MEGR recirculated through the low-pressure circuit EGRLP for each cylinder 3, can also be applied, in an equivalent manner, in case of a high-pressure exhaust gas recirculation circuit EGRHP.
Finally, the total mass MTOT of air taken in by the internal combustion engine 1 is calculated through the following formula:

MTOT=(m+MSCAV+MEXH_SCAV)*NCYL  [23]

where:

MTOT represents the total mass of air taken in by the internal combustion engine 1;
m represents the mass of air trapped in each cylinder 3;
MSCAV represents the portion of fresh air inside the intake manifold 4 directly directed towards the exhaust manifold 6 for each cylinder 3 through the respective exhaust valve 7 during the overlap phase and obtained using formula [17];
MEXH_SCAV represents the mass of exhaust gases present in the cylinder 3 from the previous cycle and expelled, upon exhaust, by the scavenging flow; and
NCYL represents the number of cylinders 3.
On the other hand, the mass of gases OFF produced by the combustion in the previous work cycle and present inside the cylinder 3, in case the pressure of the intake manifold 4 is greater than the pressure in the exhaust manifold 6, is calculated through the following equation:

OFF=PEXH*VCC/(R*TEXH)−MEXH_SCAV 

where:

PEXH represents the pressure of the gas flow in the exhaust;
TEXH represents the temperature of the gas flow in the exhaust;
VCC represents the dead volume of the combustion chamber of the cylinder 3;
MEXH_SCAV represents the residual mass of exhaust gases present inside the combustion chamber of the cylinder 3 and directly directed towards the exhaust manifold 6 through the respective exhaust valve 7; and
R represents the constant of the mixture of fresh air and/or exhaust gases.
If the internal combustion engine 1 comprises a low-pressure gas recirculation circuit, the method comprises the further steps of calculating a quantity REGR indicating the incidence of a low-pressure circuit on the gas mixture flowing in the intake duct 6:

REGR=MEGR_LP/MTOT 

where:

MTOT represents the mass of the gas mixture flowing through the intake duct 6;
MEGR_LP represents the mass of exhaust gases recirculated through the low-pressure circuit which flows in the intake duct 6; and calculating the mass of gases OFF produced by the combustion in the previous work cycle and present inside the cylinder 3 using the following equation:

OFF=PEXH*VCC/(R*TEXH)−MEXH_SCAV*(1−REGR)

The mass of gases OFF produced by the combustion in the previous work cycle and present inside the cylinder 3 is caused to be equal to zero (is saturated), in case the entire flow rate of gases produced by the combustion in the previous work cycle and present inside the cylinder 3 is directly directed towards the exhaust manifold 6 during the overlap phase through the respective exhaust valve 7.
On the other hand, the residual mass MEXH_SCAV of exhaust gases is calculated as a function of the mass MOVL flowing from the intake to the exhaust through the intake valve 5 and the exhaust valve 7. The residual mass MEXH_SCAV of exhaust gases is calculated as a function of the speed n of rotation of the internal combustion engine 1. The residual mass MEXH_SCAV of exhaust gases is advantageously calculated as a function of the pressure PEXH and of the temperature TEXH of the gas flow in the exhaust and of the dead volume VCC of the combustion chamber of the cylinder 3.
The residual mass MEXH_SCAV of exhaust gases is, in particular, calculated using the following equation:

MEXH_SCAV=f(MOVL,n)*PEXH*VCC/(R*TEXH)  [14]

where:

PEXH, TEXH represent the pressure and temperature of the gas flow in the exhaust;
VCC represents the dead volume of the combustion chamber of the cylinder 3;
n represents the speed of rotation of the internal combustion engine 1; and
MOVL represents the mass flowing from the exhaust to the intake and sucked again into the cylinder 3, during the intake stroke, through the intake valve 5.
The residual mass MEXH_SCAV of exhaust gases is calculated using the following equation:

MEXH_SCAV=MOVL*f(MOVL,n)*g1(G,n)

where:

n represents the speed of rotation of the internal combustion engine 1;
MOVL represents the mass flowing from the intake to the exhaust through the intake valve 5 and the exhaust valve 7; and
G represents the centre of gravity of the overlap phase.
Function g1 is defined in an experimental phase through a (2d) map as a function of the speed n of rotation of the internal combustion engine 1 and of the centre G of gravity of the overlap phase, respectively.
The mass MOVL is determined using the following equation:



 
  
   M
   OVL
  
  =
  
   
    S
    id
   
   *
   
    β
    ⁡
    (
    
     
      P
      
       P
       0
      
     
     ,
     n
    
    )
   
   *
   
    
     P
     0
    
    
     P
     
      0
      ⁢
      
       _
       ⁢
       REF
      
     
    
   
   *
   
    
     
      T
      
       0
       ⁢
       
        _
        ⁢
        REF
       
      
     
     
      T
      0
     
    
   
   *
   
    1
    n
   
  
 



where:

Sid represents the ideal section;
n represents the speed of the internal combustion engine (1);
P0_REF represents the reference pressure upstream of the passage section (or overlap);
T0_REF represents the reference temperature upstream of the passage section (or overlap);
T0 represents the temperature upstream of the passage section (or overlap); and
P0, P represent the pressure upstream and downstream, respectively, of the passage section (or overlap).
The development of function β is shown in FIG. 5 as a function of the compressibility factor P/P0. Function β is experimentally characterized as a function of the speed n of the internal combustion engine 1.
The ideal section S is calculated as the product between a first function A of the speed n of the internal combustion engine 1 and of the duration OVL of the overlap phase, during which each intake valve 5 and the respective exhaust valve 7 are simultaneously open, and a second function G of the speed n of the internal combustion engine 1 and of the angular difference between the top dead centre PMS and the centre of gravity G of the overlap phase.
The mass (m) of air trapped in each cylinder 3 is further calculated as a function of a number of (two) multiplying coefficients K1, K2, which take into account the angular extent VVTI of a difference relative to the reference values of the intake valve 5, the angular extent VVTE of a difference relative to the reference values of the exhaust valve 7 and the speed n of rotation of the internal combustion engine 1.
In one embodiment, the mass m of air trapped in each cylinder 3 is calculated as a function of a first multiplying coefficient K1, which takes into account the angular extent VVTI of a difference relative to the reference values of the intake valve 5 and the angular extent VVTE of a difference relative to the reference values of the exhaust valve 7, and of a second multiplying coefficient K2, which takes into account the speed n of rotation of the internal combustion engine 1 and the angular extent VVTE of a difference relative to the reference values of the exhaust valve 7.
In case the internal combustion engine 1 further comprises the exhaust gas recirculation circuit EGRLP, EGRLP, the method involves determining the mass m of air trapped in each cylinder 3 also as a function of a mass MEGR recirculated through the circuit EGRLP, EGRLP for each cylinder 3.
Hence, the mass m of air trapped in each cylinder 3 is calculated using the following formula:

m=(P*V−OFF)*f1(T,P)*f2(TH2O,P)*−MEGR  [22]

where:

f1 f2 are functions taking into account the temperature T inside the intake manifold 4, the intake pressure P and the temperature TH2O of the coolant fluid of the internal combustion engine 1;
OFF represents the mass of gases produced by the combustion in the previous work cycle and present inside the cylinder 3; and
MEGR represents the mass recirculated through the EGR circuit for each cylinder 3.
The dead volume VCC of the combustion chamber of the cylinder 3 is a function of the speed n of rotation of the internal combustion engine 1 and of a first parameter TVC, which is alternatively equal to the closing delay angle EVC of the exhaust valve 7 or to the greatest value between zero and the smallest value between the closing delay angle EVC of the exhaust valve 7 and the value of the opening advance angle IVO of the intake valve 5 multiplied by −1. The volume is determined using a map, which is a function of the speed of rotation n of the internal combustion engine 1 and of the first parameter TVC, and using a map, which is a function of the speed n of rotation of the internal combustion engine 1 and of the duration OVL of the overlap phase.
The method further comprises determining, based on a calculation model using measured and/or estimated physical quantities, the mass mobj of combustion air needed by each cylinder 3 in order to fulfil the torque request Ct*; and determining the objective pressure value POBJ inside the intake manifold 4 based on said model as a function of the mass mobj of combustion air needed by each cylinder 3 in order to fulfil the torque request Ct*, of the actual inner volume V of each cylinder 3 and of the first group of reference quantities. The method further involves controlling the throttle valve 12 to obtain the objective pressure value POBJ inside the intake manifold 4.
Finally, the method comprises detecting a first angular extent VVTI of the opening or closing difference relative to the reference values concerning the intake valve 5; acquiring the reference closing angle IVCref of the intake valve 5; and determining the closing delay angle IVC of the intake valve 5 using the respective reference angle IVCref and the first angular extent VVTI. Furthermore, the method comprises detecting a second angular extent VVTE of the opening or closing difference relative to the reference values concerning the exhaust valve 7; acquiring the reference closing angle EVCref of the exhaust valve 7; and determining the closing delay angle EVC of the exhaust valve 7 using the respective reference angle EVCref and the second angular extent VVTE.
The mass MSCAV of fresh air inside the intake manifold 4 directly directed towards the exhaust manifold 6 is calculated as the difference between the mass Move, flowing through the overlap and the residual mass MEXH_SCAV of exhaust gases present inside the combustion chamber of the cylinder 3 and directly directed towards the exhaust manifold 6 through the respective exhaust valve 7.
In case the internal combustion engine 1 comprises the low-pressure exhaust gas recirculation circuit, the method comprises calculating the quantity REGR and calculating the mass MSCAV of fresh air inside the intake manifold 4 directly directed towards the exhaust manifold 6 using the following formula:

MSCAV=(MOVL−MEXH_SCAV)*(1−REGR)

It is then possible to use the masses for each cylinder 3 and for each engine cycle in order to calculate the flow rates of the internal combustion engine 1, taking into account the number of cylinders 3 and the engine speed n (in particular, multiplying the number of cylinders 3 by the engine speed n multiplied by 1⁄2).
The description above explicitly relates to a supercharged internal combustion engine 1, but the strategy described herein can also find advantageous application in an internal combustion engine 1 which is not provided with a supercharging system.
The advantages of the model described herein are evident from the description above.
In particular, the model described herein represents a method that allows manufacturers to determine the mass m of air trapped in each cylinder 3, the total mass MTOT of air taken in by the internal combustion engine 1, the scavenging mass MSCAV and the inner EGR mass MEGRI in a manner that is deemed to be efficient (i.e. with an adequate precision), effective (i.e. quickly and without requiring an excessive calculation power for the electronic control unit 30) and economic (i.e. without requiring the installation of expensive additional components and/or sensors, such as for example the air flow meter).
The invention has been described in an illustrative manner. It is to be understood that the terminology which has been used is intended to be in the nature of words of description rather than of limitation. Many modifications and variations of the invention are possible in light of the above teachings. Therefore, within the scope of the appended claims, the invention may be practiced other than as specifically described.

## Claims

1. A method to determine a mass (m) of air trapped in at least one cylinder of an internal combustion engine; wherein the internal combustion engine includes an intake manifold and an exhaust manifold, with a piston that is moveably supported in the at least one cylinder such that the piston and the at least one cylinder define a combustion chamber and wherein the at least one cylinder includes an intake valve that opens and closes and an exhaust valve that opens and closes such that the at least one cylinder is in fluid communication with the intake manifold from which the at least one cylinder receives an air/fuel mixture through the intake valve when the intake valve is open, and in fluid communication with the exhaust manifold into which the at least one cylinder introduces exhaust gases produced by combustion through the exhaust valve when the exhaust valve is open; wherein each intake of the air/fuel mixture into the at least one cylinder, each combustion of the air/fuel mixture in the combustion chamber and the exhaust gases from the at least one cylinder together define a work cycle and wherein the opening and closing of the intake valve and the exhaust valve are controlled so as to change the timing of their respective opening and closing; wherein the internal combustion engine is controlled by an electronic control unit which is electrically connected to an intake manifold temperature sensor which senses a temperature of the air/fuel mixture in the intake manifold and is electrically connected to an intake manifold pressure sensor that senses a pressure of the air/fuel mixture present in the intake manifold and is electrically connected to a rotational speed sensor which senses an angular position and, hence, a speed (n) of rotation of the internal combustion engine; the method comprising the steps of:
detecting the pressure of the air/fuel mixture inside the intake manifold using the intake manifold pressure sensor, detecting the speed (n) of rotation of the internal combustion engine using the rotational speed sensor and detecting a closing delay angle (IVC) of the intake valve associated with the cylinder;
determining a mass of gases produced by the combustion in a previous work cycle and present inside the at least one cylinder;
determining an inner volume (V) of the at least one cylinder as a function of the speed (n) of rotation of the internal combustion engine and of the closing delay angle (IVC) of the intake valve associated with the at least one cylinder; and
determining the mass of air (m) trapped in the at least one cylinder through a product of a pressure (P) inside the intake manifold by the inner volume (V) of the at least one cylinder, from which the mass of gases is subtracted using a calculation model stored in the electronic control unit;
wherein the step of determining the mass of gases comprises the sub-steps of:
calculating a mass (MOVL) of gases which flows, during an overlap phase in which the intake valve and the exhaust valve are simultaneously open, from an exhaust to an intake and which is re-sucked into the at least one cylinder, during the following intake phase, through the intake valve associated with the at least one cylinder using the following formula, in which the combustion chamber of the at least one cylinder represents a passage section:




 
  
   M
   OVL
  
  =
  
   
    S
    id
   
   *
   
    β
    ⁡
    (
    
     
      P
      
       P
       0
      
     
     ,
     n
    
    )
   
   *
   
    
     P
     0
    
    
     P
     
      0
      ⁢
      
       _
       ⁢
       REF
      
     
    
   
   *
   
    
     
      T
      
       0
       ⁢
       
        _
        ⁢
        REF
       
      
     
     
      T
      0
     
    
   
   *
   
    1
    n
   
  
 



where:
Sid represents an area of the passage section;
n represents the speed of rotation (n) of the internal combustion engine;
P0_REF represents a reference pressure in the exhaust manifold;
T0_REF represents a reference temperature in the exhaust manifold;
T0 represents a temperature in the exhaust manifold;
P0 represents a pressure in the exhaust manifold;
P represents the pressure in the intake manifold;
wherein the area (Sid) of the passage section is calculated as a product between a first function (A) of the speed (n) of rotation of the internal combustion engine and of a duration (OVL) of the overlap phase and a second function of the speed (n) of rotation of the internal combustion engine and of an angular difference between a top dead center (PMS) and a center of gravity of the overlap phase; and
calculating the mass of gases as a function of the mass (MOVL) of gases to determine a quantity of fuel to be injected into said cylinder.

2. The method according to claim 1, wherein the electronic control unit is electrically connected to an exhaust gas pressure sensor that senses a pressure (PEXH) of a gas flowing through the exhaust manifold and electrically connected to an exhaust gas temperature sensor that senses the temperature (TEXH) of the gas flowing through the exhaust manifold, said method comprising the further steps of:
detecting the pressure (PEXH) of the gas flowing through the exhaust manifold and detecting the temperature (TEXH) of the gas flowing through the exhaust manifold;
calculating the mass OFF of gases using the following equation:


OFF=PEXH*VCC/R*TEXH+MOVL 


VCC represents a dead volume of the combustion chamber of the cylinder;
R represents a constant of a mixture of fresh air and exhaust gases; and
TEXH represents the temperature the gas flow in the exhaust manifold.

3. The method according to claim 1 and comprising the further step of calculating the mass of gases, in case the pressure (P) inside the intake manifold is greater than a pressure (PEXH) of a gas flow in the exhaust, using the following equation:


OFF=PEXH*VCC/(R*TEXH)−MEXH_SCAV 


where:
VCC represents a dead volume of the combustion chamber of the at least one cylinder;
MEXH_SCAV represents a residual mass of exhaust gases present inside the combustion chamber and directed towards the exhaust manifold through the respective exhaust valve;
R represents a constant of the mixture of fresh air and exhaust gases; and
TEXH represents a temperature of the gas flow in the exhaust manifold.

4. The method according to claim 3, wherein the internal combustion engine includes a low-pressure exhaust gas recirculation circuit; the method includes the further steps of calculating a quantity (REGR) indicating the incidence of the low-pressure exhaust gas recirculation circuit on the fresh air/gas mixture flowing in an intake duct; and calculating the mass OFF of gases also as a function of said quantity (REGR) indicating the incidence of the low-pressure exhaust gas recirculation circuit.

5. The method according to claim 1, wherein the inner volume (V) of the cylinder is further calculated using a first map, which is a function of the closing delay angle (IVC) of the intake valve and of the speed (n) of rotation of the internal combustion engine, as well as using a second map, which is a function of the pressure (P) inside the intake manifold and of the speed (n) of rotation of the internal combustion engine.

6. The method according to claim 1, wherein the mass (m) of air trapped in the at least one cylinder is multiplied by a first factor, which is a function of a temperature (T) inside the intake manifold and of the pressure (P) inside the intake manifold, and by a second factor, which is a function of a temperature (TH2O) of a coolant fluid of the internal combustion engine and of the pressure (P) inside the intake manifold.

7. The method according to claim 1, wherein the mass (m) of air trapped in the at least one cylinder is calculated as a function of a first multiplying coefficient (K1) and a second multiplying coefficient (K2), which take into account an angular extent (VVT1) of a difference relative to reference values of the intake valve, an angular extent (VVTE) of a difference relative to the reference values of the exhaust valve and the speed (n) of rotation of the internal combustion engine.

8. The method according to claim 7, wherein the mass (m) of air trapped in the at least one cylinder is calculated as a function of the first multiplying coefficient (K1), which takes into account the angular extent (VVT1) of the difference relative to the reference values of the intake valve and the angular extent (VVTE) of the difference relative to the reference values of the exhaust valve, and of the second multiplying coefficient (K2), which takes into account the speed (n) of rotation of the internal combustion engine and the angular extent (VVTE) of the difference relative to the reference values of the exhaust valve.

9. The method according to claim 1, wherein the internal combustion engine further includes an EGR exhaust gas recirculation circuit (EGRLP, EGRHP), which includes, in turn, a bypass duct; wherein along the bypass duct there is arranged an EGR valve, which adjusts a flow rate of the exhaust gases flowing through the bypass duct; the method further includes the step of determining the mass (m) of air trapped in the at least one cylinder as a function of a mass (MEGR) recirculated through the EGR exhaust gas recirculation circuit (EGRLP, EGRHP) for the at least one cylinder.

10. The method according to claim 1, wherein a dead volume (VCC) of the combustion chamber of the at least one cylinder is a function of the speed (n) of rotation of the internal combustion engine and of a first parameter (TVC), which is equal to a closing delay angle (EVC) of the exhaust valve.

11. The method according to claim 1, wherein a volume (VCC) of the combustion chamber is determined using a third map, which is a function of the speed (n) of rotation of the internal combustion engine and of a first parameter (TVC), which is equal to a closing delay angle (EVC) of the exhaust valve.

12. The method according to claim 1 and comprising the further steps of:
determining, based on a calculation model using measured and/or estimated physical quantities, a mass (mobj) of combustion air needed by the at least one cylinder in order to fulfill a torque request (Ct*); and
determining an objective pressure value (POBJ) inside the intake manifold based on said calculation model as a function of the mass (mobj) of combustion air needed by the at least one cylinder in order to fulfill the torque request (Ct*), of the inner volume (V) of the at least one cylinder and of a first group of reference quantities.

13. The method according to claim 12, wherein the internal combustion engine includes a valve, which is designed to adjust a flow rate of a gas mixture comprising both exhaust gases and fresh air through an intake duct, directed towards the intake manifold; the method further including the step of controlling said valve so as to obtain the objective pressure value (POBJ) inside the intake manifold.

14. The method according to claim 1 and comprising the further step of calculating a mass (MSCAV) of fresh air inside the intake manifold directed towards the exhaust manifold by using the difference between the mass (MOVL) of gases which flows, during the overlap phase in which the intake valve and the exhaust valve are simultaneously open, and a residual mass (MEXH_SCAV) of exhaust gases present inside the combustion chamber and directed towards the exhaust manifold through the at least one respective exhaust valve.

15. The method according to claim 1, wherein a dead volume (VCC) of the combustion chamber is determined as a function of the speed (n) of rotation of the internal combustion engine and of a first parameter (TVC), which is equal to a greatest value between zero and the smallest value between a closing delay angle (EVC) of the exhaust valve and a value of an opening advance angle (IVO) of the intake valve multiplied by −1.

16. The method according to claim 1, wherein a volume (VCC) of the combustion chamber is determined using a third map, which is a function of the speed (n) of rotation of the internal combustion engine and of a first parameter (TVC), which is equal to a greatest value between zero and the smallest value between a closing delay angle (EVC) of the exhaust valve and a value of an opening advance angle (IVO) of the intake valve multiplied by −1, and using a fourth map, which is a function of the speed (n) of rotation of the internal combustion engine and of the duration of the overlap phase (OVL).

17. The method according to claim 1 and comprising the further steps of:
determining a pressure (PEXH) of a gas flowing through the exhaust manifold using a first model and determining a temperature (TEXH) of the gas flowing through the exhaust manifold using a second model;
calculating the mass OFF of gases using the following equation:


OFF=PEXH*VCC/R*TEXH+MOVL 


VCC represents a dead volume of the combustion chamber of the cylinder;
R represents a constant of a mixture of fresh air and exhaust gases; and
TEXH represents the temperature the gas flow in the exhaust manifold.

18. A method to determine a mass (m) of air trapped in at least one cylinder of an internal combustion engine; wherein the internal combustion engine includes an intake manifold and an exhaust manifold, a piston that is moveably supported in the at least one cylinder such that the piston and the at least one cylinder define a combustion chamber and wherein the at least one cylinder includes an intake valve that opens and closes and an exhaust valve that opens and closes such that the at least one cylinder is in fluid communication with the intake manifold from which the at least one cylinder receives an air/fuel mixture through the intake valve when the intake valve is open, and in fluid communication with the exhaust manifold into which the at last one cylinder introduces exhaust gases produced by combustion through the exhaust valve when the exhaust valve is open; wherein each intake of the air/fuel mixture into the at least one cylinder, each combustion of the air/fuel mixture in the combustion chamber and the exhaust gases from the at least one cylinder together define a work cycle and wherein the opening and closing of the intake valve is controlled so as to change the timing of the opening and closing; wherein the internal combustion engine is controlled by an electronic control unit which is electrically connected to an intake manifold temperature sensor which senses a temperature of the air/fuel mixture in the intake manifold and is electrically connected to an intake manifold pressure sensor that senses a pressure of the air/fuel mixture present in the intake manifold and electrically connected to a rotational speed sensor which senses the angular position and, hence, the speed (n) of rotation of the internal combustion engine; the method comprises the steps of:
detecting the pressure of the air/fuel mixture inside the intake manifold using the intake manifold pressure sensor, detecting the speed (n) of rotation of the internal combustion engine using the rotational speed sensor and detecting a closing delay angle (IVC) of the intake valve associated with the at least one cylinder;
determining a mass of gases produced by the combustion in a previous work cycle and present inside the at least one cylinder;
determining an inner volume (V) of the cylinder as a function of the speed (n) of rotation of the internal combustion engine and of the closing delay angle (IVC) of the intake valve associated with the at least one cylinder; and
determining the mass of air (m) trapped in the at least one cylinder through a product of a pressure (P) inside the intake manifold by the inner volume (V) of the at least one cylinder, from which the mass of gases is subtracted using a calculation model stored in the electronic control unit;
wherein the step of determining the mass of gases comprises the sub-steps of:
calculating a mass (MOVL) of gases which flows, during an overlap phase in which the intake valve and the exhaust valve are simultaneously open, from an exhaust to an intake and which is re-sucked into the at least one cylinder, during the following intake phase, through the respective intake valve associated with the at least one cylinder using the following formula, in which the combustion chamber represents a passage section:




 
  
   M
   OVL
  
  =
  
   
    S
    id
   
   *
   
    β
    ⁡
    (
    
     
      P
      
       P
       0
      
     
     ,
     n
    
    )
   
   *
   
    
     P
     0
    
    
     P
     
      0
      ⁢
      
       _
       ⁢
       REF
      
     
    
   
   *
   
    
     
      T
      
       0
       ⁢
       
        _
        ⁢
        REF
       
      
     
     
      T
      0
     
    
   
   *
   
    1
    n
   
  
 



where:
Sid represents an area of the passage section;
n represents the speed of rotation (n) of the internal combustion engine;
P0_REF represents a reference pressure in the exhaust manifold;
T0_REF represents a reference temperature in the exhaust manifold;
T0 represents a temperature in the exhaust manifold;
P0 represents a pressure in the exhaust manifold;
P represents the pressure in the intake manifold;
wherein the area (Sid) of the passage section is calculated as a product between a first function (A) of the speed (n) of rotation of the internal combustion engine and of a duration (OVL) of the overlap phase and a second function of the speed (n) of rotation of the internal combustion engine and of an angular difference between a top dead center (PMS) and a center of gravity of the overlap phase; and
calculating the mass of gases as a function of the mass (MOVL) of gases to determine a quantity of fuel to be injected into said cylinder.

19. A method to determine a mass (m) of air trapped in at least one cylinder of an internal combustion engine; wherein the internal combustion engine includes an intake manifold and an exhaust manifold, a piston that is moveably supported in the at least one cylinder such that the piston and the at least one cylinder define a combustion chamber and wherein the at least one cylinder includes an intake valve that opens and closes and an exhaust valve that opens and closes such that the at least one cylinder is in fluid communication with the intake manifold from which the at least one cylinder receives an air/fuel mixture through the intake valve when the intake valve is open, and in fluid communication with the exhaust manifold into which the at least one cylinder introduces exhaust gases produced by combustion through the exhaust valve when the exhaust valve is open; wherein each intake of the air/fuel mixture into the at least one cylinder, each combustion of the air/fuel mixture in the combustion chamber and the exhaust gases from the at least one cylinder together define a work cycle and wherein the opening and closing of the exhaust valve is controlled so as to change the timing of opening and closing the exhaust valve; wherein the internal combustion engine is controlled by an electronic control unit which is electrically connected to an intake manifold temperature sensor which senses a temperature of the air/fuel mixture in the intake manifold and is electrically connected to an intake manifold pressure sensor that senses a pressure of the air/fuel mixture present in the intake manifold and electrically connected to a rotational speed sensor which senses the angular position and, hence, the speed (n) of rotation of the internal combustion engine; the method comprises the steps of:
detecting the pressure of the air/fuel mixture inside the intake manifold using the intake manifold pressure sensor, detecting the speed (n) of rotation of the internal combustion engine using the rotational speed sensor and detecting a closing delay angle (IVC) of the intake valve associated with the at least one cylinder;
determining a mass of gases produced by the combustion in a previous work cycle and present inside the at least one cylinder;
determining an inner volume (V) of the at least one cylinder as a function of the speed (n) of rotation of the internal combustion engine and of the closing delay angle (IVC) of the intake valve associated with the at least one cylinder; and
determining the mass of air (m) trapped in the at least one cylinder through a product of a pressure (P) inside the intake manifold by the inner volume (V) of the at least one cylinder, from which the mass of gases is subtracted using a calculation model stored in the electronic control unit;
wherein the step of determining the mass of gases comprises the sub-steps of:
calculating a mass (MOVL) of gases which flows, during an overlap phase in which the intake valve and the exhaust valve are simultaneously open, from an exhaust to an intake and which is re-sucked into the at least one cylinder, during the following intake phase, through the intake valve associated with the at least one cylinder using the following formula, in which the combustion chamber of the cylinder represents a passage section:




 
  
   M
   OVL
  
  =
  
   
    S
    id
   
   *
   
    β
    ⁡
    (
    
     
      P
      
       P
       0
      
     
     ,
     n
    
    )
   
   *
   
    
     P
     0
    
    
     P
     
      0
      ⁢
      
       _
       ⁢
       REF
      
     
    
   
   *
   
    
     
      T
      
       0
       ⁢
       
        _
        ⁢
        REF
       
      
     
     
      T
      0
     
    
   
   *
   
    1
    n
   
  
 



where:
Sid represents an area of the passage section;
n represents the speed of rotation (n) of the internal combustion engine;
P0_REF represents a reference pressure in the exhaust manifold;
T0_REF represents a reference temperature in the exhaust manifold;
T0 represents a temperature in the exhaust manifold;
P0 represents a pressure in the exhaust manifold;
P represents the pressure in the intake manifold;
wherein the area (Sid) of the passage section is calculated as a product between a first function (A) of the speed (n) of rotation of the internal combustion engine and of a duration (OVL) of the overlap phase and a second function of the speed (n) of rotation of the internal combustion engine and of an angular difference between a top dead center (PMS) and a center of gravity of the overlap phase; and
calculating the mass of gases as a function of the mass (MOVL) of gases to determine a quantity of fuel to be injected into the at least one cylinder.

20. A method to determine a mass (m) of air trapped in at least one cylinder of an internal combustion engine; wherein the internal combustion engine includes an intake manifold and an exhaust manifold, a piston that is moveably supported in the at least one cylinder such that the piston and the at least one cylinder define a combustion chamber and wherein the at least one cylinder includes an intake valve that opens and closes and an exhaust valve that opens and closes such that the at least one cylinder is in fluid communication with the intake manifold from which the at least one cylinder receives an air/fuel mixture through the intake valve when the intake valve is open, and in fluid communication with the exhaust manifold into which the at least one cylinder introduces exhaust gases produced by combustion through the exhaust valve when the exhaust valve is open; wherein each intake of the air/fuel mixture into the at least one cylinder, each combustion of the air/fuel mixture in the combustion chamber and the exhaust gases from the at least one cylinder together define a work cycle and wherein the opening and closing of the intake valve and/or of the exhaust valve are controlled so as to change the timing of their respective opening and closing; wherein the internal combustion engine is controlled by an electronic control unit which is electrically connected to an intake manifold temperature sensor which senses a temperature of the air/fuel mixture in the intake manifold and is electrically connected to an intake manifold pressure sensor that senses a pressure of the air/fuel mixture present in the intake manifold and electrically connected to a rotational speed sensor which senses the angular position and, hence, the speed (n) of rotation of the internal combustion engine; the method comprises the steps of:
detecting a pressure of the air/fuel mixture inside the intake manifold using the intake manifold pressure sensor, detecting the speed (n) of rotation of the internal combustion engine using the rotational speed sensor and detecting a closing delay angle (IVC) of the intake valve associated with the at least one cylinder;
determining a mass of gases produced by the combustion in a previous work cycle and present inside the at least one cylinder;
determining an inner volume (V) of the at least one cylinder as a function of the speed (n) of rotation of the internal combustion engine and of the closing delay angle (IVC) of the intake valve associated with the at least one cylinder; and
determining the mass of air (m) trapped in the at least one cylinder through a product of a pressure (P) inside the intake manifold by the inner volume (V) of the at least one cylinder, from which the mass of gases is subtracted using a calculation model stored in the electronic control unit;
wherein the step of determining the mass of gases comprises the sub-steps of:
calculating a mass (MOVL) of gases which flows, during an overlap phase in which the intake valve and the exhaust valve are simultaneously open, from an exhaust to an intake and which is re-sucked into the at least one cylinder, during the following intake phase, through the intake valve associated with the at least one cylinder using the following formula, in which the combustion chamber of the at least one cylinder represents a passage section:




 
  
   M
   OVL
  
  =
  
   
    S
    id
   
   *
   
    β
    ⁡
    (
    
     
      P
      
       P
       EXH
      
     
     ,
     n
    
    )
   
   *
   
    
     P
     EXH
    
    
     P
     
      0
      ⁢
      
       _
       ⁢
       REF
      
     
    
   
   *
   
    
     
      T
      
       0
       ⁢
       
        _
        ⁢
        REF
       
      
     
     
      T
      EXH
     
    
   
   *
   
    1
    n
   
  
 



where:
Sid represents an area of the passage section;
n represents the speed of rotation (n) of the internal combustion engine;
P0_REF represents a reference pressure in the exhaust manifold;
T0_REF represents a reference temperature in the exhaust manifold;
TEXH represents a temperature of a gas flow in the exhaust manifold;
PEXH represents a pressure of a gas flow in the exhaust manifold;
P represents the pressure in the intake manifold;
wherein the area (Sid) of the passage section is calculated as a product between a first function (A) of the speed (n) of rotation of the internal combustion engine and of a duration (OVL) of the overlap phase and a second function of the speed (n) of rotation of the internal combustion engine and of an angular difference between a top dead center (PMS) and a center of gravity of the overlap phase; and
calculating the mass of gases as a function of the mass (MOVL) of gases to determine a quantity of fuel to be injected into the at least one cylinder.