Abstract:
In one example, a method of operating an engine having a diesel particulate filter in a vehicle, the particulate filter having a length and depth, includes performing particulate filter regeneration in response to temperature variation across the length and/or depth of the particulate filter.

Description:
BACKGROUND AND SUMMARY 
       [0001]    Diesel powertrains may have a particulate filtration system referred to as a Diesel Particulate Filter (DPF), where engine generated soot may be collected. The collection, or loading, of soot leads to an increase in exhaust pressure, which may degrade engine performance. As such, collected soot can be periodically combusted (e.g., regenerated, or purged) to clean the device and reduce the performance impact. 
         [0002]    It may be advantageous to vary when a particulate filter is regenerated to reduce fuel consumption and extend filter usable life. In some examples, filter soot loading may be inferred and/or correlated to a measure of filter flow restriction, such as based on upstream and/or downstream pressures. However, the restriction over the DPF may depend heavily on the amount of flow, which in turn may vary with temperature in and around the DPF. Further, since temperature may vary both along the length of the filter and/or across the filter width, especially during transients, using a measure or estimate of DPF temperature and/or exhaust temperature may produce errors, especially during low flow conditions (e.g., idle) where errors in models may be amplified. Such errors may lead to unnecessary regeneration, thus increasing fuel usage and decreasing durability. 
         [0003]    Thus, in one approach, the restriction and/or the decision and timing of the regeneration may be correlated to loading taking into account temperature and/or flow distribution along and/or across the DPF. Further, in one embodiment, such correlation may be used during higher flow and/or higher temperature conditions to provide improved accuracy and address the problem of low flow restriction variability. 
         [0004]    The inventors herein have recognized the above issues and approaches, which will be more fully described herein with reference to the description and/or figures. 
     
    
     
       BRIEF DESCRIPTION OF THE FIGURES 
         [0005]      FIG. 1  is a schematic diagram of an engine; 
           [0006]      FIG. 2  is a schematic diagram of exemplary emission control system; 
           [0007]      FIGS. 3 and 7  are example routines for managing particulate filter regeneration; 
           [0008]      FIG. 4  shows an example flow model for a DPF; 
           [0009]      FIGS. 5-6 , and  8  show example filter flow data. 
       
    
    
     DETAILED DESCRIPTION 
       [0010]    Internal combustion engine  10 , comprising a plurality of cylinders, one cylinder of which is shown in  FIG. 1 , is controlled by electronic engine controller  12 . Engine  10  includes combustion chamber  30  and cylinder walls  32  with piston  36  positioned therein and connected to crankshaft  40 . Combustion chamber  30  is shown communicating with intake manifold  44  and exhaust manifold  48  via respective intake valve  52  and exhaust valve  54 . Engine  10  is shown as a direct injection engine with injector  80  located to inject fuel directly into cylinder  30 . Fuel is delivered to fuel injector  80  by a fuel system (not shown), including a fuel tank, fuel pump, and high pressure common rail system. Fuel injector  80  delivers fuel in proportion to the pulse width of signal FPW from controller  12 . Both fuel quantity, controlled by signal FPW and injection timing may be adjustable. Engine  10  may utilize compression ignition combustion under some conditions, for example. 
         [0011]    Controller  12  is shown in  FIG. 1  as a microcomputer including: microprocessor unit  102 , input/output ports  104 , read-only memory  106 , random access memory  108 , and a conventional data bus. Controller  12  is shown receiving various signals from sensors coupled to engine  10 , in addition to those signals previously discussed, including: engine coolant temperature (ECT) from temperature sensor  112  coupled to cooling sleeve  114 ; a measurement of manifold pressure (MAP) from pressure sensor  116  coupled to intake manifold  44 ; a measurement (AT) of manifold temperature from temperature sensor  117 ; an engine speed signal (RPM) from engine speed sensor  118  coupled to crankshaft  40 . 
         [0012]    An emission control system  20  is coupled to an exhaust manifold  48  and several exemplary embodiments of the system in accordance with the present invention are described with particular reference to  FIGS. 2A-2C . 
         [0013]    In one example, engine  10  may be a diesel fueled engine that operates with stratified charge combustion in excess oxygen conditions. Alternatively, fuel timing adjustments, and multiple fuel injections, can be utilized to obtain homogeneous charge compression ignition combustion. While lean operation may be utilized, it is also possible to adjust engine conditions to obtain stoichiometric or rich air-fuel ratio operation. 
         [0014]    In another alternative embodiment, a turbocharger can be coupled to engine  10  via the intake and exhaust manifolds. The turbocharger may include a compressor in the intake and a turbine in the exhaust coupled via a shaft. Further, the engine may include a throttle and exhaust gas recirculation. 
         [0015]    Referring now to  FIG. 2 , the emission control system  20  optionally includes a catalyst system  13  upstream of the particulate filter  15 . Various types of catalysts can be optionally used, such as, for example: a urea based Selective Catalytic Reduction (SCR) catalyst, an oxidation catalyst, and/or a NOx absorber, or these catalysts could be combined with the particulate filter. In the case of an SCR catalyst, in one example, it may include a base metal/zeolite formulation with optimum NOx conversion performance in the range of 200-500° C. Reductant, such as aqueous urea, can be stored on-board and injected in the exhaust system upstream of the SCR catalyst. Alternatively, any other structure known to those skilled in the art to deliver reductant to an exhaust gas aftertreatment device may be used, such as late injection in a direction injection type engine. 
         [0016]    Alternatively, catalyst system  13  may include (separate or in addition to the SCR catalyst) an oxidation catalyst, which may include a precious metal catalyst, preferably one containing platinum, for rapid conversion of hydrocarbons (HC), carbon monoxide (CO) and nitric oxide (NO) in the engine exhaust gas. The oxidation catalyst may also be used to supply heat in the exhaust system, wherein an exotherm is created when extra HC is reduced over the oxidation catalyst. This can be accomplished through, for example, in-cylinder injection during either or both of a power or exhaust stroke of the engine (in a direct injection engine) or any of a number of other alternatives, such as retarding injection timing, increasing EGR and intake throttling, or another approach to increase the HC concentration in the exhaust gas. Alternatively, hydrocarbons may be injected directly into the exhaust gas stream entering the oxidation catalyst. Reductant delivery system  19  may be used to deliver HC from the fuel tank or from a storage vessel to the exhaust system to generate heat for heating the particulate filter  15  for regeneration purposes. 
         [0017]    Particulate filter  15 , in one example a diesel particulate filter (DPF), may be coupled downstream of the catalyst system and may be used to trap particulate matter (e.g., soot) generated during the drive cycle of the vehicle. The DPF can be manufactured from a variety of materials including cordierite, silicon carbide, and other high temperature oxide ceramics. Once soot accumulation has reached a predetermined level, regeneration of the filter can be initiated. Filter regeneration may be accomplished by heating the filter to a temperature that will burn soot particles at a faster rate than the deposition of new soot particles, for example, 400-600° C. In one example, the DPF can be a catalyzed particulate filter containing a washcoat of precious metal, such as Platinum, to lower soot combustion temperature and also to oxidize hydrocarbons and carbon monoxide to carbon dioxide and water. 
         [0018]    Further note that a temperature sensor  21  is shown coupled to the DPF. The sensor, or additional temperature sensors, could also be located within the DPF, or upstream of the filter, or DPF temperature (or exhaust temperature) can be estimated based on operating conditions using an exhaust temperature model. In one particular example, multiple temperature sensors can be used, e.g. one upstream and one downstream of the DPF. 
         [0019]    Also, a differential pressure signal (Δp) is shown being determined from pressure sensors  124  and  126 . Note that a single differential pressure can also be used to measure the differential pressure across DPF  15 . A single port gauge pressure sensor (SPGS) may also be used. In yet another alternative embodiment, the DPF can be located in an upstream location, with an optional catalyst (or catalysts) located downstream. 
         [0020]    As will be appreciated by one skilled in the art, the specific routines described below in the flowcharts may represent one or more of any number of processing strategies such as event-driven, interrupt-driven, multi-tasking, multi-threading, and the like. As such, various acts or functions illustrated may be performed in the sequence illustrated, in parallel, or in some cases omitted. Likewise, the order of processing is not necessarily required to achieve the features and advantages, but is provided for ease of illustration and description. Although not explicitly illustrated, one or more of the illustrated acts or functions may be repeatedly performed depending on the particular strategy being used. Further, these Figures graphically represent code to be programmed into the computer readable storage medium in controller  12 . 
         [0021]    Referring now to  FIG. 3 , a routine is described for controlling particulate filter regeneration, such as based on a determined flow restriction that may be correlated to soot loading. In one example, where restriction over the DPF depends heavily on volumetric flow, which in turn depends on temperature, a distributed correlation may be used. Temperature (and restriction) may be modeled as a distributed quantity over the length of the DPF, especially during transients, rather than a single “lumped” temperature/restriction. However, this is just one approach, and various other may be alternatively used, or used in additional to distributed temperature modeling. 
         [0022]    While improved determinations of the restriction using a distributed approach can provide more appropriately timed DPF regeneration, errors may still persist in the determination. Further, in some cases, such improved estimation approaches may not be used due to timing restrictions, processing power restrictions, system degradation, etc. As such, flow variability may persist. Thus, alternatively, or in addition, the DPF scheduling routine of  FIG. 3  may impose conditions to be satisfied before the updating and/or applying an estimate of flow restriction. Specifically, the routine may limit updating the measured restriction only when post DPF temperature&gt;minimum threshold and exhaust volumetric flow rate&gt;minimum threshold. 
         [0023]    Specifically, first in  310 , the routine reads operating parameters, such as differential pressure temperatures, etc. 
         [0024]    Then, in  312 , the routine determines whether the temperature downstream of the DPF (T_postDPF) is less than a minimum temperature value (T_min) and whether exhaust volumetric flow (Q_v) is less than a minimum flow value (Q_min). Alternatively, the routine may determine whether temperature downstream of the DPF (T_postDPF) is less than a minimum temperature value (T_min) or whether exhaust volumetric flow (Q_v) is less than a minimum flow value (Q_min). 
         [0025]    If the answer to  312  is yes, the routine continues to  314  to freeze the estimated restriction value (R) at its previous vale (R_prev), which would be zero upon initialization. Otherwise, the routine continues to  316  to update the estimated restriction (R) based one or more approaches, such as using a Darcy model, and/or a distributed model as described below with regard to  FIGS. 4-6 . Further, in one particular example, a filter constant, which is a function of the temperature gradient across the particulate filter) may be used to modify a restriction value, where the restriction is based on Darcy&#39;s law, as noted below. Alternatively, the filter constant may be used with a simplified lumped model may also be used 
         [0026]    Next, in  318 , the routine regenerates the DPF based on the determination of  310 , such as by increasing exhaust temperature to a regeneration temperature. The exhaust temperature may be increased, as noted herein, by increasing throttling, late injection, etc. Additionally, the routine may further identify degradation based on the soot loading, such as degradation of the particular filter, and indicate such degradation to an operator, and/or set a code that can be communicated from the vehicle controller. 
         [0027]    While the routine of  FIG. 3  imposes the flow and temperature boundary conditions, various alternative approaches may also be used. For example, such conditions may be avoided in some examples by using alternative estimation techniques, such as described with regard to  FIGS. 7-8 . 
         [0028]    Turning now to  FIGS. 4-7 , information is provided relating to determining flow restriction in a DPF using a distributed approach. Specifically, distributed temperature/flow of the DPF may be modeled as illustrated in  FIG. 4 , and then used to obtain a more accurate restriction correlation. The figure shows an example model of the DPF flow using a circuit analogy. The DPF is approximated as several “parallel” segments exposed to different temperatures, where the temperature of each layer may be a function of inlet and outlet temperature, and each layer is exposed to a fraction of the total flow (e.g., a fraction of the total flow passes through each layer). Darcy&#39;s law may be applied to this system, in which: 
         [0000]        q =( k/ μ)( dp/dx )=( k/ μ)( Δp/x ),  q=v   w   [m/s], k=m   2    
         [0000]      Δ p=xq (μ/ k )= q ( xμ/k ), Δ p =wall pressure drop 
         [0000]      V=IR 
         [0000]        R   w,i   =x   i μ( T   i )/ k    
         [0029]    where: 
         [0030]    q=heat transfer 
         [0031]    k=area 
         [0032]    p=pressure 
         [0033]    x=distance 
         [0034]    μ is viscosity 
         [0035]    R=restriction 
         [0036]    V=flow 
         [0037]    In one example, DPF flow can be assumed to be comprised of multiple (n=number of slices&gt;=2) flow paths. Lumping channel losses with wall losses, the network of  FIG. 4  reduces to a simple parallel network of restrictions (however, in an alternative embodiment, channel losses may be separately modeled). The equivalent restriction can then be calculated as 
         [0000]      1 /R   soot,eq =1/ R   1 +1/ R   2 +1/ R   3 +1/ R   4 + . . . +1/ R   n =Σ(1/ Ri ),, 
         [0000]        R   i =( Δp−c 0− c 2*ρ( Ti )* Qî 2))/ c 1μ( Ti ) Qi,  where μ is viscosity and ρ is density 
         [0000]        Qi =α i *mex h/ ρ( Ti ), mex h =exhaust mass flow rate, α i  belongs to {α 1 , α 2 , α 3 , α 4 , . . . , α n }, Σα i =1, which then dictates the fraction of flow seeing temperature  Ti.    
         [0000]        Ti=T ( i,t ), where, assuming linear temperature drop between  T in and  T out, 
         [0000]        T ( i,t )= T in( t )+( i− 1)( T out( t )− T in( t ))/( n− 1),  i= 1,2,  . . . , n    
         [0038]    The coefficients c0, c1, and c2 can be obtained from experimental flow testing of the DPF. Further, the density and viscosity of exhaust gas can be estimated based on exhaust gas temperature and experimental test data. The restriction R i  is a monotonic function of the soot load in grams/liter. In this way, measured pressure across the DPF can be correlated to a restriction. 
         [0039]    Note that taking a linear temperature profile is without loss of generality. It can be interpreted as adopting a non-linear (e.g. logarithmic) axial grid spacing that samples the linear increments of temperature. Since the results do not (explicitly) depend on axial length, the grid spacing can change dynamically. 
         [0040]      FIG. 6  shows example simulation results for different numbers of slices (2,3) and different distribution of flow at different temperatures (α) The graph shows a normalized restriction (where the restriction is artificially normalized to 1 by dividing by the first value R(0). A comparison metric (100*σ(R)/μ(R)) can be generated for each simulation, which is independent of normalization constant of the restriction and provides a measure of the variance in the restriction. 
         [0000]    
       
         
               
               
               
               
               
             
               
               
               
               
               
             
           
               
                   
                   
               
               
                   
                 Case 
                 num_slices 
                 alpha 
                 Metric (%) 
               
               
                   
                   
               
             
             
               
                   
               
             
          
           
               
                   
                 0 
                 1 
                 N/A* 
                 7.9762 
               
               
                   
                 1 
                 2 
                 [.5 .5] 
                 6.9653 
               
               
                   
                 2 
                 2 
                 [.3 .7] 
                 5.4673 
               
               
                   
                 3 
                 3 
                 [.1 .1 .8] 
                 6.2697 
               
               
                   
                 4 
                 3 
                 [.33 .01 .66] 
                 5.7001 
               
               
                   
                 5 
                 3 
                 [.45 .1 .45] 
                 7.5119 
               
               
                   
                 6 
                 3 
                 [.01 .33 .66] 
                 6.4144 
               
               
                   
                 7 
                 3 
                 [.33 .33 .34] 
                 8.1006 
               
               
                   
                 8 
                 3 
                 [.45 .33 .22] 
                 8.7540 
               
               
                   
                 9 
                 3 
                 [.1 .45 .45] 
                 7.0568 
               
               
                   
                 10 
                 3 
                 [.33 .45 .22] 
                 8.2187 
               
               
                   
                 11 
                 3 
                 [.45 .45 .1] 
                 8.8337 
               
               
                   
                   
               
             
          
         
       
     
         [0041]    As indicated above, it may be possible to further improve the correlation between flow restriction and loading by including channel losses. As indicated below, by including the channel losses, the restriction calculation may be further stabilized. Specifically, the correlation may be modified to include the temperature offset between gas at the end of the inlet channel and the exit temperature (t_postdpf). The offset exists for heat transfer to occur at the wall towards dpf outlet, and affects the results as shown in the simulation data of  FIG. 6  and table below. 
         [0042]    Specifically, the channel losses may be modeled as: 
         [0000]        dΔp= 4 f* ( dx   —   dpf/D _cell)*(½) ρV ( x   —   dpf )̂2, where  f= 64/ Re  (Reynolds number), 
         [0043]    since channel Re&lt;2000 almost always, and 0&lt;x_dpf&lt;L_dpf, and Δpchannel_loss=∫dΔp. 
         [0044]    The temperature offset between gas at the end of the inlet channel and the exit temperature (t_postdpf) may be modeled as: 
         [0000]        T ( i,t )= T in( t )+( i− 1)( T _offset( t )+ T out( t )− T in( t ))/( n− 1), where  I= 1, 2,  . . . , n, A  is a constant offset of 75 deg C. 
         [0000]    The simulation data including channel losses used n=2 slices with alpha=[0.3 0.7], and again R artificially normalized to 1, with the same comparison metric. As indicated, including the channel losses further improved the correlation. 
         [0000]    
       
         
               
               
               
               
               
             
           
               
                   
                   
               
               
                   
                 Case 
                 t_offset_c 
                 channel loss 
                 Metric (%) 
               
               
                   
                   
               
             
             
               
                   
                 0 
                 N/A 
                 N/A 
                 7.9762 
               
               
                   
                 1 
                 0 
                 Not incl. 
                 5.4673 
               
               
                   
                 2 
                 75 
                 Not incl. 
                 4.4815 
               
               
                   
                 3 
                 75 
                 included 
                 3.8158 
               
               
                   
                   
               
             
          
         
       
     
         [0045]    Referring now to  FIGS. 7-8 , still another approach to determining filter restriction is provided. In particular, the approach may include applying a Darcy model to gas flows for lower flow soot load estimation. Specifically, variability at lower flows can be avoided during in some examples by avoiding determining the restriction estimate at such conditions (as noted in  FIG. 3 ) and using a lumped parameter estimate in the example where temperature gradients may be an in-significant noise factor for high-flow restriction variability. However, during lower flow conditions such as idle or tip-out, the restriction is frozen and as such may not be updated for a significant duration, depending on the vehicle drive cycle, such as during an extended idle. 
         [0046]    Thus, in still another approach, a first mapping approach (which may include a first estimation routine) can be used during higher flow and higher temperature conditions, and an alternative mapping may be used during lower flow (and/or lower temperature) conditions, such as described with regard to  FIG. 7 . For example, the approach of  FIGS. 4-6  may be used during the higher conditions, or other mappings may be used, as noted below. 
         [0047]    Specifically,  FIG. 7  shows an example flow chart of a routine that may be used. During higher flow conditions, a lumped approach may be used. However, during lower flow conditions, a mapping using Darcy&#39;s law and a linear transformation may be used, where: 
         [0000]    
       
         
           
             
               
                 R 
                 D 
               
               = 
               
                 
                   ( 
                   
                     1 
                     
                       
                         p 
                         1 
                       
                        
                       
                         Q 
                         exh 
                       
                        
                       μ 
                     
                   
                   ) 
                 
                  
                 
                   
                     ( 
                     
                       
                         p 
                         1 
                         2 
                       
                       - 
                       
                         p 
                         2 
                         2 
                       
                     
                     ) 
                   
                   
                     2 
                      
                     
                       R 
                       
                         D 
                          
                         
                             
                         
                          
                         0 
                       
                     
                   
                 
               
             
             , 
             
               
                 with 
                  
                 
                     
                 
                  
                 R 
                  
                 
                     
                 
                  
                 d 
                  
                 
                     
                 
                  
                 0 
               
               = 
               
                 1.44 
                  
                 
                     
                 
                  
                 e 
                  
                 
                     
                 
                  
                 5 
               
             
           
         
       
     
         [0048]    where p1 is upstream DPF pressure, p2 is downstream DPF pressure, and Qexh is exhaust flow, and: 
         [0049]    the linear transformation follows the equation, R=RD*1.79−0.78. Note that the value Rd0 may vary with system component specifications. In particular, Rd0 represents a normalization constant Rd for a clean (no soot) DPF to a fixed number, e.g. 1. This also applies to the linear transformation. 
         [0050]    Then, during higher flow conditions, the restriction may be determined as: 
         [0000]        R =(Δ p−c 0− c 2*ρ( T )* Q exĥ2))/ c 1μ( T ) Q exh. 
         [0051]    Specifically, referring to  FIG. 7 , in  710 , the routine reads operating parameters, such as differential pressure, temperatures, flow (e.g. engine flow), etc. 
         [0052]    Then, in  712 , the routine determines whether the temperature downstream of the DPF (T_postDPF) is less than a minimum temperature value (T_min). If so, the routine continues to  714  to freeze the estimated restriction value (R) at its previous vale (R_prev), which would be zero upon initialization. Otherwise, the routine continues to  716  to determine whether exhaust volumetric flow (Q_v) is less than a minimum flow value (Q_min). If so, the routine continues to  718  to determine RD using the above second mapping, and then in  720 , transforms RD to R using the linear transformation. Alternatively, when the answer to  716  is no, the routine continues to  722  to determine the restriction (R) based on a first mapping, such as illustrated above using parameters c0, c1, and c2. From either  722 ,  720 , or  714 , the routine continues to  724  to determine whether to regenerate the DPF based on R and various operating conditions, such as vehicle speed, ambient temperature, desired engine output, etc. 
         [0053]    Referring now to  FIG. 8 , data illustrates results for the approach illustrated in  FIG. 7 , where different sized restrictions are shown versus mass flow for a 6.4 L engine having a 9″×12″ sized DPF. The data illustrate that the combination of different estimation routines for different flow and/or temperature conditions can provide improved estimation results for determining flow restriction of a DPF in an engine exhaust. 
         [0054]    It will be appreciated that the configurations and routines disclosed herein are exemplary in nature, and that these specific embodiments are not to be considered in a limiting sense, because numerous variations are possible. The subject matter of the present disclosure includes all novel and nonobvious combinations and subcombinations of the various systems and configurations, and other features, functions, and/or properties disclosed herein. 
         [0055]    The following claims particularly point out certain combinations and subcombinations regarded as novel and nonobvious. These claims may refer to “an” element or “a first” element or the equivalent thereof. Such claims should be understood to include incorporation of one or more such elements, neither requiring nor excluding two or more such elements. Other combinations and subcombinations of the disclosed features, functions, elements, and/or properties may be claimed through amendment of the present claims or through presentation of new claims in this or a related application. Such claims, whether broader, narrower, equal, or different in scope to the original claims, also are regarded as included within the subject matter of the present disclosure.