Impact of Spray Cone Angle on the Performances of Methane/Diesel RCCI Engine Combustion under Low Load Operating Conditions

The behaviors of spray, in Reactivity Controlled Combustion Ignition (RCCI) dual fuel engine and subsequent emissions formation, are numerically addressed. Five spray cone angles ranging between 5° and 25° with an advanced injection timing of 22° Before Top Dead Center (BTDC) are considered. The objective of this paper is twofold: (a) to enhance engine behaviors in terms of performances and consequent emissions by adjusting spray cone angle and (b) to outcome the exergy efficiency for each case. The simulations are conducted using the Ansys-forte tool. The turbulence model is the Renormalization Group (RNG) K-epsilon, which is selected for its effectiveness in strongly sheared flows. The spray breakup is governed by the hybrid model Kelvin–Helmholtz and Rayleigh–Taylor spray models. A surrogate of n-heptane, which contains 425 species and 3128 reactions, is used for diesel combustion modeling. The obtained results for methane/diesel engine combustion, under low load operating conditions, include the distribution of heat transfer flux, pressure, temperature, Heat Release Rate (HRR), and Sauter Mean Diameter (SMD). An exergy balance analysis is conducted to quantify the engine performances. Output emissions at the outlet of the combustion chamber are also monitored in this work. Investigations show a pressure decrease for a cone angle θ = 5° of roughly 8%, compared to experimental measurement (θ = 10°). A broader cone angle produces a higher mass of NOx. The optimum spray cone angle, in terms of exergy efficiency, performance, and consequent emissions is found to lie at 15° ≤ θ ≤ 20°.


Introduction
The increase in greenhouse gases emissions demands a deeper encouragement of combustion research. Since it is heavily related to fossil fuel combustion, the transport sector has an important part in global warming and climate change [1,2]. The depletion of fossil fuels requires, also, an in-depth strategy in fuel economy [3]. Therefore, further improvement within Compression Ignition (CI) engine is always recommended. The combustion process dealing with compression ignition, specific power output, and fuel consumption, should be substantially altered. Reducing emissions is, also, required for the CI engine. This opens a new path for novel/alternative fuels [4].
To improve engine performances, various techniques were developed during the last decade. Low-Temperature Combustion (LTC), thanks to its capability of low-temperature combustion, is widely used. It serves to reduce soot, Carbon Oxides (CO), and Nitrogen hybrid Eulerian-Lagrangian Kelvin Helmholtz-Rayleigh Taylor (KH-RT) spray breakup model describes the spray atomization. The combined models stand for a detailed numerical investigation of the spray cone angle adjusting for combustion improvement of the CH 4 /diesel RCCI engine. The alternative fuel is methane. It is injected, as premixed, with the oxidizer. n-heptane represents diesel. The modeled engine is a single cylinder four stroke that operates under low load. The present work examines drop atomization by retrofitting spray cone angle. Five spray angles θ = 5 • , 10 • , 15 • , 20 • , and 25 • , are studied. The Forte software linked to the Ansys-CHEMKIN library is employed for the CFD calculation. The heat transfer flux, in-cylinder temperature, Sauter Mean Diameter (D 32 ), pressure, and Heat Release Rate (HRR) are studied. An exergy balance analysis is conducted to investigate the RCCI performances. Output emissions at Exhaust Valve Opening (EVO) are also reported.
The objective of this paper is twofold: (a) to enhance RCCI engine behaviors in terms of performances and consequent emissions by adjusting the spray cone angle and (b) to control the exergy efficiency for each case.
The present paper is organized as follows. The next section, Section 2, is adopted for the numerical method, which holds experimental configuration, presentation of the reaction mechanism for diesel surrogate, mesh and boundary condition followed by the governing equations and exergy analysis. The numerical validation is also reported in this section. The results, including heat transfer flux, pressure, HRR, Weber number, SMD, exergy efficiency, and output species are presented and discussed in Section 3. Section 4 summarizes the work in a conclusion.

Experimental Configuration
A numerical study was performed at an engine speed of 900 rpm and under 25% of engine load. The modeled engine operated with methane/diesel dual fuel strategy. It runs under lean equivalence ratio, ER = 0.41. The fuel composition contained 99% premixed airmethane mass fraction, while 1% was direct injection diesel. The Exhaust Gas Recirculation (EGR) rate equaled zero. To predict engine performances, five spray cone angles were employed. The assumed cone angles were 5 • , 10 • , 15 • , 20 • and 25 • . Both chemistry and physics, were resolved using the CFD-code Ansys-Forte. The present study utilizes a threedimensional model. Since such a configuration is axisymmetric, the modeled combustion chamber is set to a 60 • periodic sector. The main geometry components are highlighted in Figure 1. The modeled engine was a single cylinder four stroke as given by Yousfi et al. [10]. The main engine specifications are listed in Table 1.

Diesel Surrogate for the Reaction Mechanism
A unique component, N-heptane (C7H16), surrogate was applied for diesel modeling. The physical properties of the diesel spray and vaporization were represented by n-tetradecane (C14H30). The reaction mechanism contained 425 total species and 3128 reactions

Mesh and Boundary Conditions
Forte software, based on the finite volume method, was chosen for transport equations resolving. Table 2 holds the required boundary and initial conditions. A structured quadratic mesh was generated ( Figure 2). The regular mesh allows easy data management during computation. Therefore, data were smoothly transferred between cells. Furthermore, dynamic mesh made using a layering technique was fast and preserves good quality mesh. The computational model had roughly 13,500 cells at IVC. A similar number of control volumes were used by Sage L. Kokjohn and Rolf D. Reitz [30]. They showed that results were grid-independent. The turbulence K-epsilon RNG was selected for its effectiveness in strongly sheared flows and, relatively, reduced computational costs. It is also wellsuited for governing the turbulence in non-isotropic combustion. The mesh and timestepindependent spray breakup model, called the gas jet model, was applied. The atomization of sprays was governed using the hybrid Kelvin-Helmholtz Rayleigh-Taylor (KH-RT) breakup model.

Governing Equations
To describe the spray dynamics, thermodynamics, and chemistry properties as well as the turbulent multiphase flows evolving in this configuration, the mass, momentum, energy, and species transport equations were solved following an Eulerian framework while the droplets were tracked within a Lagrangian approach. The mass conservation for each control volume is expressed by the following equation:

Governing Equations
To describe the spray dynamics, thermodynamics, and chemistry properties as well as the turbulent multiphase flows evolving in this configuration, the mass, momentum, energy, and species transport equations were solved following an Eulerian framework while the droplets were tracked within a Lagrangian approach. The mass conservation for each control volume is expressed by the following equation: where ρ is the fluid density, u represents the mean velocity vector, and . ρ s is the source term due to the spray evaporation. Note that the overbar represents Reynolds averaging and the tilde is Favre averaging.
The momentum equation, which represents the motion of the fluid, is given as: where, p, σ, F s , and g are the in-cylinder pressure, viscous shear stress, spray-induced source term, and specific body force, respectively. Γ is the Reynolds stress tensor.
To calculate the effects of heat transfer, turbulent transport, turbulent dissipation, and chemical reactions, the internal energy equations have to be considered.
where, I is the Internal energy, p represents the pressure, J is derived from the total heat flux from enthalpy diffusion and heat conduction.
. Q S and . Q C represent the source terms generated during diesel injection and chemical reactions, respectively. The quantity H accounts for the effects of filtered convection term, while ε is the dissipation rate.
The K-epsilon RNG turbulence model was derived from the momentum equations using a mathematical approach called the "Renormalization group" [31]. Both k and ε, which represent the turbulent kinetic energy and the dissipation rate, respectively, are calculated as follows: In both equations Pr k , Pr ε , C s , C ε1 , C ε2 and C ε3 represent the model coefficients. They are summarized in Table 3. The quantities µ and µ T are the renormalized and the turbulence is the source term due to spray vaporization and Γ is the Reynolds stress tensor. The terms υ and υ T represent the renormalized and the turbulent kinematic viscosity, respectively. R is related to the strain tensor as presented in [31]. Table 3. Constants in the RNG k-epsilon model [32]. The gas phase, in the combustion engine, was modeled as a mixture of either gas components or species. This composition varies along with the engine cycle due to molecular diffusion, flow convection, and turbulent transport interaction during the combustion processes. The transport equation for the species, k, is represented as follows The term ρ represents the density and the subscript k denotes the species index, the capital term K denotes the number of species, → u is the velocity vector of the flow, and the quantity y k = ρ k /ρ represents the mass fraction of the transported species, k, while the term D is the mixture-averaged molecular diffusion coefficient. Note that the term Φ accounts for the effect of ensemble-averaging of the convection term, which equals Φ = ρ k u − ρ k u. The quantities . ρ k S and . ρ k C are source terms due to spray vaporization and chemical reactions, respectively.

Diesel Spray Model
A solid cone diesel spray model is chosen. The hybrid Kelvin-Helmholtz and Rayleigh-Taylor spray models govern the drop atomization. The details of these models are provided in [34]. For the self-consistency of the paper, a concise description is outlined here. The spray was processed through three steps: first, thicker film formation at the tip of the jet; second, primary breakup; and later, droplets atomization.
In the near nozzle, the liquid was assumed to be a dense blob core and the KH instability tracked the primary breakup out of the jet. The KH model proposes a parent parcel with radius, r, which equals the nozzle diameter. The child droplets with radius, r c , which is formed during the primary breakup, are modeled with: where, B 0 equals 0.61 and Λ KH represents the KH wavelength of the accurate growing wave, Ω KH . The growth rate of the fastest wave and its related wavelengths are assumed as: Here, r p is the radius of jet, σ is the surface tension, ρ l represents the density of the injected liquid, and T, Z, and We g are the dimensionless Taylor, Ohnsorge, and Weber numbers of gas, respectively. During the breakup, a new child droplet is formed due to mass loss of the parent parcel. The radius of the new droplet is calculated using Equation (10): where, τ KH represents the breakup time.
To predict the secondary break up, The RT model ensures the calculation of the frequency and wavelength of the rapid trained growing wave using: The quantity a represents the deceleration owing to drag force. The quantity ρ g is the gas density. The term Λ RT denotes the wave length of the accurate growing wave, while Ω RT represents the growth rate during the secondary breakup. The formation of the child droplets, r c , and the breakup time, τ RT , are assumed, respectively, as: B RT = 40 and C RT = 0.1 represent the size constant and the time constant of the RT breakup.
Ex in is the inlet exergy rate, which depends on the air and fuel exergies rates as written in Equation (16), Ex air + . Ex f uel (16) .
Ex out is the sum of exergy rates in form of work ( . Ex W ), exhaust gases . Ex exh , and heat transferred . Ex heat , as provided in Equation (17), Equations (15)- (17) can be put together and expressed as follows: .
Ex des (18) In this case, the system is adiabatic ( Ex des is calculated through entropy generation according to the second law of thermodynamics. Ex des are assumed as the losses quantities ( . Ex loss ) and they are written in Equation (19) . Based on Equations (18) and (19), the exergy efficiency η ex can be deduced as:
Since the intake air is atmospheric, the exergy rate of air is neglected ( . Ex air = 0). Equation (23) outlines that the exergetic work rate represents the net work [36] .
In this work, the engine operates under methane/diesel dual fuel conditions. The .
Ex f uel is calculated by adding the exergy of each fuel [37]. For diesel-methane as fuels, the input chemical exergy rate is computed by using Equation (24) as: .
For hydrocarbon fuels, which is the case in this work (N-C 7 H 16 and CH 4 ), the input fuel exergy rate is generally expressed as [38,39]: .
The parameters used for exergy efficiency (η ex ) calculation are tabulated in Tables 4 and 5.

Model Validation
Prior to the various parameter studies carried out and reported in this work, a numerical validation was conducted with the experimental results of Yousfi et al. [10]. The experimental results were obtained using a single-cylinder caterpillar 4100 engine. Details are listed in Table 2. The experiment focused on the effect of the Start Of Injection (SOI). Based on a wide range of SOI, from 10 • to 30 • BTDC, data explore that the in-cylinder pressure increases by making the injection timing earlier. Early injection also shifts the peak pressure near the Top Dead Center (TDC). Here, the validation covered both in-cylinder pressure and heat release rate. The numerical simulation took place from Intake Valve start of injection (SOI) took place at 22° BTDC. The comparison between them is pointed out in Figure 3.
Using the same operating conditions, the numerical results were in good agreement with the experimental measurements. For HRR, an error of 38% was recorded, at CA = 360°. This important error was most probably caused by the combustion mixing and important variance of the local equivalence ratio. The lean/rich zones worsen flame propagation, therefore, most of the chemical reactions were not complete, thus, leading to reduced HRR.

Results and Discussion
The simulations run from IVC to EVO. The chemistry calculation was activated from 315° CA to 400° CA. After this range, chemistry was found negligible, similar to the result in [28,33]. The maximum time step was set to 10 −5 s, while the initial simulation time step equaled 10 −7 s.
The obtained results under low load operating conditions include the distribution of heat transfer flux, pressure, temperature, Heat Release Rate (HRR), and Sauter Mean Diameter (SMD). An exergy analysis, together with the RCCI performances, is provided. Finally, the results of the emissions, captured at Exhaust Valve Opening (EVO) are presented. Using the same operating conditions, the numerical results were in good agreement with the experimental measurements. For HRR, an error of 38% was recorded, at CA = 360 • . This important error was most probably caused by the combustion mixing and important variance of the local equivalence ratio. The lean/rich zones worsen flame propagation, therefore, most of the chemical reactions were not complete, thus, leading to reduced HRR.

Results and Discussion
The simulations run from IVC to EVO. The chemistry calculation was activated from 315 • CA to 400 • CA. After this range, chemistry was found negligible, similar to the result in [28,33]. The maximum time step was set to 10 −5 s, while the initial simulation time step equaled 10 −7 s.
The obtained results under low load operating conditions include the distribution of heat transfer flux, pressure, temperature, Heat Release Rate (HRR), and Sauter Mean Diameter (SMD). An exergy analysis, together with the RCCI performances, is provided. Finally, the results of the emissions, captured at Exhaust Valve Opening (EVO) are presented.

Distribution of Heat Transfer Flux
The heat transfer flux varies considerably with droplet loading, position as well as surface covered. Heat transfer fluxes, for various spray angles, are outlined in Figure 4. The heat transfer flux affects the start of ignition, in conjunction with the ignition duration. At θ = 5 • , the start of combustion is delayed to TDC. From θ = 10 • onwards, the combustion starts earlier at θ = −5 • ATDC and the combustion duration expands as long as the spray angle increases. The maximum heat flux is registered at θ = 20 • , which means that the stratification of the mixture at the mentioned angle generates better combustion.
The heat transfer flux varies considerably with droplet loading, position as well as surface covered. Heat transfer fluxes, for various spray angles, are outlined in Figure 4. The heat transfer flux affects the start of ignition, in conjunction with the ignition duration. At θ = 5°, the start of combustion is delayed to TDC. From θ = 10° onwards, the combustion starts earlier at θ = −5° ATDC and the combustion duration expands as long as the spray angle increases. The maximum heat flux is registered at θ = 20°, which means that the stratification of the mixture at the mentioned angle generates better combustion.

Temperature
It can be observed, in Figure 5, that a significant variation in the in-cylinder temperature occurs for the different injector angles. The temperature and/or rich-mixture are, in particular, the main sources of nitrogen oxide formation. The maximum temperature is one of the most critical parameters during IC engine combustion. Data show a temperature variation between 1561 K and 1766 K. The temperature peak value is recorded for θ = 15 • and θ = 20 • . Compared to the experiment, the temperature increases by almost 6%. A longer combustion duration is registered by θ = 15 • and θ = 20 • . For cone angles smaller than 15 • , the combustion duration is reduced. A narrow spray angle affects the evaporation, due to the dense liquid blob. Therefore, the flame propagation needs more residence time to outbreak. Low temperature is a result of high-unburned hydrocarbon, thus a considerable CO emission is obtained. The subsequent increase in temperature, at CA = 352 • is thought to initiate a flame outbreak across the methane region of the combustion chamber. Methanediesel flame temperature is lower compared to that of a gasoline-diesel flame, resulting in lower NO x formation [40]. For θ = 15 • and θ = 20 • , a high in-cylinder temperature close to TDC is registered, generating considerable NO x . The same results were obtained by Poorghasemi et al. [12].
one of the most critical parameters during IC engine combustion. Data show a temperature variation between 1561 K and 1766 K. The temperature peak value is recorded for θ = 15° and θ = 20°. Compared to the experiment, the temperature increases by almost 6%. A longer combustion duration is registered by θ = 15° and θ = 20°. For cone angles smaller than 15°, the combustion duration is reduced. A narrow spray angle affects the evaporation, due to the dense liquid blob. Therefore, the flame propagation needs more residence time to outbreak. Low temperature is a result of high-unburned hydrocarbon, thus a considerable CO emission is obtained. The subsequent increase in temperature, at CA = 352° is thought to initiate a flame outbreak across the methane region of the combustion chamber. Methane-diesel flame temperature is lower compared to that of a gasoline-diesel flame, resulting in lower NOx formation [40]. For θ = 15° and θ = 20°, a high in-cylinder temperature close to TDC is registered, generating considerable NOx. The same results were obtained by Poorghasemi et al. [12].

Pressure
The impingement of fuel, injected over the cylinder liner, is a challenge for flame homogeneity and hence, for IC engine improvement. Figure 6 outlines the pressure variation for five spray cone angles. The pressure peak values for the five spray cone angles equal 65 bar, 71 bar, 79 bar, 79.5 bar, and 78 bar for θ = 5°, 10°, 15°, 20°, and 25°, respectively. It is observable that the pressure increases as the spray angle increases from 5° to 20°, however, it slightly decreases for θ =25°. This indicates that the effect of the spray angle for θ > 20° becomes weaker. The RCCI combustion duration, for θ = 15°, 20°, and 25° is longer than θ = 5° and θ = 10°. This could be explained as follows. For θ = 5° and θ = 10°, the liquid jet remains dense until impinging against the piston, which decreases the expansion work, and therefore, reduces chemical reaction rates. Here, the CO2 specie mass, which is a good indicator of a complete combustion development, is reduced. For 15° ≤ θ ≤ 20°, the pressure is extremely high because the injection takes place around the center of the piston bowl. Similar results were found by Balijepalli et al. [41]. The droplets interact with the piston bowl. The heat and mass transfer are improved, therefore, the combustion process is accelerated, which is consistent with the We-Number results.

Pressure
The impingement of fuel, injected over the cylinder liner, is a challenge for flame homogeneity and hence, for IC engine improvement. Figure 6 outlines the pressure variation for five spray cone angles. The pressure peak values for the five spray cone angles equal 65 bar, 71 bar, 79 bar, 79.5 bar, and 78 bar for θ = 5 • , 10 • , 15 • , 20 • , and 25 • , respectively. It is observable that the pressure increases as the spray angle increases from 5 • to 20 • , however, it slightly decreases for θ =25 • . This indicates that the effect of the spray angle for θ > 20 • becomes weaker. The RCCI combustion duration, for θ = 15 • , 20 • , and 25 • is longer than θ = 5 • and θ = 10 • . This could be explained as follows. For θ = 5 • and θ = 10 • , the liquid jet remains dense until impinging against the piston, which decreases the expansion work, and therefore, reduces chemical reaction rates. Here, the CO 2 specie mass, which is a good indicator of a complete combustion development, is reduced. For 15 • ≤ θ ≤ 20 • , the pressure is extremely high because the injection takes place around the center of the piston bowl. Similar results were found by Balijepalli et al. [41]. The droplets interact with the piston bowl. The heat and mass transfer are improved, therefore, the combustion process is accelerated, which is consistent with the We-Number results.
The best spray angle is θ = 20 • . Here, the pressure increases roughly 11% compared to experimental results (θ = 10 • ). It has a positive effect on air-methane/diesel mixing and on the start of combustion compared to other cone spray-angles. The best spray angle is θ = 20°. Here, the pressure increases roughly 11% compared to experimental results (θ = 10°). It has a positive effect on air-methane/diesel mixing and on the start of combustion compared to other cone spray-angles.

Heat Release Rate (HRR)
The HRR is an important parameter to pursue the combustion stratification, and therefore it is used as a parameter for thermal efficiency. The curves, in Figure 7, are calculated based on the specific heat ratio, pressure, and volume variation over the crank angle. The peak value increases proportionally to the spray angle and reaches the maximum at 352° CA. The heat release rate curves are divided into two different zones. First, the rate of heat release rises to the maximum at CA = 352°. Then, a second peak, ranging between CA = 356-362°, occurs at a lean equivalence ratio. At the combustion stroke, the energy is released and the combustion slightly drops off and carries on at a constant level [42]. Thus, RCCI produces a long combustion duration. This is due to a large difference in fuel component volatility, which results in a sequence and long duration of auto-combustion [7]. The HRR pattern shows that the combustion phasing is changed when diesel is injected into a methane-air environment. Dual fuels with various reactivities for RCCI are denoted to control combustion phasing (CP) and HRR in the engine [19,43]. The HRRvalue shows a minimum at TDC. This is understandable as the burned charge is fully released at 351°, far away from the TDC. It is worth noting that the spray cone angle decreases as the combustion phasing is retarded, resulting in more HRR at TDC.

Heat Release Rate (HRR)
The HRR is an important parameter to pursue the combustion stratification, and therefore it is used as a parameter for thermal efficiency. The curves, in Figure 7, are calculated based on the specific heat ratio, pressure, and volume variation over the crank angle. The peak value increases proportionally to the spray angle and reaches the maximum at 352 • CA. The heat release rate curves are divided into two different zones. First, the rate of heat release rises to the maximum at CA = 352 • . Then, a second peak, ranging between CA = 356-362 • , occurs at a lean equivalence ratio. At the combustion stroke, the energy is released and the combustion slightly drops off and carries on at a constant level [42]. Thus, RCCI produces a long combustion duration. This is due to a large difference in fuel component volatility, which results in a sequence and long duration of auto-combustion [7]. The HRR pattern shows that the combustion phasing is changed when diesel is injected into a methane-air environment. Dual fuels with various reactivities for RCCI are denoted to control combustion phasing (CP) and HRR in the engine [19,43]. The HRR-value shows a minimum at TDC. This is understandable as the burned charge is fully released at 351 • , far away from the TDC. It is worth noting that the spray cone angle decreases as the combustion phasing is retarded, resulting in more HRR at TDC.
The best spray angle is θ = 20°. Here, the pressure increases roughly 11% compared to experimental results (θ = 10°). It has a positive effect on air-methane/diesel mixing and on the start of combustion compared to other cone spray-angles.

Heat Release Rate (HRR)
The HRR is an important parameter to pursue the combustion stratification, and therefore it is used as a parameter for thermal efficiency. The curves, in Figure 7, are calculated based on the specific heat ratio, pressure, and volume variation over the crank angle. The peak value increases proportionally to the spray angle and reaches the maximum at 352° CA. The heat release rate curves are divided into two different zones. First, the rate of heat release rises to the maximum at CA = 352°. Then, a second peak, ranging between CA = 356-362°, occurs at a lean equivalence ratio. At the combustion stroke, the energy is released and the combustion slightly drops off and carries on at a constant level [42]. Thus, RCCI produces a long combustion duration. This is due to a large difference in fuel component volatility, which results in a sequence and long duration of auto-combustion [7]. The HRR pattern shows that the combustion phasing is changed when diesel is injected into a methane-air environment. Dual fuels with various reactivities for RCCI are denoted to control combustion phasing (CP) and HRR in the engine [19,43]. The HRRvalue shows a minimum at TDC. This is understandable as the burned charge is fully released at 351°, far away from the TDC. It is worth noting that the spray cone angle decreases as the combustion phasing is retarded, resulting in more HRR at TDC.

Weber Number (We)
The Weber number is a relevant parameter of spray atomization and therefore for droplet vaporization. It indicates whether the surface tension or the kinetic energy is dominant. Figure 8 denotes the weber number, for various spray cone angles, with respect to engine crank angles. The narrower spray cone angle (θ = 5 • ), shows a lower weber number. This indicates that the surface tensions of the gas mixture are dominated by the injected liquid inertia forces. Thus, the break-up time becomes larger and the combustion time increases. A higher We-Number is registered for θ = 15 • . It equals 1378, which increases by 12% compared to the experimental work (θ = 10 • ). A second peak is registered (with a lower We-Number), for various spray cone angle, at CA = 351 • . It is due to the sprays' (droplets') impingement against the piston bowl as reported by J.D Naber [44]. The data are consistent with the obtained HRR value.

Weber Number (We)
The Weber number is a relevant parameter of spray atomization and therefore for droplet vaporization. It indicates whether the surface tension or the kinetic energy is dominant. Figure 8 denotes the weber number, for various spray cone angles, with respect to engine crank angles. The narrower spray cone angle (θ = 5°), shows a lower weber number. This indicates that the surface tensions of the gas mixture are dominated by the injected liquid inertia forces. Thus, the break-up time becomes larger and the combustion time increases. A higher We-Number is registered for θ = 15°. It equals 1378, which increases by 12% compared to the experimental work (θ = 10°). A second peak is registered (with a lower We-Number), for various spray cone angle, at CA = 351°. It is due to the sprays' (droplets') impingement against the piston bowl as reported by J.D Naber [44]. The data are consistent with the obtained HRR value.

Sauter Mean Diameter (D32)
The Sauter mean diameter is sought to discover the droplet volume covered by the available surface. It is important for mass transfer, and therefore, combustion efficiency. Figure 9 exhibits the size distribution of droplets during the injection phase. D32 decreases considerably with increasing spray angle. The curves indicate that a broad spray angle decreases the collision between droplets. A smaller SMD (D32) promotes liquid-fuel vaporization, thus better and faster mixing is achieved. The droplet diameter increases for θ = 25°. Here, the higher droplet diameter impedes the start of combustion. This is remarkably observed by HRR results (Figure 7). The minimum D32 value is registered for θ = 20°. It equals 85.3 microns. It is 37% smaller compared to the 5° and 10° spray cone angles and 9% smaller than the 25° and 15° results. The obtained result agrees with outlined temperature and pressure behavior. Maximum pressure and temperature are also registered by θ = 20°.

Sauter Mean Diameter (D 32 )
The Sauter mean diameter is sought to discover the droplet volume covered by the available surface. It is important for mass transfer, and therefore, combustion efficiency. Figure 9 exhibits the size distribution of droplets during the injection phase. D 32 decreases considerably with increasing spray angle. The curves indicate that a broad spray angle decreases the collision between droplets. A smaller SMD (D 32 ) promotes liquid-fuel vaporization, thus better and faster mixing is achieved. The droplet diameter increases for θ = 25 • . Here, the higher droplet diameter impedes the start of combustion. This is remarkably observed by HRR results (Figure 7). The minimum D 32 value is registered for θ = 20 • . It equals 85.3 microns. It is 37% smaller compared to the 5 • and 10 • spray cone angles and 9% smaller than the 25 • and 15 • results. The obtained result agrees with outlined temperature and pressure behavior. Maximum pressure and temperature are also registered by θ = 20 • .

Exergy Efficiency
To investigate the interaction between diesel fuel consumption and output power generation, an exergetic analysis is conducted. The engine performances are studied, for various spray angles. The exergy efficiency, as a function of different spray cone angle, is presented in Figure 9. D32 over crank angle for different spray cone−angles.

Exergy Efficiency
To investigate the interaction between diesel fuel consumption and output power generation, an exergetic analysis is conducted. The engine performances are studied, for various spray angles. The exergy efficiency, as a function of different spray cone angle, is presented in Figure 10. The exergy efficiency for θ = 15 • and θ = 20 • is roughly 38%. It increases by 5% compared to the experimental work (θ = 10 • ). This is an indication that combustion occurs in better performances. It is worth noticing that a 33% exergy efficiency is recorded for the narrower spray angle (θ = 5 • ). Here, 8% exergy efficiency is lost compared to experimental results (θ = 10 • ). It is explained by the low temperature, as shown in Figure 5.

Exergy Efficiency
To investigate the interaction between diesel fuel consumption and output generation, an exergetic analysis is conducted. The engine performances are stud various spray angles. The exergy efficiency, as a function of different spray cone a presented in . The exergy efficiency for θ = 15° and θ = 20° is roughly 38%. It increases compared to the experimental work (θ = 10°). This is an indication that combustion in better performances. It is worth noticing that a 33% exergy efficiency is recorded narrower spray angle (θ = 5°). Here, 8% exergy efficiency is lost compared to experi results (θ = 10°). It is explained by the low temperature, as shown in Figure 5.

Output Species
To control the produced species during the combustion process, the masses CO, and EINOx are postprocessed at the exit of the combustion chamber. Error! Ref source not found. points out these species for the various spray cone angles.

Output Species
To control the produced species during the combustion process, the masses of CO 2 , CO, and EINO x are postprocessed at the exit of the combustion chamber. Figure 11 points out these species for the various spray cone angles. Error! Reference source not found. shows the wide change of EINOx, CO, and CO2 mass when changing theta. For 5° < θ < 20°, Emission Index (EI) NOX, which represents the mass of NOX converted by kilogram of fuel consumed, increases by 49%, 173%, 227%, and 250% as the spray angle increases. This quantity is useful to quantify the flame behavior, as compared with others flame types, for example [45]. Since NOx is strongly correlated to the temperature, this indicates that the Zeldovich NOx mechanism is the predominant way for NOx production. For θ > 20°, even though the temperature decreases, a considerable NOx formation is obtained. This is most probably due to the prompt NOx formation mechanism, which occurs at lower temperature in fuel rich regions. Compared to θ = 10°, with which validation took place, NOx species remarkably increase up to 49%, 54%, and 57% at 15°, 20°, and 25°, respectively.
The CO formation of methane/diesel RCCI engine is outlined in Error! Reference source not found.. CO emissions decrease with increasing spray angle, between θ = 5° and θ = 15°. For θ > 15°, CO slightly increases. A higher percentage of CO is observed for a narrower injector cone angle, θ = 5°. Broader cone angles lead to better vaporization, therefore, better mixing and combustion. The obtained results are conformal with the pressure registered ( Figure 6). As shown in Figure 5, the temperature at θ = 5° is the lowest, this indicates that combustion was incomplete.
The formation of CO2 depends on the combustion behavior and the total injected mass of fuel [46]. The total diesel mass injected is 0.01649 g. The CO2 specie shows an important mass at θ = 15°. This indicates an excellent evaporation characteristic obtained. For the case of θ = 5°, due to poor mixing and the presence in a rich region, CO could not be converted into CO2, that is why a lower CO2 concentration is formed.
To evaluate engine performances in terms of output emissions and the overall performance for the different spray angle, combustion at θ = 15° is best achieved leading with the lowest CO.

Conclusions
The present study explored the impact of different spray cone angles on the combustion and exergy efficiencies in the Internal Combustion Engine (ICE). The experimental data, from the literature, were performed on a modified single-cylinder caterpillar 3400 heavy-duty diesel engine. The numerical work is conducted on a single cylinder fourstroke engine, using the dual fuel RCCI strategy. The KH-RT models the solid cone sprays and the diesel jet atomization. The assessed cone angles are 5°, 10°, 15°, 20° and 25°. A Figure 11. Alteration of EINOx, CO, and CO 2 emissions for different spray cone−angles. Figure 11 shows the wide change of EINO x , CO, and CO 2 mass when changing theta. For 5 • < θ < 20 • , Emission Index (EI) NO X , which represents the mass of NO X converted by kilogram of fuel consumed, increases by 49%, 173%, 227%, and 250% as the spray angle increases. This quantity is useful to quantify the flame behavior, as compared with others flame types, for example [45]. Since NO x is strongly correlated to the temperature, this indicates that the Zeldovich NO x mechanism is the predominant way for NO x production. For θ > 20 • , even though the temperature decreases, a considerable NO x formation is obtained. This is most probably due to the prompt NO x formation mechanism, which occurs at lower temperature in fuel rich regions. Compared to θ = 10 • , with which validation took place, NO x species remarkably increase up to 49%, 54%, and 57% at 15 • , 20 • , and 25 • , respectively.
The CO formation of methane/diesel RCCI engine is outlined in Figure 11. CO emissions decrease with increasing spray angle, between θ = 5 • and θ = 15 • . For θ > 15 • , CO slightly increases. A higher percentage of CO is observed for a narrower injector cone angle, θ = 5 • . Broader cone angles lead to better vaporization, therefore, better mixing and combustion. The obtained results are conformal with the pressure registered ( Figure 6). As shown in Figure 5, the temperature at θ = 5 • is the lowest, this indicates that combustion was incomplete.
The formation of CO 2 depends on the combustion behavior and the total injected mass of fuel [46]. The total diesel mass injected is 0.01649 g. The CO 2 specie shows an important mass at θ = 15 • . This indicates an excellent evaporation characteristic obtained. For the case of θ = 5 • , due to poor mixing and the presence in a rich region, CO could not be converted into CO 2 , that is why a lower CO 2 concentration is formed.
To evaluate engine performances in terms of output emissions and the overall performance for the different spray angle, combustion at θ = 15 • is best achieved leading with the lowest CO.

Conclusions
The present study explored the impact of different spray cone angles on the combustion and exergy efficiencies in the Internal Combustion Engine (ICE). The experimental data, from the literature, were performed on a modified single-cylinder caterpillar 3400 heavyduty diesel engine. The numerical work is conducted on a single cylinder four-stroke engine, using the dual fuel RCCI strategy. The KH-RT models the solid cone sprays and the diesel jet atomization. The assessed cone angles are 5 • , 10 • , 15 • , 20 • and 25 • . A one-component n-heptane (C 7 H 16 ) surrogate is used for diesel combustion. A detailed thermodynamic analysis balance, for Hydrocarbons (CH 4 /n-C 7 H 16 ) dual fuel RCCI engine is conducted. The rate of exergy work and the exergy efficiency data, for various spray cone angle, are reported at low load operating conditions. The exergy analysis, at ambient initial temperature condition, is performed on a reduced n-heptane (425 species and 3128 reactions) reaction mechanism. The obtained pressure and heat release rate, at θ = 10 • , validated the simulation with experimental data from the literature. The finding results of the numerical investigations can be summarized as follow: • The choice of methane as a low reactive fuel serves to extend the combustion phasing, which enhances engine performances leading, thus, to a reduction of CO specie concentration.

•
The spray cone angle has a significant impact on the heat transfer flux, pressure, temperature, and subsequent pollutants formation, as well.

•
The combustion efficiency is remarkably improved for θ = 20 • . This is registered by the temperature and pressure increase. At the mentioned spray angle, the pressure and the temperature equal 79.5 bar and 1766 K, respectively. They increase roughly by 11% and 6%, compared to the experimental work, carried out for θ = 10 • . Following the outlined results, θ = 20 • is assumed to be the best spray cone angle for combustion efficiency.
• In spite of the low emissions of CO 2 and NO x for θ = 5 • , the peak pressure at this angle is low (65 bar), this could be a limitation factor for the engine. The high CO emissions at Exhaust Valve Opening (EVO) indicate that the combustion was not complete.

•
The optimal spray angle ranges between 15 • and 20 • for the best combustion performances and reduced species (CO, CO2, and NOx). For θ < 15 • , the pressure decreases, and the combustion is not complete. For θ > 20 • , NO x emissions are remarkably larger than θ = 10 • . It increases by 58% compared to the experimental work.

•
The results of exergy analysis show a maximum value of exergy efficiency for 15 which represent an increase of roughly 5%, compared to experimental work, at θ = 10 • .