# Investigation of the Effects of Steam Injection on Equilibrium Products and Thermodynamic Properties of Diesel and Biodiesel Fuels

Authors: Jean Paul Gram Shou; Marcel Obounou; Timoléon Crépin Kofané; Mahamat Hassane Babikir
Journal: Journal of Combustion
 (2020)

---

## Abstract

The effects of steam injection on combustion products and thermodynamic properties of diesel fuel, soybean oil-based biodiesel (NBD), and waste cooking oil biodiesel (WCOB) are examined in this study by considering the chemical equilibrium. The model gives equilibrium mole fractions, specific heat of the exhaust mixtures of 10 combustion products, and adiabatic flame temperatures. The results show that the mole fractions of carbon monoxide (CO) and carbon dioxide (CO2) decrease with the steam injection ratios. Nitric oxide (NO) mole fractions decrease with the steam injections ratios for lean mixtures. The specific heat of combustion products increases with the steam injection ratios. The equilibrium combustion products obtained can be used to calculate the nonequilibrium values of NO in the exhaust gases using some existing correlations of NO kinetics.


---

## Body

## 1. Introduction

Fossil-based fuels derived from crude oil are the main sources of energy for internal combustion engines worldwide in general and in Cameroon in particular. However, combustion in internal combustion engines is a major source of pollutant emissions. More severe environmental restrictions have been put into effect in many countries worldwide due to the growing concern about environmental issues. In addition, the use of fossil-based fuels is still dominant in almost all internal combustion engines despite the fact that the resources of the fossil-based fuels are depleting rapidly. Thus, researchers carry out studies on combustion in internal combustion engines and alternative fuels. In this regard, biodiesel could be considered as an alternative fuel to fossil diesel to meet the growing fuel demand and the severe environmental norms. Biodiesel is produced from wider feedstock by transesterification, microemulsion, dilution, and pyrolysis. Rapeseed, soybean, and palm oils have been largely preferred as edible raw materials for the production of biodiesel worldwide [1]. Biodiesel has gained prominence as an attractive fuel because it is biodegradable, nontoxic, environmentally friendly, and free from sulfur and aromatic contents [2]. The diesel engine operating with biodiesel can decrease carbon monoxide (CO), hydrocarbons (HC), and smoke emissions while increasing nitrogen oxides (NOx) [3–5]. The biodiesel properties are very close to the ones of diesel. However, the shortcomings of biodiesel are higher viscosity, density, cetane number, flash point, poor low-temperature properties and lower volatility, energy content, and oxidative and storage stability in comparison to the diesel fuel [6]. Higher cloud point and pour point, higher nitrogen oxide emissions, injector coking, engine compatibility, and higher engine wear are also the shortcomings of biodiesel which affect the combustion characteristics [7]. The poor low-temperature properties of biodiesel are an obstacle to its utilization in cold weather conditions [8]. The differences in chemical properties of fuels affect the combustion process development, which lead to the alteration of performance and emissions of engine [9]. Biodiesels are formed from many different classes of saturated and unsaturated fatty acids, which form their physical and chemical properties [10]. Higher viscosity results in higher spray penetration into the combustion chamber, which is very important to enhance air-fuel mixing quality especially for big and slow engines with large open combustion chambers, but there are more adverse impacts for small and fast engines such as poor fuel atomization, larger droplet size, bad mixing quality, incomplete combustion, and lower combustion efficiency, which contribute to decreasing the engine performance and increasing the exhaust emissions [11]. Biodiesel combustion is difficult to model because of the diversity of its sources and complexity in molecular structure [12–16]. Chemical equilibrium model is a better estimation for a combustion model assuming that the combustion products are at a chemical equilibrium in a high-temperature combustion. The combustion products react together to produce and remove each species at equal rates, and consequently, the net change in species composition remains constant at a given condition [17]. Chemical equilibrium models are used in thermodynamic simulations of internal combustion engines [18]. A simple model including only six species of the products, namely, carbon dioxide (CO2), water (H2O), nitrogen (N2), oxide (O2), carbon monoxide (CO), and hydrogen (H2), has been defined by Heywood [19]. Gonca [20] investigated the effects of steam injection on performance and nitric oxide (NO) emissions of a diesel engine running with ethanol-diesel blend. In his study, various ethanol-diesel blends were evaluated depending on the steam/air ratios by means of constituted solving schema based on chemical equilibrium. Ngayihi Abbe et al. [21] carried out a numerical study to compare four biodiesel surrogates using a zero-dimensional thermodynamic model involving a chemical equilibrium combustion model. A chemical equilibrium model was used by Ust and Kayadelen [22]. In their study, 10 combustion products were used to predict the equilibrium and thermodynamic properties in an H2O injected combustion system at various H2O rates. They neglected the CO and H2 formations in lean mixtures. They formed their equation systems for lean and rich mixtures separately. Rakopoulos et al. [23] performed a two-zone thermodynamic model for combustion and emissions formation in a direct injection diesel engine. They used a chemical equilibrium scheme involving 11 species of combustion products for calculating the constituents in exhaust gases. Diotallevi [24] developed a multi-zone thermodynamic model of a diesel engine for NOx formation. He created a chemical equilibrium scheme by considering 10 combustion products. The nonlinear system of equations created was solved using an iterative method developed by the author using the Matlab program. Mourya and Roy [25] performed a study on the combustion modelling of a diesel engine operating using jatropha biodiesel and diesel engine blends according to a chemical equilibrium model. The C program and Mathematica software developed by the authors were used to solve their equations. Kayadelen [26] developed a multi-featured model for estimation of the mole fractions of 14 equilibrium combustion products, thermodynamic properties and constant-pressure adiabatic flame temperature of fuels, surrogates, and fuel additives. He constituted his equation systems for lean and rich mixtures separately. CO and H2 mole fractions were neglected in lean mixtures. The nonlinear equation systems were solved using both Newton-Raphson and Gauss-Seidel methods. Some equilibrium computer programs exist like GASEQ [27] and Chemical Equilibrium with Applications (CEA) of NASA [28]. However, GASEQ software neglects the effects of dissociations on the specific heat of the combustion products. Consequently, noticeable errors may result in engine performance calculations. On the other hand, CEA of NASA necessitates adiabatic flame temperature or enthalpy of combustion to be given. Important differences are reported between the specific heat results of GASEQ and CEA of NASA particularly at high temperatures. In literature, there is a lack of studies carried out on the effects of steam injection on equilibrium products and thermodynamic properties of diesel and biodiesel fuels. Nevertheless, the equilibrium combustion products can be used to calculate the nonequilibrium NO emissions in exhaust gases. The present model takes into account the effect of combustion products dissociations with temperature on specific heat of the gas mixture. The adiabatic flame temperature or combustion enthalpy is not needed as an input, but they can be calculated by the model. The main purpose of the present study is to investigate the effects of steam injection on equilibrium products and thermodynamic properties of diesel fuel, soybean oil-based biodiesel (NBD), and waste cooking oil biodiesel (WCOB). For this purpose, the solving schema of the combustion model based on a chemical equilibrium was developed for diesel and biodiesel fuels. The nonlinear system of equations is solved using Newton-Raphson and Gauss-Seidel methods. The diesel and biodiesel with chemical formulae C14.09H24.78 [29], C18.74H34.43O2 [29], and C18H36O2 [30] were used in this study to represent diesel fuel, soybean oil-based biodiesel (NBD), and waste cooking oil biodiesel (WCOBD), respectively.


## 2. Theoretical Equilibrium Combustion Products Model

### 2.1. Equilibrium Combustion Products

The combustion products are supposed to consist of ten species which are all assumed as ideal gases and they are defined by dissociation considerations as follows [31]:(1)12H2⇌H,K1=y7p1/2y61/2,(2)12O2⇌O,K2=y8p1/2y41/2,(3)12H2+12O2⇌OH,K3=y9y41/2y61/2,(4)12O2+12N2⇌NO,K4=y10y31/2y41/2,(5)H2+12O2⇌H2O,K5=y2y61/2y6p1/2,(6)CO+12O2⇌CO2,K6=y1y41/2y5p1/2,where the unit of pressure p is atmospheres and K1 to K6 are the equilibrium constants of the reactions. Olikara and Borman [32] have curve-fitted the equilibrium constants Ki to tabulated JANAF thermochemical tables for the temperatures between 600 K and 4000 K and its values are calculated by using(7)logKi=AilnT1000+BiT+Ci+DiT+EiT2,where the equilibrium constant curve-fit coefficients Ai, Bi, Ci, Di, and Ei are listed in Table 1 and presented by Ferguson [31].Table 1
Olikara and Borman equilibrium constant curve-fit coefficients.
iAiBiCiDiEi10.432168−0.112464 × 1050.267269 × 101−0.745744 × 10−40.242484 × 10−820.310805−0.129540 × 1050.321779 × 101−0.738336 × 10−40.344645 × 10−83−0.141784−0.213308 × 1040.8534610.355015 × 10−4−0.310227 × 10−840.150879 × 10−1−0.470959 × 1040.6460960.272805 × 10−5−0.154444 × 10−85−0.7523640.124210 × 105−0.260286 × 1010.259556 × 10−3−0.162687 × 10−76−0.415302 × 10−20.148627 × 105−0.475746 × 1010.124699 × 10−3−0.900227 × 10−8


### 2.2. Equilibrium Combustion Model

#### 2.2.1. Low-Temperature Products

At low temperatures (T < 1000 K) and for α/γ ratios less than one, the overall chemical reaction, which describes the combustion (ϕ < 3) is given, as follows:(8)εϕCαHβOγNδ+0.21O2+0.79N2⟶ν1CO2+ν2H2O+ν3N2+ν4O2+ν5CO+ν6H2.For the lean and stoichiometric combustion, it is assumed that there would be enough oxygen to oxidize all CO and H2 and that means CO and H2 are negligible. It can be noticed that the equilibrium composition is independent from temperature and pressure. It only depends on equivalence ratio ϕ. The product composition is determined using atom balance equations. Products of wet combustion for the low temperature are given in Table 2. For the rich case, it is assumed that there is insufficient oxygen to oxidize all CO and H2, which means all the oxygen is consumed.Table 2
Products of wet combustion forT < 1000 K.
iSpeciesϕ ≤ 1ϕ > 11CO2αϕεαϕε − ν52H2Oβϕε + xϕε (γ − 2α) + 0.42 + ν53N20.79 +δϕε/20.79 +δφε/24O20.21 (1 −ϕ)05CO0Ν56H200.42 (ϕ − 1) − ν5When writing atom balances, five unknowns and four equations are found. One more equation (10) is used and it depends on the equilibrium constant for the water-gas reaction given by Ferguson [31] for the range of temperatures between 400 K and 3200 K:(9)lnK=2.743−1.761t−1.611t2+0.2803t3,where t = T/1000.The equilibrium constant at constant pressure K of the water-gas reaction is expressed as follows:(10)CO2+H2⇌CO+H2O,K=ν2ν5ν1ν6.


#### 2.2.2. High-Temperature Products

At higher temperatures (T ≥ 1000 K) and for α/γ ratios less than one, the overall chemical reaction, which describes the combustion (ϕ < 3), is given as follows:(11)εϕCαHβOγNδ+0.21O2+0.79N2+xH2O⟶ν1CO2+ν2H2O+ν3N2+ν4O2+ν5CO+ν6H2+ν7H+ν8O+ν9OH+ν10Nwhere νi denotes the number of moles for each product, and α, β, γ, and δ represent the atom numbers of carbon, hydrogen, oxygen, and nitrogen in the used fuel, respectively. For diesel fuel, γ and δ are zeroes. When the product mole fraction composition is known at a given temperature, equivalence ratio, and pressure, the thermodynamic properties of interest such as enthalpy, entropy, specific volume, and internal energy can be computed.The atom balance for various elements in equation (11), the constraint that the mole fraction of all the products adds up to one, the unknown total product moles N, and the six dissociation equations (1)–(6) provided by the criteria of equilibrium among combustion products yield 11 equations for 11 unknowns, namely, unknown mole fractions yi and unknown total product moles N. After writing the mole fractions of the other combustion species with respect to the mole fractions of four independent variables y3, y4, y5, and y6, the following nonlinear equations are obtained:(12)f1=C6y41/2y5+C5y41/2y6+y3+y4+y5+y6+C1y61/2+C2y41/2+C3y41/2y61/2+C4y31/2y41/2−1=0,(13)f2=2C5y41/2y6+2y6+C1y61/2+C3y41/2y61/2−d1C6y41/2y5+y5=0,(14)f3=2C6y41/2y5+C5y41/2y6+2y4+y5+C2y41/2+C3y41/2y61/2+C4y31/2y41/2−d2C6y41/2y5+y5=0,(15)f4=2y3+C4y31/2y41/2−d3C6y41/2y5+y5=0,(16)C1=K1p1/2,(17)C2=K2p1/2,(18)C3=K3.(19)C4=K4,(20)C5=K5p1/2,(21)C6=K6p1/2,(22)d1=εϕβ+2xεϕα,(23)d2=εϕγ+0.42+xεϕα,(24)d3=εϕδ+1.58εϕα.The equivalence ratioϕ and the molar fuel/air ratio ε are expressed as follows:(25)ϕ=mf/mamf/mas,(26)ε=0.21α−γ/2+β/4.By lettingY = (y3, y4, y5, y6)T be the vector of the independent variables, δY = (δy3, δy4, δy5, δy6)T the solution vector of the linear system, and F the vector containing the functions f3, f4, f5 and f6, the system of nonlinear equations can be written in a single expression using vectors:(27)FY=f3Y,f4Y,f5Y,f6Y=0.The exact solution of the problem isY=Y0+δY0, which means(28)FY=FY0+δY0=0,where Y0 is a given rough initial guess for the solution. In the neighbourhoods of Y0, the functions fi (Y) can be expanded in Taylor series and truncated after the first derivative. We obtain a set of linear equations for the vector δY0 that move each value fi (Y0) closer to zero simultaneously:(29)ℑY0δY0=−FY0.The Jacobian matrixℑ is expressed as(30)ℑY0=∂f1Y0∂y3∂f1Y0∂y4⋯∂f1Y0∂y6∂f2Y0∂y3∂f2Y0∂y4⋯∂f2Y0∂y6⋮⋮⋱⋮∂f4Y0∂y3∂f4Y0∂y4⋯∂f4Y0∂y6.The set of linear equations equation (29) is solved for δY(0) using Gauss-Seidel iterative method. The new approximated solution is calculated as Y1=Y0+δY0.IfY1 does not approximate the solution to the given tolerance, successive approximations of the solution are obtained as follows:(31)ℑYkδYk=−FYk.The procedure is repeated untilδYk reaches a specified tolerance which leads to the values for Y:(32)Y=Yk+δYk.The other dependent unknowns are expressed after determining the mole fractionsy3, y4, y5, and y6 as follows:(33)y1=C6y41/2y5,(34)y2=C5y41/2y6,(35)y7=C1y61/2,(36)y8=C2y41/2,(37)y9=C3y41/2y61/2,(38)y10=C4y31/2y41/2.


### 2.3. Thermodynamic Properties

At constant pressure, temperature variation has an effect on specific heat due to the dissociations of species at a high temperature. The effect of temperature on mole fractions should be considered during the equilibrium specific heat calculation differentiating with respect to temperature written as(39)ℑYkYTk=−FTYk,(40)FT=∂F∂T,YT=∂Y∂T.The same procedure is used to solve equation (39). The results are used in computing the specific heat of the gas mixture in equation (46).Gordon and McBride [33] proposed the expressions of molar specific heat at constant pressure c¯p,i enthalpy h¯i and entropy s¯i0 values of each species that were curve-fitted to the tabulated JANAF thermochemical tables [34].At constant pressure, enthalpy of the gas mixtureh¯ changes due to the dissociations because the mole fractions of the species change with temperature:(41)h¯=∑i=110yih¯i,(42)s¯=∑i=110Ruyis¯i0Ru−lnyipp0,(43)M=∑i=110yiMi,(44)h=1M∑i=110yih¯i,(45)s=RuM∑i=110yis¯i0Ru−lnyipp0,(46)cp=1M∑i=110yic¯p,i+h¯i∂yi∂T−MTMyih¯i,where MT=∂M/∂T=∑i=110Mi∂yi/∂T and cp=∂h/∂Tp.


### 2.4. Constant-Pressure Adiabatic Flame Temperature

The constant-pressure adiabatic flame temperature (Tad) is calculated iteratively using (47). An initial temperature is required to determine the total enthalpy of reactants hu:(47)Tadi+1=Tadi−hu−hcp,where h is the enthalpy of combustion products obtained by (44) and cp is the specific heat of the combustion products obtained by (46).For each iterationi, h and cp are reevaluated until an acceptable prescribed tolerance is reached.


### 2.5. Fuel Properties and Standard Specifications

The physical and chemical properties of the fuels are presented in Table3.Table 3
Chemical and physical properties of test fuels.
Fuel propertiesDiesel [30]NBD [30]WCOBDensity at 15°C (kg·m−3)853.7881.4865.19 [36]Kinematic viscosity at 40°C (mm2·s−1)2.82714.26913.5 [36]LHV (kJ·kg−1)426403738837000 [31]Oxygen (wt. %)0.0010.9811.27Cetane number42.651.554 [31]Chemical formulaC14.09H24.78C18.74H34.43O2C18H36O2 [31]Standard specifications for biodiesel and diesel fuels are shown in Table4 according to ASTM D6751 and EN 14214.Table 4
Standard specifications for biodiesel.
Fuel propertiesASTM D6751EN 14214 [36]Density at 15°C (kg·m−3)880 [37]860–900Kinematic viscosity at 40°C (mm2·s−1)1.9–6.0 [37]3.5–5.0Cetane numberMin. 47 [37]Min. 51Flash point (°C)Min. 93 [37]Min. 101Sulfur % (mg/kg)Max. 0.05 [30]Max. 10Oxidation stability at 110°C (hr)Min. 3 [37]Min. 6


### 2.6. Validation of the Model

The validation of the model is done by comparing the model results with the ones of software CEA of NASA [28] and software GASEQ [27], which use element potential method and minimization of Gibbs free energy approach, respectively. The simulations of methane with steam injection are carried out. The mole fraction of chemical equilibrium products and thermodynamic properties obtained with the model compared with the ones of software CEA of NASA and GASEQ at ϕ = 0.6 with steam injection 10% are mentioned in Table 5 and at ϕ = 1.2 with steam injection 10% are listed in Table 6.Table 5
Comparison of the present model results with the ones obtained from CEA of NASA and GASEQ software for equilibrium combustion products of methane atϕ = 0.6 and steam injection 10%.
ϕ = 0.6Model resultsGASEQ resultsCEA resultsGASEQ relative error (%)CEA relative error (%)CO20.0515060.051540.051510.0660.066H2O0.2339200.233990.233850.0290.029N20.6455590.645530.64553−0.005−0.005O20.0683530.068270.06827−0.122−0.122CO1.328 × 10−71.33 × 10−73.27 × 10−7−0.157−0.157H22.419 × 10−72.41 × 10−72.41 × 10−7−0.360−0.360H1.269 × 10−91.27 × 10−91.27 × 10−90.030.03O1.506 × 10−71.50 × 10−71.50 × 10−7−0.429−0.429OH3.263 × 10−53.39 × 10−53.39 × 10−53.7513.751NO6.278 × 10−46.27 × 10−46.27 × 10−4−0.126−0.126Cp (kJ·kg−1·K−1)1.464881.453831.45383−0.76−0.76Tin (K)3003003000.0000.000Tad (K)1482.31482.31482.30.0000.000Table 6
Comparison of the present model results with the ones obtained from CEA of NASA and GASEQ software for equilibrium combustion products of methane atϕ = 1.2 and steam injection 10%.
ϕ = 1.2Model resultsGASEQ resultsCEA resultsGASEQ relative error (%)CEA relative error (%)CO20.0637290.063680.06368−0.077−0.077H2O0.2781730.278220.278220.0170.017N20.5947990.594800.594800.000170.00017O25.826 × 10−85.78 × 10−85.78 × 10−8−0.790−0.790CO0.0311380.031190.031190.1680.168H20.0321060.032060.03206−0.144−0.144H2.988 × 10−52.98 × 10−52.98 × 10−5−0.260−0.260O1.535 × 10−81.53 × 10−81.53 × 10−8−0.345−0.345OH2.199 × 10−52.95 × 10−62.68 × 10−62.6792.57NO2.962 × 10−62.95 × 10−60.646096−0.394−0.394Cp (kJ·kg−1·K−1)1.645921.630421.63042−0.95−0.95Tin (K)3003003000.0000.000Tad (K)1919.31919.31919.30.0000.000


## 2.1. Equilibrium Combustion Products

The combustion products are supposed to consist of ten species which are all assumed as ideal gases and they are defined by dissociation considerations as follows [31]:(1)12H2⇌H,K1=y7p1/2y61/2,(2)12O2⇌O,K2=y8p1/2y41/2,(3)12H2+12O2⇌OH,K3=y9y41/2y61/2,(4)12O2+12N2⇌NO,K4=y10y31/2y41/2,(5)H2+12O2⇌H2O,K5=y2y61/2y6p1/2,(6)CO+12O2⇌CO2,K6=y1y41/2y5p1/2,where the unit of pressure p is atmospheres and K1 to K6 are the equilibrium constants of the reactions. Olikara and Borman [32] have curve-fitted the equilibrium constants Ki to tabulated JANAF thermochemical tables for the temperatures between 600 K and 4000 K and its values are calculated by using(7)logKi=AilnT1000+BiT+Ci+DiT+EiT2,where the equilibrium constant curve-fit coefficients Ai, Bi, Ci, Di, and Ei are listed in Table 1 and presented by Ferguson [31].Table 1
Olikara and Borman equilibrium constant curve-fit coefficients.
iAiBiCiDiEi10.432168−0.112464 × 1050.267269 × 101−0.745744 × 10−40.242484 × 10−820.310805−0.129540 × 1050.321779 × 101−0.738336 × 10−40.344645 × 10−83−0.141784−0.213308 × 1040.8534610.355015 × 10−4−0.310227 × 10−840.150879 × 10−1−0.470959 × 1040.6460960.272805 × 10−5−0.154444 × 10−85−0.7523640.124210 × 105−0.260286 × 1010.259556 × 10−3−0.162687 × 10−76−0.415302 × 10−20.148627 × 105−0.475746 × 1010.124699 × 10−3−0.900227 × 10−8


## 2.2. Equilibrium Combustion Model

### 2.2.1. Low-Temperature Products

At low temperatures (T < 1000 K) and for α/γ ratios less than one, the overall chemical reaction, which describes the combustion (ϕ < 3) is given, as follows:(8)εϕCαHβOγNδ+0.21O2+0.79N2⟶ν1CO2+ν2H2O+ν3N2+ν4O2+ν5CO+ν6H2.For the lean and stoichiometric combustion, it is assumed that there would be enough oxygen to oxidize all CO and H2 and that means CO and H2 are negligible. It can be noticed that the equilibrium composition is independent from temperature and pressure. It only depends on equivalence ratio ϕ. The product composition is determined using atom balance equations. Products of wet combustion for the low temperature are given in Table 2. For the rich case, it is assumed that there is insufficient oxygen to oxidize all CO and H2, which means all the oxygen is consumed.Table 2
Products of wet combustion forT < 1000 K.
iSpeciesϕ ≤ 1ϕ > 11CO2αϕεαϕε − ν52H2Oβϕε + xϕε (γ − 2α) + 0.42 + ν53N20.79 +δϕε/20.79 +δφε/24O20.21 (1 −ϕ)05CO0Ν56H200.42 (ϕ − 1) − ν5When writing atom balances, five unknowns and four equations are found. One more equation (10) is used and it depends on the equilibrium constant for the water-gas reaction given by Ferguson [31] for the range of temperatures between 400 K and 3200 K:(9)lnK=2.743−1.761t−1.611t2+0.2803t3,where t = T/1000.The equilibrium constant at constant pressure K of the water-gas reaction is expressed as follows:(10)CO2+H2⇌CO+H2O,K=ν2ν5ν1ν6.


### 2.2.2. High-Temperature Products

At higher temperatures (T ≥ 1000 K) and for α/γ ratios less than one, the overall chemical reaction, which describes the combustion (ϕ < 3), is given as follows:(11)εϕCαHβOγNδ+0.21O2+0.79N2+xH2O⟶ν1CO2+ν2H2O+ν3N2+ν4O2+ν5CO+ν6H2+ν7H+ν8O+ν9OH+ν10Nwhere νi denotes the number of moles for each product, and α, β, γ, and δ represent the atom numbers of carbon, hydrogen, oxygen, and nitrogen in the used fuel, respectively. For diesel fuel, γ and δ are zeroes. When the product mole fraction composition is known at a given temperature, equivalence ratio, and pressure, the thermodynamic properties of interest such as enthalpy, entropy, specific volume, and internal energy can be computed.The atom balance for various elements in equation (11), the constraint that the mole fraction of all the products adds up to one, the unknown total product moles N, and the six dissociation equations (1)–(6) provided by the criteria of equilibrium among combustion products yield 11 equations for 11 unknowns, namely, unknown mole fractions yi and unknown total product moles N. After writing the mole fractions of the other combustion species with respect to the mole fractions of four independent variables y3, y4, y5, and y6, the following nonlinear equations are obtained:(12)f1=C6y41/2y5+C5y41/2y6+y3+y4+y5+y6+C1y61/2+C2y41/2+C3y41/2y61/2+C4y31/2y41/2−1=0,(13)f2=2C5y41/2y6+2y6+C1y61/2+C3y41/2y61/2−d1C6y41/2y5+y5=0,(14)f3=2C6y41/2y5+C5y41/2y6+2y4+y5+C2y41/2+C3y41/2y61/2+C4y31/2y41/2−d2C6y41/2y5+y5=0,(15)f4=2y3+C4y31/2y41/2−d3C6y41/2y5+y5=0,(16)C1=K1p1/2,(17)C2=K2p1/2,(18)C3=K3.(19)C4=K4,(20)C5=K5p1/2,(21)C6=K6p1/2,(22)d1=εϕβ+2xεϕα,(23)d2=εϕγ+0.42+xεϕα,(24)d3=εϕδ+1.58εϕα.The equivalence ratioϕ and the molar fuel/air ratio ε are expressed as follows:(25)ϕ=mf/mamf/mas,(26)ε=0.21α−γ/2+β/4.By lettingY = (y3, y4, y5, y6)T be the vector of the independent variables, δY = (δy3, δy4, δy5, δy6)T the solution vector of the linear system, and F the vector containing the functions f3, f4, f5 and f6, the system of nonlinear equations can be written in a single expression using vectors:(27)FY=f3Y,f4Y,f5Y,f6Y=0.The exact solution of the problem isY=Y0+δY0, which means(28)FY=FY0+δY0=0,where Y0 is a given rough initial guess for the solution. In the neighbourhoods of Y0, the functions fi (Y) can be expanded in Taylor series and truncated after the first derivative. We obtain a set of linear equations for the vector δY0 that move each value fi (Y0) closer to zero simultaneously:(29)ℑY0δY0=−FY0.The Jacobian matrixℑ is expressed as(30)ℑY0=∂f1Y0∂y3∂f1Y0∂y4⋯∂f1Y0∂y6∂f2Y0∂y3∂f2Y0∂y4⋯∂f2Y0∂y6⋮⋮⋱⋮∂f4Y0∂y3∂f4Y0∂y4⋯∂f4Y0∂y6.The set of linear equations equation (29) is solved for δY(0) using Gauss-Seidel iterative method. The new approximated solution is calculated as Y1=Y0+δY0.IfY1 does not approximate the solution to the given tolerance, successive approximations of the solution are obtained as follows:(31)ℑYkδYk=−FYk.The procedure is repeated untilδYk reaches a specified tolerance which leads to the values for Y:(32)Y=Yk+δYk.The other dependent unknowns are expressed after determining the mole fractionsy3, y4, y5, and y6 as follows:(33)y1=C6y41/2y5,(34)y2=C5y41/2y6,(35)y7=C1y61/2,(36)y8=C2y41/2,(37)y9=C3y41/2y61/2,(38)y10=C4y31/2y41/2.


## 2.2.1. Low-Temperature Products

At low temperatures (T < 1000 K) and for α/γ ratios less than one, the overall chemical reaction, which describes the combustion (ϕ < 3) is given, as follows:(8)εϕCαHβOγNδ+0.21O2+0.79N2⟶ν1CO2+ν2H2O+ν3N2+ν4O2+ν5CO+ν6H2.For the lean and stoichiometric combustion, it is assumed that there would be enough oxygen to oxidize all CO and H2 and that means CO and H2 are negligible. It can be noticed that the equilibrium composition is independent from temperature and pressure. It only depends on equivalence ratio ϕ. The product composition is determined using atom balance equations. Products of wet combustion for the low temperature are given in Table 2. For the rich case, it is assumed that there is insufficient oxygen to oxidize all CO and H2, which means all the oxygen is consumed.Table 2
Products of wet combustion forT < 1000 K.
iSpeciesϕ ≤ 1ϕ > 11CO2αϕεαϕε − ν52H2Oβϕε + xϕε (γ − 2α) + 0.42 + ν53N20.79 +δϕε/20.79 +δφε/24O20.21 (1 −ϕ)05CO0Ν56H200.42 (ϕ − 1) − ν5When writing atom balances, five unknowns and four equations are found. One more equation (10) is used and it depends on the equilibrium constant for the water-gas reaction given by Ferguson [31] for the range of temperatures between 400 K and 3200 K:(9)lnK=2.743−1.761t−1.611t2+0.2803t3,where t = T/1000.The equilibrium constant at constant pressure K of the water-gas reaction is expressed as follows:(10)CO2+H2⇌CO+H2O,K=ν2ν5ν1ν6.


## 2.2.2. High-Temperature Products

At higher temperatures (T ≥ 1000 K) and for α/γ ratios less than one, the overall chemical reaction, which describes the combustion (ϕ < 3), is given as follows:(11)εϕCαHβOγNδ+0.21O2+0.79N2+xH2O⟶ν1CO2+ν2H2O+ν3N2+ν4O2+ν5CO+ν6H2+ν7H+ν8O+ν9OH+ν10Nwhere νi denotes the number of moles for each product, and α, β, γ, and δ represent the atom numbers of carbon, hydrogen, oxygen, and nitrogen in the used fuel, respectively. For diesel fuel, γ and δ are zeroes. When the product mole fraction composition is known at a given temperature, equivalence ratio, and pressure, the thermodynamic properties of interest such as enthalpy, entropy, specific volume, and internal energy can be computed.The atom balance for various elements in equation (11), the constraint that the mole fraction of all the products adds up to one, the unknown total product moles N, and the six dissociation equations (1)–(6) provided by the criteria of equilibrium among combustion products yield 11 equations for 11 unknowns, namely, unknown mole fractions yi and unknown total product moles N. After writing the mole fractions of the other combustion species with respect to the mole fractions of four independent variables y3, y4, y5, and y6, the following nonlinear equations are obtained:(12)f1=C6y41/2y5+C5y41/2y6+y3+y4+y5+y6+C1y61/2+C2y41/2+C3y41/2y61/2+C4y31/2y41/2−1=0,(13)f2=2C5y41/2y6+2y6+C1y61/2+C3y41/2y61/2−d1C6y41/2y5+y5=0,(14)f3=2C6y41/2y5+C5y41/2y6+2y4+y5+C2y41/2+C3y41/2y61/2+C4y31/2y41/2−d2C6y41/2y5+y5=0,(15)f4=2y3+C4y31/2y41/2−d3C6y41/2y5+y5=0,(16)C1=K1p1/2,(17)C2=K2p1/2,(18)C3=K3.(19)C4=K4,(20)C5=K5p1/2,(21)C6=K6p1/2,(22)d1=εϕβ+2xεϕα,(23)d2=εϕγ+0.42+xεϕα,(24)d3=εϕδ+1.58εϕα.The equivalence ratioϕ and the molar fuel/air ratio ε are expressed as follows:(25)ϕ=mf/mamf/mas,(26)ε=0.21α−γ/2+β/4.By lettingY = (y3, y4, y5, y6)T be the vector of the independent variables, δY = (δy3, δy4, δy5, δy6)T the solution vector of the linear system, and F the vector containing the functions f3, f4, f5 and f6, the system of nonlinear equations can be written in a single expression using vectors:(27)FY=f3Y,f4Y,f5Y,f6Y=0.The exact solution of the problem isY=Y0+δY0, which means(28)FY=FY0+δY0=0,where Y0 is a given rough initial guess for the solution. In the neighbourhoods of Y0, the functions fi (Y) can be expanded in Taylor series and truncated after the first derivative. We obtain a set of linear equations for the vector δY0 that move each value fi (Y0) closer to zero simultaneously:(29)ℑY0δY0=−FY0.The Jacobian matrixℑ is expressed as(30)ℑY0=∂f1Y0∂y3∂f1Y0∂y4⋯∂f1Y0∂y6∂f2Y0∂y3∂f2Y0∂y4⋯∂f2Y0∂y6⋮⋮⋱⋮∂f4Y0∂y3∂f4Y0∂y4⋯∂f4Y0∂y6.The set of linear equations equation (29) is solved for δY(0) using Gauss-Seidel iterative method. The new approximated solution is calculated as Y1=Y0+δY0.IfY1 does not approximate the solution to the given tolerance, successive approximations of the solution are obtained as follows:(31)ℑYkδYk=−FYk.The procedure is repeated untilδYk reaches a specified tolerance which leads to the values for Y:(32)Y=Yk+δYk.The other dependent unknowns are expressed after determining the mole fractionsy3, y4, y5, and y6 as follows:(33)y1=C6y41/2y5,(34)y2=C5y41/2y6,(35)y7=C1y61/2,(36)y8=C2y41/2,(37)y9=C3y41/2y61/2,(38)y10=C4y31/2y41/2.


## 2.3. Thermodynamic Properties

At constant pressure, temperature variation has an effect on specific heat due to the dissociations of species at a high temperature. The effect of temperature on mole fractions should be considered during the equilibrium specific heat calculation differentiating with respect to temperature written as(39)ℑYkYTk=−FTYk,(40)FT=∂F∂T,YT=∂Y∂T.The same procedure is used to solve equation (39). The results are used in computing the specific heat of the gas mixture in equation (46).Gordon and McBride [33] proposed the expressions of molar specific heat at constant pressure c¯p,i enthalpy h¯i and entropy s¯i0 values of each species that were curve-fitted to the tabulated JANAF thermochemical tables [34].At constant pressure, enthalpy of the gas mixtureh¯ changes due to the dissociations because the mole fractions of the species change with temperature:(41)h¯=∑i=110yih¯i,(42)s¯=∑i=110Ruyis¯i0Ru−lnyipp0,(43)M=∑i=110yiMi,(44)h=1M∑i=110yih¯i,(45)s=RuM∑i=110yis¯i0Ru−lnyipp0,(46)cp=1M∑i=110yic¯p,i+h¯i∂yi∂T−MTMyih¯i,where MT=∂M/∂T=∑i=110Mi∂yi/∂T and cp=∂h/∂Tp.


## 2.4. Constant-Pressure Adiabatic Flame Temperature

The constant-pressure adiabatic flame temperature (Tad) is calculated iteratively using (47). An initial temperature is required to determine the total enthalpy of reactants hu:(47)Tadi+1=Tadi−hu−hcp,where h is the enthalpy of combustion products obtained by (44) and cp is the specific heat of the combustion products obtained by (46).For each iterationi, h and cp are reevaluated until an acceptable prescribed tolerance is reached.


## 2.5. Fuel Properties and Standard Specifications

The physical and chemical properties of the fuels are presented in Table3.Table 3
Chemical and physical properties of test fuels.
Fuel propertiesDiesel [30]NBD [30]WCOBDensity at 15°C (kg·m−3)853.7881.4865.19 [36]Kinematic viscosity at 40°C (mm2·s−1)2.82714.26913.5 [36]LHV (kJ·kg−1)426403738837000 [31]Oxygen (wt. %)0.0010.9811.27Cetane number42.651.554 [31]Chemical formulaC14.09H24.78C18.74H34.43O2C18H36O2 [31]Standard specifications for biodiesel and diesel fuels are shown in Table4 according to ASTM D6751 and EN 14214.Table 4
Standard specifications for biodiesel.
Fuel propertiesASTM D6751EN 14214 [36]Density at 15°C (kg·m−3)880 [37]860–900Kinematic viscosity at 40°C (mm2·s−1)1.9–6.0 [37]3.5–5.0Cetane numberMin. 47 [37]Min. 51Flash point (°C)Min. 93 [37]Min. 101Sulfur % (mg/kg)Max. 0.05 [30]Max. 10Oxidation stability at 110°C (hr)Min. 3 [37]Min. 6


## 2.6. Validation of the Model

The validation of the model is done by comparing the model results with the ones of software CEA of NASA [28] and software GASEQ [27], which use element potential method and minimization of Gibbs free energy approach, respectively. The simulations of methane with steam injection are carried out. The mole fraction of chemical equilibrium products and thermodynamic properties obtained with the model compared with the ones of software CEA of NASA and GASEQ at ϕ = 0.6 with steam injection 10% are mentioned in Table 5 and at ϕ = 1.2 with steam injection 10% are listed in Table 6.Table 5
Comparison of the present model results with the ones obtained from CEA of NASA and GASEQ software for equilibrium combustion products of methane atϕ = 0.6 and steam injection 10%.
ϕ = 0.6Model resultsGASEQ resultsCEA resultsGASEQ relative error (%)CEA relative error (%)CO20.0515060.051540.051510.0660.066H2O0.2339200.233990.233850.0290.029N20.6455590.645530.64553−0.005−0.005O20.0683530.068270.06827−0.122−0.122CO1.328 × 10−71.33 × 10−73.27 × 10−7−0.157−0.157H22.419 × 10−72.41 × 10−72.41 × 10−7−0.360−0.360H1.269 × 10−91.27 × 10−91.27 × 10−90.030.03O1.506 × 10−71.50 × 10−71.50 × 10−7−0.429−0.429OH3.263 × 10−53.39 × 10−53.39 × 10−53.7513.751NO6.278 × 10−46.27 × 10−46.27 × 10−4−0.126−0.126Cp (kJ·kg−1·K−1)1.464881.453831.45383−0.76−0.76Tin (K)3003003000.0000.000Tad (K)1482.31482.31482.30.0000.000Table 6
Comparison of the present model results with the ones obtained from CEA of NASA and GASEQ software for equilibrium combustion products of methane atϕ = 1.2 and steam injection 10%.
ϕ = 1.2Model resultsGASEQ resultsCEA resultsGASEQ relative error (%)CEA relative error (%)CO20.0637290.063680.06368−0.077−0.077H2O0.2781730.278220.278220.0170.017N20.5947990.594800.594800.000170.00017O25.826 × 10−85.78 × 10−85.78 × 10−8−0.790−0.790CO0.0311380.031190.031190.1680.168H20.0321060.032060.03206−0.144−0.144H2.988 × 10−52.98 × 10−52.98 × 10−5−0.260−0.260O1.535 × 10−81.53 × 10−81.53 × 10−8−0.345−0.345OH2.199 × 10−52.95 × 10−62.68 × 10−62.6792.57NO2.962 × 10−62.95 × 10−60.646096−0.394−0.394Cp (kJ·kg−1·K−1)1.645921.630421.63042−0.95−0.95Tin (K)3003003000.0000.000Tad (K)1919.31919.31919.30.0000.000


## 3. Results and Discussion

### 3.1. Mole Fraction of Combustion Products

Adiabatic combustion simulation with steam has been carried out for the biodiesel and the diesel fuels in order to investigate the effects of steam injection on the equilibrium combustion products and thermodynamic properties. Before the combustion, the temperature of the air and the fuel is assumed to be equal to 300 K and the pressure of the air and the fuel mixture is equal to the combustion chamber pressure. The pressure inside the combustion chamber is 30 atm, which is supposed to be the pressure of injection steam. The temperature of injected steam is 300°C. The temperature of the unburned mixture is found by assuming a thermal equilibrium amongst the reactants. Initial temperature is obtained from the thermal equilibrium of fuel-air and steam mixture.The results have been comparatively presented for lean and rich combustion with increasing steam injection from 0% to 10%.Figures1(a) and 1(b) show the CO2 equilibrium mole fractions with steam injection ratios and they decrease remarkably with increasing steam injection ratios. NBD gives the highest results, while WCOBD presents the lowest results. Kayadelen [26] obtained the same behaviour in his study.Figure 1
Variation of CO2 equilibrium mole fractions with steam injection ratios.
(a)(b)Figures2(a) and 2(b) illustrate the H2O equilibrium mole fractions with steam injection ratios. It is clear that the H2O equilibrium mole fraction increases with increasing steam injection ratios. The maximum H2O equilibrium mole fraction is obtained with WCOBD and the minimum H2O equilibrium mole fraction is observed with diesel fuel. The reason is the WCOBD has the lowest α/β ratio and diesel fuel has the highest α/β ratio. The smaller the α/β ratio is, the larger the H2O equilibrium mole fraction is.Figure 2
Variation of H2O equilibrium mole fractions with steam injection ratios.
(a)(b)Figures3(a) and 3(b) illustrate the N2 equilibrium mole fractions with steam injection ratios. It is clear that the N2 equilibrium mole fraction decreases with increasing steam injection ratios. The maximum N2 equilibrium mole fraction is obtained with diesel fuel and the minimum N2 equilibrium mole fraction is observed with WCOBD. The reason is the WCOBD has the lowest α/β ratio and diesel fuel has the highest α/β ratio.Figure 3
Variation of N2 equilibrium mole fractions with steam injection ratios.
(a)(b)Figures4(a) and 4(b) shows the O2 equilibrium mole fractions with steam injection ratios. It is clear that the O2 equilibrium mole fraction decreases with increasing steam injection ratios in the case of lean combustion. The maximum O2 equilibrium mole fraction is obtained with diesel fuel and the minimum O2 equilibrium mole fraction is observed with WCOBD. The O2 equilibrium mole fraction is much lower in rich combustion due to lower oxygen concentration compared to lean combustion. The maximum O2 equilibrium mole fraction occurs with WCOBD and the minimum O2 equilibrium mole fractions occurs with diesel fuel in the case of rich combustion. The influence of the steam injection has the similar trend in Figures 5(a) and 5(b).Figure 4
Variation of O2 equilibrium mole fractions with steam injection ratios.
(a)(b)Figure 5
Variation of CO equilibrium mole fractions with steam injection ratios.
(a)(b)Figures6(a) and 6(b) show the CO equilibrium mole fractions with steam injection ratios. The CO equilibrium mole fraction decreases with increasing steam injection ratios. The maximum CO equilibrium mole fraction is observed with diesel fuel and the minimum CO equilibrium mole fraction is observed with WCOBD in the case of rich combustion. The maximum CO equilibrium mole fraction is obtained with NBD and the minimum CO equilibrium mole fraction is observed with WCOBD in the case of lean combustion.Figure 6
Variation of H2 equilibrium mole fractions with steam injection ratios.
(a)(b)Figures7(a) and 7(b) show the H2 equilibrium mole fractions with steam injection ratios which increase in both rich and lean combustions. The combustion of WCOBD gives the highest H2 equilibrium mole fraction. The influence of the steam injection has the similar trend in Figures 8(a) and 8(b).Figure 7
Variation of H equilibrium mole fractions with steam injection ratios.
(a)(b)Figure 8
Variation of O equilibrium mole fractions with steam injection ratios.
(a)(b)Figures9(a) and 9(b) demonstrate the OH equilibrium mole fractions with steam injection ratios. The OH equilibrium mole fractions increase slightly with increasing steam injection ratios in rich and lean combustion. WCOBD releases the highest OH equilibrium mole fractions while diesel fuel releases the lowest OH equilibrium mole fractions in both rich and lean combustions. Figures 10(a) and 10(b) demonstrate the NO equilibrium mole fractions with steam injection ratios. The NO equilibrium mole fractions decrease with increasing steam injection ratios in lean combustion. Diesel fuel releases the highest NO equilibrium mole fractions, while WCOBD releases the lowest NO equilibrium mole fractions in lean combustions.Figure 9
Variation of OH equilibrium mole fractions with steam injection ratios.
(a)(b)Figure 10
Variation of NO equilibrium mole fractions with steam injection ratios.
(a)(b)


### 3.2. Thermodynamic Properties of Combustion Products

The specific heat and the enthalpy of the mixture of the gases increase with the steam injection ratios due to the specific heat of superheated steam which is higher than the one of air at that temperature condition. Kayadelen [26] obtained the same behaviour in his study.Figures11(a) and 11(b) show the variation of specific heats of combustion products with respect to steam injection ratios. Higher specific heats are obtained in the rich combustions compared to the lean combustions. The highest specific heats are seen with WCOBD and the lowest specific heats are obtained with diesel fuel. In both combustions, specific heats increase with increasing steam injection ratios because the specific heat of the steam is much more than the ones of the combustion products at the same temperature except the specific heat of CO2.Figure 11
Variation of gas mixture specific heat with steam injection ratios.
(a)(b)Figures12(a) and 12(b) show the variation of enthalpy of combustion products with respect to steam injection ratios. Higher enthalpies of combustion products are obtained in the lean combustions compared to the rich combustions. The highest enthalpies of combustion products are seen with diesel fuel and the lowest enthalpies of combustion products are obtained with WCOBD. In both combustions, enthalpy of combustion products decreases with increasing steam injection ratios due to the enthalpy of the steam which is much less than the ones of the combustion products at the same temperature.Figure 12
Variation of gas mixture enthalpy with steam injection ratios.
(a)(b)Figures13(a) and 13(b) show the variation of entropy of combustion products with respect to steam injection ratios. Higher entropies of combustion products are obtained in the rich combustions compared to the lean combustions. The highest entropies of combustion products are seen with WCOBD and the lowest entropies of combustion products are obtained with diesel fuel. In both combustions, entropy of combustion products increases with increasing steam injection ratios.Figure 13
Variation of gas mixture entropy with steam injection ratios.
(a)(b)


## 3.1. Mole Fraction of Combustion Products

Adiabatic combustion simulation with steam has been carried out for the biodiesel and the diesel fuels in order to investigate the effects of steam injection on the equilibrium combustion products and thermodynamic properties. Before the combustion, the temperature of the air and the fuel is assumed to be equal to 300 K and the pressure of the air and the fuel mixture is equal to the combustion chamber pressure. The pressure inside the combustion chamber is 30 atm, which is supposed to be the pressure of injection steam. The temperature of injected steam is 300°C. The temperature of the unburned mixture is found by assuming a thermal equilibrium amongst the reactants. Initial temperature is obtained from the thermal equilibrium of fuel-air and steam mixture.The results have been comparatively presented for lean and rich combustion with increasing steam injection from 0% to 10%.Figures1(a) and 1(b) show the CO2 equilibrium mole fractions with steam injection ratios and they decrease remarkably with increasing steam injection ratios. NBD gives the highest results, while WCOBD presents the lowest results. Kayadelen [26] obtained the same behaviour in his study.Figure 1
Variation of CO2 equilibrium mole fractions with steam injection ratios.
(a)(b)Figures2(a) and 2(b) illustrate the H2O equilibrium mole fractions with steam injection ratios. It is clear that the H2O equilibrium mole fraction increases with increasing steam injection ratios. The maximum H2O equilibrium mole fraction is obtained with WCOBD and the minimum H2O equilibrium mole fraction is observed with diesel fuel. The reason is the WCOBD has the lowest α/β ratio and diesel fuel has the highest α/β ratio. The smaller the α/β ratio is, the larger the H2O equilibrium mole fraction is.Figure 2
Variation of H2O equilibrium mole fractions with steam injection ratios.
(a)(b)Figures3(a) and 3(b) illustrate the N2 equilibrium mole fractions with steam injection ratios. It is clear that the N2 equilibrium mole fraction decreases with increasing steam injection ratios. The maximum N2 equilibrium mole fraction is obtained with diesel fuel and the minimum N2 equilibrium mole fraction is observed with WCOBD. The reason is the WCOBD has the lowest α/β ratio and diesel fuel has the highest α/β ratio.Figure 3
Variation of N2 equilibrium mole fractions with steam injection ratios.
(a)(b)Figures4(a) and 4(b) shows the O2 equilibrium mole fractions with steam injection ratios. It is clear that the O2 equilibrium mole fraction decreases with increasing steam injection ratios in the case of lean combustion. The maximum O2 equilibrium mole fraction is obtained with diesel fuel and the minimum O2 equilibrium mole fraction is observed with WCOBD. The O2 equilibrium mole fraction is much lower in rich combustion due to lower oxygen concentration compared to lean combustion. The maximum O2 equilibrium mole fraction occurs with WCOBD and the minimum O2 equilibrium mole fractions occurs with diesel fuel in the case of rich combustion. The influence of the steam injection has the similar trend in Figures 5(a) and 5(b).Figure 4
Variation of O2 equilibrium mole fractions with steam injection ratios.
(a)(b)Figure 5
Variation of CO equilibrium mole fractions with steam injection ratios.
(a)(b)Figures6(a) and 6(b) show the CO equilibrium mole fractions with steam injection ratios. The CO equilibrium mole fraction decreases with increasing steam injection ratios. The maximum CO equilibrium mole fraction is observed with diesel fuel and the minimum CO equilibrium mole fraction is observed with WCOBD in the case of rich combustion. The maximum CO equilibrium mole fraction is obtained with NBD and the minimum CO equilibrium mole fraction is observed with WCOBD in the case of lean combustion.Figure 6
Variation of H2 equilibrium mole fractions with steam injection ratios.
(a)(b)Figures7(a) and 7(b) show the H2 equilibrium mole fractions with steam injection ratios which increase in both rich and lean combustions. The combustion of WCOBD gives the highest H2 equilibrium mole fraction. The influence of the steam injection has the similar trend in Figures 8(a) and 8(b).Figure 7
Variation of H equilibrium mole fractions with steam injection ratios.
(a)(b)Figure 8
Variation of O equilibrium mole fractions with steam injection ratios.
(a)(b)Figures9(a) and 9(b) demonstrate the OH equilibrium mole fractions with steam injection ratios. The OH equilibrium mole fractions increase slightly with increasing steam injection ratios in rich and lean combustion. WCOBD releases the highest OH equilibrium mole fractions while diesel fuel releases the lowest OH equilibrium mole fractions in both rich and lean combustions. Figures 10(a) and 10(b) demonstrate the NO equilibrium mole fractions with steam injection ratios. The NO equilibrium mole fractions decrease with increasing steam injection ratios in lean combustion. Diesel fuel releases the highest NO equilibrium mole fractions, while WCOBD releases the lowest NO equilibrium mole fractions in lean combustions.Figure 9
Variation of OH equilibrium mole fractions with steam injection ratios.
(a)(b)Figure 10
Variation of NO equilibrium mole fractions with steam injection ratios.
(a)(b)


## 3.2. Thermodynamic Properties of Combustion Products

The specific heat and the enthalpy of the mixture of the gases increase with the steam injection ratios due to the specific heat of superheated steam which is higher than the one of air at that temperature condition. Kayadelen [26] obtained the same behaviour in his study.Figures11(a) and 11(b) show the variation of specific heats of combustion products with respect to steam injection ratios. Higher specific heats are obtained in the rich combustions compared to the lean combustions. The highest specific heats are seen with WCOBD and the lowest specific heats are obtained with diesel fuel. In both combustions, specific heats increase with increasing steam injection ratios because the specific heat of the steam is much more than the ones of the combustion products at the same temperature except the specific heat of CO2.Figure 11
Variation of gas mixture specific heat with steam injection ratios.
(a)(b)Figures12(a) and 12(b) show the variation of enthalpy of combustion products with respect to steam injection ratios. Higher enthalpies of combustion products are obtained in the lean combustions compared to the rich combustions. The highest enthalpies of combustion products are seen with diesel fuel and the lowest enthalpies of combustion products are obtained with WCOBD. In both combustions, enthalpy of combustion products decreases with increasing steam injection ratios due to the enthalpy of the steam which is much less than the ones of the combustion products at the same temperature.Figure 12
Variation of gas mixture enthalpy with steam injection ratios.
(a)(b)Figures13(a) and 13(b) show the variation of entropy of combustion products with respect to steam injection ratios. Higher entropies of combustion products are obtained in the rich combustions compared to the lean combustions. The highest entropies of combustion products are seen with WCOBD and the lowest entropies of combustion products are obtained with diesel fuel. In both combustions, entropy of combustion products increases with increasing steam injection ratios.Figure 13
Variation of gas mixture entropy with steam injection ratios.
(a)(b)


## 4. Conclusion

The equilibrium combustion model is used to provide mole fractions and thermodynamic properties of combustion products in chemical equilibrium. Equilibrium mole fractions and the adiabatic flame temperature are necessary for estimating thermodynamic properties of exhaust gases and for providing key data to obtain nonequilibrium mole fractions. The precise calculation of thermodynamic properties is necessary for accurate performance estimations of internal combustion engines. The influences of steam injection on equilibrium combustion products and thermodynamic properties of diesel fuel and biodiesels have been modelled.The equilibrium mole fractions and thermodynamic properties of combustion products showed good agreement with the results obtained from CEA of NASA and GASEQ using methane.The highest CO2 equilibrium mole fractions are seen with NBD, and the lowest ones are formed with WCOBD. The reduction rate of NO is higher for biodiesel fuels as compared with the one of diesel fuel in lean combustion conditions. Thus, the steam injection ratios can be used for biodiesel fuels in order to lower NO equilibrium combustion products in the lean combustion conditions.The model can be a recourse tool for researchers studying fuels, surrogates, emissions, and internal combustion modelling due to its accuracy and convenience.The equilibrium mole fractions are used in combustion engineering, for example, in chemical kinetics to look for a reference point and in parametric emission monitoring systems to estimate the emissions of power systems.Simulations of any combustion process including exhaust gas recirculation, supplementary firing, steam injection, water injection, fuel emulsion, and reheat process in gas turbines can easily be conducted with higher accuracy and speed by integrating this model to the thermodynamic models.



---

*Source: 2805125-2020-09-03.xml*