Evaluation of usefulness of mass flow meter to the survey of SI engine cylinder filling in one working cycle

Similar documents
IDENTIFICATION OF FUEL INJECTION CONTROL SYSTEM IN A GDI ENGINE

Fault simulation of the sensors in gasoline engine control system

The influence of fuel injection pump malfunctions of a marine 4-stroke Diesel engine on composition of exhaust gases

ANALYSIS OF THE ENGINE FUELS IMPACT ON CARBON DIOXIDE EMISSIONS

ANALYSIS OF PERFORMANCES OF A DUAL-FUEL TURBOCHARGED COMPRESSION IGNITION ENGINE

Fuel control. The fuel injection system tasks. Starting fuel pump (FP)

Propagation of flame whirl at combustion of lean natural gas charge in a chamber of cylindrical shape

Research in use of fuel conversion adapters in automobiles running on bioethanol and gasoline mixtures

Diagnosing methods common rail fuel injectors

INFLUENCE OF RME FUELS ON CHOSEN VARIABLES OF DIESEL ENGINE SUPPLY PROCESS

MULTIPOINT SPARK IGNITION ENGINE OPERATING ON LEAN MIXTURE

UNIT 2 POWER PLANTS 2.1 INTRODUCTION 2.2 CLASSIFICATION OF IC ENGINES. Objectives. Structure. 2.1 Introduction

DETERMINATION OF THE TOTAL EFFICIENCY IN THE ENGINE WITH DIRECT INJECTION OF THE PETROL

Tips & Technology For Bosch Partners

Throttleless control of SI engine load by fully flexible inlet valve actuation system

MONO - MOTRONIC 0309 En

Experimental Investigation of Acceleration Test in Spark Ignition Engine

Preparing and programming of ESGI 2 LPG supply system manual

MULTIPLE INVESTIGATIONS OF FUME EMISSIONS OF ENGINES WITH AUTOMATIC IGNITION

Kul Internal Combustion Engine Technology. Definition & Classification, Characteristics 2015 Basshuysen 1,2,3,4,5

SAMPLE STUDY MATERIAL

NEW DIAGNOSTIC METHODS OF ALTERNATIVE LPG/CNG INJECTION SYSTEMS

Kubota Engine Training: WG1605, spark ignited

The impact of inlet channel geometry on in-cylinder swirl

INFLUENCE OF THE MARINE 4-STROKE DIESEL ENGINE MALFUNCTIONS ON THE NITRIC OXIDES EMISSION

THE DIAGNOSTIC MODEL PROPOSITION OF THE ENGINE VIBRATION SIGNAL

Installation And Programming Manual of OPTIMA Eco Tec and OPTIMA Pro Tec OBD/CAN

EXPERIMENTAL COMPARATIVE STUDIES OF INJECTION SYSTEMS

COMPUTER SIMULATION POSSIBILITIES IN MODELLING OF IGNITION ADVANCE ANGLE CONTROL IN MOTOR AND AGRICULTURAL VEHICLES

The influence of thermal regime on gasoline direct injection engine performance and emissions

Delomatic 400 (DM 400) based Gas Engine and Combined Heat and Power (CHP) Control and Management System

2.61 Internal Combustion Engines

Analysis of Effect of Throttle Shaft on a Fuel Injection System for ICES

VALVE TIMING DIAGRAM FOR SI ENGINE VALVE TIMING DIAGRAM FOR CI ENGINE

VISUALIZATION IN OF INSIDE CYLINDER PROCESSES IN GASOLINE DIRECT INJECTION ENGINE

WEEK 4 Dynamics of Machinery

Per Andersson and Lars Eriksson

Simple Finite Heat Release Model (SI Engine)

RESEARCH ON INFLUENCE OF SELECTED FAILURES ON THE EXHAUST GAS CONTENT OF SHIP DIESEL ENGINE WORKING ON HEAVY FUEL OIL

COMBUSTION AND EXHAUST EMISSION IN COMPRESSION IGNITION ENGINES WITH DUAL- FUEL SYSTEM

LAMBDA SENSOR CONTROLLER

Foundations of Thermodynamics and Chemistry. 1 Introduction Preface Model-Building Simulation... 5 References...

INFLUENCE OF INTAKE AIR TEMPERATURE AND EXHAUST GAS RECIRCULATION ON HCCI COMBUSTION PROCESS USING BIOETHANOL

CHAPTER 3 EXPERIMENTAL SET-UP AND TECHNIQUES

Simulation Method of Hydraulic Confined Piston Engine

Chapter 6. Supercharging

Error codes Diagnostic plug Read-out Reset Signal Error codes

IGNIJET DUCATI - detailed description

4. With a neat sketch explain in detail about the different types of fuel injection system used in SI engines. (May 2016)

Department of Engineering Science University of Oxford. Particulate Matter Emissions from a Highly Boosted GDI engine

Instruction of connection and programming of the VECTOR controller

The influence of non-cooled exhaust gas recirculation on the diesel engine parameters

The influence of non-cooled exhaust gas recirculation on the indicator diagrams and heat release parameters in diesel engine cylinder

LECTURE NOTES INTERNAL COMBUSTION ENGINES SI AN INTEGRATED EVALUATION

SUCCESSFUL DIESEL COLD START THROUGH PROPER PILOT INJECTION PARAMETERS SELECTION. Aleksey Marchuk, Georgiy Kuharenok, Aleksandr Petruchenko

Comparison of Swirl, Turbulence Generating Devices in Compression ignition Engine

DEUTZ Corporation 914 Gas. Customer / Event DEUTZ Corporation Presentation DATE, 2010

Effect of the boost pressure on basic operating parameters, exhaust emissions and combustion parameters in a dual-fuel compression ignition engine

UNIT IV INTERNAL COMBUSTION ENGINES

Effects of Pre-injection on Combustion Characteristics of a Single-cylinder Diesel Engine

837. Dynamics of hybrid PM/EM electromagnetic valve in SI engines

THE RELATIONSHIP BETWEEN THE FLUCTUATION OF THE INDICATED WORK, THE FLUCTUATION OF THE CRANKSHAFT SPEED AND THE ELECTRIC CURRENT FROM A GENERATING SET

IGNIJET_ Detailed description

Comparative Study Of Four Stroke Diesel And Petrol Engine.

Diagnosing the automobile starting system

EG DYNAMIC user manual

EXHAUST BRAKE SYSTEM MODEL AND TORQUE SIMULATION RESULTS ON A DIESEL SINGLE-CYLINDER ENGINE

Selected aspects of the use of gaseous fuels blends to improve efficiency and emission of SI engine

COMBUSTION in SI ENGINES

Timing is everything with internal combustion engines By: Bernie Thompson

Faraday's Law of Induction

Template for the Storyboard stage

Marcin Zacharewicz 1, Tomasz Kniaziewicz 2

Problem 1 (ECU Priority)

DEPENDENCE OF THE TOXIC COMPONENTS EXHAUST EMISSION FROM THE CAR ENGINE STARTING TEMPERATURE

Normal vs Abnormal Combustion in SI engine. SI Combustion. Turbulent Combustion

Practical Exercise: Computation of the engine output characteristics for a 4-stroke spark ignition engine

INVESTIGATION OF THE FUEL PROPERTY INFLUENCE ON NUMBER OF EMITTED PARTICLES AND THEIR SIZE DISTRIBUTION IN A GASOLINE ENGINE WITH DIRECT INJECTION

Impact of air conditioning system operation on increasing gases emissions from automobile

MIXTURE FORMATION IN SPARK IGNITION ENGINES. Chapter 5

Setup Tabs. Basic Setup: Advanced Setup:

Task 4: Read the texts, look at the illustrations and do the activities below.

Prop-Tech Vacuum Analyzer

GASOLINE DIRECT INJECTION IN SI ENGINES B. PAVAN VISWANADH P. ASHOK KUMAR. Mobile No : Mobile No:

CFD numerical simulation of the indirect cooling system of an internal combustion engine

ME 466 PERFORMANCE OF ROAD VEHICLES 2016 Spring Homework 3 Assigned on Due date:

Development of Rattle Noise Analysis Technology for Column Type Electric Power Steering Systems

ANALYSIS OF MIDDLE RANGE FREQUENCY VIBRATION OF AIRCRAFT RECIPROCATING ENGINE. Tomasz Goliasz Jacek Czarnigowski

ADDIS ABABA UNIVERSITY INSTITUTE OF TECHNOLOGY

Lecture 5. Abnormal Combustion

CHARGING SYSTEM OF SPARK IGNITION ENGINE WITH TWO TURBOCHARGERS

Computer Power. Figure 1 Power-curves from Viper and Venom bottom left and right. (Source: D Quinlan)

EEN-E2002, Gas exchange and supercharging, lecture 4a

Megasquirt EX, Installation Instructions document revision 1.4 (Includes also the version equipped with wasted-spark ignition)

4.0L CEC SYSTEM Jeep Cherokee DESCRIPTION OPERATION FUEL CONTROL DATA SENSORS & SWITCHES

IC ENGINE(4 STROKE) G.H.R.I.E&M JALGAON. Sec.(Mech) Sec.(Mech) Sec.(Mech) Sec.(Mech) Mehta chirag Shah sagar Patel jainish talele amit

IMPROVING ENERGETICAL AND ENVIRONMENTAL PERFORMANCE OF DIESEL ENGINES, BY THE EFFICIENCY SUPERCHARGE PROCESS

Parameter Setting Basic. Voltage Fuel 1 Fuel 2 Ignition 1 Ignition 2 Twin Injector COPYRIGHT 2016 HKS CO.LTD.ALLRIGHT RESERVED

CHAPTER 6 IGNITION SYSTEM

UAV engine EFI components Specifications Manual

Transcription:

Article citation info: WOŁCZYŃSKI, Z. Evaluation of usefulness of mass flow meter to the survey of SI engine cylinder filling in one working cycle. Combustion Engines. 217, 169(2), 187-192. DOI: 1.1926/CE-217-233 Zbigniew WOŁCZYŃSKI CE-217-233 Evaluation of usefulness of mass flow meter to the survey of SI engine cylinder filling in one working cycle The article presents the technique of the usage of mass flow meter to the fill survey of a cylinder in the multicylinder engine, in each working cycles. The result of correctness of the survey was grounded on the comparison of air masses: the one registered with the usage of the flow meter and the other one calculated on the basis of mixture parameters. The flow meter was located in the consolidated part of suction system in front of the throttle valve. ixture parameters were determined on the basis of time of fuel injection and measured λ zm factor of burnt mixture. The studies were carried out in the whole operating range of engine. Discrepancies between the measured and calculated values for 93% of all of the cycles fall within the ambit of ±3% of the actual value. Key words: SI engine, the fill survey, mass flow meter, the phase of fill survey 1. Introduction Filling with the air of particular cylinders of multicylinder SI engine is inimitable in further working cycles. Its reason is a string of hard- and nonmeasurable phenomena e.g. outflow of fresh air parts to the outlet manifold during the valve overlap. In the present article the following question was formed: Is it possible, with the usage of mass flow meter located in front of the throttle valve, to measure with requested correctness the air mass flowing to any cylinder in each of the working cycles. There was also a trial conducted of determining the correctness with which it is possible to measure the filling of this cylinder. Due to it the assumption was adopted that despite several nonmeasurable phenomena influencing the filling, the air mass suck to the cylinder is going to be proportional to the one measured with flow meter. It enabled the implementation of proportionality factor between measured and actual air mass, called the k correction factor. Therefore, the k correction factor ought to be used to revise measured air mass to reach its actual value, i.e. the air mass suck to the cylinder. Non-zero response time of flow meter and phenomena occurring in the suction system dictate to assume that the correction factor of measured air mass is going to have different values for each of the points of engine work. The decisive parameters from which the value of correction factor depends on are: n rotational speed and angle of the α throttle valve. Along with the change of n rotational speed wave phenomena in the suction system change as well and thereby the filling of engine cylinders. They also influence on flow meter indications located under the throttle valve. The location of the throttle valve influences on repression of flow and its pulsations. Therefore, decisive meaning of the opening angle of the α throttle valve to the k correction factor comes from it. Named parameters, i.e. n and α, also influence on the delay of signal of air mass flow m a rate in proportion to the air flow filling the engine cylinder. Therefore, these parameters decide about the range in which on the basis of m a signal the air mass is measured. In terms of above phenomena the method of air mass survey flowing to any cylinder in every filling and the way of their verification were developed. 2. The method of air mass survey and its verification Presented method is used to the air mass survey filling a cylinder of four-cylinder SI engine by way of mass flow meter located in the consolidated part of suction system in front of the throttle valve. The essence of this survey method was graphically presented on the basis of signal processes in the Fig. 1. This picture shows signals in the system of controlling SI engine necessary to determining air mass in further and successive working cycles. The survey consists of several phases. The first one determines once every turn in TDC piston: n rotational speed of engine and opening angle of the α zm throttle valve according to the dependencies: where: t time of rotation 6 n = t (1) zm α l.p. α = (2) where: l.p. number of samples of α signal. Just before the opening of inlet valve the angle of the beginning of phase of φ pn fill survey is determined. For the engine which the studies were made on it is the angle of 672 of the previous cycle. Its value is read from a map on the basis of rotational speed and opening angle of the throttle valve. In the moment of the accomplishment with crankshaft the φ pn location the Azm air mass survey flowing to the engine cylinder begins which finishes half of the rotation later. Azm survey is an integration after the time of signal of air mass flow m a rate according to the formula: φ + 18 φ = dm dt (3) Azm pn pn It is equivalent with determining a zone under the curve of m a strength (in the Fig. 1 bold with the half of the skyblue curve) in the designated borders. Simultaneously durance of injection pulses of engine is measured and after the finishing the each of them a mass of injected fuel is calculated. a COBUSTION ENGINES, 217, 169(2) 187

Evaluation of usefulness of mass flow meter to the survey of SI engine cylinder filling in one working cycle Fig. 1. Processes of chosen signals in the control system of SI engine necessary to determining air mass in each of the cycles of chosen cylinder Fig. 2. Timing of ma signal of the function of crankshaft for different rotational speeds and openings angle of the throttle valve (in the columns: small, medium, big rotation speed; in the rows: small, medium, big opening angle of the throttle valve) 188 COBUSTION ENGINES, 217, 169(2)

t w = t t (4) w2 where: t w time of injection, t w1 time of the beginning of injection, t w2 time of the finishing of injection F1;F2 w1 = 1.842 t 1.411 (5) Each of the injections is assigned to particular working cycle. The fuel injected during the opening of F2 inlet valve (in the Fig. 1 it is a range shown with the bold blue line) and during the injection on the closed inlet valve in the previous F1 cycle (in the Fig. 1 these are pulses shown with the bold pink line) creates fueled-airy mixture with the air which mass was measured in specific cycle. F F1 w = + (6) where: F mass of injected fuel per a cycle, F1 mass of the first injection of fuel in a cycle, F2 mass of the second injection of fuel in a cycle. As a result of burning of mixture exhaust gas are going to appear which analysis will allow to determine λ factor of burnt mixture. The exhaust gas will drain during the opening of exhaust valve (in the fig. 1 it is a range shown with the bold black line). The average value of λ signal from the analyzer (bold part of red curve) equals the λ zm factor of burnt mixture in particular cycle according to the formula: F2 rotation of crankshaft. Therefore, this signal was plotted in many points involving the entire zone of engine work and the angle of the beginning of the φ pn fill survey phase was determined in each of them (Fig. 2). The angle of the end of the fill phase is located 18º further. It is a result of this that flow meter signal during a cycle (two rotations four filling of cylinders) holds four characteristic parts of process. They represent the fillings of the following four engine cylinders. On the basis of φ pn value in many points of work there was created a map of this angle for tested engine cylinder as dependency from n rotational speed and α zm opening angle of the throttle valve (Fig. 3). λ zm 12 o λ o (7) 54 = l.p. where: l.p. number of λ signal samples in the range of counting the average, i.e. from 54 of particular cycle to 12 of the following one. The delay with which the average value of λ signal in proportion to the range of the opening of exhaust valve is calculated results from the time necessary for the outflow of exhaust gas and their analysis. On the basis of fuel mass in particular cycle and λ factor of mixture in the same cycle Ao air mass as an actual air mass was measured according to the formula: Fig. 3. ap of fill survey phase Similarly, i.e. for many points of engine work there was created a map of the k correction factor (Fig. 4) determining its value according to the dependency (9). Ao = F 14, 7 λzm (8) The proportion of calculated A mass to measured Azm mass equals the k correction factor: Ao k = (9) A zm The presented method of the air survey was developed as a result of the observation of influence of rotational speed and the opening angle of the throttle valve on the process of flow meter signal. It was observed that these borders submit to displacement along with the rotational speed of engine and the opening of the throttle valve. This displacement is different for the particular points of engine work and is connected with wave phenomena in the suction system. It was also noticed that for each of the points of engine work variability time of m a signal equals 18 of Fig. 4. ap of the k correction factor COBUSTION ENGINES, 217, 169(2) 189

3. Studies In order to develop the map: the angle of the beginning of the φ pn fill survey phase and the k correction factor and determining an error of air mass survey with the use of mass flow meter the studies were planned and carried out. The measurement site was built (Fig. 5). hundred thousands of samples of every signal. The following voltage signals in the function of time were registered (Fig. 7): λ composition of burnt mixture measured in the branch of outlet manifold of the fourth cylinder. m a air mass flow rate p u negative pressure in the intake manifold α angle of the throttle valve measured with TPS sensing device which is a part of the control system torque on the brake shaft U w voltage signal of controlling the power of injector of the first or the fourth cylinder U z voltage signal of controlling the current in the primary winding of ignition coil for the first and the fourth cylinder swk transformed with equipment into the bistate signal of the location of engine crankshaft. Fig. 5. easurement site anufactured in series SI engine with multipoint, indirect injection of fuel with the capacity of 16 cm 3 was the object of study. It was installed on the stationary measurement site and was managed with the factory control system. Control system was equipped with measuring-controlling interface enabling a survey of all of the signals occurring in it. Interface belongs to the equipment of measurement site. Furthermore in the suction system of engine between air filter and throttle valve additionally HF-5 mass flow meter of BOSCH company was installed. The broadband lambda sensor was the second additional sensing device installed. It was used as a sensing device of analyzer of the composition of EXA7λ burnt mixture of HORIBA company. It was installed in the branch of outlet manifold of the fourth cylinder. The PC computer equipped with the OCTOPUS CS838 measuring card with eight measuring channels and GageScope software of GAGE company was used to the registration of any signals in the time function. Each of the measuring channels is equipped fourteen-bit A/D converter. The studies were carried out on the entire zone of points of engine work determined with the rotational speed and the opening angle of the throttle valve, i.e. in 215 points (Fig. 6). n [obr/min] 58 53 48 43 38 33 28 23 18 13 8 2 4 6 8 1 α [%] Fig. 6. Points of engine work in which the studies were carried out Signals mentioned below were being registered in every point. They were sampled by 1 µs during 1 s, i.e. one Fig. 7. Exemplary process of registered signals 4. The results of studies Registered in the studies signals were submitted for initial modifications in order to prepare them for further analysis: bistate signals: U w, U z, swk were transformed into zero-one signals (Fig. 8) the remaining signals (analogue ones) were averaged with the use of method of average from 2 samples (Fig. 9) on the basis of the swk zero-one signal the angle of the location of shaft in the working cycle ( 72 ) was assigned to every sample of registered signals (Fig. 1). swk, U w, U z λ, m a, pu, α,, Uw, Uz, swk [j.k.] 9 8 7 6 5 4 3 2 1-1 -2-3 -4 57 62 67 72 77 82 57 62 67 72 77 82 Fig. 8. Bistate signals before and after transformations into zero-one signals (on the upper part of figure registered bistate signals, on the lower part corresponding to the zero-one signal) 19 COBUSTION ENGINES, 217, 169(2)

λ, m a, p u, α, [j.k.] 1 8 6 4 2-2 -4 57 62 67 72 77 82 Fig. 9. Analogue signals before and after averaging with the use of method of average from 2 samples (black lines represent averaged processes) 7 6 72 time of injection of t w fuel [ms] on the basis of U w signal and corresponding to this time dose of injected F fuel [mg] according to (4), (5), (6), mass of measured Azm air [mg] on the basis of m a signal from HF-5 flow meter in designated borders of survey according to (3), measured value of excess air α zm [-] factor on the basis of λ signal from EXA7λ meter according to (7), measured air mass k Azm rectified with the correction factor. Later on air masses were compared: the measured one with the actual one and the error of survey in every entire registered cycle was calculated according to the formula: k zm Ao δ A = 1% (1) Ao The results of this error are shown in the Fig. 11. λ, m a, α, U w, U z, swk [j.k.] λ, m a, α, Uw, Uz, swk 54 5 4 36 3 2 18 1 2 25 3 35 4 45 5 55 6 65 7 7 6 5 4 3 2 1 18 36 54 72 φ [ OWK] φ [ OWK] δa [%] 8% 7% 6% 5% 4% 3% 2% 1% % -1% -2% -3% -4% -5% -6% -7% -8% Fig. 11. Errors of air mass survey suck to the studied engine cylinder Every point in the picture represents the value of error of air mass survey in one filling cycle. The picture includes values of error of survey for 4895 filling cycles. The number of cycles for particular ranges of error was shown in the Fig. 12. Cycles in which error of survey is bigger than 3% represents less than 7% of all of them. Fig. 1. Assignment of registered in time samples of signals to temporary angular position of engine shaft (chart above presents value of angle of φ crankshaft in a cycle against the background of registered in time samples of signals, chart below present the same samples of signals in the function of φ angle) odifications mentioned above enabled carrying out of analysis of registered signals for every entire working cycle. This analysis consists in determining all of the parameters of a cycle necessary to the calculation of air mass on the basis of m a signal and error of determination of this mass. They are: n rotational speed [rpm] for a rotation (from TDC to TDC) on the basis of swk according to (1), angular location of the α zm [%] throttle valve on the basis of signal from TPS (% - entirely closed, 1% - entirely open] according to (2), the opening angle of φ pn and the end of φ pn +18 of the phase of air mass survey [ ] (retrieved from map on the basis of n and α zm ), Number of cycles 35 3 25 2 15 1 5-8% -7% -6% -5% -4% -3% -2% -1% % 1% 2% 3% 4% 5% 6% 7% 8% δ A [%] Fig. 12. Layout of the number of cycles for which error of survey falls within the ambit of particular borders 5. Conclusions 1. With the usage of mass flow meter located in the consolidated part of suction system in front of the throttle COBUSTION ENGINES, 217, 169(2) 191

valve it is possible to measure air mass filling one engine cylinder in each of the working cycles. 2. The analysis of flow meter signal have to take place in the angular range of the rotation of crankshaft ancillary to the rotational speed of engine and the opening angle of the throttle valve. 3. For all of the analyzed work cycles from the entire zone of engine work the error of air mass survey did not pass 8%. 4. easuring air mass with the usage of presented method the error of survey is less: from 3% for 93% of cycles, from 4% for 98% of cycles, from 5% for 99,4% of cycles. Nomenclature SI spark ignition n rotational speed α signal of location of the throttle valve α zm opening angle of the throttle valve moment on the brake shaft m a air mass flow rate Azm measured air mass A actual air mass φ angular position of crankshaft in a cycle φ pn opening angle of the filling survey phase t w time of injection t w1 time of the beginning of injection time of the finishing of injection t w2 F F1 F2 λ λ zm U w U z swk p u fuel mass in a cycle fuel mass injected in the first injection per a cycle fuel mass injected in the second injection per a cycle signal of mixture composition with EXA7λ analyzer measured value of mixture composition voltage signal of controlling the injector of the first or the fourth cylinder voltage signal of controlling the ignition coil for the first and the fourth cylinder signal from the sensing device of crankshaft negative pressure in the intake manifold Bibliography [1] NITA, J. Niepowtarzalność cyklowego składu mieszanki w silniku GDI. Czasopismo techniczne. 212, 5, 229-241. [2] PRZYBYŁA, G., POSTRZEDNIK S. Niepowtarzalność cykli pracy silnika przy jego wysokiej prędkości obrotowej oraz różnym obciążeniu. Czasopismo techniczne. 28, 7, 151-159. [3] DYKIER,., FLEKIEWICZ,. Analiza częstotliwości drgań czynnika roboczego w układzie dolotowym silnika SI. Czasopismo techniczne. 28, 8, 173-183. [4] NITA, J., WOŁCZYŃSKI Z. Ocena fluktuacji składu mieszanki wysterowanej z cyklu na cykl w silniku z pośrednim wtryskiem benzyny. Silniki Spalinowe. 213, 3, 12-11. [5] KOORSKA, I. (red.) Eksperymentalna ocena możliwości poprawy sprawności silnika benzynowego zasilanego LPG poprzez optymalizację jego sterowania. UTH Radom. 216, 18. Zbigniew Wołczyński, DEng. Faculty of echanical Engineering at Kazimierz Pulaski University of Technology and Humanities in Radom. e-mail: Z.Wolczynski@uthrad.pl 192 COBUSTION ENGINES, 217, 169(2)