Empirical values for burning factor
1. Introduction
Currently, one of the worldwide most used
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
The
During the
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
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
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
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
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.
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.
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
Solving Ec. 15 for
Replacing ideal gas equation in Ec. 19 and solving for
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
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
Replacing Ec. 27 in Ec. 20 and solving for
Now, considering:
And differentiating:
Differentiating ideal gas equation and solving for
Replacing Ec. 30 in Eq. 31 and solving for
Replacing Ec. 32 in Ec. 27:
Defining:
Collecting terms containing
Defining:
Replacing Ec. 36 in Eq. 35 and solving for
Since:
Time differentiating and solving for
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:
and:
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 terms
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:
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:
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]:
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.
Value | Hardenberg model | Watson model | Heywood model |
a1 | 0.746 | 0.926 | 0.80 - 0.95 |
b1 | 0.35 | 0.37 | 0.25 - 0.45 |
cc1 | 0.35 | 0.26 | 0.25 – 0.50 |
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.
Constant | Heywood | Watson | |
Expresión | Valores sugeridos | ||
C1 | 3 | ||
C2 | 5000 | 5000 | 5000 |
C3 | 6.908 | ||
C4 | 1.4 |
Constant | Miyamoto |
mpre | 4 |
mdif | 1.5 ID o 1.9 ID |
a | 0.69 |
φpre | +7° |
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
Blow by mass can be found from mass conservation:
Defining:
Net heat entering the system, considering that a part is lost by blow by, is:
with:
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 for
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]:
where
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:
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:
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:
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.
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 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.
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.
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 duration | 30º | 45º | 60º |
ηv | 0.72 | 0.72 | 0.72 |
Wi [kJ] | 0.60 | 0.73 | 0.75 |
Wnet [kJ] | 0.57 | 0.65 | 0.64 |
9.56 | 10.98 | 10.72 | |
mip [kPa] | 994.75 | 1200.21 | 1232.88 |
ηi net | 0.35 | 0.40 | 0.39 |
Pmax [kPa] | 7689.17 | 7134.93 | 6614.27 |
Tmax [K] | 2235.64 | 2400.71 | 2496.92 |
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.
Engine speed | 1500 rpm | 2000 rpm | 2500 rpm |
ηv | 0.92 | 0.69 | 0.55 |
β | 50.57 | 45.33 | 40.89 |
Wi [kJ] | 1.42 | 1.18 | 1.03 |
14.63 | 16.01 | 17.28 | |
mip [kPa] | 2314.47 | 1922.25 | 1678.68 |
ηi | 0.71 | 0.74 | 0.77 |
pmax [kPa] | 6703.23 | 6034.92 | 5506.50 |
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.
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 ratio | 18 | 20 | 22 |
ηv | 0.72 | 0.72 | 0.72 |
β | 45.33 | 42.68 | 40.24 |
Wi [kJ] | 0.73 | 0.75 | 0.77 |
10.99 | 11.36 | 11.68 | |
mip [kPa] | 1200.22 | 1233.52 | 1258.36 |
ηi | 0.45 | 0.46 | 0.47 |
Pmax [kPa] | 7134.93 | 7791.11 | 8349.81 |
Tmax [K] | 2400.71 | 2399.34 | 2392.72 |
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.
Atmospheric pressure | 100 kPa | 90 kPa | 80 kPa |
ηv | 0.72 | 0.72 | 0.72 |
β | 45.25 | 45.93 | 46.59 |
Wi [kJ] | 0.74 | 0.65 | 0.58 |
11.13 | 9.90 | 8.80 | |
mip [kPa] | 1217.03 | 1073.99 | 948.11 |
ηi | 0.45 | 0.44 | 0.44 |
pmax [kPa] | 7231.67 | 6406.81 | 5669.70 |
Tmax [K] | 2401.91 | 2391.06 | 2380.09 |
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.
Notation
A | area |
c | piston stroke |
Cd | discharge coefficient |
Cp | constant pressure specific heat |
Cv | constant volume specific heat |
d( ) | preceding term ordinary derivative |
Dp | piston diameter |
dv | valve diameter |
E | energy |
e | specific energy |
f | fuel air ratio |
h | specific enthalpy |
H | enthalpy |
hc | convection heat transfer coefficient |
mean convection heat transfer coefficient | |
Hi | lower heating value |
ID | ignition delay |
j | engine strokes |
k | conduction heat transfer coefficient |
Lv | valve lift |
m | mass |
mass flow | |
adimensional mass | |
mip | mean indicated pressure |
Nu | Nusselt number |
p | pressure |
adimensional pressure | |
pdown | downstream pressure |
pup | upstream pressure |
Q | heat |
adimensional heat | |
heat flow | |
q | heat per unit mass |
R | gas constant |
rc | compression ratio |
Re | Reynolds number |
RLA | connecting rod length to crank radius ratio |
rpm | engine speed rpm |
T | temperature |
adimensional temperature | |
t | time |
u | specific internal energy |
U | internal energy |
V | volume |
adimensional volume | |
v | specific volume |
Vmp | piston mean velocity |
W | work |
power | |
w | specific work |
Xb | burned mass fraction |
yi | i gas fraction |
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
a | air |
adm | admission |
ar | real air |
at | theoric air |
bb | blow by |
cc | combustion chamber |
cil | cylinder |
com | combustion |
com | compression |
d | displaced |
dif | diffusive phase |
e | input |
EGR | exhaust gas recirculated |
exh | exhaust |
exp | expansion |
f | fuel |
gr | residual gasses |
i | indicated |
lost | lost |
max | maximum |
net | net |
no comb | without combustion |
P | products |
pre | premixed phase |
pum | pumping |
R | reactants |
ref | reference condition |
s | output |
sto | stoichiometric |
Tot | total |
v | volumetric |
vc | control volume |
wall | cylinder wall |
References
- 1.
Automotive vehicle engines. Moscu: Editorial MIR;Jovaj M. S 1982 - 2.
andPayri F Desantes J. M Motores de combustión interna alternativos. Valencia, (Spain): Servicio de Publicaciones Universidad Politécnica de Valencia;2011 - 3.
Internal combustion engines. New Delhi: Laxmi Publications Ltd;Rajput R. K 2007 - 4.
The diesel engine. Warrendale (PA): Society of Automotive Engineers;Holt D. J 2004 - 5.
Internal combustion engines. New Delhi: Tata McGraw-Hill Publishing Company Ltd;Ganesan V 2008 - 6.
Arcoumanis C Bicen A. F Whitelaw J. H Effect of inlet parameters on the flow characteristics in a four-stroke model engine Warrendale, (PA): Society of Automotive Engineers,1982 - 7.
The thermodynamics and gas dinamics of internal combustión engines. Oxford, (England): Clarendon Press;Benson R 1986 - 8.
Depcik C Assanis D A universal heat tranfer correlation for intake and exhaust flows in a spark ignition internal combustion engine. SAE Transaction 2002-01-0372,1972 - 9.
Woschni G Computer programs to determine the relationship, and thermal load in diesel engines. SAE Technical Paper 650450, 1965 - 10.
Annand W. J Pinfold D Heat transfer in the cylinder of a motored reciprocating engine SAE Technical Paper 800457,1980 - 11.
Heat Transfer and termal loading in internal combustión engine. Budapest: Akadémiai Kiadó;Sitkei G 1974 - 12.
Computer modeling of flow and heat transfer in engines, progress and prospects. In Proceedings of COMODIAGosman A. D 1985 Tokio; 1985. - 13.
Computer simulation of intake and exhaust manifold flow and heat transfer. SAE Technical Paper 860242,Meisner S Sorenson S 1986 - 14.
Taylor C. F The internal-combustion engines in theory and practic e: combustion, fuels, materials, design. Cambridge, (MA): MIT Press;1982 - 15.
Heywood J. B Internal combustion engine fundamentals New York: McGraw-Hill;1988 - 16.
A Predictive ignition delay correlation under steady-state and transient operation of a direct injection diesel engine. Eng Gas Turbines PowerAssanis D. N Filipi Z. S Fiveland S. B Syrimis M 2003 - 17.
Liu H Simulation model for steady state and transient cold starting operation of diesel engines ETD Collection for Wayne State University. Paper AAI3037106, Michigan,2001 - 18.
Woschni G A universally applicable equation for the instantaneous heat transfer coefficient in the internal combustion engine SAE Technical Paper 670931,1967 - 19.
Hardenberg H. O Hase F 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.
Diesel starting: a mathematical model. SAE Technical Paper 880426,Gardner T Henein N 1988 - 21.
Programs for determination of composition and thermodynamic properties of combustion products for internal combustion engines calculations. Proceedings of the Inst Mech EngrsWay R. J 1977 - 22.
Fundamentos de combustión. Valencia, (Spain): Servicio de Publicaciones Universidad Politécnica de Valencia;Desantes J. M Lapuerta M 1991 - 23.
El- El-Mahallawy F Din H. S Fundamentals and technology of combustion London: Elsevier Science Ltd;2002 - 24.
Cálculo de los productos de combustión empleando métodos numéricos. Mecánica ComputacionalJarquin G Polupan G Rodríguez G. J 2003 22 2442 2452 - 25.
Mcbride B. J Gordon S Fortran IV program for calculation of thermodynamic data NASA Publication TN-D-4097,1967 - 26.
Gordon S Mcbride B. J Computer program for calculation of complex chemical equilibrium compositions, rocket performance, incident and reflected shocks, a nd Chapman-Jouguet detonations. NASA Publication SP 273,1971 - 27.
The Calculation of equilibrium flame gas composition. Journal Inst FuelHarker J. H 1967 40 206 210 - 28.
The calculation of the temperature and composition of flame gases. Journal Inst FuelHarker J. H Allen D. A 1969 42 183 187 - 29.
A computer program for calculating properties of equilibrium combustion products with some applications to IC. engines. SAE Technical Paper Nº 750468,Olikara C Borman G. L 1975 - 30.
TheAgrawal D. D Sharma S. P Gupta C. P calculation of temperature and pressure of flame gases following constant volume combustion ournal Inst Fuel1977 50 121 124 - 31.
Computer program for constant pressure or constant volume combustion calculations in hydrocarbon-air ssystems. Transactions of the ASMEAgrawal D Gupta C. P 1977 1977 246 254 - 32.
Diagnosis of DI diesel combustión from in-cylinder pressure signal by estimation of mean thernodynamic properties of the gas. Applied Thermal Engineering,Lapuerta M Armas O Hernández J 1997 - 33.
Modelado de la combustión en un MECH-CFR. In Memorias del IV Congreso Nacional de Ingeniería Mecánica, Mérida (Venezuela),Araque J Fygueroa S Martín M 2001 - 34.
Internal combustion: applied thermosciences. New York: John Wiley and Son;Ferguson C Kirkpatrick A 2001 - 35.
Benson R Annand W Baruah P A simulation model including intake and exhaust systems for a single cylinder four stroke cycle spark ignition engine 1975 17 2 97 124 - 36.
Air standard modelling for closed cycle diesel engines. Proc Instn Mech EngrsWoodward J 1995 - 37.
Hull T. E Enright W. H Jackson K. R Runge-Kutta research at Toronto Journal Applied Numerical Mathematics1996 22 225 236 - 38.
Assanis D. N Valve event optimization in a spark-ignition engine J Eng Gas Turbines Power1990 112 3 341 348 - 39.
Watson N Pilley A. D Marzouk M A Combustion correlation for diesel engine simulation SAE Technical Paper 800029,1980 - 40.
Numerical simulation of premixed lean diesel combustion in a DI engine. In Proceedings of COMODIA 98, Kioto,Miyamoto T Hayashi K Harada A Sasaki S Akagawa H Tsujimura K 1998 - 41.
Heat transfer in the cylinders of reciprocating internal ccombustion engines. Proceedings of the Institution of Mechanical EngineersAnnand W. J 1963 177 1 973 996 - 42.
Métodos numéricos para ingenieros. México: Mc Graw Hil;Chapra S Canale R 2002