Enhancement of combustion characteristics of VCR diesel engine by optimizing engine parameters

In recent years, engine emissions have been one of the important problems which are of great concern. Hence, there is a growing need to develop engines with reduced emission. In the present study, Variable Compression Ratio diesel engine model has been validated by comparing the simulation results with the experimental. The study is aimed at analyzing the effect of compression ratio, exhaust gas recirculation, fuel injection pressure and start of injection on engine performance and emission characteristics. Using composite desirability technique, the engine parameters have been optimized to achieve lower NOx, soot and ISFC. The optimum combination has been observed at Compression ratio 17.52, Start of injection −30.1 °aTDC, Fuel injection pressure 736.06 bar and Exhaust gas recirculation 28.29%. ISFC, NOx and soot are reduced by 2.37%, 29.11% and 83.81% respectively. Higher Target Fuel Distribution Index indicates the improved mixture homogeneity for the optimized parameters.


Introduction
Emission control has been a major challenge faced by the modern world. IC engine being a major source of emission have to be optimized in order to improve combustion and reduce emission. In addition, diesel engines have proved to be a great threat to the environment compared to SI engines in consideration to their high NOx and soot emission. Hence, diesel engines must be designed in order to improve their performance. These days the use of catalytic converters has come into existence. However, these are expensive. Designing an engine without incorporating any such components would be a cost effective approach [5, 16-18, 23, 26, 31].
Experimental work carried out by [1] explains the effect of EGR on performance, emission and carbon deposits in a diesel engine. It explains how the EGR displace oxygen in the intake air and reduces the oxygen content in the combustion chamber. The specific heat of the intake air mixture increases which reduces the flame temperature. This also reduces the specific fuel consumption. In a nutshell, EGR is an effective method that can be applied to a diesel engine to reduce NOx emission without sacrificing its fuel economy and efficiency. Min et al. [20] studied the effect of EGR by changing the N 2 mass fraction results show that the indicated specific NOx emissions were decreased with increased EGR fractions while increases soot emissions. The simultaneous reduction of NOx and soot was done by optimizing the engine parameters. Jayashankara and Ganesan [14] worked on diesel engine to study the effect of fuel injection timing and injection pressure on its performance. Advanced injection | https://doi.org/10.1007/s42452-021-04739-6 timing resulted in an increase in in-cylinder temperature, pressure, heat release rate and NOx emissions. Advancing the injection timing increased the soot while retarding the injection timing reduced the soot emission. But at high fuel injection pressures (1.7-2 bar) advancing the injection timing reduced the soot emission. According to the work of [28] the increase in fuel injection pressure reduced the specific fuel consumption and specific energy consumption whereas decrease in injection pressure increased fuel consumption per power output.
Song et al. [30] investigated on the effect of compression ratio on combustion characteristics for homogeneous charged compression ignition engine. It was observed that compression temperature and pressure in the cylinder vary significantly with compression ratio. The study concluded almost zero NOx emission at compression ratios 10.7 and 14. Prasad et al. [24] investigated that increasing compression ratio helps in improving efficiency and reduces soot where as NOx emissions increased. The effect of FIP (Fuel Injection Pressure) on the combustion and emission characteristics was not monotonic and there exists an optimum value of fuel injection pressure. Various studies have also been reported by using preheated alternative fuels blend with various oxide nano particles in the works of Elumalai [9,11]). It was found that inherent oxygen content improved its combustion characteristics.
From the literature it can be ascertained that improvement in the mixture homogeneity can be achieved by improvising the air-fuel mixing rate. The parameters such as EGR, compression ratio, SOI, injection pressure, etc. have shown a significant effect on air-fuel mixing rate. This can be done by optimizing the engine parameters such as fuel injection pressure, compression ratio, injection timing, EGR, etc. [4].
The primary objective of this work is to improve the mixture homogeneity and analyze the effects of various parameters and their interactions on the ISFC, NOx, and soot. The simulations were performed based on design of experiments (DOE) using response surface methodology. The work also aims to find the optimal combination of compression ratio, exhaust gas recirculation, fuel injection pressure, and start of injection, which would result in better performance with reduced NOx, soot emissions.

Computational domain and methodology
CONVERGE™ [27] has been used for modeling the axi-symmetric engine sector and for simulation purpose. The computational domain of the VCR engine is shown in Fig. 1. The fluid flow properties are governed by the following equations: Conservation of mass

Conservation of momentum
Conservation ofenergy Equation of state presented below is used to couple density, temperature and pressure Optimization of engine parameters has been done using Design of Experiments (DOE) method. This method is effective in optimizing the engine parameters by evaluating the combined effect of input factors on output responses. The DOE technique considered in the present study is Response Surface Methodology (RSM). RSM is an efficient way of reducing the computational time with reasonable accuracy. The main aim of RSM is optimizing an output response which is effected by various input factors [6,22]. The following 2nd order polynomial equation represents a simple approximating function where, 'n' denotes the number of variables, 0 denotes the constant term, i denotes the coefficient of linear terms, ij denotes the coefficient of interaction terms and ε denotes the error.

Validation of VCR engine model (Diesel)
The simulation and experimental results are compared to validate numerical results obtained for VCR diesel case from CFD code CONVERGE™. The VCR engine specifications are given in Table 1. Validation of the predicted numerical results with that of the experimental data was carried out by comparing the pressure variation with crank angle. Figure 2 shows the variation of pressure with crank angle for experimental and simulation results. The trends of simulation are similar to that of experiment. However, a slight difference of around 1.2% in peak pressure is observed. Experimental pressure shows lesser as compared to simulation. This may be due the fact that during the compression stroke the gases may escape from the crevice region and would end up with lower in-cylinder pressure. This shows that the present model is in good accordance with the experimental results of the VCR diesel engine.
Validation of emissions were also performed by comparing the experimental soot and NOx emissions with the simulation values and given in Table 2. It is observed that they are in good harmony.

Results and discussions
Box-Behnken method was applied for designing the numerical experiments. In this method each variable is maintained at 3 equal intervals. For example, if the CR range is from 12 to 18, it is assigned as 12, 15 and 18 and considered for numerical experiment design. Since there are 4 variables and 3 center points, it becomes a set of 27 experiments with each block having 2 variations (means two independent variables are varied per experiment). The design of experiments and their responses are given the Table 3.

Interaction effects of VCR Diesel engine
The interactive effect graphs are generated through mod-eFrontier [12] software. Figure 3 shows the interaction plots of EGR and CR at different levels of SOI and FIP VCR diesel engine. It can be observed from figure that in the interaction of EGR and CR, as the EGR is increased from 0  to 40% ISFC increased from the level 195 to 265 g/kWh. It is observed from the figure that the ISFC is favorable in the mid to high range of FIP and SOI. The operating range of CR between high and mid (fig (a), (b) and (d)) seems to be quite favorable in reducing ISFC. Further, a detailed view can reveal that operating range was good enough in the case of high FIP and mid SOI. Figure 4 represents the interaction of CR and FIP on NOx emissions at different levels of EGR and SOI. It is observed from the figure that the NOx emissions are favorable in mid to high range of EGR and low to mid-range of SOI. It can also be interpreted from the figure that the increase in the fuel injection pressure increases NOx emissions [19]. Low EGR and high SOI combination is not fruitful at any stage of the interaction of CR and FIP. The effect of fuel injection pressure and compression ratio on NOx emissions mainly depends on the state of the mixture and ignition delay [29]. As the EGR is increased, the NOx decreases. This is due to the fact that EGR increases the specific heat of intake mixture, thus reduces the combustion temperature. It may also lower the amount of oxygen in the intake mixtures.
This combination of reduced combustion temperatures and reduced oxygen content lowers the NOx emissions. Figure 5 shows contour plots of soot for SOI and CR at different levels of FIP and EGR. It is seen from the figure that the interaction of CR and SOI is dominant in the region of high EGR and mid to high FIP. The variation of interaction is less aggressive in the low EGR range. The interaction causes more fruitful region when the other two parameters (EGR and FIP) are set at low and high levels respectively [2,15].
It is observed from the Table 4 that the difference between predicted and adjusted R 2 values is less than 0.2 for all the three responses, which implies that the model was able to fit the data with reasonably good accuracy. Analysis of variance (ANOVA) is used to determine the significant and insignificant parameters and their interactions. ANOVA gives the p value for the different response parameters, such as ISFC, NOx, and soot emissions, and are given in Table 5. All the individual parameters (CR, SOI, FIP, and EGR) are significant for ISFC, NOx, and soot [13].

Optimization of VCR diesel
The objective of any engine design is to minimize emissions and maximize the performance. It can be observed from the above discussions that the tradeoff is observed between the NOx and soot emissions. So there is a need for optimization of the design parameters aimed to reduce both NOx and soot without losing its performance. Optimization was performed using the desirability approach which is search based optimization technique [6,7]. The whole approach is to modify all the responses (X i ) as a separate desirability function (D i ) and then will be varied in the entire range.
If the response is at its mark value or objective, then (Di) = 1 on the other hand, if the mark value is in the infeasible region then the (D i ) = 0 [21]. Each of the desirability function is established depending upon the goal of optimization. Table 6 shows the criteria of optimization in the desirability approach. where U i is the upper limit of the responses (X i ), T i is the target of the response. The optimum combination of the VCR diesel case was found out to be the 17.52 CR, 30.1°bTDC SOI, FIP 736.06 bar and EGR 28.29% with a composite desirability of 0.95.The compression ratio did not alter much from the baseline case which is also beneficial in terms of the design aspect.

Comparison of optimized and baseline configuration of VCR (Diesel)
The optimized and baseline cases of the VCR diesel engine were compared for pressure, temperature, NOx, soot, HC and CO vs crankangle and are shown in the Figs. 6, 7, 8, 9, 10 and 11 respectively. It was noticed from the Fig. 6 that the optimized case is superior in performance as the peak pressure of optimized is more as compared to the baseline case this is also manifested in Table 7 with reduced ISFC for optimized case. It is noticed from the Fig. 7 that the average temperature of the optimized case reduces during the combustion. This is because of the presence of EGR which could suppress the combustion temperature [8]. Figures 8 and 9 shows comparision of the NOx and soot formation for the optimised and baseline cases. It is evident that the both NOx and soot levels were deceased for the optimized case. It is also observed from Figs. 10 and 11 that the HC and CO emissions are decreased in the optimum case than the baseline. Table 7 gives the comparison of baseline and optimized cases of VCR diesel based on performance and emissions. It is noticed that the ISFC decreased by 2.37%

Mixture homogeneity comparison of the baseline and optimized cases
TFDI (Target Fuel Distribution Index) determines the homogeneity of the mixture [25] and is compared for the baseline and optimized cases at different crank angle instances during the combustion as shown in Fig. 12.The target range of the equivalence ratio for the practical diesel engine applications is from 0.3 to 1.2 [3]. It is observed that the optimized case has more TFDI than the baseline case in all the instances. This is a clear indication of the better homogeneous mixture formation in the optimized case than the baseline case. Better homogeneity was attained due to advanced injection, higher fuel injection pressure and presence of EGR in the optimized case. Better homogeneity favors the lesser soot formation in optimized case. Figure 13 shows the temperature distribution plotted in the target region. It is also observed that the target region has lesser in-cylinder temperatures for optimum configuration than the baseline case which is responsible for lesser NOx.

Conclusions
The VCR diesel engine was analyzed using DOE based study and optimization has also been done using desirability technique. Start of Injection is the second most influential parameter in view of ISFC. This implies that the CR and SOI are two potential parameters in deciding the engine performance. However, there is a limit in advancing SOI due to the possibility of negative work. Overall the individual effects play a major role in deciding the ISFC. It is crucial to mention that the optimal solutions are only for a particular load and speed other load and speed shall have different optimum conditions. Hence, in future numerical experiments can be carried out to obtain optimal conditions for various loads and speeds. The parametric study may be analyzed further and can be optimized using Genetic Algorithm based optimization techniques such as NSGA II, ARMOGA etc. to obtain Pareto optimal solutions.
Author contributions Prabhakara Rao Ganji: Collected the data from the literature survey, Conceived and design the analysis, Graphical work. Rajesh Khana Raju Vysyaraju: Supervisor 1, Contributed data analysis, Numerical analysis, Optimization. Srinivasa Rao Surapaneni: Supervisor 2, Reference work, Performed the result analysis, Wrote the paper. B Karuna Kumar: Reference work, Performed the result analysis, Proof read the paper.

Conflict of interest
The authors declare that they have no conflict of interests.

Ethical approval
The authors confirm that the final version of the manuscript has been reviewed, approved, and consented for publication by all author.
Ethical Standards This study was as per the ethical standards approved by the Animal Use Ethics Committee of the University of Oeste Paulista (UNOESTE), Presidente Prudente, São Paulo, Brazil, by Protocol Number 4396.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.