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°.

1. 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 Oxides (NOx) simultaneously. This option is commercially presented in three different approaches: 1/Homogeneous Charge Compression Ignition (HCCI) [5,6,7], which provides the ignition of a lean homogeneous air-fuel mixture; 2/Premixed Charge Compression Ignition (PCCI), is an adaptative form of HCCI, it controls combustion instability by a second injection to enrich wherever flame should be started; 3/Reactivity Controlled Compression Ignition (RCCI) [7,8], based on high reactivity fuel, that is, diesel and n-heptane as well as a low reactive fuel such as gasoline, iso-octane, and natural gas. The LTC is based on a combination of injection timing, mixture homogeneity, and dual fuel mode. It leads to a higher thermal efficiency, which allows the combustion of the lean mixture. A homogeneous LTC mixture results, also, at a lower and uniform flame temperature.

Within this work, the RCCI methane/diesel dual-fuel is used. Here, combustion deals with two main stages. First, the injection of low reactive fuel, which is air-methane based. Second, at a high pressure level of the compression stroke, diesel is injected and combustion takes place [9].

In fact, RCCI dual fuel engine experiences several problems. Looking behind at the conventional diesel engine, the RCCI technique suffers from unstable combustion performances, high fraction of unburned fuel, which delivers considerable CO emissions as well as low thermal efficiency [10].

Accordingly, valuable studies have been addressed to encourage RCCI engine efficiency. Bo yang et al. [11] evaluated the chronology timing of diesel and methane using Low Pressure Dual-fuel Direct Injection (LPDDI). They found a better compromise by retarding injection timing of CH4 at −112° CA ATDC and early diesel injection at −250° ATDC. The injection pressure and diesel–natural gas mixture fraction impacts on the combustion phasing (CP) were also reported experimentally by Poorghasemi et al. [12]. Combustion chamber geometry and bowl shape have the potential to alternate combustion and emissions [13]. A variety range of inlet valve closing temperature (TIVC) and exhaust gas recirculation at RCCI dual fuel have been reported [14]. Peng Jiang et al. [15] compared gasoline/hydrogenated-catalytic-biodiesel (HCB) RCCI with conventional gasoline/diesel port fuel injection. They tested the injection timings for each case. As result, direct injection effectively controlled mixture homogeneity and improved the combustion process. The injection delays have the effect of increasing CO and unburned hydrocarbon (UHC), although the delay of HCB direct injection improves combustion efficiency and lowers the output emissions. Zhu et al. [16] applied a direct injection of n-heptane combined with ethanol, gasoline, and butanol for every time. They found that injection timing advancement postpones the ignition delay, and the ethanol/n-heptane mixture has the capability for CO and soot reduction. Dempsey et al. [17] studied the effect of cetane number improvements on ethanol, gasoline, and methanol. They found that the mixture containing methanol and cetane number improvements could be similar to diesel. The effects of spray cone angle and swirl ratio [18], piston bowl and compresssion ratio [19,20], initial temperature [21], biodiesel-gasoline RCCI [22], have been also examined. Regarding the entropy production in related configurations, Ries et al. [23] studied the generated entropy in a turbulent impinging jet. Ganjehkaviri et al. [24] conducted a study to improve IC engine exergy using a heat recovery system. The thermodynamic cycle for the IC engine was also investigated [25]. The exergy losses [26] and the entropy generation [27], in a Detailed, Reduced, and Skeletal n-heptane combustion, are processed.

Previous studies point out a wide parametric variety that could affect RCCI performances. They have even recorded much progress. In spite of the mentioned advancement, very scarce research studies have addressed the effect of injection atomization and droplet pulverization on the emissions and engine efficiency of a light-duty methane/diesel dual fuel engine. Due to high output in Unburned Hydrocarbon (UHC), operating at low load mode remains challenging.

Therefore, an adapted RANS-based combustion model, as developed by “Kong-Reitz”, is used [28,29] in this paper to especially isolate the effect of spray cone angle on methane/diesel RCCI engine performances under low load operating conditions. The 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 CH4/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 (D32), 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.

2. Numerical Method

2.1. 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 air-methane 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 three-dimensional 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.

Figure 1

3D Computational model.

Table 1

Diesel engine specification.

Parameter	Value
Engine type	Caterpillar 3400
Bore × stroke [mm]	137.2 × 165.1
Connecting rod length [mm]	261.62
Displacement volume [L]	2.44
Compression ratio [-]	16.25
Nozzle type (hole × diameter) [mm]	6 × 0.23
Diesel fuel injection type	Direct injection
Natural gas injection type	Direct injection

2.2. 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 [30].

2.3. 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 well-suited for governing the turbulence in non-isotropic combustion. The mesh and timestep-independent 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.

Table 2

Boundary and initial conditions.

Boundary and Initial Condition	Value
Combustion chamber temperature at IVC [K]	360
Combustion chamber pressure at IVC [bar]	1.02
Kinetic energy [J]	10
Turbulence dissipation rate [m2/s3]	1732
Cylinder head wall temperature [K]	400
Piston temperature [K]	400
Liner wall temperature [K]	400
Turbulence model	k-ε RNG
Spray model	KH-RT
Injector type	Solid cone

Figure 2

Mesh representation: (a) Front view; (b) Top view.

2.4. 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, represents the mean velocity vector, and 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, , , 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. and 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 and , which represent the turbulent kinetic energy and the dissipation rate, respectively, are calculated as follows: In both equations , , , , and represent the model coefficients. They are summarized in Table 3. The quantities and are the renormalized and the turbulence dynamic viscosity, respectively. is the source term due to spray vaporization and is the Reynolds stress tensor. The terms and 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].

	C_ε1	C_s	C_ε2	C_ε3	1/Pr_k	1/Pr_ε
k-ε RNG	1.42	1.5	1.68	−0.9 to 1.726 [33]	1.39	1.39

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, is the velocity vector of the flow, and the quantity 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 . The quantities and are source terms due to spray vaporization and chemical reactions, respectively.

2.5. 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, which is formed during the primary breakup, are modeled with: where, equals 0.61 and represents the KH wavelength of the accurate growing wave, . The growth rate of the fastest wave and its related wavelengths are assumed as: Here, r is the radius of jet, is the surface tension, represents the density of the injected liquid, and T, Z, and We 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, τ 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 represents the deceleration owing to drag force. The quantity is the gas density. The term denotes the wave length of the accurate growing wave, while represents the growth rate during the secondary breakup. The formation of the child droplets, r, and the breakup time, , are assumed, respectively, as:

= 40 and  = 0.1 represent the size constant and the time constant of the RT breakup.

2.6. Exergy Balance Analysis

According to [35], the exergy balance is assessed following Equation (15): where is the inlet exergy rate, which depends on the air and fuel exergies rates as written in Equation (16), is the sum of exergy rates in form of work (), exhaust gases , and heat transferred , as provided in Equation (17),

Equations (15)–(17) can be put together and expressed as follows:In this case, the system is adiabatic , while is calculated through entropy generation according to the second law of thermodynamics.

, , and are assumed as the losses quantities and they are written in Equation (19) Based on Equations (18) and (19), the exergy efficiency can be deduced as: written using Equation (21) [36,37] yields:Note that the exergy of air can be calculated using Equation (22) [36]:Since the intake air is atmospheric, the exergy rate of air is neglected (). 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 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-C7H16 and CH4), the input fuel exergy rate is generally expressed as [38,39]:The parameters used for exergy efficiency () calculation are tabulated in Table 4 and Table 5.

Table 4

Specification of the used fuels.

	E˙x_fuel(kW)	LHV (MJ/kg)	m˙ (kg/h)	Inlet Temperature (K)
Diesel	5.7	44.643	0.4503	298.15
Methane	18.075	50	1.259

Table 5

Work exergy rate for various injector angle.

	θ = 5^°	θ = 10^°	θ = 15^°	θ = 20^°	θ = 25^°
E˙x_w(kW)	7.911	8.573	9.003	9.013	8.816

2.7. 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 Closing (IVC) to Exhaust Valve Opening (EVO). The start of injection (SOI) took place at 22° BTDC. The comparison between them is pointed out in Figure 3.

Figure 3

Numerical validation of pressure and heat release rate at DIT = 22° BTDC.

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.

3. 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.

3.1. 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.

Figure 4

View of heat transfer flux for various cone angles.

3.2. 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. 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].

Figure 5

Cylinder temperature variation over crank angle for different spray cone angles.

3.3. 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.

Figure 6

Cylinder pressure variation over crank angle for different spray cone 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.

3.4. 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.

Figure 7

HRR over crank angle for different cone angles.

3.5. 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.

Figure 8

Weber number for various spray cone−angles.

3.6. 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°.

Figure 9

D32 over crank angle for different spray cone−angles.

3.7. 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.

Figure 10

RCCI exergy efficiency for various spray angles.

3.8. Output Species

To control the produced species during the combustion process, the masses of CO2, CO, and EINOx are postprocessed at the exit of the combustion chamber. Figure 11 points out these species for the various spray cone angles.

Figure 11

Alteration of EINOx, CO, and CO2 emissions for different spray cone−angles.

Figure 11 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 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 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.

4. 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 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 (C7H16) surrogate is used for diesel combustion. A detailed thermodynamic analysis balance, for Hydrocarbons (CH4/n-C7H16) 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 CO2 and NOx 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°, NOx 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 , which represent an increase of roughly 5%, compared to experimental work, at θ = 10°.