InTechOpen uses cookies to offer you the best online experience. By continuing to use our site, you agree to our Privacy Policy.

Engineering » Mechanical Engineering » "Internal Combustion Engines", book edited by Kazimierz Lejda and Pawel Wos, ISBN 978-953-51-0856-6, Published: November 14, 2012 under CC BY 3.0 license. © The Author(s).

Chapter 4

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

By Simón Fygueroa, Carlos Villamar and Olga Fygueroa
DOI: 10.5772/50028

Article top


Engine control volume
Figure 1. Engine control volume
Block diagram and routines used by MECID computer program
Figure 2. Block diagram and routines used by MECID computer program
p vs φ and T vs φ diagrams for various combustion beginning angles
Figure 3. p vs φ and T vs φ diagrams for various combustion beginning angles
Accumulated and instant burned mass fraction for various combustion beginning angles
Figure 4. Accumulated and instant burned mass fraction for various combustion beginning angles
Pressure and accumulated burned fraction for different combustion duration
Figure 5. Pressure and accumulated burned fraction for different combustion duration
In-cylinder pressure for various engine speeds
Figure 6. In-cylinder pressure for various engine speeds
Pressure and accumulated burned fraction for different compression ratios
Figure 7. Pressure and accumulated burned fraction for different compression ratios
Pressure and accumulated burned fraction for various atmospheric pressures
Figure 8. Pressure and accumulated burned fraction for various atmospheric pressures

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

Simón Fygueroa1, 2, 3, 4, 5, Carlos Villamar6 and Olga Fygueroa3, 6, 7

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.

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:


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:


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

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


From equivalence ratio (mixture richness) definition:


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


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:


Since the work due to a volume change is:

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


or in terms of enthalpy:


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


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


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


From Eq. 14:


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


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


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


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


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


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:


and defining:


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

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


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


Now, considering:

And differentiating:


Differentiating ideal gas equation and solving for dRdt


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


Replacing Ec. 32 in Ec. 27:




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




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



Time differentiating and solving for dρdt:


Replacing Ec. 39 in Ec. 37:


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:




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]:


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

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


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

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


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.

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:


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:


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:

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 φ:


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


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:


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:


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,


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:


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:


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:


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


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


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


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


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

ExpresiónValores sugeridos

Table 2.

Constants for Heywood and Watson expressions

mdif1.5 ID o 1.9 ID

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:


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


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


Blow by mass can be found from mass conservation:



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




and defining the following dimensionless variables:


dimensionless heat losses for unit crankshaft rotation angle are:


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


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:


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:

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


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]:


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]:


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


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:


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:


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:


The theoretical incoming air mass is calculated by:

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

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:


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:


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

H22H K1=y42y2(pp0)
O22O K2=y42y9(pp0)
H2+O22OH K3=y32y2y9
H2+12O2H2O K5=y1y2y912(pp0)12
CO+12O2CO2 K6=y7y8y812(pp0)12

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:

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:


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:


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:


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


and Eqs. 114 to 116 are transformed into:


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.


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 + %EGRProducts

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]:


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


Net power is obtained through the following expression:


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


Efficiencies are determined with:


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

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º
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
Wi [kJ]1.421.181.03
Wnet [kW]14.6316.0117.28
mip [kPa]2314.471922.251678.68
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
Wi [kJ]0.730.750.77
Wnet [kW]10.9911.3611.68
mip [kPa]1200.221233.521258.36
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
Wi [kJ]0.740.650.58
Wnet [kW]11.139.908.80
mip [kPa]1217.031073.99948.11
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.

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.


cpiston stroke
Cddischarge coefficient
Cpconstant pressure specific heat
Cvconstant volume specific heat
d( )preceding term ordinary derivative
Dppiston diameter
dvvalve diameter
especific energy
ffuel air ratio
hspecific enthalpy
hcconvection heat transfer coefficient
hcmean convection heat transfer coefficient
Hilower heating value
IDignition delay
jengine strokes
kconduction heat transfer coefficient
Lvvalve lift
m˙mass flow
m~adimensional mass
mipmean indicated pressure
NuNusselt number
p~adimensional pressure
pdowndownstream pressure
pupupstream pressure
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
T~adimensional temperature
uspecific internal energy
Uinternal energy
V~adimensional volume
vspecific volume
Vmppiston mean velocity
wspecific work
Xbburned mass fraction
yii gas fraction

Greek. letters

βpremixed phase burning factor
γpolytropic coefficient
∂( )preceding term partial derivative
φcrankshaft rotation angle
ωangular velocity
ϕequivalence ratio (mixture richness)


arreal air
attheoric air
bbblow by
cccombustion chamber
difdiffusive phase
EGRexhaust gas recirculated
grresidual gasses
no combwithout combustion
prepremixed phase
refreference condition
vccontrol volume
wallcylinder wall


1 - M. S Jovaj, Automotive vehicle engines. Moscu: Editorial MIR; 1982
2 - F Payri, and J. M Desantes, Motores de combustión interna alternativos.Valencia, (Spain): Servicio de Publicaciones Universidad Politécnica de Valencia; 2011
3 - R. K Rajput, Internal combustion engines. New Delhi: Laxmi Publications Ltd; 2007
4 - D. J Holt, The diesel engine. Warrendale (PA): Society of Automotive Engineers; 2004
5 - V Ganesan, Internal combustion engines. New Delhi: Tata McGraw-Hill Publishing Company Ltd; 2008
6 - C Arcoumanis, A. F Bicen, J. H Whitelaw, Effect of inlet parameters on the flow characteristics in a four-stroke model engineWarrendale, (PA): Society of Automotive Engineers, 1982
7 - R Benson, The thermodynamics and gas dinamics of internal combustión engines. Oxford, (England): Clarendon Press; 1986
8 - C Depcik, D Assanis, A universal heat tranfer correlation for intake and exhaust flows in a spark ignition internal combustion engine.SAE Transaction 2002-01-0372, 1972
9 - G Woschni, Computer programs to determine the relationship, and thermal load in diesel engines. SAE Technical Paper 650450, 1965
10 - W. J Annand, D Pinfold, Heat transfer in the cylinder of a motored reciprocating engineSAE Technical Paper 800457, 1980
11 - G Sitkei, Heat Transfer and termal loading in internal combustión engine. Budapest: Akadémiai Kiadó; 1974
12 - A. D Gosman, Computer modeling of flow and heat transfer in engines, progress and prospects. In Proceedings of COMODIA 1985Tokio; 1985.
13 - S Meisner, S Sorenson, Computer simulation of intake and exhaust manifold flow and heat transfer. SAE Technical Paper 860242, 1986
14 - C. F Taylor, The internal-combustion engines in theory and practice: combustion, fuels, materials, design. Cambridge, (MA): MIT Press; 1982
15 - J. B Heywood, Internal combustion engine fundamentalsNew York: McGraw-Hill; 1988
16 - D. N Assanis, Z. S Filipi, S. B Fiveland, M Syrimis, A Predictive ignition delay correlation under steady-state and transient operation of a direct injection diesel engine. Eng Gas Turbines Power 2003
17 - H Liu, Simulation model for steady state and transient cold starting operation of diesel enginesETD Collection for Wayne State University. Paper AAI3037106, Michigan, 2001
18 - G Woschni, A universally applicable equation for the instantaneous heat transfer coefficient in the internal combustion engineSAE Technical Paper 670931, 1967
19 - H. O Hardenberg, F Hase, An 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 - T Gardner, N Henein, Diesel starting: a mathematical model. SAE Technical Paper 880426, 1988
21 - R. J Way, Programs for determination of composition and thermodynamic properties of combustion products for internal combustion engines calculations. Proceedings of the Inst Mech Engrs 1977
22 - J. M Desantes, M Lapuerta, Fundamentos de combustión. Valencia, (Spain): Servicio de Publicaciones Universidad Politécnica de Valencia; 1991
23 - El-F Mahallawy, El-H. S Din, Fundamentals and technology of combustionLondon: Elsevier Science Ltd; 2002
24 - G Jarquin, G Polupan, G. J Rodríguez, Cálculo de los productos de combustión empleando métodos numéricos. Mecánica Computacional 20032224422452
25 - B. J Mcbride, S Gordon, Fortran IV program for calculation of thermodynamic dataNASA Publication TN-D-4097, 1967
26 - S Gordon, B. J Mcbride, Computer program for calculation of complex chemical equilibrium compositions, rocket performance, incident and reflected shocks, and Chapman-Jouguet detonations. NASA Publication SP 273, 1971
27 - J. H Harker, The Calculation of equilibrium flame gas composition. Journal Inst Fuel 196740206210
28 - J. H Harker, D. A Allen, The calculation of the temperature and composition of flame gases. Journal Inst Fuel 196942183187
29 - C Olikara, G. L Borman, A computer program for calculating properties of equilibrium combustion products with some applications to IC. engines. SAE Technical Paper Nº 750468, 1975
30 - D. D Agrawal, S. P Sharma, C. P Gupta, The calculation of temperature and pressure of flame gases following constant volume combustionJournal Inst Fuel 197750121124
31 - D Agrawal, C. P Gupta, Computer program for constant pressure or constant volume combustion calculations in hydrocarbon-air ssystems. Transactions of the ASME 19771977246254
32 - M Lapuerta, O Armas, J Hernández, Diagnosis of DI diesel combustión from in-cylinder pressure signal by estimation of mean thernodynamic properties of the gas. Applied Thermal Engineering, 1997
33 - J Araque, S Fygueroa, M Martín, Modelado de la combustión en un MECH-CFR. In Memorias del IV Congreso Nacional de Ingeniería Mecánica, Mérida (Venezuela), 2001
34 - C Ferguson, A Kirkpatrick, Internal combustion: applied thermosciences. New York: John Wiley and Son; 2001
35 - R Benson, W Annand, P Baruah, A simulation model including intake and exhaust systems for a single cylinder four stroke cycle spark ignition engineInternational Journal of Mechanical Sciences197517297124
36 - J Woodward, Air standard modelling for closed cycle diesel engines. Proc Instn Mech Engrs 1995
37 - T. E Hull, W. H Enright, K. R Jackson, Runge-Kutta research at TorontoJournal Applied Numerical Mathematics 199622225236
38 - D. N Assanis, Valve event optimization in a spark-ignition engineJ Eng Gas Turbines Power 19901123341348
39 - N Watson, A. D Pilley, M Marzouk, A Combustion correlation for diesel engine simulationSAE Technical Paper 800029, 1980
40 - T Miyamoto, K Hayashi, A Harada, S Sasaki, H Akagawa, K Tsujimura, Numerical simulation of premixed lean diesel combustion in a DI engine. In Proceedings of COMODIA 98, Kioto, 1998
41 - W. J Annand, Heat transfer in the cylinders of reciprocating internal ccombustion engines. Proceedings of the Institution of Mechanical Engineers 19631771973996
42 - S Chapra, R Canale, Métodos numéricos para ingenieros. México: Mc Graw Hil; 2002