# Patent

## Title
System and method for model based boost control of turbo-charged engines

## Abstract

A system and method for controlling boost pressure in various turbo-charged engine configurations as well as variable geometry turbine (VGT) arrangements includes an electronic controller programmed to receive a predetermined desired boost pressure PBoostdes. A desired pressure delta ΔPWGdes across a waste-gate valve is determined using the desired boost pressure PBoostdes. A control signal is generated for controlling the waste-gate valve so as to achieve the desired pressure delta ΔPWGdes. In boost pressure and vacuum pneumatically-actuated waste-gate valve arrangements, the respective solenoid duty cycles are obtained through use of various data structures. Where a waste-gate valve position is controlled by an electrical motor, the valve position is determined using a data structure as a function of desired waste-gate valve flow at sonic standard conditions.

## Background
RELATED APPLICATIONS
This application is related to co-pending application entitled “SYSTEM AND METHOD FOR MODELING TURBO-CHARGED ENGINES AND INDIRECT MEASUREMENT OF TURBINE AND WASTE-GATE FLOW AND TURBINE EFFICIENCY,” Ser. No. 11/867,422 owned by the common assignee of the present invention and herein incorporated by reference in its entirety.
TECHNICAL FIELD
The present invention relates to a system and method for model based boost control of turbo-charged engines.
BACKGROUND OF THE INVENTION
A turbo-charged internal combustion engine includes additional components and physical processes in both the intake and exhaust stream. On the intake side of the engine, a centrifugal compressor and intercooler are provided and are located between the air cleaner and a throttle valve. On the exhaust side, a turbine and a waste-gate—which defines a parallel exhaust stream path with the turbine—are both located between the exhaust manifold and the catalyst/muffler. It is known to provide an engine management system (EMS) configured to control the operation of a turbo-charged engine, including boost control. However, such an EMS is conventionally configured to perform its functions with only a minimal amount of additional information, notwithstanding the increased system complexity, in order to maintain reduced costs (i.e., by reducing the number of sensors). Conventionally, the additional sensors added when an engine is turbo-charged are all located on the intake side (e.g., a boost pressure sensor and boost temperature sensor).
As to boost control, conventional systems do not adequately address efficiency considerations. More specifically, conventional controls do not adequately address the desire for producing the optimum amount of boost for a needed amount of engine torque. As a result, the turbine presents an unnecessarily large restriction to the engine exhaust, which decreases efficiency. Also, the compressor ends up producing too much boost, which must thereafter be decreased downstream by the intake throttle, also resulting in operating inefficiency. This unnecessarily decreases fuel economy. Also, conventional control systems are specific to a particular implementation type, which reduces its applicability to other waste-gate mechanization approaches.
There is therefore a need for a system and method for controlling the boost pressure to a turbo-charged engine that minimizes or eliminates one or more of the problems set forth above.
SUMMARY OF THE INVENTION
The present invention provides for increases engine efficiency. The invention enables accurate boost control. This in turn allows setting the desired boost to an optimum, desired boost needed to satisfy requested engine torque without any intake throttling. This results in increased efficiency. Moreover, the invention incorporates the flexibility to work with a variety of waste-gate mechanizations.
A method is provided of controlling the boost pressure to an internal combustion engine having a turbo-charger with a compressor and an exhaust driven turbine with a parallel-path waste-gate adjusted by a waste-gate valve. The method includes a number of steps. In a preferred embodiment, the first step involves determining a desired waste-gate mass flow rate at standard sonic conditions {dot over (m)}WGdes@SonicStdCond across the waste-gate valve using a predetermined (i.e., given) desired boost pressure PBoostdes. This step is common across a wide variety of waste-gate mechanizations, resulting in increased flexibility. The method also involves the step of generating a control signal for controlling the waste-gate valve so as to achieve the desired flow rate {dot over (m)}WGdes@SonicStdCond. In alternate embodiments, the control signal may be customized to various pneumatically-actuated and electrically-actuated waste-gate valve mechanizations.
In an alternate embodiment, one use for the method involves determining a desired boost pressure PBoostdes and a desired engine air mass flow rate {dot over (m)}eng,airdes based on a requested engine torque. As above, this allows for increased engine operating efficiency. Of course, the method for controlling boost pressure may be used for many other applications, and is dependent on just a given value for the desired boost pressure PBoostdes as an input.
Other features, object and advantages of the present invention are also presented.



BRIEF DESCRIPTION OF THE DRAWINGS
The present invention will now be described by way of example, with reference to the accompanying drawings:
FIG. 1 is simplified diagrammatic and block diagram of a turbo-charged engine system having a controller configured for model-based boost control.
FIG. 2 is a flow chart diagram showing the boost control method of the present invention.
FIGS. 3-5 are diagrammatic views of boost-based pneumatic, vacuum-based pneumatic, and electrically actuated waste-gate valve mechanizations.
FIG. 6 is a block diagram of one embodiment of model-based boost control.
FIG. 7 is a diagram showing turbine isentropic efficiency as a function of turbine pressure ratio.
FIG. 8 is a diagram showing turbine pressure ratio as a function of corrected engine flow.
FIG. 9 is a chart of data in a data table showing the desired waste-gate diaphragm delta pressure.
FIG. 10 is a diagrammatic view showing a calibration setup for obtaining data for a duty cycle data table.



DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT
Referring now to the drawings wherein like reference numerals are used to identify identical components in the various views, FIG. 1 is a diagrammatic view of a turbo-charged internal combustion engine system 10 configured in accordance with the present invention. The system 10 includes an internal combustion engine 12 controlled by an electronic engine controller 14 all in accordance with the present invention.
Engine 12 may be a spark-ignition engine that includes a number of base engine components, sensing devices, output systems and devices, and a control system. Alternatively, the present invention may be used with compression-ignition engines, such as diesel or the like.
Generally, electronic controller 14 is configured via suitable programming to contain various software algorithms and calibrations, electrically connected and responsive to a plurality of engine and vehicle sensors, and operably connected to a plurality of output devices. Controller 14 includes at least one microprocessor or other processing unit, associated memory devices such as read only memory (ROM) 14a and random access memory (RAM) 14b, input devices for monitoring input from external analog and digital devices, and output drivers for controlling output devices. In general, controller 14 is operable to monitor engine operating conditions and operator inputs using the plurality of sensors, and control engine operations with the plurality of output systems and actuators, using pre-established algorithms and calibrations that integrate information from monitored conditions and inputs. The software algorithms and calibrations which are executed in electronic controller 14 may generally comprise conventional strategies known to those of ordinary skill in the art. These programmed algorithms and calibrations are configured, when executed, to monitor the engine operating conditions and operator demands using the plurality of sensors, and control the plurality of engine actuators accordingly. The software algorithms and calibrations are preferably embodied in pre-programmed data stored for use by controller 14.
System 10 further includes a turbo-charger 15 having a compressor 16, which may include a compressor recirculation path 18, and an exhaust gas driven turbine 20, which includes a parallel waste-gate flow path 22. As known, the compressor is driven by the turbine, and the amount of boost is controlled principally by a waste-gate control mechanism (e.g., valve) shown schematically as a waste-gate valve 24. For purposes of the present invention, the waste-gate valve 24 may actuated through any of several waste-gate actuation mechanizations, including but not limited to (1) pneumatic actuation of the waste-gate valve using boost pressure as an actuation source, which will have a solenoid under the control of the controller 14 adjusting the amount of boost pressure routed to the actuator (see FIG. 2); (2) pneumatic actuation of the waste-gate valve using vacuum from a vacuum pump as the actuation source, which will have a solenoid under the control of the controller 14 adjusting the amount of vacuum routed to the actuator (see FIG. 3); and (3) an electrically controlled waste-gate valve (i.e., where the waste-gate valve is directly moved by an electric motor or the like under the control of controller 14—see FIG. 4).
On the air intake side of the engine 12, FIG. 1 shows an air intake port 26, an air filter 28, an intercooler 30 configured to cooperate with and complement compressor 16, a throttle valve 32, and an intake manifold 34. These features are well known and understood in the art. In the context of the present invention, these features may comprise conventional implementations.
On the exhaust side of the engine 12, FIG. 1 shows an exhaust gas manifold 36. Additionally, various downstream exhaust components are conventionally included in system 10, such as a catalytic converter and a muffler, and are shown schematically as a single exhaust restriction block 38, which feeds into exhaust gas outlet 40. These features are well known and understood in the art. In the context of the present invention, these features may comprise conventional implementations.
Conventionally, a variety of feedback paths are provided in system 10. For example, FIG. 1 shows an exhaust gas recirculation (EGR) tube or the like coupled between the exhaust manifold 36 and the intake manifold 34, and whose flow path is adjusted by way of an EGR valve 44. As known, the EGR valve 44 may be controlled by the electronic controller 14 in accordance with conventional EGR algorithms configured to achieve predetermined performance criteria. Generally, varying the position of the valve 44 alters the amount of exhaust gas that is provided to the intake manifold 34 for mixing with intake air, fuel and the like destined for combustion in engine 12.
With continued reference to FIG. 1, additional feeds may also be provided. For example, evaporative emissions control and diagnostics generally call for an evaporative (“evap”) emissions canister (not shown) be provided in an automotive vehicle that includes system 10. The evap canister is coupled to a fuel tank (not shown) as well as to inlets 46 and 48 by a combination of vent, purge and check valves, all as known in the art. For purposes of the present invention, however, the impact of these features may be ignored while the evap emissions and control system is not in operation (i.e., when not performing a purge cycle or a diagnostic cycle).
FIG. 1 also shows a variety of sensors deployed on the intake side of the engine 12, including an ambient or barometric pressure sensor 50 configured to produce a barometric pressure signal 52, an ambient air temperature sensor such as an intake air temperature (IAT) sensor 54 configured to generate an IAT signal 56, a boost air temperature sensor 58 configured to generate a boost air temperature signal 60, a boost pressure sensor 62 configured to generate a boost pressure signal 64, and an intake manifold pressure sensor such as a manifold absolute pressure (MAP) sensor 66 configured to generate a MAP signal 68. These sensors and their functioning are all well known and understood in the art. For purposes of the present invention, these sensors may all comprise conventional components.
Additionally, system 10 includes capabilities for determining a value for the mass air flow {dot over (m)}C, which may be obtained either via measurement by an air meter (e.g., mass air flow sensor or MAF sensor-not shown) typically placed just upstream of the compressor 16, or, in an alternate embodiment, calculated by the well known speed-density equation, for example as set forth in U.S. Pat. No. 6,393,903 entitled VOLUMETRIC EFFICIENCY COMPENSATION FOR DUAL INDEPENDENT CONTINUOUSLY VARIABLE CAM PHASING to Reed et al., assigned to the common assignee of the present invention, and incorporated herein by reference in its entirety.
As described in the Background, in systems that include a turbo-charger 15, as compared to a naturally aspirated engine, the engine management system (EMS) implemented in electronic controller 14 must perform all its normal calculations and control functions, but is also configured to control, among other things, the amount of boost provided by the turbo-charger.
FIG. 1 shows in block form various desired exhaust states 70, as described in co-pending application entitled “SYSTEM AND METHOD FOR MODELING TURBO-CHARGED ENGINES AND INDIRECT MEASUREMENT OF TURBINE AND WASTE-GATE FLOW AND TURBINE EFFICIENCY.” It should be understood that this representation is not meant to mean necessarily that these exhaust parameters are communicated physically outside of the controller 14.
As will be described in greater detail below, the present invention employs a first set of various data tables or the like to store required data for the purpose of exhaust state estimation, including but not limited to a compressor enthalpy delta data table 72 (sometimes referred to herein as “Table1”), a turbo-charger corrected rotational speed data table 74 (sometimes referred to herein as “Table2”), a turbine isentropic efficiency data table 76 (sometimes referred to herein as “Table3”), a turbine pressure ratio/corrected engine flow data table 78 (sometimes referred to herein as “Table4”), a turbine power and boundary term (A*B) data table 80 (sometimes referred to herein as “Table5”), an exhaust pressure drop data table 82 (sometimes referred to herein as “Table6”), and a heat transfer calibration data 84 (sometimes referred to herein as “Table7”).
FIG. 2 is a simplified flowchart of a method according to the invention. The method begins in step 86.
In step 86, a desired boost pressure is provided to the method for controlling boost pressure. The desired boost pressure is determined based on operator input and the like which collectively correspond to a desired produced engine torque. Typically, an engine management system (i.e., EMS implemented in controller 14) in which model-based boost control is implemented is of the torque based control type. In this case, a torque control structure will translate the driver torque request into the desired engine airflow, and the air control structure will translate that into a desired intake manifold pressure, and ultimately decide throttle position. An advantage of turbo-charged engines with actively controlled boost pressure is improved fuel economy by operating un-throttled as much as possible. This can be achieved by setting the desired boost pressure PBoostdes (upstream of the throttle) equal to the desired intake manifold pressure mentioned above. It should be understood, however, that trade-offs with regard to other operating parameter may result in a desired boost pressure that is different than what might be considered optimum for fuel efficiency purposes (e.g., the EMS may decide to keep the boost pressure higher than optimum and throttle it down using the throttle so as to keep the turbine rotating, which can improve responsiveness-many other situations are possible). Furthermore, since the maximum boost pressure that a turbo-charged engine is capable of producing will exceed the engines structural limitation, and knock limit for a spark ignition engine, the desired boost will be limited to a calibratable maximum value. The result of this conventional torque-based control logic is the input to the boost control logic of the present invention, namely; desired boost pressure PBoostdes (and optionally desired engine air flow {dot over (m)}eng,airdes which as noted above is calculated and available per a torque based control). Alternatively, however, in a non torque-based control configuration, where a desired engine flow is not normally calculated, the desired boost pressure can nonetheless be translated to a desired engine flow {dot over (m)}eng,air@des Boost by the well known speed-density calculation referred to elsewhere herein by replacing actual intake manifold pressure (e.g., MAP signal 68 in FIG. 1) with desired boost PBoostdes in the speed-density calculations. In either configuration, now having a given desired boost pressure PBoostdes, the method proceeds to step 88.
In step 88, the method involves determining a desired waste-gate mass flow rate at sonic standard conditions {dot over (m)}WGdes@SonicStdCond from the inputs desired Boost PBoostdes. The present invention uses a model of the exhaust states as the basis for model-based boost control. Ultimately, this allows the conversion of a desired boost into (1) the desired waste-gate valve flow and (2) desired waste-gate valve delta pressure (which equals desired turbine delta pressure), for all three waste-gate mechanizations described herein. It is significant that the exhaust state estimation model allows these calculations using only predictive (desired) variables as opposed to actual states, which avoids control instability. Furthermore, these calculations are identical for all waste-gate mechanizations to be described, and use various predetermined data tables (i.e., calibrations that may already be available for the estimation logic). The method proceeds to step 90.
In step 90, the method converts the desired waste-gate flow rate at standard sonic conditions {dot over (m)}WGdes@SonicStdCond into an actuator control command, and is unique for each mechanization type. In other words, a control signal is generated for controlling the waste-gate valve so as to achieve the desired waste-gate flow rate at standard sonic conditions {dot over (m)}WGdes@SonicStdCond. The further conversion of the desired waste-gate valve flow and delta pressure into the actuator command is specific to each mechanization. Generally, for the pneumatic actuators (i.e., boost pressure or vacuum type), it is recognized that the waste-gate valve position is the result of the force balance across the waste-gate valve linkage. The force on one side is generated by the delta pressure across the waste-gate valve and the counteracting force is generated by the delta pressure across the waste-gate actuation diaphragm. A given force balance results in a given waste-gate valve position, which is equivalent to an effective flow area which can be converted to a waste-gate valve flow. This means that a unique calibration can be created which describes the desired waste-gate diaphragm delta pressure as a function of desired waste-gate valve flow at standard sonic conditions and desired waste-gate delta pressure. This can be converted into the desired solenoid duty-cycle using the tabulated solenoid characteristic. For electrically-controlled type mechanization, the desired waste-gate valve flow can be converted into a desired waste-gate valve position by inverting the waste-gate flow characteristic just like converting desired throttle flow into desired throttle position.
FIG. 3 shows, in greater detail, a first pneumatic actuation mechanization 100 employing boost pressure. Turbine 20 includes a turbine flow path and a waste-gate flow path 22. The waste-gate valve 24 is connected via a linkage 102 to a waste-gate diaphragm 104. The waste-gate diaphragm 104, in turn, is pneumatically controlled via a waste-gate solenoid 106 coupled to a source of boost pressure. Typically, the source is the boost pressure generated by the turbo compressor.
FIG. 4 shows, in greater detail, a second pneumatic actuation mechanization 108 including a waste-gate diaphragm 110 and a waste-gate solenoid 112 coupled to a source of vacuum. Typically, the vacuum source is from an engine driven vacuum pump.
FIG. 5 shows, in greater detail, a third mechanization 114 where the waste-gate valve 24 is controlled, either directly or through use of a linkage 102, by an electrically-controlled device such as an electrical motor 116.
FIG. 6 is a block diagram view showing electronic controller 14 including additional tables specific to converting desired boost (see block 118) into specific control signals (see block 120) for all three waste-gate mechanizations 100, 108 and 114. Additional tables include a desired waste-gate diaphragm delta pressure data table 122 (sometimes referred to herein as “TableC1”), a transient response improvement term data table 124 (sometimes referred to herein as “TableC2”), a desired waste-gate solenoid (for boost pressure configurations) duty cycle data table 126 (sometimes referred to herein as “TableC3”), a desired waste-gate solenoid (for vacuum configurations) duty cycle data table 128 (sometimes referred to herein as “TableC4”) and a desired waste-gate valve position data table 130 (sometimes referred to herein as “TableC5”). The purpose of these tables will be described more fully below. Also, a description of how to populate the data in these tables will also be described below.
In the case of torque based control, the desired engine air flow calculated in the torque control structure of controller 14 may under high load request conditions exceed the actual engine airflow while operating at the maximum desired boost level. Therefore, the well known speed-density calculation is used to estimate the engine airflow at the maximum desired boost level {dot over (m)}eng,air@Max des Boost by replacing actual intake manifold pressure (e.g., MAP signal 68 in FIG. 1) with the maximum desired boost PBoostdes,Max in the speed-density calculations.
In the case of non-torque based control, where a desired engine flow is not normally calculated, the desired boost can be translated to a desired engine flow {dot over (m)}eng,air@des Boost by the well known speed-density calculation by replacing actual intake manifold pressure (e.g., MAP signal 68 in FIG. 1) with desired boost PBoostdes in the speed-density calculations.
The desired compressor flow {dot over (m)}Cdes can now be calculated as:
(1)




  
    
      
        
          
            m
            .
          
          C
          des
        
        =
        
          {
          
            
              
                
                  
                    m
                    .
                  
                  
                    eng
                    ,
                    air
                  
                  
                    
                      @
                      
                        
                      
                      ⁢
                      des
                    
                    ⁢
                    
                      
                    
                    ⁢
                    Boost
                  
                
              
              
                ,
              
              
                
                  non
                  ⁢
                  
                    -
                  
                  ⁢
                  torque
                  ⁢
                  
                    
                  
                  ⁢
                  based
                
              
            
            
              
                
                  MIN
                  ⁡
                  
                    (
                    
                      
                        
                          m
                          .
                        
                        
                          eng
                          ,
                          air
                        
                        
                          
                            
                          
                          ⁢
                          des
                        
                      
                      ,
                      
                        
                          m
                          .
                        
                        
                          eng
                          ,
                          air
                        
                        
                          
                            @
                            
                              
                            
                            ⁢
                            Max
                          
                          ⁢
                          
                            
                          
                          ⁢
                          des
                          ⁢
                          
                            
                          
                          ⁢
                          Boost
                        
                      
                    
                    )
                  
                
              
              
                ,
              
              
                
                  
                    torque
                    ⁢
                    
                      
                    
                    ⁢
                    based
                  
                  ⁢
                  
                    
                  
                
              
            
          
        
      
    
    
      
        
      
    
  




where {dot over (m)}Cdes is the desired compressor mass flow rate; {dot over (m)}eng,airdes is the desired engine air mass flow rate from the torque control structure; and {dot over (m)}eng,airdes Boost is the desired engine air mass flow rate at the desired boost pressure level; and {dot over (m)}eng,air@Max des Boost is the desired engine air mass flow rate at the maximum desired boost pressure level.

It is noted that the above calculation assumes that all engine airflow will flow through the compressor. If there are other air sources, e.g. IAC valve airflow or purge airflow, they should naturally be accounted for in equation (1).
The desired compressor pressure ratio is calculated as:
(2a) PC,outdes=PBoostdes+ΔPIntercooler where PC,outdes is the desired compressor outlet stagnation pressure; PBoostdes is the desired boost pressure (after intercooler, before throttle); and ΔPIntercooler is the intercooler pressure drop.
(2b) PC,in=Pamb−ΔPAirfilter where PC,in is the compressor inlet stagnation pressure; Pamb is the ambient pressure; and ΔPAirfilter is the air filter pressure drop.
(2c)




  
    
      
        
          Pr
          C
          des
        
        =
        
          
            P
            
              C
              ,
              out
            
            des
          
          
            P
            
              C
              ,
              
                i
                ⁢
                
                  
                
                ⁢
                n
              
            
          
        
      
    
    
      
        
      
    
  




where PrCdes is the desired compressor pressure ratio; PC,outdes is the desired compressor outlet stagnation pressure; and PC,in is the compressor inlet stagnation pressure.

The desired compressor power PCdes is calculated as:
(3a)




  
    
      
        
          
            V
            .
          
          
            C
            ,
            cor
          
          des
        
        =
        
          
            
              m
              .
            
            C
            des
          
          *
          
            
              
                T
                
                  C
                  ,
                  
                    i
                    ⁢
                    
                      
                    
                    ⁢
                    n
                  
                
              
              ⁢
              R
            
            
              P
              
                C
                ,
                
                  i
                  ⁢
                  
                    
                  
                  ⁢
                  n
                
              
            
          
          *
          
            
              
                T
                
                  C
                  ,
                  
                    i
                    ⁢
                    
                      
                    
                    ⁢
                    n
                  
                
                reference
              
              
                T
                
                  C
                  ,
                  
                    i
                    ⁢
                    
                      
                    
                    ⁢
                    n
                  
                
              
            
          
        
      
    
    
      
        
      
    
  




where {dot over (V)}C,cordes is the desired corrected compressor volume flow rate; {dot over (m)}Cdes is the desired compressor mass flow rate; TC,in is the compressor inlet air stagnation temperature; R is the gas constant; PC,in is the compressor inlet stagnation pressure; and TC,inreference is the reference compressor inlet air stagnation temperature (from turbo manufacturers data).

(3b) {dot over (N)}T,cordes=Table2({dot over (V)}C,cordes, PrCdes) where {dot over (N)}T,cordes is the desired corrected turbo-charger rotational speed; {dot over (V)}C,cordes is the desired corrected compressor volume flow rate; and PrCdes is the desired compressor pressure ratio.
(3c) ΔhCdes=Table1({dot over (V)}C,cordes, {dot over (N)}T,cordes) where ΔhCdes is the desired compressor enthalpy delta; {dot over (V)}C,cordes is the desired corrected compressor volume flow rate; and {dot over (N)}T,cordes is the desired corrected turbo-charger rotational speed.
(3d) PCdes={dot over (m)}Cdes*ΔhCdes where PCdes is the desired power absorbed by compressor; {dot over (m)}Cdes is the desired compressor mass flow rate; and ΔhCdes is the desired compressor enthalpy delta.
During real-time execution by the electronic controller 14, the dependencies in the above equations (3a) through (3d) are evaluated.
First, equation (3a) is evaluated to obtain a value for the desired corrected volume flow rate {dot over (V)}C,cordes. The dependencies in the right hand side of equation (3a) can be met either through direct sensor measurement or via estimation. For example, TC,inreference and R will be known, TC,in≈Tamb will be known via measurement by intake air temperature (IAT) sensor 54, PC,in will be known via measurement by ambient pressure sensor 50 (Pamb) as modified by ΔPAirfilter (i.e., using equation (2b)), and the desired compressor flow {dot over (m)}Cdes will be known via prior evaluation of equation (1). It should be understood that the pressure drops across the air filter 28 and intercooler 30, ΔPAirfilter and ΔPIntercooler respectively, may be empirically determined by characterizing such pressure drops as a function of engine flow.
Next, equation (3b) is evaluated to obtain a value for the desired corrected turbo rotational speed ({dot over (N)}T,cordes). The inputs needed for use of the corrected rotational speed data table 74 (“Table2”) are the desired corrected volume flow rate {dot over (V)}C,cordes, which can be calculated from equation (3a) and the desired compressor pressure ratio PrCdes, which can be calculated from equation (2c). The compressor inlet pressure PC,in in has already been calculated in evaluating equation (2b). The desired compressor outlet pressure PC,outdes may be determined using equation (2a). The desired compressor pressure ratio PrCdes is then calculated, thus allowing equation (3b) to be evaluated to obtain a value for the desired, corrected rotational speed {dot over (N)}T,cordes.
Next, equation (3c) is evaluated to obtain a value for the desired compressor enthalpy delta ΔhCdes, using the just-determined values for {dot over (V)}C,cordes and {dot over (N)}T,cordes.
Finally, equation (3d) is evaluated, using the value of the desired compressor enthalpy delta ΔhCdes and the desired mass flow rate {dot over (m)}Cdes, to obtain a value of the desired power absorbed by the compressor PCdes.
The desired turbine flow {dot over (m)}Tdes is calculated as:
(4) ηT=Table3(PrT) where ηT is the turbine isentropic efficiency and PrT is the turbine pressure ratio.
(5)




  
    
      
        
          
            A
            des
          
          =
          
            (
            
              
                P
                C
                des
              
              
                
                  η
                  T
                
                ⁢
                
                  c
                  p
                  e
                
                ⁢
                
                  T
                  
                    T
                    ,
                    
                      i
                      ⁢
                      
                        
                      
                      ⁢
                      n
                    
                  
                
              
            
            )
          
        
        ,
        
          
            B
            des
          
          =
          
            B
            =
            
              
                
                  T
                  
                    T
                    ,
                    
                      i
                      ⁢
                      
                        
                      
                      ⁢
                      n
                    
                  
                
              
              
                P
                
                  T
                  ,
                  out
                
              
            
          
        
      
    
    
      
        
      
    
  




where Ades is the desired turbo power term and





  A
  =
  
    (
    
      
        P
        C
      
      
        
          η
          T
        
        ⁢
        
          c
          p
          e
        
        ⁢
        
          T
          
            T
            ,
            
              i
              ⁢
              
                
              
              ⁢
              n
            
          
        
      
    
    )
  




is the turbo power term; Bdes is the desired turbine boundary term; and





  B
  =
  
    
      
        T
        
          T
          ,
          
            i
            ⁢
            
              
            
            ⁢
            n
          
        
      
    
    
      P
      
        T
        ,
        out
      
    
  




is the turbine boundary term.

(6) {dot over (m)}T,cordes=Table5(Ades*Bdes) where {dot over (m)}T,cordes is the desired corrected turbine mass flow rate; Ades is the desired turbo power term; and Bdes is the desired turbine boundary term.
(7)




  
    
      
        
          
            m
            .
          
          T
          des
        
        =
        
          
            
              m
              .
            
            
              T
              ,
              cor
            
            des
          
          ⁢
          
            
              P
              
                T
                ,
                out
              
            
            
              
                T
                
                  T
                  ,
                  
                    i
                    ⁢
                    
                      
                    
                    ⁢
                    n
                  
                
              
            
          
        
      
    
    
      
        
      
    
  




where {dot over (m)}Tdes is the desired turbine mass flow rate; {dot over (m)}T,cordes is the desired corrected turbine mass flow rate; PT,out is the turbine outlet stagnation pressure; and TT,in is the turbine inlet stagnation temperature.

The desired turbine pressure ratio PrTdes and turbine inlet PEMdes and outlet PT,outdes pressures are calculated as:
(8) PrTdes=Table4({dot over (m)}T,cordes) where PrTdes is the desired turbine pressure ratio and {dot over (m)}T,cordes is the desired corrected turbine mass flow rate.
(9a)




  
    
      
        
          tmp
          des
        
        =
        
          
            
              
                m
                .
              
              exh
              des
            
            ⁢
            
              
                R
                *
                
                  T
                  
                    T
                    ,
                    out
                  
                
              
            
          
          
            
              P
              amb
            
            *
            Calibration
            ⁢
            
              
            
            ⁢
            1
          
        
      
    
    
      
        
      
    
  




where tmpdes is a temporary variable used in equation (9a) to allow easier substitution into equation 9(b); {dot over (m)}exhdes is the desired exhaust system mass flow rate; R is the gas constant; TT,out is the turbine outlet stagnation temperature; and Pamb is the ambient pressure.

(9b) PT,outdes=Pamb*Table4(tmpdes) where PT,outdes is the desired turbine outlet stagnation pressure; Pamb is the ambient pressure; and tmpdes is the temporary variable referred to above in connection with equation (9a).
(10) PEMdes=PT,indes=PrTdes*PT,outdes where PEMdes is the desired turbine outlet stagnation pressure; PT,indes is the desired turbine inlet stagnation pressure; PrTdes is the desired turbine pressure ratio; and PT,outdes is the desired turbine outlet stagnation pressure.
The desired waste-gate mass flow {dot over (m)}WGdes and desired waste-gate mass flow at sonic standard conditions {dot over (m)}WGdes@SonicStdCond are calculated as:
(11) {dot over (m)}WGdes={dot over (m)}Cdes−{dot over (m)}Tdes where {dot over (m)}WGdes is the desired waste-gate mass flow rate; {dot over (m)}Cdes is the desired compressor mass flow rate; and {dot over (m)}Tdes is the desired turbine mass flow rate.
(12a)




  
    
      
        
          
            m
            .
          
          WG
          
            des
            ⁢
            
              
            
            @
            
              
            
            ⁢
            SonicStdCond
          
        
        =
        
          
            
              m
              .
            
            WG
            des
          
          ⁢
          
            
              P
              std
            
            
              P
              EM
              des
            
          
          ⁢
          
            
              
                T
                EM
              
              
                T
                std
              
            
          
          ⁢
          
            1
            
              
                β
                2
              
              ⁡
              
                (
                
                  1
                  /
                  
                    Pr
                    T
                    des
                  
                
                )
              
            
          
        
      
    
    
      
        
      
    
  




where {dot over (m)}WGdes@ SonicStdCond is desired waste-gate mass flow rate at standard sonic conditions; {dot over (m)}WGdes is the desired waste-gate mass flow rate; Pstd is a Standard Pressure (a defined reference); Tstd is a Standard Temperature (a defined reference); PEMdes is the desired exhaust manifold stagnation pressure; TEM is the exhaust manifold gas stagnation temperature; β2 is the normalized valve flow dependency on pressure ratio for compressible fluid; and PrTdes is the desired turbine pressure ratio.

(12b)




  
    
      
        
          tmp
          ⁡
          
            (
            Pr
            )
          
        
        =
        
          
            
              2
              
                k
                -
                1
              
            
          
          ⁢
          
            
              
                
                  (
                  Pr
                  )
                
                
                  2
                  /
                  k
                
              
              -
              
                
                  (
                  Pr
                  )
                
                
                  
                    (
                    
                      k
                      +
                      1
                    
                    )
                  
                  /
                  k
                
              
            
          
        
      
    
    
      
        
      
    
  




where tmp(Pr) is a temporary variable used in equation (12b) to allow for easier substitution in equation (12d);

(12c)




  
    
      
        
          Pr
          crit
        
        =
        
          
            (
            
              2
              
                k
                +
                1
              
            
            )
          
          
            k
            /
            
              (
              
                k
                -
                1
              
              )
            
          
        
      
    
    
      
        
      
    
  




where Prcrit is the critical pressure ratio;

(12d)




  
    
      
        
          β
          2
        
        =
        
          {
          
            
              
                1
              
              
                ,
              
              
                
                  
                    for
                    ⁢
                    
                      
                    
                    ⁢
                    Pr
                  
                  <=
                  
                    Pr
                    crit
                  
                
              
            
            
              
                
                  
                    tmp
                    ⁡
                    
                      (
                      Pr
                      )
                    
                  
                  
                    tmp
                    ⁡
                    
                      (
                      
                        Pr
                        crit
                      
                      )
                    
                  
                
              
              
                ,
              
              
                
                  
                    for
                    ⁢
                    
                      
                    
                    ⁢
                    Pr
                  
                  >
                  
                    Pr
                    crit
                  
                
              
            
          
        
      
    
    
      
        
      
    
  




where β2 is the normalized valve flow dependency on pressure ratio for compressible fluid and Prcrit is critical pressure ratio;

and where Beta2 (β2) is the normalized effect of pressure ratio (Pr) on valve flow of a compressible fluid (gas), and is detailed in equations (12b)-(12d) above. Its deduction is e.g., given in Appendix C, John B. Heywood, “Internal Combustion Engine Fundamentals,” ISBN P-070100499-8, herein incorporated by reference in its entirety. For its use for waste-gate flow, the ratio of specific heats for exhaust gas should be used, k=kexh=1.33. For real time implementations, the Beta2 (β2) parameter should be tabulated in advance (i.e., as opposed to real time calculation).
For the pneumatic waste-gate mechanizations 100 and 108, the desired delta pressure across the waste-gate is the same as the desired delta pressure across the turbine:
(13) ΔPWGdes=ΔPTdes=PEMdes−PT,outdes where ΔPWGdes is the desired waste-gate valve delta pressure; ΔPTdes is the desired turbine delta pressure; PEMdes is the desired exhaust manifold stagnation pressure; and PT,outdes is the desired turbine outlet stagnation pressure.
FIG. 3 shows the physical linkage between the waste-gate valve delta pressure ΔPWG and the waste-gate diaphragm delta pressure ΔPWG,Dphr. The position of the waste-gate valve 24 is the result of the force balance across the waste-gate valve linkage. The force on one side is generated by the delta pressure across the waste-gate valve ΔPWG. The counteracting force is generated by the delta pressure across the waste-gate actuation diaphragm ΔPWG,Dphr. A given force balance results in a given position of the waste-gate valve 24, which position is equivalent to an effective flow area. The effective flow area in turn can be converted to a waste-gate valve flow at sonic standard conditions. This means that a unique calibration can be created which describes the desired waste-gate diaphragm delta pressure ΔPWG,Dphrdes as a function of desired waste-gate valve flow at sonic standard conditions {dot over (m)}WGdes@SonicStdCond and desired waste-gate delta pressure ΔPWGdes.
(14) ΔPWG,Dphrdes,base=TableC1({dot over (m)}WGdes@SonicStdCond, ΔPWG) where ΔPWG,Dphrdes,base is the desired waste-gate diaphragm delta pressure, base look-up; {dot over (m)}WGdes@SonicStdCond is the desired waste-gate mass flow rate at standard sonic conditions; and ΔPWGdes is the desired waste-gate valve delta pressure.
Furthermore, boost control transient response may be improved further by adding a transient term calibrated as a function of boost control error:
(15) ΔPWG,Dphrdes,TransientTerm=TableC2(PBoost−PBoost) where ΔPWG,Dphrdes,TransientTerm is the desired waste-gate diaphragm delta pressure, transient term; PBoost is the boost pressure (after intercooler, before throttle); and PBoostdes is the desired Boost Pressure (after intercooler, before throttle).
Where transient response improvement term table 124 (“TableC2”) should be calibrated to zero for zero boost control error (PBoost−PBoostdes=0).
When adding a closed-loop correction term ΔPWG,Dphrdes,CL for robustness to calibration and model inaccuracies, it should be applied to ΔPWG,Dphrdes.
(16) ΔPWG,Dphrdes=ΔPWG,Dphrdes,CL+ΔPWG,Dphrdes,CL+ΔPWG,Dphrdes,TransientTerm where ΔPWG,Dphrdes is the desired waste-gate diaphragm delta pressure; ΔPWG,Dphrdes,base is the desired waste-gate diaphragm delta pressure, base look-up; ΔPWG,Dphrdes,CL is the desired waste-gate diaphragm delta pressure, closed-loop term; and ΔPWG,Dphrdes,TransientTerm is the desired waste-gate diaphragm delta pressure, transient term.
As to the waste-gate mechanization 100 (FIG. 3), it can be seen that the desired waste-gate diaphragm delta pressure ΔPWG,Dphrdes can be converted into the desired waste-gate solenoid delta pressure ΔPWG,S In ddes by:
(17) ΔPWG,S In ddes=PBoostdes−Pamb−ΔPWG,Dphrdes where ΔPWG,S In ddes is the desired waste-gate solenoid delta pressure; PBoostdes is the desired Boost Pressure (after intercooler, before throttle); Pamb is the ambient pressure; and ΔPWG,Dphrdes is the desired waste-gate diaphragm delta pressure.
The relationship between the solenoid duty cycle, its upstream pressure and delta pressure can be tabulated, and thus used to look-up the waste-gate solenoid duty cycle DutyCycleWG,S In d:
(18) DutyCycleWG,S In d=TableC3(ΔPWG,S In ddes, PBoost) where DutyCycleWG,S In d is the waste-gate solenoid commanded duty-cycle; ΔPWG,S In ddes is the desired waste-gate solenoid delta pressure; and PBoost is the boost pressure (after intercooler, before throttle).
As to the waste-gate mechanization 108 (FIG. 4), when using a vacuum source (pump) to actuate the waste-gate diaphragm, equation (19a) defines the relationship:

ΔPWG,S In ddes=Pvac−Pamb−ΔPWG,Dphrdes  (19a)

where ΔPWG,S In ddes is the desired waste-gate solenoid delta pressure; Pvac is the pressure of the vacuum source; Pvac is the ambient pressure; and ΔPWG,Dphrdes is the desired waste-gate diaphragm delta pressure.
The relationship between the solenoid duty cycle, its upstream pressure and delta pressure can be tabulated, and thus used to look-up the waste-gate solenoid duty cycle DutyCycleWG,S In d as set forth in equation (19b).

DutyCycleWG,S In d=TableC4(ΔPWG,S In ddes,Pvac)  (19b)

where DutyCycleWG,S In d is the waste-gate solenoid commanded duty-cycle; ΔPWG,Dphrdes is the desired waste-gate diaphragm delta pressure; and Pvac is the pressure of the vacuum source.
For a constant vacuum, equation (19b) may be simplified into equation (19c):
(19c) DutyCycleWG,S In d=TableC4(ΔPWG,Dphrdes) where DutyCycleWG,S In d is the waste-gate solenoid commanded duty-cycle and ΔPWG,Dphrdes is the desired waste-gate diaphragm delta pressure.
As to the waste-gate mechanization 114 (FIG. 5) based on an electric motor-controlled waste-gate valve, a tabulated waste-gate valve flow characteristic is used to convert the desired waste-gate valve flow at sonic standard conditions {dot over (m)}WGdes@SonicStdCond into a desired valve position:
(20) ValvePosWGdes=TableC5({dot over (m)}WGdes@SonicStdCond) where ValvePosWGdes is the desired waste-gate valve position and {dot over (m)}WGdes@SonicStdCond is the desired waste-gate mass flow rate at standard sonic conditions.
Finally, it is noted that the model based controls approach also works for engine with both waste-gate and VGT by simply adding the VGT position dependency to tables 76, 78 and 80 (i.e., “Table3”, “Table4” and “Table5”).
It should be appreciated that a model form as described above allows for an explicit one-step through calculation of the desired waste-gate flow and delta pressure from the desired boost pressure. This one-step through (“loop”) solvability makes it suitable for real-time systems, as opposed to iterative approaches which are too resource intensive to be implemented in practical, real time systems. Additionally, this model form ensures good transient boost response, because it is based on predictive states (as opposed to actual).
A description detailing how to populate estimation model data tables 74-84 will now be set forth.
It is common practice for manufacturers of turbo-chargers to make the following data available: {dot over (V)}C,cor, PrC, {dot over (N)}T,cor, ηC at a specified TC,inreference, PC,inreference, where {dot over (V)}C,cor is the corrected compressor volume flow rate, {dot over (N)}T,cor is the corrected turbo-charger rotational speed and TC,inreference, PC,inreference are the reference compressor inlet air stagnation temperature and reference compressor inlet air stagnation pressure, respectively. This data is what is referred to herein as the compressor characteristics map. While one of ordinary skill in the art will understand the form of the data commonly available, for the sake of clarity, and for example purposes only, a turbo-charger manufacturer may provide a map/diagram that graphically illustrates the compressor isentropic efficiency ηC (e.g., as rings or partial rings indicating various efficiency levels such as 50%, 60%, 70%, 72%, 75%, etc.) on an X-Y chart where the X-axis is the corrected volume (air) flow rate and the Y-axis is the compressor pressure ratio. Overlaying the efficiency level rings may be a series of traces, generally offset from each other, with each indicating a respective turbo-charger corrected speed. It should be understood, however, that other forms of data or information may be available, including without limitation data in electronic format. Moreover, while it is preferred, due to convenience primarily, that such manufacturer-provided data concerning the compressor be used, it should be understood that independent measurements and characterization of the compressor of a turbo-charger may be employed to obtain the same information as described above, without any limitation of the present invention.
The values in tables 72 and 74 (“Table1” and “Table2” in the equations) are preferably calculated and stored in advance of real-time execution by electronic controller 14. Specifically, both tables 72 and 74 may be populated in advance using the compressor characteristic map.
The compressor enthalpy delta (ΔhC) table 72 (“Table1”) is a table that takes as inputs the corrected compressor volume flow rate ({dot over (V)}C,cor) and the corrected turbo-charger rotational speed ({dot over (N)}T,cor) and provides as an output a value for the compressor enthalpy delta (i.e., change in enthalpy of the gas across the compressor), in the form of equation (T1) below:

ΔhC=Table1({dot over (V)}C,cor,{dot over (N)}T,cor)  (T1)

Table 72 may be populated using the data conventionally available from the turbo manufacturers characterizing the compressor, as well as the following equation (T2):




  
    
      
        
          Δ
          ⁢
          
            
          
          ⁢
          
            h
            C
          
        
        =
        
          
            1
            
              η
              C
            
          
          ⁢
          
            c
            
              p
              ,
              air
            
          
          ⁢
          
            
              T
              
                C
                ,
                
                  i
                  ⁢
                  
                    
                  
                  ⁢
                  n
                
              
            
            ⁡
            
              [
              
                
                  
                    (
                    
                      Pr
                      C
                    
                    )
                  
                  
                    (
                    
                      
                        (
                        
                          
                            k
                            air
                          
                          -
                          1
                        
                        )
                      
                      /
                      
                        k
                        air
                      
                    
                    )
                  
                
                -
                1
              
              ]
            
          
        
      
    
    
      
        (
        T2
        )
      
    
  




The corrected turbo-charger rotational speed ({dot over (N)}T,cor) table 74 is a table that takes as inputs the corrected compressor volume flow rate ({dot over (V)}C,cor) and the compressor pressure ratio (PrC), and provides as an output the corrected turbo-charger rotational speed ({dot over (N)}T,cor), as in the form of equation (T3) below:

{dot over (N)}T,cor=Table2({dot over (V)}C,cor,PrC)  (T3)

Table 74 may be populated by using the data conventionally available from the turbo manufacturer characterizing the compressor.
FIG. 7 is a chart showing turbine isentropic efficiency as a function of turbine pressure ratio, as observed in several turbo-charged engines. FIG. 7 shows that a simple model of the turbine isentropic efficiency ηT is adequate for the present invention. Table 76(“Table3”) is preferably calculated and stored in advance of real-time execution by electronic controller 14, and takes the form as set forth below in equation (T4):

ηT=Table3(PrT)  (T4)

The data for table 76 (“Table3”) is preferably populated as follows. In one embodiment, the data for table 76 (“Table3”) may be obtained by the indirect measurement of turbine isentropic efficiency ηT from equation (T5) below, by obtaining 1) measurements of TC,in, TT,in, PC,in, PC,out, PT,in, PT,out, {dot over (m)}C; 2) {dot over (m)}T from the indirect method described below in the section “Turbine and waste-gate flow indirect measurement methods”; and 3) ηC from reference to the compressor characteristics data provided by the turbo manufacturer. Additionally, the turbine pressure ratio PrT may be calculated from measurements of PT,in, PT,out. Trace 132 is exemplary of the data that may populate the data table.




  
    
      
        
          η
          T
        
        =
        
          
            1
            
              η
              C
            
          
          ⁢
          
            
              
                m
                .
              
              C
            
            
              
                m
                .
              
              T
            
          
          ⁢
          
            
              c
              
                p
                ,
                air
              
            
            
              c
              
                p
                ,
                exh
              
            
          
          ⁢
          
            
              T
              
                C
                ,
                
                  i
                  ⁢
                  
                    
                  
                  ⁢
                  n
                
              
            
            
              T
              
                T
                ,
                
                  i
                  ⁢
                  
                    
                  
                  ⁢
                  n
                
              
            
          
          ⁢
          
            
              [
              
                
                  
                    (
                    
                      Pr
                      C
                    
                    )
                  
                  
                    (
                    
                      
                        (
                        
                          
                            k
                            air
                          
                          -
                          1
                        
                        )
                      
                      /
                      
                        k
                        air
                      
                    
                    )
                  
                
                -
                1
              
              ]
            
            
              [
              
                1
                -
                
                  
                    (
                    
                      Pr
                      T
                    
                    )
                  
                  
                    (
                    
                      
                        (
                        
                          
                            k
                            exh
                          
                          -
                          1
                        
                        )
                      
                      /
                      
                        k
                        exh
                      
                    
                    )
                  
                
              
              ]
            
          
        
      
    
    
      
        (
        T5
        )
      
    
  




FIG. 8 is a chart showing turbine pressure ratio PrT as a function of corrected engine flow. The following empirical relationship (i.e., equations (T6) and (T7)) for a turbine shown in FIG. 8 was found in Müller, M., et al., “Mean Value Modeling of Turbocharged Spark Ignition Engines”, Society of Automotive Engineers (SAE) Technical Paper, no. 980784, International Congress and Exposition, Detroit, Mich., USA, Feb. 23-26, 1998, incorporated by reference herein in its entirety.




  
    
      
        
          
            Pr
            T
          
          =
          
            Table
            ⁢
            
              
            
            ⁢
            4
            ⁢
            
              (
              
                
                  m
                  .
                
                
                  T
                  ,
                  cor
                
              
              )
            
          
        
        ⁢
        
          

        
        ⁢
        Where
      
    
    
      
        (
        T6
        )
      
    
  
  
    
      
        
          
            m
            .
          
          
            T
            ,
            cor
          
        
        =
        
          
            
              m
              .
            
            T
          
          ⁢
          
            
              
                T
                
                  T
                  ,
                  
                    i
                    ⁢
                    
                      
                    
                    ⁢
                    n
                  
                
              
            
            
              P
              
                T
                ,
                out
              
            
          
        
      
    
    
      
        (
        T7
        )
      
    
  




The pressure ratio/corrected engine flow data table 78 (i.e., the “Table4” in equation (T6)) may be populated with values in advance of real-time execution by electronic controller 14 using measurements of PT,in, PT,out. TT,in where {dot over (m)}T={dot over (m)}eng is used for all data points where the waste-gate and EGR valve are closed. More specifically, first measuring PT,in, PT,out, TT,in and {dot over (m)}eng test data. Then generating a first set of data-pairs for turbine pressure ratio (PrT) and corrected engine mass flow rate ({dot over (m)}eng,cor) using the measured test data. This would involve converting measured {dot over (m)}eng to {dot over (m)}eng,cor using the form of equation (T7) but substituting engine flow rate for turbine flow rate. Then, generating a second set of data-pairs for turbine pressure ratio (PrT) and corrected turbine mass flow rate ({dot over (m)}T,cor) from the first set of data pairs wherein the second set of data-pairs is indicative of conditions where the waste-gate valve would be closed. And finally, populating the data table 78 (“Table4”) with the second set of data-pairs. Such conditions would include the values for the turbine pressure ratio (PrT), for a given corrected turbine mass flow rate ({dot over (m)}T,cor) would always be equal to or greater than the turbine pressure ratio (PrT) for the same value {dot over (m)}eng,cor (i.e., from the first set of data pairs). One would understand that this data in table 78 represent a smooth and continuous shape (e.g., as shown as trace 134 in FIG. 8).
Equation (T8) is set forth below:




  
    
      
        
          
            (
            
              
                P
                C
              
              
                
                  η
                  T
                
                ⁢
                
                  c
                  p
                  e
                
                ⁢
                
                  T
                  
                    T
                    ,
                    
                      i
                      ⁢
                      
                        
                      
                      ⁢
                      n
                    
                  
                
              
            
            )
          
          *
          
            (
            
              
                
                  T
                  
                    T
                    ,
                    
                      i
                      ⁢
                      
                        
                      
                      ⁢
                      n
                    
                  
                
              
              
                P
                
                  T
                  ,
                  out
                
              
            
            )
          
        
        =
        
          
            
              m
              .
            
            
              T
              ,
              cor
            
          
          ⁡
          
            [
            
              1
              -
              
                
                  (
                  
                    Pr
                    T
                  
                  )
                
                
                  (
                  
                    
                      (
                      
                        
                          k
                          exh
                        
                        -
                        1
                      
                      )
                    
                    /
                    
                      k
                      exh
                    
                  
                  )
                
              
            
            ]
          
        
      
    
    
      
        (
        T8
        )
      
    
  




Two new variables A, B are introduced, as set forth in equation (T9) below.




  
    
      
        
          A
          =
          
            (
            
              
                P
                C
              
              
                
                  η
                  T
                
                ⁢
                
                  c
                  p
                  e
                
                ⁢
                
                  T
                  
                    T
                    ,
                    
                      i
                      ⁢
                      
                        
                      
                      ⁢
                      n
                    
                  
                
              
            
            )
          
        
        ,
        
          
        
        ⁢
        
          B
          =
          
            
              
                T
                
                  T
                  ,
                  
                    i
                    ⁢
                    
                      
                    
                    ⁢
                    n
                  
                
              
            
            
              P
              
                T
                ,
                out
              
            
          
        
      
    
    
      
        (
        T9
        )
      
    
  




Where A is a Turbo Power Term; and
B is a Turbine Boundary Term.
Substituting new variables A and B of equation (T9) into equation (T8) results in a rewritten equation, designated as equation (T10) below:

A*B={dot over (m)}T,cor[1−(PrT)((kexh−1)/kexh)]  (T10)

Equation (T6), making use of data table 78 (“Table4”), shows that there is a unique relationship between {dot over (m)}T,cor and PrT, and equation (T10) therefore stands for the proposition that a given product (A*B) describes one such unique combination of {dot over (m)}T,cor and PrT. It is therefore possible to directly tabulate the relationship between (A*B) and {dot over (m)}T,cor in a further table, namely turbine power and boundary term data table 80 (“Table5”), the form of which is set forth in equation (T11) below.

{dot over (m)}T,cor=Table5(A*B).  (T11)

Table 80 (“Table5”) may be populated with values in advance of real-time execution by controller 14 (i.e., can be calculated off-line) using data table 78 (“Table4”), by calculating PrT from equation (T6) using an array of {dot over (m)}T,cor and then calculating the corresponding (A*B) values from equation (T10).
Turbine and Waste-Gate Flow Indirect Measurement Methods
FIG. 8 and equation (T6) describe the relationship between PrT and {dot over (m)}T,cor, which is calibrated for the data points when the waste-gate valve (e.g., waste-gate valve 24 in FIG. 1) is closed. The data points for an open waste-gate valve, however, will deviate from this calibration. Inverting equation (T6), gives equation (T12) below.

{dot over (m)}Tcor=Table4−1(PrT)  (T12)

Equation (T7) may be used to convert {dot over (m)}T,cor to {dot over (m)}T. A value for {dot over (m)}T may then be used to obtain the waste-gate valve flow rate {dot over (m)}WG={dot over (m)}exh−{dot over (m)}T where {dot over (m)}exh={dot over (m)}eng−{dot over (m)}EGR.
The pressure drop data table 82 (“Table6”), contains predetermined data characterizing the pressure drop characteristics of catalyst/muffler restriction 38 as a function of flow. PT,out is assumed known as this may be described by conventional models of the pressure drop across the catalyst and muffler, which opens to atmosphere (i.e., the ambient pressure at tailpipe exit opening 40—best shown in FIG. 1). For completeness, equations (T13a) and (T13b) below may be used by electronic controller 14 to produce a value for the turbine outlet pressure PT,out.




  
    
      
        tmp
        =
        
          
            
              
                m
                .
              
              exh
            
            ⁢
            
              
                R
                *
                
                  T
                  
                    T
                    ,
                    out
                  
                
              
            
          
          
            
              P
              amb
            
            *
            Calibration
            ⁢
            
              
            
            ⁢
            1
          
        
      
    
    
      
        (
        T13a
        )
      
    
  
  
    
      
        
          P
          
            T
            ,
            out
          
        
        =
        
          
            P
            amb
          
          *
          Table
          ⁢
          
            
          
          ⁢
          6
          ⁢
          
            
          
          ⁢
          
            (
            tmp
            )
          
        
      
    
    
      
        (
        T13b
        )
      
    
  




Where tmp is a temporary variable, and
Calibration1 is a scaling calibration value selected so that the variable “tmp” is within the bounds or limits established for Table6.
Part of the exhaust gas flows through the turbine across which the temperature drops as the expansion process is converted into absorbed power of the turbine. The other part of the exhaust gas which flows through the waste-gate path does not experience a temperature drop due to expansion-work. Furthermore, due to the high flows and turbulence levels in the turbine and waste-gate path, a temperature drop is caused by convective heat transfer to the turbine and waste-gate housing. Therefore the model for the turbine outlet temperature is a mixing model of the outlet temperatures from the turbine and waste-gate plus the heat transfer model, expressed in equations (T14a), (T14b) and (T14c) below:




  
    
      
        tmp
        =
        
          1
          -
          
            
              η
              T
            
            ⁡
            
              (
              
                1
                -
                
                  
                    (
                    
                      
                        P
                        
                          T
                          ,
                          out
                        
                      
                      
                        P
                        
                          T
                          ,
                          
                            i
                            ⁢
                            
                              
                            
                            ⁢
                            n
                          
                        
                      
                    
                    )
                  
                  
                    (
                    
                      
                        (
                        
                          
                            k
                            exh
                          
                          -
                          1
                        
                        )
                      
                      /
                      
                        k
                        exh
                      
                    
                    )
                  
                
              
              )
            
          
        
      
    
    
      
        (
        T14a
        )
      
    
  
  
    
      
        
          T
          
            T
            ,
            out
            ,
            
              
                w
                /
                o
              
              ⁢
              
                
              
              ⁢
              HT
            
          
        
        =
        
          
            T
            
              T
              ,
              
                i
                ⁢
                
                  
                
                ⁢
                n
              
            
          
          *
          
            (
            
              
                tmp
                ⁢
                
                  
                
                ⁢
                
                  
                    
                      m
                      .
                    
                    T
                  
                  
                    
                      m
                      .
                    
                    exh
                  
                
              
              +
              
                
                  
                    m
                    .
                  
                  WG
                
                
                  
                    m
                    .
                  
                  exh
                
              
            
            )
          
        
      
    
    
      
        (
        T14b
        )
      
    
  
  
    
      
        
          T
          
            T
            ,
            out
          
        
        =
        
          
            T
            amb
          
          -
          
            
              (
              
                
                  T
                  amb
                
                -
                
                  T
                  
                    T
                    ,
                    out
                    ,
                    
                      
                        w
                        /
                        o
                      
                      ⁢
                      
                        
                      
                      ⁢
                      HT
                    
                  
                
              
              )
            
            ⁢
            Table
            ⁢
            
              
            
            ⁢
            7
            ⁢
            
              
            
            ⁢
            
              (
              
                
                  m
                  .
                
                exh
              
              )
            
          
        
      
    
    
      
        (
        T14c
        )
      
    
  




Where tmp is a temporary variable, and
“Table7” corresponds to heat transfer calibration data table 84.
Heat transfer calibration data table 84 is preferably populated with data in advance of real time execution by controller 14. Table 84 (“Table7”) may be populated by obtaining; 1) measurements of: Tamb, TT,in, TT,out, PT,in, PT,out, {dot over (m)}exh; 2) {dot over (m)}T, {dot over (m)}WG from the methods described in the section “Turbine and waste-gate flow indirect measurement methods”; 3) ηT either from method described in the section “Turbine isentropic efficiency indirect measurement method” or set as a reasonable constant value, e.g. 55%. Rearranging equations (T14a)-(T14c) in order to isolate “Table7”, the value of Table7 can be calculated for each data point based on the data gathered, as described above in this paragraph in 1), 2) and 3). Table 84 (“Table7”) can now be tabulated with that data. In sum, the heat transfer data in table 84 (“Table7”) may be calibrated by fitting the model to the measured TT,out or a reasonable, constant value of ηT may be used.
A description detailing how to populate data tables 122-130 (“TableC1” through “TableC5”) will now be set forth.
FIG. 9 shows graphically an example of the data table 122 (37 TableC1”). The desired waste-gate diaphragm delta pressure data table 122 (“TableC1”) may be populated in either of two ways. The first uses steady-state data that may be collected already for the calibration of the estimation model data tables 74-84 as described above. This is done by additionally measuring ΔPWG,Dphr and ΔPWG=ΔPT=PEM−PT,out, and by indirectly measuring {dot over (m)}WG by the method described above in the section “Turbine and waste-gate flow indirect measurement methods”. The conversion from {dot over (m)}WG to {dot over (m)}WG@SonicStdCond is described above by equation (12). With the foregoing data measured (directly or indirectly), the table 122 (“TableC1”) may now be tabulated in the form of ΔPWG,Dphr=TableC1({dot over (m)}WG@SonicStdCond, ΔPWG).
The other way of calibrating table 122 (“TableC1”) is during actual engine operation, and is to change the values of ΔPWG,Dphr in table 122 (“TableC1”) until {dot over (m)}WG matches {dot over (m)}WGdes.
Populating the table 124 (“TableC2”) with data may be done during calibration fine tuning while performing transient maneuvers to obtain improved response without overshoot.
The data table 126 (“TableC3”), which is the calibration of the solenoid characteristics, may be performed without the engine running. The form is given as follows: DutyCycleWG,S In d=TableC3(ΔPWG,S In ddes, PBoost). The setup to characterize the solenoid valve as shown in FIG. 10. The operation of the valve may cause substantial pressure fluctuations, making accurate pressure measurement difficult. If so, one may add simple damper volumes. By adjusting the manual valve, one may reduce the compressed air source to a desired PBoost breakpoint (i.e., this is Pin in the FIG. 10). Then, change the waste-gate control solenoid duty cycle (i.e., this is DC in the FIG. 10) and log the delta pressure (dP). Data should be taken with a sufficiently fine solution since the surface is fairly non-linear. Finally, the captured data is inverted to fit the calibration format of table 126 (“TableC3”). Thus, the first step involves measuring data in the format as follows: dP=f(Pin, DC). The next step involves inverting the captured data into the data format of the data table 126, as follows: DC=f(dP, Pin).
The data table 128 (“TableC4”) may be populated with data as follows. First, for a constant vacuum, measure ΔPWG,Dphr versus DutyCycleWG,S In d and record the data. Then use the recorded data to populate the data table 128 (“TableC4”). Alternatively, if a variable pressure vacuum source is used, measure ΔPWG,Dphr versus DutyCycleWG,S In d and Pvac, and then record the data. A setup and approach similar to that used for data table 126 (“TableC3”) where compressed air is replaced with vacuum in FIG. 10.
The data table 130 (“TableC5”) may be populated with data as follows. First, the section describing the calibration of table 122 (“TableC1”) described the measurement of {dot over (m)}WG and how it is converted into {dot over (m)}WG@SonicStdCond One may, having collected the data for various waste-gate valve positions, ValvePosWGdes, may simply populate the data table 130 (“TableC5”).
It should be understood that electronic controller 14 as described above may include conventional processing apparatus known in the art, capable of executing pre-programmed instructions stored in an associated memory, all performing in accordance with the functionality described herein. That is, it is contemplated that the processes described herein will be programmed in a preferred embodiment, with the resulting software code being stored in the associated memory. Implementation of the present invention, in software, in view of the foregoing enabling description, would require no more than routine application of programming skills by one of ordinary skill in the art. Such an electronic controller may further be of the type having both ROM, RAM, a combination of non-volatile and volatile (modifiable) memory so that the software can be stored and yet allow storage and processing of dynamically produced data and/or signals.
It is to be understood that the above description is merely exemplary rather than limiting in nature, the invention being limited only by the appended claims. Various modifications and changes may be made thereto by one of ordinary skill in the art, which embody the principles of the invention and fall within the spirit and scope thereof.
Chart 1 correlates the various symbols/terms with its description:











Symbol
Description







R
Gas constant


kair
Ratio of specific heats, air


kexh
Ratio of specific heats, exhaust






  
    c
    
      p
      ,
      air
    
  
  =
  
    c
    p
    a
  




Specific heat at constant pressure, air






  
    c
    
      p
      ,
      exh
    
  
  =
  
    c
    p
    c
  




Specific heat at constant pressure, exhaust


Baro
Barometric pressure


Pamb
Ambient pressure


PC,in
Compressor inlet stagnation pressure






  P
  
    C
    ,
    in
  
  reference




Reference compressor inlet stagnation pressure (from turbo manufacturers data)


Pstd
Standard Pressure (a defined reference).


PC,out
Compressor outlet stagnation pressure






  P
  
    C
    ,
    out
  
  des




Desired compressor outlet stagnation pressure


PBoost
Boost pressure (after intercooler, before throttle)






  P
  Boost
  des




Desired Boost Pressure (after intercooler, before throttle)


MAP
Intake Manifold Absolute Pressure


PEM (= PT,in)
Exhaust manifold stagnation pressure






  P
  EM
  des




Desired exhaust manifold stagnation pressure


PT,in (= PEM)
Turbine inlet stagnation pressure


PT,out
Turbine outlet stagnation pressure






  P
  
    T
    ,
    out
  
  des




Desired turbine outlet stagnation pressure


ΔPAirfilter
Airfilter pressure drop


ΔPIntercooler
Intercooler pressure drop


Pvac
Pressure of vacuum source used for waste-gate



actuation.


ΔPWG
Waste-gate valve delta pressure






  
    ΔP
    WG
    des
  
  ⁢
  
    
  
  ⁢
  
    (
    
      =
      
        ΔP
        T
        des
      
    
    )
  




Desired waste-gate valve delta pressure


ΔPWG,Dphr
Waste-gate diaphragm delta pressure






  ΔP
  
    WG
    ,
    Dphr
  
  des




Desired waste-gate diaphragm delta pressure






  ΔP
  
    WG
    ,
    Dphr
  
  
    des
    ,
    base
  




Desired waste-gate diaphragm delta pressure, base look-up






  ΔP
  
    WG
    ,
    Dphr
  
  
    des
    ,
    TransientTerm
  




Desired waste-gate diaphragm delta pressure, transient term






  ΔP
  
    WG
    ,
    Dphr
  
  
    des
    ,
    CL
  




Desired waste-gate diaphragm delta pressure, closed-loop term






  ΔP
  
    WG
    ,
    
      S
      ⁢
      
        
      
      ⁢
      ln
      ⁢
      
        
      
      ⁢
      d
    
  
  des




Desired waste-gate solenoid delta pressure






  
    ΔP
    T
    des
  
  ⁢
  
    
  
  ⁢
  
    (
    
      =
      
        ΔP
        WG
        des
      
    
    )
  




Desired turbine delta pressure


PrC
Compressor pressure ratio






  Pr
  C
  des




Desired compressor pressure ratio


PrT
Turbine pressure ratio






  Pr
  T
  des




Desired turbine pressure ratio


Prexh
Exhaust system pressure ratio (catalyst, muffler)


Prcrit
Critical pressure ratio


Tamb
Ambient air temperature


TC,in (≅ Tamb)
Compressor inlet air stagnation temperature






  T
  
    C
    ,
    in
  
  reference




Reference compressor inlet air stagnation temperature (from turbo manufacturers data)


Tstd
Standard Temperature (a defined reference).


TEM (= TT,in)
Exhaust manifold gas stagnation temperature


TT,in (= TEM)
Turbine inlet stagnation temperature


TT,out
Turbine outlet stagnation temperature


TT,out,w/o HT
Turbine outlet stagnation temperature not accounting



for heat loss


{dot over (m)}C
Compressor mass flow rate






  
    m
    .
  
  C
  des




Desired compressor mass flow rate


{dot over (m)}T
Turbine mass flow rate


{dot over (m)}T,cor
Corrected turbine mass flow rate






  
    m
    .
  
  T
  des




Desired turbine mass flow rate






  
    m
    .
  
  
    T
    ,
    cor
  
  des




Desired corrected turbine mass flow rate


{dot over (m)}eng
Engine mass flow rate


{dot over (m)}eng,cor
Corrected engine mass flow rate






  
    m
    .
  
  
    cng
    ,
    air
  
  des




Desired engine air mass flow rate from torque control






  
    m
    .
  
  
    
      e
      ⁢
      ng
    
    ,
    air
  
  
    
      @
      
        
      
      ⁢
      Max
    
    ⁢
    
      
    
    ⁢
    des
    ⁢
    
      
    
    ⁢
    Boost
  




Engine air mass flow rate at maximum desired boost






  
    m
    .
  
  
    
      e
      ⁢
      ng
    
    ,
    air
  
  
    
      @
      des
    
    ⁢
    
      
    
    ⁢
    Boost
  




Engine air mass flow rate at desired boost


{dot over (m)}WG
Waste-gate mass flow rate






  
    m
    .
  
  WG
  des




Desired waste-gate mass flow rate






  
    m
    .
  
  WG
  
    des
    ⁢
    
      
    
    @
    
      
    
    ⁢
    SonicStdCond
  




Desired waste-gate mass flow rate at standard sonic conditions


{dot over (m)}exh
Exhaust system mass flow rate






  
    m
    .
  
  exh
  des




Desired exhaust system mass flow rate


{dot over (m)}EGR
EGR mass flow rate


{dot over (V)}C,cor
Corrected compressor volume flow rate






  
    V
    .
  
  
    C
    ,
    cor
  
  des




Desired corrected compressor volume flow rate


ηC
Compressor isentropic efficiency


ηT
Turbine isentropic efficiency


PC
Power absorbed by compressor






  P
  C
  des




Desired power absorbed by compressor


PT
Power produced by turbine


ΔhC
Compressor enthalpy delta






  Δh
  C
  des




Desired compressor enthalpy delta


{dot over (N)}T,cor
Corrected turbo-charger rotational speed






  
    N
    .
  
  
    T
    ,
    cor
  
  des




Desired corrected turbo-charger rotational speed


VE
Volumetric efficiency


 






  A
  =
  
    (
    
      
        P
        C
      
      
        
          η
          T
        
        ⁢
        
          c
          p
          e
        
        ⁢
        
          T
          
            T
            ,
            in
          
        
      
    
    )
  




Turbo Power Term (home made word)


 


Ades
Desired Turbo Power Term (home made word)


 






  B
  =
  
    
      
        T
        
          T
          ,
          in
        
      
    
    
      P
      
        T
        ,
        out
      
    
  




Turbine Boundary Term (home made word)


 


Bdes
Desired Turbine Boundary Term (home made word)


β2
Normalized valve flow dependency on pressure ratio



for compressible fluid


DutyCycleWG,S ln d
Waste-gate solenoid commanded duty-cycle






  ValvePos
  WG
  des




Desired waste-gate valve position

## Claims

1. A method of optimizing boost pressure of an internal combustion engine having a turbo-charger with a compressor and an exhaust driven turbine with a waste-gate flow path that bypasses the turbine in a parallel path, and a waste-gate mass flow rate {dot over (m)}WG through the waste-gate flow path being adiustable by a waste-gate valve, comprising the steps of:
providing an engine controller being in electrical communication with the engine; providing a desired boost pressure PBoostdes from the engine to the engine controller:
calculating a desired waste-gate mass flow rate {dot over (m)}WGdes in the waste-gate flow path with the engine controller, said calculated desired waste-gate mass flow rate {dot over (m)}WGdes being a function of the provided desired boost pressure PBoostdes ; and
generating the control signal from the controller being received by the waste-gate valve, wherein the engine controller is adapted to operatively control the waste-gate valve using the control signal such that the waste-gate mass flow rate {dot over (m)}WG through the waste-gate flow path attains the calculated desired waste-gate mass flow rate {dot over (m)}WGdes ,
wherein said step of calculating the desired waste-gate mass flow rate {dot over (m)}WGdes further includes the desired waste-gate mass flow rate {dot over (m)}WGdes being a desired waste-gate mass flow rate at standard sonic conditions {dot over (m)}WGdes@SonicStdCond across the waste-gate valve, and calculating the desired waste-gate mass flow rate at standard sonic conditions {dot over (m)}WGdes@SonicStdCond using a desired pressure ratio PrTdes includes the sub-steps of,
determining a desired compressor mass flow rate {dot over (m)}Cdes with the engine controller;
determining a desired compressor pressure ratio PrCdes using the provided determined, desired boost pressure PBoostdes with the engine controller;
determining a desired compressor power PCdes with the engine controller using the determined desired compressor mass flow rate {dot over (m)}Cdes and the determined desired compressor pressure ratio PrCdes and predetermined compressor characteristics data;
determining a desired corrected turbine mass flow rate {dot over (m)}T,cordes with the engine controller using the determined desired compressor power PCdes , the desired corrected turbine mass flow rate {dot over (m)}T,cordes being a function of a product term (A*B), and the product term (A*B) including a turbine power term (A) having the determined desired compressor power PCdes and a turbine boundary term (B) having a turbine outlet pressure PT,out ;
determining a desired turbine pressure ratio PrTdes with the engine controller using the determined desired corrected turbine mass flow rate {dot over (m)}T,cordes ; and
determining the desired waste-gate mass flow rate at standard sonic conditions {dot over (m)}WGdes@SonicStdCond with the engine controller using the desired turbine pressure ratio PrTdes.

2. The method of claim 1 wherein said step of providing the desired boost pressure PBoostdes further includes the desired boost pressure PBoostdes is a function of a desired engine torque.

3. The method of claim 1 wherein a position of the waste-gate valve is controlled by an electrical motor, said step of generating the control signal comprises the sub-steps of,
determining a waste-gate valve position ValvePosWGdes using the desired waste-gate mass flow rate at standard sonic conditions {dot over (m)}WGdes@SonicStdCond;
producing the control signal in accordance with said determined waste gate valve position; and
applying the control signal from the engine controller to the electrical motor.

4. The method of claim 3 further including the step of,
establishing a position-versus-flow mathematical model that correlates desired waste-gate mass flow rate at standard sonic conditions {dot over (m)}WGdes@SonicStdCond with waste-gate valve position ValvePosWGdes; and
obtaining a value for the waste-gate valve position ValvePosWGdes using the position-versus-flow mathematical model.

5. The method of claim 1 wherein the turbo-charger includes a variable geometry turbine (VGT) turbo-charger, and said step of generating the control signal further include controlling the waste-gate valve so as to achieve the desired waste-gate mass flow rate at standard sonic conditions {dot over (m)}WGdes@SonicStdCond being performed as a function of VGT position.

6. The method of claim 1 wherein the step of providing the engine controller further includes providing data tables disposed in a memory of the engine controller, and said data tables including data values measured for a plurality of engine system operating conditions associated with the engine, and the data tables comprising a predetermined compressor enthalpy delta table, a predetermined turbo corrected speed data table, a predetermined turbine efficiency (ηT) data table, a predetermined turbine pressure ratio/corrected engine flow data table, a predetermined product term (A*B) data table, a predetermined pressure data drop-exhaust system data table, and a predetermined heat transfer calibration table, wherein the engine controller uses the data tables for the calculating the desired waste-gate flow rate {dot over (m)}WGdes and the engine controller accesses at least one of the data tables after the step of providing the engine controller and prior to the step of generating the control signal.

7. The method of claim 1 further including,
configuring a control signal with the engine controller based on the calculated desired waste-gate mass flow rate {dot over (m)}WGdes to allow operative communication with the waste-gate valve before the step of generating the control signal,
wherein the step of configuring the control signal with the engine controller further includes providing data tables disposed in a memory of the engine controller, and said data tables including data values measured for a plurality of waste-gate value operating conditions associated with the waste-gate value in the engine, and the data tables comprising at least one of,
(i) a desired waste-gate valve position table,
(ii) a desired waste-gate diaphragm delta pressure ΔPwG,Dphr data table,
(iii) a transient response improvement term data table,
(iv) a desired waste-gate solenoid duty cycle data table where the waste-gate valve is pneumatically actuated using boost pressure adjusted through a solenoid, and
(v) a desired waste-gate solenoid duty cycle data table where the waste-gate valve is pneumatically actuated using a vacuum source adjusted through a solenoid, wherein the engine controller uses at least one of the data tables to further configure the control signal to operatively communicate with the waste-gate valve, and the engine controller accesses at least one of the data tables after the step of providing the engine controller.

8. The method of claim 1 wherein the calculated desired waste-gate mass flow rate {dot over (m)}WGdes is a calculated desired waste-gate mass flow rate at standard sonic conditions {dot over (m)}WGdes@SonicStdCond defined by the mathematical relationship




  
    
      m
      .
    
    WG
    
      des
      @
      SonicStdCond
    
  
  =
  
    
      
        m
        .
      
      WG
      des
    
    ⁢
    
      
        P
        std
      
      
        P
        EM
        des
      
    
    ⁢
    
      
        
          T
          EM
        
        
          T
          std
        
      
    
    ⁢
    
      1
      
        
          β
          2
        
        ⁡
        
          (
          
            1
            /
            
              Pr
              T
              des
            
          
          )
        
      
    
  



where
{dot over (m)}WGdes@SonicStdCond = the calculated desired waste-gate mass flow rate at standard sonic conditions,
{dot over (m)}WGdes@SonicStdCond = the desired waste-gate mass flow rate,
Pstd= a Standard Pressure (a defined reference),
PEMdes= a desired exhaust manifold stagnation pressure,
TEM= a Exhaust manifold gas stagnation temperature,
Tstd= a Standard Temperature (a defined reference),
β2= a normalized valve flow dependency on pressure ratio for compressible fluid, and
PrTdes= the desired turbine pressure ratio.

9. A method of optimizing boost pressure of an internal combustion engine having a turbo-charger with a compressor and an exhaust driven turbine with a waste-gate flow path that bypasses the turbine in a parallel path, and a waste-gate mass flow rate {dot over (m)}WG through the waste-gate flow path being adiustable by a waste-gate valve, comprising the steps of:
providing an engine controller being in electrical communication with the engine; providing a desired boost pressure PBoostdes from the engine to the engine controller:
calculating a desired waste-gate mass flow rate {dot over (m)}WGdes in the waste-gate flow path with the engine controller, said calculated desired waste-gate mass flow rate {dot over (m)}WGdes being a function of the provided desired boost pressure PBoostdes ; and
generating the control signal from the controller being received by the waste-gate valve, wherein the engine controller is adapted to operatively control the waste-gate valve using the control signal such that the waste-gate mass flow rate {dot over (m)}WG through the waste-gate flow path attains the calculated desired waste-gate mass flow rate {dot over (m)}WGdes,
wherein the waste-gate valve is pneumatically actuated using boost pressure adiusted through a solenoid, said step of generating the control signal with the engine controller comprises the sub-steps of,
determining a duty cycle DutyCycleWG,S In d with the engine controller using the desired waste-gate mass flow rate at standard sonic conditions {dot over (m)}WGdes@SonicStdCond ;
producing the control signal in accordance with said determined duty cycle; and
applying the control signal from the engine controller to the solenoid,
wherein the waste-gate valve is linked to and controlled by a waste-gate diaphragm and where the solenoid adiusts the boost pressure applied to the diaphragm, said sub-step of determining a duty cycle further including includes the sub-step of,
determining a desired pressure delta across the waste-gate diaphragm ΔPWG,Dphrdes,
wherein said step of determining the desired pressure delta across the waste-gate diaphragm ΔPWG,Dphrdes further includes,
establishing a first mathematical model that correlates the desired waste-gate mass flow rate at standard sonic conditions {dot over (m)}WGdes@SonicStdCond and desired waste gate valve delta pressure ΔPWGdes being a function of a desired base pressure delta across the waste-gate diaphragm ΔPWG,Dphrdes,base; and
obtaining a value for the desired base pressure delta across the waste-gate diaphragm ΔPWG,Dphrdes,base using the first mathematical model,
wherein said step of obtaining the value for the determining a desired base pressure delta across the waste-gate diaphragm ΔPWG,Dphrdes,base further includes,
establishing a second mathematical model that correlates a difference between boost ressure PBoost and the desired boost pressure PBoostdes to a transient response pressure delta across the waste-gate diaphragm ΔPWG,Dphrdes,TransientTerm; and
obtaining a value for the transient response pressure delta across the waste-gate diaphragm ΔPWG,Dphrdes,TransientTerm using the second mathematical model,
wherein said step of obtaining the value for the desired pressure delta across the waste-gate diaphragm PWG,Dphrdes is an arithmetic summation of ΔPWG,Dphrdes,base and ΔPWG,Dphrdes,TransientTerm and ΔPWG,Dphrdes,CL, where ΔPWG,Dphrdes,CL is a closed loop correction term to correct model inaccuracies.

10. The method of claim 9 wherein said step of determining the duty cycle DutyCycleWG,S In d further includes,
determining a desired waste-gate solenoid pressure delta ΔPWG,S In ddes.

11. The method of claim 10 wherein said step of determining a desired waste-gate solenoid pressure delta ΔPWG,S In ddes further includes,
subtracting an ambient pressure Pamb and the waste-gate diaphragm pressure delta ΔPWG,Dphrdes from the desired boost pressure PBoostdes to obtain the desired waste-gate solenoid pressure delta ΔPWG,S In ddes.

12. The method of claim 11 wherein said step of determining the duty cycle DutyCycleWG,S In d further includes,
establishing a third data structure correlating boost pressure PBoost and the desired waste-gate solenoid pressure delta ΔPWG,S In ddes to the duty cycle DutyCycleWG,S In d; and
obtaining a value for the duty cycle DutyCycleWG,S In d using the third mathematical model.

13. A method of optimizing boost pressure of an internal combustion engine having a turbo-charger with a compressor and an exhaust driven turbine with a waste-gate flow path that bypasses the turbine in a parallel path, and a waste-gate mass flow rate {dot over (m)}WG through the waste-gate flow path being adiustable by a waste-gate valve, comprising the steps of:
providing an engine controller being in electrical communication with the engine; providing a desired boost pressure PBoostdes from the engine to the engine controller;
calculating a desired waste-gate mass flow rate {dot over (m)}WGdes in the waste-gate flow path with the engine controller, said calculated desired waste-gate mass flow rate {dot over (m)}WGdes being a function of the provided desired boost pressure PBoost des; and
generating the control signal from the controller being received by the waste-gate valve, wherein the engine controller is adapted to operatively control the waste-gate valve using the control signal such that the waste-gate mass flow rate {dot over (m)}WG through the waste-gate flow path attains the calculated desired waste-gate mass flow rate {dot over (m)}WGdes,
wherein the waste-gate valve is pneumatically actuated using a vacuum source adiusted through a solenoid, said step of generating the control signal comprises the sub-steps of,
determining a duty cycle DutyCycleWG,S In d using the desired waste-gate mass flow rate at standard sonic conditions {dot over (m)}WGdes@SonicStdCond ;
producing the control signal in accordance with the determined duty cycle;
applying the control signal from the engine controller to the solenoid,
wherein the waste-gate valve is linked to and controlled by a waste-gate diaphragm and where the solenoid adiusts the vacuum applied to the diaphragm, said sub-step of determining the duty cycle includes the sub-step of,
determining a desired pressure delta across the waste-gate diaphragm ΔPWG,Dphrdes,
wherein said step of determining the desired pressure delta across the waste-gate diaphragm ΔPWG,Dphrdes includes,
establishing a first mathematical model that correlates desired waste-gate mass flow rate at standard sonic conditions {dot over (m)}WGdes@SonicStdCond and desired waste-gate valve delta pressure ΔPWGdes to a desired base pressure delta across the waste-gate diaphragm ΔPWG,Dphrdes,base;and
obtaining a value for the desired base pressure delta across the waste-gate diaphrag ΔPWG,Dphrdes,base using the first mathematical model,
wherein said step of obtaining the value for the determining a desired base pressure delta across the waste-gate diaphragm ΔPWG,Dphrdes,base includes,
establishing a second mathematical model that correlates a difference between boost ressure PBoost and the desired boost pressure PBoostdes to a transient response pressure delta across the waste-gate diaphragm ΔPWG,Dphrdes,TransientTerm ; and
obtaining a value for the transient response pressure delta across the waste-gate diaphragm ΔPWG,Dphrdes,TransientTerm using the second mathematical model, desired pressure delta across the waste-gate ΔPWG,Dphrdes is an arithmetic summation of ΔPWG,Dphrdes,base and ΔPWG,Dphrdes,TransientTerm and ΔPWG,Dphrdes,CL, where ΔPWG,Dphrdes,CL is a closed loop correction term to correct model inaccuracies.

14. The method of claim 13 wherein said step of determining the duty cycle DutyCycleWG,S In d includes,
determining a desired waste-gate solenoid pressure delta ΔPWG,S In ddes.

15. The method of claim 14 wherein said step of determining the desired waste-gate solenoid pressure delta ΔPWG,S In ddes includes,
subtracting an ambient pressure Pamb and the waste-gate diaphragm pressure delta ΔPWG,Dphrdes from a vacuum source pressure Pvac to obtain the desired waste-gate solenoid pressure delta ΔPWG,S In ddes.

16. The method of claim 15 wherein said step of obtaining the duty cycle DutyCycleWG,S In d further includes,
establishing a vacuum solenoid mathematical model correlating desired waste-gate solenoid pressure delta ΔPWG,S In ddes and vacuum source pressure Pvac with the duty cycle DutyCycleWG,S In d; and
obtaining a value for the duty cycle DutyCycleWG,S In d using the vacuum solenoid mathematical model.

17. The method of claim 15 wherein said step of obtaining the duty cycle DutyCycleWG,S In d further includes,
establishing a vacuum solenoid mathematical model Correlating desired waste-gate diaphragm delta pressure ΔPWG, Dphrdes with the duty cycle DutyCycleWG,S In d; and
obtaining a value for the duty cycle DutyCycleWG,S In d using the vacuum solenoid mathematical model.