Investigation on an Injection Strategy Optimization for Diesel Engines Using a One-Dimensional Spray Model

Common rail systems have been widely used in diesel engines due to the stricter emission regulations. The advances in injector technology and ultrahigh injection pressure greatly promote the development of multiple-injection strategy, leading to the shorter injection duration and more variable injection rate shape, which makes the mixing process more significant for the formation of pollutant emission. In order to study the mixing process of diesel sprays under variable injection rate shapes and find the optimized injection strategy, a one-dimensional spray model was modified in this paper. The model was validated by the measured spray penetrations based on shadowgraphy experiments with the varying injection rate. The simulations were performed with five injection rate shapes, triangle, ramping-up, ramping-down, rectangle and trapezoid. Their spray penetrations, entrainment rates and equivalence ratios along spray axial distance are compared. The potentials of multiple-injection and gas-jet after end-of-injection (EOI) to improve mixing process and emission reduction are discussed finally. The results indicated that ramping-up injection rate obtains the highest entrainment rate after EOI, and it needs 1.5 times of injection duration for the entrainment wave to arrive at the spray tip. For the other four injection rates, the sprays can be treated as a steady-like state, needing twice of injection duration from EOI to the time the entrainment wave reaches the spray tip. The multiple-injection with proper injection rate shape enhanced the entrainment rate, and the gas-jet after EOI affected the mixture distribution and entrainment rate in spray tail under ramping-down injection rate.


Introduction
With the implementation of more stringent emission regulations, such as US Tier 3 and post-Euro 6, diesel engines are facing unprecedented challenges to reduce their adverse environmental impact in terms of pollutant and greenhouse emissions. To cope with the emission problem and improve the performance of diesel engines, the common rail technology has flourished around the world. The advanced fuel injection technologies in a common rail system, such as the ultrahigh injection pressure [1][2][3][4] and the multiple-injection strategy [5][6][7] as well as the high control accuracy [8], have greatly improved the performance of diesel engines. Under the action of ultrahigh injection pressures, the fuel can be injected into the cylinder in a short time, while the multiple-injection strategy divides the whole injection duration into many parts, creating shorter injection durations. Short injection duration often causes EOI much earlier than ignition under certain conditions, such as premixed charge combustion ignition (PCCI) and low temperature combustion (LTC) cases, resulting in the more non-combusting diesel sprays in a constant volume chamber. Then the mixing process was studied by calculating and analyzing the spray penetrations, entrainment rates and equivalence ratio distributed along the spray axial distance. After that, the improved mixing characteristics of the multiple-injection strategy that consists of ramping-down and ramping-up injection rates were discussed in depth. Finally, this paper also investigated the effect of the gas-jet after EOI on mixture distribution and emission reduction under the ramping-down injection rate shape.

Spray Model
After EOI, the deceleration state of turbulent sprays has been found to travel downstream from the injector nozzle, meanwhile, the entrainment rate of the passing region is greatly enhanced. The propagation of the increased entrainment region is commonly known as the "entrainment wave".
The original one-dimensional discrete spray model, developed by Musculus et al. [23], could not only predict the entrainment wave in diesel spray, but also calculate the spray penetration. The computation results showed that the spray penetration after the entrainment wave arrived at the spray tip was smaller when compared with that of a steady jet. Moreover, under the action of the entrainment wave, the dependence of spray penetration on time transited from the square-root to the fourth-root. This result was validated by several experimental data. It indicates that this model has a significant advantage in mimicking the diesel spray propagation after EOI. However, the original model mainly focused on the constant injection rate case and didn't show any application for the varying injection rate. To overcome this limitation, the original model is improved in this paper to introduce the effect of varying injection rate on the spray propagation by calculating the velocity at the nozzle exit under varying fuel mass flow rate. Because based on the one-dimensional discrete spray model that proposed by Musculus et al. [23], the improved one has the same assumptions as the former. Furthermore, taking the coking and cavitation phenomenon as well as the liquid diesel spray contraction into account, the nozzle diameter is used with a discharge coefficient in computation. Figure 1 shows the schematic diagram of one-dimensional discrete diesel spray model. The fuel mass and axial momentum in each control volume are solved by the transport Equations (1) and (2). M out (2) where m f is fuel mass and M is axial momentum in a control volume. Overdot means integral flux crossing the inlet or outlet of each control volume. Where m f is fuel mass and M is axial momentum in a control volume. Using an upwind differencing scheme, the discretized equations for fuel mass and momentum are: where = X f,i (t), = u i (t) are the fuel volume fraction and fluid velocity averaged over the spray cross section at i − th control volume respectively, ρ f is the fuel density, β i is the coefficient which determines the radial velocity profile at i − th control volume, A i is the cross-sectional area at the downstream surface where ∆Z is the length of a control volume, θ is the spreading angle of diesel spray, ρ a (t) is the density of surrounding gas. A i is evaluated at the axial midpoint of each control volume in the spray model of Musculus et al. [23], however, A i defined at the downstream surface area of a control volume gives almost the same results when a sufficiently small ∆Z is used.

Varying Injection Rate Calculation
Boundary conditions for solving these equations above are according to Equations (9) and (10): where the subscript "0" means the inlet surface of the first control volume, i.e., the nozzle exit. u f (t) is a cross-section averaged fuel velocity at the nozzle exit and is calculated by Equation (11): where . m f,0 (t) is the injection mass rate, µ is the flow coefficient of the nozzle orifice, d n is the geometric orifice diameter. m f,0 (t) is assumed to have a trapezoidal pattern, as shown in Figure 2, and it can be calculated according to Equation (12).
where the injection mass rate peak . m f,max is obtained by Equation (13): .
is the injection mass rate, μ is the flow coefficient of the nozzle orifice, d n is the is assumed to have a trapezoidal pattern, as shown in Figure 2, and it can be calculated according to Equation (12).
where the injection mass rate peak m  f,max is obtained by Equation (13):

Validation of the Spray Model with Short Injection Duration
For model validation, the diesel spray tip penetration even after EOI with short injection duration was calculated and compared with the shadowgraph data of non-combusting diesel sprays in a constant volume chamber. The experimental setup and measuring equipment are the same as that in our previous study [11]. The constant volume chamber is in cylindrical shape with a volume of 150 cm 3 . Its diameter and depth are 80 mm and 30 mm, respectively. In the shadowgraph imaging system, a Xeon lamp is used as the light source. The focal lengths of the two same concave mirrors used are 1910 mm. The high-speed camera has a shooting speed of 50,000 fps and the exposure time was 5 ms. With the aid of a common-rail injection system, an advanced injector with a discharge coefficient of 0.8 is employed to inject diesel fuel into the chamber. The operating conditions are listed in Table 1, and the calculation was also performed based on these conditions.

Validation of the Spray Model with Short Injection Duration
For model validation, the diesel spray tip penetration even after EOI with short injection duration was calculated and compared with the shadowgraph data of non-combusting diesel sprays in a constant volume chamber. The experimental setup and measuring equipment are the same as that in our previous study [11]. The constant volume chamber is in cylindrical shape with a volume of 150 cm 3 . Its diameter and depth are 80 mm and 30 mm, respectively. In the shadowgraph imaging system, a Xeon lamp is used as the light source. The focal lengths of the two same concave mirrors used are 1910 mm. The high-speed camera has a shooting speed of 50,000 fps and the exposure time was 5 ms. With the aid of a common-rail injection system, an advanced injector with a discharge coefficient of 0.8 is employed to inject diesel fuel into the chamber. The operating conditions are listed in Table 1, and the calculation was also performed based on these conditions.  Figure 3 shows the comparison of spray tip penetrations between the simulation result and the experiment data from shadowgraph with the injection quantity of 0.417 mg and 0.833 mg for single nozzle hole under 180 MPa injection pressure, respectively. The spray tip penetrations and time are expressed on a logarithmic scale. As the injection rates in Figure 3 show, the actual injection rates have a shape like a triangle, therefore the calculations were based on the triangle injection rate and the injection durations were kept the same with that of the measurements. The results reveal that the simulation results are slightly smaller than the experiment data, which can be ascribed to the slightly lower measurement accuracy of injection rate and spray tip penetration under a very small amount of fuel injected, especially in the initial injection period. The simplification of the triangle injection rate for simulation further increases such a difference. However, the simulation results of spray tip penetration show similar tendencies to the measured data. Furthermore, the spray penetrations are able to capture the relation of square-root dependence on time around the moment of twice the injection duration. Then the spray penetrations change to the relation of fourth-root dependence on time. This phenomenon is coherent with that affected by the entrainment wave [24]. Thus, it is reasonable to apply this spray model to analyze the diesel spray mixing process.   Figure 3 shows the comparison of spray tip penetrations between the simulation result and the experiment data from shadowgraph with the injection quantity of 0.417 mg and 0.833 mg for single nozzle hole under 180 MPa injection pressure, respectively. The spray tip penetrations and time are expressed on a logarithmic scale. As the injection rates in Figure 3 show, the actual injection rates have a shape like a triangle, therefore the calculations were based on the triangle injection rate and the injection durations were kept the same with that of the measurements. The results reveal that the simulation results are slightly smaller than the experiment data, which can be ascribed to the slightly lower measurement accuracy of injection rate and spray tip penetration under a very small amount of fuel injected, especially in the initial injection period. The simplification of the triangle injection rate for simulation further increases such a difference. However, the simulation results of spray tip penetration show similar tendencies to the measured data. Furthermore, the spray penetrations are able to capture the relation of square-root dependence on time around the moment of twice the injection duration. Then the spray penetrations change to the relation of fourth-root dependence on time. This phenomenon is coherent with that affected by the entrainment wave [24]. Thus, it is reasonable to apply this spray model to analyze the diesel spray mixing process.

Calculation of Spray Tip Penetrations Under Varying Injection Rate
As has been found in a previous study [19], the ramping-up and ramping-down injection rate shapes both have a significant influence on FPEG. Due to the similarity between FPEG and the conventional diesel engine, these two injection rate shapes must also play an important role in the latter. However, there is little relative research about this. Inversely, as the common used and typical injection rate shapes in diesel engines, scholars have devoted a great deal of effort to the triangle and rectangle as well as trapezoid shapes [9][10][11][12][13]. In order to compare the potential of varying injection rate shapes to lower the pollutant emissions of diesel engines more objectively and comprehensively and find the optimized injection strategy, five injection rate shapes (including triangle, ramping-up,

Calculation of Spray Tip Penetrations Under Varying Injection Rate
As has been found in a previous study [19], the ramping-up and ramping-down injection rate shapes both have a significant influence on FPEG. Due to the similarity between FPEG and the conventional diesel engine, these two injection rate shapes must also play an important role in the latter. However, there is little relative research about this. Inversely, as the common used and typical injection rate shapes in diesel engines, scholars have devoted a great deal of effort to the triangle and rectangle as well as trapezoid shapes [9][10][11][12][13]. In order to compare the potential of varying injection rate shapes to lower the pollutant emissions of diesel engines more objectively and comprehensively and find the optimized injection strategy, five injection rate shapes (including triangle, ramping-up, ramping-down, rectangle and trapezoid) were selected in this paper, as shown in Figure 4. To simulate the diesel spray propagation under the conditions extant in engines' cylinders, the ambient gas conditions are set as the actual in-cylinder pressure and temperature, so the calculation conditions are the same as given in Table 1, and the injection rates have been shown in a previous study [11]. ramping-down, rectangle and trapezoid) were selected in this paper, as shown in Figure 4. To simulate the diesel spray propagation under the conditions extant in engines' cylinders, the ambient gas conditions are set as the actual in-cylinder pressure and temperature, so the calculation conditions are the same as given in Table 1, and the injection rates have been shown in a previous study [11].  Figure 5 shows the spray tip penetrations with the injection quantity of 0.417 mg and 0.833 mg for a single nozzle hole, respectively. Both Figure 5a,b reveal the same tendency that the initial spray tip penetration depends on the injection rate. A higher injection rate leads to a larger spray tip penetration during the initial injection. A higher maximum injection rate finally leads to a larger spray tip penetration, which is different from the previous study [19] and can be seen from two aspects. On the one hand, the  Figure 5 shows the spray tip penetrations with the injection quantity of 0.417 mg and 0.833 mg for a single nozzle hole, respectively. Both Figure 5a,b reveal the same tendency that the initial spray tip penetration depends on the injection rate. A higher injection rate leads to a larger spray tip penetration during the initial injection. ramping-down, rectangle and trapezoid) were selected in this paper, as shown in Figure 4. To simulate the diesel spray propagation under the conditions extant in engines' cylinders, the ambient gas conditions are set as the actual in-cylinder pressure and temperature, so the calculation conditions are the same as given in Table 1, and the injection rates have been shown in a previous study [11].  Figure 5 shows the spray tip penetrations with the injection quantity of 0.417 mg and 0.833 mg for a single nozzle hole, respectively. Both Figure 5a,b reveal the same tendency that the initial spray tip penetration depends on the injection rate. A higher injection rate leads to a larger spray tip penetration during the initial injection. A higher maximum injection rate finally leads to a larger spray tip penetration, which is different from the previous study [19] and can be seen from two aspects. On the one hand, the A higher maximum injection rate finally leads to a larger spray tip penetration, which is different from the previous study [19] and can be seen from two aspects. On the one hand, the triangle, ramping-up and ramping-down injection rate shapes have the same maximum injection rate, although their spray tip penetrations evolve with different growth rates in the early injection duration, they all gradually tend to the same penetration value and slope over time. On the other hand, compared with the rectangle and trapezoid injection rate shapes, the triangle injection rate reaches the larger spray tip penetration finally due to its higher maximum injection rate though its injection rate is smaller than that of the other two during the initial injection. However, in spite of the fact that the triangle and ramping-up as well as ramping-down injection rate shapes reach the final same spray penetration, the ramping-down injection rate has been keeping the greater spray tip penetration until the differences between them can be ignored, meaning that its overall spray gets more time to mix with air. Therefore, the ramping-down injection rate shape could obtain the better fuel-air mixing if the ignition delay is long enough.

Analysis of Entrainment Rates Under Varying Injection Rate
One of the indicators that allow for the insightful analysis on mixing process of diesel spray is the entrainment rate. Figure 6 shows the entrainment rates with varying injection rate shapes. Comparing the entrainment rates for the same injection rate shape under different fuel injection quantity, the larger fuel injection quantity case ( Figure 6b) visibly gets much higher entrainment rate, with the variation tendency extremely similar to the smaller fuel injection quantity case ( Figure 6a). For the same fuel injection quantity, it can be found that the higher injection rate causes the higher entrainment rate during the initial injection. Each injection rate shape has the special entrainment rate pattern and the entrainment rate comes to the peak soon after EOI. The entrainment rate peaks (from large to small) are in the order: ramping-up, triangle, trapezoid, ramping-down and rectangle. It is the entrainment wave after EOI that should be responsible for such difference. triangle, ramping-up and ramping-down injection rate shapes have the same maximum injection rate, although their spray tip penetrations evolve with different growth rates in the early injection duration, they all gradually tend to the same penetration value and slope over time. On the other hand, compared with the rectangle and trapezoid injection rate shapes, the triangle injection rate reaches the larger spray tip penetration finally due to its higher maximum injection rate though its injection rate is smaller than that of the other two during the initial injection. However, in spite of the fact that the triangle and ramping-up as well as ramping-down injection rate shapes reach the final same spray penetration, the ramping-down injection rate has been keeping the greater spray tip penetration until the differences between them can be ignored, meaning that its overall spray gets more time to mix with air. Therefore, the ramping-down injection rate shape could obtain the better fuel-air mixing if the ignition delay is long enough.

Analysis of Entrainment Rates Under Varying Injection Rate
One of the indicators that allow for the insightful analysis on mixing process of diesel spray is the entrainment rate. Figure 6 shows the entrainment rates with varying injection rate shapes. Comparing the entrainment rates for the same injection rate shape under different fuel injection quantity, the larger fuel injection quantity case (Figure 6b) visibly gets much higher entrainment rate, with the variation tendency extremely similar to the smaller fuel injection quantity case ( Figure  6a). For the same fuel injection quantity, it can be found that the higher injection rate causes the higher entrainment rate during the initial injection. Each injection rate shape has the special entrainment rate pattern and the entrainment rate comes to the peak soon after EOI. The entrainment rate peaks (from large to small) are in the order: ramping-up, triangle, trapezoid, ramping-down and rectangle. It is the entrainment wave after EOI that should be responsible for such difference.   As discussed in the previous study [2], the velocity at the nozzle exit is high enough to affect the spray when the injection rate just begins to decrease for the ramping-down, triangle and trapezoid injection rates. Once the velocity at the nozzle exit is too low to affect the spray, the real time of the EOI (tEOI) can be defined. tEOI is the start of entrainment wave and the velocity at the nozzle exit at tEOI (UEOI) determines the spray penetration level. In fact, the peak of entrainment rate occurs when the entrainment wave arrives at the spray tip. Thus, the entrainment rate peak is decided by tEOI and UEOI. To calculate tEOI and UEOI under varying injection rate shapes, the Equation (14) proposed by Liu et al. [11] is used in this paper: As discussed in the previous study [2], the velocity at the nozzle exit is high enough to affect the spray when the injection rate just begins to decrease for the ramping-down, triangle and trapezoid injection rates. Once the velocity at the nozzle exit is too low to affect the spray, the real time of the EOI (t EOI ) can be defined. t EOI is the start of entrainment wave and the velocity at the nozzle exit at t EOI (U EOI ) determines the spray penetration level. In fact, the peak of entrainment rate occurs when the entrainment wave arrives at the spray tip. Thus, the entrainment rate peak is decided by t EOI and U EOI . To calculate t EOI and U EOI under varying injection rate shapes, the Equation (14) proposed by Liu et al. [11] is used in this paper: where R IV is a dimensionless ratio and is selected to be 0.65 in this paper, t decay is the period that injection velocity decays, t inj is the injection duration, u max is the maximum injection velocity. The calculation results of t EOI and U EOI for the injection quantity of 0.417 mg and 0.833 mg are shown in Tables 2 and 3 respectively. t peak denotes the time at the entrainment rate peak. A same law can be discovered from Tables 2 and 3, that is, t peak in the other four injection rate shapes, except the ramping-up, is almost twice that of t EOI , which is because the entrainment wave travels downstream at twice the initial spray penetration rate. As reference [24] demonstrated, the typical phenomenon is revealed in the case of rectangle injection rate, the diesel jet with rectangle injection rate can be treated as a steady jet, and the entrainment wave front arrives at the jet tip at twice of injection duration. Figure 7 shows the mean velocity over the cross-sectional area against the axial position in the diesel spray under five injection rate shapes with time varying. As shown in Figure 7d (rectangle injection rate), the same axial position has the same velocity during injection due to the steady state. There is an envelope velocity curve for the steady spray, and the velocity at the position that belong to the spray steady part after EOI is on the envelope velocity curve. After the entrainment wave front arriving at the spray tip, the total spray transfers to the decelerating state, which causes the velocities at any axial position are lower than that of steady state.
For the trapezoid injection rate, the constant injection rate takes main injection process, so that the spray can be treated as the steady state. According to the previous study [11], the twice of t EOI is close to the t peak . As shown in Figure 7e, the envelope velocity curve starts from the steady part and it extends to the t peak . After t peak , the spray gets into the decelerating state. In Figure 7c, although the spray is not in the steady state due to the ramping-down injection rate, the spray tip velocities at different time before t peak are on the same envelope velocity curve. After t peak , the spray tip velocities depart from the envelope velocity curve, and the values are lower than that on the envelope curve. It might be reasonable to use the twice of U EOI to represent the average entrainment wave speed, and the U EOI can be treated as the average effective fuel injection velocity that considers the fuel penetration as the spray penetration. Then the t peak is similar with twice of t EOI even for the unsteady spray with ramping-down injection rate.
However, because the momentum flux at each axial position of spray needs a certain response time to the change of the injected velocity, the velocities at some positions still increase even the injection velocity turns to decrease or stop, which can be found in Figure 7a,b. In Figure 7a, the envelope velocity curve can cover the spray tip velocity up to t peak from 0.16 ms. That means the injection velocity increase effect has been reached spray tip, the decrease of injection velocity effect mainly control the spray. Meanwhile, the t peak is almost equal to the twice of t EOI , which is indicated the effectiveness of the average entrainment wave and spray penetration speed concept. (e) (f) Regarding the ramping-up injection rate, there is no steady-like state as shown in Figure 7b, because the injection rate is reduced to zero sharply. At the tEOI, the UEOI is the maximum injection Regarding the ramping-up injection rate, there is no steady-like state as shown in Figure 7b, because the injection rate is reduced to zero sharply. At the t EOI , the U EOI is the maximum injection velocity, the average entrainment wave speed can be calculated by twice of U EOI . However, the increase of the injection velocity effect has not traveled to spray tip, and the U EOI cannot represent the average spray penetration speed. Similar to the assumption of Liu et al. [11], a same parameter R IV (0.65) is selected in this paper as a factor of U EOI to calculate the average fuel penetration speed, which represents the average spray penetration speed and one third of average entrainment wave speed. As shown in Equation (15), t 1 is the duration from the t EOI to the time of entrainment wave reaching the spray tip: Equation (16) can be simplified as: then the following expression is obtained: Thus, the entrainment wave arrives at the spray tip at t = 1.5 t EOI , which is consistent with the calculation result of t peak and t EOI under the ramping-up injection rate, t peak (0.25 ms) is about 1.5 times (0.2625 ms) of t EOI for the injection quantity of 0.417 mg, while t peak (0.42 ms) is also about 1.5 times (0.4455 ms) of t EOI for the injection quantity of 0.833 mg. In addition, as shown in Figure 7b, the entrainment wave takes 0.075 ms to arrive at the spray tip with the spray penetration of 15.8 mm at 0.25 ms. The average entrainment wave traveling speed is 210.67 mm/ms, and the average spray penetration rate is about 63.2 mm/ms. The average entrainment wave traveling speed is 3.33 times of the average spray penetration rate, which is similar with the result calculated using the U EOI .
Moreover, due to the largest U EOI, the spray of ramping-up injection rate penetrates fast after EOI, compared with the triangle and ramping-down injection rates ( Figure 5), so the terminal spray penetration is a combined effect result of injection rate before EOI and U EOI . The triangle, ramping-up and ramping-down injection rate shapes tend to the same spray penetration finally just because an opportune tradeoff is reached.
As shown in Figure 7, there are some similarities between five injection rate shapes: the mean axial velocity decreases along the axial distance before EOI and increases along the axial distance after EOI, but the velocity gradients are different. When fuel injection supply stops after EOI, the mass conservation and continuity theorem require a fast gas entrainment to compensate for the decreasing fuel mass flux. Therefore, the fuel mass flux, or speaking directly, the velocity gradient determines the intensity of entrainment wave. As mentioned above, the entrainment rate peak occurs when the entrainment wave arrives at the spray tip and t peak is decided by t EOI as well as the speed of entrainment wave. In order to explain the difference of entrainment rate peak for five injection rate shapes, the velocity gradient at t peak should be taken into account. Shown in Figure 7f, it is obvious that five injection rate shapes have different velocity gradients and velocity magnitude at t peak . Because the entire spray has been enhanced by the entrainment wave at this moment, associated with the cone structure of spray, the entrainment effect of downstream spray accounts for a larger proportion of the entire spray due to its larger volume. Thus, the ramping-up injection rate has the highest velocity gradient and the fast mean axial velocity in downstream, which leads to the biggest entrainment rate peak in Figure 6, followed by the triangle injection rate. The trapezoid injection rate gets the slightly larger entrainment rate peak than the ramping-down on account of its higher velocity gradient, too. While for the rectangle injection rate, its velocity gradient and velocity magnitude are both the lowest, resulting in a much weaker entrainment effect than the other four injection rate shapes, so its entrainment rate peak is the lowest.

Analysis of Mixing Process Based on Equivalence Ratios
Spray tip penetration and entrainment rate can estimate the diesel spray mixing process in the macroscopic view. However, the microscopic mixing process is also a very significant factor for diesel spray combustion. To analyze the microscopic mixing process, the equivalence ratio along the spray axial position with varying time was calculated using the spray model. Figure 8 shows the equivalence ratio distributions under five injection rate shapes. Figure 8a,c,e show the similar equivalence ratio level after EOI in the cases of triangle, ramping-down as well as trapezoid injection rates, the equivalence ratio is almost evenly distributed along the spray axial distance after EOI.
Especially at the time of 0.36 ms, the equivalence ratio is around 1.0. The fuel could be burned completely if the ignition occurred at that time, as well as there is no enough O 2 remained for NO x formation. Therefore, the fuel consumption and NO x emission can decrease simultaneously. But compared with the triangle and trapezoid injection rates, the ramping-down is able to reach much longer spray penetration before EOI, while the downstream spray has lower equivalence ratio, indicating a stronger mixing process. So the ramping-down injection rate has an advantage over the other two in creating homogeneous mixture. However, due to the low velocity near the nozzle exit, the fuel-rich mixture formation near the nozzle exit, it usually leads to the high soot formation. Figure 8b,d show the similar results of equivalence ratio along the spray axial distance with varying time in the cases of ramping-up and rectangle injection rates, the mixture near the nozzle exit quickly gets leaner after EOI, and the equivalence ratio increases with the axial distance increasing. It means that the both two injection rates with longer ignition delay can decrease the soot formation near the nozzle exit on certain conditions. It can be observed that the ramping-up injection rate obtains much leaner mixture than that of rectangle injection rate after EOI. This phenomenon indicates that the ramping-up injection rate is a better choice to realize LTC, because the ignition is occurred after EOI (with the overlong ignition delay not considered). In summary, the spray characteristics among the penetration and entrainment rate as well as equivalence ratio under varying injection rate shapes can be synoptically shown in Table 4, where I to V denote the maximum to the minimum of each parameter. From the practical point of view, these spray characteristics can provide some important guidance and reference. For example, the long penetration is always expected in the low-speed diesel engine with big bore and long stroke, so the ramping-down injection rate is an optimal choice. Instead, the spray-wall impingement process is tried to be avoided sometimes, hence a relatively short penetration is needed and the rectangle injection rate should be taken into account first. In some advanced combustion modes, such as PCCI and LTC, the ramping-up injection rate can be used to greatly promote the mixing of fuel and the cold ambient air or exhaust gas, then quickly decrease the temperature of mixture as well as the formation of NO x . As for homogeneous charge compression ignition (HCCI), the uniform charge in cylinder is indispensable, so that all points in diesel spray are ignited at the same time and the heat is evenly distributed throughout the cylinder. In this way, the triangle, ramping-down as well as trapezoid injection rates are the theoretically feasible alternatives because the equivalence ratio is almost evenly distributed along the spray axial distance after EOI in these cases.  Triangle  II  II  II  Ramping-up  III  I  V  Ramping-down  I  IV  I  Rectangle  V  V  IV  Trapezoid  IV  III  III macroscopic view. However, the microscopic mixing process is also a very significant factor for diesel spray combustion. To analyze the microscopic mixing process, the equivalence ratio along the spray axial position with varying time was calculated using the spray model. Figure 8 shows the equivalence ratio distributions under five injection rate shapes. Figure 8a,c,e show the similar equivalence ratio level after EOI in the cases of triangle, ramping-down as well as trapezoid injection rates, the equivalence ratio is almost evenly distributed along the spray axial distance after EOI.  (e)

Effect of the Multiple-Injection
According to the information in Table 4, the ramping-down injection rate could reach the longest penetration but also the highest equivalence ratio near the nozzle, while the ramping-up injection rate obtains the largest entrainment rate peak and the lowest equivalence ratio near the nozzle. So it is a good idea to develop a new multiple-injection strategy that both contains the ramping-down and ramping-up injection rates. Figure 9 shows the comparisons of spray tip penetration and entrainment rate between the ramping-down single injection and the multiple-injection. In Figure 9a, the spray tip penetration of the multiple-injection is equal to that of the ramping-down single injection before t = 0.1 ms. After that, the spray tip penetration of the multiple-injection becomes smaller than that of the ramping-down single injection due to the sharply decrease of injection rate. As shown in Figure 9b, although the entrainment rate of the multiple-injection is smaller than that of the ramping-down single injection from t = 0.05 ms to t = 0.3 ms, it increases rapidly over the latter after t = 0.3 ms because the ramping-up sub-injection begins to play a leading role in the entrainment process. multiple-injection is smaller than that of the ramping-down single injection from t = 0.05 ms to t = 0.3 ms, it increases rapidly over the latter after t = 0.3 ms because the ramping-up sub-injection begins to play a leading role in the entrainment process.  Figure 10 shows the comparison of entrainment rate after EOI between the multiple-injection and the ramping-down single injection. The spray head of the multiple-injection is a little shorter than that of the ramping-down single injection, which is in agreement with Figure 9a. When it is at t = 0.18 ms and t = 0.24 ms, since the sub-spray generated by the ramping-up sub-injection has not overtaken the sub-spray generated by the ramping-down sub-injection, and much fuel accumulates in the former spray head, so the high-concentration fuel causes little gas to be entrained. Due to the slower injection velocity, the entrainment rate of the ramping-down sub-injection in the multiple-injection is smaller than that of the ramping-down single injection even though they have the similar injection rate shape. However, the upstream sub-spray generated by the ramping-up sub-injection of the multiple-injection has much higher entrainment rate than that of the ramping-down single injection. After t = 0.24 ms, the upstream sub-spray overtakes the downstream sub-spray. Because the entrainment wave of the ramping-up sub-injection plays a dominant role in the whole spray, the entrainment rate of the multiple-injection is much higher than that of the ramping-down single injection, indicating that the multiple-injection can get better mixing effect.   Figure 10 shows the comparison of entrainment rate after EOI between the multiple-injection and the ramping-down single injection. The spray head of the multiple-injection is a little shorter than that of the ramping-down single injection, which is in agreement with Figure 9a. When it is at t = 0.18 ms and t = 0.24 ms, since the sub-spray generated by the ramping-up sub-injection has not overtaken the sub-spray generated by the ramping-down sub-injection, and much fuel accumulates in the former spray head, so the high-concentration fuel causes little gas to be entrained. Due to the slower injection velocity, the entrainment rate of the ramping-down sub-injection in the multiple-injection is smaller than that of the ramping-down single injection even though they have the similar injection rate shape. However, the upstream sub-spray generated by the ramping-up sub-injection of the multiple-injection has much higher entrainment rate than that of the ramping-down single injection. After t = 0.24 ms, the upstream sub-spray overtakes the downstream sub-spray. Because the entrainment wave of the ramping-up sub-injection plays a dominant role in the whole spray, the entrainment rate of the multiple-injection is much higher than that of the ramping-down single injection, indicating that the multiple-injection can get better mixing effect. Figure 11 shows the comparisons of equivalence ratio after EOI between the multiple-injection and the ramping-down single injection. In the multiple-injection case, the equivalence ratio after EOI remains the characteristic of ramping-up injection rate shape, resulting from the major influence of the ramping-up sub-injection. The variation of equivalence ratio of the multiple-injection is seen to be quite different from the ramping-down single injection. The equivalence ratio of the multiple-injection is much lower near the nozzle exit and increases gradually along the axial distance. In addition, it is obvious that the spray tail of the multiple-injection moves downstream quickly, while the ramping-down single injection consistently keeps the high equivalence ratio near the nozzle exit. Therefore, the multiple-injection is beneficial to improve the working environment of the nozzle and decrease the soot emission. the similar injection rate shape. However, the upstream sub-spray generated by the ramping-up sub-injection of the multiple-injection has much higher entrainment rate than that of the ramping-down single injection. After t = 0.24 ms, the upstream sub-spray overtakes the downstream sub-spray. Because the entrainment wave of the ramping-up sub-injection plays a dominant role in the whole spray, the entrainment rate of the multiple-injection is much higher than that of the ramping-down single injection, indicating that the multiple-injection can get better mixing effect.  Figure 10. Comparisons of entrainment rate after EOI between the ramping-down single injection and the multiple-injection (solid lines: ramping-down single injection; dotted lines: multiple-injection).   Figure 11 shows the comparisons of equivalence ratio after EOI between the multiple-injection and the ramping-down single injection. In the multiple-injection case, the equivalence ratio after EOI remains the characteristic of ramping-up injection rate shape, resulting from the major influence of the ramping-up sub-injection. The variation of equivalence ratio of the multiple-injection is seen to be quite different from the ramping-down single injection. The equivalence ratio of the multiple-injection is much lower near the nozzle exit and increases gradually along the axial distance. In addition, it is obvious that the spray tail of the multiple-injection moves downstream quickly, while the ramping-down single injection consistently keeps the high equivalence ratio near the nozzle exit. Therefore, the multiple-injection is beneficial to improve the working environment of the nozzle and decrease the soot emission.

Effect of the Gas-Jet After EOI
However, incomplete combustion will occur in the fuel-lean zone if the equivalence ratio is too low, so the originality of combining the ramping-down and ramping-up injection rates to optimize the working conditions of the nozzle and emission performance still has some limitations. To overcome the problem, another improvement approach is presented in this paper.
The main idea of the new method is to use the gas-jet after EOI to blow the mixture near the nozzle exit downstream under the ramping-down injection rate shape, as shown in Figure 12. In this way, the stagnant fuel near the nozzle exit due to the influence of entrainment wave could be carried by the gas-jet to join in the combustion of downstream spray, which can not only improve the working conditions of the nozzle, but also decrease the soot emission.

Effect of the Gas-Jet After EOI
However, incomplete combustion will occur in the fuel-lean zone if the equivalence ratio is too low, so the originality of combining the ramping-down and ramping-up injection rates to optimize the working conditions of the nozzle and emission performance still has some limitations. To overcome the problem, another improvement approach is presented in this paper.
The main idea of the new method is to use the gas-jet after EOI to blow the mixture near the nozzle exit downstream under the ramping-down injection rate shape, as shown in Figure 12. In this way, the stagnant fuel near the nozzle exit due to the influence of entrainment wave could be carried by the gas-jet to join in the combustion of downstream spray, which can not only improve the working conditions of the nozzle, but also decrease the soot emission.
This paper studied the influence of gas jet with two different velocities on the spray evolution, the selected gas-jet velocities (100 m/s, 200 m/s) are both lower than the local speed of sound (648.89 m/s), so the influence of shock wave is inexistent. The calculation conditions for gas-jet are listed in Table 5. The comparisons of spray tip penetration and entrainment rate between the ramping-down injection without gas-jet and the ramping-down injection with gas-jet are shown in Figure 13. It seems that whether there is gas-jet after EOI or not has no effect on the spray tip penetration and entrainment rate, which is related to the limited effect range of gas-jet. overcome the problem, another improvement approach is presented in this paper.
The main idea of the new method is to use the gas-jet after EOI to blow the mixture near the nozzle exit downstream under the ramping-down injection rate shape, as shown in Figure 12. In this way, the stagnant fuel near the nozzle exit due to the influence of entrainment wave could be carried by the gas-jet to join in the combustion of downstream spray, which can not only improve the working conditions of the nozzle, but also decrease the soot emission.    Table 5. The comparisons of spray tip penetration and entrainment rate between the ramping-down injection without gas-jet and the ramping-down injection with gas-jet are shown in Figure 13. It seems that whether there is gas-jet after EOI or not has no effect on the spray tip penetration and entrainment rate, which is related to the limited effect range of gas-jet.   Figure 14 shows the variations of equivalence ratio of the ramping-down injection without gas-jet and the ramping-down injection with gas-jet. It is obvious that the gas-jet after EOI has a great impact on the equivalence ratio of spray tail mixture. On the one hand, the gas-jet makes the spray tail move downstream, and the higher the velocity of gas-jet, the farther the spray tail leaves from the nozzle exit and the larger the zone affected by the gas-jet. For instance, when it is at t = 0.48 ms, the spray tail is around 4 mm away from the nozzle exit and the zone from the nozzle exit to 7.5 mm downstream is affected by the gas-jet with velocity of 100 m/s, while the spray tail is near 6 mm away from the nozzle exit and the zone from the nozzle exit to 10 mm downstream is affected in the case of gas-jet with velocity of 200 m/s. Under the action of the gas-jet, the stagnant fuel near the nozzle exit due to the influence of entrainment wave can flow downstream quickly, so the formation of fuel-lean zone near the nozzle exit is reduced and the following combustion environment is improved, decreasing the generation of soot and HC. Nevertheless, because of the rapid momentum exchange with low velocity mixture in spray tail, the gas-jet is not strong enough to affect the downstream zone any more. Thus, the ramping-down injection with gas-jet keeps the same equivalence ratio in the spray downstream as the ramping-down injection without gas-jet. In addition, owing to the limited effect range of gas-jet, there is almost no difference in spray tip Figure 13. Spray tip penetrations (a) and entrainment rates (b) of the ramping-down injection with gas-jet and the ramping-down injection without gas-jet (black solid lines: ramping-down injection without gas-jet; red chain lines: ramping-down injection with gas-jet of 100 m/s; blue dotted lines: ramping-down injection with gas-jet of 200 m/s). Figure 14 shows the variations of equivalence ratio of the ramping-down injection without gas-jet and the ramping-down injection with gas-jet. It is obvious that the gas-jet after EOI has a great impact on the equivalence ratio of spray tail mixture. On the one hand, the gas-jet makes the spray tail move downstream, and the higher the velocity of gas-jet, the farther the spray tail leaves from the nozzle exit and the larger the zone affected by the gas-jet. For instance, when it is at t = 0.48 ms, the spray tail is around 4 mm away from the nozzle exit and the zone from the nozzle exit to 7.5 mm downstream is affected by the gas-jet with velocity of 100 m/s, while the spray tail is near 6 mm away from the nozzle exit and the zone from the nozzle exit to 10 mm downstream is affected in the case of gas-jet with velocity of 200 m/s. Under the action of the gas-jet, the stagnant fuel near the nozzle exit due to the influence of entrainment wave can flow downstream quickly, so the formation of fuel-lean zone near the nozzle exit is reduced and the following combustion environment is improved, decreasing the generation of soot and HC. Nevertheless, because of the rapid momentum exchange with low velocity mixture in spray tail, the gas-jet is not strong enough to affect the downstream zone any more. Thus, the ramping-down injection with gas-jet keeps the same equivalence ratio in the spray downstream as the ramping-down injection without gas-jet. In addition, owing to the limited effect range of gas-jet, there is almost no difference in spray tip penetration and overall entrainment rates between the ramping-down injection without gas-jet and the ramping-down injection with gas-jet as Figure 13 shows. penetration and overall entrainment rates between the ramping-down injection without gas-jet and the ramping-down injection with gas-jet as Figure 13 shows. On the other hand, since the fuel in the spray tail flows downstream quickly due to the gas-jet, so more fuel accumulates in the new spray tail, which makes the equivalence ratio peak increase slightly. But the velocity of gas-jet has little impact on the equivalence ratio peak. Therefore, compared with the ramping-down injection without gas-jet, the one with gas-jet has more potential to achieve complete combustion and reduce the emission of soot and HC by organizing the ideal fuel-air distribution.

Conclusions
To find the optimized injection strategy for lower pollutant emissions, the mixing process of diesel spray under varying injection rate shapes has been studied in this paper by using a modified one-dimensional discrete spray model. The calculations were performed with five injection rate shapes: triangle, ramping-up, ramping-down, rectangle and trapezoid. Then the spray penetrations, entrainment rates and equivalence ratio distributed along the spray axial distance were analyzed. And this paper also discussed the ameliorative mixing characteristic of the multiple-injection strategy that was combined by the ramping-down and ramping-up injection rates, with the impact of the gas-jet after EOI on mixture distribution and pollutant emission reduction under the ramping-down injection rate shape investigated too. The conclusions are summarized as follows： (1) The higher injection rate leads to the larger spray tip penetration during the initial injection and the higher maximum injection rate leads to the larger spray tip penetration finally. (2) The ramping-up injection rate needs 1.5 times the injection duration for the entrainment wave to arrive at the spray tip. For the other four injection rates, the sprays can be treated as steady-like state, needing twice the injection duration from EOI to the time the entrainment wave reaches the spray tip. The concept using UEOI to represent the average entrainment wave and spray penetration speed is easy and effective to estimate the time when the entrainment wave arrives at the spray tip. (3) The ramping-down injection rate enhances the spray penetration and mixing rate during the initial period, and obtains an almost even distribution of the equivalence ratio along the spray axis after EOI. The ramping-up injection rate leads to the leaner mixture after EOI and even much leaner near the nozzle exit, which are useful to control the combustion process with longer ignition delay and decrease the soot emissions. On the other hand, since the fuel in the spray tail flows downstream quickly due to the gas-jet, so more fuel accumulates in the new spray tail, which makes the equivalence ratio peak increase slightly. But the velocity of gas-jet has little impact on the equivalence ratio peak. Therefore, compared with the ramping-down injection without gas-jet, the one with gas-jet has more potential to achieve complete combustion and reduce the emission of soot and HC by organizing the ideal fuel-air distribution.

Conclusions
To find the optimized injection strategy for lower pollutant emissions, the mixing process of diesel spray under varying injection rate shapes has been studied in this paper by using a modified one-dimensional discrete spray model. The calculations were performed with five injection rate shapes: triangle, ramping-up, ramping-down, rectangle and trapezoid. Then the spray penetrations, entrainment rates and equivalence ratio distributed along the spray axial distance were analyzed. And this paper also discussed the ameliorative mixing characteristic of the multiple-injection strategy that was combined by the ramping-down and ramping-up injection rates, with the impact of the gas-jet after EOI on mixture distribution and pollutant emission reduction under the ramping-down injection rate shape investigated too. The conclusions are summarized as follows: (1) The higher injection rate leads to the larger spray tip penetration during the initial injection and the higher maximum injection rate leads to the larger spray tip penetration finally. (2) The ramping-up injection rate needs 1.5 times the injection duration for the entrainment wave to arrive at the spray tip. For the other four injection rates, the sprays can be treated as steady-like state, needing twice the injection duration from EOI to the time the entrainment wave reaches the spray tip. The concept using U EOI to represent the average entrainment wave and spray penetration speed is easy and effective to estimate the time when the entrainment wave arrives at the spray tip.
(3) The ramping-down injection rate enhances the spray penetration and mixing rate during the initial period, and obtains an almost even distribution of the equivalence ratio along the spray axis after EOI. The ramping-up injection rate leads to the leaner mixture after EOI and even much leaner near the nozzle exit, which are useful to control the combustion process with longer ignition delay and decrease the soot emissions. (4) The spray characteristic summaries about the penetration and entrainment rate as well as equivalence ratio can provide a certain guidance and reference for practical application. (5) The multiple-injection consisting of the ramping-down and ramping-up injection rates can obtain much higher entrainment rates than the ramping-down single injection, and the equivalence ratio of the former is much lower near the nozzle exit than that of the latter. (6) The gas-jet after EOI has no obvious effect on the spray tip penetration and overall entrainment rate because the momentum of gas-jet is limited and the affected zone is small. However, the gas-jet after EOI has a great impact on the equivalence ratio of spray tail mixture. The spray tail will move downstream quickly under the action of the gas-jet. The higher the velocity of gas-jet, the farther the spray tail leaves from the nozzle exit and the larger the zone affected by the gas-jet. With the gas-jet after EOI, the equivalence ratio in spray tail will increase slightly due to the more fuel accumulated.