Research Article  Open Access
Estimation of Individual Cylinder AirFuel Ratio in Gasoline Engine with Output Delay
Abstract
The estimation of the individual cylinder airfuel ratio (AFR) with a single universal exhaust gas oxygen (UEGO) sensor installed in the exhaust pipe is an important issue for the cylindertocylinder AFR balancing control, which can provide highquality torque generation and reduce emissions in multicylinder engine. In this paper, the system dynamic for the gas in exhaust pipe including the gas mixing, gas transport, and sensor dynamics is described as an output delay system, and a new method using the output delay system observer is developed to estimate the individual cylinder AFR. With the AFR at confluence point augmented as a system state, an observer for the augmented discrete system with output delay is designed to estimate the AFR at confluence point. Using the gas mixing model, a method with the designed observer to estimate the individual cylinder AFR is presented. The validity of the proposed method is verified by the simulation results from a spark ignition gasoline engine from engine software enDYNA by Tesis.
1. Introduction
Cylindertocylinder airfuel ratio (AFR) balancing control in internal combustion engines with multiple cylinders is one of the technology trends to satisfy the increasingly stringent emission regulations, which can also improve the engine performance, such as thermal efficiency and fuel economy. The AFR of each cylinder is decided by the aspirated air mass, the injected fuel mass, and the residual gas from the prior cycle, in which the combustion stroke of each cylinder sequentially occurs along the rotational angle of the crankshaft. Due to the air breathing variability and injector variability, there exists AFR imbalance between cylinders, leading to adverse impacts on emission performance using the conventional controllers [1, 2].
In order to improve the AFR control accuracy, there has been a great deal of research that focuses on the AFR control of individual cylinders [3–10]. In fact, the estimation of the individual cylinder AFR with a single universal exhaust gas oxygen (UEGO) sensor installed in the exhaust pipe is one of the key technology trends for the individual cylinder AFR control. The digital filtering techniques are employed to extract the AFR imbalance signals from oxygen sensor voltage signals in [11], in which the oxygen sensor voltage signal is processed to determine imbalanced cylinder identification and AFR cylinder imbalance levels. In [12], a modeling method to estimate the individual fuelgas ratio is proposed to estimate AFR, which is used for an adaptive generalized predictive control approach to balance the individual cylinder characteristics in the static engine operation mode. A static steady state observer based on the individual cylinder AFR model along the air mass flow, gas mixing, and sensor dynamics in an exhaust manifold can be found in [13]. In the diesel engines, a nonlinear observer is proposed to estimate the individual cylinder AFR [14]. In [15], a PI compensator is designed to compensate the cylinderbycylinder variations, in which an input observer is proposed to estimate individual cylinder AFR.
However, the transport delay and sensor delay from the exhaust confluence point to UEGO sensor output are ignored in the proposed observers from the above papers, which may reduce the accuracy of the AFR estimation of each cylinder. In order to improve the individual cylinder AFR estimation accuracy, the system dynamics in the exhaust pipe including gas transport and sensor dynamics is described as an augmented discrete system with output delay in this paper, in which the AFR at confluence point is augmented as a system state. Then, an observer for the augmented discrete system with output delay is designed. With the combination of the designed observer and the gas mixing model at confluence point, the method to estimate the individual cylinder AFR is presented. The performance of the proposed method is validated against the simulation result from engine software enDYNA provided by Tesis, and a comparison with existing method is given during an urban driving cycle, which demonstrates that the proposed method can improve the accuracy of the individual cylinder AFR estimation.
This paper is organized as follows. In Section 2, the system dynamics in the exhaust pipe including gas transport and UEGO sensor dynamics is described as an augmented system with output delay. In Section 3, an observer for the output delay system is designed, and the method to estimate the individual cylinder AFR is presented. Simulation results from enDYNA are presented in Section 4, and the conclusions are summarized in Section 5.
2. Problem Formulation
A schematic diagram of a 4cylinder sparkignited (SI) gasoline engine is shown in Figure 1, where the fuel injectors equipped at the inlet port near to the intake valve are controlled individually. The fuel mass burnt in each cylinder is injected by the corresponding injector, and the fuel injection command is delivered to the injector of each cylinder serially along the crank angle. The AFR of each cylinder is , where is the fuel mass into the cylinder and is the air mass into the cylinder.
After combustion, the combusted gas of each cylinder is exhausted into the corresponding runner during the exhaust stroke of each cylinder and passes through their runners and confluence in the public exhaust manifold. Then, the AFR of the mixed gas is measured by a UEGO sensor, and the mixed gas runs to the outside passing through the catalyst. The system dynamics in the exhaust pipe includes the gas mixing, gas transport, and sensor dynamics, in which the transfer function from the confluence point to sensor output can be given by [13]where is the AFR at the exhaust confluence point, is the measured AFR of the UEGO sensor, is the time delay including the transport delay and sensor delay , is the time constant of the mixing process, and is the sensor time constant.
For an engine with 4 cylinders, the combined UEGO sensor signal is sampled at the exhaust top dead center of each cylinder, and the sampling period is related to the engine cycle period as , in which the engine cycle period is , where is the engine speed in rpm. Consisting of the zeroorder holder and the mixing and sensor dynamics, the discrete form of model (1) with the sampling period can be given by [13]whereFurthermore, (2) can be written as a statespace form with input delay [15]:where
Define the new system state as , and (4) with input delay can also be rewritten as the following discrete system with output delay:
in (6) is unknown, which can be considered as a system state due to the small rate of change of . Then, an augmented discrete system with output delay from (6) can be obtained:where
Equation (7) indicates that the estimation of the AFR at the exhaust confluence point becomes the state estimation of the discrete system with output delay (7).
3. Observer Design for Discrete System with Output Delay
The observer for the output delay system (7) is given:where is the state estimate and is the feedback gain matrix. The asymptotical stability of the proposed observer (9) is analyzed in the following theorem.
Theorem 1. There exist matrices , , , and , such that the following linear matrix inequality (LMI) is feasible: then observer (9) is asymptotically stable.
Proof. Set the estimation error , and the error dynamic system between (7) and (9) is obtained:Denote ; then we haveChoose a Lyapunov functional candidate asDefine ; then along the solution of (11) and (12) we haveAccording to Lemma in [16], for any constant matrix , , and , the following inequality holds: With the combination of (14) and (15), we havewhereBy Schur complement [17], the following LMI (18) guarantees , which can guarantee and the asymptotical stability of observer (9):Now, condition (10) guaranteeing (18) must be proved. Define , and we haveBecause of the fact that , we have . Therefore, condition (10) can guarantee ; then (18) holds.
The estimation of the AFR at the exhaust confluence point can be obtained according to observer (9). In order to obtain the AFR of each cylinder, the relationship between the AFR of each cylinder and the AFR at the exhaust confluence point is analyzed in the following.
The combusted gas of each cylinder is discharged into the corresponding exhaust port and flows to the exhaust confluence point in the exhaust manifold, in which we assume that exhaust gas mixing in the individual exhaust runner can be neglected. Hence, the AFR in the exhaust runner is constant during one engine cycle, and the AFR at the exhaust confluence point in the th engine cycle can be given by [15]where is the exhaust air mass flow in the th exhaust manifold passing through the confluence point and is the AFR in the th cylinder. Furthermore, under the assumption that air mass aspirated into cylinders in each cycle is constant, exhaust air flow has the same shape between successive cycles; then a periodic function can be obtained:Therefore, the gas mixing behavior (20) can be rewritten in the domain aswhere . The relationship between the AFR of each cylinder and the AFR at the exhaust confluence point can be obtained by (22).
The algorithm to estimate the individual cylinder AFR is as follows: First, the AFR at the exhaust confluence point is obtained according to observer (9). Then, the individual cylinder AFR can be calculated through (22). With the combination of (9) and (22), the method for the estimation of each cylinder AFR can be given as follows:
4. Simulation Studies
In this section, the simulation study of the estimation of the individual cylinder AFR is presented in the environment of a 2.0 L 4cylinder SI gasoline engine from enDYNA [18, 19]. The enDYNA is a professional software tool for the realtime simulation of internal combustion engines, providing readytouse models for all common engine types comprising crank angle synchronous combustion, gas path, fuel system, cooling system, drivetrain, driver, and softECU. The R4cylinder SIengine is an example in enDYNA to simulate a 4cylinder SI gasoline engine, whose specifications are given in Table 1. The observer architecture is illustrated in Figure 2.

The parameters of the system dynamics in the exhaust pipe (2) and the gas mixing (22) are presented as follows [15]:
Then, the system matrix can be obtained:
According to the inequality (10), the gain matrix can be given by .
Here, the input of the throttle angle in enDYNA is designed as a step signal presented in Figure 3. Accordingly, the estimation results of the individual cylinder AFR by the proposed method are shown in Figure 3, and the estimation errors are plotted in Figure 4. When the throttle changes suddenly, the estimation errors of the individual cylinder AFR are about 0.06%, and then the steady state errors tend to decay within 0.01%.
In order to verify the effectiveness of the proposed method under driving cycle condition, one segment of the urban driving cycle ECE (Economic Commission for Europe) is used [19], under which the throttle angle , engine speed , intake manifold pressure , and vehicle velocity are plotted in Figure 5. Accordingly, the comparison of the individual cylinder AFR estimation between the proposed method and the input observer in [15] are presented in Figure 6, and the estimation errors are plotted in Figure 7. Clearly, the error of the proposed method is smaller than the input observer when the AFR changes slowly, in which the steady state error of the input observer is 0.03%. When the AFR changes severely, there exist fluctuations of the AFR estimation error from both the proposed method and input observer in [15]. However, the estimation error from the proposed method is smaller. It is demonstrating that the proposed method considering time delay can improve the accuracy of the individual cylinder AFR estimation.
5. Conclusion
An efficient method for the estimation of the individual cylinder AFR with a single UEGO sensor was developed to improve the estimation accuracy. The system dynamics in the exhaust pipe was described as an augmented discrete system with output delay, in which the AFR at confluence point was augmented as a system state and beneficial to be estimated comparing the system with input delay. Then, an observer for the augmented system with output delay was designed to estimate the AFR at confluence point, which can avoid accurately inverting the engine model including delays. Using the gas mixing model, a method to estimate the individual cylinder AFR based on the proposed observer was presented. The performance of the proposed method was validated by the simulation data from engine software enDYNA provided by Tesis, and a comparison with existing method was obtained during ECE cycle, demonstrating that the proposed method considering time delay from exhaust gas transport and UEGO sensor dynamics can improve the accuracy of the individual cylinder AFR estimation.
Competing Interests
The authors declare that they have no competing interests.
Acknowledgments
This work was supported by China Automobile Industry Innovation and Development Joint Fund (no. U1564213).
References
 J. B. Heywood, Internal Combustion Engine Fundamentals, vol. 930, McGrawHill, New York, NY, USA, 1988.
 J. W. Grizzle, K. L. Dobbins, and J. A. Cook, “Individual cylinder airfuel ratio control with a single EGO sensor,” IEEE Transactions on Vehicular Technology, vol. 40, no. 1, pp. 280–286, 1991. View at: Publisher Site  Google Scholar
 J. Chauvin, P. Moulin, G. Corde, N. Petit, and P. Rouchon, “Realtime nonlinear individual cylinder air fuel ratio observer on a diesel engine test bench,” in Proceedings of the 16th Triennial World Congress of International Federation of Automatic Control (IFAC '05), pp. 194–199, Prague, Czech Republic, July 2005. View at: Google Scholar
 J. Chauvin, P. Moulin, G. Corde, N. Petit, and P. Rouchon, “Kalman filtering for realtime individual cylinder air fuel ratio observer on a diesel engine test bench,” in Proceedings of the American Control Conference, pp. 1886–1891, IEEE, Minneapolis, Minn, USA, June 2006. View at: Google Scholar
 J. Chauvin, N. Petit, and P. Rouchon, “Six degrees crankshaft individual air fuel ratio estimation of diesel engines for cylinder balancing purpose,” Sae Technical Papers, 2006. View at: Google Scholar
 Y. Liu and T. Shen, “Modeling and experimental validation of airfuel ratio under individual cylinder fuel injection in gasoline engines,” IEEJ Journal of Industry Applications, vol. 1, no. 3, pp. 155–163, 2012. View at: Publisher Site  Google Scholar
 J. F. Burkhard, “Individual cylinder fuel control for a turbocharged engine,” SAE Technical Paper 01487191, SAE International, 2014. View at: Google Scholar
 H. Li, Y. Huang, G. Li, and Y. Yang, “Research on the cylinderbycylinder variations detection and control algorithm of diesel engine,” SAE Technical Paper 01487191, 2015. View at: Google Scholar
 S. Nakagawa, A. Numata, and T. Hori, “Individual cylinder control for airfuel ratio cylinder imbalance,” SAE Technical Paper 2015011624, 2015. View at: Publisher Site  Google Scholar
 M. Kassa, C. Hall, A. Ickes, and T. Wallner, “Incylinder oxygen mass fraction estimation method for minimizing cylindertocylinder variations,” SAE Technical Paper 01487191, 2015. View at: Google Scholar
 N. Qiao, C. Krishnamurthy, and N. Moore, “Determine airfuel ratio imbalance cylinder identification with an oxygen sensor,” SAE International Journal of Engines, vol. 8, no. 3, pp. 1005–1011, 2015. View at: Publisher Site  Google Scholar
 K. Suzuki, T. Shen, J. Kako, and S. Yoshida, “Individual A/F estimation and control with the fuelgas ratio for multicylinder IC engines,” IEEE Transactions on Vehicular Technology, vol. 58, no. 9, pp. 4757–4768, 2009. View at: Publisher Site  Google Scholar
 L. Benvenuti, M. D. Di Benedetto, S. Di Gennaro, and A. SangiovanniVincentelli, “Individual cylinder characteristic estimation for a spark injection engine,” Automatica, vol. 39, no. 7, pp. 1157–1169, 2003. View at: Publisher Site  Google Scholar  MathSciNet
 J. Chauvin, G. Corde, P. Moulin, N. Petit, and P. Rouchon, “High frequency individual cylinder estimation for control of diesel engines,” Oil and Gas Science & Technology, vol. 61, no. 1, pp. 57–72, 2006. View at: Publisher Site  Google Scholar
 B. He, T. Shen, J. Kako, and M. Ouyang, “Input observerbased individual cylinder airfuel ratio control: modelling, design and validation,” IEEE Transactions on Control Systems Technology, vol. 16, no. 5, pp. 1057–1065, 2008. View at: Publisher Site  Google Scholar
 W.A. Zhang and L. Yu, “Stability analysis for discretetime switched timedelay systems,” Automatica, vol. 45, no. 10, pp. 2265–2271, 2009. View at: Publisher Site  Google Scholar  MathSciNet
 S. P. Boyd, L. El Ghaoui, E. Feron, and V. Balakrishnan, Linear Matrix Inequalities in System and Control Theory, vol. 15, SIAM, 1994.
 TESIS DYNAware, enDYNA® THERMOS® 2.0 Block Reference Manual, 2006.
 T. DYNAware, enDYNAⓇ THERMOSⓇ 2.0 User Manual, 2006.
Copyright
Copyright © 2016 Changhui Wang and Zhiyuan Liu. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.