Modeling and Multi-Objective Optimization of NO x Conversion Efficiency and NH 3 Slip for a Diesel Engine

The objective of the study is to present the modeling and multi-objective optimization of NOx conversion efficiency and NH3 slip in the Selective Catalytic Reduction (SCR) catalytic converter for a diesel engine. A novel ensemble method based on a support vector machine (SVM) and genetic algorithm (GA) is proposed to establish the models for the prediction of upstream and downstream NOx emissions and NH3 slip. The data for modeling were collected from a steady-state diesel engine bench calibration test. After obtaining the two conflicting objective functions concerned in this study, the non-dominated sorting genetic algorithm (NSGA-II) was implemented to solve the multi-objective optimization problem of maximizing NOx conversion efficiency while minimizing NH3 slip under certain operating points. The optimized SVM models showed great accuracy for the estimation of actual outputs with the Root Mean Squared Error (RMSE) of upstream and downstream NOx emissions and NH3 slip being 44.01 ˆ 10 ́6, 21.87 ˆ 10 ́6 and 2.22 ˆ 10 ́6, respectively. The multi-objective optimization and subsequent decisions for optimal performance have also been presented.


Introduction
Diesel engines are widely used in almost all kinds of commercial vehicles and some passenger vehicles because of their good dynamic performance, fuel economy, and durability.However, one of the main challenges regarding the use of diesel engines is the reduction of NO x emissions, which is a common air pollutant and may cause photochemical smog.This is especially true against the background that there is a growing global concern over environment pollution and climate change issues resulting from the use of automobiles.Generally, the urea based selective catalytic reduction (urea-SCR) is widely adopted as a promising and efficient technique for decreasing the NO x emissions [1][2][3][4] through redox reactions between NH 3 and NO x .This technology was first used in stationary applications and has become pretty popular for vehicles over the last decade [5][6][7].
The main advantages of SCR technology include the fact that it has a high NO x conversion rate (up to 90% or more) and can improve fuel economy.Furthermore, the reliability and durability of the SCR system are also superb.However, this technique requires external supply of reductant (urea) and will generate additional NH 3 slip because the non-reactive NH 3 will be emitted into the environment.NH 3 is also a pollutant and may cause severe respiratory diseases.Therefore, the main challenge of the SCR is the control of urea injection amount to deal with the optimal trade-off problem between the NO x conversion efficiency and NH 3 slip [8].
Researchers have proposed various control strategies and methods for urea injection, including open-loop control [9], closed-loop control [10][11][12] and model predictive and adaptive control [13][14][15].Another way of dealing with this problem is the use of multi-objective genetic algorithms.The first step of the proposed method is to build engine and SCR models for the prediction of the upstream and downstream NO x emissions and NH 3 slip to obtain the two conflicting objective functions (NO x conversion efficiency and NH 3 slip) for the multi-objective optimization.The multi-objective optimization method is then used to optimize the decision variable (urea injection amount) to maximize NO x conversion efficiency while minimizing the NH 3 slip under certain operating points.
Figure 1 shows a schematic of a catalytic converter illustrating how NO x emissions are converted into nitrogen gas and water.Figure 1 also indicates the possible emissions of non-reactive NH 3 (NH 3 slip) and non-reactive NO x .
Sustainability 2016, 8, 478 2 of 13 challenge of the SCR is the control of urea injection amount to deal with the optimal trade-off problem between the NOx conversion efficiency and NH3 slip [8].
Researchers have proposed various control strategies and methods for urea injection, including open-loop control [9], closed-loop control [10][11][12] and model predictive and adaptive control [13][14][15].Another way of dealing with this problem is the use of multi-objective genetic algorithms.The first step of the proposed method is to build engine and SCR models for the prediction of the upstream and downstream NOx emissions and NH3 slip to obtain the two conflicting objective functions (NOx conversion efficiency and NH3 slip) for the multi-objective optimization.The multi-objective optimization method is then used to optimize the decision variable (urea injection amount) to maximize NOx conversion efficiency while minimizing the NH3 slip under certain operating points.
Figure 1 shows a schematic of a catalytic converter illustrating how NOx emissions are converted into nitrogen gas and water.Figure 1   A good number of research activities have been conducted and have proposed a variety of ways for predicting NOx emissions [16][17][18][19].A study [20] used artificial neural networks to build a prediction model for CO2, soot, and NOx.In references [19,21], an adaptive least squares support vector machine model was built for the prediction of NOx emissions with a novel update to tackle process variations.The support vector machine (SVM), which is established based on the structural risk minimization principle, has proven to exhibit better generalization performance than neural networks and other methods [22].In the work of Martinez-Moraleset [23], artificial neural networks were used for the non-linear identification of a gasoline engine for the purpose of evaluating objective functions used within an optimization framework involving the use of the NSGA-II genetic algorithm and the Multi-Objective Particle Swarm Optimization (MOPSO) algorithm.It was observed that the NSGA-II algorithm performed better than the MOPSO algorithm in the Pareto based optimization process.In the study [24], it was found, again, that the multi-objective optimization solution for the NSGA-II algorithm was better than that of the ε-constraint method.
Against the background given above, the current study proposed a novel method based on an SVM (Support Vector Machine) and GA (Genetic algorithm) for modeling a target diesel engine and SCR system to predict upstream and downstream NOx emissions and NH3 slip.The NSGA-II genetic algorithm was used for the multi-objective optimization involving the maximization of NOx conversion efficiency while minimizing NH3 slip because it has proven to give better optimal results compared with other optimization methods as observed in references [23,24].

Experimental Setup and Methods
As discussed in Section 1, the first step of the current study is to build prediction models for estimating engine and SCR system output as close as possible to the actual output y for an input vector , , , … , using optimal SVM (support vector machine) models.The SVM is the most appropriate method for obtaining the prediction model of a stochastic system from a small A good number of research activities have been conducted and have proposed a variety of ways for predicting NO x emissions [16][17][18][19].A study [20] used artificial neural networks to build a prediction model for CO 2 , soot, and NO x .In references [19,21], an adaptive least squares support vector machine model was built for the prediction of NO x emissions with a novel update to tackle process variations.The support vector machine (SVM), which is established based on the structural risk minimization principle, has proven to exhibit better generalization performance than neural networks and other methods [22].In the work of Martinez-Moraleset [23], artificial neural networks were used for the non-linear identification of a gasoline engine for the purpose of evaluating objective functions used within an optimization framework involving the use of the NSGA-II genetic algorithm and the Multi-Objective Particle Swarm Optimization (MOPSO) algorithm.It was observed that the NSGA-II algorithm performed better than the MOPSO algorithm in the Pareto based optimization process.In the study [24], it was found, again, that the multi-objective optimization solution for the NSGA-II algorithm was better than that of the ε-constraint method.
Against the background given above, the current study proposed a novel method based on an SVM (Support Vector Machine) and GA (Genetic algorithm) for modeling a target diesel engine and SCR system to predict upstream and downstream NO x emissions and NH 3 slip.The NSGA-II genetic algorithm was used for the multi-objective optimization involving the maximization of NO x conversion efficiency while minimizing NH 3 slip because it has proven to give better optimal results compared with other optimization methods as observed in references [23,24].

Experimental Setup and Methods
As discussed in Section 1, the first step of the current study is to build prediction models for estimating engine and SCR system output ŷ as close as possible to the actual output y for an input vector X " px 1 , x 2 , x 3 , . . ., x n q using optimal SVM (support vector machine) models.The SVM is the most appropriate method for obtaining the prediction model of a stochastic system from a small amount of experimental data [25].After obtaining the prediction models, non-dominated sorting genetic algorithm (NSGA-II) was implemented to deal with the trade-off between the two conflicting objective functions of NO x conversion efficiency and NH 3 slip under certain working points.The whole procedures of modeling and multi-objective optimization were implemented in MATLAB ® (MathWorks, Inc., Natick, MA, USA).

Data Collection
In this study, data was collected by running a four-stroke, six-cylinder YUCHAI YC6L-42 diesel engine on a dynamometer.The main specifications of the experimental engine as well as the dynamometer and associated instrumentation are presented in Tables 1 and 2 respectively.The Eddy-current dynamometer and associated instrumentation were used in the experiments to measure engine speed and torque with the fuel supply per cycle being obtained from the CAN bus directly.The AVL 4000 and the LDS6 instrumentation were used for the measuring of NO x emissions and NH 3 slip, respectively.In order to cover as many common operating points of the diesel engine as possible, 180 equally spaced operating points were selected, with speed ranging from 900 rpm to 2600 rpm (with a step size of 100 rpm) and engine load ranging from 10% to 100% (with a step size of 10%).Parts of these are shown in Table 3.

Data Preprocessing
Generally, there are some numerical discrepancies among input variables x 1 , x 2 , x 3 , . . ., x n , and, as a consequence, some variables may not have a significant influence on the system's output y.Therefore, in order to eliminate the negative impact caused by the huge numerical discrepancies of input variables, there was the need to normalize each variable to vary from 0 to 1.For each column of data sample x i , the normalized data sample is given by: where x i is the original data, x i 1 is the normalized data, and x i (max) and x i (min) are the maximum and minimum in x i , respectively.Subsequently, the normalized data x i 1 were used as the inputs for the SVM models.

Building and Optimizing the SVM Models
The Support Vector Machine (SVM) [26] was developed by Vapnik and Cortes in 1995.As a novel kind of machine learning method, the SVM is gaining increasing popularity because of its many attractive features and promising empirical performance.A detailed description of SVM theory can be found in references [27][28][29].
Here, only a brief description is given.A support vector machine takes advantage of the kernel function to map the input data onto a high-dimensional feature space.Linear regression is then performed in the high-dimensional feature space.As a result, non-linear problems can be addressed in a linear space through non-linear feature mapping.After training on the input data sample, the SVM model can be used to predict variables whose values are unknown.
The final prediction function used by an SVM is as follows: where a i is Lagrange multiplier, and x i is a feature vector corresponding to a training variable.The components of vector a i and the constant b are optimized during training.C is a penalty factor, which indicates the degree of attention paid to outliers and determines the range of a i .The larger the value of C is the more attention is paid to the outliers.The kernel function is represented by K px i , xq, and is one of the most important parts of the SVM model.There are four common kinds of kernel functions.Among them, the Gaussian radial basis function kernel is most commonly used because of its effectiveness and speed in the modeling process [30].
The Gaussian function takes the form: where g is the parameter of the kernel function, which is as important as penalty factor C, with x and x i representing independent variables.

Model Parameter Optimization with Grid Search and GA
Based on the preceding discussion, it is clear that the kernel function parameters C and g are considerably vital in the SVM model.The selection of parameters C and g directly influences the performance of the SVM model.Thus, it could be deduced that the performance of the SVM model lies in the choice of the kernel function parameter g and penalty factor C. In this study, a multi-algorithm combined method is proposed for optimizing the model parameters.
A Genetic Algorithm (GA) [31] is a kind of optimization algorithm, which uses genetics to simulate the natural evolution process.Three basic operators of selection, crossover and mutation are used to generate better offspring populations in order to find exact or approximate solutions to optimization problems.GAs have been successfully used to solve a broad spectrum of optimization problems owing to its high effectiveness and lower time consumption [32].
Therefore, in this study, a genetic algorithm was used to optimize the C and g parameters of the SVM model.Furthermore, the k-fold cross-validated root mean squared error (RMSE) [33] was selected as the fitness of the objective function for the GA to reflect the prediction accuracy of the SVM model.

RMSE "
where N is the number of samples, with y and ŷ representing the actual and predicted values respectively.The main objective of the GA was to yield the smallest k-fold cross-validated RMSE by searching for the best combinations of the C and g parameters for the SVM model.In a GA, prior to the initial random population generation, the ranges of each parameter need to be given, which are almost given empirically at present.If the ranges were too wide, the optimization process would be time-consuming.On the other hand, if the ranges were too small, the best parameters may not be captured in the ranges.Based on this, a grid search is proposed for determining the rough scope in this study.In addition, existing studies have shown that identifying a good parameter pair by searching exponentially in the sequences is more practical and less time-consuming.For example, C = 2 ´2, 2 ´1, ¨¨¨, 2 10 ; g = 2 ´1, 2 0 , ¨¨¨, 2 10 [34].
The procedure of optimizing the C and g parameters of the SVM with grid search and GA is shown in Figure 2 as follows: Sustainability 2016, 8, 478 5 of 13 Therefore, in this study, a genetic algorithm was used to optimize the C and g parameters of the SVM model.Furthermore, the k-fold cross-validated root mean squared error (RMSE) [33] was selected as the fitness of the objective function for the GA to reflect the prediction accuracy of the SVM model.
where N is the number of samples, with y and y representing the actual and predicted values respectively.The main objective of the GA was to yield the smallest k-fold cross-validated RMSE by searching for the best combinations of the C and g parameters for the SVM model.In a GA, prior to the initial random population generation, the ranges of each parameter need to be given, which are almost given empirically at present.If the ranges were too wide, the optimization process would be time-consuming.On the other hand, if the ranges were too small, the best parameters may not be captured in the ranges.Based on this, a grid search is proposed for determining the rough scope in this study.In addition, existing studies have shown that identifying a good parameter pair by searching exponentially in the sequences is more practical and less timeconsuming.For example, C = 2 −2 , 2 −1 , ⋯, 2 10 ; g = 2 −1 , 2 0 , ⋯, 2 10 [34].
The procedure of optimizing the C and g parameters of the SVM with grid search and GA is shown in Figure 2 as follows: Step 1: Involves the use of a grid search in the range of C ∈ [2 −10 , 2 15 ] and g ∈ [2 −5 , 2 10 ] according to the exponential sequences mentioned earlier.The k-fold cross-validated RMSE for all pairs of C and g parameters is evaluated for the SVM model in the process of finding the best pair.
Step 2: Based on the best pair of C and g parameters determined in Step 1, the ranges of each parameter in the GA are established.Subsequently, a two-dimensional random initial population, which is binary encoded, is created.Each dimension (chromosome) represents C and g, respectively.
Step 3: Involves the calculation of the fitness function of all populations given by the k-fold crossvalidated RMSE for the SVM model.
Step 4: Based on fitness level (k-fold cross-validated RMSE), a better offspring population is generated through the use of the three basic genetic operators of selection, crossover and mutation.In this study, the probabilities for the selection, crossover and mutation operators were set to 0.9, 0.8 and 0.05, respectively.
Step 5: In this step, Step 4 was repeated until the stopping criterion (100 generations for the current study) was satisfied.Step 1: Involves the use of a grid search in the range of C P [2 ´10 , 2 15 ] and g P [2 ´5, 2 10 ] according to the exponential sequences mentioned earlier.The k-fold cross-validated RMSE for all pairs of C and g parameters is evaluated for the SVM model in the process of finding the best pair.
Step 2: Based on the best pair of C and g parameters determined in Step 1, the ranges of each parameter in the GA are established.Subsequently, a two-dimensional random initial population, which is binary encoded, is created.Each dimension (chromosome) represents C and g, respectively.
Step 3: Involves the calculation of the fitness function of all populations given by the k-fold cross-validated RMSE for the SVM model.
Step 4: Based on fitness level (k-fold cross-validated RMSE), a better offspring population is generated through the use of the three basic genetic operators of selection, crossover and mutation.In this study, the probabilities for the selection, crossover and mutation operators were set to 0.9, 0.8 and 0.05, respectively.
Step 5: In this step, Step 4 was repeated until the stopping criterion (100 generations for the current study) was satisfied.
After these steps, the best pair of C and g was obtained, ending the optimization process for the SVM models.

Multi-Objective Optimization
After building and optimizing the three prediction models for the upstream and downstream NO x emissions as well as NH 3 slip, the ultimate goal of this study was to optimize the urea injection amount to maximize the NO x conversion efficiency while minimizing the NH 3 slip, a situation that characterizes a typical multi-objective optimization problem.
A multi-objective optimization problem is usually concerned with maximizing or minimizing a number of objective functions in the presence of certain inequality and equality constraints, as well as other constraints in the form of lower and upper bounds defining the decision variable space.Generally, multi-objective optimization problems do not have a single optimal solution but a Pareto optimal set for the conflict between the various goals, and, as a consequence, it is therefore virtually impossible for multiple objectives to simultaneously achieve optimal results.
The NSGA-II genetic algorithm was proposed by Kalyanmoy et al. [35] to help reduce the computational complexity based on a certain number of decision variables and a given population of solutions, preserve the elite members of a population of solutions and eliminate the need for a sharing parameter associated with other multi-objective evolutionary algorithms like the Pareto-archived evolutionary strategy (PAES) [36] and the strength-Pareto evolutionary algorithm (SPEA) [37].NSGA-II takes advantage of the non-dominated sorting and crowding distance so that the algorithm has the capacity to approximate the best Pareto frontier and ensure that the obtained Pareto optimal solution has a good spreading.The specific description of the NSGA-II can be found in references [35,38,39].
The basic procedure for executing NSGA-II for the current multi-objective problem is summarized in Figure 3 and the steps are as follows: Step 1: Involves fixing the parameters and the range of the decision variable (urea injection amount).In the current study, an initial solutions population of 100 and a maximum number of generations of 200 were used for the multi-objective optimization.The lower and upper bounds for the urea injection amount variable were 0 and 2000 mL/h, respectively.
Step 2: Based on the parameters and ranges determined in Step 1, a random initial population, which is binary encoded, was created.
Step 3: Involves the calculation of the ranks and crowding distance of the population, followed by the application of the tournament selection method for selecting the best solutions in a particular population for creating a mating pool for producing child solutions based on non-dominated sorting.
Step 4: Involves the generation of offspring and parent population through crossover and mutation operators to ensure that the diversity within different generations of solutions is preserved.The best individuals were selected as the new population from the population incorporating offspring and parent population based on non-dominated sorting.In this study, probabilities of 0.9 and 0.1 were used for the crossover and mutation operators for the multi-optimization framework, respectively.
Step 5: In this step, a repetition of Step 3 and Step 4 was carried out until the stopping criterion (200 generations for the current study) was satisfied.
After these steps, the optimal Pareto set of the decision variable (urea injection amount) was obtained.

Results for SVM Model Simulation
Based on the modeling and optimization methods discussed in Section 2, three prediction models were obtained.One was an engine model for predicting upstream NOx emissions.The other two were an SCR model for predicting downstream NOx emissions and another for NH3 slip under different operating points.Half (1/2) of the equally spaced data samples were selected as the training dataset for these optimal SVM models.After the models were trained, the rest 1/2 samples were applied to the model to test the generalization ability of the optimal models.The flow chart for designing the optimal SVM models is shown in Figure 4. Figure 5a shows the comparison between the actual outputs and model prediction for the upstream NOx emissions.Figure 5b also shows the quality of fit between the actual outputs and the model prediction.It is evident from Figure 5a that the curves for actual and predicted values almost coincide, and the errors between them are considerably small.It can also be seen from Figure 5b that all the points are evenly distributed along the line where the predicted values tracked the actual values fairly well, implying that the upstream NOx emissions were predicted with good accuracy for all samples.Similar statistics are shown graphically for the downstream NOx emissions and NH3 slip in Figures 6 and 7, which also show good prediction accuracy.

Results for SVM Model Simulation
Based on the modeling and optimization methods discussed in Section 2, three prediction models were obtained.One was an engine model for predicting upstream NO x emissions.The other two were an SCR model for predicting downstream NO x emissions and another for NH 3 slip under different operating points.Half (1/2) of the equally spaced data samples were selected as the training dataset for these optimal SVM models.After the models were trained, the rest 1/2 samples were applied to the model to test the generalization ability of the optimal models.The flow chart for designing the optimal SVM models is shown in Figure 4.

Results for SVM Model Simulation
Based on the modeling and optimization methods discussed in Section 2, three prediction models were obtained.One was an engine model for predicting upstream NOx emissions.The other two were an SCR model for predicting downstream NOx emissions and another for NH3 slip under different operating points.Half (1/2) of the equally spaced data samples were selected as the training dataset for these optimal SVM models.After the models were trained, the rest 1/2 samples were applied to the model to test the generalization ability of the optimal models.The flow chart for designing the optimal SVM models is shown in Figure 4. Figure 5a shows the comparison between the actual outputs and model prediction for the upstream NOx emissions.Figure 5b also shows the quality of fit between the actual outputs and the model prediction.It is evident from Figure 5a that the curves for actual and predicted values almost coincide, and the errors between them are considerably small.It can also be seen from Figure 5b that all the points are evenly distributed along the line where the predicted values tracked the actual values fairly well, implying that the upstream NOx emissions were predicted with good accuracy for all samples.Similar statistics are shown graphically for the downstream NOx emissions and NH3 slip in Figures 6 and 7, which also show good prediction accuracy.Figure 5a shows the comparison between the actual outputs and model prediction for the upstream NO x emissions.Figure 5b also shows the quality of fit between the actual outputs and the model prediction.It is evident from Figure 5a that the curves for actual and predicted values almost coincide, and the errors between them are considerably small.It can also be seen from Figure 5b that all the points are evenly distributed along the line where the predicted values tracked the actual values fairly well, implying that the upstream NO x emissions were predicted with good accuracy for all samples.Similar statistics are shown graphically for the downstream NO x emissions and NH 3 slip in Figures 6 and 7 which also show good prediction accuracy.The RMSE (root mean square error), MAPE (mean absolute percent error) and R 2 (correlation coefficient) values as well as the fits of the three models are given in Table 4.The RMSE (root mean square error), MAPE (mean absolute percent error) and R 2 (correlation coefficient) values as well as the fits of the three models are given in Table 4.The RMSE (root mean square error), MAPE (mean absolute percent error) and R 2 (correlation coefficient) values as well as the fits of the three models are given in Table 4.The RMSE (root mean square error), MAPE (mean absolute percent error) and R 2 (correlation coefficient) values as well as the fits of the three models are given in Table 4.The MAPE for the upstream and downstream NO x emission prediction models were under 5% for all datasets except for the test dataset for the downstream NO x emissions, which were considered good enough for estimating the outputs.It was meaningless to calculate the MAPE for the NH 3 slip for the reason that the actual value of NH 3 slip may be 0 for certain particular working points, and the MAPE may not be defined mathematically (that is a constant divided by 0 cannot be defined mathematically).For example, it is difficult to achieve catalytic reduction when the temperature of SCR catalyst is under 200 ˝C.Under such a condition the urea injection amount is 0 and there is no NH 3 slip.The RMSE for three models of upstream NO x emissions, downstream NO x emissions and NH 3 slip were 44.01 ˆ10 ´6, 21.87 ˆ10 ´6 and 2.22 ˆ10 ´6 respectively, which were considerably small, while the squared correlation coefficients for the three models were 0.99, 0.99 and 0.98 respectively.This implied that the accuracies of the three models were good enough for the prediction of the upstream and downstream NO x emissions and NH 3 slip.The quality of fit between actual values and prediction values for the three models are also shown in Table 4.
It can be also seen from Table 4 that the RMSE and MAPE for the test datasets are a little larger than for the training datasets.The main reasons for this are as follows: Firstly, it was impossible to cover all the operating points of the diesel engine when choosing the training datasets for the model, which had a negative impact on the prediction accuracy.Secondly, some subjective factors existed in the parameter selection and data processing.Thirdly, there could be some outliers from such a large experimental data.Furthermore, there was no absolute guarantee that the experimental data would to be obtained in a completely steady state, and could lead to prediction errors as a consequence.

NSGA-II Multi-Objective Optimization Results
The prediction models for upstream NO x emissions, downstream NO x emissions and NH 3 slip are coded M1, M2 and M3, respectively, and are shown together with the inputs and outputs of the three models in Table 5.
Table 5.The inputs and outputs of the three prediction models.Equations ( 5) and ( 6) were used as objective functions for the optimization framework using the NSGA-II genetic algorithm, with Equations ( 7) and ( 8) as inequality constraints.The aim of the current research was to maximize Yε (NO x conversion efficiency) while minimizing the Y NH3 (NH 3 slip) using the NSGA-II genetic algorithm for the multi-objective optimization procedure discussed in Section 3.
To observe how the whole process functions, mode 4 in Table 3 was taken as an example, and the optimization was run and terminated after 200 iterations.Under certain operating speed, torque, fuel supply, temperature and upstream NO x emissions were fixed.The only decision variable that could be changed by the researcher to deal with the trade-off between Yε (NO x conversion efficiency) and Y NH3 (NH 3 slip) was the urea injection amount.Figure 8 shows the final 100 Pareto optimal solutions for the objective functions Yε (NO x conversion efficiency) and Y NH3 (NH 3 slip), and their corresponding decision variables (urea injection amount) used for the multi-objective optimization.It could be seen from the graph that it is not possible to improve (increase) NO x conversion efficiency without worsening (increasing) NH 3 slip.This trend generally characterizes a conflicting multi-objective optimization problem.
Among the 100 Pareto optimal solutions, in the initial stages there was a rapid increase in the NO x conversion efficiency with the increase of the urea injection amount, while the value for NH 3 slip increased slightly.On the other hand, the NO x conversion efficiency gradually decreased, while the rate of NH 3 slip gradually increased.Finally, NH 3 slip goes up rapidly with the increase of urea injection amount, while the NO x conversion efficiency remains constant at 0.95.
For this particular multi-objective optimization problem, related emissions legislation (WHTC for Euro VI) requires that the real-time value of NH 3 slip should be less than 10 ppm.Therefore, a vertical line was drawn through the 10 ppm point on the horizontal axis, with the two points A(9.7, 0.88) and B(9.7, 940) which are closest to the left of vertical line (NH 3 slip less than 10 ppm) on the curves for the Pareto optimal frontier and the urea injection amount were taken as the optimal solution for the multi-objective problem under operating mode 4.
Point A represents the final optimal solution for the multi-objective optimization problem and point B is the corresponding urea injection amount under operating mode 4. In other words, for mode 4, with an upstream NO x concentration of 684 ppm and a catalytic converter temperature of 250 ˝C, the optimal amount of urea which results in the optimal performance of the SCR catalytic converter (88% reduction of NO x emissions and 9.7 ppm NH 3 slip) is 940 mL/h.Similar results for other working conditions can also be achieved using the same method.The optimization results for modes 4-9 (engine speed = 1500 rpm and 2000 rmp) of Table 3 are summarized and presented in Table 6.The designer can also obtain other optimal solutions to the multi-objective optimization problem from the Pareto front by setting other restrictions in the form of constraints.
To observe how the whole process functions, mode 4 in Table 3 was taken as an example, and the optimization was run and terminated after 200 iterations.Under certain operating speed, torque, fuel supply, temperature and upstream NOx emissions were fixed.The only decision variable that could be changed by the researcher to deal with the trade-off between Yε (NOx conversion efficiency) and YNH3 (NH3 slip) was the urea injection amount.Figure 8 shows the final 100 Pareto optimal solutions for the objective functions Yε (NOx conversion efficiency) and YNH3 (NH3 slip), and their corresponding decision variables (urea injection amount) used for the multi-objective optimization.It could be seen from the graph that it is not possible to improve (increase) NOx conversion efficiency without worsening (increasing) NH3 slip.This trend generally characterizes a conflicting multiobjective optimization problem.
Among the 100 Pareto optimal solutions, in the initial stages there was a rapid increase in the NOx conversion efficiency with the increase of the urea injection amount, while the value for NH3 slip increased slightly.On the other hand, the NOx conversion efficiency gradually decreased, while the rate of NH3 slip gradually increased.Finally, NH3 slip goes up rapidly with the increase of urea injection amount, while the NOx conversion efficiency remains constant at 0.95.
For this particular multi-objective optimization problem, related emissions legislation (WHTC for Euro VI) requires that the real-time value of NH3 slip should be less than 10 ppm.Therefore, a vertical line was drawn through the 10 ppm point on the horizontal axis, with the two points A(9.7, 0.88) and B(9.7, 940) which are closest to the left of vertical line (NH3 slip less than 10 ppm) on the curves for the Pareto optimal frontier and the urea injection amount were taken as the optimal solution for the multi-objective problem under operating mode 4.
Point A represents the final optimal solution for the multi-objective optimization problem and point B is the corresponding urea injection amount under operating mode 4. In other words, for mode 4, with an upstream NOx concentration of 684 ppm and a catalytic converter temperature of 250 °C, the optimal amount of urea which results in the optimal performance of the SCR catalytic converter (88% reduction of NOx emissions and 9.7 ppm NH3 slip) is 940 mL/h.Similar results for other working conditions can also be achieved using the same method.The optimization results for modes 4-9 (engine speed = 1500 rpm and 2000 rmp) of Table 3 are summarized and presented in Table 6.The designer can also obtain other optimal solutions to the multi-objective optimization problem from the Pareto front by setting other restrictions in the form of constraints.

Conclusions
A novel ensemble methodology with GA and SVM has been presented for establishing the models for the prediction of upstream and downstream NO x emissions and the NH 3 slip.The various aspects of the modeling and optimization processes have been discussed in detail.The NSGA-II was used to solve the multi-objective optimization problem of maximizing NO x conversion efficiency while minimizing NH 3 slip based on the decision variable of urea injection amount under certain operating points.Based on the current study, the following conclusions could be drawn: (1) The prediction accuracy of the engine and SCR models could be improved by using an SVM, the parameters of which were optimized using a GA.The RMSE of upstream and downstream NO x emissions and NH 3 slip for the all datasets was 44.01 ˆ10 ´6, 21.87 ˆ10 ´6 and 2.22 ˆ10 ´6, respectively.The MAPE of the models were all under 5%, and was good enough for estimating the actual outputs.(2) The optimized urea injection amounts under certain operating points were obtained through multi-objective genetic algorithm for maximizing NO x conversion efficiency while minimizing NH 3 slip.
also indicates the possible emissions of non-reactive NH3 (NH3 slip) and non-reactive NOx.

Figure 1 .
Figure 1.Schematic for the conversion of NOx to nitrogen and water.

Figure 1 .
Figure 1.Schematic for the conversion of NO x to nitrogen and water.

Figure 3 .
Figure 3. General procedure for executing the Non-dominated Sorting Genetic Algorithm (NSGA-II) multi-objective optimization in the current study.

Figure 4 .
Figure 4. Processes for designing the proposed SVM models.

Figure 3 .
Figure 3. General procedure for executing the Non-dominated Sorting Genetic Algorithm (NSGA-II) multi-objective optimization in the current study.

Figure 3 .
Figure 3. General procedure for executing the Non-dominated Sorting Genetic Algorithm (NSGA-II) multi-objective optimization in the current study.

Figure 4 .
Figure 4. Processes for designing the proposed SVM models.

Figure 4 .
Figure 4. Processes for designing the proposed SVM models.

Figure 5 .
Figure 5.Comparison between experimental and predicted values for upstream NOx emissions.(a) Actual versus predicted output; (b) Quality of fit.

Figure 6 .Figure 7 .
Figure 6.Comparison between experimental and predicted values for downstream NOx emissions.(a) Actual versus predicted output; (b) Quality of fit.

Figure 5 .
Figure 5.Comparison between experimental and predicted values for upstream NO x emissions.(a) Actual versus predicted output; (b) Quality of fit.

Figure 5 .Figure 6 .Figure 7 .
Figure 5.Comparison between experimental and predicted values for upstream NOx emissions.(a) Actual versus predicted output; (b) Quality of fit.

Figure 6 .
Figure 6.Comparison between experimental and predicted values for downstream NO x emissions.(a) Actual versus predicted output; (b) Quality of fit.

Figure 5 .Figure 6 .Figure 7 .
Figure 5.Comparison between experimental and predicted values for upstream NOx emissions.(a) Actual versus predicted output; (b) Quality of fit.

Figure 7 .
Figure 7.Comparison between experimental and predicted values for NH 3 slip.(a) Actual versus predicted output; (b) Quality of fit.
NO x emissions, urea injection amount Downstream NO x emissions M3 Temperature, Upstream NO x emissions, urea injection amount NH 3 slip

Figure 8 .
Figure 8. Pareto-based optimal solutions to the multi-objective optimization problem.

Figure 8 .
Figure 8. Pareto-based optimal solutions to the multi-objective optimization problem.

Table 1 .
Technical specifications of the experimental engine.

Table 2 .
The parameter measurement methods.

Table 3 .
Engine operating conditions.

Table 4 .
The fit statistics for the three prediction models.

Table 6 .
Optimization results for the engine speeds of 1500 rpm and 2000 rpm.