Open access

Thermodynamic Study of the Working Cycle of a Direct Injection Compression Ignition Engine

Written By

Simón Fygueroa, Carlos Villamar and Olga Fygueroa

Submitted: 21 November 2011 Published: 14 November 2012

DOI: 10.5772/50028

From the Edited Volume

Internal Combustion Engines

Edited by Kazimierz Lejda and Pawel Wos

Chapter metrics overview

2,985 Chapter Downloads

View Full Metrics

1. Introduction

Currently, one of the worldwide most used energy sources are fuels derived from oil, such as hydrocarbons, that burn with oxygen releasing large amount of thermal energy. This energy can be transformed into mechanical work by mean of internal combustion engines [1]. Internal combustion engine is a device that allows obtaining mechanical energy from the thermal energy stored in a fluid due to a combustion process [2].

It should be noted that in reciprocating internal combustion engines (RICE) the combustion products constitute the working fluid; this simplifies their design and produces high thermal efficiency. For this reason these engines are one of the lighter weight generating units known and thus actually are the most commonly used transport engines [3].

The RICE operate following a mechanical cycle consisting of two main parts: the first one is the closed cycle, where the compression, combustion and expansion processes are carried out, and the second one is the open cycle where the working fluid is renewed, known as the gas exchange process and constituted by the intake and exhaust processes [4]. When RICE is study, it is mandatory to determine the working fluid thermodynamic properties, as well as the mixture amount that enters and leaves the cylinder [5].

The flow characteristics in spark-ignition engines (SIE) or compression ignition engines (CIE) can be summarized according to [6] as follows: transient as a result of piston movement, fully turbulent for all cylinders due to engine velocities and admission duct dimensions, and three-dimensional due to the engine geometry that also varies during the cycle (contours varying with time) producing different local velocity fields.

During the gas exchange ondulatory and inertial phenomena processes occur, as well as instability in the processes that occur within the engine. The variation of in-cylinder pressure during the intake and the exhaust has a complex pattern, for this reason the analytical calculation of gas exchange considering the above-mentioned phenomena is quite complicated and requires the use of specialized computer programs that use coefficients obtained experimentally [1].

The basis for calculation of non-stationary non-isentropic flow characteristic of RICE inlet and outlet ducts, and NO emissions was established in [7]. Various empirical correlations to take into account heat transfer during gas exchange process and to adjust the Reynolds number exponential factor in such a way that reduce to only one the adjustment coefficients were considered in[1,2,4 and 8].

A procedure widely used in both experimental and theoretical study of flow in engines, is to analyze the engine cycle in absence of combustion, simulating the process of compression-expansion and making measurements to the engine operating at this condition [5].

In all RICE working cycle processes, there is heat transfer to the cylinder walls, which occurs with greater intensity during the combustion and expansion due to the high temperature gradients reached. Woschni [9] proposed equations to determine turbulent convective heat transfer considering average speed of in-cylinder gases and Annand [10] find correlations to calculate instant average coefficients for turbulent convection heat transfer using gas average temperature and proposed correlations to evaluate flame radiation emitted during combustion. Correlations for convection heat transfer taking into consideration surface change and cylinder enclosed volume as piston moves was established in [11]. A computer program to calculate the heat transfer in a RICE combustion chamber using models to consider turbulence was presented in [12]. A universal correlation for mixture flow in the admission and exhaust process, correcting the coefficients of Nusselt, Reynolds and Prandtl numbers was proposed in [13].

In present study the compression process is considered adiabatic and reversible, but in real engines there is heat transfer between working fluid, valves and cylinder walls. At the beginning of compression the fluid temperature is lower than the temperature of the surfaces that surround the cylinder volume, causing an increase in the fluid temperature, some instant later temperatures get equal and latter on, heat is transferred from the working fluid to the walls, therefore the politropic coefficient varies during the process [1].

The complexity of the combustion process in RICE because of untimely and incomplete combustion, dissociation and heat transfer, has encouraged the development of special techniques for carrying out studies. Adequate realization of this process is decisive in terms of the produced power and its efficiency having great influence on the engine life and reliability [14].

Various models have been proposed to study the combustion process such as the Wiebe burning law applicable to SIE and the Watson law, applicable to CIE [15]. These laws determine burned mass fraction and heat released depending on the crankshaft rotation angle. These models used physical constants obtained experimentally. The Rasselier and Withrow relationship, along with burning laws allow obtaining combustion pressure per crankshaft rotation degree. To quantify ignition delay there are many correlations as the one proposed by [16], [17], [18] or [19]. Models proposed by [15], [20] and [17] are used to calculate the burning factor, which is representative of the burned mass fraction in the pre-mixed phase and the diffusive phase.

Volume changes may be evaluated using an expression proposed by [15], which correlates the engine dimensions: compression ratio, displaced volume, combustion chamber volume, crank radius, connecting rod length, and crankshaft rotation angle. Average temperature during combustion process can be determined using in-cylinder pressure and ideal gas equation [4].

Calculation methods used to obtain equilibrium composition and final state of chemical species present in the combustion products of a fuel air mixture are well known and are referenced in the literature [21, 22, 23 and 24]. One of the most complete programs is perhaps the NASA-Lewis code CEC [25 and 26] which considers liquid and gaseous chemical species, is extremely versatile and can be used to calculate thermodynamic state, chemical equilibrium, rockets theoretical behavior and even Chapman-Jouguet detonation properties.

Computer programs for calculating CHO and CHON constant pressure combustion systems assuming that combustion products were composed of eight and ten chemical species were presented respectively in [27 and 28]. A code less general than the NASA code, limited to twelve chemical species CHON combustion systems specifically designed to be applied to the analysis of internal combustion engines processes was published in [29]. A program for calculating twelve species CHON constant volume combustion systems, applicable to temperatures up to 3400 K was presented in [30]. A program, valid for temperatures up to 6000 K, which can be calculate, both constant pressure and constant volume combustion for an eighteen chemical species CHON system is available in [31].

The working fluid properties function of its temperature, pressure and richness can be determined by applying thermodynamic basic equations for ideal gas mixtures, considering the mass fractions of each component in the mixture [32]. Also can be determined through routines such as FARG and ECP [33 and 34]. To complement the study of the combustion process, models to determine NO emissions as the extended Zeldovich mechanism, have been considered. The reason for the use of these models is because the specific constants of the reaction rate for NO are very small compared to the combustion rate, for this reason it is supposed that all the species present in the products with the exception of NO, are in chemical equilibrium.

The expansion process produces mechanical work from energy released during combustion and ends when exhaust valve opens. At this moment products are expelled from the cylinder initially at critical speed ranging between 500 and 700 m/s, and then are pushed by piston movement towards the upper dead center [4 and 15]. Towards the end of the exhaust during the valve overlap, part of the fresh mixture escapes contributing to the emission of unburned hydrocarbons and reducing the engine efficiency.

To investigate the gas exchange process, using gas dynamics to analyze the gas flow in transient processes with variable composition and variable specific heats, models such as [35] have been used.

To improve the gas exchange process we must advance inlet valve opening (AIVO) and delay exhaust valve closure (DEVC). Because of this, there is a period in which the two valves remain open simultaneously, this period is known as valve overlap, which helps to remove as much gas and admit as much air or fresh mixture. This is due to the depression originated in the inlet valve vicinity, due to the ejection effect produced by burned gas movement through the exhaust valve; this will contribute to increase efficiency and power produced by the RICE [1].

Two research methods are employed to study the working cycle of RICE. The first one is based on data acquisition from experimental tests and the second one is based only on mathematical simulation. The latter method is more versatile and reduces the required research empirical data depending on the employed calculating method and imposed simplifications. However, to validate mathematical simulation results, experimental parameters obtained in laboratory are required [5]. Use of numerical analysis methods has now greatly developed and increased due to increase of velocity and calculation capacity of modern computers. These methods provide faster performance, versatile and can handle more information than can be measured in an experimental test. However, results accuracy obtained by applying the models depends on made assumptions.

Modeling is a research technique employed in RICE, its use has grown in last two decades due to the cost decrease obtained by eliminating or reducing the laboratory tests, since these require a large amount of repetitive tests to obtain appropriate results, bringing time and money losses in preparing, calibrating, measuring, repairing and replacing testing engines. RICE designers must build more efficient engines due to higher fuel cost and new regulations on combustion emissions produced in the process that occurs inside engine. In order to optimize these designs, numerous trial and error tests are required. Implementing the tests implies expensive construction and testing of several prototypes. Modeling is a procedure that allows realizing numerous tests with relatively low cost.

To determine engine p vs. V diagram, working fluid is considered as an ideal gas, the mass entering the cylinder is calculated using a filling model that takes into account valve rise and discharge coefficient. Initial in-cylinder mass are residual gases, same quantity that was used as a reference value to control the expelled mass during exhaust. Instant volume was determined using an equation in terms of crankshaft rotation angle [15]. Final compression temperature was found from the first law of thermodynamic considering a uniform flow process and convection heat transfer. Power, mean indicated pressure and maximum pressure and temperature were calculated using the methods proposed in [1], [4] and [15]. Cyclic dispersion were studied using mean indicated pressure variation coefficient and pressure variation as a function of main combustion phase angle in the range of 10° before top dead center (TDC) and 10° after TDC, [1]. Calculations for exhaust process were similarly to those of the admissions process. A model to study the closed loop of a CIE limited pressure cycle, replacing the constant volume heat rejection process by an isentropic expansion process followed by a constant pressure heat rejection is proposed in [36].

There are commercial packages which constitute a very useful tool in the field of research and development of RICE being employed by different companies in the automotive sector. These include the ECARD (Engine Computer Aided Research & Development), developed by the IMST group, a global model that allows simulating engine operation throughout its working cycle, using similar complexity models for the different processes involved. The OpenWAM is a free, open source 1-dimensional gas-dynamics code produced by CMT group that can be used to predict the flow movement through the elements of an internal combustion engine. NEUROPART, uses neural networks to determine the product properties and composition influence on the exhaust emissions and particles formation. CHEMKIN, uses the chemical kinetic concepts to analyze fluids in gas phase through fluid dynamic simulation. EQUIL, calculates the composition at equilibrium of combustion products. PREMIX, calculates combustion speed for different fuels. SENKIN, allows to determine the time delay for different fuels and the combustion kinetic evolution depending on the species involved in a process.

Advertisement

2. Mathematical model

Present paragraph will develop fundamentals and mathematical equations that govern the phenomena occurring in CIE. For this purpose, volume control in Figure 1, which shows the mass and energy interactions with surroundings will be considered.

Figure 1.

Engine control volume

It should be noted that the control volume during gas exchange processes works as an open system. During compression, combustion and expansion processes works as a closed system, that is why corrections should be made to take into account the exchanged mass due to leakage and the supply of fuel.

2.1. Mass conservation

The mass conservation principle establishes that the total mass change in a control volume is:

mvc=memsE1

The summation is used when there are several inputs and/or output flows. Expressing Ec. 1 in differential form and dividing by a time differential we obtain the mass time rate of change:

dmvcdt=dmedtdmsdtE2

To express last equation in terms of air and fuel mass entering the control volume, we define:

f=mfmE3

Time differentiating and rearranging we obtain the fuel time rate of change:

f=dfdt=(memsm)[(fefs)]E4

From equivalence ratio (mixture richness) definition:

ϕ=mfma(mfma)sto=mfma(mfma)stoE5

replacing Ec. 3 in Ec. 5 and time deriving:

ϕ=dϕdt=1(mfma)stof(1f)2E6

2.2. Energy conservation

The first law of thermodynamics for an open system, disregarding changes in kinetic and potential energy can be written in differential form as:

dEdt=dQdtdWdt+mehemshsE7

Since the work due to a volume change is:

dWdt=W=PdVdtE8

and the first term on the left-hand side of Eq. 7 can be evaluated in terms of the internal energy:

dEdt=ddt(mu)=(mdudt)vc+(udmdt)vcE9

or in terms of enthalpy:

dEdt=ddt(mh)ddt(pV)E10

Substituting Eqs. 8 and 10 in Eq. 7 we have:

(mdudt)vc+(udmdt)vc=QW+mehemshsE11

Since internal energy, enthalpy and density are T, p and ϕ functions, their time rate of change are:

dudt=(uT)dTdt+(up)dpdt+(uϕ)dϕdtE12
dhdt=(hT)dTdt+(hp)dpdt+(hϕ)dϕdtE13
dρdt=(ρT)dTdt+(ρp)dpdt+(ρϕ)dϕdtE14

Assuming the working fluid is an ideal gas, differentiating ideal gas equation and rearranging we have:

pdVdT+VdpdT=mRdTdt+mTdRdt+RTdmdtE15
dρdt=dpdtρRdTdtρTdRdtRTE16

From Eq. 14:

dpdt=dρdt(ρT)dTdt(ρϕ)dϕdt(ρp)E17

substituting Eq. 17 in Eq. 16, rearranging and solving for dpdT:

dpdt=ρTdTdtρRdRdt(ρT)dTdt(ρϕ)dϕdt(ρp)1RTE18

Solving Ec. 15 for dRdt, simplifying and substituting in Ec. 18:

dpdt[(ρp)1RT]=ρTdTdtρR[pmTdVdt+VmTdpdtRTdTdtRmdmdt](ρT)dTdt(ρϕ)dϕdtE19

Replacing ideal gas equation in Ec. 19 and solving for dpdT:

dpdt=ρ(ρp)[dVdtV+dmdtm1ρ(ρT)dTdt1ρ(ρϕ)dϕdt]E20

Differentiating ideal gas equation with respect to p and T, we obtain:

(ρp)=1RTE21
(ρT)=pRT2E22

and substituting Ecs. 21 and 22 in Ec. 20:

dpdt=p[dVdtV+dmdtm+(dTdt)1TRTp(ρϕ)dϕdt]E23

Later equation, function of density, volume, mass and mixture richness will be used to obtain the in-cylinder pressure when time varies (indicator diagram).

The procedure to obtain a similar expression for temperature variation over time will be illustrated below. Solving Eq. 11 for dudt:

dudt=QmvcpmvcdVdt+1mvc(mehemshs(udmdt)vc)E24

and defining:

B=RTdVdtV+1m(Q+mehemshs(udmdt)vc)E25

On the other hand, introducing Ecs. 8 in Ec. 11, get the following expression:

dudt=BE26

Substituting Ec. 26 in Ec. 12 and solving for dpdt:

dpdt=B(udt)dTdt(udϕ)dϕdt(up)E27

Replacing Ec. 27 in Ec. 20 and solving for dTdtgives:

dTdt=(up)[ρdVdtV+ρdmdtm(ρϕ)dϕdtB+(uϕ)dϕdt](up)(ρT)(ρp)(uT)E28

Now, considering:

R=R(T,p,ϕ,t)E29

And differentiating:

dRdt=(RT)dTdt+(Rp)dpdt+(Rϕ)dϕdtE30

Differentiating ideal gas equation and solving for dRdt

dRdt=dpdtpRdTdtTRdρdtρRE31

Replacing Ec. 30 in Eq. 31 and solving fordpdt:

dpdt=(RT)dTdt+(Rϕ)dϕdt+dTdtRT+dρdtRρRp(Rp)E32

Replacing Ec. 32 in Ec. 27:

(uT)dTdt+(up)[(RT)dTdt+(Rϕ)dϕdt+dTdtRT+dρdtRρRp(Rp)]+(uϕ)dϕdt=BE33

Defining:

D=1pR(Rp)E34

Collecting terms containing dTdt and substituting Eq. 34 in Eq. 33:

(dTdt)[(uT)+(up)pDR{(RT+RT)}]+(up)1D[(Rϕ)dϕdt+dϕdtRρ]+(uϕ)dϕdt=BE35

Defining:

C=1+TR(RT)E36

Replacing Ec. 36 in Eq. 35 and solving fordTdt:

dTdt=B(up)pD[1R(Rϕ)dϕdt+dρdt1ρ](uϕ)dϕdt(uT)+CDpT(up)E37

Since:

m=ρVE38

Time differentiating and solving for dρdt:

dρdt=dmdtρdVdtVE39

Replacing Ec. 39 in Ec. 37:

dTdt=B(up)pD[(Rϕ)dϕdt1R+dmdt1mdVdt1V](uϕ)dϕdt(uT)+CDpT(up)E40

This equation will be used to determine the in-cylinder temperature when time varies.

If Eqs. 25, 34, 36 and 37 are replaced in Eq. 23 and collecting terms we obtained:

dpdt=QmbbhbbdVdt[mCV+p]dmdt[(D(uϕ)hcil+u)C(DR(Rϕ+1))]m[C(1p1R(Rϕ))+(up)]E41

and:

dmdt=QmbbhbbdVdt[mCV+p]m[C(1p1R(Rϕ))+(up)]dpdt(D(uϕ)hcil+u)C(DR(Rϕ+1))E42

Eqs. 41 and 42 will be used to obtain the indicator diagram (p vs. V or p vs. φ diagram) and burned mass fraction diagram (m vs. t diagram) respectively.

2.3. Instant in-cylinder volume

The instant volume inside the control volume in terms of the displaced volume, compression ratio, connecting rod length to crank radius ratio and crankshaft rotation angle, can be obtained through the following expression [15]:

V(φ)=Vd[1rc1+12[RLA+1cosφ(RLA2sen2φ)12]]E43

In this expression RLA is the connecting rod length (l) to crank radius (a) ratio.

RLA=laE44

Deriving Ec. 43 with respect to crankshaft rotation angle, we obtain:

dVdφ=Vd2[senφ+senφcosφ(RLA2sen2φ)12]E45

The time in seconds it takes to describe some crankshaft rotation angle can be calculated with the following expression:

t=φ6rpmE46

Solving previous expression for φ and replacing in Ec. 45 to make the corresponding conversion from degrees to radians we obtained:

dVdt=3Vd(rpm)[senπrpm30t+senπrpm30tcosπrpm30t(RLA2sen2πrpm30t)12]E47

Previous expression allows determining the in-cylinder volume variation with respect to time, while Ec. 45 will be used to calculate the volume variation with respect to crankshaft rotation angle.

Advertisement

3. Equations, models and calculations

Models and assumptions used to analyze each of the thermodynamic processes that are carried out in a CIE will be presented in this paragraph. Routines commonly used in RICE are employed to calculate the thermodynamic properties of the chemical species formed during combustion. FARG and ECP routines [34] are used to determine properties depending on gas temperature. PER routine [29] is used to obtain same properties depending on mixture richness. DVERK routine [37] found in the International Mathematics and Statistics Library is used for solving differential equations systems by the Runge - Kutta Verner fifth and sixth order method.

3.1. Admission process

The parameter characterizing the admissions process is the volumetric efficiency defined as:

ηv=marmat=marρ0iVdrpm30jE48

It takes into account the losses in the inlet valve and all the admission system if atmospheric density value is used for ρ0.

Real mass air flow entering the cylinder is determined by the following equations [15] function of the ratio pdown/pup:

pdownpup<1m=CdArefp0RaT0(pdownpup)(2γγ1[1(pdownpup)γ+1γ])E49
pdownpup1m=CdArefpup(γRT0)(2γ1)γ+12(γ+1)E50

where pdown is the downstream pressure and pup is the upstream pressure. Although the discharge coefficient Cd varies during the process, in present study, we assume that it is constant and equal to its average value. The reference area Aref, usually called curtain area, since it depends on the valve lifting Lv, is taken as:

Aref=πdvLvE51

A model proposed in [38], was used to theoretically determine the lifting valve profile which is function of the maximum lifting, Lv max and the crankshaft rotation angle φ:

Lv(φ)=Lvmax+C2φ2+Cpφp+Cqφq+Crφr+CsφsE52

Coefficients C2, Cp, Cq, Cr y Cs are determined with the following equations:

C2=pqrsh[(p2)(q2)(r2)(s2)cmed2]E53
Cp=2qrsh[(p2)(qp)(rp)(sp)cmedp]E54
Cq=2prsh[(q2)(qp)(rq)(sq)cmedq]E55
Cr=2pqsh[(r2)(rp)(rq)(sr)cmedr]E56
Cs=2pqrh[(s2)(sp)(sq)(sr)cmeds]E57

Recommended values for p, q, r, and s are: p = 6; q = 8; r = 10; s = 12.

Gas pressure and temperature variation over time in this process is calculated from Eqs. 23 and 40. Since a CIE only compresses air, the term corresponding to mixture richness variation with time is zero. For this reason the above equations are:

dpdt=p[dVdtV+dmdtm+(dTdt)1T]E58
dTdt=B(up)pD[dmdt1mdVdt1V](uT)+CDpT(up)E59

To solve these equations heat release, heat transfer, blow by, ignition delay and chemical species formation models are required. In addition the termsuT,up, RT y Rp must be determined by using routines FARG and ECP.

With Eqs. 43 and 47 we calculate the volume and its time derivative, respectively, while with Eqs. 49 or 50 depending on the case determine the flow mass. The in-cylinder accumulated mass is obtained by summation of the mass flows times the values obtained by Eq. 46.

3.2. Closed loop cycle

Closed loop cycle corresponds to compression, combustion and expansion processes. Compression starts when inlet valve closes. Variation in pressure and temperature with time during this process is determined taking into account that only air is compressed (Ecs 58 and 59) and there are mass losses due to blow by. When fuel injection begins the mixture composition varies; therefore the expressions used to determine temperature and pressure time variation during the closed loop cycle are Eq. 23 and 40. When exhaust valve opens, begins the exhaust process.

3.3. Exhaust process

Equations used in this process are the same used during the admissions process, but noting the working fluid is a mixture of burned gases and heat transfer is higher than during admission because of the high temperature present.

During valve overlap we want to extract as much as possible of burned gases and taking advantage of the dynamic effects, increase the amount of fresh charge entering the cylinder. Equations used during this process are the same used during the intake and exhaust, but considering that there is simultaneously fresh charge entry and burned gases exit.

3.4. Ignition delay model

Ignition delay in CIE characterizes the heat quantity will release immediately occur the fuel auto ignition and has a direct influence on engine rumble and pollutants formation. The model presented in [19] points out that ignition delay depends on in-cylinder temperature and pressure, engine speed and accumulated fuel amount and may be calculated in degrees and in milliseconds with the following expressions:

ID[deg]=Apcn(EA)E60
ID[ms]=ID[deg]0.006(rpm)E61

where: A=0.36+0.22Vmp,Vmp=c(rpm)30,n=0, EA=exp[Ru(1RTc-117190)(21.2pc-12.4)],E=618840NC+25,pc=pambrcnc,Tc=Tambrcnc-1,

nc=1.30a1.37,Ru=8.3143[J/molK].E62

Other models based on experimental data suggest correlations that use an Arrhenius expression similar to that proposed in [15], in which the constants estimated by [39] are as follows: A=3.45,n=1.02,EA=exp[EaRuTc],ERu=2100.

Another model whose constants are the same as in previous case use an A term function of richness as shown by the following expression [16]: A=2.4ϕ-0.2.

3.5. Heat release model

Considering the fourth term numerator in Eq. 41, which represents the heat released during the combustion process and applying the Watson relationship we obtain the following equation:

mcHidXb=dmdt[(D(uϕ)hcil+u)C(DR(Rϕ+1))]E63

Model of fuel apparent burning will be used to represent the combustion process. It uses two empirical equations, one for the pre-mixed combustion phase and another for the diffusive combustion phase. The instantaneous total amount of heat released by crankshaft rotation degree is given by the sum of the two components:

(dmcdφ)Tot=(dmcdφ)pre+(dmcdφ)difE64

3.6. Burning factor

Heat release model requires defining, depending on the process physical condition, the initial fuel amount burned during the pre-mixed phase. An initial fuel burning factor is used for this purpose [15] [17]. This factor estimated, depending on initial richness and delay period, what part of the injected fuel is burned during the pre-mixed phase. The difference is burned during the diffusive phase. The burning factor is defined as:

β=(mc)pre(mc)TotE65

and may be calculated by the following expression [15]:

β=1a1ϕ0b1IDScc1E66

The a1, b1 and cc1 values which are shown in Table 1 [39, 19 and 15] depend on the used model.

ValueHardenberg modelWatson modelHeywood model
a10.7460.9260.80 - 0.95
b10.350.370.25 - 0.45
cc10.350.260.25 – 0.50

Table 1.

Empirical values for burning factor

Taking into account the heat released during each phase Ec. 63 becomes

(dmcdφ)Tot=β(dmcdφ)pre+(1β)(dmcdφ)difE67

Heat released during each phase is evaluated by the empirical expressions proposed in [39] and [40]. Equations proposed in [39] are:

(dXbdφ)pre=C1C2(φφ0Δφ)C11(1(φφ0Δφ)C1)C21E68
(dXbdφ)dif=C3C4(φφ0Δφ)C41exp(C3(φφ0Δφ)C4)E69

Equation proposed in [40] that use the duration and heat released percentage in each phase, unlike [39] equations, is:

dQdφ=a(Qpreφpre)mpre(φφpre)mpre1exp[a(φφpre)mpre]+(Qdifφdif)mdif(φφdif)mdif1exp[a(φφdif)mdif]E70

Table 2 shows the constants for Heywood [15] and Watson [39] equations and Table 3 shows the constants for Miyamoto [40] equation.

ConstantHeywoodWatson
ExpresiónValores sugeridos
C12+1.25E8(ID[ms]rpm)2.42+0.002703*IDa2.43
C2500050005000
C314.2ϕ0.64410ϕ1.5056.908
C40.79C30.250.48C30.4231.4

Table 2.

Constants for Heywood and Watson expressions

ConstantMiyamoto
mpre4
mdif1.5 ID o 1.9 ID
a0.69
φpre+7°

Table 3.

Constants for Miyamoto expression

3.7. Blow by model

Due to the in-cylinder pressure increase during compression and combustion processes, a part of the gasses (mbb) is lost through the rings resulting in a produced power reduction. The model that will be described below takes into account such losses.

From logarithmic derivative of the ideal gas equation:

1pdpdφ+1VdVdφ=1mdmdφ+1TdTdφE71

Applying the First Law of Thermodynamics in differential form to an open system, we obtain:

mCvdTdφ+CvTdmdφ=dQdφpdVdφmbbCpTωE72

Replacing dT/dφ from Ec. 70 and rearranging:

dpdφ=γpVdVdφ+γ1VdQdφγmbbωmpE73

Blow by mass can be found from mass conservation:

dmdφ=mbbωE74

Defining:

C=mbbmE75

Net heat entering the system, considering that a part is lost by blow by, is:

dQ=QedxdQlostE76

with:

dQlostdt=hcA(TTwall)E77

and defining the following dimensionless variables:

p˜=pp1E78
V˜=VV1E79
T˜=TT1E80
Q˜=QeP1V1E81
Q˜lost=QlostP1V1E82
h˜c=hT1(Acc4VccDp)P1V1ωE83
β=4V1Dp(Acc4VccDp)E84

dimensionless heat losses for unit crankshaft rotation angle are:

dQ˜lostdφ=h˜c(1+βV˜)(p˜V˜T˜wall)E85

Replacing Ecs. 77 to 83 in Ecs. 72, 73 and 84 and applying expansion work definition, we have the system of equations:

dp˜dφ=γp˜V˜dV˜dφ+γ1V˜[Q˜dxdφh˜c(1+βV˜)(p˜V˜m˜T˜wall)]γCp˜ωE86
dQ˜lostdφ=h˜c(1+βV˜)(p˜V˜m˜T˜wall)E87
dmdφ=CmωE88
dW˜dφ=p˜dV˜dφE89

Solving this ordinary differential equations system one can calculate the mass lost by leaks, lost heat and power, and pressure variation.

In-cylinder mass varies with time and can be calculated depending on the speed of the engine. Replacing Ec. 73 in Eq. 87 and solving form¯we obtain:

m¯=exp[Cbb(φ+π)π(rpm)30]E90

3.8. Heat transfer models

Calculation of instant heat transfer in RICE is a complex problem. Expressions for calculating the total flow heat by combined conduction, convection and radiation requires to use empirical correlations [41]. One of the more used is based on the relationship between the Nusselt number and the Reynolds number for forced convection:

Nu=aRebE91

Replacing the Nusselt and Reynolds numbers in Eq. 112, we have:

h¯cLk=a(ρVLμ)bE92

In this expression L represents the characteristic length, which in RICE is the piston diameter, a magnitude that does not vary.

Considering also separately the radiation heat transfer, gets the expression to evaluate the total heat flow [41]:

qA=akDpRebΔT+c(Tgas4Twall4)E93

In these expressions a, b and c are constants that usually take the following values: a = 0.35 – 0.80, b = 0.70 – 0.80 and c = 0.576 σ, where σ is the Stefan – Boltzmann constant.

An expression used to calculate convection heat transferred is based on the model proposed in [18]:

Qlost=h¯c*A*(T¯gasT¯wall)E94

whereh¯crepresents the overall convection heat transfer coefficient that can be determined from the expression:

h¯c=a*Dpm1pmTgas0.751.62mwmE95

The constant values in this equation are: a = 0.13 and m = 0.8. Pressure must be in bar and the temperature in kelvin.

w is the average gas speed calculated by the following expression:

w=[C1wVmp+C2wVdTrefprefVref(ppnocomb)]E96

Subscript ref is used for a reference state that may be the compression process beginning. The empirical constant values depend on the process are: C1w = 6.18 during gas exchange process, C1w = 2.28 during compression, combustion and expansion processes, C2w = 0.0 during gas exchange and compression processes and C2w = 3.24E-3 during combustion and expansion processes.

Another expression used to consider the convection heat transfer process of RICE proposed in [19] is:

h¯c=C1HVd0.06p0.8T0.4(Vmp+C2H)0.8E97

where C1H and C2H are empirical constants which take into account local variations due to intake turbulence or combustion chamber geometry of the, their values are: C1H = 130E-4 and C2H = 1.40. Eq. 96 considers the increase in gas speed due to engine velocity and uses as characteristic length the volume rather than the piston diameter, as it is the case when using expressions proposed in [18].

3.9. In-cylinder mass

Total mass that fills the cylinder is composed of air, fuel and residual gases. Fuel enters the cylinder when compression process end and fuel injection begins. Residual gases have two components: gases remaining in the cylinder at the end of the exhaust process and recirculated gases entering the cylinder with admitted air as a pollution control measure. There is also a mass lost by blow by. Therefore we have:

mtotal=maire+mgr+mf+mEGRmbbE98

The theoretical incoming air mass is calculated by:

maire=ρVdE99

While the real mass entering the cylinder is obtained by the expression:

maire=madmtE100
madm is calculated with Eqs. 49 or 50 and time t is obtained from Eq. 46.

3.10. Combustion products model

This model allows determining theoretically the composition at equilibrium and combustion product thermodynamic properties of a fuel-air mixture whereas reactive products consist of ten chemical species. The combustion global equation for a ten chemical species system is:

mcCmHnOl+n+m4+l2ϕ(O2+3.7N2)y1H2O+y2H2+y3OH+y4H+y5N2+y6NO+y7CO2+y8CO+y9O2+y10OE101

Applying conditions of mass conservation to elements C, H, O and N and reactants we obtain five equations with eleven unknowns (product molar fractions and fuel mass) as shown in following expressions:

C:y1+y8=nE102
H:2y1+2y2+y3+y4=mE103
O:y1+y3+y6+2y7+y8+2y9+y10=l+2(n+m412)ϕE104
N2:2y5+y6=7.428(n+m412)ϕE105
i=110yi=1E106

Applying the chemical equilibrium to combustion reaction yields six additional algebraic equations to close the system:

H22H K1=y42y2(pp0)E107
O22O K2=y42y9(pp0)E108
H2+O22OH K3=y32y2y9E109
12O2+12N22OHK4=y62y912y512E110
H2+12O2H2O K5=y1y2y912(pp0)12E111
CO+12O2CO2 K6=y7y8y812(pp0)12E112

In Eqs. 101 to 111 p is the products pressure, y1 to y10 are the species molar fractions and K1 to K6 are the equilibrium constants which are function of temperature. The system is solved by calculating the partial derivatives with respect to temperature, pressure, and mixture richness through FARG, ECP y PER routines. Applying the first law of thermodynamics to reactants and products depending on the type of combustion process we can find their thermodynamic properties. For constant volume combustion and constant pressure combustion respectively:

UR=UPE113
HR=HPE114

3.11. NOx formation model

The nitrous oxides (NOx) produced by a RICE are NO and NO2 and can be calculated by Zeldovich mechanism. This theory postulates that the production of oxides during the combustion process can be explained through the following chemical reactions:

R1:O+N2NO+NE115
R2:N+O2NO+OE116
R3:N+OHNO+HE117

A differential equation which allows finding the NO concentration by unit time during combustion, using Zeldovich reactions and considering the possibility of occurrence in both directions can be obtained using the basic theory of chemical kinetics:

d[NO]dt=k1+[O][N2]+k2+[N][O2]+k3+[N][OH]k1[NO][N]k2[NO][O]k3[N][OH]E118

Considering the steady state condition for [N], which is equivalent to assume its change is small compared to other interesting species and proceeding in similar mode to the previous case, we obtained:

d[N]dt=k1+[O][N2]k2+[N][O2]k3+[N][OH]k1[NO][N]+k2[NO][O]+k3[NO][H]E119
d[N]dt=0E120

Substituting Ec. 118 in Ec.117 and using the equilibrium reactions (Ecs. 114 to 116) gets the following expression:

d[NO]dt=2R1{1([NO]/[NO]e)2}1+([NO]/[NO]e)R1(R2+R3)E121

and Eqs. 114 to 116 are transformed into:

R1:k1+[O]e[N2]eE122
R2:k2[NO]e[O]eE123
R3:k3[NO]e[H]eE124

Terms enclosed in square brackets represent the chemical equilibrium concentration for corresponding species obtained using routines FARG and ECP. Constants for Zeldovich reactions, depending on temperature, are expressed in the following way.

K1+=7.6E13exp(38000/T)E125
K2=1.5E09exp(19500/T)E126
K3=2.0E14exp(23650/T)E127

Positive sign indicates reaction tendency to form products, while negative sign indicates tendency to form reactants.

3.12. Exhaust gas recirculation (EGR) model

Exhaust gas recirculation (EGR) is a very effective technique employed to reduce nitrogen oxide emissions (Lapuerta, 2000). The method involves replacing a part of the intake air with exhaust gases during the admissions process. In this way a decrease in the amount of air available for combustion reduces the final temperature of combustion, and therefore lowers the production of nitrous oxides emissions. This process can be expressed in the following way:

Fuel+ Air + %EGRProductsE128

Recirculated gases will be taken at engine exit and will be introduced into the cylinder at a temperature greater than air admitted. Composition and properties calculation of involved species will be made with routines FARG and ECP.

3.13. Work, mean effective pressure and efficiency

To determine the work produced by closed loop cycle, we consider that it is represented by the loop enclosed area in a p-V diagram. The area was calculated with the trapezoidal rule [42]:

A=W={p0ΔV+p1ΔV+p2ΔV+p3ΔV........pn1ΔV+pnΔV}E129

Gas exchange work, indicated work and net work are obtained with the following equations:

Wpum=WexhWadmE130
Wi=WexpWcomE131
Wnet=WiWpumE132

Net power is obtained through the following expression:

Wnet=Wnetrpm30jE133

Mean indicated pressure, mean net indicated pressure and mean indicated gas exchange pressure are calculated with equations:

mip=Wi/VdE134
mipnet=Wnet/VdE135
mippum=Wpum/VdE136

Efficiencies are determined with:

ηi=WimcHiE137
ηinet=WnetmcHiE138

3.14. MECID computer program

MECID program consist of a routines group for determining: optimum lifting of intake and exhaust valves, gas leakage across piston rings (blow by), composition variation of chemical species that constitute the working fluid, effect of heat transfer, ignition delay, heat release, residual gases and recirculated gases, as well as nitrous oxides formation and pressure drop during intake and exhaust processes due to valve contraction. To calculate all parameters it uses the mathematical relationships and models proposed in preceding paragraphs. Figure 2 shows the MECID program flow diagram.

Using the ENGINE routine the program requested information on engine works site conditions (temperature and pressure) and on engine geometric characteristics (piston bore and stroke; valves diameter, maximum lift and average discharge coefficient, compression ratio, connecting rod length to crank radius ratio and engine speed) With VALVE routine requested information on calculation conditions (calculation angular interval, advance and delay of intake and exhaust valve closure, richness, residual gases ratio, recirculated gas percentage and cylinder wall average temperature). Selected heat transfer, ignition delay and heat release models which will be used in the application, as well as the angles when combustion begins and ends, the program allows to perform studies to determine what should be the working conditions to ensure maximum fuel energy utilization the influence of different variables on engine operation. Results of more important studies will be presented in next paragraph.

Figure 2.

Block diagram and routines used by MECID computer program

Advertisement

4. Results and discussion

4.1. Effect of burning process beginning on CIE operating parameters

For combustion duration constant and varying combustion start we obtained the p and T vs. φ diagrams shown in Figure 3. We note that if the combustion process is advanced too much, from 360° to 310°, it reached a high maximum pressure while the piston is at some point in the compression stroke. This brings as a consequence, greater work consumption and reduced cycle efficiency. If combustion beginning delays too much, 360° to 410°, the maximum pressure is very low. If combustion beginning advance is such that the process occurs around TDC we obtain intermediate maximum pressures during expansion process and a higher maximum temperature, representing greater thermal energy available to do work. In view of the above, for subsequent studies we will take a 345° angle as combustion beginning angle. Maximum temperature is obtained for a greater angle than the maximum pressure one because the burning model uses constant combustion duration.

Figure 3.

p vs φ and T vs φ diagrams for various combustion beginning angles

Figure 4 shows how varies the accumulated burned mass fraction and the instant burned mass fraction depending on combustion start angle. Burned mass is responsible for pressure and temperature rapid increase during combustion process. It is observed that peak values are the same for all three cases considered. This is because richness and volumetric efficiency remain unchanged.

Figure 4.

Accumulated and instant burned mass fraction for various combustion beginning angles

4.2. Effect of combustion process duration on CIE operating parameters

In present case, the start of combustion process is kept constant. Figure 4 shows that at lower combustion duration the greater the maximum pressure reached since a great amount of heat is released in a short time period, as seen by the burned curves slope. In these curves can be noted that in considered cases, the maximum burned fractions do not change because the richness was kept constant and the obtained volumetric efficiency remained almost constant, therefore the total in-cylinder mixture mass is the same for all referred cases.

Figure 5.

Pressure and accumulated burned fraction for different combustion duration

Table 4 presents the principal results obtained for operation parameters. From Table 4 it can be concluded that the more efficient cycle producing higher power, is one in which combustion takes 45° and therefore this duration will be used in future studies.

Combustion process duration30º45º60º
ηv0.720.720.72
Wi [kJ]0.600.730.75
Wnet [kJ]0.570.650.64
Wnet [kW]9.5610.9810.72
mip [kPa]994.751200.211232.88
ηi net0.350.400.39
Pmax [kPa]7689.177134.936614.27
Tmax [K]2235.642400.712496.92

Table 4.

Summary of main results obtained when varying combustion process duration.

4.3. Effect of angular velocity on CIE operating parameters

Pressure variation depending on crankshaft rotation angle when engine speed change, is shown in Figure 6. It can be noted from this figure that higher pressures are reached when the engine operates at lower speeds. This is because while the engine works at lower speeds, greater will be its volumetric efficiency. If the richness is kept constant, the greater the mass which is burned during combustion process and thus more energy will be released.

Figure 6.

In-cylinder pressure for various engine speeds

Engine speed1500 rpm2000 rpm2500 rpm
ηv0.920.690.55
β50.5745.3340.89
Wi [kJ]1.421.181.03
Wnet [kW]14.6316.0117.28
mip [kPa]2314.471922.251678.68
ηi0.710.740.77
pmax [kPa]6703.236034.925506.50

Table 5.

Summary of main results obtained varying engine speed.

Table 5 presents the principal results obtained for operating parameters. From this table it can be concluded that at higher rpm volumetric efficiency is reduced due to increased air velocity in the intake system and increased pressure friction losses. Additionally, as there is less time to fill the cylinder, which results in lower air intake, the maximum pressure reached is reduced as well as the indicated work. Net power increases due to the influence of increasing rpm. Indicated efficiency tends to increase with increasing rpm, because burning the same fuel amount at higher rpm and net power is higher. Figure 6 shows that maximum pressure is obtained in all cases close to the TDC. One can also appreciate that the higher the rpm, the mass fraction burned during the premixed phase and the mean indicated pressure are reduced

4.4. Effect of compression ratio on CIE operating parameters

Now we will consider the effect of varying the compression ratio on the main engine operating parameters assuming valve inlet pressure equal to atmospheric pressure. Figure 7 shows the in-cylinder pressure variation and the accumulated burned mass fraction for various compression ratios and in Table 6 we can see the summary of main results obtained for different compression ratios.

Figure 7.

Pressure and accumulated burned fraction for different compression ratios

By analyzing Figure 7 and Table 6 we can conclude that: i) volumetric efficiency is unaffected by compression ratio since is not dependent on it, ii) at higher compression ratios there is a reduction of burned mass fraction in the premixed phase, iii) there is an indicated work increase because maximum pressure and temperature increase, resulting in increased efficiency and mean pressures, iv) pumping work increases when compression ratio increases as pressures are greatest during exhaust, v) cycle net work increases because net work and power increase at higher rates than pumping work.

Compression ratio182022
ηv0.720.720.72
β45.3342.6840.24
Wi [kJ]0.730.750.77
Wnet [kW]10.9911.3611.68
mip [kPa]1200.221233.521258.36
ηi0.450.460.47
Pmax [kPa]7134.937791.118349.81
Tmax [K]2400.712399.342392.72

Table 6.

Summary of main results obtained varying compression ratio.

4.5. Effect of atmospheric pressure on CIE operating parameters

It is well known that increasing the engine speed the start of combustion must anticipate in order that maximum pressures remain nearly constant and occur near TDC. For this reason in present study the combustion duration will remain constant and, its beginning will be advanced. Figure 8 shows the indicator diagram for various engine speeds. In Figure 8 which shows the indicator diagram for various engine speeds, it is observed that the maximum pressure in all studied cases occurs near TDC and its value is approximately constant. Since optimum injection angle for each speed is not known it is not easy to superimpose the curves which is what is desired in an engine in order to obtain maximum power and efficiency in a wide speeds range.

Figure 8.

Pressure and accumulated burned fraction for various atmospheric pressures

Atmospheric pressure100 kPa90 kPa80 kPa
ηv0.720.720.72
β45.2545.9346.59
Wi [kJ]0.740.650.58
Wnet [kW]11.139.908.80
mip [kPa]1217.031073.99948.11
ηi0.450.440.44
pmax [kPa]7231.676406.815669.70
Tmax [K]2401.912391.062380.09

Table 7.

Summary of main results obtained varying engine speed.

The summary of results in Table 7 shows that increasing the engine speed i) volumetric efficiency is reduced which is the expected behavior for CIE, ii) mass fraction burned during the premixed phase is reduced and iii) mean pressure decreases possibly due to a combustion beginning advance not enough for the increase in speed considered.

Advertisement

5. Conclusions

A methodology to theoretically determine the thermodynamic working cycle of direct injection compression ignition engines was presented.

A computer program based on the application of the above methodology was introduced.

Results of several studies performed with the program was presented and discussed. From some results we can conclude:

  • As the engine speed increases the start of combustion should anticipate maintaining nearly constant maximum cycle pressure as well as the angle at which this occurs.

  • At higher compression ratios, the greater the potency and mean effective pressure, this coincides with the theoretical internal combustion engines behavior.

  • At higher atmospheric pressure the engine volumetric efficiency will increase which ensures a greater air intake and thus an increase in produced power.

Notation

Aarea
cpiston stroke
Cddischarge coefficient
Cpconstant pressure specific heat
Cvconstant volume specific heat
d( )preceding term ordinary derivative
Dppiston diameter
dvvalve diameter
Eenergy
especific energy
ffuel air ratio
hspecific enthalpy
Henthalpy
hcconvection heat transfer coefficient
hcmean convection heat transfer coefficient
Hilower heating value
IDignition delay
jengine strokes
kconduction heat transfer coefficient
Lvvalve lift
mmass
m˙mass flow
m~adimensional mass
mipmean indicated pressure
NuNusselt number
ppressure
p~adimensional pressure
pdowndownstream pressure
pupupstream pressure
Qheat
Q~adimensional heat
Q˙heat flow
qheat per unit mass
Rgas constant
rccompression ratio
ReReynolds number
RLAconnecting rod length to crank radius ratio
rpmengine speed rpm
Ttemperature
T~adimensional temperature
ttime
uspecific internal energy
Uinternal energy
Vvolume
V~adimensional volume
vspecific volume
Vmppiston mean velocity
Wwork
W˙power
wspecific work
Xbburned mass fraction
yii gas fraction

Advertisement

Greek. letters

βpremixed phase burning factor
γpolytropic coefficient
∂( )preceding term partial derivative
ηefficiency
μviscosity
π3.1416
ρdensity
φcrankshaft rotation angle
ωangular velocity
ϕequivalence ratio (mixture richness)

Subindex

aair
admadmission
arreal air
attheoric air
bbblow by
cccombustion chamber
cilcylinder
comcombustion
comcompression
ddisplaced
difdiffusive phase
einput
EGRexhaust gas recirculated
exhexhaust
expexpansion
ffuel
grresidual gasses
iindicated
lostlost
maxmaximum
netnet
no combwithout combustion
Pproducts
prepremixed phase
pumpumping
Rreactants
refreference condition
soutput
stostoichiometric
Tottotal
vvolumetric
vccontrol volume
wallcylinder wall

References

  1. 1. JovajM. SAutomotive vehicle engines. Moscu: Editorial MIR; 1982
  2. 2. PayriFand DesantesJ. MMotores de combustión interna alternativos.Valencia, (Spain): Servicio de Publicaciones Universidad Politécnica de Valencia; 2011
  3. 3. RajputR. KInternal combustion engines. New Delhi: Laxmi Publications Ltd; 2007
  4. 4. HoltD. JThe diesel engine. Warrendale (PA): Society of Automotive Engineers; 2004
  5. 5. GanesanVInternal combustion engines. New Delhi: Tata McGraw-Hill Publishing Company Ltd; 2008
  6. 6. ArcoumanisCBicenA. FWhitelawJ. HEffect of inlet parameters on the flow characteristics in a four-stroke model engineWarrendale, (PA): Society of Automotive Engineers, 1982
  7. 7. BensonRThe thermodynamics and gas dinamics of internal combustión engines. Oxford, (England): Clarendon Press; 1986
  8. 8. DepcikCAssanisDA universal heat tranfer correlation for intake and exhaust flows in a spark ignition internal combustion engine.SAE Transaction 2002-01-0372, 1972
  9. 9. WoschniGComputer programs to determine the relationship, and thermal load in diesel engines. SAE Technical Paper 650450, 1965
  10. 10. AnnandW. JPinfoldDHeat transfer in the cylinder of a motored reciprocating engineSAE Technical Paper 800457, 1980
  11. 11. SitkeiGHeat Transfer and termal loading in internal combustión engine. Budapest: Akadémiai Kiadó; 1974
  12. 12. GosmanA. DComputer modeling of flow and heat transfer in engines, progress and prospects. In Proceedings of COMODIA 1985Tokio; 1985.
  13. 13. MeisnerSSorensonSComputer simulation of intake and exhaust manifold flow and heat transfer. SAE Technical Paper 860242, 1986
  14. 14. TaylorC. FThe internal-combustion engines in theory and practice: combustion, fuels, materials, design. Cambridge, (MA): MIT Press; 1982
  15. 15. HeywoodJ. BInternal combustion engine fundamentalsNew York: McGraw-Hill; 1988
  16. 16. AssanisD. NFilipiZ. SFivelandS. BSyrimisMA Predictive ignition delay correlation under steady-state and transient operation of a direct injection diesel engine. Eng Gas Turbines Power 2003
  17. 17. LiuHSimulation model for steady state and transient cold starting operation of diesel enginesETD Collection for Wayne State University. Paper AAI3037106, Michigan, 2001
  18. 18. WoschniGA universally applicable equation for the instantaneous heat transfer coefficient in the internal combustion engineSAE Technical Paper 670931, 1967
  19. 19. HardenbergH. OHaseFAn empirical formula for computing the pressure rise delay of a fuel from its cetane number and from the relevant parameters of direct-injection diesel engines.SAE Technical Paper 790493, 1979
  20. 20. GardnerTHeneinNDiesel starting: a mathematical model. SAE Technical Paper 880426, 1988
  21. 21. WayR. JPrograms for determination of composition and thermodynamic properties of combustion products for internal combustion engines calculations. Proceedings of the Inst Mech Engrs 1977
  22. 22. DesantesJ. MLapuertaMFundamentos de combustión. Valencia, (Spain): Servicio de Publicaciones Universidad Politécnica de Valencia; 1991
  23. 23. El-MahallawyFEl-DinH. SFundamentals and technology of combustionLondon: Elsevier Science Ltd; 2002
  24. 24. JarquinGPolupanGRodríguezG. JCálculo de los productos de combustión empleando métodos numéricos. Mecánica Computacional 20032224422452
  25. 25. McbrideB. JGordonSFortran IV program for calculation of thermodynamic dataNASA Publication TN-D-4097, 1967
  26. 26. GordonSMcbrideB. JComputer program for calculation of complex chemical equilibrium compositions, rocket performance, incident and reflected shocks, and Chapman-Jouguet detonations. NASA Publication SP 273, 1971
  27. 27. HarkerJ. HThe Calculation of equilibrium flame gas composition. Journal Inst Fuel 196740206210
  28. 28. HarkerJ. HAllenD. AThe calculation of the temperature and composition of flame gases. Journal Inst Fuel 196942183187
  29. 29. OlikaraCBormanG. LA computer program for calculating properties of equilibrium combustion products with some applications to IC. engines. SAE Technical Paper Nº 750468, 1975
  30. 30. AgrawalD. DSharmaS. PGuptaC. PThe calculation of temperature and pressure of flame gases following constant volume combustionJournal Inst Fuel 197750121124
  31. 31. AgrawalDGuptaC. PComputer program for constant pressure or constant volume combustion calculations in hydrocarbon-air ssystems. Transactions of the ASME 19771977246254
  32. 32. LapuertaMArmasOHernándezJDiagnosis of DI diesel combustión from in-cylinder pressure signal by estimation of mean thernodynamic properties of the gas. Applied Thermal Engineering, 1997
  33. 33. AraqueJFygueroaSMartínMModelado de la combustión en un MECH-CFR. In Memorias del IV Congreso Nacional de Ingeniería Mecánica, Mérida (Venezuela), 2001
  34. 34. FergusonCKirkpatrickAInternal combustion: applied thermosciences. New York: John Wiley and Son; 2001
  35. 35. BensonRAnnandWBaruahPA simulation model including intake and exhaust systems for a single cylinder four stroke cycle spark ignition engineInternational Journal of Mechanical Sciences197517297124
  36. 36. WoodwardJAir standard modelling for closed cycle diesel engines. Proc Instn Mech Engrs 1995
  37. 37. HullT. EEnrightW. HJacksonK. RRunge-Kutta research at TorontoJournal Applied Numerical Mathematics 199622225236
  38. 38. AssanisD. NValve event optimization in a spark-ignition engineJ Eng Gas Turbines Power 19901123341348
  39. 39. WatsonNPilleyA. DMarzoukMA Combustion correlation for diesel engine simulationSAE Technical Paper 800029, 1980
  40. 40. MiyamotoTHayashiKHaradaASasakiSAkagawaHTsujimuraKNumerical simulation of premixed lean diesel combustion in a DI engine. In Proceedings of COMODIA 98, Kioto, 1998
  41. 41. AnnandW. JHeat transfer in the cylinders of reciprocating internal ccombustion engines. Proceedings of the Institution of Mechanical Engineers 19631771973996
  42. 42. ChapraSCanaleRMétodos numéricos para ingenieros. México: Mc Graw Hil; 2002

Written By

Simón Fygueroa, Carlos Villamar and Olga Fygueroa

Submitted: 21 November 2011 Published: 14 November 2012