Threedimensional (3D) computational code was implemented to solve conservation equations based on finite volume method as to simulate 1.8 L Ford diesel engine. Velocity and pressure of each computational cell is achieved by SIMPLE (semiimplicit method for pressurelinked equations) algorithm. For the exergetic aspect, the initial condition is set at 0.1 MPa and 300 K. The engine modeling is performed with 130 °, 140 °, and 150 ° with respect to xaxis under 1500 and 2500 rpm engine speeds. The results, however, indicate better air/fuel mixture (near stoichiometric equivalence ratio) for 130 ° of injection angle, albeit smaller spray droplets (lower sauter mean diameter) were introduced with 140 °. It is seen that higher soot and NOx mass fraction is attributed to 1500 rpm engine speed. The highest NOx and soot are exhausted at 130 ° and 150 ° of injection, respectively. Second law efficiency was calculated for different spray angle and engine speed schemes such that 36.62%, 30.2%, and 32.07% are associated with 130 °, 140 °, and 150 ° of injection angle under 1500 rpm, respectively. In terms of engine performance, that is, indicated mean effective pressure, indicated specific fuel consumption, and temperature, the best performance metrics are of 130 ° equal to 15.4 bar, 0.3856 kg/kWh, and 2074.97 K under 1500 rpm, respectively. Instant irreversibility rate is the highest amount with peak value of 17.48 J/deg for 130 deg1500 rpm, while 140 ° shows higher mean irreversibility rate over crank angle (CA) period.
Energy utilization is a vital issue; hence, presentation of a framework to serve the purpose of both frugality and high efficiency seems inevitable. Energy evaluation of internal combustion engines (ICEs) alone, fails to determine the key reasons governing the inefficiency and combustion deficiency since energy balance equation gives the amount of energy either as work delivery or heat generation/dissipation. The study in quality of energy (exergy) sheds light on the sources, which can assist to reduce the losses from a defined control volume. Nowadays special attention is paid for the estimation of operational parameters in energyrelated devices through numerical methods. Several studies were devoted to the casesensitive investigation of the engines and their performance with various operational conditions [15]. The effect of spray injection angle on combustion and exhaust emissions were conducted experimentally [6]. The tests were performed by 60 ° and 70 ° injection angle and compared with conventional angle of 156 °. The results, however, showed that the pressure for narrowangle injection (60 ° or 70 °) is increasing compared to wideangle strategy (156 °). Park and Reitz [7] applied 2sprayangle of grouphole nozzle aimed at squish zone (150 ° and 170 °) and bowl (80 ° spray angle). Their work revealed better air/fuel mixture, higher incylinder pressure, and reduced emissions with sprayguided injection compared to the baseline conventional case. The authors' previous work covers availability and irreversibility analysis with engine and combustion process based on different operative parameters such as various fuels, inlet charge temperature, and equivalence ratio [811].
Literature survey revealed the following results for the consideration of the first and second law in the ICEs; Ghazikhani et al. [12] performed an experimental investigation of exhaust gas recirculation (EGR) impact on irreversibility and brakespecific fuel consumption of diesel engine. They concluded that with EGR introduction, the total incylinder irreversibility increases due to the extension of flame region. The effect of nheptane and natural gas blends were surveyed in homogeneous charge compression ignition (HCCI) engines from the second thermodynamic law viewpoint [13]. A complete chemical kinetic mechanism was utilized in order to analyze the availability of system for varied EGR and natural gas percentage addition. They proved that the natural gas addition accounts for lowering exergy destruction accompanied by the second law efficiency rise. Ozkan et al. [14] discussed the effect of preinjection timing on diesel engines exergetic and energetic performance elaborately. Their results divulged that with adopting appropriate preinjection strategy, no perceptible influence could be noticed on thermal and exergetic efficiency; however, NOx emissions were reduced significantly by 7.4%. A new study has been dedicated to determine the optimal operating condition of HCCI engine that is based on crankangle resolved method [15]. It was observed that the equivalence ratio should be maintained at high values for majority of engine conditions while pressure was used to regulate engine load. Moreover, it was indicated that combustion timing has to be adjusted just before the sharp increase in unburned species losses. Jafarmadar and Zehni [16] carried out exergy analysis in combustion chamber of the IDI diesel engine by threedimentional modeling. They investigated dwell time duration effect on exergy terms, which resulted in work exergy and exergy efficiency decrease when dwell time increased from 5 ° to 30 °CA.
This work deals with different spray injection angle into the bowl section of CI diesel engine while 1500 and 2500 rpm engine speeds were employed. Spray structure, combustion, and emissions were inspected from the first law perspective. Then the second law analyses were conducted in various forms such as the indicated work, heat release, and chemical exergy as well as the irreversibility term. Finally, regression analysis was performed to develop a fourth order polynomial at six specified cases. Finally, among six cases with different speed and injection angle, 130 deg1500 rpm shows optimal exergetic efficiency, although, higher emission was issued.
The simulation is according to 1.8L Ford diesel engine. The numerical modeling in the simulation is detailed in ref. [17]. AVL FIRE code was employed for multidimensional simulation of the respective diesel engine. The swirl/tumble motion in combustion chamber of diesel engine has positive characteristics such as lowNOx emission, low noise, and high engine speed [1820]. The schematic of combustion chamber with guided spray injection is illustrated in Figure 1. Three sprays of 130 °, 140 °, and 150 ° with respect to xaxis are shown. As seen, 140 ° spray angle collides at central spot of the bowl, 130 ° at near pedestal section inducing the counter clockwise spray diffusion, and 150 ° at near squish zone with clockwise swirl motion of the spray. Boundary conditions, dead state condition, flow and emission simulation submodels are listed in Table 1. The engine specifications are summarized in Table 2 (for the sake of brevity further explanation on numerical implementation was omitted, more information are detailed in [8]). 3dimensional gridding of the chamber with sufficient 168,421 mostly unstructured cells was constructed (Fig. 2). The veracity of simulated results is validated according to pressure courses of 1500 and 2500 rpm engine speeds. Figure 3 shows excellent agreement between the modeling and experimental results. The obtained results by computational fluid dynamics (CFD) numerical method were compared with experimental results in terms of the incylinder pressure at 1500 and 2500 rpm engine speeds [21]. The highest discrepancies at peak values for 1500 and 2500 engine speeds are lower than 0.1% between experimental and modeling curves.
Head temperature  550.15 K 
Piston temperature  575.15 K 
Cylinder temperature  475.15 K 
Dead state pressure  1.01325 bar 
Dead state temperature  290.15 K 
Spray breakup  Modified KHRT 
Combustion model  Eddy BreakUp 
Turbulence model  kζf 
Evaporation  Dukowicz 
NO  Extended Zeldovich 
Soot  KennedyHiroyasuMagnussen 
Wall treatment  Hybridized wall 
Heat transfer wall mode  Standard wall function 
Residual gas ratio  0.5 
Fuel injection quantity (mg/cycle)  31.3 
Bore× stroke  82.5 × 82 mm 
Displacement  438 cm^{3}/cylinder 
Compression ratio  19.5:1 
Swirl ratio @ IVC  3 
Rail pressure  540–1255 bar (based on engine speed) 
Nozzle geometry  5 × 0.15 mm 
Bowl diameter  0.0483 m 
Bowl depth  0.0144 m 
Clearance  0.86 mm 
Number of nozzle holes  4 
Injection start timing  3 °CA BTDC 
Injection spray angle  130 °, 140 °, 150 ° 

Figure 1. Schematic illustration of guided spray scheme.


Figure 2. CFD meshing model of combustion chamber at TDC.


Figure 3. Pressure history for experimental and simulation for (A) 1500 rpm, (B) 2500 rpm [20].

Exergy is defined as a feature representing the capability of a system in transforming a given amount of energy at prescribed condition to useful work. In other words, exergy characterizes the maximum useful work attained from a system, which is plausible when a reversible process can be established between the initial and final point. The maximum extractable work depends on both ambient condition and the characteristics of the system. In general, total exergy is classified to thermo mechanical exergy and chemical exergy. The former governs the case when no chemical reaction tendency between the working fluid and the environment exists, whereas the latter is related to thermal and physical equilibrium with environment. In the engine applications, all the components of the working medium have to be either oxidized (e.g. fuel, CO, H), or reduced (e.g. NO, OH) in a reversible manner as the system approaches the dead state. The only components of the system, which cannot react chemically with the atmosphere are O_{2}, N_{2}, CO_{2} and H_{2}O [22, 23]. The thermo mechanical (Ex_{tm}) and chemical (Ex_{ch}) exergies are formulated as follows:

(1) 

(2) 
where m_{i}, x_{i}, and μ_{i0} are the mass flow, mass fraction, and chemical potential of species i obtained at restricted dead state. Here, μ_{i0} is chemical potential associated with restricted dead state, whereas superscript 0 as in, μ_{i}^{0} and x_{i}^{0} denote chemical potential and mass fraction associated with actual dead state condition for the species i. At the restricted dead state, there are thermal and mechanical equilibrium, but a chemical potential exists due to chemical reaction tendency with environment. At the actual dead state, this chemical potential is obviated. P_{0}, T_{0} are the pressure and temperature at dead state conditions, accordingly. The total exergy (Ex) is the sum of chemical and thermomechanical exergies:

(3) 
The rate of indicated work exergy is quantified as given:

(4) 
dV/dθ denotes volume change of cylinder volume with crank angle (CA) and P is the instantaneous pressure of the cylinder obtained by first law analyzed data. Ex_{q} represents the exergy of heat release within cylinder. The exergy of heat release rate (HRR) can be presented by (dQ/dθ is HRR):

(5) 
Irreversibility is denoted by I term that is the availability destruction concerning the combustion process in the chamber given as:

(6) 
The chemical availability of the C_{x}H_{y}O_{z} liquid fuel is calculated by the following formula [24]:

(7) 
where x, y, and z are the number of C, H, and O atoms in a general hydrocarbon fuel. LHV in the equation (7) defines the low heat value of the fuel. Note that w is the water content of the fuel. A_{f,ch} expresses the chemical availability (exergy) of a specific fuel.

(8) 
In order to quantify the second law efficiency of the engine the following formula is applicable, where A_{w} is the work transfer availability (exergy):

(9) 
Injection angle variation can influence the fuel flow behavior in the nozzle section due to the pressure and velocity field alteration. Subsequently, spray structure in terms of sauter mean diameter (SMD) and penetration was influenced; thereby air/fuel mixture uniformity can undergo drastic changes. SMD and liquid kinetic energy (injected spray jet energy) of the spray were taken into account as to evaluate the spray characteristics injected into the bowl of the combustion chamber. Figure 4 shows the liquid energy of injected spray at different injection angles with respect to crankangle under 1500 and 2500 rpm engine speeds. It is apparent that at 2500 rpm engine speed, higher liquid kinetic energy was accomplished. In fact, in most light duty diesel engines, smaller jet penetration is anticipated at higher speed since both swirl and tumble flows are stronger so that the axial momentum of the jet is lower while the radial one is larger. This phenomenon leads to a better atomization of sprays at higher engine speeds. No significant change in liquid energy can be noticed at 1500 rpm engine speed with injection angle variation, however, considering 2500 rpm, the spray direction of 150 ° shows the highest liquid energy with peak value of 0.1627 J. Although spray energy in the case of 140 ° is lower compared to the other injection angles, the prolonged liquid energy duration together with higher liquid energy after 7 °CA ATDC contributes to the highest spray penetration. Figure 5 demonstrates SMD of different sprays injected from different angles against crankangle. As shown in Figure 5, deviation of spray injection angle from bowl center towards squish and piston bottom wall (130 ° or 150 °) would increase spray droplet diameter. The overall higher liquid jet energy at 140 ° brings about intensified liquid jet disintegration.

Figure 4. Liquid kinetic energy as a function of CA (crank angle) for various injection angles under 1500, 2500 rpm speeds.


Figure 5. SMD (sauter mean diameter) evolution over CA (crank angle) for various injection angles under 2500 rpm speed.

Figure 6 illustrates equivalence ratio of air/fuel mixture as a function of crankangle for different spray injection angles at 1500 rpm engine speed. As fuel injection starts, fuel to air proportion and equivalence ratio starts to increase. According to Figure 6, lower equivalence ratio variance from stoichiometric ratio (~1) was achieved using 130 ° injection angle that signifies better mixing process, fuel evaporation, and more complete combustion. There are two reasons to explain this trend: first, it goes to the spray flow motion, which induces counter clockwise recirculation, therefore diffusing fuel into squish and top chamber wall. This makes better use of oxygen content in squish area. Another reason is attributed to reduced injectionwall distance and wall wetting for 150 ° spray angle. The uniform mixture distribution yields more obtainable work and heat generation, which will be addressed in subsequent sections.

Figure 6. Equivalence ratio as a function of CA (crank angle) for various injection angles at 2500 rpm speed.

Figure 7 depicts soot mass fraction distribution with crankangle for injection angles of 130 °, 140 °, and 150 ° at engine speeds of 1500 and 2500 rpm. Figure 7 indicates that higher soot concentration in combustion chamber is associated with low engine speed since higher engine speed increases the tumble flow motion and gas resistance in the chamber. Moreover, deflecting spray direction towards squish or pedestal edge demonstrates bigger droplet diameter and higher probability of spraywall impingement (see Fig. 5) that can increase soot substantially. The results of soot curves confirms the results of Figure 6 denoting that the highest soot mass fraction is attributed to the highest equivalence ratio for 140 ° case. According to Figure 8, the high intensity NOx concentration is produced with nozzle angle of 130 °, corresponding to the uniform air/fuel mixture within the piston volume (counter clockwise spray motion contributes in well mixing). The homogeneity of fuelair mixing is conducive to higher burning rate of fuel that can increase the incylinder temperature and pressure. NOx emission is expected to rise when the temperature increases. With engine speed increase, NOx emission decreases due to lesser air/fuel mixing time, thus lower evaporation, heat generation, and temperature is expected. Note that the lowest NOx for 1500 and 2500 rpm speeds attributes to 140 ° and 150 ° of the injection angles, respectively (according to maximum temperature presented in Table 3). Figure 9 clearly shows that in order to have the lowest NOx and soot emissions 140 ° and 150 ° of spray angle have to be adopted at 2500 rpm, respectively. This is mainly originated from the role of wallspray impingement and swirl motion of the flow in the swirl chamber. Irrespective of the spray injection angle, in all cases, 2500 rpm produced lesser emission magnitude.
deg  1500 rpm  2500 rpm  

T_{max} (K)  IMEP (bar)  ISFC (kg/kWh)  A_{f,ch} (J)  η_{II} (%)  Average (dI/dθ) (J/deg)  T_{max} (K)  IMEP (bar)  ISFC (kg/kWh)  η_{II} (%)  Average (dI/dθ) (J/deg)  
 
130  2074.97  15.4  0.385  1381.7  36.62  2.642  1951.18  12.5  0.540  34.03  4.166 
140  1815.57  12.1  0.574  1368.7  30.2  2.708  1695.98  10.4  0.721  28.12  2.887 
150  1911.72  13.3  0.491  1376.6  32.07  2.481  1822.77  11.2  0.639  31.12  2.746 

Figure 7. Soot mass fraction as a function of CA (crank angle) for various injection angles under 1500, 2500 rpm speeds.


Figure 8. NOx mass fraction as a function of CA (crank angle) for various injection angles under 1500, 2500 rpm speeds.


Figure 9. Emission comparison at different injection angles under 1500, 2500 rpm speeds.

Figure 10 shows accumulated work exergy as a function of CA (time) for various spray injection angles at specified 1500 and 2500 rpm speeds. The compression stroke up to 720 °CA demonstrates negative work exergy values and no perceptible difference in terms of work exergy can be observed between spraying orientation schemes. The accumulated work exergy for expansion and combustion period produce work and starts to increase up to 20 °ATDC, thereafter decreases accordingly. The increase in work exergy in lower engine speed (1500 rpm) is higher than that of higher engine speed. It was also demonstrated that 130 ° of spray injection angle generates higher work exergy amount due to higher uniformity of air/fuel mixture and air entrainment (closer equivalence ratio to 1), therefore higher evaporation and stronger incylinder pressure was achieved. Note that 130 ° injection has longer injector – wall distance that prevents fuel deposit on the bowl wall and availability destruction as a result of flame quenching and combustion deterioration. Although higher engine speed expedites spray breakup process as stronger tumble occurs in the chamber, due to irreversibility nature of higher speed and less mixing time the work availability content tends to decrease at high engine speed. The peak values of accumulative work exergy equals 495.1, 460.1 J corresponding to 130 ° of spray angle at 1500 and 2500 rpm, respectively. Higher engine speed reduces the systems capability to fully exploit the provided energy resulted in combustion process, thereby exergy destruction occurs. The results of Table 3 regarding IMEP (indicated mean effective pressure) are in accordance with the work exergy. The highest mean effective pressure is at 130 ° reaching to 15.4 and 12.5 bar at 1500 and 2500 rpm, respectively. With increasing engine speed from 1500 to 2500 rpm, 18.8%, 14.04%, and 15.78% IMEP reduction occurs while 69.8%, 63.1%, and 29.5% peak work exergy reduction happens for 130 °, 140 °, and 150 °, respectively. Capturing heat potential of combustion with increasing engine speed is dropped since increasing speed as a source of irreversibility decreases the heatwork transformation.

Figure 10. Accumulated work availability as a function of CA (crank angle) for various injection angles under 1500, 2500 rpm speeds.

Rate of heat release exergy was depicted in Figure 11 as a function of CA for various spray angles under 1500 and 2500 rpm engine speeds. As shown in Figure 11, longer heat release exergy pertains to low engine speed, meanwhile, heat exergy is delayed in the case of 2500 rpm. The highest and lowest peak values of heat release exergy are associated with 150 ° and 140 ° that happens at 2500 rpm engine speeds. Considering the highest second law efficiency and maximum temperature for the case of 130 °, it can be concluded that 150 ° is incapable of converting useful heat potential into work production, denoting much higher rate of the heat transfer from the chamber wall boundary. Nonetheless, 140 ° possesses the lowest temperature and HRR exergy due to high equivalence ratio and ISFC (indicated specific fuel consumption) and shortage of sufficient oxygen content for mixture evaporation and fuel burning. Increasing engine speed is a source of exergy destruction, therefore lower exergy content is found under higher engine speeds in all exergy forms. Applying 130 ° spray injection is an effective method for creating strong squish since it can induce upward spray flow motion and tumble along the compression stroke, hence powerful pressure was resulted. Higher pressure before combustion initiation makes stronger turbulence and homogeneously distributed fuel air mixture. With upward movement of the spray injection angle, combustion site gets closer to squish zone, hence the majority of the released heat is trapped within the squish zone giving rise to deficient combustion. Increasing temperature and pressure resulted from high quality of spray structure (indepth penetration of spray and near stoichiometric equivalence ratio) and combustion process of 130 ° spray injection into the bowl section of piston. Higher temperature and pressure amplitude of the combustion chamber for the case of 130 ° of injection angle explains higher overall availability, according to Figure 12. One can find that the total availability rise was delayed a few CAs for 130 ° and 150 ° which implies more opportunity for effective air/fuel mixing before combustion initiates. Consequently, exergy peak values of 130 ° and 150 ° injection angle are 29.8%, 14.6% higher than 140 ° under 1500 rpm, respectively. Table 3 reveals higher maximum temperature and effective pressure (IMEP) when 130 ° of injection angle was applied. As a result, higher total availability for the case of 130 ° once again was proved since availability is correlated with temperature and pressure of a system. Total chemical availability is depicted in Figure 13 with respect to CA for various injection angle and engine speed schemes. Total chemical availability is consisted of fuel chemical availability (reduction or oxidation) and diffusivity availability. As offered by Figure 13, a monotonous trend is presented during compression stroke for chemical availability since no fuel was injected and chemical reaction would not happen. The mole fraction increase and dramatic pressure and temperature rise with combustion initiation induce exergy increase. With regard to 130 °, more fuel was involved in chemical reaction due to fuel diffusion in squish area and higher oxygen content accessibility. However, the rate of chemical exergy increase was reduced after 740 °CA due to combustion chamber expansion and reduced diffusivity as mole concentration of combustion products is stabilized. 130 ° of injection angle prevails in total chemical availability than those of other cases that stems from better combustion implementation and effectual fuel burning reaction, together with higher fuel diffusion with residual gas. Second law efficiency is represented in detail in Table 3, where 130 ° showed superior amount equal to 36.62%, 34.03% under 1500, and 2500 rpm speeds, respectively. On the other hand, the lowest second law efficiency is associated with 140 ° such that 30.2%, 28.12% was achieved regarding 1500 and 2500 rpm, respectively. The results of second law efficiency are consistent with IMEP, ISFC, and T_{max} of each case. Irreversibility rate as a function of CA for different spray aiming scheme was depicted in Figure 14 under 1500 rpm engine speed. Based on Table 3 and Figure 14, mean irreversibility rate per CA for 140 ° is dominant (2.708 J/deg over CA domain) while instantaneous peak irreversibility is dominant for 130 ° at 724 °CA (17.48 J/deg). This is because the system cannot seize the whole heat generated over crankangle period. The mean irreversibility shows the highest amount for the case of 140 ° under both 1500 and 2500 rpm.

Figure 11. HRR (heat release rate) availability as a function of CA (crank angle) for various injection angles under 1500, 2500 rpm speeds.


Figure 12. Total availability as a function of CA (crank angle) for various injection angles under 1500, 2500 rpm speeds.


Figure 13. Chemical availability as a function of CA (crank angle) for various injection angles under 1500 and 2500 rpm speeds.


Figure 14. Irreversibility rate as a function of CA (crank angle) for various injection angles at 1500 rpm.

A correlation between total availability (A) and CA was developed based on the presented equation (10) (see Table 4). The best curve that data points can map onto is attained with fourth order polynomial. A linear regression analysis for 130 °, 140 °, and 150 ° injection orientation was performed for 1500 and 2500 rpm speed scenarios. Based on the obtained data for a_{0}, a_{1}, a_{2}, a_{3}, and a_{4} coefficients for the best curve fitting, the most precise formula according to correlation coefficient is presented for 130 deg1500 rpm case with R^{2} = 91.93%. However, the lowest RMSE (root mean square error) is the lowest for 140 ° of injection angle to the bowl segment. The formula is presented as follows:

(10) 
1500 rpm  2500 rpm  

a_{0} ×(10^{7})  a_{1} (10^{5})  a_{2} ×(10^{2})  a_{3} ×(10^{−1})  a_{4} ×(10^{−4})  RMSE –  R^{2} (%)  a_{0} ×(10^{7})  a_{1} ×(10^{5})  a_{2} ×(10^{2})  a_{3} ×(10^{−1})  a_{4} ×(10^{−5})  RMSE –  R^{2} (%)  
 
130  4.06  −2.2  4.46  −4.02  1.35  3.43  91.93  1.94  −1.04  2.09  −1.86  6.21  3.25  91.78 
140  3.05  −1.65  3.36  −3.03  1.02  2.5  91.86  2.37  −1.28  2.61  −2.36  7.98  2.2  89.38 
150  3.19  −1.73  3.51  −3.51  1.06  2.79  91.87  1.67  −2.33  1.8  −1.61  5.39  2.87  91.1 
The performance and accuracy of developed correlation is gauged by statistical factors such as coefficient of determination (R^{2}), RMSE, defined as:

(11) 

(12) 

(13) 
where Y_{actual} and Y_{predicted} are computed and predicted values of objective parameters by the developed models, respectively.
Energy analysis was performed along exergy assessment for different injection strategy under different engine speeds. A correlation was observed between the fluid flow and the combustion process with exergy terms. The results can be classified as follows:
None declared.
Published on 01/06/17
Submitted on 01/06/17
Licence: Other
Are you one of the authors of this document?