Similar articles

Abstract: In this paper, a multi-objective particle swarm optimization (MOPSO) algorithm and a nondominated sorting genetic algorithm II (NSGA-II) are used to optimize the operating parameters of a 1.6 L, spark ignition (SI) gasoline engine. The aim of this optimization is to reduce engine emissions in terms of carbon monoxide (CO), hydrocarbons (HC), and nitrogen oxides (NO _x ), which are the causes of diverse environmental problems such as air pollution and global warming. Stationary engine tests were performed for data generation, covering 60 operating conditions. Artificial neural networks (ANNs) were used to predict exhaust emissions, whose inputs were from six engine operating parameters, and the outputs were three resulting exhaust emissions. The outputs of ANNs were used to evaluate objective functions within the optimization algorithms: NSGA-II and MOPSO. Then a decision-making process was conducted, using a fuzzy method to select a Pareto solution with which the best emission reductions can be achieved. The NSGA-II algorithm achieved reductions of at least 9.84%, 82.44%, and 13.78% for CO, HC, and NO _x , respectively. With a MOPSO algorithm the reached reductions were at least 13.68%, 83.80%, and 7.67% for CO, HC, and NO _x , respectively.

1.  Introduction

 Emissions of air pollutants are known to have a variety of negative effects on public health and natural environment. For several years, emissions of nitrogen oxides (NO _x ), hydrocarbons (HC), carbon monoxide (CO), and particulate matter (PM) have been regulated in different regions of the globe. It has been identified that hydrocarbons are among the most relevant pollutants as they are the major contributor to smog. They are the result of partially burned fossil fuels occurring during the operation of internal combustion engines; this is also known as incomplete combustion. One of the adverse effects to human health of this pollutant is the reduction of the blood’s ability to carry oxygen, which can lead to fatal results in cases of overexposure. Other emissions of interest are NO _x , which are generated when the nitrogen contained in the air reacts with oxygen at high temperature and pressure inside the engine. While NO _x are precursors to smog and acid rain, they are also part of greenhouse gases. This means that they trap the heat from solar radiation within our atmosphere provoking a rise in the temperature of the Earth’s surface. Thus, engine emissions play a significant role in global warming (Mori, 1997).

 In Europe, as in other regions, European emissions standards are created to define the acceptable limits for exhaust emissions for new vehicles sold in the European Union (EU). The compliance is determined by running the engine through a standardized test cycle. However, the new standards do not apply to vehicles already in use. The engines of these vehicles should be checked and recalibrated to meet current or planned standards.

 Spark ignition (SI) engines have multiple relevant inputs, such as the ignition advance timing, injected fuel mass flow, injection time, engine speed, and in addition air/fuel ratio, and intake manifold absolute pressure (IMAP), among many others that are feedforward or feedback controlled. The outputs of the engine (emissions, drivability) are influenced in a nonlinear way (Guerrier and Cawsey, 2004). Therefore, in the production of a new car, the design of control systems is a very time-consuming task. This requires recycling those engines that have been utilized in the past to power-up new vehicles. An interesting engineering paradox therefore takes place: in order to decrease the design time an old engine is used, but it is necessary to invest a lot of time in modifying the old engine in order to meet the new emission regulations.

 In order to optimize the engine operating parameters, a parameter search method should be utilized for many of the engine operating conditions to determine the ideal exhaust emissions. However, when this parameter search method is executed experimentally, it requires a great deal of resources, time, and labor. The optimization of engine parameters by simulation on a computer-based environment is very useful and is replacing the traditional tune-up method. This process is based on modeling the engine responses at a limited set of operating points obtained from experimental measurements using an engine test bench. Artificial neural network (ANN) models have been widely used in several engineering fields due to their ability to find highly nonlinear relationships. This ability makes ANNs appropriate to solve complex nonlinear problems. Many studies of engine calibration have investigated the application of ANNs to model emissions in combustion engines necessary to minimize engine responses (Langouet et al., 2011; Yap et al., 2012). Ismail et al. (2012) reported on an ANN model to predict nine different diesel engine responses including CO and carbon dioxide (CO₂). Ghobadian et al. (2009) used a multilayer perceptron network (MLP) to model the brake power, output torque, specific fuel consumption, and exhaust emissions of a diesel engine. Canakci et al. (2009) used an ANN model to predict the performance and exhaust emissions of a diesel engine including CO, NO _x and unburned hydrocarbons. However, Saerens et al. (2009) used a mean value model of an internal combustion engine instead of ANN to minimize the fuel consumption by acting on the throttle valve.

 Many researchers are conducting tests for the optimization of engine parameters with different types of optimization methods. Particle swarm optimization (PSO) (Kennedy and Heberhart, 1995) is one of the newest techniques within the family of evolutionary optimization algorithms. PSO is a swarm intelligence method that models social behavior to guide swarms of particles towards the most promising regions of the search space and has been extended in order to solve multiple objective optimization problems (Coello and Lechuga, 2002). A multiple objective particle swarm optimization (MOPSO) algorithm is efficient at solving a broad variety of engineering optimization problems due to its fast convergence (Zhao and Cao, 2005; Abido, 2009). On the other hand, genetic algorithms have been used to solve multiobjective optimization problems (MOPs) including minimization of exhaust emissions and fuel consumption. Langouet et al. (2011) used a local linear model tree (LOLIMOT) and multiobjective covariance matrix adaptation evolution strategy (MO-CMA-ES) algorithm to carry out the optimization problem of minimizing emissions of NO _x , HC, CO, CO₂, also for a diesel engine, utilizing six engine parameters for many engine operating conditions. D’Errico et al. (2011) used the nondominated sorting genetic algorithm II (NSGA-II) in order to find parameters that minimize the break specific fuel consumption and NO _x emissions, as well as maximize the torque of a single cylinder SI engine. Kesgin (2004) increased the efficiency keeping NO _x emissions under a constraint value of a turbocharged natural gas engine. Shi and Reitz (2010) minimized CO, unburned hydrocarbon (UHC), NO _x emissions of a heavy-duty compression ignition engine fueled with diesel and gasoline by using NSGA-II to find optimal combinations of eight optimization parameters. Atashkari et al. (2007) preferred NSGA-II to run a Pareto-based optimization of a variable valve-timing engine considering conflicting objectives such as fuel consumption and torque. Alonso et al. (2007) combined ANNs with genetic algorithms to optimize the diesel engine settings reaching important reductions in emissions as CO, HC, NO _x , break specific fuel consumption (BSFC), and PM for two engine operating conditions.

 Although the minimization of engine emissions via computerized optimization processes is not a novel proposal, only a few studies have been reported for gasoline engines (Atashkari et al., 2007; Saerens et al., 2009), while most of the work reported is for diesel engines. In this research project an ascendant approach is performed. In the first instance, ANN models are obtained for each exhaust emission of interest, i.e., CO, HC, and NO _x . The inputs for each model are the engine parameters found in the electronic control unit (ECU) of the engine, which can be modified to change the operating point, and eventually to minimize the corresponding emissions. Then, a minimization process is run using two multi-objective optimization algorithms (MOPSO and NSGA-II). A set of optimal solutions, called Pareto front, is found. The aim of this paper is to show the comparison of emissions reductions obtained by two evolutionary algorithms, and then to make an informed decision on the new values of the engine parameters. A final solution is obtained from the Pareto front by applying a fuzzy method.

 It is of a key importance to note that this work intends to contribute in three particular areas. In the first place, having a model that allows predicting the behavior of the internal combustion engine of a vehicle in terms of identified significant pollutants becomes substantial for shortening the time to market for a new engine. This study will show that it has the necessary accuracy to facilitate this conclusion. In the second place, a model that can be run to make a new set point of circulating vehicles in order to make them fulfill the strictest environmental regulations is necessary. It is estimated that testing old vehicles takes a lot of resources while keeping a high level of incertitude in the results. Finally, the adaptation process of a vehicle to a local market can take into account local environmental conditions that affect engine performance, which becomes critical when trying to keep emissions low without compromising performance.

2.  Gasoline engine emissions models

 Using real data obtained from an engine test bench, three neural network models are generated in order to predict engine emissions.

2.1.  Experimental setup

 A four cylinder SI gasoline engine coupled to a hydraulic dynamometer controlling the engine speed in real time is used to conduct the experiments. Exhaust emissions are measured by a gas analyzer type FGA4000XDS from InfraredIndustries® (USA), for 60 different values of the engine angular speed (n _eng) in revolutions per minute and the angle of the admission throttle valve (α _th) in %. There are four additional operating parameters whose values are determined for each one of 60 operating conditions. The additional operating parameters are the time of the open state of an injector (t _inj) in ms measured by a Tektronix® oscilloscope from USA. The injected fuel mass flow (m _fuel) in lb/h (1 lb=0.454 kg) is measured by the test bench sensor and conditioned by the software of the SF-902 hydraulic dynamometer from SuperFlow® Corporation in USA. Moreover, the ignition advance timing (θ _ign) in degrees (°) relative to #1 cylinder and the IMAP in kPa are taken by an on-board diagnostics II (OBD II) system.

 The ECU controls all actuators and regulators of the engine according to the measured signals obtained from the sensors. There are prerecorded integrated look up tables in memory for each operating condition that impose the values of the engine operating parameters. These look up tables, called cartography, can be adjusted in order to change the engine response in performance as well as in exhaust emissions. Tables 1 and 2 show the characteristics of the studied internal combustion engine and the hydraulic dynamometer, respectively.







 Table 3 shows the taken values for a total of 60 operating points (four angles of the admission throttle valve for 15 engine speed values), which are evaluated in the experimental setup to collect the values of CO, HC, and NO _x emissions. These 60 engine conditions are selected covering approximately the New European Driving Cycle (NEDC).




 In Fig. 1, the equipment to measure engine operating parameters and the engine test-bench used in this study is shown. It consists of a dynamometer connected with the crankshaft of the engine. This dynamometer has an electrically controlled throttle actuator to vary the angle of the throttle valve. A fuel canister is used to measure the fuel consumption of the engine in the range from 20 to 270 lb/h.




 The ignition advance timing is obtained with a crankshaft position sensor. This sensor indicates that when a notch on the crankshaft goes by and sends a pulse to the ECU, as it is referenced to the top dead point of the piston 1, it is possible to calculate the angular position of the crankshaft. A manifold absolute pressure (MAP) sensor provides manifold pressure information to the engine ECU. This MAP sensor is also used with an OBD-II PC. A tachometer (magnetic pickup on a 60-tooth gear) measures the engine speed and it is transformed to r/min. Fig. 2 shows the schematic diagram of the experimental apparatus. Measurements are performed with the engine operating at steady state; this data are recorded 10 s after setting an operating point to avoid the transient phase.




 Before carrying out the first measurement the engine is warmed for 10 min to ensure that the engine operates in an optimal temperature condition to establish a reference scenario. Table 4 shows some sample data of the studied gasoline engine for four operating points defined at an engine speed of n _eng=2500 r/min.

2.2.  Neural models of the engine exhaust emissions

 The LOLIMOT approximates a nonlinear function with piece-wise linear models (Hafner et al., 2000; Nelles, 2001). The predicted emission by LOLIMOT neural network y _pred is calculated by adding the contributions of all M linear models (1) , where x =[x ₁, x ₂, ···, x_Q ]^T are the Q inputs of the network, each neuron represents a local linear model (LLM) and an associated validity function that determines the validity region of the LLM. The outputs of each LLM are given by (2) , where p_ij and p_{i 0} denote the LLM parameters of the ith hidden neuron. The validity function of the ith hidden neuron corresponds to standard Gaussian function defined as (3) , where (4) , where the centers ρ_ij and standard deviations σ_ij define the position of the center of the peak and the width of the Gaussian functions, respectively. LOLIMOT is an incremental construction algorithm that divides the input space in hyper-rectangles. The validity function of the corresponding LLM is located in their center. The standard deviations σ_ij are proportional to the size of the hyper-rectangle.

 A new LLM is attached to a model partition during every iteration, which is derived from the LLM with the largest local error until the number of iterations is exceeded. The whole model is a network with one hidden layer (Nelles, 2001), and each neuron in this layer produces a linear part of the model. The output layer calculates the weighted sum of the outputs of local linear models. The approximating function generated by the neural network is used to transform the nonlinear behavior of the combustion engine inputs x into the outputs y _pred which approximates the actual output of the engine for each operating condition.

 Fig. 3 shows the structure of the used LOLIMOT neural network, whose engine input variables are the engine speed, the angle of the admission throttle valve, the injection time, the injected fuel mass flow, the ignition advance timing, and the IMAP. The output is one of the selected exhaust emissions. Due to the nonlinear relationship of the data, one ANN is constructed for each exhaust emission.

2.3.  Predicted emissions and performance criteria

 In this study, a total of 60 operating conditions are under consideration. Forty of them are used for training of the ANN model, while the other 20 cases are used for the validation and comparison of the ANN model. It is relevant to mention that a total of nine parameters are monitored during the process. Each operating condition is represented by a dataset. With the aim of conserving the statistical validity of the experiment, it was decided that the datasets for each phase, training as well as testing, would be randomly chosen with a uniform distribution. Thus, all data have the same probability of being selected, ensuring a different selection for each phase. Figs. 4 and 5 present the input and the output training and test datasets for each operating parameter and engine emission, respectively.







 Three different criteria are considered in order to evaluate the ANN performance: (1) the mean absolute percentage error (MAPE) that expresses accuracy as a percentage; (2) the correlation coefficient R ² used to evaluate the linear relationship between predicted and measured output; and (3) the root mean square errors (RMSE) that are taken into account in order to estimate the ANNs prediction qualities for many architectures of the models. MAPE, RMSE, and R ² performances are calculated as (Canakci et al., 2009) (5) , (6) , (7) , where is a measured value and a predicted value by the neural model over all test datasets.

 To choose the size of ANNs, MAPE, RMSE and R ² are calculated for the testing data for different numbers of hidden neurons as shown in Fig. 6 and Table 5. The RMSE and MAPE are smaller and R ² is close to one, with (6, 10, 1), (6, 9, 1), and (6, 15, 1) LOLIMOT neural network architectures for CO, HC, and NO _x prediction, respectively. The corresponding correlations between recorded and predicted exhaust emissions for these architectures are shown in Figs. 7–9. The obtained values in each model are very close to the experimental data. The values of R ² are 0.982, 0.972, and 0.998 for testing datasets in the LOLIMOT neural networks models in predicting CO, HC, and NO _x , respectively.

3.  Multi-objective optimization problem

3.1.  Definition

 The optimization approach described here uses ANNs as the basis for the evaluation of the objective functions. These are used in the optimization process. Mathematically, the multi-objective optimization problem can be written as (8) , where Ω is the parameter space, x is a point in this space, N is the number of objective functions and ζ is the vector-valued constraint function with upper and lower limits ζ _U and ζ _L. The three engine responses considered for minimizing are CO, HC, and NO _x emissions of the studied gasoline engine. This minimization is done by the variables of injection time t _inj, the injected fuel mass flow m _fuel, the ignition advance timing θ _ign, and the IMAP. The optimization process is performed for each operating condition defined by engine speed n _eng and the angle of the admission throttle valve α _th. The complete details of the multi-objective optimization problem are specified in Table 6.

3.2.  Objective functions

 The definition of the objective function is very important for solving an optimization problem, considering that the main objective of the optimization is to reduce CO, HC, and NO _x emissions, subject to the restrictions of the operating parameters of the studied gasoline engine. The objective functions are defined by (9) , where e( x ) _i represents the value of the ith emission for the operating parameters x , ME _i represents the maximum value allowed for the ith emission, and γ_i is a penalty factor. Here i=1, 2, 3, for CO, HC, and NO _x objective function f( x ) _i respectively. The larger the value of γ_i , the greater the value of the objective function, causing the algorithm to explore other regions in the search space. However, the values of γ_i are set to 1 for the multi-objective optimization problems treated here by using MOPSO and NSGA-II methods.

3.3.  MOPSO optimization method

 A MOP requires the simultaneous satisfaction of a number of different objective functions, under a number of inequality or equality constraints. MOPs generally give a possibly uncountable set of solutions, whose evaluated vectors represent the best possible trade-offs in the objective function space. In this sense, Pareto optimality is the procedure to establish a hierarchy among the solutions in order to determine whether a solution is really one of the best possible trades-off.

 The solution x ^* is said to be Pareto optimal (minimal) if no other solution can be found to dominate x ^* using the definition of Pareto dominance. Given two solutions x ^* and x , vector F ( x ^*) is said to dominate vector F ( x ) (denoted by ) if and only if (10) , (11) .

 In MOPSO, a population of potential solutions randomly distributed called particles is created over the search space with random velocity. The positions of the particles are moving towards the positions of the local best element X _j ^* (t) and the global best element X _j ^** (t). These are selected from the updated set of recorded nondominated solutions. Then, the particles change their positions during generations until a termination criterion is met (Coello and Lechuga, 2002).

 The following gives a brief introduction to the operation of the proposed MOPSO algorithm. Consider a swarm of n particles, and each particle is an m-dimensional vector, where m is the number of optimized parameters.

1. Initialize the number of iterations until a given limit is reached.

2. Generate randomly n particles, X _j , with j=1, 2, ···, n, where X _j =[x_{j ,1}, x_{j ,2}, ···, x_j,m ].

3. Generate randomly initial velocities of all particles V _j , with j=1, 2, ···, n, where V _j =[v_{j ,1}, v_{j ,2}, ···, v_{j ,m} ].

4. All initial particles are evaluated by using the objective functions.

5. Define the local best element for each particle as X _j ^* (0)= X _j (0), j=1, 2, ···, n.

6. Define a global set S ^** (0), with all nondominated solutions.

7. Calculate the distances between X _j (0) and the members in S ^**(0) by using the Euclidian distance. Define the global best element X _i ^** (t) as the nearest member in S ^**(t) to the ith particle.

8. Set the initial value of the inertia weight w(0).

9. Update the inertia weight w(t) according to (12) , where iter_max is the number of generations, and iter is the current number of iteration.

10. The jth particle velocity in the kth dimension is updated by (13) , where c ₁ and c ₂ are positive constants, j=1, 2, ···, n, k=1, 2, ···, m, and r ₁ and r ₂ are uniformly distributed random numbers in [0, 1].

11. Based on the updated velocities, each particle changes its position according to (14) .

12. The updated position of each particle X _i (t+1) is added to the local set S _i ^* (t+1).

13. Update the global set S ^**(t+1) with all the nondominated solutions of .

14. The nearest member in S _j ^* (t+1) to the jth particle is defined as X _j ^* (t+1). Similarly, the nearest member in the global set S ^**(t+1) to the ith particle is defined as the global best element X _j ^** (t+1).

15. If the number of iterations exceeds its limit then stop, else go to step 9.

3.4.  NSGA-II optimization method

 The NSGA-II method has a diversity preservation mechanism that involves calculation of the so-called crowding distance, which evaluates the density of individuals surrounding a particular individual in the population. Details on NSGA-II can be found in (Deb et al., 2002). It can be summarized as follows:

1. Generate a random population P ₀ of l individuals.

2. Evaluate objective functions for each individual, and sort the population using nondomination order and their crowding distances.

3. Initialize current population P _l = P ₀.

4. Apply binary tournament in order to create a mating pool of parents of size l/2 from P _l based on Pareto front ranking and crowding distance.

5. Randomly select couples from mating pool and apply bimodal crossover and polynomial mutation operators in order to generate offspring population R _l called children of size l.

6. Parents and children populations are merged, thus creat a single intermediate population made of 2l individuals; sort them using nondomination order relation.

7. The population is ranked according to the fitness parameters and only the best l individuals are chosen to create the new parents population for the subsequent generation P _{l +1} of size l.

8. Go to step 4 with P _l = P _{l +1}.

 In this study, the minimization process finds the optimal solution for each operating condition, defined by engine speed and the angle of the admission throttle valve. For each engine condition the operating parameters are constrained by maximum and minimum values defined by their cartography nominal levels.

3.5.  Fuzzy decision-making

 The multi-objective evolutionary algorithms provide a set of nondominated solutions which form the Pareto front; here a fuzzy-based mechanism is used to find an acceptable solution on the Pareto front. A membership value for the ith objective of the jth solution in the Pareto front is calculated using the membership function given by (15) where f_i is the value of the ith objective function, f_i ^max and f_i ^min are the maximum and minimum values of the ith objective function, respectively, and mf_i ^j indicates how well the jth nondominated solution is able to satisfy the ith objective. The accomplishment of each Pareto solution with respect to all the D nondominated solutions can be obtained as (16) .

 The solution with the maximum value of mf ^j is the solution that can be accepted by the decision-maker.

4.  Results

 For the MOPSO, the inertia weight factor w decreases linearly during the optimization run, w _max and w _min are set to 1.2 and 0.1, respectively, the amount of particles n is set at 50, iter_max=10, and constants c ₁=c ₂=0.5. Whereas for NSGA-II a population of l=80 individuals and G=20 generations, with crossover probability p _c=0.9 and mutation probability p _m=0.033 are considered here. The simulation programs are written in MATLAB and run in a 2.20 GHz Pentium-4 processor with 1 GB of RAM. To evaluate the convergence of the algorithm, a map of the objective space is created. The Pareto solutions indicating the selection made by the fuzzy decision method are presented in Figs. 10a–10c. For the optimization of the engine operating condition defined by n _eng= 3500 r/min and α _th=75%, reductions of 27.56%, 99.67%, and 83.72% for CO, HC, and NO _x emissions with NSGA-II are achieved, respectively. Utilizing MOPSO, reductions of 26.6%, 99.96%, and 69.97% for CO, HC, and NO _x emissions are achieved, respectively.




 After the NSGA-II and MOPSO converged, the values for the emissions corresponding to the obtained Pareto solutions and fuzzy decision-maker (for 10 engine operating conditions) are analyzed as shown in Tables 7 and 8, respectively. As a result of the optimization processes and the reduction percentages, the values of emissions were obtained with the nominal values of the operating parameters (Table 7 and 8). Table 7 shows the reduction for the operating point defined by n _eng=2900 r/min and α _th=25%. This is about 91.44% with NSGA-II for CO emissions, which is greater than that achieved with MOPSO that reduces emissions only 64.17%. Also for the engine operating condition n _eng=3100 r/min and α _th=75%, a reduction with NSGA-II (72.19% for CO) is achieved, which is higher than that with MOPSO which reduces 42.93% CO emissions. However, with the MOPSO algorithm, greater reductions (41.27%) for NO _x at the operating point defined by n _eng=3400 and α _th=75% are achieved, whereas the NSGA-II reduces emissions by only 13.78%. Similarly for the engine operating condition n _eng= 3750 r/min and α _th=75%, a reduction of 41.02% for NO _x emission with MOPSO is reached, while NSGA-II only reduces this emission by 26.64%.







 While n _eng=3750 r/min and α _th=75%, the reductions achieved are 96.00% and 97.33% for HC by employing NSGA-II and MOPSO, respectively.

 Tables 7 and 8 present reductions of at least 9.84%, 82.44%, and 13.78% achieved for CO, HC, and NO _x emissions respectively with NSGA-II, while with MOPSO algorithm reductions of at least 13.68%, 83.8%, and 7.67% for CO, HC and NOx emissions are reached, respectively. On the other hand, similar reductions are obtained for HC emissions with both algorithms in most engine operating conditions, for example reductions of 88.247% and 87.580% using NSGA-II and MOPSO respectively for the operating point n _eng=2000 r/min and α _th=50% are obtained.

 It is possible to use the NSGA-II to diminish emissions in some engine operating conditions and the MOPSO method in the rest of the engine conditions, such that the engine is calibrated with parameter values that result in a greater reduction of emissions. Tables 9 and 10 show the optimal configurations of parameters obtained with NSGA-II and MOPSO, respectively, for the 10 operating points of the studied gasoline engine which resulted in the identified emissions reductions.







 Thus, the emissions of the gasoline engine can be reduced by updating the engine operating parameters. As a comparison, Alonso et al. (2007) found reductions of 25.8%, 36.2%, and 4.4% for CO, HC, and NO _x emissions; however, this was for a diesel engine at full load operating condition, using the NSGA-II algorithm. Although the case study is not the same, similar reductions are obtained here by using NSGA-II algorithm.

5.  Conclusions

 In this paper, neural network models for exhaust emissions of an SI gasoline engine are used as the basis for the evaluation of the objective functions which are evaluated within the MOPSO and NSGA-II optimization methods. These methods are applied to find a solution to the problem of reducing the emissions of a gasoline engine. LOLIMOT neural networks are utilized to determine suitable predictions for three exhaust emissions under stationary conditions. The carbon monoxide (CO), hydrocarbons (HC) and nitrogen oxides (NO _x ) emissions are modeled and subsequently reduced. Performance of the ANNs has been evaluated by calculating the coefficient of correlation (R ²) and mean absolute percentage error (MAPE). LOLIMOT networks with architectures (6, 10, 1), (6, 9, 1), and (6, 15, 1) are found to be capable to approximate the CO, HC, and NO _x emissions of the studied gasoline engine, with R ² of 0.982, 0.972, and 0.998, respectively. The MOPSO and NSGA-II methods are performed for 10 engine operating conditions, where the limits in the engine operating parameters are established according to the values in the nominal engine cartography. Important improvements with respect to the nominal manufacturer engine cartographies are observed, reaching reductions of at least 9.84%, 82.44%, and 13.78% for CO, HC, and NO _x emissions respectively with NSGA-II, although with MOPSO algorithm reductions are reached at least 13.68%, 83.8%, and 7.67% for CO, HC, and NO _x emissions, respectively. The main limitation of this study is that the determined parameters are only valid for the particular engine used in this case study. However, the methodology can be used to perform the minimization of emissions on hydrogen or biodiesel fueled internal combustion engines.