Printer Friendly

The modelling of operation of the compression ignition engine powered with diesel fuel with LPG admixture.

1. Introduction

The gaseous fuels are used in transportation due to their lower price However, the car manufacturers rarely offer engine units powered with this type of fuel. That is why the row of research courses is run aiming the adaptation different sorts of engines for powering with gaseous fuels, what is connected with the cost increase, thus longer investment return [1-3]. Application the engines to operate on gaseous fuel is related to many problems [4, 5]. They are mainly: adaptation costs, refilling, engine durability, fuel tank placing, service and operation safety. As fuels capable of transportation usage can be named: natural gas (CNG--compressed, LNG--liquefied), liquid mixture of propane and butane (LPG) and biogas. Basic parameters of these fuels are listed in Table 1.

In the case of petrol engines, the gaseous fuel can completely replace petrol fuel due to external ignition. In diesel engines, the diesel fuel is essential to initiate the ignition (dual-fueled engine) or applying the spark plug, then the engine can be powered only with gaseous fuel [6, 7]. Due to the fact, that the diesel engine operates at the air excess, the delivering the gaseous fuel is possible in cylinder filling process at the expense of reduced air mass, containing oxygen, essential in combustion process. In this case the engine uses advantages of both fuels and is capable of operating with minimal interference in engine construction. The second case requires higher range of interfering with engine infrastructure, where the additional installation of subsystems responsible for ignition initialization (spark plug) and reducing the compression ratio are required.

That is why the attempt was made to create the simulation model of the diesel engine powered with admixture of LPG of different degrees. Modern, commercial solutions of LPG powering systems dedicated for diesel engines are meeting several operational problems. That is why the calculations were needed, basing on which, the conclusion about source of some engine issues was possible to draw.

Internal combustion engines used in road transportation are mostly compression ignition engines powered with diesel fuel. The delay of injection start, which initiates the ignition of dual-fuel engine (diesel/LPG) effects in increasing the efficiency and decreasing the hard work effect [8]. On the other hand, the knocking combustion occurs resulting in lowered durability, including failure wear [9]. The effect of the knock can be lowered by means of proper selection of diesel/LPG mixing ratio, lowering by this hydrocarbons and carbon monoxide emissions [10]. Increasing the preliminary dosage mass results in thermal efficiency increase, simultaneously increasing the possibility of knock, thus should it not be used in dual-fuel powering at maximal loads [11]. The negative influence of dualfuel powering on CO emission can be reduced by increasing the preliminary dosage as a result of lengthening the injection, however influences this the NOx emission increase [12].

Unfortunately, the dominant feature with LPG powering is the variability of chemical composition and pollution content, what influences on the diversity of toxic compounds emission to the atmosphere. In the paper [13] it was shown, that at proper LPG composition, the engine works smoothly with up to 90 % ON dosage replaced by LPG. In general, engines powered with gaseous fuel are characterized by lower toxic compounds emission, as it is in the case of petrol or diesel engines, however in the results of some technical imperfections, the total emission increases compared to base fuel powering [14].

The amount of heat transferred with the exhaust gases to the diesel engine exhaust system is lower for the dual-fuel solution. Faster gas combustion in smaller combustion chamber volume causes the increase of maximum temperature and cycle pressure [15].

On the other hand, the CNG injection to the combustion chamber results in higher mechanical load as an effect of maximum pressure increase, comparing to CNG brought via intake manifold, where as bringing CNG via intake manifold results in increased thermal loads [16].

2. The characteristics of research object and necessary conditions from simulation initialization

The object of research is the compression ignition engine with parameters listed in Table 2. The simulation research concerns the engine powered only by diesel fuel and powered by mixture of diesel and LPG in gaseous phase injected indirectly to the inlet duct during intake stroke.

For dual-fuel powering the system presented in Fig. 1 was applied, where the nozzle was inserted to support gaseous fuel to the intake manifold. The engine in intake stroke is partially powered with LPG, making the combustible mixture. The injection of diesel fuel at the final phase of compression stroke is responsible for ignition initiation.

In the next step of research, the following assumptions were applied:

a) constant wall temperature: 360°K;

b) density of flowing gas is pressure and temperature dependent;

c) air with following properties:

--density in normal conditions--1.205 kg/[m.sup.3],

--heating value--1.005 kJ/(kg °K),

--thermal conductivity--0.0259 W/(m °K),

--dynamic viscosity--1.81e-05 kg/ms,

--kinematic viscosity--1.506e-05 [m.sup.2]/s,

--molecular weight--28.96 kg/kmol;

d) gravity omitted;

e) turbulence model k-epsilon.

The simulation was conducted by using the transient option, which means that the parameters of medium flow are variable in time. The boundary conditions of the engine operation were applied as for driving with high loads.

Thus assumed:

--engine speed--2000 rpm,

--intake pressure--0.1 MPa,

--back pressure--0.1 MPa.

For simulation initialization purposes, it was necessary to determine the valves timing and put it into the software workspace (Fig. 2).

As the initial conditions for the whole fluid computational space the constant pressure at value of 100 kPa, temperature at value of 290°K and percentage composition of air with 21% content of oxygen, water steam of 0.5%, carbon dioxide of 0.03% were applied.

The diesel fuel injection was assumed, in which the stream is cone-shaped with the apex angle of 120 deg, drop diameter of 40 µm, the mass flow of 4.5 gm/s and the stream velocity equal to 250 m/s. The LPG injection war realized by assuming the constant injection pressure of 70 kPa and injector nozzle diameter of 2 mm. The LPG composition was assumed as 50 % of propane and 50% of butane. Injection times for specific variants are shown in Table 3.

3. The simulation research methodology

For simulation research purpose the solid model was created, which corresponds to research object geometry, in which one full power cycle is possible to perform [17]. It consists of the volume, which is inside intake duct, combustion chamber and exhaust duct, this is the space in which flow processes occur (Fig. 3).

As an environment for model research the software was chosen, which bases on computational fluid dynamics (CFD) and uses the finite volume elements method --ANSYS Fluent. The model changes during simulation as a result of piston and valves displacement according to crankshaft movement and timing. The change of the computational space is possible due to dynamic mesh application. The simulation begins from the TDC of the piston, which corresponds to the angular crankshaft position at 0° value.

The finite volume mesh is built of about 300000 elements in TDC position and about 500000 in BDC position. The quality of the mesh is defined by means of statistical tool called Skewness, available in ANSYS Mesh module, of which the maximum value was 0.869 and the average 0.24.

Skewness has an influence on the mesh structure during rebuilding at consecutive time steps, the value of which stays on the assumed level, i.e. in this case 0.6. Skewness contained in proper range determines high quality of given results. For 3D models permissible values range from 0 to 0.95. The ways of computing the parameter are shown in Figs. 4 and 5 and are represented accordingly by dependencies (1) and (2) for tetrahedronical cells and cells of other shapes.

Skewness = optimal cell size - cell size/ optimal cell size, (1)

Skewness = MAX [[?.sub.max] - [?.sub.e]/180-[?.sub.e], [?.sub.e] - [?.sub.min]/ [?.sub.e]], (2)

The model has multisolid structure, it consists of ten separate solids. Multisolidness of the model is determined by ANSYS Fluent software requirements regarding variable geometry during simulation. The simplifications were not applied in the model, in the form of section of the elements or of symmetrical volume, which significantly shorten the calculation time, but are able to be applied in the case of predictible steady state research. In the analyzed case they are not applicable due to the turbulent flow character, strong air swirl during load exchange and uneven flame spread during combustion.

The model consists mainly of tetrahedral elements, which occur in solids marked in Fig. 6 by numbers 1, 2, 4, 7, 9. It is needed to say that solid 9 is subjected to rebuild during the whole simulation, the other are constant. In solids 3, 5, 6, 8, 10 the hexaheddral elements occur with prismatic shape with quadrangular or triangular base, which are used accordingly to software requirements for getting in these areas cell extension or division, thus creation of new layers during movable elements movement.

This solution allows valve outline to displace without necessity of rebuilding solids 1 and 2, which would be related to lengthening the calculation time. Between consecutive solids the interfaces are created letting maintain the continuity of the flow between displacing areas of various structure, or creating the wall out of the interface denying further flow in the time range corresponding to the valves closure.

For simulating processes occuring inside the compression ignition engine the models were used, which are supplied by ANSYS Fluent software. The tool was used allowing initiating the combustion process as a result of compression ignition. It takes into account the amount of time passing from injecting the fuel to self-ignition occurence, i.e. the ignition lag, which is dependent on the parameters due to following formula [19]:


where [C.sub.1] is pro-exponential; CN is engine speed; [[bar.S].sub.p] is mean piston speed; [E.sub.a] is effective activation energy; R is universal gas constant; T is temperature; p is absolute pressure; [e.sub.p] is pressure exponent.

[E.sub.a] = [E.sub.hh]/CN + 25, (4)

where [E.sub.hh] is activation energy; CN is cetane number.

The combustion process was realized by using Eddy-dissipation model for turbulent models. The fuels are oxidated according to the reactions [20]:

2[C.sub.10] [H.sub.22] +31[O.sub.2] =20C[O.sub.2] +22[H.sub.2]O ; (5)

[C.sub.3][H.sub.8] + 5[O.sub.2] = 3C[O.sub.2] + 4[H.sub.2]O; (6)

2[C.sub.4][H.sub.10] + 13[O.sub.2] = 8C[O.sub.2] = 10[H.sub.2] O. (7)

4. Research results and discussion

The speed distribution during cylinder filling (Fig. 7--the case when cylinder is filled only with air) showed the correctness of dividing the fluid to components depending on prognosed work parameters gradients. The swirls can be observed in almost whole analyzed space (stream lines with velocity values). Maximum velocity values were noted at level of 17 m/s in valve slot necking.

The average pressure and temperature courses show differences in the case of applying LPG admixture (Fig. 8 and 9). After occurence of ignition the pressure reaches greater gradient in the case of diesel+LPG, maximum value is higher by about 0.25 MPa as well. This implies that the engine powered by LPG admixture is capable of reaching higher values of working indicators, i.e. torque.

The increase of temperature by about 50 °K is related to that (Fig. 9).

That is why it requires monitoring the temperature in exhaust systems of engine of such type, what is practiced in commercial LPG system solutions.

Static pressure and temperature distributions at TDC for Simulation 1 (diesel) are presented in Fig. 10. The biggest pressure values occur in lower part of combustion chamber what is the result of proper directing of fuel streams. (Fig. 10, a and d). The temperature reaches the biggest values at the edges of the stream as well, according to the flame spread direction. The diesel fuel injection to the space filled with air admixtured by LPG show increased pressure in the area of the top part of combustion chamber, the similar thing occurs in the case of temperature (Simulation 2--diesel+LPG--Fig. 10, b and e). Simulation 3 (diesel initial injection + LPG) shows the lowest pressure and temperature values, what causes the difficulties of reaching the indicators values at a level comparable with diesel powering.

The maximum value of reaction heat (Fig. 11) in beginning combustion phase in the case of Simulation 2 is significantly lower than in diesel (Simulation 1), the contrary situation occurs in the case of Simulation 3, where lowered amount of reaction heat value can be observed, what let conclude about lower thermal loads occuences for this kind of solution.

6. Conclusions

Basing on the conducted research and analysis can be stated that:

1. The compression ignition engine works with air excess and is qualitative controlled, it is possible to supplement the supplying air with LPG admixture without negative influence on oxygen content necessary in combustion process.

2. Using the software based on computational fluid dynamics with finite volume method allows process analyzing of the engines in pre-prototype research.

3. Supplementing the air with LPG admixture during compression ignition engine filling results in increase both pressure and temperature during combustion process.

4. The temperature distribution at TDC after combustion initiation shows that in the case of diesel powering the combustion chamber walls (piston) are more heat loaded than in the case of powering LPG admixture or diesel initial injection + LPG.

5. The course of heat development in cylinder shows lengthening the process in the case of diesel powering, what can be the result of fuel cracking, the decay of the higher hydrocarbons to lower, which flame slower.

6. In the following phase the charged engine simulations are planned to be conducted, because modern compression ignition power units are equipped with air supplying systems of such type., and commercial attempts of LPG applying as the admixture to the base fuel cause wear problems.


The investigations described in this paper are a part of the research project No S/WM/2/13 realized at Bialystok University of Technology.


[1.] Lejda, K.; Jaworski, A. 2006. Gas fueling of agricultural engines, Motrol--Motorization and Power Industry in Agriculture 8: 131-138 (in Polish).

[2.] Gis, W.; Menes, W.; Waskiewicz, J. 2012. The analysis of the potential for the use of gaseous fuels in the municipal bus transport in Poland. Part 1, Combustion Engines 1(148): 17-26. PTNSS-2012-111.

[3.] Gis, W.; Menes, W.; Waskiewicz, J. 2012. The analysis of the potential for the use of gaseous fuels in the municipal bus transport in Poland. Part 2, Combustion Engines 3(150): 82-88. PTNSS-2012-SS3-302.

[4.] Lejda, K.; Jaworski, A. 2000. Gas fuelling problems of automotive engines, Proceedings of Vth International Scientific Conference Gas fueled engines 2000. Czestochowa, 229-246 (in Polish).

[5.] Lejda, K.; Jaworski, A. 2003. Gas fuelling problems of agricultural engines, Proceedings of International Conference Motrol Automotive and energetics in agriculture 5: 92-100.

[6.] Flekiewicz, M. 1999. Gas fuelling systems in trucks and buses, Fuels, Oils, Lubricants in Operation 66: 3740.

[7.] Zoltkowski, A.; Jarczewski, M. 2000. Adaptating diesel engines to gaseous fuel, Proceedings of Vth International Scientific Conference Gas fueled engines 2000, Czestochowa, 387-394 (in Polish).

[8.] Luft, S. 2001. The choice of control parameters dual fuel compression ignition engine fueled mainly liquid propane-butane, Journal of Kones Internal Combustion Engines 8(3-4): 53-59.

[9.] Luft, S.; Koziol, S. 2003. Selected results of dual fuel diesel engine fueled mainly fuel LPG, Journal of Kones Internal Combustion Engines 10(1-2): 163-171.

[10.] Kowalewicz, A. 2008. Adaptation of a diesel engine to the natural gas supply, Technical Transactions 7-M: 6778 (in Polish).

[11.] Abd Alla, G.H.; Soliman, H.A.; Badr, O.A.; Abd Rabbo, M.F. 2000. Effect of pilot fuel quantity on the performance of a dual fuel engine, Energy Conversion & Management 41: 559-572.

[12.] Papagiannakis, R.G.; Hountalas, D.T.; Rakopoulos, C.D. 2007. Theoretical study of the effects of pilot fuel quantity and its injection timing on the performance and emissions of a dual fuel diesel engine, Energy Conversion and Management 48: 2951-2961.

[13.] Saleh, H.E. 2008. Effect of variation in LPG composition on emissions and performance in a dual fuel diesel engine, Fuel 87: 3031-3039.

[14.] Zoltowski, A. 2002. Emission problems of the car engines fuelled with LPG, CNG, Journal of Kones Internal Combustion Engines 1-2: 303-310.

[15.] Luft, S.; Michalczewski, A. 2002. Trying to explain to improve the overall efficiency of the diesel engine dually fueled mainly LPG, Journal of Kones Internal Combustion Engines 1-2: 164-172.

[16.] Michalowski R.; Miksiewicz, K.; Tkaczyk, M. 2012. Numerical analysis of comparative combustion engines CNG, Combustion Engines 1(148): 62-71. PTNSS2012-SS1-109.

[17.] Jemni, M.A.; Kantchev, G.; Abid, M.S. 2011. Influence of intake manifold design on in-cylinder flow and engine performances in a bus diesel engine converted to LPG gas fuelled, using CFD analyses and experimental investigations, Energy 36: 2701-2715.

[18.] Ansys Fluent 12.1 Documentation 2009.

[19.] Hardenburg, O.F.; Hase, F.W. 1979. 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, SAE.

[20.] Rychter, T. 1990. Mathematical modeling of the working cycle piston engine, PWN, Warsaw (in Polish).

Received March 27, 2015

Accepted September 17, 2015

G. Pulawski *, D. Szpica **

* Bialystok University of Technology, 45C Wiejska Str., 15-351 Bialystok, Poland, E -mail:

** Bialystok University of Technology, 45C Wiejska Str., 15-351 Bialystok, Poland, E -mail:

Table 1
Properties of fuels [1]

Property Diesel LPG
 (70% of butane)

Heating value, MJ/kg 42.50 45.78

Heating value of 2.47 2.77
stoichiometric mixture,

Liquid/gas density, 815-855 540
kg/[m.sup.3] 2.06

Octane number - 105

Cetane number 50-85 -

Flammability limits, % 0.6 2
 6.5 9
Theoretical air request 14.5 15.5
for combustion, kg/kg

Property Natural gas Biogas Hydrogen

Heating value, MJ/kg 50.00 17.50 120.00

Heating value of 2.74 2.46 3.42
stoichiometric mixture,

Liquid/gas density, 424 71
kg/[m.sup.3] 0.72 1.2 0.09

Octane number 110 134 70

Cetane number - - -

Flammability limits, % 5 5 4
 15 20 75
Theoretical air request 17.2 6.1 34.0
for combustion, kg/kg

Table 2
Engine technical data

Engine type CI, natural aspirated
Number, configure of cylinders 4, in-line
Timing system 2 valves per cylinder
Displacement, [dm.sup.3] 1.896
Bore, mm 79.5
Stroke, mm 89.5
Connecting rod length, mm 150
Compression ratio 18
Maximum power, kW 77 at 4400 rpm
Maximum torque, Nm 250 at 1900 rpm

Table 3
Fuel injection parameters

Parameter Simulation 1 Simulation 2

Fuel type diesel diesel+LPG

Injection angle diesel 20° 10°
[°CA] LPG - 20

Injection start diesel 18° 18°
angle [°CA] before TDC before TDC

 LPG - 60°
 after TDC

Parameter Simulation 3

Fuel type diesel initial
 injection +LPG

Injection angle diesel 1°
[°CA] LPG 40°

Injection start diesel 18° before TDC
angle [°CA] LPG 60° after TDC
COPYRIGHT 2015 Kauno Technologijos Universitetas
No portion of this article can be reproduced without the express written permission from the copyright holder.
Copyright 2021 Gale, Cengage Learning. All rights reserved.

Article Details
Printer friendly Cite/link Email Feedback
Author:Pulawski, G.; Szpica, D.
Article Type:Report
Date:Nov 1, 2015
Previous Article:The castellated beams deflections calculated with theory of composed bars.
Next Article:Interpretation of parameters in strain energy density bone adaptation equation when applied to topology optimization of inert structures.

Terms of use | Privacy policy | Copyright © 2022 Farlex, Inc. | Feedback | For webmasters |