Combustion Model and Control Parameter Optimization Methods for Single Cylinder Diesel Engine

This research presents a method to construct a combustion model and a method to optimize some control parameters of diesel engine in order to develop a model-based control system. The construction purpose of the model is to appropriately manage some control parameters to obtain the values of fuel consumption and emission as the engine output objectives. Stepwise method considering multicollinearity was applied to construct combustion model with the polynomial model. Using the experimental data of a single cylinder diesel engine, the model of power, BSFC, NO x , and soot on multiple injection diesel engines was built. The proposed method succesfully developed the model that describes control parameters in relation to the engine outputs. Although many control devices can be mounted to diesel engine, optimization technique is required to utilize this method in finding optimal engine operating conditions efficiently beside the existing development of individual emission control methods. Particle swarm optimization (PSO) was used to calculate control parameters to optimize fuel consumption and emission based on the model. The proposed method is able to calculate control parameters efficiently to optimize evaluation item based on the model. Finally, the model which added PSO then was compiled in a microcontroller.


Introduction
Currently, energy and environmental pollution have become hot topics.Public and private transportations system plays an important role in people's life.Vehicle has become one of the most important factors in this system.In Jakarta, it is reported that 70% exhaust gas pollution is from vehicles [1].
Internal combustion engine is the main power of the vehicle.The global challenge in internal combustion engine is the reduction of exhaust gas such as soot, NO  [2,3], CO, HC, and particulate matter (PM) emissions from vehicles by improving fuel consumption without sacrificing the vehicle performance.On the other side, most internal combustion engine methods that reduce either NO  or soot emissions cause an increase in the other emissions [4].Recently, with the environmental restrictions and sustainable development, pollution standards have been more and more stringent.For example, in Japan, the standard of diesel engine emission regulation is very high.In 2003 the MOE finalized very stringent 2005 emission standards for both light and heavy vehicles although it remained relaxed through the 1990s.In 2005 heavy-duty emission standards (NO  = 2 g/kWh, PM = 0.027 g/kWh) were the most stringent diesel emission regulations in the world.In 2009, this standards were more tightened (NO  = 0.7 g/kWh, PM = 0.01 g/kWh) to a level in-between the US 2010 and Euro V requirements [5].So engine pollutant emission reduction became a major interest for engine development.
Additionally, a technology to decrease the exhaust gas emission often increases brake specific fuel consumption (BSFC) and causes deprivation in fuel economy, which in turn increases the emissions of CO 2 .Thus, it is of importance to reduce the emissions while improving BSFC.
There are some methods that have been proposed to reduce diesel engine exhaust gas emissions.One of them is method for in-cylinder control, that is, exhaust gas recirculation (EGR) which can reduce NO  emissions by reducing the combustion temperatures [6,7].The other method is ultrahigh injection pressure which can reduce soot emissions [8].However, due to high combustion temperatures associated with high injection pressures, NO  emissions are likely to increase.
The current diesel engine has been equipped with a number of control devices such as multiple injection equipment with common-rail system and turbocharger [9].Appropriate configurations of multiple injections can reduce soot emissions without a significant increase in NO  emissions [10].
The previous methods individually have been proven to work toward the reduction of a particular pollutant.Evaluation of combinations of these methods, however, is a heavy task considering the total number of experiments that need to be performed.
Although many control devices can be mounted in diesel engine, a technology which sets the multiple control parameter of these control devices optimally is needed.In order to control the large number of control parameters appropriately by considering the fuel consumption and exhaust gas components as the engine output objectives, a combustion model construction which reproduces the characteristic value of fuel consumption and exhaust gas components from control parameter is needed.In order to know the characteristics of the combustion, the construction of a combustion model to reproduce the characteristic values of fuel consumption and exhaust gas components from the control parameters is called for.Furthermore, in order to improve the fuel efficiency and to reduce the exhaust gas emission, the control parameter optimization is needed.
This research has target to construct a combustion model, to get the optimal control parameter, and to validate in the real diesel engine experimental device.In this study, we applied the stepwise method considering multicollinearity to construct a combustion model.In the construction of a combustion model, the method of approximating by a polynomial model based on experimental data was used.We used the experimental data of single cylinder diesel engine to build a predictive model of fuel consumption and exhaust gas in multiple injections.Furthermore, PSO (particle swarm optimization) method was used to optimize the control parameters and validated in diesel engine.Finally, we evaluated its performance with validation in diesel engine.

Methodology
To increase the fuel efficiency and decrease exhaust emission, diesel engine optimization on model-based control system is proposed.In model-based control system, engine conditions are directly calculated by models, which achieve a time benefit and improve accuracy in comparison with conventional maps.
The methodology of this research consists of several steps.First, we get the data from the diesel engine laboratory to build the combustion model.Second, we build the combustion model by stepwise method considering multicollinearity using diesel engine data.Third, we get the optimal value of control parameter by the PSO.Fourth, we validate the control performance in diesel engine laboratory.Finally, we compile the combustion model which added PSO in a microcontroller, and control performance is evaluated.For more details, we can see the schematic of diesel engine system in Figure 1.

Single Cylinder Diesel Engine Experiment
3.1.Specification of Diesel Engine.In this research, we developed an experimental device for studying the control technology to improve the fuel consumption and reduce the diesel engine emission.The diesel engine experimental device is in Figure 2.
In this research, we choose a four-cycle one-cylinder diesel engine.This research used Yanmar TF70 V-E diesel engine with 4 cycle horizontal type water-cooling and equipped with a turbocharger (in Figure 2).The specification of single cylinder diesel engine is as in Table 1.

Experiment Condition of Multistage Injection Diesel
Engine.In this research, the diesel engine is set with three stage injection, that is, pilot1 injection, pilot2 injection, and main injection as shown in Figure 3.Then, the rotation speeds were 1000 rpm, 1500 rpm, and 2000 rpm; EGR rate was 0% and the engine temperature was set at about 90 ∘ C degrees.The engine control parameters were set as in Table 2 and the engine optimization objectives are listed as in Table 3.
The experiment was done under the condition mentioned above and the experiment results will be discussed in the next part.

Stepwise Method.
Stepwise method is a systematic method for adding and removing terms from a multilinear model based on their statistical significance in a regression [11].It uses statistical methods to remove the redundant variables.In this study, it is to identify good model with considerably less computing than is required for all possible regressions.The main approaches in this method are forward selection, backward elimination, and bidirectional elimination.Forward selection is as follows: starting with no variables in the model, testing the addition of each variable using a chosen model comparison criterion, adding the variable that improves the model the most, and repeating this process until none improves the model.Backward elimination is as follows: starting with all candidate variables, testing the deletion of each variable using a chosen model comparison criterion, deleting the variable that improves the model the most by being deleted, and repeating this process until no further improvement is possible.Bidirectional elimination is as follows: a combination of the above, testing at each step  for variables to be included or excluded.In this study, we use bidirectional elimination to identify the model.
The model expresses the value of a response variable as a multilinear function of one or more predictor variables in order 1 and order 2: where   is response variable in observation ,  0 is coefficient constant,   is coefficient on the   predictor,   is coefficient on the   predictor and   predictor,   is coefficient on the   predictor 2 order,  is the total number of predictors, and   is error term.The model ( 1) is estimated by least squares method, which yields parameter estimates such that the sum of squares of errors is minimized.In order to select a clause effectively in presumption, we redefine one explaining variable about all the order 2, respectively.For example, we redefine  1 =  1 2 and  2 =  1  2 .Thus, the polynomial model selects the combination of an explaining variable effectively in presumption and constituted by the stepwise method.In this research, we use the multiple correlation coefficient  that indicated the matching level of the calculation datum by the regression equation and the original datum; the result is better when  is closer to 1. Statistic values  indicate the significance of the regression equation, whose values obey  distribution.

The Stepwise Procedure with Consideration of Multicollinearity.
The general procedure of the stepwise method consists of three steps.First, we select an initial regression model.Second, the procedure repeatedly alters the model by adding or removing a predictor variable in accordance with the -test [12].Finally, the search is terminated when the response variables which satisfy the stepping criterion do not exist anymore, or when the iteration step has reached a specified maximum number.The flow chart of the stepwise model-building procedures is shown in Figure 4.
Multicollinearity is a term reserved to describe the case when the intercorrelation of predictor variables is high.Multicollinearity does not invalidate the regression model in the sense that the predictive value of the equation may still be good as long as the prediction is based on combinations of predictors within the same multivariate space used to calibrate the equation.But there are several negative effects of multicollinearity such as the variance of the regression coefficients which can be inflated so much that the individual coefficients are not statistically significant even though the overall regression equation is strong and the predictive ability is good.
The variance inflation factor (VIF) is a statistic that can be used to identify multicollinearity [13] in a matrix of predictor variables.Variance inflation refers here to the mentioned effect of multicollinearity on the variance of estimated regression coefficients.Multicollinearity depends not just on the bivariate correlations between pairs of predictors, but on the multivariate predictability of any one predictor from the other predictors as well.Accordingly, the VIF is based on the multiple correlation coefficients in regression of each predictor in multivariate regression on all the other predictors: (2)

Particle Swarm Optimization
PSO algorithm is an adaptive algorithm based on a socialpsychological metaphor; a population of individuals (referred to as particles) adapts by returning stochastically toward previously successful regions [14].Particle swarm has two primary operators: velocity update and position update.During each generation, each particle is accelerated toward the particles previous best position and the global best position.
For each iteration a new velocity value for each particle is calculated based on its current velocity, the distance from its previous best position, and the distance from the global best position.The new velocity value is then used to calculate the next position of the particle in the search space.This process is then iterated for a set number of times, or until a minimum error is achieved.The flow chart of particle swarm optimization algorithm is shown in Figure 5.

Combustion Model Compilation in Microcontroller
In this research, we want to compile the combustion model which added PSO result in microcontroller and install in diesel engine to do real experiment/validation in diesel engine.First, we have to create a Simulink model using MAT-LAB and Simulink.Instead of programming C code manually, we can implement the control algorithm graphically using Simulink blocks.The models are saved as MDL files.
In this research we built combustion model with stepwise method and optimized by PSO.All of this work was done by Simulink and MATLAB.Finally, we get the comparison of validation result and the estimation result of an engine control parameter.

Yes
Return the q = 0 q < p the variable and all variable of Q and then calculate the F value q = 0 in Q to P. q = q − 1, p = p + 1 Calculate the minimum F value and then define without variable and then calculate decreased F value by this removal Is the all VIF less than threshold ?
Were all variable selected from P?
Make VIF between the variable of F max and each variable of Q Add variable of F max in P to Q. q = q + 1, p = p − 1 Were all variable selected from Q? (3) In the same way we obtain the combustion model of BSFC, NO  , and soot constructed by stepwise method considering multicollinearity: BSFC = 37.745 − 57.0318 5 + 4.9796 − 002 ( = 0.95147 and 7 variables) ( = 0.99603 and 6 variables) .
In this research, we have reported our predicted model of power, BSFC, NO  , and soot in multiple injection diesel engines by stepwise method considering multicollinearity.It shows that the predictive accuracy by stepwise method considering multicollinearity is high.This is proved by the multiple correlation coefficient 0.9000 or more and  value is very high.It can be regarded that the stepwise method can effectively estimate the objectives.In order to improve exhaust emissions and fuel efficiency in a diesel engine, in next step, we use particle swarm optimization (PSO), one of optimization techniques to find the optimal engine operating condition efficiently.

The Optimization of an Engine Control
Parameter.Before getting the optimal value of engine control parameter, we have to create a Simulink model using MATLAB and Simulink.Instead of programming C code manually, we can implement the control algorithm graphically using Simulink blocks.The models are saved as MDL files.In this research we built combustion model with stepwise method and optimized by PSO.All of this work was done using Simulink and MATLAB.In Figure 8, we can see the combustion model with PSO built by Simulink.
Simulation with the PSO was repeated 10 times based on the engine speed (rpm) from 1000 rpm to 2000 rpm, and 10 simulated engine optimal control parameters were obtained.The engine optimal control parameters are listed in Table 4.

Comparison of Validation Result and the Estimation
Result of an Engine Control Parameter.The engine optimal control parameters were validated in an engine test bench, and the results of validation based on the simulated optimal control parameters were compared with the calculated engine optimal objective values.The experiment result as the validation result and the calculation result are shown in Figures 9,10,11,and 12.The comparison of the experiment and the calculated engine optimal value of power, BSFC, NO  , and soot are illustrated in Figures 9-12, respectively.Based on Figures 9-12, it is observed that the calculated optimal value of power is similar to the validated value from engine test bench.The validation optimal value of BSFC is similar to the calculated optimal value except in 1600 rpm that increases 17.22%.The validated optimal value of NO  increases 13.77% in 1500 rpm and 13.59% in 2000 rpm but the validated optimal value of soot decrease 15.38% in 1300 rpm.Based on this result, the PSO method shows the effective methods in this optimization problem.Based on the observation in engine test bench, result of the optimization simulation, and the validation, conclusions for this study are the optimal control input parameters obtained by PSO were tested and analyzed.The result proved that the PSO is an effective method for engine optimization problem.

Conclusions
Based on the experiment data, in order to control the large number of control parameters appropriately in consideration of power, BSFC, NO  , and soot as the engine output objectives, the model construction which reproduces the characteristic value of power, BSFC, NO  , and soot from control parameter is needed.In this study, the stepwise method considering multicollinearity was applied to construct the polynomial model order 1 and order 2.
The accuracy of predictions made using stepwise method models considering multicollinearity depends on how well the regression function fits the data; there should be regular checks to see how well a regression function fits a given data set.This can be done through regular updates to ensure that the error values are always below a prespecified error threshold.
In this research, we have reported our predicting model of power, BSFC, NO  , and soot in multiple injection diesel engines by stepwise method considering multicollinearity.This paper shows that the predictive accuracy by stepwise method considering multicollinearity is high.This is proved by the multiple correlation coefficient 0.9000 or more and  value is very high.It can be regarded that the stepwise method can effectively estimate the objectives.
In order to improve exhaust emissions and fuel efficiency in a diesel engine, we use particle swarm optimization (PSO), one of optimization techniques to find the optimal engine operating condition efficiently.
Based on the observation in engine test bench, result of the optimization simulation, and the validation, conclusions for this research are the optimal control input parameters obtained by PSO were tested and analyzed.The result proved that the PSO is an effective method for engine optimization problem.

Figure 1 :
Figure 1: Schematic of diesel engine system.

Figure 7 :
Figure 7: Estimate value based on stepwise considering multicollinearity.

Figure 9 :Figure 10 :
Figure 9: Comparison of estimation and calculation of optimal value of power.

Figure 11 : 1 )Figure 12 :
Figure 11: Comparison of estimation and calculation of optimal value of NO  .

Table 1 :
Specification of diesel engine.

Table 2 :
Diesel engine control parameters.
The dotted line is experiment value and solid line is estimate value.The multiple correlation coefficient  is 0.98009 but this model uses 17 variables and  value is only 692.2417.The estimate value and actual measurement of power by stepwise method considering multicollinearity with 2 order polynomial equations are shown in Figure7.If compared with Figure6, the multiple correlation coefficient  is only Xj and the velocity Vj of each particle are initialized by random numbers Evaluate fitness value J = f(Xj) of each objective function Comparison with the old best value f(XPbest) of each particle Comparison with the best value f(XGbest) of whole particle Update the position Xj and velocity Vj by XGbest and XPbest min < F out variable of F min F value as F max maximum F value as F max the minimum F value as F min Select variable from P. Make a regression analysis calculation using Select variable from Q. Make a regression analysis calculation using Q Arrange the F value in descending order.Define P = group of unadded variable Q = group of added variable p = number of variable in P q = number of variable in Q Figure4: Model-building procedure of stepwise regression.

Table 4 :
The engine optimal control parameters.