Some previous research results have shown that EGR (exhaust gas recirculation) rate, pilot fuel quantity, and main injection timing closely associated with engine emissions and fuel consumption. In order to understand the combined effect of EGR rate, pilot fuel quantity, and main injection timing on the NOx (oxides of nitrogen), soot, and ISFC (indicated specific fuel consumption), in this study, CFD (computational fluid dynamics) simulation together with the Taguchi method and the ANOVA (analysis of variance) technique was applied as an effective research tool. At first, simulation model on combustion and emissions of a light duty diesel engine at original baseline condition was developed and the model was validated by test. At last, a confirmation experiment with the best combination of factors and levels was implemented. The study results indicated that EGR is the most influencing factor on NOx. In case of soot emission and ISFC, the greatest influence parameter is main injection timing. For all objectives, pilot fuel quantity is an insignificant factor. Furthermore, the engine with optimized combination reduces by at least 70% for NOx, 20% in soot formation, and 1% for ISFC, in contrast to original baseline engine.
1. Introduction
Because of robust design and high thermal efficiency, diesel engines increasingly and widely are applied to many power equipments, including light duty, medium duty, and heavy duty vehicles [1, 2]. However, under the pressure of environmental pollution and energy crisis, diesel engine performance is being continuously enhanced in order to pursue a higher power density, improve fuel economy, and reduce emissions characteristics to meet more stringent emissions regulations [1, 3–6].
There are many factors that will affect the combustion and emissions of diesel engines, due to the fact that the combustion process and emission formation are a process of complex physical and chemical action. Taguchi method is one of the efficient and cost-effective methods to parametric analysis and optimization performance of objectives with a remarkable decrease in number of trials and time involved [7–9]. Thus, the combination of test and Taguchi method is an effective way to investigate the effect of different factors in controlling the emissions and fuel economy [1–3, 10, 11]. However, the decreased experimental numbers still spend a lot of time and money.
On the other hand, current multidimensional CFD tools have become sufficiently mature to guide the design and development of more efficient and cleaner internal combustion engines [4]. Indeed, simulation models can provide insight into the fundamental processes occurring during combustion [12]. For example, Mobasheri et al. [13] successfully used a CFD simulation to explore the combined effects of pilot, post- and multiple-fuel injection strategies and EGR on engine performance and emission formation in a heavy duty DI-diesel engine. They found that injecting adequate fuel in postinjection at an appropriate EGR allows significant soot reduction without a NOx penalty rate.
In addition, numerous studies have shown that EGR rate, pilot fuel quantity, and main injection timing associated with engine NOx, soot emissions and fuel consumption [1–3, 10, 11, 13]. It is necessary to conduct an in-depth study in order to understand these three factors on the relationship with emissions and fuel consumption. Hence, in this study, CFD simulation together with the Taguchi method was used to investigate the combined effect of EGR rate, pilot fuel quantity, and main injection timing on the NOx, soot, and ISFC in a diesel engine at 2000 r/min, 50% load, with minimum number of experimental work. Moreover, for purpose of determining the effects of influencing factors on the objectives, the ANOVA technique was applied to analyze simulation results.
2. Computational Procedure, Model Validation, and Methodology2.1. Computational Procedure
The multidimensional CFD calculations were carried out with the commercial CFD package AVL FIRE. The turbulent combustion model is ECFM-3Z (extended coherent flame model-3 zones), which includes a description of the local mixture stratification by considering three mixing zones (see Figure 1): the unmixed fuel zone F, the mixed zone M containing fuel and air, and the unmixed air zone A. Each region is separated into two parts: unburned gases region and burned gases region [14, 15]. For turbulent combustion phenomena, the ECFM-3Z model uses a 2-step chemistry mechanism for the fuel conversion like [16]
(1)CαHβOγ+(α+β4-γ2)O2⟶αCO2+β2H2O(2)CαHβOγ+(α2-γ2)O2⟶αCO+β2H2
where α, β, and γ are the atomic number of C, H, and O elements, respectively.
Principle of the ECFM3Z model.
The k-ζ-f turbulent model is used for modeling turbulent flow in the combustion chamber. Hanjalić et al. developed this model. They proposed an eddy viscosity model based on the concept of Durbin’s elliptic relaxation, which solves a transport equation for the velocity scales ratio ζ=v2¯/k instead of v2¯, thus increasing the robustness of the model, reducing the sensitivity to nonuniform grids [17].
The extended Zeldovich mechanism is widely accepted as a successful method for predicting NO concentration. Chemical reactions are given by [6, 17]
(3)O+N2⟷κ1NO+N(4)N+O2⟷κ2NO+N(5)N+OH⟷κ3NO+H
where κ1, κ2, and κ3 are the rate of reaction; the unit is cm3/mol·s.
The Soot emission is calculated using the kinetic model, which contains a reduced number of species and reactions and has been developed in order to provide a computationally efficient kinetic overall soot model [16]. Table 1 presents several submodels used in present simulation model [13, 18].
Computational submodels.
Breakup model
KH-RT model
Evaporation model
Dukowicz model
Collision model
Nordin model
Wall interaction model
Bai Gosman model
Because of the symmetry of the combustion chamber and the injector at the center of the combustion chamber has 6 nozzle holes, the simulations were performed on one sixth of the combustion chamber for reducing calculation time. The mesh of the simulation model on the top dead center was shown in Figure 2.
Mesh of the simulation model on the TDC.
2.2. Model Validation
The engine used in this study was a light duty, 4 cylinders, four stroke, turbocharging, water-cooled diesel engine. Specifications of simulated engine were listed in Table 2. The experimental system was shown in Figure 3. The in-cylinder pressure was measured using a Kistler 6125B quartz crystal pressure sensor and an ONO SOKKI DS-9100 Combustion analyzer was used to analyze the in-cylinder pressure. The NOx emissions were measured with a HORIBA MEXA-7100D exhaust gas analyzer. Smoke was measured by smoke meter.
Engine specifications.
Bore (mm)
93
Stroke (mm)
102
Compression ratio
17.2
Swirl ratio
2
Number of nozzle holes
6
Number of cylinders
4
Displacement (L)
2.771
Intake valve closure (CAD)
595
Exhaust valve opening (CAD)
846
Schematic of experimental setup. (1) Dynamometer controller; (2) dynamometer; (3) charge amplifier; (4) PC control; (5) pressure sensor; (6) engine; (7) encoder; (8) intake cooler; (9) Compressor; (10) gas flow meter; (11) air cleaner; (12) exhaust gas analyzer; (13) smoke meter; (14) turbine; (15) EGR valve; (16) EGR cooler; (17) ECU; (18) PC control; (19) common rail; (20) fuel flow meter; (21) fuel tank.
The simulated original baseline condition of light duty diesel engine was actuated at 2000 r/min, 50% load with a pilot injection included in injection strategy. Table 3 lists initial conditions and input parameters for simulation.
Calculation conditions.
Intake temperature (IVC) (K)
330
Intake pressure (IVC) (MPa)
0.151
Engine speed (r/min)
2000
EGR rate
7%
Start of pilot injection (°BTDC)
19
Start of main injection (°ATDC)
1
Fuel injection quantity (mg/cycle)
25
TDC (top dead center) on the compression stroke was defined as 720 CAD (crank angle degree). To ensure the intake air mass amount was equal to the actual quantity of inlet air, the intake air temperature obtained after several adjustments was set to 330 K at IVC (intake valve closure). The well-justified assumptions were applied such that an entire piston surface, cylinder head, and cylinder wall are at constant temperature. Constant temperature of 500 K was used for the piston surface and cylinder head, and 400 K was used for the cylinder wall. From Figures 4 and 5, it can be seen that good agreement was achieved between experimental and simulation results. Although there are still some differences between simulation and test, Mobasheri et al. [13] considered that these variances related to testing uncertainties in input parameters for the simulations, for instance, the precise injection duration, start of injection timing, and gas temperature at intake valve closure.
Comparison of calculated and measured on cylinder pressure.
Comparison of calculated and measured on heat release rate.
2.3. Taguchi Method
Instead of having to test all possible combinations of design parameters, the Taguchi method tests the special combinations according to arrangements of orthogonal arrays (OA), in order to obtain the optimal combination and save time and expenses [9, 19, 20]. Another important feature of Taguchi method is the signal-to-noise (S/N) ratio employed as a measure of the impact of noise factors on performance. In this study, EGR rate, pilot fuel quantity, and main injection timing were considered as the influencing factors. In addition, NOx emission, soot emission, and ISFC are the response variables of the influencing factors. Therefore, the smaller-the-better characteristics of S/N ratio are used in the present study. S/N ratio is calculated as depicted in the following formula:
(6)S/N=-10log1n(∑i=1nyi2),
where yi is the value of the quality characteristic measured from the test and n is the number of experiments. The unit of S/N ratio is dB [2, 9, 19, 20].
In present study, a L9 orthogonal array for three factors and three levels was chosen and given in Table 4. The details of three levels of the chosen factors are given in Table 5. In Table 4 columns 1, 2, and 3 represent the EGR rate, the pilot fuel quantity, and main injection timing, respectively.
Taguchi’s L9 orthogonal arrays.
Trial number
Column 1 (EGR rate)
Column 2 (pilot fuel quantity)
Column 3(main injection timing)
1
1
1
1
2
1
2
2
3
1
3
3
4
2
1
2
5
2
2
3
6
2
3
1
7
3
1
3
8
3
2
1
9
3
3
2
Factors with chosen levels.
Factors
Levels of factors
1
2
3
EGR rate
7%
17%
27%
Pilot fuel quantity
7.2%
17.2%
27.2%
Main injection timing
4°CA BTDC
1°CA ATDC
6°CA ATDC
EGR rate used in this work was calculated by the concentration of carbon dioxide in intake and exhaust gas through the following equation:
(7)EGR=[CO2]inlet[CO2]exhaust×100%.
Definition of pilot fuel quantity was the proportion of an existing fuel quantity in pilot injection, which accounted for the total cycle fuel injection quantity.
2.4. ANOVA Technique
ANOVA is an analytical method that identifies the factor which shows a large dispersion by analyzing the dispersion of the characteristic value of each factor [20]. In this study, the data from CFD simulation results based on the orthogonal arrays can be analyzed by performing analysis of variance to illustrate the degree of importance of each factor that prominently influenced the response variables [11, 19, 20]. The corresponding mathematical statements established with ANOVA for present study are given in (8), (9), (10), and (11). Equation (8) expresses the mean value of the overall S/N, where k is the number of the trials. The sums of squares due to the variations of the overall mean (SS) is calculated by (9). The mean of the influencing factors (SSi) is expressed as (10). The percentage contribution of the individual factors on the chosen response variables will be considered to decide the optimized combination and it is calculated through (11) as follows [2, 9, 19, 20]:
(8)S/N¯=19∑k=19(S/N)k,(9)SS=∑i=19((S/N)ij-S/N¯)2,(10)SSi=∑j=13((S/N)ij-S/N¯)2,(11)Contribution%=SSiSS×100%.
3. Results and Discussion3.1. CFD Computational Results
After conducting the CFD simulation with the combination of influencing factors from L9 OA, the computational results of response variables were shown in Table 1. The results derived from the calculated mass formation divided by the calculated IMEP (indicated mean effective pressure), then they were normalized as discussed by Hajireza et al. [21], in respect to the corresponding value for original baseline case in order to see a comparable result.
3.2. Analysis of S/N Response Curves
Table 7 shows the S/N ratios for each factor computed in the Taguchi process and based on the results of Table 6. From the Table 7, the combination which has the maximum S/N ratio was considered as the best combination in achieving the optimum objective since as the assessment of quality characteristics, the larger S/N value, and the better results of response variables [11]. The S/N response curves are graphical expressions of change in response variables with the variation in influencing factor level. The Figures 6, 7, and 8 show the response curves of Taguchi experiment. In fact, S/N ratio in Figures 6, 7, and 8 is the value of corresponding (S/N)ij, which have considered the mean effect of each levels in each factor. U, V, and W in Tables 7 and 8 represent EGR rate, pilot fuel quantity, and main injection timing, respectively.
Calculated normalized results of response variables.
Trial number
NOx
Soot
ISFC
1
1.0876
0.5770
0.9822
2
0.6593
0.8876
1.0000
3
0.3268
1.8455
1.0251
4
0.5788
1.0485
1.0042
5
0.3426
1.8547
1.0356
6
0.5966
0.6473
0.9942
7
0.2403
2.2587
1.0592
8
0.2637
0.7622
0.9906
9
0.1129
1.5820
1.0079
S/N ratio results.
Factors
Levels
Trial number
S/N ratio
(S/N)ij
(i)
(j)
(k)
NOx
Soot
ISFC
NOx
Soot
ISFC
U
7%
1
-0.7297
4.7766
0.1561
4.2013
0.1633
-0.0199
2
3.6189
1.0356
0.0000
3
9.7148
-5.3225
-0.2157
17%
4
4.7500
-0.4114
-0.0363
6.1800
-0.6665
-0.0967
5
9.3041
-5.3657
-0.3040
6
4.4859
3.7776
0.0502
27%
7
12.3834
-7.0773
-0.4995
14.3029
-2.9007
-0.1617
8
11.5794
2.3591
0.0823
9
18.9458
-3.9839
-0.0680
V
7.20%
1
-0.7297
4.7766
0.1561
5.4679
-0.9040
-0.1266
4
4.7500
-0.4114
-0.0363
7
12.3834
-7.0773
-0.4995
17.20%
2
3.6189
1.0356
0.0000
8.1674
-0.6570
-0.0739
5
9.3041
-5.3657
-0.3040
8
11.5794
2.3591
0.0823
27.20%
3
9.7148
-5.3225
-0.2157
11.0488
-1.8429
-0.0778
6
4.4859
3.7776
0.0502
9
18.9458
-3.9839
-0.0680
W
4°CA BTDC
1
-0.7297
4.7766
0.1561
5.1118
3.6378
0.0962
6
4.4859
3.7776
0.0502
8
11.5794
2.3591
0.0823
1°CA ATDC
2
3.6189
1.0356
0.0000
9.1049
-1.1199
-0.0348
4
4.7500
-0.4114
-0.0363
9
18.9458
-3.9839
-0.0680
4°CA ATDC
3
9.7148
-5.3225
-0.2157
10.4674
-5.9218
-0.3397
5
9.3041
-5.3657
-0.3040
7
12.3834
-7.0773
-0.4995
Results of ANOVA.
Factors
S/N¯
SS
SSi
Contribution %
NOx
Soot
ISFC
NOx
Soot
ISFC
U
8.228
274.32
57.3127
5.0226
0.0101
20.9
3.23
2.88
V
-1.1347
155.39
15.5789
0.783
0.0017
5.7
0.5
0.49
W
-0.0928
0.3506
15.4943
45.693
0.1001
5.6
29.41
28.54
S/N response curves for NOx.
S/N response curves for soot.
S/N response curves for ISFC.
As it can be seen from Figure 6, the peak point of each factor occurs at the EGR rate of 27%, pilot fuel quantity of 27.2%, and main injection timing of 6°CA ATDC (after top dead center), respectively. It is obvious that the S/N ratio of factor EGR rate varies largely and three factors show a consistent trend in affecting the NOx varying with the level. The reason is that EGR has been proved to be a very effective NOx reduction technique and retarded main injection timing decreases in-cylinder peak heat release rate and peak temperature, which lead to a reduction in NOx [22]. As for pilot injection, it can shorten the ignition delay and raise the in-cylinder temperature before main combustion, thus avoiding the rapid heat release rate and increasing the gas temperature too quickly. Accordingly, the NOx emissions, which have a close relationship with temperature, can be reduced [23]. For further analysis and discussion we will draw on ANOVA technique in the following sections.
It is observed from Figure 7 that the optimum factors are at 7% EGR rate, 17.2% pilot fuel quantity and 4°CA BTDC (Before Top Dead Center) main injection timing. For pilot fuel quantity, the maximum S/N ratio in Figures 7 and 8 appears in corresponding point that at level 2. This illustrates that it is necessary to choose proper pilot fuel quantity for optimizing combustion process in-cylinder. Since the change in pilot fuel quantity changes the fuel concentration distribution, this may result in two different effects, either enhancing or deteriorating fuel-air mixing [24]. In addition, it is found that EGR rate and main injection timing have the same effect for soot emissions. In other words, with EGR rate increasing and main injection timing retarding, the soot emissions increased. It is also illustrated that EGR and main injection timing have a contrary influence for soot emissions as compared to NOx emissions. This effect is attributed to the introduced high level EGR in combustion chamber giving rise to the reduction in oxygen availability for fuel combustion, and the retarded main injection would weaken atomization and vaporization of the injected fuel; furthermore, those will lead to incomplete combustion and increase in soot emission [25].
Figure 8 shows plots of the S/N ratio of ISFC against EGR rate, pilot fuel quantity, and main injection timing. As is seen, the results show the same trend as Figure 7. This means that the soot emissions and ISFC can be reduced simultaneously by selecting the same factors and levels.
3.3. Application of ANOVA
In order to determine the effects of influencing factors (pilot fuel quantity, EGR rate, and main injection timing) on the NOx emissions, soot emissions, and ISFC of a light duty diesel engine, the ANOVA was carried out and the results were shown in Table 8. Note that the calculation results of Table 8 depended on the results in Table 7.
As shown in the Table 8, percentage contributions of EGR rate, pilot fuel quantity, and main injection timing were 20.9%, 5.7%, and 5.6%, respectively, for NOx emissions. In the case of soot emissions, the percentage contributions of EGR rate, pilot fuel quantity, and main injection timing were 3.23%, 0.5%, and 29.41% respectively. In addition, with regard to ISFC, percentage contributions of EGR rate, pilot fuel quantity, and main injection timing were 2.88%, 0.49%, and 28.54%, respectively.
The main objective of percentage contributions is to analyze engine influencing parameters which significantly affect the performance characteristics. The larger value of percentage contributions means that the factor has a greater influence on the response variable. It can be observed that EGR is the most influencing factor on NOx. When compared with EGR, the other influencing parameters had slight influence for NOx control, though retarded injection has been recognized as an effective way to reduce NOx emissions by retarding the combustion phase which reduces the peak temperature in-cylinder resulting in suppressing the molecules of oxygen and nitrogen combine. Temperature and oxygen concentration play a crucial role in the formation of oxides of nitrogen. On the other hand, EGR can reduce the peak flame temperature and oxygen concentration at the same time. Those factors can further illustrate that, as compared to other influencing parameters, EGR is a very effective control technique for decreasing NOx [1, 26]. In addition, the two parameters of pilot fuel quantity and main injection timing have approximately the same level of influence on the NOx emissions, as shown in Table 8.
Based on Table 8, it also can be seen that the greatest influence of the soot and ISFC is exhibited main injection timing, then EGR rate, and pilot fuel quantity. In view of the fact that the percentage contributions of pilot fuel quantity is almost close to zero, it can be considered that pilot fuel quantity is almost insignificant in case of soot and ISFC, compared to main injection timing and EGR. This is because the pilot fuel quantity is smaller than the injected fuel in main injection, so the emissions caused by the pilot injection can be ignored in most cases and the total emissions simply relied on the main combustion [27].
3.4. Optimum Combination
Based on the analysis of S/N response curves and results of ANOVA, it is concluded that in the case of NOx, EGR is the most influencing factor, the other influencing parameters have slight influence, and the best level is at level 3. However, main injection timing is the greatest effect in regard to soot and ISFC; the peak point S/N ratio of this factor occurs at 4°CA BTDC. On the other hand, pilot fuel quantity is an insignificant influencing factor for all response variables, especially in terms of soot and ISFC. Hence, the optimum combination of factors and levels for the chosen response variables are at 27% EGR rate, 27.2% pilot fuel quantity, and 4°CA BTDC main injection timing.
3.5. Confirmation Experiment
The last step in Taguchi optimization method is to conduct confirmation experiments with the optimum combination of influencing factors for validating the improvement in the response variable by comparing it with the original baseline condition. After conducting the confirmation experiment, optimized engine reduces by at least 70% for NOx and 20% in soot formation. As for ISFC, a slight reduction of 1% compared to original baseline engine is shown.
Figure 9 shows in-cylinder equivalence ratio and emission distributions with respect to optimized and original baseline case by simulations. As is seen, in the optimized case, after the spray impinges on the top of the piston bowl at end of the main injection, equivalence ratio distributions indicate a more homogeneous mixture in-cylinder than fuel concentration distributions in original baseline case. It is also illustrated that the fuel-air mixing was enhanced as a result of the application of optimum parameter combination. In the case of comparison of emission distributions between optimized and original case, as shown in Figure 9, it is found that significant concentrations of NOx and soot emissions are found in regions corresponding to emission formation zones in original baseline case. However, most of NOx and soot emissions are eliminated in the optimized case, because of improved mixture formation and combustion processes [28].
Comparison of equivalence ratio and emission distributions between optimized and original case.
4. Conclusions
In the present work, CFD simulation together with the Taguchi and the ANOVA technique was used to investigate the combined effect of EGR rate, pilot fuel quantity, and main injection timing on the NOx, soot, and ISFC in a diesel engine at 2000 r/min, 50% load, with minimum number of experimental work. The following conclusions are drawn from the investigation.
According to the analysis of S/N curves, it is concluded that the S/N ratio of factor EGR rate varies largely and three influencing factors show a consistent trend in affecting the NOx varying with the levels. In addition, EGR and main injection timing have a contrary influence for soot as compared to NOx. In the case of soot emissions and ISFC, they can be reduced simultaneously by choosing the same factors and levels for each other.
The ANOVA results illustrated that EGR is the most influencing factor on NOx. For soot emission and ISFC, the greatest influence factor is main injection timing. However, pilot fuel quantity does not have important effects on all objectives.
The optimum combination of factors and levels for the chosen response variables are at 27% EGR rate, 27.2% pilot fuel quantity, and 4°CA BTDC main injection timing. Confirmation experiment was carried out to validate the accuracy of the CFD simulation results. And the engine with optimized combination reduces by at least 70% for NOx and 20% in soot formation. As for ISFC, a slight reduction of 1% compared to original baseline engine was shown.
The combination of the CFD simulation, the Taguchi method, and the ANOVA technique is found to be a good method, which can significantly save time and cost, to find optimum combination for low NOx, soot, and ISFC in a light duty diesel engine.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgment
This work is supported by the National Natural Science Foundation of China (Grant no. 51176064).
SaravananS.NagarajanG.SampathS.Multi response optimization of NOx emission of a stationary diesel engine201089113235324010.1016/j.fuel.2010.05.0152-s2.0-77955473670WuH.WuZ.Combustion characteristics and optimal factors determination with Taguchi method for diesel engines port-injecting hydrogen201247141142010.1016/j.energy.2012.09.0272-s2.0-84868582128KrishnanC.VinodY.SudhindraC.Optimization of fuel injection parameters for meeting euro III exhaust emission norms on a heavy duty diesel engine using Taguchi technique20082008-01-1610ReitzR. D.Directions in internal combustion engine research201316011810.1016/j.combustflame.2012.11.0022-s2.0-84869062217ParkW.LeeJ. Y.MinK.YuJ.ParkS.ChoS. H.Prediction of real-time NO based on the in-cylinder pressure in Diesel engines201334230753082KumarS.ChauhanM. K.Numerical modeling of compression ignition engine: a review20131951753010.1016/j.rser.2012.11.0432-s2.0-84871723632RezaeiS. Z.ZhangF.XuH.GhafourianA.HerrerosJ. M.ShuaiS.Investigation of two-stage split-injection strategies for a Dieseline fuelled PPCI engine201310729930810.1016/j.fuel.2012.11.0482-s2.0-84878356912HsuQ.-C.DoA. T.Minimum porosity formation in pressure die casting by Taguchi method20132013992086510.1155/2013/920865ThipprakmasS.Application of Taguchi technique to investigation of geometry and position of V-ring indenter in fine-blanking process20103152496250010.1016/j.matdes.2009.11.0462-s2.0-76549121869DiezA.ZhaoH.Investigation of split injection in a single cylinder optical diesel engine20102010-01-0605SaravananS.NagarajanG.SampathS.Combined effect of injection timing, EGR and injection pressure in NOx control of a stationary diesel engine fuelled with crude rice bran oil methyl ester201310440941610.1016/j.fuel.2012.10.0382-s2.0-84870430503KomninosN. P.RakopoulosC. D.Modeling HCCI combustion of biofuels: a review20121631588161010.1016/j.rser.2011.11.0262-s2.0-84855780976MobasheriR.PengZ.MirsalimS. M.Analysis the effect of advanced injection strategies on engine performance and pollutant emissions in a heavy duty DI-diesel engine by CFD modeling2012331596910.1016/j.ijheatfluidflow.2011.10.0042-s2.0-84857233826ColinO.BenkenidaA.The 3-zones Extended Coherent Flame Model (ECFM3Z) for computing premixed/diffusion combustion200459659360910.2516/ogst:20040432-s2.0-13844275723TruffinK.ColinO.Auto-ignition model based on tabulated detailed kinetics and presumed temperature PDF—application to internal combustion engine controlled by thermal stratifications20115423-244885489410.1016/j.ijheatmasstransfer.2011.06.0442-s2.0-80052178561AVL LIST GmbH2013Graz, AustriaAVL LIST GmbHHanjalićK.PopovacM.HadžiabdićM.A robust near-wall elliptic-relaxation eddy-viscosity turbulence model for CFD20042561047105110.1016/j.ijheatfluidflow.2004.07.0052-s2.0-11344267319PengZ.LiuB.WangW.LuL.CFD investigation into diesel PCCI combustion with optimized fuel injection20114351753110.3390/en40305172-s2.0-79953200871ThipprakmasS.PhanitwongW.Process parameter design of spring-back and spring-go in V-bending process using Taguchi technique2011328-94430443610.1016/j.matdes.2011.03.0692-s2.0-79957838304LeeD. H.ParkJ. S.RyuM. R.ParkJ. H.Development of a highly efficient low-emission diesel engine-powered co-generation system and its optimization using Taguchi method201350149149510.1016/j.applthermaleng.2012.06.0402-s2.0-84866550023HajirezaS.RegnerG.ChristieA.EgertM.MittermaierH.Application of CFD modeling in combustion bowl assessment of diesel engines using DoE methodology20062006-01-3330SalehH. E.Effect of exhaust gas recirculation on diesel engine nitrogen oxide reduction operating with jojoba methyl ester200934102178218610.1016/j.renene.2009.03.0242-s2.0-66149164176ZhangJ.2011Birmingham, UKUniversity of BirminghamOkudeK.MoriK.ShiinoS.YamadaK.MatsumotoY.Effects of multiple injections on diesel emission and combustion characteristics20072007-01-4178GomaaM.AliminA. J.KamaruddinK. A.Trade-off between NOx, soot and EGR rates for an IDI diesel engine fuelled with JB5201111111987199310.3923/jas.2011.1987.19932-s2.0-79955563969AgarwalD.SinghS. K.AgarwalA. K.Effect of Exhaust Gas Recirculation (EGR) on performance, emissions, deposits and durability of a constant speed compression ignition engine20118882900290710.1016/j.apenergy.2011.01.0662-s2.0-79953066746SoidS. N.ZainalZ. A.Spray and combustion characterization for internal combustion engines using optical measuring techniques—a review201136272474110.1016/j.energy.2010.11.0222-s2.0-79551513971KimD.ParkS.Optimization of injection strategy to reduce fuel consumption for stoichiometric diesel combustion20129322923710.1016/j.fuel.2011.08.0672-s2.0-84855963020