# Patent

## Title
On-board aftertreatment device tail pipe hydrocarbon slip calculation

## Abstract

A method to determine hydrocarbon slip in a compression ignition engine equipped with an electronic control unit having memory and coupled to an exhaust gas aftertreatment system having a diesel oxidation catalyst, a diesel particulate filter, a selective catalyst reducer and a clean up diesel oxidation catalyst.

## Background
TECHNICAL FIELD OF THE INVENTION
Aftertreatment systems on diesel engines often contain a Diesel Oxidation Catalyst (DOC) which in many applications is utilized to oxidize a supplied hydrocarbon (via in cylinder or external dosing measures) to increase the temperature of the Diesel Particulate Filter (DPF) to regenerate or remove the soot that has been trapped within the DPF. This DOC does not necessarily burn 100% of the hydrocarbons (HC), and some of it slips the catalyst component substrates downstream of the DOC.
There is a need to develop an on board aftertreatment diagnostic to determine tail pipe hydrocarbon slip. In one embodiment, there is disclosed a method to calculate determine the hydrocarbon slip after each of these components (DPF and DOC), and also convert this determination to parts per million (ppm) all in a manner that is simple enough to be done on board the vehicle, generally within an electronic control unit (ECU), control logic is utilized to determine the ppm slip out of the DPF+DOC system to track the ability to comply with On Board Diagnostic (OBD) requirements for vehicles. This diagnostic can be utilized to track the efficiency of an included Ammonia Oxidation (AMOX) clean up catalyst should one be employed with a Selective Catalyst Reduction (SCR) system yielding tailpipe out HC slip in ppm.
BRIEF SUMMARY OF THE INVENTION
In one embodiment there is disclosed a method to determine hydrocarbon slip in an internal combustion engine, and more particularly, in a compression ignition engine, equipped with an electronic control unit having memory and coupled to an exhaust gas aftertreatment system having a diesel oxidation catalyst, a diesel particulate filter, a selective catalyst reducer and a clean up diesel oxidation catalyst. The method comprises determining a change in temperature of exhaust gas flow through the diesel oxidation catalyst component; determining a change in temperature of exhaust gas flow through the diesel particulate filter; determining a change in temperature through the selective catalyst reducer at the clean-up diesel oxidation catalyst outlet; and determining the efficiency of each exhaust gas aftertreatment component.
In another embodiment, there is disclosed a method to determine hydrocarbon slip in an compression ignition engine equipped with an electronic control unit having memory and coupled to an exhaust gas aftertreatment system having a diesel oxidation catalyst, a diesel particulate filter, a selective catalyst reducer and a clean up diesel oxidation catalyst comprising determining a change in temperature of exhaust gas flow through the diesel oxidation catalyst component; determining a change in temperature of exhaust gas flow through the diesel particulate filter; determining a change in temperature through the selective catalyst reducer at the clean-up diesel oxidation catalyst outlet; and determining the efficiency of each exhaust gas aftertreatment component.
In another embodiment, the change in exhaust gas flow temperature through the diesel catalyst component is determined according to:

ΔTdoc=Sensed DOC out Temperature−Virtual DOC Out Temperature.

In another embodiment, the change in temperature of the DPF is determined according to:

ΔTDPF=Sensed DPF Out Temperature−Virtual DPF Out Temperature.

In another embodiment, the change in temperature of the exhaust gas flow at the clean-up DOC outlet is determined according to:

ΔTSCR+Clean−upDOC=Sensed Clean-Up DOC Out Temperature−Virtual Clean-Up Out DOC Temperature.

In another embodiment, the energy released at each catalyst due to hydrocarbon burning is determined according to:

{dot over (Q)}=Cp·{dot over (m)}exh·ΔT=η·{dot over (m)}fuel·ΔH  (1)

Where η—catalyst efficiency

    
    
        {dot over (m)}fuel—HC Flow rate into the catalyst, kg/s
        ΔH—heat content of fuel, 43 MJ/kg for diesel fuel
        ΔT—Temperature increase due to HC burning (explained above)
        Cp—exhaust specific heat, kJ/kg·K
        {dot over (m)}exh—Exhaust gas flow rate through the catalyst, kg/s wherein the HC flow in mass faction, (1) becomes
    
    






  
    
      
        
          HC
          massfraction
        
        =
        
          
            
              
                m
                .
              
              fuel
            
            
              
                m
                .
              
              exh
            
          
          =
          
            
              
                
                  C
                  p
                
                ·
                Δ
              
              ⁢
              
                
              
              ⁢
              T
            
            
              
                η
                ·
                Δ
              
              ⁢
              
                
              
              ⁢
              H
            
          
        
      
    
    
      
        (
        2
        )
      
    
  




    
    
        Converting the mass fraction to volume fraction (ppm) as C14H30 (as in normal diesel fuel), (2) changes to
    
    






  
    
      
        
          HC
          
            ppm
            
              (
              
                C
                ⁢
                
                  
                
                ⁢
                14
              
              )
            
          
        
        =
        
          
            
              
                
                  m
                  .
                
                fuel
              
              /
              198
            
            
              
                
                  m
                  .
                
                exh
              
              /
              29
            
          
          =
          
            
              
                
                  29
                  ·
                  
                    C
                    p
                  
                  ·
                  Δ
                
                ⁢
                
                  
                
                ⁢
                T
              
              
                
                  198
                  ·
                  η
                  ·
                  Δ
                
                ⁢
                
                  
                
                ⁢
                H
              
            
            ·
            
              10
              6
            
          
        
      
    
    
      
        (
        3
        )
      
    
  




Since the HC concentration (ppm) is normally corrected to C1, the ppm in C1 will be




  
    
      
        
          HC
          ppm
        
        =
        
          
            HC
            
              ppm
              
                (
                
                  C
                  ⁢
                  
                    
                  
                  ⁢
                  1
                
                )
              
            
          
          =
          
            
              14
              ·
              
                
                  
                    29
                    ·
                    
                      C
                      p
                    
                    ·
                    Δ
                  
                  ⁢
                  
                    
                  
                  ⁢
                  T
                
                
                  
                    198
                    ·
                    η
                    ·
                    Δ
                  
                  ⁢
                  
                    
                  
                  ⁢
                  H
                
              
              ·
              
                10
                6
              
            
            =
            
              
                constant
                ·
                
                  
                    Δ
                    ⁢
                    
                      
                    
                    ⁢
                    T
                  
                  η
                
              
              =
              
                k
                ·
                
                  
                    Δ
                    ⁢
                    
                      
                    
                    ⁢
                    T
                  
                  η
                
              
            
          
        
      
    
    
      
        (
        4
        )
      
    
  




Apply Equ. (4) to DOC, DPF and SCR+clean up DOC respectively, we have:




  
    
      
        
          HC
          
            ppm
            
              ,
              doc
            
          
        
        =
        
          k
          ·
          
            
              Δ
              ⁢
              
                
              
              ⁢
              
                T
                doc
              
            
            
              η
              doc
            
          
        
      
    
    
      
        (
        5
        )
      
    
  
  
    
      
        
          HC
          
            ppm
            
              ,
              dpf
            
          
        
        =
        
          k
          ·
          
            
              Δ
              ⁢
              
                
              
              ⁢
              
                T
                dpf
              
            
            
              η
              dpf
            
          
        
      
    
    
      
        (
        6
        )
      
    
  
  
    
      
        
          HC
          
            ppm
            
              ,
              
                scr
                +
                cleanupdoc
              
            
          
        
        =
        
          k
          ·
          
            
              Δ
              ⁢
              
                
              
              ⁢
              
                T
                
                  scr
                  +
                  cleanupdoc
                
              
            
            
              η
              
                scr
                +
                cleanupdoc
              
            
          
        
      
    
    
      
        (
        7
        )
      
    
  







BRIEF DESCRIPTION OF THE DRAWINGS
FIG. 1 is a schematic representation of an aftertreatment device or system for use with an internal combustion engine;
FIG. 2 is a schematic representation of a software flow diagram representing one method according to the present disclosure.



DETAILED DESCRIPTION OF THE INVENTION
Turning now to the drawings wherein like numbers refer to like structures, and particularly to FIG. 1, there is disclosed a schematic representation of an aftertreatment device or system 10 for use with an internal combustion engine, including a compression ignition engine. The System 10 is connected at one end to an internal combustion engine exhaust manifold, schematically represented at 27, to vent exhaust from the engine through the system 10. The system 10 is comprised of a DOC inlet conduit 12, in fluid communication with exhaust manifold 27 at a first end and a DOC 14 at an opposite end. The DOC 14 is in close proximity to a DPF 16. DPF outlet conduit 18 is in fluid communication with SCR 20. SCR 20 is in a proximal position with a clean up DOC 22, from which exhaust outlet 24 permits the treated exhaust gas stream to vent to the atmosphere.
Temperature sensors are positioned along the system 10 to sense the temperature of the exhaust gas stream as it travels through the system. DOC inlet temperature sensor 26 and DPF inlet sensor 28 sense the temperature of the exhaust gas as it passes through the DOC 14. As one of ordinary skill will understand, the exhaust gas stream through the DOC is very hot. As the exhaust gas stream enters and leaves the DOC 14, DPF inlet temperature sensor acts as a DOC outlet temperature sensor. DPF outlet temperature sensor 30 is in close proximity to the conduit 18. DPF outlet temperature sensor 30 and SCR inlet temperature sensor 32 sense the temperature of the exhaust gas as it passes through the DPF and enters the SCR 20. The SCR outlet temperature sensor 34 senses the temperature of the exhaust gas stream as it exits the clean up DOC 22. Each of the sensors is electronically connected to an electronic control unit (ECU) 13 through electronic connections 15, 17, 19, 21, and 23, respectively. The ECU is usually on the engine and controls the operation of the engine. Specifically, the ECU contains the operating instruction and other software and algorithms necessary to the operation of the system and the engine. The ECU may be a single processor or a two processor controller, depending upon design considerations. If the ECU is a two processor controller, it may be comprised of a Motor Control Module (MCM) and a component control module (CPC), that communicate with each other. The MCM has the operating instructions for the engine, and the CPC has operating instructions for the components of the vehicle. The MCM has a memory with memory with tables or profiles that map or otherwise store values for the amount of hydrocarbons permitted in the exhaust stream at any point along the system 10. The memory may be PROM, EPROM, EEPROM, FLASH or any other memory. The MCM and the CPC communicate with the Engine and the other components. Without limitation, one suitable operating software program is DDEC available from Detroit Diesel Corporation.
FIG. 2 is a schematic representation of a software flow diagram of one method 38 of determining hydrocarbon slip in an exhaust system. By determining the change in temperature of the exhaust gas between the various components in the system 10 (FIG. 1), the hydrocarbon slip can be determined. Specifically, Step 40 is determining the change of tempareature of the exhaust gas stream through a DOC. This may be accomplished by determining the temperature difference between the DOC inlet and the DOC outlet, to determine the ΔTdoc. This permits the determination of the HCslip.doc. Step 42 is determining the change of temperature of the exhaust gas stream through the diesel particulate filter. This may be accomplished by measuring the temperature difference from the DOC outlet and the DPF outlet. The temperature difference between the DOC out and the DPF out determines the ΔTdpf. The HCslip dpf is then determined and converted to ppm. Step 44 is determining the change in temperature of the exhaust gas stream through the Selective Catalyst Reducer (SCR) and a clean up DOC outlet. The temperature difference between the SCR inlet temperature and the clean up DOC outlet temperature is used to determine the ΔTscr+clean up doc. The exhaust gas stream then passes through the clean up DOC and the HCslip tail pipe is the HCslip scr+clean up doc.
An exhaust pipe fuel injection system upstream of the DOC or an in cylinder late post injection event provides the HC species, and the DOC is designed to burn most of this fuel to increase the exhaust temperature to a level wherein the loaded soot on the DPF can be burned clean, via regeneration of the DPF or otherwise. Since the HC species may not always be oxidized 100% at each catalyst, especially as the catalysts age, an on board HC slip prediction method is necessary to permit modifications in the operation of the engine to ensure exhaust gas stream composition meets governmental regulations, or to warn the operator that the exhaust system components are in need of service.
The ΔT described earlier in this application is defined as the temperature increase across a given catalyst device due solely to the HC oxidation. The temperature rise (ΔT) is equal to the measured (sensed) catalyst out temperature minus the virtual catalyst out temperature. The virtual catalyst out temperature is determined from the measured catalyst in temperature with correction of heat exchange between the exhaust gas stream and the catalyst substrate, as well as the temperature sensor delay. Since the insulation matt material between the catalyst substrate and the metal skin of the catalyst canning plus any other insulation methods prevent significant heat transfer from the substrate to the ambient, minimal heat loss is assumed. However, it is further contemplated that this heat loss could be determined and compensated for. The heat generated from the controlled soot burning in the DPF can be neglected because of its small scale, however, it is further contemplated that this potential source of error could be compensated for. It is further contemplated that any check for the ppm of the HC in the exhaust gas stream could be done only at the end of an active regeneration where soot burning is zero, but it is also contemplated to spot check the ppm of the HC over the entire system 10.
By way of example, starting with the DOC

ΔT doc=sensed DOC out Temperature−virtual DOC out temperature.

The virtual DOC out temperature calculation needs the exhaust flow rate as input, and it can be obtained based upon engine configuration, operating condition, intake air input, and some correction tables in memory of the ECU. All other inputs are either constraints or experimental calibrated data, therefore the virtual DOC out temperature can be determined on board, as well as the ΔT doc. The same logic can be applied to the DPF and the clean up DOC catalyst to obtain ΔT dpf and ΔT scr+ clean up doc.
Assuming the energy released at each catalyst is due to hydrocarbon burning, it may be understood that:

{dot over (Q)}=Cp·{dot over (m)}exh·ΔT=η·{dot over (m)}fuel·ΔH  (1)

Where η—catalyst efficiency

    
    
        {dot over (m)}fuel—HC Flow rate into the catalyst, kg/s
        ΔH—heat content of fuel, 43 MJ/kg for diesel fuel
        ΔT—Temperature increase due to HC burning (explained above)
        Cp—exhaust specific heat, kJ/kg·K
        {dot over (m)}exh—Exhaust gas flow rate through the catalyst, kg/s
    
    


Expressing the HC flow in mass faction, (1) becomes




  
    
      
        
          HC
          massfraction
        
        =
        
          
            
              
                m
                .
              
              fuel
            
            
              
                m
                .
              
              exh
            
          
          =
          
            
              
                
                  C
                  p
                
                ·
                Δ
              
              ⁢
              
                
              
              ⁢
              T
            
            
              
                η
                ·
                Δ
              
              ⁢
              
                
              
              ⁢
              H
            
          
        
      
    
    
      
        (
        2
        )
      
    
  




Converting the mass fraction to volume fraction (ppm) as C14H30 (as in normal diesel fuel), (2) changes to




  
    
      
        
          HC
          
            ppm
            
              (
              
                C
                ⁢
                
                  
                
                ⁢
                14
              
              )
            
          
        
        =
        
          
            
              
                
                  m
                  .
                
                fuel
              
              /
              198
            
            
              
                
                  m
                  .
                
                exh
              
              /
              29
            
          
          =
          
            
              
                
                  29
                  ·
                  
                    C
                    p
                  
                  ·
                  Δ
                
                ⁢
                
                  
                
                ⁢
                T
              
              
                
                  198
                  ·
                  η
                  ·
                  Δ
                
                ⁢
                
                  
                
                ⁢
                H
              
            
            ·
            
              10
              6
            
          
        
      
    
    
      
        (
        3
        )
      
    
  




Since the HC concentration (ppm) is normally corrected to C1, the ppm in C1 will be




   
  
    
      
        
          
            
              
                
                  HC
                  ppm
                
                =
                
                ⁢
                
                  HC
                  
                    ppm
                    
                      (
                      
                        C
                        ⁢
                        
                          
                        
                        ⁢
                        1
                      
                      )
                    
                  
                
              
            
          
          
            
              
                =
                
                ⁢
                
                  14
                  ·
                  
                    
                      
                        29
                        ·
                        
                          C
                          p
                        
                        ·
                        Δ
                      
                      ⁢
                      
                        
                      
                      ⁢
                      T
                    
                    
                      
                        198
                        ·
                        η
                        ·
                        Δ
                      
                      ⁢
                      
                        
                      
                      ⁢
                      H
                    
                  
                  ·
                  
                    10
                    6
                  
                
              
            
          
          
            
              
                =
                
                ⁢
                
                  constant
                  ·
                  
                    
                      Δ
                      ⁢
                      
                        
                      
                      ⁢
                      T
                    
                    η
                  
                
              
            
          
          
            
              
                =
                
                ⁢
                
                  k
                  ·
                  
                    
                      Δ
                      ⁢
                      
                        
                      
                      ⁢
                      T
                    
                    η
                  
                
              
            
          
        
      
      
        
          (
          4
          )
        
      
    
  




If one applies Equ. (4) to DOC, DPF and SCR+clean up DOC respectively, we have:




  
    
      
        
          HC
          
            ppm
            
              ,
              doc
            
          
        
        =
        
          k
          ·
          
            
              Δ
              ⁢
              
                
              
              ⁢
              
                T
                doc
              
            
            
              η
              doc
            
          
        
      
    
    
      
        (
        5
        )
      
    
  
  
    
      
        
          HC
          
            ppm
            
              ,
              dpf
            
          
        
        =
        
          k
          ·
          
            
              Δ
              ⁢
              
                
              
              ⁢
              
                T
                dpf
              
            
            
              η
              dpf
            
          
        
      
    
    
      
        (
        6
        )
      
    
  
  
    
      
        
          HC
          
            ppm
            
              ,
              
                scr
                +
                cleanupdoc
              
            
          
        
        =
        
          k
          ·
          
            
              Δ
              ⁢
              
                
              
              ⁢
              
                T
                
                  scr
                  +
                  cleanupdoc
                
              
            
            
              η
              
                scr
                +
                cleanupdoc
              
            
          
        
      
    
    
      
        (
        7
        )
      
    
  




The HC slip after each catalyst will be:




  
    
      
        
          HC
          
            slip
            ,
            doc
          
        
        =
        
          
            
              (
              
                1
                -
                
                  η
                  doc
                
              
              )
            
            ·
            
              HC
              
                ppm
                ,
                doc
              
            
          
          =
          
            
              (
              
                1
                -
                
                  η
                  doc
                
              
              )
            
            ·
            k
            ·
            
              
                Δ
                ⁢
                
                  
                
                ⁢
                
                  T
                  doc
                
              
              
                η
                doc
              
            
          
        
      
    
    
      
        (
        8
        )
      
    
  
  
    
      
        
          HC
          
            slip
            ,
            dpf
          
        
        =
        
          
            
              (
              
                1
                -
                
                  η
                  dpf
                
              
              )
            
            ·
            
              HC
              
                ppm
                ,
                dpf
              
            
          
          =
          
            
              (
              
                1
                -
                
                  η
                  dpf
                
              
              )
            
            ·
            k
            ·
            
              
                Δ
                ⁢
                
                  
                
                ⁢
                
                  T
                  dpf
                
              
              
                η
                dpf
              
            
          
        
      
    
    
      
        (
        9
        )
      
    
  
  
    
      
        
          
            
              
                HC
                
                  slip
                  ,
                  
                    scr
                    +
                    cleanupdoc
                  
                
              
              =
              
              ⁢
              
                
                  (
                  
                    1
                    -
                    
                      η
                      
                        scr
                        +
                        cleanupdoc
                      
                    
                  
                  )
                
                ·
                
                  HC
                  
                    ppm
                    ,
                    
                      scr
                      +
                      cleanupdoc
                    
                  
                
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  (
                  
                    1
                    -
                    
                      η
                      
                        scr
                        +
                        cleanupdoc
                      
                    
                  
                  )
                
                ·
                k
                ·
                
                  
                    Δ
                    ⁢
                    
                      
                    
                    ⁢
                    
                      T
                      
                        scr
                        +
                        cleanupdoc
                      
                    
                  
                  
                    η
                    
                      scr
                      +
                      cleanupdoc
                    
                  
                
              
            
          
        
      
    
    
      
        (
        10
        )
      
    
  




Since the HC slip from upstream catalyst is the HC input for the catalyst downstream, so Equ. (6) equals to Equ. (8), and Equ. (7) equals to Equ. (9). Therefore we get:




  
    
      
        
          
            η
            dpf
          
          =
          
            
              
                
                  η
                  doc
                
                ·
                Δ
              
              ⁢
              
                
              
              ⁢
              
                T
                dpf
              
            
            
              
                
                  (
                  
                    1
                    -
                    
                      η
                      doc
                    
                  
                  )
                
                ·
                Δ
              
              ⁢
              
                
              
              ⁢
              
                T
                doc
              
            
          
        
        ⁢
        
          

        
        ⁢
        and
      
    
    
      
        (
        11
        )
      
    
  
  
    
      
        
          η
          
            scr
            +
            cleanupdoc
          
        
        =
        
          
            
              
                η
                dpf
              
              ·
              Δ
            
            ⁢
            
              
            
            ⁢
            
              T
              
                scr
                +
                cleanupdoc
              
            
          
          
            
              
                (
                
                  1
                  -
                  
                    η
                    dpf
                  
                
                )
              
              ·
              Δ
            
            ⁢
            
              
            
            ⁢
            
              T
              dpf
            
          
        
      
    
    
      
        (
        12
        )
      
    
  




Substituting Equ. (11) to Equ. (12), then Equ. (12 to Equ. (10):





   
  
    
      
        
          
            
              
                
                  HC
                  
                    slip
                    ,
                    
                      scr
                      +
                      cleanupdoc
                    
                  
                
                =
                
                ⁢
                
                  HC
                  
                    slip
                    ,
                    tailpipec
                  
                
              
            
          
          
            
              
                =
                
                ⁢
                
                  k
                  ·
                  
                    (
                    
                      
                        
                          
                            
                              (
                              
                                1
                                -
                                
                                  η
                                  doc
                                
                              
                              )
                            
                            
                              η
                              doc
                            
                          
                          ·
                          Δ
                        
                        ⁢
                        
                          
                        
                        ⁢
                        
                          T
                          doc
                        
                      
                      -
                      
                        Δ
                        ⁢
                        
                          
                        
                        ⁢
                        
                          T
                          dpf
                        
                      
                      -
                      
                        Δ
                        ⁢
                        
                          
                        
                        ⁢
                        
                          T
                          
                            scr
                            +
                            cleanupdoc
                          
                        
                      
                    
                    )
                  
                
              
            
          
        
      
      
        
          (
          13
          )
        
      
    
  




When an aftertreatment device has no clean up DOC, the ΔT scr+clean up doc=0. If there is no coating on the DPF, ΔT dpf=0. Regardless of either contingency, equation 13 is operative in both situations.
The words used in the discussion of the disclosure are words of description, not words of limitation. Those skilled in the art recognize that many variations and modifications are possible without departing from the scope and spirit of the invention as set forth in the appended claims.

## Claims

1. A method to determine hydrocarbon slip in a compression ignition engine equipped with an electronic control unit having memory; said engine coupled to an exhaust gas aftertreatment system having a diesel oxidation catalyst (DOC), a diesel particulate filter (DPF), a selective catalyst reducer (SCR) and a clean up diesel oxidation catalyst comprising:
determining a change in temperature of exhaust gas flow through the diesel oxidation catalyst component;
determining a change in temperature of exhaust gas flow through the diesel particulate filter;
determining a change in temperature through the selective catalyst reducer at the clean-up diesel oxidation catalyst outlet; and
determining efficiency of each exhaust gas aftertreatment component by determining whether the component has a coating of excess hydrocarbon.

2. The method of claim 1, wherein the change in exhaust gas flow temperature through the diesel catalyst component is determined according to:


ΔTdoc=Sensed DOC out Temperature−Virtual DOC Out Temperature.

3. The method of claim 1, wherein the change in temperature of the DPF is determined according to:


ΔTDPF=Sensed DPF Out Temperature−Virtual DPF Out Temperature.

4. The method of claim 1, wherein the change in temperature of the exhaust gas flow at the clean-up DOC outlet is determined according to:


ΔTSCR−Clean-upDOC=Sensed Clean-Up DOC Out Temperature−Virtual Clean-Up Out DOC Temperature.

5. The method of claim 1, wherein the energy released at each catalyst due to hydrocarbon burning is determined according to:


{dot over (Q)}=Cp·{dot over (m)}exh·ΔT=ƒ·{dot over (m)}fuel·ΔH  (1)

Wherein
η—catalyst efficiency
{dot over (m)}fuel—Hydrocarbon (HC) Flow rate into the catalyst, kg/s
ΔH—heat content of fuel, 43 MJ/kg for diesel fuel
ΔT—Temperature increase due to HC burning (explained above)
Cp—exhaust specific heat, kJ/kg·K
{dot over (m)}exh—Exhaust gas flow rate through the catalyst, kg/s.

6. The method of claim 5, wherein the HC flow in mass faction is determined according to:




  
    
      
        
          HC
          massfraction
        
        =
        
          
            
              
                m
                .
              
              fuel
            
            
              
                m
                .
              
              exh
            
          
          =
          
            
              
                
                  
                    C
                    p
                  
                  ·
                  Δ
                
                ⁢
                
                  
                
                ⁢
                T
              
              
                
                  η
                  ·
                  Δ
                
                ⁢
                
                  
                
                ⁢
                H
              
            
            .
          
        
      
    
    
      
        (
        2
        )

7. The method of claim 6, wherein the mass fraction may be is converted to a volume fraction according to:




  
    
      
        
          HC
          
            ppm
            
              (
              
                C
                ⁢
                
                  
                
                ⁢
                14
              
              )
            
          
        
        =
        
          
            
              
                
                  m
                  .
                
                fuel
              
              /
              198
            
            
              
                
                  m
                  .
                
                exh
              
              /
              29
            
          
          =
          
            
              
                
                  29
                  ·
                  
                    C
                    p
                  
                  ·
                  Δ
                
                ⁢
                
                  
                
                ⁢
                T
              
              
                
                  198
                  ·
                  η
                  ·
                  Δ
                
                ⁢
                
                  
                
                ⁢
                H
              
            
            ·
            
              
                10
                6
              
              .
            
          
        
      
    
    
      
        (
        3
        )

8. The method of claim 7, wherein the HC concentration is determined in ppm and corrected to C1 according to:




   
  
    
      
        
          
            
              
                
                  HC
                  ppm
                
                =
                
                ⁢
                
                  HC
                  
                    ppm
                    
                      (
                      
                        C
                        ⁢
                        
                          
                        
                        ⁢
                        1
                      
                      )
                    
                  
                
              
            
          
          
            
              
                =
                
                ⁢
                
                  14
                  ·
                  
                    
                      
                        29
                        ·
                        
                          C
                          p
                        
                        ·
                        Δ
                      
                      ⁢
                      
                        
                      
                      ⁢
                      T
                    
                    
                      
                        198
                        ·
                        η
                        ·
                        Δ
                      
                      ⁢
                      
                        
                      
                      ⁢
                      H
                    
                  
                  ·
                  
                    10
                    6
                  
                
              
            
          
          
            
              
                =
                
                ⁢
                
                  constant
                  ·
                  
                    
                      Δ
                      ⁢
                      
                        
                      
                      ⁢
                      T
                    
                    η
                  
                
              
            
          
          
            
              
                =
                
                ⁢
                
                  k
                  ·
                  
                    
                      
                        Δ
                        ⁢
                        
                          
                        
                        ⁢
                        T
                      
                      η
                    
                    .
                  
                
              
            
          
        
      
      
        
          (
          4
          )

9. The method of claim 8, wherein the HC concentration expressed in ppm at the DOC, DPF and SCR+clean up DOC respectively, is determined according to:




  
    
      
        
          HC
          
            ppm
            
              ,
              doc
            
          
        
        =
        
          k
          ·
          
            
              Δ
              ⁢
              
                
              
              ⁢
              
                T
                doc
              
            
            
              η
              doc
            
          
        
      
    
    
      
        (
        5
        )
      
    
  
  
    
      
        
          HC
          
            ppm
            
              ,
              dpf
            
          
        
        =
        
          k
          ·
          
            
              Δ
              ⁢
              
                
              
              ⁢
              
                T
                dpf
              
            
            
              η
              dpf
            
          
        
      
    
    
      
        (
        6
        )
      
    
  
  
    
      
        
          HC
          
            ppm
            
              ,
              
                scr
                +
                cleanupdoc
              
            
          
        
        =
        
          k
          ·
          
            
              
                Δ
                ⁢
                
                  
                
                ⁢
                
                  T
                  
                    scr
                    +
                    cleanupdoc
                  
                
              
              
                η
                
                  scr
                  +
                  cleanupdoc
                
              
            
            .
          
        
      
    
    
      
        (
        7
        )

10. The method of claim 9, wherein the HC slip after each catalyst is determined according to:




  
    
      
        
          HC
          
            slip
            ,
            doc
          
        
        =
        
          
            
              (
              
                1
                -
                
                  η
                  doc
                
              
              )
            
            ·
            
              HC
              
                ppm
                ,
                doc
              
            
          
          =
          
            
              (
              
                1
                -
                
                  η
                  doc
                
              
              )
            
            ·
            k
            ·
            
              
                Δ
                ⁢
                
                  
                
                ⁢
                
                  T
                  doc
                
              
              
                η
                doc
              
            
          
        
      
    
    
      
        (
        8
        )
      
    
  
  
    
      
        
          HC
          
            slip
            ,
            dpf
          
        
        =
        
          
            
              (
              
                1
                -
                
                  η
                  dpf
                
              
              )
            
            ·
            
              HC
              
                ppm
                ,
                dpf
              
            
          
          =
          
            
              (
              
                1
                -
                
                  η
                  dpf
                
              
              )
            
            ·
            k
            ·
            
              
                Δ
                ⁢
                
                  
                
                ⁢
                
                  T
                  dpf
                
              
              
                η
                dpf
              
            
          
        
      
    
    
      
        (
        9
        )
      
    
  
  
    
      
        
          
            
              
                HC
                
                  slip
                  ,
                  
                    scr
                    +
                    cleanupdoc
                  
                
              
              =
              
              ⁢
              
                
                  (
                  
                    1
                    -
                    
                      η
                      
                        scr
                        +
                        cleanupdoc
                      
                    
                  
                  )
                
                ·
                
                  HC
                  
                    ppm
                    ,
                    
                      scr
                      +
                      cleanupdoc
                    
                  
                
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  (
                  
                    1
                    -
                    
                      η
                      
                        scr
                        +
                        cleanupdoc
                      
                    
                  
                  )
                
                ·
                k
                ·
                
                  
                    
                      Δ
                      ⁢
                      
                        
                      
                      ⁢
                      
                        T
                        
                          scr
                          +
                          cleanupdoc
                        
                      
                    
                    
                      η
                      
                        scr
                        +
                        cleanupdoc
                      
                    
                  
                  .
                
              
            
          
        
      
    
    
      
        (
        10
        )

11. The method of claim 10, wherein the HC slip from an upstream catalyst is the HC input for the catalyst downstream, Equ. (6) equals to Equ. (8), and Equ. (7) equals Equ. (9), expressed according to:




  
    
      
        
          
            η
            dpf
          
          =
          
            
              
                
                  η
                  doc
                
                ·
                Δ
              
              ⁢
              
                
              
              ⁢
              
                T
                dpf
              
            
            
              
                
                  (
                  
                    1
                    -
                    
                      η
                      doc
                    
                  
                  )
                
                ·
                Δ
              
              ⁢
              
                
              
              ⁢
              
                T
                doc
              
            
          
        
        ⁢
        
          

        
        ⁢
        and
      
    
    
      
        (
        11
        )
      
    
  
  
    
      
        
          η
          
            scr
            +
            cleanupdoc
          
        
        =
        
          
            
              
                
                  η
                  dpf
                
                ·
                Δ
              
              ⁢
              
                
              
              ⁢
              
                T
                
                  scr
                  +
                  cleanupdoc
                
              
            
            
              
                
                  (
                  
                    1
                    -
                    
                      η
                      dpf
                    
                  
                  )
                
                ·
                Δ
              
              ⁢
              
                
              
              ⁢
              
                T
                dpf
              
            
          
          .
        
      
    
    
      
        (
        12
        )

12. The method of claim 11, wherein substituting Equ. (11) into Equ. (12) then substituting Equ. (12) into Equ. (10), expressed according to:




   
  
    
      
        
          
            
              
                
                  HC
                  
                    slip
                    ,
                    
                      scr
                      +
                      cleanupdoc
                    
                  
                
                =
                
                ⁢
                
                  HC
                  
                    slip
                    ,
                    tailpipec
                  
                
              
            
          
          
            
              
                =
                
                ⁢
                
                  k
                  ·
                  
                    
                      (
                      
                        
                          
                            
                              
                                (
                                
                                  1
                                  -
                                  
                                    η
                                    doc
                                  
                                
                                )
                              
                              
                                η
                                doc
                              
                            
                            ·
                            Δ
                          
                          ⁢
                          
                            
                          
                          ⁢
                          
                            T
                            doc
                          
                        
                        -
                        
                          Δ
                          ⁢
                          
                            
                          
                          ⁢
                          
                            T
                            dpf
                          
                        
                        -
                        
                          Δ
                          ⁢
                          
                            
                          
                          ⁢
                          
                            T
                            
                              scr
                              +
                              cleanupdoc
                            
                          
                        
                      
                      )
                    
                    .
                  
                
              
            
          
        
      
      
        
          (
          13
          )