28InternalCombustionEnginespressure (about 180 bars)after burning of the whole dose (Figure 20).Velocity of incrementofthemeancharge temperature inside the caloric chamber depends on thevalueofthe fueldose(Figure21)andforbiggerdosethequickerincrementofthetemperatureisobservedburningof.fueldose30002500Im,=0,035g-m,=0,04g1500m,=0,045g10005000246time [ms]Figure21.Variation ofthetemperatureinthecaloricchamberfordifferent doseof CNG1081341e+01252.196221e+0155en6Se-015e+003.0100D00tei02702.1eioetoe2.12e00z.12eim181e00-81e-0015001Se0sLii.190-00800(b)(a)764.4196116-011691.4439110-02CethsSeneee-02.790-0Leve2.120-18e1anedSettLSena1.1e*e)ad0(c)(d)Figure 22. Temperature in the caloric chamber after initiation of combustion: a) 0.5 ms, b) 0.6ms,c)1.85ms,d)7.4ms
28 Internal Combustion Engines pressure (about 180 bars) after burning of the whole dose (Figure 20). Velocity of increment of the mean charge temperature inside the caloric chamber depends on the value of the fuel dose (Figure 21) and for bigger dose the quicker increment of the temperature is observed. Figure 21. Variation of the temperature in the caloric chamber for different dose of CNG Figure 22. Temperature in the caloric chamber after initiation of combustion: a) 0.5 ms, b) 0.6 ms, c) 1.85 ms, d) 7.4 ms This image cannot currently be displayed. (a) (b) (c) (d)
Factors Determing Ignition and Efficient Combustion in Modern Engines Operating on Gaseous Fuels29Thedose of fuel influenceson variationof allthermodynamicparameters.Theinitiationofthecombustionprocess lasted about 0.5ms for all doseof the fuel.The completecombustionof all doses of the fuel without swirl and tumblefollows after 4mswithassumptionofheattransfertothewalls.FourslidesinFigure22showthespreadingoftheflameinthecaloricchamber from the spark plugto the walls almost spherically.Themaximum of temperaturenear sparkplugamounts almost 3600K and after combustion process decreasesto2700K.10.3.VerificationofignitionmodellingThe initial simulations oftheCNG combustion was carried out on themodel of the chamberused for the experimental tests on the Schlieren stand in a steady state initial conditions.Thechamber had the volume equalled 100 cm3 with diameter D-80 mm and width B=20 mm.Theinitiation of theignitionfollowedinthecentreofthechamberbytwothinelectrodes.The chamber was filled by natural gas at 5 bars and =1.4.The initial temperature of thecharge amounted 300 K, so this required much more electrical energy than forfiring engineThe ignition energy was simulated as additional internal energy in the centre of thecombustionchamber.TheLESmodelforfullypremixedchargewasusedintheCFDopensourceprogramOpenFOAM.The classical idea is to use a filter whichallowsfortheseparation of large and small length scales in the flow-field. Applying the filtering operatorto the Navier-Stokes equations provides a new equation governing the large scales exceptfor oneterminvolvingthe small velocity scale.Themodel of combustion chamberwascreated by hexahedron cells and contained 68x68x32 cells.Calculations of combustionprocess were carried out in 64-bit Linux system with visualisation of results by use Paraviewsoftware.Thecombustionprocess inthechamberlasted a longtime (above50ms),becauseof absenceof thegas motion.The oxidation ofmethanewas simulatedby the OpenFOAMcombustion procedure in Xoodles module.Thermodynamicproperties of the chargewerecalculated byusingJANAFtables.Increaseof pressurein the flatcombustionchamberwithout initial swirl or"tumble" of the charge is shown in Figure 23.CNGignition inSchlierenchamberinitialpressure-5bars5040Jvegl anssed302010002040501030time [ms]Figure 23. Increase of pressure in the chamber after ignition
Factors Determing Ignition and Efficient Combustion in Modern Engines Operating on Gaseous Fuels 29 The dose of fuel influences on variation of all thermodynamic parameters. The initiation of the combustion process lasted about 0.5 ms for all dose of the fuel. The complete combustion of all doses of the fuel without swirl and tumble follows after 4 ms with assumption of heat transfer to the walls. Four slides in Figure 22 show the spreading of the flame in the caloric chamber from the spark plug to the walls almost spherically. The maximum of temperature near spark plug amounts almost 3600 K and after combustion process decreases to 2700 K. 10.3. Verification of ignition modelling The initial simulations of the CNG combustion was carried out on the model of the chamber used for the experimental tests on the Schlieren stand in a steady state initial conditions. The chamber had the volume equalled 100 cm3 with diameter D=80 mm and width B=20 mm. The initiation of the ignition followed in the centre of the chamber by two thin electrodes. The chamber was filled by natural gas at 5 bars and =1.4. The initial temperature of the charge amounted 300 K, so this required much more electrical energy than for firing engine. The ignition energy was simulated as additional internal energy in the centre of the combustion chamber. The LES model for fully premixed charge was used in the CFD open source program OpenFOAM. The classical idea is to use a filter which allows for the separation of large and small length scales in the flow-field. Applying the filtering operator to the Navier-Stokes equations provides a new equation governing the large scales except for one term involving the small velocity scale. The model of combustion chamber was created by hexahedron cells and contained 68x68x32 cells. Calculations of combustion process were carried out in 64-bit Linux system with visualisation of results by use Paraview software. The combustion process in the chamber lasted a long time (above 50 ms), because of absence of the gas motion. The oxidation of methane was simulated by the OpenFOAM combustion procedure in Xoodles module. Thermodynamic properties of the charge were calculated by using JANAF tables. Increase of pressure in the flat combustion chamber without initial swirl or “tumble” of the charge is shown in Figure 23. Figure 23. Increase of pressure in the chamber after ignition
30 Internal Combustion EnginesThe combustion process involves the change of thermodynamic parameters of the gas,whichcanbeobserved bymoving offlame with different temperature,pressure and densityin burned and unburned spaces.Full combustion of the methane-airmixture lasts longerthan in the real engine combustion chamber at the same geometry of the combustionchamber.The propagation of chemical reactions is radial and the thick boundary of thecombustion (about8mm)isobservedbecauseof theleanmixture.Propagationoftheflamecausestheradial compression ofthegasbetween unburned andburned regionsandthinarea of twice higher density is formed.Figure 9 shows distribution of gas density in thechamberafter18msfromstartofignition.Redcolourindicatesdensityonthelevel0.0118g/cm and bluecolouronly 0.005g/cm?.Figure 24. Gas density and absolute gas velocity after 18 ms from beginning of ignitionCombustion process in the narrow area takes place with turbulent velocity.Turbulencecauses penetration of the flame into the unburned mixture with velocity higher thanlaminar combustion speed.For themethane-air stoichiometric mixture the combustionlaminarspeedamountsonly40cm/s.Fortheconsideredcasetheabsolutevelocityofcombustionintheflameregionamountsabout80m/sasoneisshowninFig.10.Howevertotal combustion speed isvery low and is close to the laminar speed of methane-airmixture0.4m/s.Experimental tests on the Schlieren stand done by Sendyka and Noga [11] showed alsoradial propagation of the flame defined by the change of the charge density.Figure 25shows thefilms of theflamepropagation in thechamber at3,7,40and 54ms after start ofthe ignition, respectively.The ignition of the CNG and air mixture with initial pressure 5bars and initial temperature 300K was initiated by twothin electrodes in centre of thecombustion chamber.The chargewasfully premixed with airexcessratio X=1.4.Theflameisdistortedbytouchingintothequartzglass inthechamber,whichis observedbyhellcircleinsidetheblackcircle.Thechange of gasdensity influencesonthedistortionofthelaserbeam and photos showdevelopment of the flame during combustion process.Theexperimental testprovestheresult obtained fromsimulationbyusingLEScombustionmodel intheOpenFOAMprogram
30 Internal Combustion Engines The combustion process involves the change of thermodynamic parameters of the gas, which can be observed by moving of flame with different temperature, pressure and density in burned and unburned spaces. Full combustion of the methane-air mixture lasts longer than in the real engine combustion chamber at the same geometry of the combustion chamber. The propagation of chemical reactions is radial and the thick boundary of the combustion (about 8 mm) is observed because of the lean mixture. Propagation of the flame causes the radial compression of the gas between unburned and burned regions and thin area of twice higher density is formed. Figure 9 shows distribution of gas density in the chamber after 18 ms from start of ignition. Red colour indicates density on the level 0.0118 g/cm3 and blue colour only 0.005 g/cm3. Figure 24. Gas density and absolute gas velocity after 18 ms from beginning of ignition Combustion process in the narrow area takes place with turbulent velocity. Turbulence causes penetration of the flame into the unburned mixture with velocity higher than laminar combustion speed. For the methane-air stoichiometric mixture the combustion laminar speed amounts only 40 cm/s. For the considered case the absolute velocity of combustion in the flame region amounts about 80 m/s as one is shown in Fig.10. However, total combustion speed is very low and is close to the laminar speed of methane-air mixture 0.4 m/s. Experimental tests on the Schlieren stand done by Sendyka and Noga [11] showed also radial propagation of the flame defined by the change of the charge density. Figure 25 shows the films of the flame propagation in the chamber at 3, 7, 40 and 54 ms after start of the ignition, respectively. The ignition of the CNG and air mixture with initial pressure 5 bars and initial temperature 300 K was initiated by two thin electrodes in centre of the combustion chamber. The charge was fully premixed with air excess ratio =1.4. The flame is distorted by touching into the quartz glass in the chamber, which is observed by hell circle inside the black circle. The change of gas density influences on the distortion of the laser beam and photos show development of the flame during combustion process. The experimental test proves the result obtained from simulation by using LES combustion model in the OpenFOAM program
Factors Determing Ignition and Efficient Combustion in Modern Engines Operating on Gaseous Fuels31Figure 25. Schlieren stand - combustion boundary of the flame after 3, 7,40 and 54 ms [15]Bothsimulationandexperimentdonotshowdeviation of thespherical combustionflameThe experiment demonstrated velocity ofcombustion in radial direction of value 40 cm/s.10.4.MixturemotionandignitionThe most important factor influencing on the ignition is the charge motion through thesparkplug.Twokindsof motions were considered:swirland tumblecaused byvalveandinlet profile, combustion chamber and squish.The combustion process is strongly connectedwith turbulence of the charge and only small part is the laminar speed of the totalcombustion velocity.Simulation was carried out in the rectangular space with centrallocationofthe spark plug.Themesh ofthe combustion chambermodel with lengthandwidth 5 cm and height 3 cm was divided into 288000 cells with rectangular prism (NX=80,NY=80 and NZ-45).The calculations were carried out in transient conditions(initial timestep le-6 s in time t=5 ms).The spark plug was located in the centre of the calculation spaceandtheobjectoftheelectrodeswascreatedbyCADsystem.Themeshintheregionofthespark plug electrodes contains fine grids with cell length equal 0.3 mm in x and y axis.At the first the ignition of CNG was simulated with,initial tumble"y-250 rad/s and p=20bars.The chargewithvelocityabout15m/s flewthrough thegapofthe sparkplug causingthe propagation of the flame inside the chamber.The simulation of combustion and gasmovementwascarried outalsobyPhoenics,whichtakes intoaccountturbulencemodelandsimple combustion of compressible fluid. The charge motion is connected with highturbulenceand this causes also the higher combustion rate.DistributionofthecombustionproductsinthemodelledspaceisshowninFigure26at0.5ms and 1.2ms after start of the ignition,respectively.After short time (about 1 ms)thewhole chargeisburned inthe calculation space.Thehigher flowvelocityisbetween theelectrodes of the spark plug. The other simulation was carried out for the central swirlaround thesparkplugwith swirl velocity15m/s on themean radius1.5cm.In this case theinteraction of the electrode shape is seen -the propagation of the flame is faster in theopened site of the electrodes.Figure 27 presents development of combustion process after1and4msfrombeginningoftheignition.The swirl in the chamber influences on the irregular propagation on the flame and extendsthe combustion process.Even after 4ms the combustion of themethane is not full Velocity
Factors Determing Ignition and Efficient Combustion in Modern Engines Operating on Gaseous Fuels 31 Figure 25. Schlieren stand – combustion boundary of the flame after 3, 7, 40 and 54 ms [15] Both simulation and experiment do not show deviation of the spherical combustion flame. The experiment demonstrated velocity of combustion in radial direction of value 40 cm/s. 10.4. Mixture motion and ignition The most important factor influencing on the ignition is the charge motion through the spark plug. Two kinds of motions were considered: swirl and tumble caused by valve and inlet profile, combustion chamber and squish. The combustion process is strongly connected with turbulence of the charge and only small part is the laminar speed of the total combustion velocity. Simulation was carried out in the rectangular space with central location of the spark plug. The mesh of the combustion chamber model with length and width 5 cm and height 3 cm was divided into 288000 cells with rectangular prism (NX=80, NY=80 and NZ=45). The calculations were carried out in transient conditions (initial time step 1e-6 s in time t=5 ms). The spark plug was located in the centre of the calculation space and the object of the electrodes was created by CAD system. The mesh in the region of the spark plug electrodes contains fine grids with cell length equal 0.3 mm in x and y axis. At the first the ignition of CNG was simulated with „initial tumble” y =250 rad/s and p=20 bars. The charge with velocity about 15 m/s flew through the gap of the spark plug causing the propagation of the flame inside the chamber. The simulation of combustion and gas movement was carried out also by Phoenics, which takes into account turbulence model and simple combustion of compressible fluid. The charge motion is connected with high turbulence and this causes also the higher combustion rate. Distribution of the combustion products in the modelled space is shown in Figure 26 at 0.5 ms and 1.2 ms after start of the ignition, respectively. After short time (about 1 ms) the whole charge is burned in the calculation space. The higher flow velocity is between the electrodes of the spark plug. The other simulation was carried out for the central swirl around the spark plug with swirl velocity 15 m/s on the mean radius 1.5 cm. In this case the interaction of the electrode shape is seen – the propagation of the flame is faster in the opened site of the electrodes. Figure 27 presents development of combustion process after 1 and 4 ms from beginning of the ignition. The swirl in the chamber influences on the irregular propagation on the flame and extends the combustion process. Even after 4 ms the combustion of the methane is not full. Velocity
32InternalCombustionEnginesof the gas flow in the spark plug gap is smaller than in the "tumble" case. For this reason thepropagationofthecombustionproducts andflameisnotuniform.IBEFigure26.Combustionproducts with initial"tumble"chargemotion after0.5ms and after1.2msLUREA21E4IE-E野5Figure27.Combustionproducts withinitial"swirl"chargemotionafter1ms and after4ms11.ConclusionsThe chapter contains results of theoretical, modelling and experimental work considered tofactors, which have verybig impact on the ignition of gaseous fuels in combustion engines.Onthefact of more and more important roleofgaseous engines,particularlythose fuelledbynatural gas, definition ofgood conditionsfor ignition of gaseous fuels is one of thetaskof development of modern spark ignition gaseous engines,particularly with high chargingratio.Experimental works with CNG ignition were done in the caloric chamber, however inconditions closed to real conditions of engine work. On the presented considerations onecandrawsomeconclusionsandremarks:Gaseous fuels, such as CNG requires higher electric energy delivered by the ignitionsystem. Higher pressure in the combustion chamber increases internal energy near thespark plug and requires also higher secondaryvoltage of the ignition coil.Forgaseousleaner mixtures an ignition system with higher energy is needed (above 60m))
32 Internal Combustion Engines of the gas flow in the spark plug gap is smaller than in the “tumble” case. For this reason the propagation of the combustion products and flame is not uniform. Figure 26. Combustion products with initial “tumble” charge motion after 0.5 ms and after 1.2 ms Figure 27. Combustion products with initial “swirl” charge motion after 1 ms and after 4 ms 11. Conclusions The chapter contains results of theoretical, modelling and experimental work considered to factors, which have very big impact on the ignition of gaseous fuels in combustion engines. On the fact of more and more important role of gaseous engines, particularly those fuelled by natural gas, definition of good conditions for ignition of gaseous fuels is one of the task of development of modern spark ignition gaseous engines, particularly with high charging ratio. Experimental works with CNG ignition were done in the caloric chamber, however in conditions closed to real conditions of engine work. On the presented considerations one can draw some conclusions and remarks: 1. Gaseous fuels, such as CNG requires higher electric energy delivered by the ignition system. Higher pressure in the combustion chamber increases internal energy near the spark plug and requires also higher secondary voltage of the ignition coil. For gaseous leaner mixtures an ignition system with higher energy is needed (above 60 mJ)