Standard specifications of test engine.
Abstract
In order to satisfy the growing demand for the reduction of fuel consumption and pollutant emissions, various technologies have been employed in diesel engines. Consequently, to determine the optimal combustion control strategy, many parameters such as injection pressure, nozzle diameter, injection timing, injection quantity, and exhaust gas recirculation (EGR) rate should be selected properly corresponding to the engine operating conditions. It is difficult to obtain the appropriate strategies without understanding the change in combustion process when varying these parameters. To realize parametric studies on combustion control strategy of modern diesel engines, a phenomenological combustion model based on stochastic method was developed. In this model, the modeling of the spray tip and tail penetration after the end of injection, and interaction between the sprays of sequent injection stages were focused on to modify the stochastic combustion model for combustion simulation with multiple injection. The effects of swirl, wall impingement, and adjacent spray interaction are formulated simply to make the combustion model more accurate and computationally efficient. The simulation results were compared with experimental data from a single-cylinder test engine for pilot/main two-stage injection. The results reveal that the model has capability to accurately predict the combustion characteristics and emissions of diesel engine with pilot/main two-stage injection.
Keywords
- diesel engine
- combustion modeling
- phenomenological model
- possibility density function
1. Introduction
Recently, continuous reduction of the oxides of nitrogen (NOx) and PM emissions for confronting increasingly stringent emission standards becomes a significant task for combustion improvement of diesel engines; especially the reduction of NOx emissions is quite difficult for aftertreatment for low-load operating conditions due to low exhaust gas temperature. Hence, in-cylinder combustion improvement is required for low-load operation, subsequently the advanced combustion concepts have been proposed for dealing with this situation, such as premixed charge compression ignition (PCCI) and low-temperature combustion (LTC). No matter whether it is PCCI or LTC, the main aim is to obtain more homogeneous mixture and suppress locally high combustion temperature. This aim is typically achieved with injection timing adjustment with high exhaust gas recirculation (EGR) rate. However, accompanying with reduction of NOx and PM, the CO and HC emissions are observed to increase and even exceed the emission standards [1–4]. Meanwhile, long ignition delay promotes the fuel-air mixing before ignition; it leads to the strong dependence of ignition timing on chemical reaction, which makes the ignition timing control difficult. And fast combustion of the mixture premixed before ignition increases the combustion noise associated with rapid rise of in-cylinder pressure. Because of these barriers, the advanced combustion concepts are limited in low-load operation.
To extend the advanced combustion concepts to high-load operation and obtain the best reduction of emissions, many technologies, including high-injection pressure, small nozzle diameter, multistage injection, EGR, and high-intake boost, have to be integrated to realize in-cylinder combustion improvement. Consequently, to determine the optimal combustion control strategy, many parameters such as injection pressure, nozzle diameter, injection timing, injection quantity, and EGR rate should be selected properly corresponding to the engine operating conditions. It brings a heavy work and high cost on experiments, moreover it is difficult to obtain the appropriate strategies without understanding the change in combustion process when varying these parameters. Thus, an accurate and computationally efficient combustion model is needed for these parameters study. Considering the trade-off between accuracy and computational efficiency, the phenomenological combustion models are often adopted for the parametric study [5].
In phenomenological combustion models, the combustion process is mimicked mathematically based on the results of observed phenomena with focus on predominant and rate-determining processes in combustion. To describe the heterogeneous distribution of fuel concentration and temperature and low-temperature combustion phenomena, some detailed phenomenological diesel combustion models are developed based on the spray behavior, which incorporate the calculation of chemical reaction in the spray region. According to the method for description of fuel concentration distribution, these models can be divided into two types. One type is to divide the spray into different zones based on the spray structure [6–8]. In the other type, fuel concentration distribution is represented by probability density function (PDF) of fuel mass fraction [9, 10]. PDF method is dominant to couple the turbulent mixing and chemical kinetics to calculate the turbulent combustion process in diesel engines, and it is easy to blend with the emission models, such as NOx and soot models. Whence the pollutants production process can be predicted and analyzed with varying relevant parameters.
In modern diesel engines, multiple-injection has played a significant role in the latest engines to reduce NOx, smoke emissions, and combustion noise simultaneously [11–13]. Proper selection of injection strategy from numerous injection parameters satisfying different operating conditions is a challenging work. Thus, to estimate the modern diesel engines combustion, the combustion process with multiple injections needs to be included in the phenomenological combustion model. Recently, some models had been proposed based on PDF method to calculate the combustion with multiple-injection [14–16]. In these models, each injection is customarily treated to form a spray region, and relevant calculations including spray propagation, fuel evaporation, turbulent mixing, ignition, and chemical reaction are performed in each spray. According to the sprays propagation phenomena under multiple-injection situation, the combustion of the latter spray is significantly dependent on the thermodynamic states of the former spray. Therefore, the interaction between former spray and latter spray has to be considered in the model. Normally, the latter injection starts after the end of the former injection, so that the evolution of the spray after the end of injection (EOI) has to be imitated. In fact, the spray behaviors after EOI are indicated quite different from that during injection by the fundamental studies on turbulent gas and diesel jets [17–21]. However, there is no simple model available in phenomenological models to describe the spray development and air entrainment after EOI. Currently, frequently used models of spray propagation are Wakuri’s model, Siebers’s model, and Hiroyasu’s model that only satisfy the spray propagation during injection with steady state and do not involve any information for that after EOI.
According to the modeling concept from the developed models [14–16], the latter spray is assumed to be injected into the former spray region, so that the entrained gas of the latter spray only comes from the former spray zone. It should be noted that all of these interaction treatments assume that later spray entrains the gas from the former spray zone immediately when the later injection starts, whereas some fundamental studies discover that the mixture near the injector tip is diluted very fast after EOI [20, 22]. This implies that the gas near the injector tip should be very lean when the later injection starts, particularly for the early pilot injection case. Thereby, it is necessary to consider this phenomenon in appropriate form to imitate the interaction between sprays from the sequent injections. Meanwhile, the swirl flow deflects the spray from the initial injection direction [23]. In multiple-injection case, the magnitudes of deflection are different for former and latter sprays, which results in that the swirl effects has to be involved for spray-to spray interaction modeling.
In addition, the phenomenological combustion models usually do not consider the spatial information. However, the early injection is a typical strategy for advanced combustion concepts, especially quite early pilot injection for multiple-injection strategy. The early injection causes the wall impingement, consequently the effects of different impinging positions and the spray flames propagation along the cylinder wall on fuel-air mixing need to be involved in the model.
Thus, in this study, a phenomenological combustion model was developed based on a stochastic mixing model. The heterogeneity of fuel concentration and temperature is described using the PDF. Evolution of the mixing and combustion is represented by the temporal change of PDF due to the mixing, fuel injection, fuel evaporation, air entrainment into a spray, heat loss to the combustion chamber walls, and chemical reactions. To focus on the simulation of new combustion concepts, such as LTC and PCCI, the effects of EGR and multiple-injection was involved in the model. The model employs a quasi-global chemical kinetics model combined with stochastic mixing model to describe the low-combustion process caused by EGR effect. The key points of the study concentrated on are as follows: (1) Spray model developing to describe the spray evolution after end of injection (EOI). (2) Defining and calculating the presumed spray tail to consider the quite lean mixture near the nozzle exit after EOI. (3) Modeling spray-to-spray interaction between sequent injections and the swirl flow effects on it. (4) Introducing the effects of spray-wall interaction on fuel-air mixing. Finally, the simulations are performed in pilot/main two-stage injection cases, and are validated by experimental data from a single-cylinder diesel engine, including in-cylinder pressures, heat release rates, and NOx and PM emissions.
2. Stochastic combustion model
The stochastic combustion model in this study was first proposed in [10] for conventional combustion with single-stage injection in DI diesel engines. And then this model was modified in [24–26] to analyze the PCCI combustion process. Up to Ref. [26], this model has employed a stochastic turbulent mixing model coupled with a quasi-global chemical kinetic model to describe the heterogeneous combustion process in diesel sprays, and the specific submodels includes air entrainment, turbulent mixing, fuel vaporization, heat loss, and chemical kinetic model.
2.1. Model description
The essential of the stochastic combustion model focus on the turbulence controlled fuel-air mixing process. In this model, the heterogeneity in fuel concentration and temperature is described using probability density functions
where (∙) expresses the time derivative,
where
The
Based on this concept, the combustion process of diesel spray is described as shown in Figure 1. In this process, the combustion chamber is divided into spray zone and surrounding air zone after the injection start. The spray zone is formed by injected fuel particles and entrained air particles whose amount is determined by momentum theory [28]. In the spray zone, the fuel-air mixing is taken into account with modified Curl’s collision and dispersion model [27]. The mixing rate is controlled by stochastic collision frequency that is dependent on turbulence kinetic energy in spray. Meanwhile, liquid fuel contained in each fluid particle is vaporized according to an average evaporation rate, sequentially chemical reaction calculation using Schreiber’s model [29] is performed in each fluid particles for combustion calculation. In addition, the heat loss from wall is also considered using Woschni’s equation [30]. Thus, in the stochastic combustion model, fuel-air mixing and chemical reactions are simultaneously calculated.
3. Spray model
Figure 2 shows the shadowgraph images of two-stage injection in a constant volume combustion chamber. In Figure 2, the part “A” is the first spray flame region. After the second injection starts, the first spray flame region and the second spray region are separated at first as shown in the left graph. Once the second spray tip penetration is long enough to overtake the boundary of the first spray flame region, the second spray subsequently penetrates into the first spray flame region and causes the interaction between the first spray flame and the second spray as shown in the right graph.
To simulate the combustion with multiple-injection, the situations shown in Figure 2 can be abstracted as that shown in Figure 3. The volume in the cylinder is divided into three zones to represent the sprays and ambient gas for two-stage injection situation, including the first spray zone, the second spray zone, and the ambient air zone. Before injection, the air and EGR gas in the cylinder is treated as ambient air zone that is presumed to be composed of a great number of fluid elements. Then after the start of the first injection, the fuel is injected into the cylinder as fluid elements, and entrained fluid elements from ambient air zone to form the first spray zone. After the end of the injection, the first spray tail is assumed to depart from nozzle exit to the downstream. When the second injection starts, the first and second spray zones exist in the cylinder and the fluid elements of ambient air zone are entrained into these two spray zones, respectively. Once the second spray tip overtakes the first spray tail, both of the fluid elements from the ambient air zone and the first spray zone are entrained into the second spray zone. Finally, these two spray zones are combined into one spray zone after all of the fluid elements of the first spray zone are entrained into the second spray zone. Obviously, the spray tip and tail penetration and air entrainment rate are primary parts to carry out the concept, and the spray tip penetration should include the spray evolution after EOI as the discussion in Section 1. Thereby, a new zero-dimensional spray penetration calculation was developed, and the air entrainment was calculated based on the developed model.
3.1. Spray penetration including the spray evolution after EOI
Indeed, because of the termination of injection, the spray is transformed from steady state into deceleration state that leads to the spray propagation different from that in the steady state. Recently, entrainment wave theory was proposed by Musculus [18] based on analytical research on the turbulent jets behaviors during the deceleration state. This theory describes the propagation of the increased entrainment region from nozzle exit to jet downstream after EOI, which also presents the propagation of information of the fuel injection termination. Once the entrainment wave arrives at the tip of the jet, the jet penetration shifts from a relation of square root dependence on time to a relation of fourth root dependence on time, the former relation is widely proved in steady jet and the later relation is same as the result from the experimental research on unsteady turbulent jets [17]. The entrainment wave was also observed in diesel jets, and a one-dimensional discrete model derived in [19] was successful to predict the entrainment wave in diesel jets and the penetration of diesel jets even after the entrainment wave reaches the jet tip. These indicate that the one-dimensional discrete spray model in [19] has potential to predict the diesel spray propagation after EOI. Thus, a zero-dimensional model of diesel spray propagation was developed referring to this one-dimensional discrete spray model to involve the spray information of penetration and air entrainment rate after EOI in this study. To conduct this work, the modeling method of the one-dimensional discrete model is only used for the treatment of the spray tip part, and the same assumptions as in the one-dimensional discrete spray model were employed for the zero-dimensional spray model. To simplify the spray modeling, the spray is treated to be formed with a constant injection rate.
According to the entrainment wave phenomenon, the spray propagation process can be divided into two stages. As shown in Figure 4, the time when the entrainment wave front arrives at the spray tip is defined as transition timing, and the diesel spray tip penetration calculation is separated in two stages by the transition timing. During injection, the spray is in steady state; after the EOI, the entrainment wave propagates from nozzle exit to the spray tip, and the spray tip keeps in the steady state before the transition timing. After the transition timing, the total spray transforms into decelerating state. In this spray propagation process, the spray tip maintains in steady state before transition timing and then decelerates from the steady state after transition timing. Thus, the spray tip penetration is calculated in two stages, before and after transition timing. In the case of constant injection rate, the transition timing is twice of injection duration.
The spray tip penetration (
where
where
Meanwhile
Substituting Eqs. (19) and (20) into (21), the spray tip velocity is represented as follows:
Before the transition timing, the spray tip is in steady state, so that the momentum flux and fuel mass flow rate integrated over the tip cross-sectional area
where
Therefore, after substituting Eq. (6) into (9), the
where
So far the spray tip penetration (
Eq. (12) cannot be solved as a continuous equation. Thus, the discrete method is used to calculate the spray tip penetration. The spray tip penetration at any time of
And the average velocity over the tip cross-section at time
After the transition timing, the total momentum flux and fuel mass flow rate over the tip cross-section decelerate from the values of steady state because the information of fuel injection termination already arrived at the spray tip [18]. However, Eq. (9) is valid also for the period after the transition timing, which had been demonstrated in [31]. Thus, the same form as Eq. (10) are available for the average velocity over the tip cross-section, where
where
where
To obtain
where
Whence the spray tip penetration after the transition timing is able to be obtained, and the method for solving Eq. (16) is same as that before transition timing.
3.2. Spray tail penetration
After EOI, the mixture near the nozzle becomes very lean due to the termination of the fuel supply. Therefore, it is assumed that the part of the spray near the nozzle acts as ambient air zone and a spray tail exists. Knowing the position of the spray tail, the start time of the interaction with the subsequent spray can be determined.
In this study, the spray tail position is determined as a cross-section of a spray where 10% of total fuel is contained up to the nozzle. It can be imagined that the tip of the spray formed by this 10% fuel equates the total spray tail. In other words, the tail penetration of the spray containing total injected fuel can be calculated by the tip penetration of the spray formed by an injection whose quantity is equal to 10% of total injection quantity. Figure 5 shows the examples of calculated tail penetrations compared to the tip penetrations for the injection quantity of 2.5 and 25 mg cases.
4. Spray-to-spray interaction modeling
Another primary respect of simulation for diesel combustion with multiple-injection is to describe the interaction between the sprays from sequent injections. Figure 6 shows the spray-to-spray interaction, the interaction occurs after the arrival of the second spray tip at the first spray tail, and the entrainment from the first spray to the second spray is involved to represent the interaction.
The entrainment behavior can be considered as that the spray entrains the ambient gas through the spray boundary that has a conical surface, therefore the entrainment rate can be represent as the product of air density, entrainment velocity over the spray boundary surface, and the area of spray boundary surface. After the second spray tip touches the first spray tail, the entrainment area of the second spray is divided into two parts by the first spray tail as shown in Figure 6. The ratio of the entrainment rate of these two parts (
where
Meanwhile, the swirl effects on the spray-to-spray interaction cannot be ignored. The swirl flow does not only decrease the spray penetration [32] but also deviates the spray path [23]. As that shown in Figure 7, the swirl effects causes the second spray tip to overtake the first spray tail and tip earlier than the case without swirl flow effects, and reduce the overlap region between the first and second sprays. Thus, based on the assumption that the entrainment rate is proportional to the spray surface area, the ratio (
where
To calculate
where
To calculate
The ratio (
5. Wall impingement effects introduced in stochastic combustion model
5.1. Wall impingement effects on the turbulent mixing
Adjusting the injection time earlier than top dead center (TDC) is often used to realize the PCCI or LTC, especially for the multiple-injection case in which the pilot injection time is usually advanced to the middle even at the early stage of compression stroke. The early injection timing makes the spray flow into the squish region and impinges on the cylinder liner or piston top as shown in Figure 9. Due to the low temperature of the walls and/or the adherence of fuel on the piston top surface, oxidation reaction and mixing in the mixture are attenuated [35]. To involve such effects in the stochastic combustion model, the reduction of fuel-air turbulent mixing rate was considered according to the volume ratio of the spray flowed into the squish area and total spray. The temperature effect mentioned above was not considered; however, the overall oxidation reaction rate is lowered as a result of reduced mixing rate.
As shown in Figure 9, if the spray tip cross-section impinges on the bowl lip edge, the spray can be divided into two parts, squish part and bowl part. Thus the volume ratio between squish part and incremental spray (
where
where
where
5.2. Wall impingement effects on the air entrainment
Figure 10 shows abstracted diagram of the situation of spray impinging on the wall in the cylinder. The upper spray represents the initial stage of wall impingement. In this stage, it is reasonable to consider that the air entrainment is enhanced because the surface area of spray is enlarged by wall impingement [36–38]. Thus, to improve the stochastic combustion model, a constant (
where
6. Emission models
6.1. NOx model
Based on the assumption that most of the NOx is NO, the production of NO during combustion process is computed by the extended Zeldovich mechanism [39]. Meanwhile, the NO normally produces at a high rate in the mixture with equivalence ratio around 1.0 accompanying with high temperature. Thus, the NO concentration is estimated in the fluid elements whose temperature is over 1200 K with the equilibrium species including C, CO, CO2, O2, O, OH, H2, H, H2O, N2, and NO.
6.2. Soot model
The soot model refers to the Moss’s soot model [40]. Moss’s soot model is a semi-empirical soot model derived based on laminar diffusion flame. The soot particles inception and coagulation are considered for the calculation of soot particles number density, and the soot particles surface growth and oxidation are calculated to obtain the soot volume fraction. The soot oxidation rates per unit area by O2 and hydroxyl radical (OH) are introduced to calculate the soot oxidation rate in this model. The soot oxidation rate per unit area by O2 (
7. Validation of the developed combustion model
7.1. Test engine
The test engine is a water-cooled single-cylinder four-stroke-cycle direct-injection diesel engine equipped common-rail injection system. The standard specifications are given in Table 1.
Engine type | Single-cylinder, DI-diesel engine |
Bore × Stroke | 85 mm × 96.9 mm |
Displacement | 550 cc |
Compression ratio | 16.3 |
Combustion chamber | Reentrant type (51.6 cavity) |
Injection system | Common-rail system 0.125 mm × 7 holes nozzle |
All of the experiments were performed at thermally steady states of the engine at a fixed speed of 1500 rpm, an inlet coolant temperature of 80°C, and a lubricating oil temperature of 80°C. The intake pressure was 0.1 MPa and the intake temperature was 35°C. Exhaust back pressure valve was fully open. The fuel was JIS No. 2 diesel fuel (density at 15°C = 820.2 kg/m3 and cetane index = 54.7). The averaged in-cylinder pressure of 50 cycles was used to calculate the heat release rate, which was measured using a piezoelectric pressure transducer (Kistler 6052A).
7.2. Calculation conditions
With the aim for validating the simulation of combustion with multistage injection, at first, the pilot/main two-stage injection strategies were conducted in experiment and calculation, and the calculation ran from intake valve closure (IVC) of –145°ATDC to exhaust valve open (EVO) of 125°ATDC. Since the start of main injection is normally set near TDC in the pilot/main injection strategy, the main injection timing (
Injection pressure | 125 MPa |
Total injection quantity | 32 mm3 per cycle |
Pilot injection quantity ( |
2, 4, and 6 mm3 per cycle |
Pilot injection timing ( |
–9, –19, and –24°ATDC |
Main injection timing | 1°ATDC |
EGR ratio | 20% |
Swirl ratio | 2.0 |
7.3. Pressure and heat release rate
To valid the combustion model, the in-cylinder pressures and the heat release rates were calculated and compared to the experimental data at first. Figures 11–13 show the in-cylinder pressures and heat release rates for different pilot injection timings and quantities. The simulation results obtain the similar levels of pressures and heat release rates to those of the experiments. It is also observed that the model is able to capture the tendencies of the pressure and heat release rate when varying the pilot injection timing and injection quantity.
7.4. Emissions calculation
Based on the good agreement in the pressures and the heat release rates, the NOx and soot emissions were calculated. Figure 14 shows the NOx (left) and soot (right) emissions against the pilot injection timing. The results reveal that the NOx calculation obtains the similar emission level and the variation caused by pilot injection conditions changing with the measured data. Regarding to the soot emission, the soot emissions level obtained by the model is comparable with the measured data, and the calculated soot emissions are able to reproduce the increase with the increase in pilot injection quantity at fixed pilot injection timings, which is observed in the experimental data. Although the tendency of soot emissions variation with pilot injection timing retarded does not completely coincide with that of every experiment with different pilot injection quantity, the general tendency can be captured correctly.
8. Conclusion
In this study, a stochastic combustion model was introduced to develop a phenomenological combustion model for modern diesel engines. In order to be able to describe the combustion process of advanced combustion mode, the spray propagations after EOI, spray-to-spray interaction, swirl effects, and wall impingement effects were modeled based on their physical phenomena in appropriate ways for the stochastic combustion model. Then, the developed combustion model was validated based on the experimental data from a single-cylinder diesel engine with pilot/main two-stage injection. The results revealed that the model is able to accurately predict the combustion of the diesel engine with pilot/main two-stage injection, and reasonable prediction of NOx and soot emissions can be obtained by this model. Specific conclusions are as follows:
(1) A zero-dimensional spray propagation model was derived. The spray model is able to predict the spray evolution including spray tip penetration and overall air entrainment after EOI. It is capable to catch the spray propagation tendency after EOI as that well recognized by the fundamental study. Thus, based on introducing this model, the accuracy of the combustion model can be improved especially for the short injection duration case in which the ignition is later than the EOI. As for the multiple-injection case, the thermodynamic states of the gas in the cylinder before the later injection start can be predict more realistically.
(2) Presumed spray tail was proposed for zero-dimensional spray model, which allows to take account of the fast dilution of the mixture near the injector tip after EOI, which resulted by the terminated fuel supply and increased air entrainment rate near the injector tip after EOI. And thanks to the presumed spray tail, the interaction between sprays from sequent injections was formulated by the rate of entrainment into the later spray from the former spray and/or the surrounding gas. Meanwhile, spray deflection by swirl flow cannot be neglected when imitating the interaction between sprays from sequent injection. And the interaction with swirl flow effect was described based on a simple geometrical consideration and formulated using the momentum from injected fuel and entrained gas. In this way, the stochastic combustion model can be used for simulation of combustion with multistage injection.
(3) To consider the wall impingement effect in early injection timing case, after wall impinging, the volume ratio of the spray flowed into the squish area, and total spray was introduced to reduce the fuel-air mixing rate that can involve the combustion chamber shape effects in the phenomenological combustion model.
(4) Air entrainment rate enhancement caused by wall impingement was considered by a factor with air entrainment rate of free spray during initial stage of wall impingement. And the reduction of main spray entrainment rate by interaction between adjacent sprays was formulated by the ratio between spray volume and chamber volume to balance the effect of enhanced air entrainment rate in the late combustion period. According to the results, effects of wall impingement and adjacent spray interactions on the entrainment rate are helpful to well predict the heat release rate in initial and late combustion periods, respectively.
Acknowledgments
The experiments and modeling work presented here had been performed in Combustion and Power Engineering Laboratory at Kyoto University. The research is supported by National Natural Science Foundation of China (Grant No. 51509051) and Natural Science Foundation of Heilongjiang Province of China (Grant No. LC2015017).
References
- 1.
Colban, W. F., Miles, P. C., and Oh, S. Effect of intake pressure on emissions from an automotive diesel engine operating in low temperature combustion regimes, SAE Paper 2007-01-4063, 2007. - 2.
Lee, S. and Reitz, R. D. Spray targeting to minimize soot and CO formation in premixed charge compression ignition (PCCI) combustion with a HSDI diesel engine, SAE Paper 2006-01-0918, 2006. - 3.
Cho, K., Han, M., Sluder, C. S., Wagner, R. M., and Lilik, G. K. Experimental investigation of the effects of fuel characteristics on high efficiency clean combustion in a light-duty diesel engine, SAE Paper 2009-01-2669, 2009. - 4.
Ojeda, W., Zoldak, P., Espinosa, R., and Kumar, R. Development of a fuel injection strategy for diesel LTC, SAE Paper 2008-01-0057, 2008. - 5.
Lakshminarayanan, P. A. and Aghav, Y. V. Modeling diesel combustion, Springer, Dordrecht/Heidelberg/London/New York, 2009. - 6.
Shahed, S. M., Chiu, W. S., and Yumlu, V. S. A preliminary model for the formation of nitric oxides in DI diesel engine and its application in parametric studies, SAE Paper 730083, 1973. - 7.
Dent, J. C. and Mehta, P. S. Phenomenological combustion model for a quiescent chamber diesel engine, SAE Paper 811235, 1981. - 8.
Hiroyasu, H., Kadota, T., and Arai, M. Development and use of a spray combustion modeling to predict diesel engine efficiency and pollutant emissions: Part 1. Combustion modeling, Bulletin of JSME, 26 (1983), pp. 569–575. - 9.
Kaminoto, T., Chang, Y. J., and Kobayashi, H. Rate of heat release and its prediction of a diesel flame in a rapid compression machine, SAE Paper 841076, 1984. - 10.
Ikegami, M., Shioji, M., and Koike, M. A stochastic approach to model the combustion process in direct-injection diesel engines. In: Twentieth International Symposium on Combustion, (1984), pp. 217–224. - 11.
Hotta, Y., Inayoshi, M., Nakakita, K., Fujiwara, K., and Sakata, I. Achieving lower exhaust emissions and better performance in an HSDI diesel engine with multiple injection, SAE Paper 2005-01-0928, 2005. - 12.
Badami, M., Mallamo, F., Millo, F., and Rossi, E. E. Influence of multiple injection strategies on emissions, combustion noise and BSFC of a DI common rail diesel engine, SAE Paper 2002-01-0503, 2002. - 13.
Vanegas, A., Won, H., Felsch, C., et al. Experimental investigation of the effect of multiple injections on pollutant formation in a common-rail DI diesel engine, SAE Paper 2008-01-1911, 2008. - 14.
Inagaki, K., Ueda, M., Mizuta, J., et al. Universal diesel engine simulator (UniDES): 1st report: Phenomenological multi-zone PDF model for predicting the transient behavior of diesel engine combustion, SAE Paper 2008-01-0843, 2008. - 15.
Dulbecco, A. and Lafossas, F. A. A 0D phenomenological approach to model diesel HCCI combustion with multi-injection strategies using probability density functions and detailed tabulated chemistry, SAE Paper 2009-01-0678, 2009. - 16.
Nicolas, B., Christian, C., and Pascal, H. A physical 0D combustion model using tabulated chemistry with presumed probability density function approach for multi-injection diesel engines, SAE Paper 2010-01-1493, 2010. - 17.
Sangras, R., Kwon, O. C., and Faeth, G. M. Self-preserving properties of unsteady round nonbuoyant turbulent starting jets and puffs in still fluids. Journal of Heat Transfer, 124 (2002), pp. 460–469. - 18.
Musculus, M.P.B. Entrainment waves in decelerating transient turbulent jets, Journal of Fluid Mechanics, 638 (2009), pp. 117–140. - 19.
Musculus, M.P.B. and Kattke, K. Entrainment waves in diesel jets, SAE Paper 2009-01-1355, 2009. - 20.
Musculus, M. P. B., Lachaux, T., Pickett, L. M., et al. End-of-injection over-mixing and unburned hydrocarbon emissions in low-temperature-combustion diesel engines, SAE Paper 2007-01-0907, 2007. - 21.
Moon, S., Matsumoto, Y., and Nishida, K. Entrainment evaporation and mixing characteristics of diesel sprays around end-of-injection, SAE Paper 2009-01-0849, 2009. - 22.
Kastengren, A., Powell, C. F., Tilocco, F. Z., et al. End-of-injection behavior of diesel sprays measured with X-ray radiography, Journal of Engineering for Gas Turbines and Power, 134 (2012), 094501 (1–4). - 23.
Yamane, K., Ikegami, M., and Shioji, M. Fuel injection pressure and nozzle orifice diameter in direct-injection diesel engine, Proceedings of International Symposium COMODIA, 94 (1994), pp. 225–230. - 24.
Ishiyama, T., Shioji, M., and Ihara, T. Analysis of ignition processes in a fuel spray using an ignition model including turbulent mixing and reduced chemical kinetics, International Journal of Engine Research, 4 (2003), pp. 155–162. - 25.
Ishiyama, T., Kee, S. S., and Kitamura, Y. Modeling and experiment of NOx formation in DI-PCCI combustion, SAE Paper 2007-01-0194, 2007. - 26.
Horibe, N. and Ishiyama, T. Study on selection of injection conditions for diesel PCCI combustion: Numerical investigation using an ignition-combustion model based on stochastic approach. Transactions of JSME(B), 76 (2010), pp. 699–706, (In Japanese). - 27.
Janicka, J., Kolbe, W., and Kollmann, W. Closure of the transport equation for the probability density function of turbulent scalar fields, Journal of Non-Equilibrium Thermodynamics, 4 (1979), pp. 47–66. - 28.
Wakuri, Y., Fujii, M., Amitani, T., and Tsuneya, R. Studies on the penetration of fuel spray in a diesel engine, Bulletin of JSME, 3 (1960), pp. 123–130. - 29.
Schreiber, M., Sakak, A. S., and Lingens, A. A reduced thermokinetic model for the autoignition of fuels with variable octane ratings, In: Twenty-fifth Symposium. (International) on Combustion, (1994), pp. 933–940. - 30.
Woschni, G. Universally applicable equation for the instantaneous heat transfer coefficient in the internal combustion engine, SAE Paper 670931, 1967. - 31.
Long Liu, Naoto Horibe, and Takuji Ishiyama, Combustion Modeling for a Diesel Engine with Multi-Stage Injection Using a Stochastic Combustion Model, Procs. IMech. E., Part D: J. Auto. Eng., 228 (5), pp. 516–532, 2014. - 32.
Hiroyasu, H., Kadota, T., Arai, M. Supplementary Comments: Fuel spray characterization in diesel engines. In: James N. Mattavi and Charles A. Amann (eds) Combustion modeling in reciprocating engines. Plenum Press., New York, (1980), pp. 369–408. - 33.
Kau, C. J., Heap, M. P., Tyson, T. J., and Wilson, R. P. The prediction of nitric oxide formation in a direct injection diesel engine, Proceedings of 16th Symposium (Internatinal) on Combustion, (1977), pp.337–350. - 34.
Long Liu, Naoto Horibe, Issei Tamura, Tatsuya Komizo, and Takuji Ishiyama, Study on NOx and Soot Emissions in Diesel Engines with Multi-Stage Injection Based on a Stochastic Combustion Model, International Journal of Automotive Engineering, Vol.5(1), pp. 7–13, 2014. - 35.
Liu, L., Horibe, N., and Ishiyama, T. Combustion modeling for a diesel engine with multi-stage injection using a stochastic combustion model, Proceedings of the Institution of Mechanical Engineers, Part D: Journal of Automobile Engineering, 228 (5), (2014), pp. 516–532. - 36.
Hiroyasu, H., Kadota, T., and Arai, M. Supplementary comments: Fuel spray characterization in diesel engines. In: Mattavi, J.N., and Amann, C.A. (eds), Combustion Modeling in Reciprocating Engines, Plenum Press., New York, (1980), pp. 369–408. - 37.
Kook, S., Park, S., and Bae, C. Influence of early fuel injection timings on premixing and combustion in a diesel engine, Energy & Fuel, 22 (2008), pp. 331–337. - 38.
Hiroyasu, H., Kadota, T., and Arai, M. Development and use of a spray combustion modeling to predict diesel engine efficiency and pollutant emissions (Part 1 Combustion Modeling), Bulletin of the JSME, 26(1983), pp.569-575. - 39.
Bruneaux, G. Mixing process in high pressure diesel jets by normalized laser induced exciplex fluorescence. Part II: Wall impinging versus free jet, SAE Paper 2005-01-2097, 2005. - 40.
Bruneaux, G. and Causse, M. Air entrainment in diesel-like gas jet by simultaneous flow velocity and fuel concentration measurements, comparison of free and wall impinging jet configurations, SAE Paper 2011-01-1828, 2011. - 41.
Lavoie, G. A., Heywood, J. B., and Keck, J.C. Experimental and theoretical investigation of nitric oxide formation in internal combustion engines, Combustion Science and Technology, 1 (1970), pp. 313–326. - 42.
Moss, J. B., Stewart, C. D., and Young, K. J. Modeling soot formation and burnout in a high-temperature laminar diffusion flame burning under oxygen-enriched conditions, Combustion and Flame, 101 (1995), pp. 491–500. - 43.
Nagle, J. and Strickland-Constable, R. F. Oxidation of carbon between 1000–2000°C, Proceedings of 5th Carbon Conference, 1 (1962), pp. 154–164. - 44.
Neoh, K.G. Soot burnout in flames, Ph.D. Thesis, Massachusetts Institute of Technology, (1980).