Abstract:
A method for controlling an internal combustion engine using a controller that controls an air flow path by adjusting at least one of a variable geometry turbine (VGT) and an exhaust gas recirculation (EGR) flow rate during engine operation. The method determines inputs, such as engine speed and fuel rate from the sensor data, and employs a switch based gain-scheduled explicit model predictive controller (MPC) responsive to the inputs to determine the air flow path.

Description:
BACKGROUND 
       [0001]    1. Field of the Disclosure 
         [0002]    This application relates to an airflow control into an internal combustion engine. More particularly, the application relates to performing an optimal control of an internal combustion engine using an explicit model predictive controller in combination with a switch gain scheduler. 
         [0003]    2. Description of the Related Art 
         [0004]    Air flow control of an internal combustion engine is an important strategy to increase engine performance and fuel economy. This is particularly true with diesel engines. Traditionally, an air flow control strategy is developed using a proportional (P), integral (I), and differential (D) control design. PID control design is focused on nominal engine operating conditions and consequently does not perform effectively under off-nominal conditions. PID is a well-known simple control approach utilized for single input single output (SISO) systems. PID approaches have been extended to multiple-input multiple output (MIMO) systems; however more advanced control techniques are able to deliver better performance, e.g. Model Predictive Controller (MPC). MPCs are popular in the process industry such as in chemical processing technologies and are currently being adopted with limited use in the automotive industry. Internal combustion engines operate under different conditions (e.g., high speed or load vs low speed or load) during their use. Each engine operating condition has an air flow dynamic that necessitates design of multiple controllers. Design and use of multiple controllers trigger problems related to the limited memory of an engine control unit (ECU) and the cost of controller calibration. 
       SUMMARY 
       [0005]    In order to reduce the memory usage and processing time, a Switch Gain Scheduled controller is designed and implemented in conjunction with an explicit MPC. 
         [0006]    According to an embodiment, a switch gain scheduler used in conjunction with an explicit MPC makes an off-nominal operating condition resemble a nominal operating condition. The switch gain scheduler changes a controller gain depending on operating conditions. 
         [0007]    The forgoing general description of the illustrative implementations and the following detailed description thereof are merely exemplary aspects of the teachings of this disclosure, and are not restrictive. 
     
    
     
       BRIEF DESCRIPTION OF THE DRAWING 
         [0008]    A more complete appreciation of this disclosure and many of the attendant advantages thereof will be readily obtained as the same becomes better understood by reference to the following detailed description when considered in connection with the accompanying drawings, wherein: 
           [0009]      FIG. 1  is a block diagram of the switch gain scheduled explicit model predictive control (MPC) showing an off-nominal plant model coupled with the switch gain scheduler to resemble a nominal plant model; 
           [0010]      FIG. 2  is an exemplary process of the switch gain scheduled explicit model predictive control; 
           [0011]      FIG. 3  is block diagram of a standard observer-based feedback controller, where the plant can either be under the nominal condition or under off-nominal operating conditions. 
           [0012]      FIG. 4  is an exemplary illustration of the constraints of applying gain scheduling directly; 
           [0013]      FIG. 5A and 5B  are an exemplary illustration of the solution proposed in this patent to the constraints of switch gain scheduled problem of  FIG. 4 ; 
           [0014]      FIG. 6  is schematic representation of an ECU embedded with the switch gain scheduled explicit MPC controller, whose control actions are send to the engine via actuator. 
           [0015]      FIG. 7  is a graph depicting the difference in the engine parameter response in case of full gain scheduling and “switched” gain scheduling. The explicit MPC is built using 1600 rpm, and simulated against 4000 rpm model. 
           [0016]      FIG. 8  is a graph depicting the difference in the engine parameter response in case of full gain scheduling, “switched” gain scheduling, and diagonal gain scheduling. 
           [0017]      FIG. 9  is a graph depicting the non-linear model simulation of intake manifold pressure (MAP) using “switched” gain scheduling. 
           [0018]      FIG. 10  is a graph depicting the non-linear model simulation of exhaust gas recirculation rate (EGR) using “switched” gain scheduling. 
           [0019]      FIG. 11  is a graph illustrating the sample switch gain scheduling matrix to be used at each time-step during the example simulation. 
       
    
    
     DETAILED DESCRIPTION 
       [0020]    In the drawings, like reference numerals designate identical or corresponding parts throughout the several views. Further, as used herein, the words “a”, “an” and the like generally carry a meaning of “one or more”, unless stated otherwise. The drawings are generally drawn to scale unless specified otherwise or illustrating schematic structures or flowcharts. 
         [0021]    Furthermore, the terms “approximately,” “proximate,” and similar terms generally refer to ranges that include the identified value within a margin of 20%, 10% or preferably 5%, and any values there between. 
         [0022]      FIG. 1  is a conceptual block diagram of a switch gain scheduled explicit MPC. The switch gain scheduled block  102  is coupled with an off-nominal plant  103  and together they represent a nominal operating plant block  104 . The explicit MPC 101  is the primary controller, which provides an optimal control action to the plant in off-nominal condition via the switch gain scheduled block  102 . The optimal control action is determined based on the difference between a reference state x ref  and an estimated state {circumflex over (x)} from an observer  105 . The observer  105  implements an algorithm that calculates an estimated plant state from the observed output y of the off-nominal plant and the commanded plant input u*. 
         [0023]    The advantage of this embodiment is only one explicit model predictive controller (MPC) may be designed. Hence, the ECU memory usage is significantly reduced in addition to the calibration requirement reduction, while maintaining the performance of the engine such as the diesel engine. 
         [0024]    Referring to  FIG. 2 , the switch gain scheduled explicit MPC design and implementation and can involve several offline and online processes. First process step  201  is an offline calibration of an explicit MPC near a nominal operating condition. Examples of nominal operating conditions include an engine operating at idle speed, say 800 rpm. Second process step  202  is also offline and involves determination of a gain scheduling matrix for different operating conditions. The first two steps are independent and commonly performed during the control system design. The process steps after  202  determine the optimal control for different operating conditions by applying gain scheduling along with an explicit MPC, defined for nominal operating conditions. In process step  203 , an offline process, the gain scheduling matrices H i  are split into diagonal and off-diagonal parts H i   d  and H i   o . 
         [0025]    In process step  204 , constraint parameters for the scaled plant with H i   d  are computed in real-time online, where i is the operating condition. Next, in process step  205 , the diagonal matrix H i   d  is used to determine the optimal control strategy. The optimal control strategy is obtained using the using explicit MPC solver and scaled constraints parameters. In process step  206 , the diagonal matrix H i   d  is multiplied to the optimal control action taken by the explicit MPC to calculate a new controlled plant input u. Then the switch gain scheduled explicit MPC waits for the next half step sample time  207 . 
         [0026]    In the process step  208 , constraint parameters for the scaled plant with H i   o  are computed in real-time, where i is the operating condition. In the process step  209 , the diagonal matrix H i   o  is used to determine the optimal control strategy. The optimal control strategy is obtained using the explicit MPC solver and scaled constraints parameters. In the process step  210 , the off-diagonal matrix H i   o  is multiplied to the optimal control action taken by the explicit MPC to calculate a new controlled plant input u. Then the switch gain scheduled explicit MPC waits for the next half step sample time  211 . The processes  204 - 211  are repeated as the time advances. 
         [0027]    The switch gain scheduled explicit MPC is general enough to be applied to any internal combustion engine control system application. In this embodiment, diesel engine air flow control is used as a sample application. 
         [0028]    Referring to  FIG. 3 , consider the nominal plant, P 0 (s)  302 , which is implemented as a set of linearized equations and off-nominal plant, P i (s)  305 , which is implemented as a set of linearized equations. P 0 (s) and P i (s) result from linearization of a non-linear plant model operating at different operating conditions. The observers  303  and  306  implement the same algorithm, but depending on the control action u* and the state of the plant, the observers provide different output to their respective explicit MPC controller. 
         [0029]    Each of the nominal and off-nominal plants may be controlled by different explicit MPCs  301  and  304 , respectively. In case of off-nominal conditions, if there are i different operating condition, then i different MPCs may be implemented, which in turn increases the memory storage and processing requirement of the ECU  601  used for the engine  604  in  FIG. 6 . Examples of off-nominal operating conditions include low speed or high speed engine condition like 800 rpm or 4000 rpm respectively. Implementing a plurality of explicit MPC in the ECU  601  of an automobile is generally discouraged due to the limited memory and processing capacity of the ECU. 
         [0030]    In the present embodiment, the plurality of MPC issue is resolved by designing a single explicit MPC at nominal operating condition and inserting the switching gain scheduler block  102 . The switch gain scheduler  102  is coupled with the off-nominal plant model  103  as illustrated in  FIG. 1 . The coupled blocks within  104 , together imitate the behavior of a plant operating under nominal conditions. Hence, only a single explicit MPC operating at a nominal condition needs to be designed. 
         [0031]    The design and implementation of a gain scheduler itself is dependent on the type of controller and system dynamics, which include different operating conditions, that render the gain scheduler a specialized module. The present embodiment designs and implements the gain scheduler for an explicit MPC type of controller. In addition, the switch gain scheduler is designed based on splitting the gain matrix H i  into H i   d  and H i   o , where i refers to the operating condition, d denotes the diagonal matrix and o denotes the off-diagonal matrix. The switch gain scheduler defines a switching variable which dynamically switches between the H i   d  and the H i   o  gain matrices. 
         [0032]    The basis of splitting the gain matrix H i  is illustrated in  FIG. 4 . Referring to  FIG. 4 , the axis u 1  and u 2  represent the two inputs (e.g., amount of fuel, engine speed) to the controller. The bounded region  401  is an exemplary true control constraint, which must be satisfied by any plant controller. The optimal control (u*) applied by the explicit MPC  401  to a plant falls within the true control constraint  401 . Thus the optimal control actions (u*) issued by the explicit MPC controller to the plant should satisfy the true control constraints. 
         [0033]    Applying the gain scheduling (GS) matrix H i  in the traditional way transforms the control issued by the explicit MPC. An exemplary transformation  402  shows that the original control constraints  401  are violated (shaded regions). 
         [0034]    The mathematical form of the exemplary  FIG. 4  is equations (1) and (2). Equation (1) represents the situation when only explicit MPC controls the plant. In equation (1), V is the coefficient matrix of the true control constraint and cannot be changed; u=u* is the plant input matrix; W is the real control constraint matrix as seen by the explicit MPC. 
         [0000]        V u≦W    (1)
 
         [0035]    Where, 
         [0000]    
       
         
           
             
               V 
               = 
               
                 [ 
                 
                   
                     
                       1 
                     
                     
                       0 
                     
                   
                   
                     
                       
                         - 
                         1 
                       
                     
                     
                       0 
                     
                   
                   
                     
                       0 
                     
                     
                       1 
                     
                   
                   
                     
                       0 
                     
                     
                       
                         - 
                         1 
                       
                     
                   
                 
                 ] 
               
             
             ; 
             
               u 
               = 
               
                 [ 
                 
                   
                     
                       
                         u 
                         1 
                       
                     
                   
                   
                     
                       
                         u 
                         2 
                       
                     
                   
                 
                 ] 
               
             
             ; 
             
               W 
               = 
               
                 [ 
                 
                   
                     
                       
                         u 
                         
                           max 
                            
                           
                               
                           
                            
                           1 
                         
                       
                     
                   
                   
                     
                       
                         - 
                         
                           u 
                           
                             min 
                              
                             
                                 
                             
                              
                             1 
                           
                         
                       
                     
                   
                   
                     
                       
                         u 
                         
                           max 
                            
                           
                               
                           
                            
                           2 
                         
                       
                     
                   
                   
                     
                       
                         - 
                         
                           u 
                           
                             min 
                              
                             
                                 
                             
                              
                             2 
                           
                         
                       
                     
                   
                 
                 ] 
               
             
           
         
       
     
         [0036]    On applying the gain scheduling after the explicit MPC, equation (1) is transformed to equation (2). The resulting matrix does not fully satisfy the constraints W. 
         [0000]        V H   i   u*≦W    (2)
 
         [0037]    Where, 
         [0000]    
       
         
           
             
               H 
               i 
             
             = 
             
               [ 
               
                 
                   
                     
                       α 
                       11 
                     
                   
                   
                     
                       α 
                       12 
                     
                   
                 
                 
                   
                     
                       α 
                       21 
                     
                   
                   
                     
                       α 
                       22 
                     
                   
                 
               
               ] 
             
           
         
       
     
         [0038]    When working with explicit MPC, the control constraints restrict one to use special cases of H i  such that the transformation u=H i  u* (refer  FIG. 1 ). In this embodiment, the gain scheduling (GS) is designed such that the original control constraints that are satisfied by explicit MPC are not violated. 
         [0039]    Referring to  FIGS. 5A and 5B , the gain scheduling matrix H i  is spilt into two matrices H i   d —composed of the diagonal element of H i  and H i   o —composed of the diagonal element of H i . The gain scheduling algorithm switches between H i   d  and H i   o  based on the sampling time T s . The H i   d  matrix is applied at an even time step, while H i   d  matrix is applied at an odd time step. If the sampling time of the explicit MPC is assumed as 2T s , then on average the desired scheduling H i  is achieved. 
         [0040]    Referring to  FIG. 5A , the control constraints as seen by the explicit MPC (W′)  501  are modified and are different from the real control constraint (W)  401 . On applying the diagonal gain scheduling H i   d , the real constraints are satisfied. 
         [0041]    Referring again to  FIG. 5A , the mathematical form of constraints seen by explicit MPC can be written as equation (3): 
         [0000]        V u≦W′    (3)
 
         [0042]    Where, 
         [0000]    
       
         
           
             
               
                 W 
                 ′ 
               
               = 
               
                 [ 
                 
                   
                     
                       
                         
                           α 
                           11 
                           
                             - 
                             1 
                           
                         
                          
                         
                           u 
                           
                             max 
                              
                             
                                 
                             
                              
                             1 
                           
                         
                       
                     
                   
                   
                     
                       
                         
                           - 
                           
                             α 
                             11 
                             
                               - 
                               1 
                             
                           
                         
                          
                         
                           u 
                           
                             min 
                              
                             
                                 
                             
                              
                             1 
                           
                         
                       
                     
                   
                   
                     
                       
                         
                           α 
                           22 
                           
                             - 
                             1 
                           
                         
                          
                         
                           u 
                           
                             max 
                              
                             
                                 
                             
                              
                             2 
                           
                         
                       
                     
                   
                   
                     
                       
                         
                           - 
                           
                             α 
                             22 
                             
                               - 
                               1 
                             
                           
                         
                          
                         
                           u 
                           
                             min 
                              
                             
                                 
                             
                              
                             2 
                           
                         
                       
                     
                   
                 
                 ] 
               
             
             ; 
           
         
       
     
         [0000]    V and u are same as in equation (1) 
         [0043]    On applying the gain scheduling after the explicit MPC, equation (3) is transformed to equation (4). The resulting matrix now fully satisfies the constraints W  401 . 
         [0000]        V H   i   d    u*≦W    (4)
 
         [0044]    Where, 
         [0000]    
       
         
           
             
               
                 H 
                 i 
                 d 
               
               = 
               
                 [ 
                 
                   
                     
                       
                         α 
                         11 
                       
                     
                     
                       0 
                     
                   
                   
                     
                       0 
                     
                     
                       
                         α 
                         22 
                       
                     
                   
                 
                 ] 
               
             
             ; 
           
         
       
     
         [0000]    V and u* are same as in equation (1) 
         [0045]    Referring to  FIG. 5B , the mathematical form of constraints seen by explicit MPC can be written as equation (5): 
         [0000]        V u≦W″   (5)
 
         [0046]    Where, 
         [0000]    
       
         
           
             
               
                 W 
                 ″ 
               
               = 
               
                 [ 
                 
                   
                     
                       
                         
                           α 
                           12 
                           
                             - 
                             1 
                           
                         
                          
                         
                           u 
                           
                             max 
                              
                             
                                 
                             
                              
                             1 
                           
                         
                       
                     
                   
                   
                     
                       
                         
                           - 
                           
                             α 
                             12 
                             
                               - 
                               1 
                             
                           
                         
                          
                         
                           u 
                           
                             min 
                              
                             
                                 
                             
                              
                             1 
                           
                         
                       
                     
                   
                   
                     
                       
                         
                           α 
                           21 
                           
                             - 
                             1 
                           
                         
                          
                         
                           u 
                           
                             max 
                              
                             
                                 
                             
                              
                             2 
                           
                         
                       
                     
                   
                   
                     
                       
                         
                           - 
                           
                             α 
                             21 
                             
                               - 
                               1 
                             
                           
                         
                          
                         
                           u 
                           
                             min 
                              
                             
                                 
                             
                              
                             2 
                           
                         
                       
                     
                   
                 
                 ] 
               
             
             ; 
           
         
       
     
         [0000]    V and u are same as in equation (1) 
         [0047]    On applying the gain scheduling after the explicit MPC, equation (5) is transformed to equation (6). The resulting matrix now fully satisfies the constraints W  401 . 
         [0000]        V H   i   o    u*≦W    (6)
 
         [0048]    Where, 
         [0000]    
       
         
           
             
               H 
               i 
               o 
             
             = 
             
               [ 
               
                 
                   
                     0 
                   
                   
                     
                       α 
                       12 
                     
                   
                 
                 
                   
                     
                       α 
                       21 
                     
                   
                   
                     0 
                   
                 
               
               ] 
             
           
         
       
     
         [0049]    A mathematical proof illustrating why a gain scheduling method that switches between a gain matrix H i   d  and H i   o  works as well as a full gain scheduling, in which there is no switching between gain matrix, is discussed hereafter. Recall, a switched gain scheduled matrix also satisfies the real control constraints as opposed to a full gain scheduler (see  FIG. 4 ). 
         [0050]    A discretized dynamic system employing a full gain scheduling can be represented mathematically as equations 7 and 8. 
         [0000]        x   k+1   =A x   k   +B H u   k    (7)
 
         [0000]        x   k+2   =A   2   x   k   +A B H u   k    (8)
 
         [0000]    Where, A and B are system dynamics matrices, H is a full gain scheduling matrix, x k  is the system input at time-step k, and u k  is the controlled plant input at time-step k. 
         [0051]    A dynamic system employing a switched gain scheduling can be represented mathematically as equations 9 and 10. 
         [0000]          x     k+1   =A  x     k   +B H   d    u   k    (9)
 
         [0000]          x     k+2   =A   2      x     k   +A B H   d    u   k   +B H   o    u   k    (10)
 
         [0052]    An error system may be defined as e k =x k −  x   k . Then, the error dynamics between a dynamic system employing full matrix scheduling and a dynamic system employing a switched system can be written as equation 11: 
         [0000]        e   k+2   =A   2    e   k +( A B H−A B H   d   −B H   o )  u   k    (11)
 
         [0053]    As long as the eigenvalues of A are inside the unit circle, and the control u k  is bounded, the error can also be bounded. Now consider a first order Taylor series expansion of an exponential matrix as given by equation 12. 
         [0000]        e   A     c     ΔT   =I+A   c   ΔT+ higher order terms   (12)
 
         [0054]    Applying (12) to the discretized equation 11 and assuming H=H d +H o  we get, 
         [0000]        e   k+2   =A   2    e   k   +ΔT   2    A B H   o    u   k +higher order terms of ΔT   (13)
 
         [0055]    Note in equation 13, for small sampling period ΔT, the term ΔT 2  A B H o  u k  dominates over the higher order terms of ΔT. As ΔT tends to zero, a switch gain scheduled system approaches the dynamic system under consideration. 
         [0056]    The results of sample implementation of the switch gain scheduled explicit MPC for a diesel engine air flow control are shown in  FIG. 7-11 .  FIG. 7  shows a linear simulation showing the use of standard full matrix gain scheduling (GS) versus the switched gain scheduling strategy (SGS). As can be seen, the switching strategy is able to approximately recover the input-output response of using full matrix gain scheduling. 
         [0057]      FIG. 8  shows the zoomed in version of simulation in  FIG. 7 , in addition the purely diagonal gain scheduling is depicted for comparison.  FIG. 8  clearly depicts the difference in responses between purely diagonal gain scheduling (DGS) and switched gain scheduling. The switched gain scheduling strategy not only guarantees that the control action is within the real control limits, but also closely follows the full matrix gain scheduling behavior, which does not necessarily satisfy the real control constraints. 
         [0058]      FIGS. 9-11  illustrate the intake pressure, EGR rate responses, and gain scheduling parameters over a portion of the NEDC drive cycle. In particular,  FIG. 9  illustrates the effect of the switched gain scheduling strategy on the intake pressure. The response time of the intake pressure is faster by approximately 100 ms on tip-ins compared to the full gain scheduler. A faster intake pressure response may be desirable for faster acceleration and better drivability in general. 
         [0059]      FIG. 10  illustrates the effect of a switched gain scheduling strategy on exhaust gas recirculation (EGR) rate measured in gram/sec. The EGR rate closely follows the desired reference, and output is comparable to one obtained from the full gain scheduled controller. In general, a faster EGR response is one of the desired characteristics of internal combustion engine control in order to meet emissions requirements. A higher EGR rate may lead to a lower NOx percentage. 
         [0060]      FIG. 11  depicts the diagonal and off-diagonal gain matrix used at each time step during the sample simulation discussed earlier. Referring back to  FIG. 6 , once the switch gain scheduled explicit MPC controller  602  is designed as summarized in  FIG. 2  and detailed herein, all the related algorithms are implemented in the engine control unit (ECU)  601  of engine  604 . 
         [0061]    An exemplary engine control unit (ECU) such as ECU  601  contains at least one micro-processor or the equivalent, such as a central processing unit (CPU) or application specific processor ASP (not shown), input and output interfaces  607 , memory circuit  606  (e.g., ROM, EPROM, EEPROM, flash memory, static memory, DRAM, SDRAM, and their equivalent), power circuitry, and other supporting components. The microprocessor is circuitry that utilizes a computer readable storage medium, such as the memory circuit  606 , configured to control the microprocessor to perform and/or control the processes discussed in this embodiment. 
         [0062]    The microprocessor or aspects thereof, in alternate implementations, can include or exclusively include a logic device for augmenting or fully implementing this disclosure. Such a logic device includes, but is not limited to, an application-specific integrated circuit (ASIC), a field programmable gate array (FPGA), a generic-array of logic (GAL), and their equivalent. The microprocessor can be a separate device or a single processing mechanism. Further, this disclosure can benefit from parallel processing capabilities of a multi-cored CPU. Control circuitry provided by one or more processors in multi-processing arrangement may also be employed to execute sequences of instructions contained in memory. Alternatively, hard-wired circuitry may be used in place of or in combination with software instructions. The exemplary implementations discussed herein are not limited to any specific combination of hardware circuitry and software.