Experimental Research on the Macroscopic and Microscopic Spray Characteristics of Diesel-PODE 3-4 Blends

: Polyoxymethylene dimethyl ether (PODE) is a low-viscosity oxygenated fuel that can improve the volatility of blended fuels. In this work, the macroscopic and microscopic spray characteristics of diesel-PODE 3-4 under different ambient temperatures and injection pressures (IP) are studied. The studied blends consisted of pure diesel (P0), two diesel blend fuels of 20% (P20) and 50% (P50) by volume fraction of PODE 3-4 . The Mie scattering and Schlieren imaging techniques are used in the experiment. The results show that with the increase in IP, the vapor phase penetration distance and the average cone angle of the three fuels increased, and the Sauter mean diameter (SMD) of the three fuels decreased. When the ambient temperature increased, the vapor phase projection area and the average vapor phase cone angle of P20 and P50 increased, and the SMD decreased, but the vapor phase projection area of pure diesel did not change signiﬁcantly. The results indicate that the blended fuel with PODE 3-4 has better spray characteristics than P0 at low temperature, and the SMD hierarchy between the three fuels is P0 > P20 > P50. Through the visualization experiment, it is helpful to further understand the evaporation characteristics of different fuel properties and develop appropriate alternative diesel fuel.


Introduction
Diesel engines are widely used in ships, generators and heavy trucks due to their high power output and high thermal efficiency [1][2][3]. However, due to the related environmental problems, the shortage of oil resources and the requirements of national laws and regulations in recent years, the internal combustion engines, which depend on traditional oil as the power source, have been severely challenged, and the quest to find a novel fuel mixed with diesel fuel is seen as a solution to these problems [4,5]. At present, the most common method is to add biomass oxygenated fuel to diesel to form a suitable blend.
Biodiesel, alcohols and ethers are the most suitable additives for compression ignition internal combustion engines [6,7]. Many countries are committed to developing biofuels. In Europe, bioeconomy has been well developed especially in Ireland, Denmark, Portugal and Austria [8]. D'Adamo et al. [9] researched a circular economy model and found that applying biomethane in the transport system of Rome leads to a reduction of emissions. Compared with diesel, biodiesel is renewable, non-toxic and has higher hexadecane value. Biodiesel contains no aromatic hydrocarbons, which is capable for reduction of unburned hydrocarbons (HC), particulate matter (PM) and carbon monoxide (CO) emissions. Thus, Colombia's government has implemented laws about promoting biodiesel production industry [10]. Hassaan et al. [11] also suggested that the Egyptian government should pay more attention to the construction of biogas and biomethane production plants in the future. PODE n /natural gas blend resulted in fewer hydrocarbon, CO and soot emissions compared to the diesel/natural gas blend, and significantly improved the thermal efficiency.
These works have mainly focused on the in-cylinder combustion process of the engine. In fact, the fuel injection, atomization and mixing with air are critical to the entire combustion process, which is of great significance to energy saving and emissions reduction. Li et al. [29] researched the macroscopic and microscopic spray characteristics of diesel/PODE blended fuels and found that, with the increase in the proportion of PODE in diesel-PODE blend, the spray penetration distance decreased, though the average spray cone angle increased. In addition, both the characteristic diameters and the Sauter mean diameter decreased with the increase in the proportion of PODE in the diesel/PODE blend, though the range of relative size of droplets showed little change. However, Li et al. did not study the evaporation characteristics of PODE, whereas the test temperature was only the room temperature. Currently, the techniques of Mie scattering and diffuse back-illumination (DBI) images are widely used methods for measuring liquid spray characteristics, while schlieren image technology is suitable for studying spray evaporation characteristics [30][31][32]. Huang et al. [33] researched the spray, evaporation and combustion characteristics of ethanol/diesel blends under low temperature combustion (LTC) using DBI and schlieren image techniques. The results showed that the evaporation characteristics of the fuel blend increased the spray spreading angle and projected area, though they had little influence on the spray penetration distance, whereas the combustion phenomenon spread from the periphery behind the spray tip to the forward of the spray and then to the nozzle. Ma et al. [34] investigated the evaporation and spray characteristics of n-pentanol/diesel blends using DBI and schlieren image techniques, and found that pure diesel has a longer spray penetration distance and a smaller spray cone angle than n-pentanol in the absence of evaporation. When the ambient temperature exceeds 800 K under evaporation conditions, the spray penetration distance of pure diesel is reduced, although the spray penetration distance of n-pentanol increases and this trend is particularly remarkable for higher proportions of n-pentanol in diesel. Payri et al. [35] used Mie scattering and schlieren images techniques to study the effects of cylindrical and conical nozzles on the liquid and vapor phase spray characteristics of n-heptane, n-hexane and diesel alternative fuels (composed of n-tetradecane, n-decane and methylnaphthalene), and found that, for the same fuel under the same working conditions, the cylindrical nozzle has a smaller vapor phase penetration distance and liquid phase length than the conical nozzle. However, the spray spreading angle shows the opposite trend. Subsequently, Payri et al. [32] compared the experimental and calculated values of the vapor and liquid phase spray penetration distances of pure diesel oil at different IPs, ambient density and ambient temperature, and found that the experimental and calculated values can be in good agreement. Therefore, it is concluded that the IP affects the vapor phase length and the ambient temperature affects the liquid phase length.
However, the research on the spray (liquid penetration distance, liquid phase cone angle and liquid phase projection area) and evaporation characteristics (vapor phase cone angle, vapor phase penetration distance and vapor phase projection area) of a diesel/PODE n blend (chain length of PODE n is: n = 3-4) is still lacking, especially regarding the droplet size distribution characteristics of blended fuel. Many studies [20,36] have shown that the mixing of diesel and PODE 3-4 is the most suitable for application in diesel engines. Therefore, in this work, the effect of different ambient temperatures and IPs on the macroscopic and microscopic spray characteristics of diesel/PODE 3-4 blend is analyzed. From the results, the evaporation characteristics of different fuel properties and appropriate alternative diesel fuels can be further understood.  Figure 1 shows the schematic of the experimental setup, which mainly consists of a constant volume combustion bomb (CVCB), fuel supply system and image acquisition system. The specific parameters of CVCB are described in detail in a previous work [37]. In order to ensure the optical path and camera shooting, three quartz windows with the diameter of 110 mm were located on the side of the CVCB. A single hole electromagnetic valve injector was used for the experiments, and the fuel supply system was procured from Bosch's third-generation high-pressure common rail test rig. Each set of experiments was repeated three times and the average of the three test data was taken. The uncertainties of the apparatus are shown in Table 1.
R PEER REVIEW 4 of 23 Figure 1 shows the schematic of the experimental setup, which mainly consists of a constant volume combustion bomb (CVCB), fuel supply system and image acquisition system. The specific parameters of CVCB are described in detail in a previous work [37]. In order to ensure the optical path and camera shooting, three quartz windows with the diameter of 110 mm were located on the side of the CVCB. A single hole electromagnetic valve injector was used for the experiments, and the fuel supply system was procured from Bosch's third-generation high-pressure common rail test rig. Each set of experiments was repeated three times and the average of the three test data was taken. The uncertainties of the apparatus are shown in Table 1.  The images were taken using a FASTCAM-SA7 high-speed camera, which was manufactured by PhotronCorp, Tokyo, Japan. During the experiments, the external computer issued an injection command, and then the electronic control unit (ECU) in a high-pressure common rail system drove the injector and high-speed camera to work synchronously according to the detected fuel injection signal. The schlieren device adopted a Zshaped arrangement, which is mainly composed of a light source slit system and a knifeedge camera system. Figure 2 shows the schematic of the Mie scattering device, which is mainly used to  The images were taken using a FASTCAM-SA7 high-speed camera, which was manufactured by PhotronCorp, Tokyo, Japan. During the experiments, the external computer issued an injection command, and then the electronic control unit (ECU) in a high-pressure common rail system drove the injector and high-speed camera to work synchronously according to the detected fuel injection signal. The schlieren device adopted a Z-shaped arrangement, which is mainly composed of a light source slit system and a knife-edge camera system. Figure 2 shows the schematic of the Mie scattering device, which is mainly used to measure the liquid spray characteristics. During the experiments, two tungsten halogen lamps were placed in each of the two quartz windows. During operation, the tungsten halogen lamp emitted a light to illuminate the spray. Then, the scattered light, reflected by the spray, was received by the high-speed camera. Finally, the liquid phase spray boundary image was displayed on the computer screen. A micro-nano particle size analyzer (Winner 318A) was used to measure the microscopic characteristics of the spray and the micro experimental test structure diagram as shown in Figure 3a. The working principle of Winner 318A is shown in Figure 3b, which obtains particles' size by using a laser beam to test the intensity of scattering spectrum of particles.

Experimental Setup
FOR PEER REVIEW 5 of 23 by the spray, was received by the high-speed camera. Finally, the liquid phase spray boundary image was displayed on the computer screen. A micro-nano particle size analyzer (Winner 318A) was used to measure the microscopic characteristics of the spray and the micro experimental test structure diagram as shown in Figure 3a. The working principle of Winner 318A is shown in Figure 3b, which obtains particles' size by using a laser beam to test the intensity of scattering spectrum of particles.

Image Processing
The Mie scattering method can capture the liquid phase spray boundary image, while the schlieren method can capture the spray boundary image of the vapor phase. Both the methods use MATLAB's self-programming process for image processing. The steps of processing the image are similar for both methods. Firstly, the background image and the spray image are subtracted by using the MATLAB program to achieve the removal of background. Then, the edge pixels of the spray are detected using two thresholds of the Sobel and the cany operator, and the boundary curve of the spray is determined to obtain the liquid phase spray characteristic parameters. The Sobel operator and cany operator are used to edge the spray image detection, obtaining the coordinate values for the edge point of the image. Finally, the spray characteristic parameters of the liquid or vapor phase are obtained. Figure 4 shows the processed images from Mie scattering and schlieren

Image Processing
The Mie scattering method can capture the liquid phase spray boundary image, while the schlieren method can capture the spray boundary image of the vapor phase. Both the methods use MATLAB's self-programming process for image processing. The steps of processing the image are similar for both methods. Firstly, the background image and the spray image are subtracted by using the MATLAB program to achieve the removal of background. Then, the edge pixels of the spray are detected using two thresholds of the Sobel and the cany operator, and the boundary curve of the spray is determined to obtain the liquid phase spray characteristic parameters. The Sobel operator and cany operator are used to edge the spray image detection, obtaining the coordinate values for the edge point  Figure 4 shows the processed images from Mie scattering and schlieren methods.

Experimental Procedure
The experimental conditions are presented in Table 2. The fuel sample is based on No. 0 diesel fuel sold in China's petroleum market. The blended PODEn is obtained from Qingdao Tong Chuan Petrochemical Engineering Company, Qingdao, China, and is mainly composed of PODE2, PODE3 and PODE4 with the mass fractions of 2.6%, 88.9% and 8.5%, respectively. Since the main components are PODE3 and PODE4, this article uses PODE3-4 to represent PODEn. There were three different volume ratios of diesel/PODE3-4 blends used in the experiments, 0%, 20% and 50% with respect to PODE3-4, and were abbreviated as P0, P20 and P50, respectively. The relevant physical and chemical properties of the blends are presented in Tables 3 and 4.

Experimental Procedure
The experimental conditions are presented in Table 2. The fuel sample is based on No. 0 diesel fuel sold in China's petroleum market. The blended PODE n is obtained from Qingdao Tong Chuan Petrochemical Engineering Company, Qingdao, China, and is mainly composed of PODE 2 , PODE 3 and PODE 4 with the mass fractions of 2.6%, 88.9% and 8.5%, respectively. Since the main components are PODE 3 and PODE 4 , this article uses PODE [3][4] to represent PODE n . There were three different volume ratios of diesel/PODE 3-4 blends used in the experiments, 0%, 20% and 50% with respect to PODE [3][4] , and were abbreviated as P0, P20 and P50, respectively. The relevant physical and chemical properties of the blends are presented in Tables 3 and 4.

Definition of Spray Parameters
The evaporation characteristics of the spray influence the droplet size, whereas the size of the vapor penetration distance determines the rate of fuel impingement on the combustion chamber. The definition of the projected area is the calculated area inside the spray boundary. The projected area of vapor and liquid phase reflects the range of fuel spray diffusion, which is the combined effect of cone angle and penetration distance, while its size can reflect the quality of mixing of fuel and the surrounding environment gas [43]. The Mie scattering method is used to obtain the liquid phase spray boundary image, while the schlieren method is used to obtain the spray boundary image of the vapor phase. Therefore, it is necessary to define and distinguish the spray image parameters obtained by the two shooting methods. Figure 5 illustrates the way the spray characteristic parameters are defined in this work. The subscripts L and V represent the liquid phase and vapor phase spray characteristic parameters, respectively. The penetration distance for each phase is defined as the axial distance, which the fuel can reach farthest from the nozzle to the spray front, and is denoted by S V and S L for vapor and liquid phase, respectively. The spray cone angle for each phase is defined based on the method of Naber et al. [44], and is defined as the angle between the line from fuel injection from the nozzle to the half of the spray penetration distance and the tangent line along the spray contour. In Figure 5, the spray cone angle is represented by θ V and θ L for vapor and liquid phase, respectively. The sum of the values of the pixel area included in the entire spray image is the spray projection area, whereas the vapor and liquid phase projected areas are represented by A V and A L , respectively.   Figure 6 shows the evolution of vapor and liquid spray patterns for different fuel blends over time at different ambient temperatures (573 K, 623 K and 673 K) and the IP of 160 MPa. It can be known from Figure 6a,c,e that the vapor phase spray for each fuel blend at different ambient temperatures gradually increases with time. Meanwhile, the spray gradually changes from "dark black" (liquid phase) to "transparent color" (vapor phase) with the increase in temperature. This trend is particularly pronounced as the proportion of PODE3-4 in the blend increases. It can be seen from Figure 6b,d,f that the liquid phase spray of the three fuels gradually increases with time at 573 K. However, when the ambient temperature is 623 K, and both the P0 and P20 are in the middle and late spray (after 0.7 ms), the front end of the liquid phase spray gradually becomes blurred, producing a "mist", which is composed of a lot of small droplets. On the other hand, the sample P50 shows a rapid reduction of spray at the same temperature (623 K). As the ambient temperature increases to 673 K, the mist area of the three blends decreases. This is because the rate of evaporation of droplets becomes faster due to an increase in the ambient temperature, which results in rapid evaporation of small droplets formed at the edge of the liquid core and the front end of spray during the spraying process.

Effect of Ambient Temperature on Fuel Spray Characteristics
(a) Vapor phase spray image for P0.
(b) Liquid phase oil beam image for P0.  Figure 6 shows the evolution of vapor and liquid spray patterns for different fuel blends over time at different ambient temperatures (573 K, 623 K and 673 K) and the IP of 160 MPa. It can be known from Figure 6a,c,e that the vapor phase spray for each fuel blend at different ambient temperatures gradually increases with time. Meanwhile, the spray gradually changes from "dark black" (liquid phase) to "transparent color" (vapor phase) with the increase in temperature. This trend is particularly pronounced as the proportion of PODE 3-4 in the blend increases. It can be seen from Figure 6b,d,f that the liquid phase spray of the three fuels gradually increases with time at 573 K. However, when the ambient temperature is 623 K, and both the P0 and P20 are in the middle and late spray (after 0.7 ms), the front end of the liquid phase spray gradually becomes blurred, producing a "mist", which is composed of a lot of small droplets. On the other hand, the sample P50 shows a rapid reduction of spray at the same temperature (623 K). As the ambient temperature increases to 673 K, the mist area of the three blends decreases. This is because the rate of evaporation of droplets becomes faster due to an increase in the ambient temperature, which results in rapid evaporation of small droplets formed at the edge of the liquid core and the front end of spray during the spraying process.   Figure 6 shows the evolution of vapor and liquid spray patterns for different fuel blends over time at different ambient temperatures (573 K, 623 K and 673 K) and the IP of 160 MPa. It can be known from Figure 6a,c,e that the vapor phase spray for each fuel blend at different ambient temperatures gradually increases with time. Meanwhile, the spray gradually changes from "dark black" (liquid phase) to "transparent color" (vapor phase) with the increase in temperature. This trend is particularly pronounced as the proportion of PODE3-4 in the blend increases. It can be seen from Figure 6b,d,f that the liquid phase spray of the three fuels gradually increases with time at 573 K. However, when the ambient temperature is 623 K, and both the P0 and P20 are in the middle and late spray (after 0.7 ms), the front end of the liquid phase spray gradually becomes blurred, producing a "mist", which is composed of a lot of small droplets. On the other hand, the sample P50 shows a rapid reduction of spray at the same temperature (623 K). As the ambient temperature increases to 673 K, the mist area of the three blends decreases. This is because the rate of evaporation of droplets becomes faster due to an increase in the ambient temperature, which results in rapid evaporation of small droplets formed at the edge of the liquid core and the front end of spray during the spraying process.    Figure 7 shows the vapor and liquid phase penetration distances for the three blends at different ambient temperatures. It can be clearly seen from the picture that the vapor phase penetration distances of P0, P20 and P50 increase with time. However, when the ambient temperature is 673 K, the vapor phase penetration distance of the three fuels is significantly lower than those for 573 K and 623 K. The reason is that, at lower temperatures, the spray mainly develops in liquid form, while at a high temperature, the fuel evaporates in a large amount, causing the spray mainly to develop in vapor form. Furthermore, the vapor spray develops at a lower rate than that in liquid phase. Therefore, when the ambient temperature is 673 K, the vapor phase penetration distances of the three blends are smaller than the vapor phase penetration distances at low temperatures. For the liquid phase penetration distance, the fuel did not undergo significant evaporation at the temperatures of 573 K and 623 K. When the temperature was increased to 673 K, the (e)Vapor phase spray image for P50.

Effect of Ambient Temperature on Fuel Spray Characteristics
(f) Liquid phase spray image for P50.  Figure 7 shows the vapor and liquid phase penetration distances for the three blends at different ambient temperatures. It can be clearly seen from the picture that the vapor phase penetration distances of P0, P20 and P50 increase with time. However, when the ambient temperature is 673 K, the vapor phase penetration distance of the three fuels is significantly lower than those for 573 K and 623 K. The reason is that, at lower temperatures, the spray mainly develops in liquid form, while at a high temperature, the fuel evaporates in a large amount, causing the spray mainly to develop in vapor form. Furthermore, the vapor spray develops at a lower rate than that in liquid phase. Therefore, when the ambient temperature is 673 K, the vapor phase penetration distances of the three blends are smaller than the vapor phase penetration distances at low temperatures. For the liquid phase penetration distance, the fuel did not undergo significant evaporation at the temperatures of 573 K and 623 K.
When the temperature was increased to 673 K, the liquid phase penetration distance for P0 was substantially similar to the vapor phase penetration distance during the initial stage of the spray. However, the passage of time gradually separated both the values. The liquid penetration distance for each of P20 and P50 is expressed as the situation, in which the initial value of the spray quickly reaches a certain value and keeps fluctuating around this value. Therefore, it is indicated that diesel will slowly evaporate after the background temperature reaches a certain value until a certain period of time. However, when the diesel blend contains a large proportion of PODE 3-4 , the liquid phase and vapor phase spray are separated at the initial moment of the spray. The liquid phase penetration distance reaches the maximum value and becomes stable. Meanwhile, the amount of evaporation and the amount of fuel injection achieve certain equilibrium. In addition, compared with the ambient temperature of 573 K and 623 K, when the ambient temperature is 673 K, the vapor and liquid phase penetration distances of the three fuels are all lower. This is consistent with the results reported by Gimeno et al. [30], who stated that "the ambient temperature rises, and the fuel vapor and liquid phase penetration distances decrease". liquid phase penetration distance for P0 was substantially similar to the vapor phase penetration distance during the initial stage of the spray. However, the passage of time gradually separated both the values. The liquid penetration distance for each of P20 and P50 is expressed as the situation, in which the initial value of the spray quickly reaches a certain value and keeps fluctuating around this value. Therefore, it is indicated that diesel will slowly evaporate after the background temperature reaches a certain value until a certain period of time. However, when the diesel blend contains a large proportion of PODE3-4, the liquid phase and vapor phase spray are separated at the initial moment of the spray. The liquid phase penetration distance reaches the maximum value and becomes stable. Meanwhile, the amount of evaporation and the amount of fuel injection achieve certain equilibrium. In addition, compared with the ambient temperature of 573 K and 623 K, when the ambient temperature is 673 K, the vapor and liquid phase penetration distances of the three fuels are all lower. This is consistent with the results reported by Gimeno et al. [30], who stated that "the ambient temperature rises, and the fuel vapor and liquid phase penetration distances decrease".  Figure 8 shows the vapor and liquid cone angles of the three fuels at different ambient temperatures. It can be known that the vapor phase cone angles of P0, P20 and P50 increase with the increase in ambient temperature. However, the change in vapor phase cone angle does not agree with the increase in vapor phase conge angle of P0 before 0.6 ms, whereas the temperature increases gradually after 0.6 ms. The reason is that the viscosity of the diesel oil is high, which initially causes the droplets to stick together after the fuel  Figure 8 shows the vapor and liquid cone angles of the three fuels at different ambient temperatures. It can be known that the vapor phase cone angles of P0, P20 and P50 increase with the increase in ambient temperature. However, the change in vapor phase cone angle does not agree with the increase in vapor phase conge angle of P0 before 0.6 ms, whereas the temperature increases gradually after 0.6 ms. The reason is that the viscosity of the diesel oil is high, which initially causes the droplets to stick together after the fuel is sprayed from the nozzle. Due to this, it is not easy for the droplets to diffuse in the radial direction. As time and temperature increase, the spray sharply evaporates, and the situation improves. In addition, when the temperature is 673 K, compared with the pure diesel, the average vapor cone angles of P20 and P50 increased by 1.9 • and 3.3 • , respectively, and the average liquid cone angle reduced by 1 • and 2.2 • , respectively, indicating that the addition of PODE [3][4] to diesel in a certain proportion can improve the evaporation and diffusion of the spray. It can also be seen that the liquid phase cone angles of P0, P20 and P50 decrease with the ambient temperature increase. The reason is that the ambient temperature increases, causing the heat exchange between the liquid phase fuel in the spray front and the edge region and the high-temperature ambient gas to increase. Therefore, the rate of evaporation of liquid fuel greatly increases, resulting in the decrease in liquid phase spray cone angle. is sprayed from the nozzle. Due to this, it is not easy for the droplets to diffuse in the radial direction. As time and temperature increase, the spray sharply evaporates, and the situation improves. In addition, when the temperature is 673 K, compared with the pure diesel, the average vapor cone angles of P20 and P50 increased by 1.9°and 3.3°, respectively, and the average liquid cone angle reduced by 1° and 2.2°, respectively, indicating that the addition of PODE3-4 to diesel in a certain proportion can improve the evaporation and diffusion of the spray. It can also be seen that the liquid phase cone angles of P0, P20 and P50 decrease with the ambient temperature increase. The reason is that the ambient temperature increases, causing the heat exchange between the liquid phase fuel in the spray front and the edge region and the high-temperature ambient gas to increase. Therefore, the rate of evaporation of liquid fuel greatly increases, resulting in the decrease in liquid phase spray cone angle.  Figure 9 shows the vapor and liquid phase projected areas of the three tested fuels at different ambient temperatures. It can be seen that the vapor and liquid phase projection areas of P0 are not much different for the ambient temperatures of 573 K and 623 K and before 0.9 ms. However, the vapor phase projection area of P0 after 0.5 ms at ambient temperature of 673 K is larger than the other two temperatures. The results for P20 and P50 show that the higher the temperature, the larger the vapor phase projection area and smaller the liquid phase projection area. It shows that the increase in pure diesel in the low temperature range (within 623 K) has little influence on the evaporation characteris-  Figure 9 shows the vapor and liquid phase projected areas of the three tested fuels at different ambient temperatures. It can be seen that the vapor and liquid phase projection areas of P0 are not much different for the ambient temperatures of 573 K and 623 K and before 0.9 ms. However, the vapor phase projection area of P0 after 0.5 ms at ambient temperature of 673 K is larger than the other two temperatures. The results for P20 and P50 show that the higher the temperature, the larger the vapor phase projection area and smaller the liquid phase projection area. It shows that the increase in pure diesel in the low temperature range (within 623 K) has little influence on the evaporation characteristics and spray diffusion, while the fuel blend with PODE 3-4 still has good diffusion capability and evaporability at low temperatures. For the liquid fuel projected area of the mixed fuel, with the ambient temperature, the projected area decreases. When the temperature is 673 K, the liquid projection area of the three blends tends to flatten with time. The reason is that the higher the ambient temperature, the more evaporation on both sides of the liquid fuel spray and the faster the evaporation rate, which results in a smaller projected area. When the evaporation speed of the droplets is equal to the diffusion speed, the liquid projected area reaches a stable state. tics and spray diffusion, while the fuel blend with PODE3-4 still has good diffusion capability and evaporability at low temperatures. For the liquid fuel projected area of the mixed fuel, with the ambient temperature, the projected area decreases. When the temperature is 673 K, the liquid projection area of the three blends tends to flatten with time. The reason is that the higher the ambient temperature, the more evaporation on both sides of the liquid fuel spray and the faster the evaporation rate, which results in a smaller projected area. When the evaporation speed of the droplets is equal to the diffusion speed, the liquid projected area reaches a stable state.  Figure 10 shows the development of vapor and liquid spray patterns at different temperatures for different IPs (80, 120 and 160 MPa) at the ambient temperature of 673 K. It can be seen that the vapor phase and the liquid phase spray lengths for P0, P20 and P50 increase with the IP increase. For the vapor phase oil bundle, the vapor phase spray profile grows evenly larger under high IPs. Before 0.3 ms, the vapor phase fuel spray regions of the three fuels are dark black, whereas the vapor phase and liquid phase spray regions substantially coincide. However, the liquid phase of the fuel is relatively large. After 0.3 ms, the vapor phase oil jets of the three blends showed obvious vapor and liquid phase separation with time. This situation is particularly remarkable under the condition of high IP and high mixing ratio of PODE3-4. For liquid phase spray, the liquid phase spray length of pure diesel oil under the same IP is significantly larger than the liquid phase spray lengths of P20 and P50. In addition, at a higher IP, the front end of the liquid fuel spray of  Figure 10 shows the development of vapor and liquid spray patterns at different temperatures for different IPs (80, 120 and 160 MPa) at the ambient temperature of 673 K. It can be seen that the vapor phase and the liquid phase spray lengths for P0, P20 and P50 increase with the IP increase. For the vapor phase oil bundle, the vapor phase spray profile grows evenly larger under high IPs. Before 0.3 ms, the vapor phase fuel spray regions of the three fuels are dark black, whereas the vapor phase and liquid phase spray regions substantially coincide. However, the liquid phase of the fuel is relatively large. After 0.3 ms, the vapor phase oil jets of the three blends showed obvious vapor and liquid phase separation with time. This situation is particularly remarkable under the condition of high IP and high mixing ratio of PODE [3][4] . For liquid phase spray, the liquid phase spray length of pure diesel oil under the same IP is significantly larger than the liquid phase spray lengths of P20 and P50. In addition, at a higher IP, the front end of the liquid fuel spray of the three fuels produces a "mist", which is composed of broken atomized droplets, indicating that increasing the IP promotes the atomization of the fuel. the three fuels produces a "mist", which is composed of broken atomized droplets, indicating that increasing the IP promotes the atomization of the fuel.   Figure 11 shows the penetration distance for P0, P20 and P50 over time for different IPs. It can be seen that the vapor and liquid phase penetration distances of the three fuels increase with the increase in IP. This is because when the IP increases, the pressure difference between the inside and outside of the nozzle increases, and the initial kinetic energy of the oil droplets becomes higher [45]. This increased kinetic energy makes the spray penetrate to a longer distance. In addition, the liquid phase penetration distances of P0 and P20 change significantly with the increase in IP. However, the liquid phase penetration distance of P50 does not show much change with the increase in injection pressure. This is because the P50 blend has a lower viscosity relative to P0 and P20, and the lower viscosity will make the fuel break into small droplets more easily after being sprayed from the nozzle [46]. Additionally, at high temperatures, small droplets at the liquid spray front are rapidly evaporated, which cause the liquid penetration distance of P50 to change little with the increase in IP. This process is beneficial to reduce the proportion of liquid phase fuel in the flame region after ignition and improve the combustion efficiency in the cylinder.  Figure 11 shows the penetration distance for P0, P20 and P50 over time for different IPs. It can be seen that the vapor and liquid phase penetration distances of the three fuels increase with the increase in IP. This is because when the IP increases, the pressure difference between the inside and outside of the nozzle increases, and the initial kinetic energy of the oil droplets becomes higher [45]. This increased kinetic energy makes the spray penetrate to a longer distance. In addition, the liquid phase penetration distances of P0 and P20 change significantly with the increase in IP. However, the liquid phase penetration distance of P50 does not show much change with the increase in injection pressure. This is because the P50 blend has a lower viscosity relative to P0 and P20, and the lower viscosity will make the fuel break into small droplets more easily after being sprayed from the nozzle [46]. Additionally, at high temperatures, small droplets at the liquid spray front are rapidly evaporated, which cause the liquid penetration distance of P50 to change little with the increase in IP. This process is beneficial to reduce the proportion of liquid phase fuel in the flame region after ignition and improve the combustion efficiency in the cylinder.  Figure 12 shows the vapor and liquid cone angles for P0, P20 and P50 under different IPs. The results show that it is not intuitive to see the regularity of the vapor and liquid phase cone angles of the three fuels as a function of IP. Therefore, the average values of the vapor and liquid phase cone angles are shown in Figure 13. It can be seen that the average vapor phase cone angles of the three fuels increase with the increase in IP. When  Figure 12 shows the vapor and liquid cone angles for P0, P20 and P50 under different IPs. The results show that it is not intuitive to see the regularity of the vapor and liquid phase cone angles of the three fuels as a function of IP. Therefore, the average values of the vapor and liquid phase cone angles are shown in Figure 13. It can be seen that the average vapor phase cone angles of the three fuels increase with the increase in IP. When the IP is increased from 80 MPa to 160 MPa, the average vapor phase cone angles for P0, P20 and P50 increase by 1.7 • , 1.8 • and 2.5 • , respectively. This is because when the IP is increased, the kinetic energy of the fuel from the nozzle outlet increases, which in turn leads to enhanced mixing of the spray with the surrounding gas [47], thus resulting in an increase in the vapor cone angle. In addition, regardless of the IP, the average vapor phase cone angles for the three samples are found in the following ascending order: P0 < P20 < P50. This is because as the proportion of PODE [3][4] in the blend increases, the viscosity of the mixed fuel decreases compared to pure diesel. This leads to an increase in the resistance of spray to ambient gas [30], causing the average vapor cone angle to increase. These results are consistent with those reported by Valentino et al. [48], who concluded that "reducing the viscosity of the fuel will increase the spray cone angle". Interestingly, the average liquid phase cone angle decreases as the proportion of PODE 3-4 in the blend increases. This is because the viscosity of the mixed fuel is lowered, which causes the fuel to break more easily and hence, atomize after being ejected from the nozzle. Therefore, it becomes easier to evaporate under high temperature conditions, and results in a decrease in average liquid phase cone angle. It can also be seen from Figure 13b that the average liquid phase cone angle of the three fuels does not differ much under different IPs. This because although increasing the IP will cause the spray to spread to both sides, at the same time, the higher IP will cause the fuel to break more severely. Due to this reason, the spray evaporation speed will increase, and therefore, the change in average liquid cone angle will not be obvious.

Effect of Injection Pressure on Spray Characteristics
the IP is increased from 80 MPa to 160 MPa, the average vapor phase cone angles for P0, P20 and P50 increase by 1.7°, 1.8° and 2.5°, respectively. This is because when the IP is increased, the kinetic energy of the fuel from the nozzle outlet increases, which in turn leads to enhanced mixing of the spray with the surrounding gas [47], thus resulting in an increase in the vapor cone angle. In addition, regardless of the IP, the average vapor phase cone angles for the three samples are found in the following ascending order: P0 < P20 < P50. This is because as the proportion of PODE3-4 in the blend increases, the viscosity of the mixed fuel decreases compared to pure diesel. This leads to an increase in the resistance of spray to ambient gas [30], causing the average vapor cone angle to increase. These results are consistent with those reported by Valentino et al. [48], who concluded that "reducing the viscosity of the fuel will increase the spray cone angle". Interestingly, the average liquid phase cone angle decreases as the proportion of PODE3-4 in the blend increases. This is because the viscosity of the mixed fuel is lowered, which causes the fuel to break more easily and hence, atomize after being ejected from the nozzle. Therefore, it becomes easier to evaporate under high temperature conditions, and results in a decrease in average liquid phase cone angle. It can also be seen from Figure 13b that the average liquid phase cone angle of the three fuels does not differ much under different IPs. This because although increasing the IP will cause the spray to spread to both sides, at the same time, the higher IP will cause the fuel to break more severely. Due to this reason, the spray evaporation speed will increase, and therefore, the change in average liquid cone angle will not be obvious.   Figure 14 shows the vapor and liquid phase projection areas of P0, P20 and P50 under different IPs. It can be seen that the vapor phase projected areas of P0, P20 and P50 increase with the IP increase, indicating that increasing the IP will increase the propagation speed of the spray. This will also improve the spread of spray to surrounding [49], thereby improving the utilization of air in the cylinder. It can be seen from the liquid projection area that the liquid projection area increases slightly with the increase in injection pressure. This is because with higher injection pressure, the axial velocity and radial momentum of the fuel are greater after it is ejected from the nozzle. Due to this, the range of oil beam space diffusion becomes wider. However, higher IP will aggravate the degree of breakage of droplets, which leads to the speed of partial evaporation of oil beam outline, due to which the overall liquid projection area does not increase much. In addition, it can be seen that the liquid phase projected area curve of P0 is substantially coincidental with the vapor phase projected area before 0.5 ms and for the IPs of 80 MPa and 120 MPa. However, for P20 and P50, these timings are 0.3 ms and before 0.2 ms, respectively, and after these time intervals, curves for both the samples gradually start to separate from each other. It shows that, compared with the pure diesel, the mixed fuel can produce diffusion and evaporation earlier, which is very important for the fuel and gas mixing in the cylinder when the engine is actually working.  Figure 14 shows the vapor and liquid phase projection areas of P0, P20 and P50 under different IPs. It can be seen that the vapor phase projected areas of P0, P20 and P50 increase with the IP increase, indicating that increasing the IP will increase the propagation speed of the spray. This will also improve the spread of spray to surrounding [49], thereby improving the utilization of air in the cylinder. It can be seen from the liquid projection area that the liquid projection area increases slightly with the increase in injection pressure. This is because with higher injection pressure, the axial velocity and radial momentum of the fuel are greater after it is ejected from the nozzle. Due to this, the range of oil beam space diffusion becomes wider. However, higher IP will aggravate the degree of breakage of droplets, which leads to the speed of partial evaporation of oil beam outline, due to which the overall liquid projection area does not increase much. In addition, it can be seen that the liquid phase projected area curve of P0 is substantially coincidental with the vapor phase projected area before 0.5 ms and for the IPs of 80 MPa and 120 MPa. However, for P20 and P50, these timings are 0.3 ms and before 0.2 ms, respectively, and after these time intervals, curves for both the samples gradually start to separate from each other. It shows that, compared with the pure diesel, the mixed fuel can produce diffusion and evaporation earlier, which is very important for the fuel and gas mixing in the cylinder when the engine is actually working.   Figure 14 shows the vapor and liquid phase projection areas of P0, P20 and P50 under different IPs. It can be seen that the vapor phase projected areas of P0, P20 and P50 increase with the IP increase, indicating that increasing the IP will increase the propagation speed of the spray. This will also improve the spread of spray to surrounding [49], thereby improving the utilization of air in the cylinder. It can be seen from the liquid projection area that the liquid projection area increases slightly with the increase in injection pressure. This is because with higher injection pressure, the axial velocity and radial momentum of the fuel are greater after it is ejected from the nozzle. Due to this, the range of oil beam space diffusion becomes wider. However, higher IP will aggravate the degree of breakage of droplets, which leads to the speed of partial evaporation of oil beam outline, due to which the overall liquid projection area does not increase much. In addition, it can be seen that the liquid phase projected area curve of P0 is substantially coincidental with the vapor phase projected area before 0.5 ms and for the IPs of 80 MPa and 120 MPa. However, for P20 and P50, these timings are 0.3 ms and before 0.2 ms, respectively, and after these time intervals, curves for both the samples gradually start to separate from each other. It shows that, compared with the pure diesel, the mixed fuel can produce diffusion and evaporation earlier, which is very important for the fuel and gas mixing in the cylinder when the engine is actually working.

Microscopic Spray Characteristics
To further understand the atomization characteristics of blended fuel, in this section, the microscopic spray characteristics of blended fuels at different temperatures (303 K, 363 K) and injection pressures are studied.

Microscopic Spray Characteristics
To further understand the atomization characteristics of blended fuel, in this section, the microscopic spray characteristics of blended fuels at different temperatures (303 K, 363 K) and injection pressures are studied. Figure 15 shows the influence of different temperature and IPs on the SMD of test fuels' droplets. As can been seen from the figure, with the IP and temperature increases, the SMD of the test fuels' droplets are reduced. This is due to the increase in IP, which increases the degree of oil bundle breakage and tends to produce more small droplets; the higher of temperature, the more volatile the fuel is, and more tiny droplets can be added, so the SMD of the test fuel is reduced. It can also be seen from the figure that as the proportion of PODE 3-4 added increases, the fuel's SMD decreases. This is because of the low viscosity of PODE 3-4 (see Table 4), as the proportion of PODE 3-4 increases, the viscosity of blended fuel decreases, which is conducive to the improvement of volatility and the SMD reduced.

Microscopic Spray Characteristics
To further understand the atomization characteristics of blended fuel, in t the microscopic spray characteristics of blended fuels at different temperatures K) and injection pressures are studied. Figure 15 shows the influence of different temperature and IPs on the S fuels' droplets. As can been seen from the figure, with the IP and temperatur the SMD of the test fuels' droplets are reduced. This is due to the increase in increases the degree of oil bundle breakage and tends to produce more small d higher of temperature, the more volatile the fuel is, and more tiny droplets can so the SMD of the test fuel is reduced. It can also be seen from the figure that portion of PODE3-4 added increases, the fuel's SMD decreases. This is because viscosity of PODE3-4 (see Table 4), as the proportion of PODE3-4 increases, the blended fuel decreases, which is conducive to the improvement of volatility an reduced.  Figure 16 shows the effects of different temperature and IPs on the drop tribution. It can be seen from the figure, when the IP is 80 MPa, the P20 and size curves shift to the left with the temperature rises, while P0 shifts to the r because as the temperature rises, the volatility of the blended fuel increases, Figure 15. The influence of different temperature and IPs on the SMD of test fuels' droplets. Figure 16 shows the effects of different temperature and IPs on the droplet size distribution. It can be seen from the figure, when the IP is 80 MPa, the P20 and P50 particle size curves shift to the left with the temperature rises, while P0 shifts to the right. This is because as the temperature rises, the volatility of the blended fuel increases, and a large number of small diameter droplets are generated. As the viscosity of P0 is relatively high, the oil bundle is easy to be broken and form large diameter droplets, while the kinetic energy of the larger droplets is relatively small and easy to be absorbed and merged into larger droplets, so more large diameter droplets are generated. With the IP increases, the droplets' size distribution curve of the fuel shifts to the left at 303 K. This is due to the higher IP, which improves the atomization and generates more small droplets, especially P20 and P50, which are highly volatile, generating more droplets of small size than P0. When the IPs is 160 MPa, with the temperature rises and the P0 curve shifts to the left, while the P20 and P50 curves shift to the right and with a bimodal distribution. This is because under high temperature and high IP, P0 can have a higher degree of atomization and smaller droplet size, but for P20 and P50, under high IP, it is easier to generate more smaller droplets (as shown in Figure 16d). However, as the temperature rises, the movement speed of the droplets increases and the probability of mutual adsorption between droplets increases, so the probability of large-size droplets appearing increases.
(b) 363 K Figure 15. The influence of different temperature and IPs on the SMD of test fuels' droplets. Figure 16 shows the effects of different temperature and IPs on the droplet size distribution. It can be seen from the figure, when the IP is 80 MPa, the P20 and P50 particle size curves shift to the left with the temperature rises, while P0 shifts to the right. This is because as the temperature rises, the volatility of the blended fuel increases, and a large number of small diameter droplets are generated. As the viscosity of P0 is relatively high, the oil bundle is easy to be broken and form large diameter droplets, while the kinetic energy of the larger droplets is relatively small and easy to be absorbed and merged into larger droplets, so more large diameter droplets are generated. With the IP increases, the droplets' size distribution curve of the fuel shifts to the left at 303 K. This is due to the higher IP, which improves the atomization and generates more small droplets, especially P20 and P50, which are highly volatile, generating more droplets of small size than P0. When the IPs is 160 MPa, with the temperature rises and the P0 curve shifts to the left, while the P20 and P50 curves shift to the right and with a bimodal distribution. This is because under high temperature and high IP, P0 can have a higher degree of atomization and smaller droplet size, but for P20 and P50, under high IP, it is easier to generate more smaller droplets (as shown in Figure 16d). However, as the temperature rises, the movement speed of the droplets increases and the probability of mutual adsorption between droplets increases, so the probability of large-size droplets appearing increases.

Conclusions
In the CVCB premixed combustion device, the macroscopic and microscopic spray characteristics of the three different diesel and PODE blends (P0, P20 and P50) were studied under different temperatures and IPs. Based on the results, the following conclusions are drawn.
1. When the ambient temperature increases, the liquid cone angles decrease and the vapor cone angles increase for test fuels. When the ambient temperature is 673 K, P0 has the smallest average vapor phase cone angles, but the average liquid cone angle is largest. 2. For the ambient temperatures of 573 K and 623 K, the P0, P20 and P50 vapor and liquid phase penetrations do not change much. When the temperature increases to 673 K, the vapor and liquid phase penetration distances of the three fuels decrease,

Conclusions
In the CVCB premixed combustion device, the macroscopic and microscopic spray characteristics of the three different diesel and PODE blends (P0, P20 and P50) were studied under different temperatures and IPs. Based on the results, the following conclusions are drawn.

1.
When the ambient temperature increases, the liquid cone angles decrease and the vapor cone angles increase for test fuels. When the ambient temperature is 673 K, P0 has the smallest average vapor phase cone angles, but the average liquid cone angle is largest.

2.
For the ambient temperatures of 573 K and 623 K, the P0, P20 and P50 vapor and liquid phase penetrations do not change much. When the temperature increases to 673 K, the vapor and liquid phase penetration distances of the three fuels decrease, while the liquid phase penetration distance decreases the most. 3.
The vapor phase projection areas of P20 and P50 show that the higher the ambient temperature, the larger the projected area, whereas the liquid phase projected area shows the opposite trend. The vapor and liquid phase projections of P0 at ambient temperatures of 573 K and 623 K do not change much, and the situation improved when the temperature increased to 673 K.

4.
With the IP increase, the vapor phase penetration distance and the vapor phase cone angles of P0, P20 and P50 increase. Meanwhile, the average liquid phase cone angles of the three fuels will decrease. The liquid phase penetration distances increase in P0 and P20 with the IP increase, although the change in the liquid phase penetration distance of P50 is not obvious.

5.
As the IP increases, the vapor and liquid phase projection areas of the test fuels increase. Compared with the P0, the vapor and liquid phase projection area curves of the blends can be distinguished at an early stage, indicating that the blended fuel can undergo earlier diffusion and evaporation, which are critical to the mixing of oil and gas in the cylinder during actual engine operation. 6.
For the vapor phase, both the smallest cone angle and penetration distance is at 673 K and 80 MPa, while the largest projected area is at 673 K and 160 MPa of P50. For the liquid phase, it has the largest cone angle at 673 K and 160 MPa, and smallest penetration at 673 K and 80 MPa. At 573 K, 160 MPa, the liquid phase projected area of P0 is the largest of all tested points. 7.
When the temperature and IPs increase, the SMD of the three fuels decrease, and the SMD hierarchy between the three fuels is P0 > P20 > P50. As the IPs and temperature increase, the droplet size decreases, especially when the IP is 160 MPa and the temperature is 363 K, the droplet size distribution of P50 and P20 is bimodal, and the droplet size is smaller.
This article accomplished research on the characteristics of the three different diesel and PODE blends, but more situations should be considered in the next experiments. In the future, the microscopic spray characteristics in high ambient temperature or ambient pressure should become the key point of discussion.