The effect of injection timing and intake valve close timing on performance and emissions of diesel PCCI engine with a full engine cycle CFD simulation
Ming Jia,Maozhao Xie,Tianyou Wang,Zhijun Peng
DOI: https://doi.org/10.1016/j.apenergy.2011.03.024
IF: 11.2
2011-01-01
Applied Energy
Abstract:A full-cycle computational fluid dynamics (CFD) simulation coupled with detailed chemical kinetics mechanism has been used to investigate the effect of start of injection (SOI) timing and intake valve close (IVC) timing on performance and emissions of diesel premixed charge compression ignition (PCCI) engine. By sweeping SOI timing from −35 to −5 °CA ATDC and IVC timing from −140 to −80 °CA ATDC with fixed 50% exhaust gas recirculation (EGR) and 1.8 bar intake pressure, the contour plots for ignition timing, nitric oxides (NO x ), soot, hydrocarbon (HC), carbon monoxide (CO), indicated specific fuel consumption (ISFC), and ringing intensity have been developed. The results indicate that the operating range can be divided into kinetically controlled region and mixing-controlled region, in which the ignition timing is solely controlled by IVC timing and SOI timing respectively. To Minimize HC, CO, NO x and soot emissions, SOI timing must be carefully adjusted within a limited range. With the retarded IVC timing, the operating range of SOI becomes wider for clean combustion. The IVC timing should be optimized with consideration of ignition timing and combustion efficiency at different SOI timing in order to improve fuel economy. For purpose of avoiding engine knock, the SOI timing around −20 °CA ATDC and early IVC timing are pursued. Keywords Premixed charge compression ignition (PCCI) Computational fluid dynamics (CFD) Variable valve actuation (VVA) Intake valve closing (IVC) 1 Introduction Premixed charge compression ignition (PCCI) combustion for diesel engines is receiving increasing attentions in recent years, because it is capable of providing both high efficiency and extremely low emissions of nitric oxides (NO x ) and particulate matter (PM). Unlike conventional diesel combustion, PCCI combustion mode uses early in-cylinder fuel injection and high exhaust gas recirculation (EGR) ratio to avoid NO x and PM emissions simultaneously. However, due to the low volatility and high ignitability of diesel fuel, several technical barriers, such as mixture preparation and control of combustion phasing, must be overcome before the commercialization of diesel PCCI engine [1] . In PCCI combustion mode, since the fuel is injected well before top dead center (TDC), obtaining good mixing in the available time and preventing wall-wetting due to spray over-penetration can be challenging [2] . A large number of injection strategies have been proposed for diesel PCCI combustion in recent years. The concept of narrow angle direct-injection (NADI) was suggested by Walter and Gatellier [3,4] to keep the fuel target within the combustion bowl and avoid the interaction between the spray the cylinder liner at advanced injection timing. Lechner et al. [5] evaluated the potential of using the injector nozzle with narrow spray cone angle for diesel PCCI combustion by varying the spray cone angle and advanced injection timing. It is found that carbon monoxide (CO) emission rises dramatically and combustion efficiency drops significantly as the fuel starts to miss the bowl. The Mie-scattering technique was used to investigate in-cylinder combustion and spray evolution processes by Fang et al. in a diesel PCCI engine with a 70-degree spray angle injector [6] . The results indicate that liquid fuel impingement on the bowl wall leads to fuel film combustion which is called “pool fire”. Because of rich air–fuel mixture and low temperature on the wall surface, the pool fire results in incomplete combustion and high soot formation for all early injection cases. Siewert [7] studied the effects of varying nozzle spray angle and rail pressure on emissions and thermal efficiency in a diesel PCCI engine with early injection timing. It is indicated that the amount of liquid spray that misses the piston bowl is directly linked to the measured increases in hydrocarbon (HC), CO and smoke emissions and reduction in thermal efficiency. It should be noted that injection timing does not provide an effective means of controlling ignition timing as in the conventional diesel engine, so that controlling the combustion phasing is still a critical issue for PCCI engine [1] . As one of the most promising strategies for controlling diesel PCCI combustion, variable valve actuation (VVA) attracts increasing attentions recently because of its fast response time. By using late intake valve closing (IVC), the effective compression ratio is decreased and advanced combustion can be avoided, thus it is helpful for the increase of thermal efficiency and reduction of NO x emissions. Helmantel and Denbratt [8] tested several valve timing strategies for controlling diesel PCCI combustion with a Lotus adjustable valve timing system. It is found that lower fuel consumption is achieved by late IVC due to the reduction of pumping losses, however, VVA cannot be used as a replacement for EGR due to the limit of knock intensity. Nevin et al. [9] incorporated a hydraulically actuated variable IVC system into a single cylinder diesel engine to explore the potential of intake valve actuation on heavy-duty diesel engine. It is shown that late IVC results in a 90% NO x reduction at high speed and intermediate load with constant airflow and no EGR. With the help of late IVC, less EGR is required to achieve the same NO x emissions, while CO and PM could be dramatically reduced by nearly 70%. He et al. [10] experimentally investigated the feasibility of using late IVC to improve the performance and emissions of diesel engine with a fully flexible valve actuation system. They found that late IVC reduces approximately 25–50% NO x emissions and more than 95% PM at some operating conditions. Walter et al. [11] evaluated the potential of VVA for diesel homogeneous charge compression ignition (HCCI) engine in a single cylinder engine with several intake and exhaust VVA configurations. It is found that the reduction of compression ratio by late IVC is not beneficial for the expansion of near-zero NO x operating range due to the decrease of total mass trapped inside the cylinder and limited equivalent ratio at a given intake pressure. Recently, Kim et al. [12] indicated that late IVC is helpful for the reduction in PM and NO x emissions without increasing fuel consumption, but it can also expand the applicable operating range of the PCCI engine. With the help of computational modeling, deeper insights into the influence of late IVC on diesel PCCI could be obtained. Genzale et al. [13] used a multi-dimensional computational fluid dynamics (CFD) model to explore the effectiveness of late IVC in controlling diesel PCCI combustion phasing over various loads, engine speeds and boost pressures. It is indicated that the application of late IVC could optimize combustion phasing and increase output work, thus allowing the possibility to employ PCCI operation at higher loads. A more detailed study based on equivalence ratio–temperature map and three-dimensional (3D) CFD model was conducted by Murata et al. [14–16] . They found that the use of late IVC combined with a large amount of EGR simultaneously reduces NO x and soot emissions while maintaining the fuel consumption performance, since local maximum combustion temperature is lowered and local over-rich regions is avoided during the combustion processes. All previous studies show the potential of late IVC for future diesel PCCI engine with the reduction of pollutant emissions and/or the improvement of fuel economy. However, the impact of late IVC combined with advanced injection timing on fuel–air mixing and combustion processes is still not fully understood. Although multi-dimensional models have been used previously [13–15] , the gas exchange processes were not considered in their simulations and the initial conditions were simplified by the assumption that the in-cylinder pressure and temperature conditions at IVC are the same as those in the intake manifold regardless of what time the intake valve closes. This simplification could lead to underestimation of the effective compression ratio and neglect of the in-cylinder flow induced by the valve motion. Since the in-cylinder tumble ratio, swirl ratio and turbulence intensity are significantly affected by the IVC timing, which will further influence the mixing and combustion processes of diesel PCCI engines, it is necessary to take account of the intake processes in the simulations. In this paper, a full engine cycle simulation was performed by using a 3D CFD model coupled with detailed chemical kinetics to understand the impact of late IVC and injection timing on fuel–air mixing and combustion characteristics in a diesel PCCI engine. 2 Test engine The engine simulated in this study is derived from a Ford four-cylinder, high-speed direct-injection (HSDI) diesel engine with a displacement of 0.5 liters/cylinder. A major modification to the original engine is decreasing the compression ratio from 18.2 to 16.0 in order to realize PCCI combustion. A new piston with an open-crater-type combustion bowl and a six-hole nozzle with narrow spray included angle (120°) were used in this study to give better mixture preparation. The geometry of the combustion bowl was optimized by Boyarski and Reitz [17] for diesel PCCI combustion with early fuel injection and narrow spray cone angle. The engine specifications are summarized in Table 1 . This study focuses on the part-load operating conditions with the torque of around 106.3 Nm and the engine speed of 1600 RPM. Our previous study indicates that boosting intake pressure is extremely important for the reduction of soot, HC and CO emissions for the strategies with late IVC due to the improved homogeneity of the fuel–air mixture and increased oxygen masses [18] . Meanwhile, sufficient EGR is needed to realize the ultra-low NO x combustion, even when late IVC timing is employed in the system. The corresponding operating conditions tested in this study are shown in Table 2 . 3 Numerical model The CFD code used in the present study is KIVA-3V [19] . To gain better predictions, various physical and chemistry submodels in KIVA-3V have been updated as listed in Table 3 . The turbulent flows were modeled using the RNG k – ε turbulence model [20] , modified for variable-density engine flows. Spray is represented by a particle method, several submodels are included in the code to account for spray atomization, breakup, collision, vaporization and spray/wall interaction [19] . Kelvin–Helmholtz Rayleigh–Taylor (KH-RT) model [21] was used to simulate the breakup process of injected droplets. The spray collision model developed by Nordin [22] with the improvement of grid independence was integrated into the KIVA-3V code, and the spray/wall interaction model used here is developed by Han et al. [23] with the consideration of the gas density in computing the size of secondary droplets in the splash processes. Heat transfer from the wall was simulated by the Han and Reitz model [24] which accounts for the variations in the gas density and the turbulent Prandtl number in the boundary layer. Diesel fuel model DF2 (C 12 H 26 Cummins model [19] ) was used to model the non-chemical-kinetic processes of diesel fuel, while a reduced n-heptane reaction mechanism with 29 species and 52 reactions [26] was chosen to simulate the diesel fuel chemistry due to the similar ignition characteristics of n-heptane to those of diesel fuel. In addition, the extended Zeldovich mechanism and an improved phenomenological soot model [27] were employed to calculate NO x and soot emissions, respectively. Based on the recent findings of Singh et al. [28] , the CHEMKIN solver [25] was integrated into KIVA-3V code to solve the chemical reaction mechanism without considering the subgrid-scale turbulence-chemistry interaction in this study. The computational model was validated by the experiments in diesel PCCI engine from the work of Lee [29] . The test engine were operated at part-load with indicated mean effective pressure (IMEP) around 5 bar. In order to realize low-temperature combustion, 55% EGR and the compression ratio of 16.0 have been used in the experiment. Fig. 1 shows the comparisons between prediction and measurement for pressure and heat release rate as SOI timing varies from −43 ATDC to −2 ATDC and IVC timing is fixed at −142 °CA ATDC. It can been observed from Fig. 1 that the predicted pressure and heat release rate match the experimental data very well. Further validations on the emissions of HC, CO, NO x and soot were indicated in Fig. 2 a and b. The predicted HC and NO x emissions shows good agreement with the measurements. However, there still exist some discrepancies between the simulation and experiment, and it can be found from Fig. 2 a that the predicted CO emissions are three times smaller than the measured value. Overall, the trends of pressure, heat release rate, as well as HC, CO, NO x and soot with the variation of SOI timings were well reproduced by the simulation, thus the model was applied to investigate the effect of IVC timing and SOI timing on diesel PCCI combustion in the following study. A full 360° mesh comprised of 132,000 cells at bottom dead center (BDC) was used to model the engine investigated in this study, and the computational grid is shown in Fig. 3 . A full 720 °CA engine cycle from exhaust valve opening (EVO) of previous cycle to the EVO of current cycle including the intake and exhaust processes has been calculated, since the time of IVC have an important impact on in-cylinder swirl ratio and tumble ratio development, which further affects the combustion and emissions of PCCI engine [18] . By sweeping the start of injection (SOI) timing from −35 to −5 °CA ATDC, as well as the IVC timing from −140 to −80 °CA ATDC, total 169 cases have been computed to investigate the interactive effects between SOI and IVC. As shown in Fig. 4 , IVC timing was retarded from 220 to 280 °CA (i.e., −140 to −80 °CA ATDC), while the exhaust valve lift profiles and maximum intake valve life were fixed in the tested cases. To focus on the impact of SOI and IVC timing, all other tested conditions were kept constant in this study, including 1.8 bar intake pressure and 50% EGR rate. The simulations were performed on an Intel quad-core Xeon (2.5 GHz) cluster. It spends averaged 58 computational hours for each case, and around 9800 computational hours have been cost for all the calculations. 4 Results and discussions 4.1 Ignition timing Adjusting the IVC timing can significantly affect the in-cylinder pressure and ignition timing due to the variation of effective compression ratio. Fig. 5 shows the effect of late IVC timing on in-cylinder pressure and heat release rate at fixed SOI timing of −20 °CA ATDC. It can be found that ignition timing is efficiently controlled by IVC timing. With late IVC, ignition timing is delayed and the peak in-cylinder pressure is decreased. When the IVC timing of −80 °CA TADC is employed, the combustion phasing locates at TDC. Fig. 6 shows the variation in the in-cylinder pressure and heat release rate traces for various SOI timings at a fixed IVC timing of −120 °CA ATDC. For the cases with SOI timing before −20 °CA ATDC, the pressure traces are quite similar and are independent of SOI timing. This is mainly due to the fact that the premixed combustion occurs in the cylinder as the advanced injection timing is used, and the kinetic reaction rates dominate the combustion processes. As SOI timing is retarded later than −15 °CA ATDC, the air–fuel mixing plays a more important role, thus the ignition timing is delayed as SOI timing is retarded. A 50% burn point map has been developed in Fig. 7 by sweeping SOI timing and IVC timing simultaneously. It is easy to divide the 50% burn points into two regions, i.e., kinetically controlled and mixing-controlled regions as shown in the figure. In the kinetically controlled region, 50% burn point is independent of SOI timing, and is only the function of IVC timing. When IVC timing is varied, burn point is effectively controlled, which is consistent with the results shown in Fig. 5 . As IVC is retarded beyond −115 °CA ATDC, the kinetically controlled region begins to expand. For the case with IVC of −80 °CA ATDC, 50% burn point is not affected by the SOI timing until SOI timing is later than −10 °CA ATDC. However, in the mixing-controlled region, SOI timing dominates the autoignition processes. With the delayed SOI timing, 50% burn point dramatically retarded close to TDC. Further deep understanding of the interesting phenomenon can be found from the distribution of ignition delay in Fig. 8 . It is easy to observe that, the ignition delay is less than 8 °CA when IVC timing and SOI timing lie in the mixing-controlled region. Due to the limited mixing timing, the ignition and combustion processes of PCCI engine are determined by the fuel–air mixing, so that 50% burn point is the only function of SOI timing. However, as the ignition delay is longer than 8 °CA, the time is enough to form the premixed charge in the cylinder, thus IVC timing plays a decisive role in the autoignition processes. 4.2 NO x and soot emissions The NO x emissions for IVC timing and SOI timing sweep are shown in Fig. 9 a. It is found that, at fixed SOI timing, NO x emission decreases with retarding IVC timing, especially for the cases with advanced injection timing. NO x emission shows a more complicated behavior with the varied injection timing due to the spray-bowl interaction. The in-cylinder temperature and equivalence ratio distributions dramatically affects NO x formation, which will be further discussed in the following contexts. As indicated in Fig. 9 a, it can be seen that the delayed SOI timing generally leads to the decrease of NO x emission. The maximum amount of NO x is formed in the region where SOI timing is earlier than −25 °CA ATDC and IVC timing is earlier than −130 °CA ATDC. Fig. 9 b shows the distribution of soot emissions. Two regions with serious soot emissions can be observed. For the cases with very advanced SOI timing and retarded IVC timing, the heavy liquid-spray impingement on the piston wall results in significant increase of soot emissions. When SOI timing is later than −10 °CA ATDC, poor fuel–air mixing in the cylinder also increases the soot emissions. The in-cylinder equivalence ratio and temperature distributions at the crank angle with peak pressure for nine cases indicated in Fig. 9 are shown in Fig. 10 . For further deep understanding of the source of NO x and soot emissions, the in-cylinder equivalence ratio and temperature distributions are plotted on the ϕ – T (equivalence ratio–temperature) map in Fig. 11 . The ϕ – T map constructed in this study is based on the work of Akihama et al. [30] and Bergman et al. [31] by performing calculations in a zero-dimensional, closed perfectly stirred reactor under various temperature and equivalence ratio conditions. According to the conditions of diesel PCCI combustion, all data were recalculated at constant pressure and temperature conditions with 50% EGR and pressure of 10 MPa, which corresponds to in-cylinder peak pressure. The reaction time was considered to be 1.0 ms because it is close to the combustion duration in the diesel PCCI engine investigated in this study. The equivalence ratio ϕ used in this study is calculated based on the intermediate and product composition. The definition for ϕ is the number of oxygen atoms to fully oxidize the available carbon and hydrogen atoms divided by the number of available oxygen atoms as follows: (1) ϕ = ∑ i = 1 N species 2 n C i [ C i ] + 1 2 n H i [ C i ] ∑ i = 1 N species 2 n O i [ C i ] where n C i , n H i and n O i are the number of carbon, hydrogen and oxygen atoms in the i th species, respectively. [ C i ] is the concentration of i th species with the unit of mole/m 3 . As shown in Fig. 11 , NO x volume fraction and soot yield are shown as contours in the ϕ – T map with equivalence ratio ranging from 0 to 4, and temperature ranging from 1500 to 2500 K. The soot yield is defined as percentage of the carbon present as soot referred to the total carbon available. It is evident that most of soot forms in the region where temperature is between 1600 and 2200 K and equivalence ratio is higher than 2.0. The source of NO x mainly lies in the region where temperature exceeds 2150 K and equivalence ratio is lower than 2.0. The equivalence ratio and temperature conditions of each gas cell in the cylinder from 3D-CFD calculation can be plotted on the ϕ – T map to gain understanding of the source of NO x and soot emissions. In order to give a better presentation and a clear comparison for the distributions of equivalence ratio and temperature under various SOI and IVC timings in the ϕ – T map, only the boundaries of the respective cell clouds were plotted in Fig. 11 . It can be observed from Fig. 11 a that, for the cases of SOI = −30 °CA ATDC with IVC timing retarding from −130 to −90 °CA ATDC, a significant reduction of maximum combustion temperature, resulting in the decrease of NO x emissions. However, as late IVC is employed, the lower ambient pressure and temperature produce more wall-wetting, so obvious increase of the region with higher equivalence ratio for IVC = −90 °CA ATDC can be found in Fig. 11 a. Consequently, soot emissions increase dramatically as IVC timing is retarded. For the cases with SOI = −20 °CA ATDC in Fig. 11 b, the combustion temperature is kept below approximately 2200 K which could avoid high NO x concentrations compared to those of SOI = −30 °CA ATDC. The peak combustion temperature slightly decreases as IVC timing is delayed, so NO x emissions could be further reduced. However, the equivalence ratio in local region within the cylinder is beyond 2.0 for all the IVC timings, leading to significant soot emissions, even for the case with IVC = −130 °CA ATDC. The late IVC timing plays a complicated role on the soot emissions. On the one hand, the retarded IVC timing prolongs the ignition delay by decreasing the efficient compression ratio, which helps the mixing of the fuel and air. However, the low ambient pressure caused by late IVC will also increase the wall-wetting, which should be avoided in PCCI combustion. Overall, at the fixed SOI timing of −20 °CA ATDC, the case with IVC = −90 °CA ATDC predicts the lowest soot emissions by providing a relatively best fuel–air mixing before autoignition, but the case with IVC = −110 °CA ATDC leads to the highest soot emissions due to the serious wall-wetting. As SOI timing is retarded to −10 °CA ATDC, incomplete combustion occurring in the cylinder results in considerable local region with rich mixture and consequently the rapid increase of soot emissions as shown in Fig. 11 c. Retarding IVC timing reduces NO x emissions which is consistent with the cases of earlier SOI timing, but the benefit of late IVC with SOI = −10 °CA ATDC on NO x emissions is not as significant as those of SOI = −30 °CA ATDC due to the relatively low combustion temperature. 4.3 HC and CO emissions The contour plots of HC and CO emissions over SOI and IVC timing are shown in Fig. 12 . As SOI timing is retarded later than −10 °CA ATDC, HC and CO emissions increase rapidly due to the incomplete combustion resulting from the low compression ratio and limited mixing time. Peak HC and CO also occur in the region where SOI timing is earlier than −30 °CA ATDC and IVC timing is later than −110 °CA ATDC. This can be explained by the fact that serious wall-wetting happens at the end of combustion phasing when the ambient pressure is low. It is seen from Fig. 13 that the liquid fuel impinging the piston wall increases rapidly when advanced SOI timing and late IVC timing are employed. For the case with SOI timing of −35 °CA ATDC and IVC timing of −80 °CA ATDC, as high as 11.4% liquid fuel still leaves on the piston surface at the end of combustion stage. By summarizing the information about NO x , soot, HC and CO emissions from Figs. 9 and 12 , four major regions can be identified for various SOI and IVC timings in Fig. 14 . The region with SOI timing later than −10 °CA ATDC is characterized by the incomplete combustion. As IVC timing is retarded to −80 °CA ATDC, the incomplete region expands to SOI timing of −17 °CA ATDC due to the low combustion temperature. The region with advanced SOI timing and early IVC timing is characterized by NO x region because of the highest NO x emissions. Using retarded IVC timing and late SOI timing can significantly reduce NO x emission by decreasing the combustion temperature. However, the cases with retarded IVC timing and early SOI timing lie in the serious wall-wetting region due to the low ambient pressure. Therefore, the final clean combustion is confined in a small region where SOI timing and IVC timing must be carefully adjusted. Noted that the clean combustion shown in Fig. 14 is defined as the region where NO x and soot emissions are less than 1.0 g/kg-fuel and 3.0 g/kg-fuel, respectively, and HC and CO emissions are kept at lower levels simultaneously. For the cases with IVC timing earlier than −120 °CA ATDC, SOI timing should be restricted in a narrow range around −10 °CA ATDC to avoid all of NO x , soot, HC and CO emissions. When IVC timing is retarded later than −110 °CA, the operating range of SOI becomes wider. It is beneficial for the control of PCCI combustion to keep the stable operation. 4.4 Fuel consumption The contour plot for indicated specific fuel consumption (ISFC) is shown in Fig. 15 . The ISFC is calculated as: (2) ISFC = m f / W i = m f / ∮ p · dV where m f is the injected fuel mass per cycle, W i is the indicated work per cycle, which is obtained by integrating in-cylinder pressure with respect to corresponding in-cylinder volume throughout the entire four-stroke cycle. As can be expected, the regions with very high CO and HC emissions lead to increased ISFC, including the serious wall-wetting and incomplete combustion regions. Two pursued regions with as low as 185 g-fuel/kWh ISFC can be observed in the figure. For advanced SOI timing around −25 °CA ATDC, retarded IVC timing improves the fuel efficiency. This is mainly due to the fact that the early ignition timing always happens when early SOI timing is employed, which increases the compression work and decreases ISFC. When late IVC timing is used, the ignition timing can be efficiently delayed to TDC which helps the decrease of the compression work, thus ISFC is reduced. However, for the conditions with SOI timing of −15 °CA ATDC, early IVC timing is helpful for the decrease of ISFC. This is because the combustion efficiency dominates the ISFC, so the increased combustion temperature with early IVC is beneficial for the improvement of fuel efficiency. 4.5 Engine knock Another issue should be overcome for PCCI engine is the engine knock. With the early injecting timing, the high cetane number of diesel fuel can lead to early ignition and serious engine knock. In order to evaluate the influence of SOI and IVC timing on engine knock, the ringing intensity developed by Eng [32] was used in this study, which is calculated as: (3) Ringing Intensity ≈ 1 2 γ · ( 0.05 · dP dt max ) 2 P max · γ RT max where P max is the maximum cylinder pressure, ( dP / dt ) max is the maximum rate of pressure rise, γ is the ratio of specific heat, R is the ideal gas constant, and T max is the maximum cylinder temperature. The variation of ringing intensity with SOI and IVC timing is shown in Fig. 16 . Compared to the trends of emissions and ISFC, the ringing intensity shows more complicated behavior. As SOI timing is advanced beyond −20 °CA ATDC in the cases of −120 °CA ATDC IVC, the ignition timing is fixed at −10 °CA ATDC as shown in Fig. 6 . Therefore, with early SOI timing, the more homogenous charge is formed and simultaneous combustion occurs, which increases the ringing intensity. However, as SOI timing is further retarded later than −17 °CA ATDC, the ignition timing is near TDC, part of the local rich mixture autoignition also raises the ringing intensity. The optimized SOI timing for the minimum ringing intensity is around −20 °CA ATDC. The delayed IVC timing does not help the reduction of ringing intensity, although the combustion temperature is efficiently decreased. This is mainly due to the fact that less air is induced in the cylinder as intake valve is closed for the cases with late IVC timing. Therefore, the peak heat release rate increases slightly with retarded IVC timing in Fig. 5 , although the ignition timing is delayed. Boosting the intake pressure could effectively solve this problem, which is beyond the scope of this study. 5 Conclusions The effect of injection timing and intake valve close timing on combustion and emission characteristics of a diesel PCCI engine has been studied by a full-cycle three-dimensional CFD model coupled with detailed chemical kinetics. By sweeping SOI timing from −35 to −5 °CA ATDC, as well as the IVC timing from −140 to −80 °CA ATDC, the contour plots for NO x , soot, HC, CO, ISFC and ringing intensity have been developed. The results show that the kinetically controlled region and mixing-controlled regions can be divided based on the 50% burn point distributions for various SOI and IVC timings, in which the ignition timing are controlled by IVC timing and SOI timing respectively. Using retarded IVC timing and optimized SOI timing can significantly reduce NO x emission by decreasing the combustion temperature. However, the cases with retarded IVC timing and early SOI timing lie in the serious wall-wetting region due to the low ambient pressure. The region with very late SOI timing is characterized by incomplete combustion, and the incomplete combustion region expands with retarded IVC. Overall, the operating range of SOI becomes wider for clean combustion with the retarded IVC timing. The cases with advanced SOI timing around −25 °CA ATDC and retarded IVC timing, as well as SOI timing of −15 °CA ATDC and early IVC timing are helpful for the decrease of ISFC due to the optimization of ignition timing and combustion efficiency. In order to avoid engine knock, the SOI timing should be kept around −20 °CA ATDC to minimize the ringing intensity. The delayed IVC timing does not help the reduction of ringing intensity, although the combustion temperature is efficiently decreased. Acknowledgements This work is supported by the National Natural Science Foundation of China (Grant No. 50806008 ), National Basic Research Project of China (Grant No. 2007CB210002 ) and State Key Laboratory of Engines, Tianjin University (Grant No. K2010-03). References [1] F. Zhao T.W. Asmus D.N. Assanis J.E. Dec J.A. Eng P.M. Najt Homogeneous charge compression ignition (HCCI) engines: key research and development issues 2003 Society of Automotive Engineers Warrendate, PA [2] J.E. Dec Advanced compression-ignition engines – understanding the in-cylinder processes Proc Combust Inst 32 2 2009 2727 2742 [3] Walter B, Gatellier B. Development of the high-power nadist concept using dual-mode diesel combustion to achieve zero nox and particulate emissions. SAE paper 2002-01-1744; 2002. [4] B. Gatellier A. Ranini M. Castagné New developments of the nadi concept to improve operating range, exhaust emissions and noise Oil Gas Sci Technol 61 1 2006 7 23 [5] Lechner GA, Jacobs T, Chryssakis C, Assanis DN, Siewert RM. Evaluation of a narrow spray cone angle, advanced injection timing strategy to achieve partially premixed compression ignition combustion in a diesel engine. SAE paper 2005-01-0167; 2005. [6] Fang T, Lee CF, Coverdill RE, White RA. Smokeless combustion within a small-bore hsdi diesel engine using a narrow angle injector. SAE paper 2007-01-0203; 2007. [7] Siewert RM. Spray angle and rail pressure study for low nox diesel combustion. SAE paper 2007-01-0122; 2007. [8] Helmantel A, Denbratt I. HCCI operation of a passenger car di diesel engine with an adjustable valve train. SAE paper 2006-01-0029; 2006. [9] Nevin RM, Sun Y, Gonzalez MA, Reitz RD. PCCI investigation using variable intake valve closing in a heavy duty diesel engine. SAE paper 2007-01-0903; 2007. [10] He X, Durrett R, Sun Z. Late intake valve closing as an emissions control strategy at tier 2 bin 5 engine-out nox level. SAE paper 2008-01-0637; 2008. [11] B. Walter P. Pacaud B. Gatellier Variable valve actuation systems for homogeneous diesel combustion: how interesting are they? Oil Gas Sci Technol 63 4 2008 517 534 [12] H. Kim J. Lee K. Kim K. Lee Effect of the atkinson cycle combined with calibration factors on a two-stage injection-type premixed charge compression ignition engine Energy Fuels 23 10 2009 4908 4916 [13] C.L. Genzale S.C. Kong R.D. Reitz Modeling the effects of variable intake valve timing on diesel HCCI combustion at varying load, speed, and boost pressures J Eng Gas Turb Power 130 5 2008 0528061 0528068 [14] Murata Y, Kawano D, Kusaka J, Daisho Y, Suzuki. H, Ishii H, Goto Y, Odaka M. Achievement of medium engine speed and load premixed diesel combustion with variable valve timing. SAE paper 2006-01-0203; 2006. [15] Murata Y, Kusaka J, Daisho Y, Kawano D, Suzuki H, Ishii H, Goto Y. Miller-PCCI combustion in an hsdi diesel engine with VVT. SAE paper 2008-01-0644; 2008. [16] Y. Murata Y. Nishio J. Kusaka Y. Daisho D. Kawano H. Suzuki Numerical analysis of miller-premixed charge compression ignition combustion on a dynamic ϕ – t map Int J Eng Res 11 2 2010 89 98 [17] Boyarski NJ, Reitz RD. Premixed compression ignition (PCI) combustion with modeling-generated piston bowl geometry in a diesel engine. SAE paper 2006-01-0198; 2006. [18] Z. Peng M. Jia Full engine cycle cfd investigation of effects of variable intake valve closing on diesel PCCI combustion and emissions Energy Fuels 23 12 2009 5855 5864 [19] Amsden AA. KIVA-3v: A block-structured KIVA program for engines with vertical or canted values. Los Alamos National Laboratory Technical Report, LA-13313-MS; 1997. [20] Z.Y. Han R.D. Reitz Turbulence modeling of internal combustion engines using rng κ–ε models Combust Sci Technol 106 4–6 1995 267 295 [21] L.M. Ricart R.D. Reitz J.E. Dec Comparisons of diesel spray liquid penetration and vapor fuel distributions with in-cylinder optical measurements J Eng Gas Turb Power 122 4 2000 588 595 [22] Nordin N. Complex chemistry modeling of diesel spray combustion. PhD thesis, Chalmers University; 2001. [23] Z.Y. Han Z. Xu N. Trigui Spray/wall interaction models for multidimensional engine simulation Int J Eng Res 1 1 2000 127 146 [24] Z.Y. Han R.D. Reitz A temperature wall function formulation for variable-density turbulence flows with application to engines convective heat transfer modeling Int J Heat Mass Transfer 40 3 1997 613 625 [25] Kee RJ, Rupley FM, Meeks E, Miller JA. Chemkin-iii: A fortran chemical kinetics package for the analysis of gas phase chemical and plasma kinetics. Sandia National Laboratories, SAND96-8216; 1996. [26] Patel A, Kong SC, Reitz RD. Development and validation of reduced reaction mechanism for HCCI engine simulation. SAE paper 2004-01-0558; 2004. [27] M. Jia Z.J. Peng M.Z. Xie Numerical investigation of soot reduction potentials with diesel homogeneous charge compression ignition combustion by an improved phenomenological soot model Proc Inst Mech Eng Part D J Automob Eng 223 3 2009 395 412 [28] S. Singh R.D. Reitz M.P.B. Musculus T. Lachaux Validation of engine combustion models against detailed in-cylinder optical diagnostics data for a heavy-duty compression-ignition engine Int J Eng Res 8 1 2007 97 126 [29] Lee S-S. Investigation of two low emissions strategies for diesel engines: Premixed charge compression ignition (PCCI) and stoichiometric combustion. PhD Thesis, University of Wisconsin-Madison; 2006. [30] Akihama K, Takatori Y, Inagaki K, Sasaki S, Dean AM. Mechanism of the smokeless rich diesel combustion by reducing temperature. SAE paper 2001-01-0655; 2001. [31] Bergman M, Golovitchev V. Application of transient temperature vs. Equivalence ratio emission maps to engine simulations. SAE paper 2007-01-1086; 2007. [32] Eng JA. Characterization of pressure waves in HCCI combustion. SAE paper 2002-01-2859; 2002.