 Open Access
 Authors : Lister. Munodawafa. Dzikiti , Patrick Mukumba
 Paper ID : IJERTV9IS080295
 Volume & Issue : Volume 09, Issue 08 (August 2020)
 Published (First Online): 03092020
 ISSN (Online) : 22780181
 Publisher Name : IJERT
 License: This work is licensed under a Creative Commons Attribution 4.0 International License
CFD Contextual Modelling of Biogas Combustion in Internal Combustion Engine: A Review
Lister. Munodawafa. Dzikiti 1 Patrick Mukumba1
1Department of Physics, University of Fort Hare, Alice South Africa;
* Correspondence: Lister Munodawafa Dzikiti.
Abstract: The perpetual use of petroleum products as fossil fuel nowadays is the primary cause of many economical, ecological and environmental problems. The exhaust emissions from fossilfuelled engines world over are a major cause of hazardous emissions such as NOx (oxides of nitrogen), CO (Carbon monoxide), HC (Hydro Carbons) and PM (Particulate Matter). Furthermore high temperature combustion in fossilfuelled IC engines has been identified as the major contributor of greenhouse gases characterised by high pollution emissions. Hence this review paper focuses on findings from contemporary renewable energy researchers based on attempting to substitute fossilfuelled IC engines with biogas fuelled ones. Recent studies have been focusing on combustion modelling and technology that would lead to the improvement in biogas combustion efficiency. Confirmation from literature indicates that biogas usage as an engine fuel has the greatest potential of substituting the environmental unfriendly fuel derived from petroleum products due to its very low CO2 emissions. CFD simulations based on biogas combustion have further proved that high energetic contents of biogas similar to that of natural gas can be combusted for the generation of mechanical energy. The SA automotive industry is growing fast but still based on fossilfuelled IC engines. Hence, the findings from literature suggest the need to develop CFD modelling of biogas combustion in internal combustion engine in the context of SA automotive industry
Keywords: Computational Fluid Dynamics, Compression Ignition, internal Combustion, Navier Stokes Equations, Spark Ignition, Combustion Modelling
 INTRODUCTION
Energy is universally recognized as a prime agent of economic development [1]. Contemporary energy researchers have already verified the existence of a relationship between energy availability and economic activity [2]. Rapid growth in world economies together with associated exponential increase in population has led to an abrupt increase in energy demand [3]. About 80% of overall energy demand is mainly derived from fossil fuel [4].
In the recent years fossil fuel prices have been on the rise due to limited supply and deposits especially of crude oil as well as significant increase in demand of petroleum fuels [5]. The reality of risks associated with environmental degradation and climate change due to petroleum fuel usage are now more apparent [2]. Carbon monoxide, sulphur dioxide, nitrogen oxides and particles are undesirable emissions associated with burning fossil fuels. These compounds are toxic, contribute to acid rain and smog and can ultimately cause respiratory problems. The automotive industry is worlds primary consumer of energy extracted from fossil fuels [6]. The current rate of energy consumption from fossil fuel sources means that all reservoirs will be depleted by 2042 [7]. Furthermore high temperature combustion of fossil fuels has been identified as the major contributor of greenhouse gases due to its high pollution emissions [8]. Currently, a fossil fuel derivative like diesel (a mixture of hydrocarbons with C15 – C18 carbon atoms and an approximate calorific value of about 11,000 kcals/kg) is used as fuel in compression ignition engines. However, the perpetual use of petroleum products like diesel as fossil fuel nowadays is the primary cause of many economical, ecological and environmental problems. For instance, the exhaust emissions from the diesel engines world over are a major cause of respiratory problems and, heavy pollutions with hazardous emissions such as (oxides of nitrogen (NOx), Carbon monoxide (CO), Hydro Carbon (HC) and Particulate Matter (PM) [9].
The primary reason for this pollution is due to the heterogeneous mixture of air and diesel in the combustion chamber. This fuel will be exhausted soon because of its excessive usages and nonrenewable nature. The most critical issue in the present is the replacement of fossil fuels with renewable sources. Renewable energy resources are potentially among the most effective and efficient solutions to nonrenewable sources such as fossil fuel. Suggesting that the use of alternative energy source in internal combustion (IC) engines will have greater impact on energy generation and consumption [10]. It is against such a catastrophic background that an alternative automotive fuel resource is required to meet the requirements of the current environment to reduce hazardous emissions without compromising CI engine efficiency. Hence, biogas is a better alternative fuel, which can be used to circumvent the above crisis [11]. Furthermore, reducing consumption of oilderived energy products will greatly minimize greenhouse gas emissions [12]. A renewable energy source like biogas is what the future holds even for the automotive industry. It is therefore imperative to simulate and analyse the fundamental impacts of firing biogas in an IC engine using CFD analysis methods. This will help to accurately predict the physical conditions required before engineers can modify the engine into accepting biogas as fuel.
 COMBUSTION ENGINE SYSTEMS
There are two main types of combustion engines namely internal combustion engine (IC) and external combustion (EC) engines [11]. The compression ignition (CI) engine and Spark ignition (SI) engine exist as the two types of IC engines in the autoindustry. The combustion process in a SI and CI engine can easily be described through use of Wiebe equation shown [13].
m1
b
b
x 1 exp a 0
t
(1)
Where = crank angle, 0 = the position of the crankshaft angle at which the spark was ignited,
t = the total combustion
duration, a and m are adjustable parameters with a representing the maximum fraction of energy released/burned mass, ( xb 0
to xb 1). However in mechanical engineering, the CI engine is preferred to SI due to its high thermal efficiency, low fuel
consumption and very low emissions of uncombusted hydrocarbons [14].
In an IC engine ignition the combustion process is generally characterised by spark timing as well as duration and completeness of combustion [11].

 Internal Combustion Engine
Internal combustion (IC) engine is widely used for fuel combustion particularly in the automotive industry [15]. It is important to improve the engine efficiency through reduction in fuel consumption, hazardous emissions, noise pollution and the negative effect these leaves on the environment [16].
The IC engine exists in many different types and sizes but can all be categorised as twostroke cycle or four stroke cycle. The following are examples of internal combustion engines:
 Compression ignition (CI) also called single fuel nonpremixed compression ignition engine. (SFNPCI)
 Spark ignition (SI) also called the single fuel premixed spark ignition engine (SFPSI)
 Homogenous charge compression ignition (HCCI)
 Reactivity controlled compression ignition (RCCI) also referred to as the dual fuel premixed compression ignition engine. (DFPCI)
 Variable compression ratio (VCR)
 Opposed piston internal combustion engine (OPICE)
The IC engine can either be run by Otto or theDiesel principle. For Otto IC engines, a spark from a spark plug [17] ignites the premixed airfuel mixture. For the Diesel principle, the air is compressed beforehand in the cylinder and the incoming fuel spray is ignited by the high pressure and temperature [18]. In sparkignited (SI/Otto) engine, the fourstroke cycle according to Chiodi [19], starts with the piston positioned in top dead centre (TDC) and as the piston travels down towards bottom dead centre (BDC) the intake valve opens, letting the fresh charge of airfuel mixture to enter the cylinder during the intake stroke, as shown in Figure 1a. With the piston at BDC the intake valve closes and as the piston travels towards TDC again the fresh charge is compressed during the compression stroke, see Figure 1b. As the piston reaches TDC the compressed airfuel mixture is ignited by the spark plug and the chemical energy of the fuel is converted to heat during combustion. This increases the cylinder gas temperature and pressure, thereby adding work to the crankshaft during the power stroke or expansion stroke, see Figure 1c. When the piston reaches BDC again the exhaust valve opens and as the piston travels towards TDC the exhaust gas is pushed out from the cylinder during the exhaust stroke, see Figure 1d. When the piston has reached TDC the cycle is restarted again. In one power cycle the crankshaft has done two full revolutions (720 crank angle degrees or CAD) and the piston has travelled up and down the cylinder four times, therefore the name four stroke cycles.
Figure 1: The four stroke cycle (adapted from Semin, 2008)
 Internal Combustion Engine

 POTENTIAL OF BIOGAS AS FUEL FOR IC ENGINES
Biogas is a renewable energy fuel derived from anaerobic digestion of organic matter such as municipal and agricultural wastes [13]. Biogas is mainly composed of a mixture of methane (CH4) and carbon dioxide (CO2) [12]. The heating value of biogas depends primarily on its CH4 content [19]. The methane composition of biogas produced in most anaerobic digesters ranges between 40% 75% [20]. Under standard conditions of temperature and pressure, methane has a density of approximately 0.75kg/m3 [15]. Adding carbon dioxide that is slightly heavier, the biogas density reaches about 1.15kg/m3 [12]. Pure methane has an upper calorific value of 39.8MJ/m3, which corresponds to 11.06Kwh/m3 [19]. The quality of biogas from a biogas digester plant (BDP) varies with period of digestion (adapted from Reddy & Reddy, 2014). Biogas of higher quality burns with a blue flame in the presence of excess oxygen (O2) [21].
The use of biogas as an engine fuel is advantageous due to its low CO2 emissions [13]. Biogas has the capability of lowering the ignition limit and can increase the power output by augmenting the energy density at lean mixtures hence, increasing the biogastocarbon ratio [22] Biogas cannot be used directly in a compression ignition engine because of its high auto ignition temperature [23]. Huang & Crookes [14] recommended pilot ignition techniques in which diesel engines is fitted with spark plugs or a use of small amount of diesel fuel to ignite the biogas.
The combustion of a mixture of biogas and air produces low gas temperature inside the combustor compared to natural gas (Mihic, 2004). The biogas combustion temperature depends on the carbon dioxide to methane ratio (CO2/CH4) (Razbani, et al., 2011). The biogas temperature decrease by about 37% for every 22% increase in CO2/CH4 ratio.
Huang & Crookes [14] observed that the NOx emissions inside a combustor decreases when diesel or gasoline fuel is replaced with biogas due to lower heating value. Noting that thermal NOx emission is a function of gas temperature, thus high gas temperature produces high NOx (mainly NO) emissions [24]. The presence of noncombustible constituents in the biogas such as carbon dioxides and nitrogen also reduces the temperature of the flame and consequently the NO mass fractions [15].
The NOx emissions at the exit point of the combustor decreases with increase in biogas CO2/CH4 ratio but the emissions of CO2 at the exit of the combustor increase with an increase of CO2 fractions in the biogas fuels (higher CO2 input) [22]. The presence of nonreactive gas (CO2) with methane cools down the reaction by absorbing energy from the combustion and modifies the reaction zone by reducing flame burning rate and velocity [15]. The emissions from biogas combustion vary with the composition of biogas. The higher CO2 percentage in the biogas fuel, the lower the NOx emissions and the higher the CO2 emissions at the exit point of the combustor. The CO emissions also decrease by increasing the CO2 /CH4 ratio of the biogas fuel [24].
 COMPUTATIONAL FLUID DYNAMICS
Many computational tools aid in the design and development of IC chambers and analysis of combustion process with the aim of improving combustion efficiency [25]. Computational fluid dynamics (CFD) is a computational tool that can be used in calculating combustion dependent parameters as well as combustion process. CFD uses a numerical method in solving a set of nonlinear partial differential equations to predict fluid flow. The complex engineering problems solved using CFD approaches involve turbulence [11], hot spots analysis, compression ratios, ignition limits, ignition timing [26], emission pollutants and performance of combustors. [15].
The fundamental governing equations of computational fluid dynamics include continuity, momentum and energy equations [17]. These equations are focused on three basic parameters:


 Mass is conserved. (continuity equation)
 Fnet = ma (Newtons second law), momentum is conserved
 Energy is conserved.

The NavierStokes (NS) equations are a set of nonlinear partial differential equations used to describe the flow of a fluid whose stress depends linearly on velocity gradients and pressures [17]. The NS equation ignores the fact that the fluid is made up of discrete molecules. To quantify the transition between laminar and turbulent flows then Reynolds number is preferably used [27]. The conservation laws are a statement of the fact that the rate of change of mass, momentum, or energy in a certain volume is equal to the rate at which it enters the borders of the volume plus the rate at which it is created inside [3].

 COMPUTATIONAL MODELLING
The main concept of CFD is the discretized solution of a set of partial differential equations commonly known as the NavierStokes equations [21]. These equations are discretized in time, by solving the equations in small time steps, and in space by dividing the domain in a large number of small computational cells or control volumes (CV) [15]. The space discretized solution domain is commonly referred to as the computational grid or mesh [28].
The method commonly used for space discretization is called the Finite Volume Method (FVM), and for a static mesh (no moving boundaries), it is based on the integral form of the conservation equation over a control volume fixed in space [18]. The simulation of a fluid dynamics in a combustion engine requires moving boundaries of the solution domain, because the structure of the flow field is very much dependent on the moving piston and valves [15].. When the mesh has moving boundaries the integral form of the conservation equation for a tensorial property defined per unit mass in an arbitrary moving volume V bounded by a closed surface S is given in equation 2.
d dV
dS (u u
)
dS q S dV
(2)
V
V
dt S
b S
V
where is the density, u is the fluid velocity, ub
sources/sinks of respectively [19].
is the boundary velocity and q
and
S are the surface and volume
The CFD equations govern the flow of fluid [2]. The conservation form of CFD governing equations is obtained directly from a control volume that is fixed in space rather han moving with the fluid [17]. When the volume is fixed in space the concern is on i) the flux of mass into and out of volume. ii) momentum into and out of the volume. iii) energy into and out of the volume [12].
 Governing Equations
The governing equations used for biogas combustion include equations of conservation of mass, momentum, energy, equations of
the turbulent kinetic energy , and the dissipation rate of the turbulent kinetic energy the standard turbulence
model), the mixture fraction equation (mixture fraction/PDF (probability density function)) model for nonpremixed combustion modeling), and the equation of state. For steady turbulent non premixed combustion [29] the time averaged gas phase equations are given below.
THE CONTINUITY EQUATION
The continuity equation governs the conservation of mass, which means that the rate of change of mass in an arbitrary control volume must be equal to the total mass flow over the control volume boundaries:
(U ) 0
t
CONSERVATION OF MOMENTUM
(3)
The momentum equation governs the conservation of linear and angular momentum. According to Newtons second law, the rate of change of momentum on a fluid particle equals the sum of forces acting on that particle:
u (uu) g
t (4)
Where the two right hand terms denote the gravitational body force and viscous stress tensor respectively. CONSERVATION OF ENERGY
This equation is based on the first law of thermodynamics which states that the rate of change of total energy of a fluid element (internal, kinetic, mechanical, chemical energy) must equal the net rate of heat and work flow over its boundaries [30]. However, an alternative formulation of the energy equation, used in the ANSYS FLUENT code, is the enthalpy equation:
h (uh)
t t
)h Q Dp
dt
(5)
where is the density, u is the fluid flow vector, g is the body force, is the stress tensor, h is the enthalpy, is the
thermal diffusivity for enthalpy,
t is the turbulent thermal diffusivity, Q is the energy source term, and p is the fluid
pressure.
MIXTURE FRACTION EQUATION
t f
fu
j
j
x
x
t
j S
m
m
(6)
The mixture fraction f is given by
x j
x j
K
K
Z
f
 ZK ,O
(7)
ZK ,F ZK ,O
Where Z K is the element mass fraction of element k . Subscripts F and O denote fuel and oxidizer inlet stream values,
respectively.
Sm is the term source that includes the mixture of chemical species.
EQUATION OF STATE
Assuming that the biogas behaves ideally, then the equation for the temperature and the pressure can be given as
T 1 M h (T ) i E(T )
And
R0 i
i i
(8)
p R T i
0
0
i Mi
(9)
where
Mi is the molar mass of species i and
R0 is the ideal gas constant
Thermodynamic relation between state variables is given by
e e(T , p)
and for a calorically perfect gas this relation would be
(10)
where
e cvT cv = specific heat at constant volume.
(11)
THE VELOCITY EQUATION
V 0.2gd w ,
b
m s1
(12)
THE PRESSURE EQUATION
P wgh ,
N m2
(13)
THE REYNOLDS NUMBER
R bVD , dimensionless (14)
E
g
The NS equations are very difficult to solve. As the Reynolds number is increased, the scale of the interesting dynamics gets smaller so that most solutions of the full NS equations are done at Reynolds numbers of 1 to 10,000 for simple geometries.
COMPRESSION RATIO ( CR )
V
V
R
R
C VBDC
TDC
(15)
whereVBDC is the total volume of combustion chamber when piston is at the bottom of its stroke, VTDC is the total volume of
the combustion chamber when the piston is at the top of its stroke.
THERMAL EFFICIENCY ( T )
C
C
1
1
1
T k 1
R
(16)
Where k represents ratio of specific heat, a property of air.
 ZK ,O
 Governing Equations
 Turbulence Modelling
The flow in the engine room is time dependent, unsteady, compressible and turbulent [28]. Turbulent flow is a flow that is irregular, random and chaotic [17]. It can be described as a state of continuous instability in the flow, where it is still possible to separate the fluctuations from the mean flow properties [21]. Turbulence is associated with eddies which affects the mean flow, where the range of the scales is very large, from the smallest turbulent eddies characterized by Kolmogorov microscales, to the flow features comparable with the size of the geometry. The largest eddies extract their energy from the mean flow [31]. These eddies transfer energy to smaller eddies in a cascade process. The large scale eddies have an orientation imposed by the mean
flow, but the smaller eddies will notremember their origin and orientation and will behave isotropic that is independent of direction. The scale of these eddies, referred to as Kolmogorovs microscale, are small and dissipative forces prevails [19]. In this regime, the kinetic energy is destroyed by viscous forces, and viscosity and dissipation affect the length scales [18]. The cascade process, from the biggest to the smallest eddies, occurs over a wide length scale spectrum and a very fine mesh is needed to resolve the smallest Kolmogorov eddies.
 Averaging methods for NS equation
A very fine resolution in time is also needed, since turbulent flow is always unsteady. The requirements on mesh resolution and timestep size for an exact solution of the flow, puts very high demands on computer resources, rendering it unsuitable for engineering applications [17]. This can be solved by simulating turbulent flow using statistical methods.
Using statistical or averaging methods, the local value of the variable can be separated into the mean and the fluctuation around the mean which then makes it possible to derive the equations for the mean properties themselves [19]. The averaging methods can appear in two versions, Reynolds averaging and densityweighted Favre averaging.
Most CFD packages (including ANSYS FLUENT) are to solve the Reynolds Averaged Navier Stokes (RANS) equations [18].
RANS equations govern the mean velocity u( y)
and pressure but are primarily used to solve mean velocity only. If we
formally average the NS equations and simplify for this geometry we arrive at the following
u`v` 1 p
2 u( y)
y x
v y 2
(17)
Where
u`v`is known as the Reynolds stress tensor. Applying the averaging procedures to the NS equations will introduce the
unknown term Reynolds stress tensor. This generates a NS system of equations with more unknowns, situation commonly referred to as the closure problem. To close the system, further modeling is necessary to solve the Reynolds stress tensor (Alias, 2008). A widely used method to handle the closure problem is the Boussinesq approximation, which states that the Reynold stresses are
linked to mean rates of deformation by introducing a turbulent viscosity vt
[19].The turbulent viscosity vt can be evaluated through use of the k – turbulence model.
 The k – Turbulence Model
This twoequation model makes use of two partial differential equations for the determination of the turbulent viscosity (Alias, 2008). k and are related to the primitive variables through
~ 1 u~“2
Where k is the turbulent kinetic energy And
2
~ u“ : u“
(18)
(19)
Where is the turbulent dissipation rate and
v
L
is the laminar viscosity.
The turbulent viscosity vt can then be calculated using the PrandtlKolmogorov equation [17] This equation expresses vt as
a function of the turbulent kinetic energy k and its dissipation rate is given as
2
2
k
k
~
vt Cv ~
(20)
Where C is an empirically determined model constant (C = 0.09) while k and are obtained by the solution of their respective transport equations.
 Averaging methods for NS equation
 Combustion Modelling
 COMPUTATIONAL MODELLING
The combustion model used for the simulations in this study is the
b
Weller combustion model because it is the standard
combustion model incorporated in the ANSYS FLUENT CFD code [28]. This model is based on the evaluation of the dimensionless variable b that describes the species concentration of the combustion reactants in each computational cell [19]. The value of b ranges from 0 to 1. Where a value of 0 means that the combustion process is completed in the whole cell and a value of 1 means that no combustion process has taken place.
The transport equation to describe the evolution of b in time and space states:
~
~~ ~
~
b
t
ub bb
( u
 ) min , eq Su b
(21)
Where is the density, u
is the density in the unburnt mixture, u is the fluid flow vector,
t
b = , (22)
h
Where t
is the turbulent viscosity, h
is the turbulent Prandtl coefficient, Su
is the laminar flame speed. The sign Â¯
indicates the average and ~ indicates the Favre average based on the density (i.e. ~ u /
). The subscript eq signifies the
equilibrium value which serves as an upper limit of . Lastly, is defined as the local flame wrinkling factor and it is defined
as the ratio between the turbulent flame speed
St and the laminar flame speed Su :
St
Su
(23)
In solving for one can use a transport equation, which requires reproducing all the details of the flame structure or the transient response of the turbulent flame speed [12]. Under these conditions the transport equation need not be solved at all, instead a simple analytical (or algebraic) expression for can be used in the b equation.
Hence, the algebraic expression of is the given as:
u`
u
u
eq 1 A S
(24)
where A is a free coefficient,
u` is the turbulence intensity and
Su is the laminar flame speed.
The general hydrocarbon stoichiometric equation in a combustion process can be expressed as shown equation 25 [28].
C H (n m)(O 3.76N ) nCO m H O 3.76(n m)N
(25)
n m 4 2 2 2 2 2 4 2
For biogas stoichiometry the general form of combustion equation is given as
(iCH 4 jH 2 kCO2 ) (2i 0.5 j)(1.0O2 3.8N2 )
(i k)CO2 (2i j)(H 2O) (7.6i 1.9 j)(N2 )
(26)
Where i, j, and k represent the mole fraction of CH4, hydrogen (H2) and CO2 respectively. The detailed study carried out on models used for internal combustion, spray, and autoignition to using CFD are presented below.

 BIOGAS COMBUSTION SIMULATION STUDIES
The literature presented here shall predicate more on combustion modelling which have been studied, simulated and experimentally validated world over.
Porpatham [32] examined and modelled the effect of methane concentration on combustion in biogas fuelled SI engine. The observations from the study indicated that raising methane concentration in biogas led to the significant increase in SI engine performance and sudden reduction in total hydrocarbons (THC) emissions. Furthermore, the abrupt decrease in CO2 concentration from biogas resulted in greater methane/ oxygen content in the charge leading to a faster combustion rate and tremendous power output at specific equivalence ratio. The simulated outcome of this study was validated and its results were found to be in agreement with experimental.
Another study reported that an increase in the concentration of CO2 in biogas fuel led to a decrease in engine power output thus compromising engine performance. This also resulted in more production of pollutant compared to methane rich biogas. The efficiency of the IC was compromised with an increased fraction of CO2, which led to more consumption of biogas fuel for similar power output. Although the mole fraction of NOx decreased significantly, the exhaust CO and THC concentration increased with increase in mole fraction of CO2 [33].
Intensive studies on the BiogasDiesel dual fuel combustion system was done by Shaik et al [34] using CFD analysis aided with ANSYS FLUENT Software package. In this study the benefits of biogas substitution in an IC engine in terms of combustion flame velocity (CFV) and probable formation of NOx was simulated for five compression ratios. The turbulence model was modified to aid in simulating turbulent kinetic energy (TKE) and turbulent dissipation rate (TDR) for the dual fuel system. The effect of compression ratio in ignition engine system in which biogas was premixed with air and diesel as primary fuel in this duel fuel set up was simulated.
Guessab et al [35] carried out a turbulent swirling flow analysis while applying a three dimensional RANSRSM turbulence procedures for diffusion chamber, using biogas as fuels. The study was conducted using a cantype gas turbine combustor model. The flame was stabilised by gradually adding CO2 at lean fuel/air equivalence ratios needed in the reduction of NOx emissions. There was a significant change in aerodynamic and thermal properties when methane/natural gas was substituted with biogas. This was due to increase in fuel mass flow variations and lower heating value associated with biogas. This is because the two gaseous fuels differ in terms of their chemical composition. The results of the study were considered costeffective as they significantly proved the possibilities of using methane or biogas in turbulent combustion systems.
In a study based on combustion of synthetic gases in an IC engine, Arroyo, [10] concluded that CO /CO2 percentage is the major contributory factor of increased exhaust gas pollutants if compared with biogas and other related fuels. It was also observed in this study that the exhaust gas concentration of THC was relatively low due to the smallunburnt methane content. Although the presence of such chemical species like excess air and CO2 provided mitigation for NOx emissions. However greater mole fraction of H2 in the synthetic gas led to an increased NOx emissions being execrated by high combustion flame temperatures associated with the presence of H2.
Huang [36] studied a variable compression ratio (VCR) single cylinder SI engine fuelled with biogas. The simulation of the biogas operated VCR SI engine was achieved with different mixtures of domestic natural gas and CO2. The content of CO2 was gradually increased from 0% to 40% in order to replicate a biogas fuel. A range of air/fuel ratios were also studied from the lean to the rich operating limit at various speeds and variable (CR). The results suggested that CO2 in biogas fuel greatly influenced the reduction of the toxic NOx emissions and also enhance the SI engine in operating at a high CR. However, it is unfortunate to observe that incomplete combustion as evidenced by high emissions of THC was prominent in the presence of excessive amounts of CO2. Meanwhile the emissions of CO remained unaffected by CR and crankshaft speed but CO correlation existed with the relative airtofuel ratio. Jung et al [37] studies suggested that biogas combustion in engines could be improved by increasing the methane percentage composition of the biogas. The study invetigated the extension of lean operational limit up to a relative air/fuel ratio of 1:5 with variable biogas compositions and for CH4: CO2 volume ratios of up to 70%: 30 % with no engine knocking. The emissions of NOx were reduced to nearzero due to the decrease in the combustion temperature, which was aided by use of the lean combustion mixture.
Bedoya et al [38] studied the effects of biogas composition on homogeneous charge compression ignition (HCCI) combustion system. They adopted a numerical methodological approach which combines chemical kinetics and CFD numerical simulation to analyse the trends in emissions and combustion parameters. They numerically analysed the effect of biogas fuel composition using the HCCI engine in terms of its combustion and emission parameters. Obtained results showed that NOx emissions were not primarily influenced by biogas composition however, CO and noncombusted THC emissions increased with increase in the percentage composition of CO2. The numerical results also indicated that biogas composition could be used as an effective way of controlling the onset of combustion as well as combustion phasing in biogas fuelled HCCI engines thus allowing for combustion stability and operational safety with low and high equivalence ratios. Abade [39] study numerically simulated CH4/ CO2 fuel mixture in order to investigate its effect on emissions and catalyst performance of a biogasfuelled SI engine. The study revealed that increasing the CO2 content led to a sharp increase in THC and a decrease in NOx at the exhaust manifold. This concurred with the findings of Bedoya et al [38] study. Furthermore, it was noted in another study that doubling catalyst residence time will cause a rise in THC and CO emissions by a factor of 2. It was also noted that effective catalytic conversion of THC emissions depended more on the relative proportions of THC, NOx and CO emissions at the level of the exhaust manifold [39]. Toppo [40] study presented a CFD model for an advanced combustion processes in CI engine. Temperature profiling inside the combustion chamber and NOx emission at the exhaust manifold was compared for Jathropha fuel with conventional diesel fuel. The simulated results found were subjected to experimental validation for both Jathropha and diesel fuels in a CI engine system. The simulation was done in ANSYS FLUENT code using nonpremixed combustion modelling to imitate a real incylinder combustion phenomenon. Twodimensional combustion system with mesh deformation was adopted in the simulation. The dimensions of the computational domain were aligned with the actual CI test engine. The findings of this study only revealed that the NOx emitted for jatropha was more compared to that emitted by diesel. However, the study made no mention of a comparative study of biogas and diesel/petrol in terms NOx emissions.
Noor et al [41] study presented a stepwise CFD simulation method for nonpremixed moderate and intense low and diluted (MILD) oxygen biogas fuelled combustion furnace. A Nonpremixed combustion model with a turbulent realizable standard kepsilon was adopted in their simulation. The CFD simulation was achieved perfectly for the MILD regime and the maximumtoaverage temperature ratio was found to be less than the expected 23%. Furthermore, it was also noted that meshing quality had optimum influence on convergence or divergence of calculations. Indicating that good CFD results is a product of high quality meshing achievable if and only if calculations converge. VitÃ¡zek et al [42] study was aimed at developing a new biogas mixture thermodynamics methodology using tabular exact parameters to calculate operating values of biogas composition. A biogas combustion mathematical model was developed and elaborated in addition to a simplified simulation diagram that was developed to depict the combustion device. The simulation diagram with indications of mass and energy movements was the kind, which would support the combustion of gaseous fuel as the limiting reactant. The findings revealed that using biogas directly proved difficult due to the presence of some H2O, CO2, N2, H2S, THC, organosilicon compounds impurities. The presence of such impurities would make the potential of biogas usage as fuel to be costly due to additional purification cost.
Palaniswamy et al [11] used CFD software to investigate the maximum distance and optimal injector angle needed for a better biogas/air stoichiometric mixture at the intake manifold. Their findings indicated that the biogas introduced at 45 injector angle and 336mm distance from the axis of the inlet valve led to a better mixture. The experimentally validated results proved that 50% of blended biogas with a supercharged system was 25% more mechanically efficient compared to the natural aspiration system. The results also reported less CO emission values for a supercharged system.
Mameri et al [43] CFD simulation study revealed that the presence of CO2 in biogas has a thermal dissociation effects. It absorbs energy from the combustion of methane, which leads to high exhaust temperature. The simulation study also showed convincingly that the useful energetic part of the biogas is methane. Therefore, suffice to say that biogas can be a better fossil fuel substitute provided CFD simulation feasibility studies are to be carried out prior to industrial implementation. KuÂ´znia et al [44] used a Chemked II program to model the biogas combustion kinetic process. The results of the study suggested that biogas combustion with greater CH4 content does not necessarily lead to an increase in exhaust CO2 gas concentration but rather improves engine performance.
Pablo [45] conducted an investigation using a diesel CI engine with a maximum output power 8.5 kW. The CI was converted into SI engine to allow gaseous fuel conditions and usage to be appropriate. In this study three fuel types were considered namely: simulated biogas, biogas (25% CH4) and biogas (50% CH4) by volume. The findings of this study suggested that CO2 presence allowed the SI engine to operate at very high compression ratios (CR) even under usual combustion conditions. The study also found that enrichment of biogas with CH4 resulted in an abrupt decrease in NOx, CO and unburnt CH4/THC emissions at the level of the exhaust manifold. This also concurred with findings from other studies [10, 37 39].
 CONCLUSION
Carbon monoxide, sulphur dioxide, nitrogen oxides and particles are undesirable emissions associated with burning fossil fuels. The growing need to substitute fossil fuels with biogas for internal combustion engine particularly in automotive industry motivated the present study. This study successfully simulated and analysed previous studies focused on the fundamental impacts of firing biogas in an IC engine using CFD. From the study, it was gathered that the presence of CO2 in biogas has a thermal dissociation effects during combustion in an IC engine. It absorbs energy from the combustion of methane, which leads to high exhaust temperature. The simulation study also showed convincingly that the useful energetic part of the biogas is methane. Also the present study revealed that using biogas directly proved difficult due to the presence of some H2O, CO2, N2, H2S, THC, organosilicon compounds impurities. The presence of such impurities would make the potential of biogas usage as fuel to be costly due to additional purification cost. Therefore, suffice to say that biogas can be a better fossil fuel substitute provided CFD simulation feasibility studies are to be carried out prior to industrial implementation.
 RECOMMENDATION
Computational work specifically using CFD simulation is now becoming more and more important due to its lower cost and acceptable accuracy with minimum error. A purely theoretical and numerical CFD simulation was employed in this study to investigate the combustion properties of biogas in an IC engine. The present study has shown that achieving unique operating conditions for biogasuelled engines, particularly in IC engines is not that simple. This is mainly due to the significant difference in biogas stoichiometric composition. Hence, an experimental investigation of biogas composition and its direct effect on emissions in IC engines should be carried out prior to its sole use as substitute for fossil fuel. From a South African perspective, the country has a large potential of biogas reservoirs whose energy can be exploited in the country`s growing automotive industry. Hence, South Africa is perhaps better suited to explore the possibilities of utilising biogas energy as an alternative energy source especially for its automotive industry. In South Africa`s automotive industry a drastic reduction in exhaust CO2 concentration will go a long way in lowering greenhouse gas emissions.
ACKNOWLEDGMENTS REFERENCES
 Giles. D. E., Som, S. and Aggarwal. S. K. (2006). NOx Emission Characteristics of Counterflow Syngas Diffusion Flames with Airstream Dilution.
Fuel Vol.85 (1213): 172942.
 Barik, D., Sah, S., and Murugan, S. (2013). Biogas Production and Storage for Fueling Internal Combustion Engines. International Journal of Emerging Technology and Advanced Engineering 3 (3): 193202.
 Noor, M. M., Wandel, A.P. and Yusaf, T. (2013). Analysis of recirculation zone and ignition position of nonpremixed bluffbody for biogas MILD combustion. International Journal of Automotive and Mechanical Engineering, 8, p.1176.
 Maczulak, A. (2010). Renewable Energy: Sources and Methods. New York: Facts on File Inc.
 Hairuddin. A. A, Yusaf. T. and Wandel, A. P. (2016). Singlezone zerodimensional model study for dieselfuelled homogeneous charge compression ignition (HCCI) engines using Cantera, International Journal of Automotive and Mechanical Engineering, Vol. 13 (2), 3309 3328.
 IEA. (2011). World energy outlook. Paris: International Energy Agency.
 Shafiee, S., & Topal, E. (2009). When will fossil fuel reserves be diminished? Energy Policy, 37(1), 181189.
 Talibi M, (2017). Combustion and exhaust emission characteristics, and incylinder gas composition, of hydrogen enriched biogas mixtures in a diesel engine Energy. International Journal of Engineering Research & Technology Vol.124: 397e412.
 Huang. H, (2017). The potentials for improving combustion performance and emissions in diesel engines by fuelling butanol/diesel/PODE34 blends Energy Procedia 105 914 920 ScienceDirect, The 8th International Conference on Applied Energy.
 Arroyo J, (2014). Combustion behavior of a spark ignition engine fueled with synthetic gases derived from biogas. Fuel. International journal of engineering research & technology Vol.117, 508
 Palaniswamy1. D, Ramesh. G, Sivasankaran. S, Sooryaprakash. K, (2016). CFD Analysis for Homogenous Effect of Biogas and Air in the Intake Manifold of Dual Fuel CI Engine International Journal of Advanced Engineering Technology eissn 09763945.
 Suzuki.A.B.P, Fernandes.D.M, Faria.R.A.P & Vidal.T.C.M, (2011). Use of biogas in internal combustion engines, Brazilian Journal of Applied Technology for Agricultural Science, Vol. 4 (1), 221237.
 Carrera. J. (2013). Numerical study on the combustion process of a biogas spark. Ignition engine. Thermal science, Vol.17, 24154.
 Huang, J. and Crookes, R. J. (1998). Assessment of simulated biogas as a fuel for the spark ignition engine. Fuel, London, v. 77, n. 15, p. 17931801, 1998.
 Mihic, S., (2004). Biogas Fuel for internal combustion Engines, Annals of Faculty Engineering, Hunedoara. Vol 2 (3), 411.
 Yusaf, T., Baker, P., Hamawand, I. and Noor, M.M., (2013). Effect of compressed natural gas mixing on the engine performance and emissions.
International Journal of Automotive and Mechanical Engineering, 8, p.1416.
 Hussain, S.M., BSP, D.K. and KVK, D.R., (2012). CFD analysis of combustion and emissions to study the effect of compression ratio and biogas substitution in a diesel engine with experimental verification. International Journal of Engineering Science and Technology, 4(2), pp.473492.
 Jemni, M.A., Kantchev, G. and Abid, M.S., (2011). Influence of intake manifold design on incylinder flow and engine performances in a bus diesel engine converted to LPG gas fuelled, using CFD analyses and experimental investigations. Energy, 36(5), pp.27012715.
 Chiodi, M., (2011). An innovative 3DCFDapproach towards virtual development of internal combustion engines. Braunschweig: Vieweg+ Teubner Verlag
 Mukumba P, Makaka G, Mamphweli S & Misi S. 2013. A possible design and justification for a biogas plant at Nyazura Adventist High School, Rusape, Zimbabwe. Journal of Energy in Southern Africa, Vol 24 (4).
 Semin, R.A.B., (2008). A technical review of compressed natural gas as an alternative fuel for internal combustion engines. American J. of Engineering and Applied Sciences, 1(4), pp.302311
 Anggono, W., Wardana, I., Lawes, K., Hughes, K. J.,Wahyudi, S., and Hamidi, N., (2012). Laminar Burning Characteristics of BiogasAir Mixtures in Spark Ignited Premix Combustion. Journal of Applied Sciences Research 8 (8): 412632
 Marchaim, U. (1992), Biogas processes for sustainable development FAO Agricultural Services Bull, pp. 95.
 Razbani, O., Mirzamohammad, N., and Assadi, M., (2011). Literature Review and Road Map for Using Biogas in Internal Combustion Engines. 3rd Internal Conference on Applied Energy, Perugia, Italy.
 Rodrigues, S. M, Saslow, L. R., Garcia, N., John, O. P., Keltner D. (2009) Oxytocin receptor genetic variation relates to empathy and stress reactivity in humans, Proceedings of the National Academy of Sciences 106(50):2143741.
 Koten, H., Mustafa, Y., Zaafer Gul, M. (2014). Compressed BiogasDiesel DualFuel Engine Optimization Study for Ultralow Emission, Advances in Mechanical Engineering, vol 2014, 18.
 Mukaro. R, (2014). Digital, statistical and wavelet study of turbulence flow structure in laboratory plunging water waves, PhD thesis. University of KwaZulu Natal, Durban.
 Noor, M.M., Wandel, A.P. and Yusaf, T., (2013, July). Detail guide for CFD on the simulation of biogas combustion in bluffbody mild burner. In
Proceedings of the 2nd International Conference of Mechanical Engineering Research (ICMER 2013) (pp. 125). Universiti Malaysia Pahang.
 Mare, F., Jones, W. P., Menzies, K. R. (2004) Large Eddy Simulation of a Model Gas Turbine Combustor, Combustion and Flame 137(3): 278294.
 Reddy, R. and Reddy, P., (2014). Analysis of producer gas carburetor for different airfuel ratios using CFD. International Journal of Research in Engineering and Technology, 3, pp.470474.
 Bicsak, G., Hornyak, A., Veress, A. (2012) Numerical Simulation of Combustion Processes in a Gas Turbine, AIP Conference Proceedings 1493, 140 (2012); https://doi.org/10.1063/1.4765482.
 Porpatham.E, Ramesh.A, Nagalingam. B, (2012). Effect of compression ratio on the performance and combustion of a biogas fuelled spark ignition engine. Fuel. International journal of engineering research & technology,Vol. 95, 24756.
 Lee J, (2010). A study on performance and emissions of a 4stroke IC engine operating on landfill gas with the addition of H2, co and syngas. New york: columbia university.
 Shaik, S. Dewir, Y.H. Singh, N. Nicholas A. Micropropagation and bioreactor studies of the medicinally important plant Lessertia (Sutherlandia) frutescens LSouth African Journal of Botany, 76 (2010), pp. 180186.
 Guessab.A, Baki.T, Mansour.C, (2016). Combustion of methane and biogas fuels in gas turbine cantype combustor model, journal of applied fluid mechanics, vol. 9 (5), 22292238.
 Huang.J, (1999). Sparkignition engine erformance with simulated biogas: a comparison with gasoline and natural gas. Fuel and energy abstracts. Vol.40, 2839.
 Jung. C, Park J, Song S, (2015). Performance and NOx emissions of a biogasfueled turbocharged internal combustion engine. Energy (oxford). Vol.86, 18695.
 Bedoya.I.D, Saxena. S, Cadavid.F.J, Dibble. R.W, (2013). Numerical analysis of biogas composition effects on combustion parameters and emissions in biogas fueled HCCI engines for power. Journal of engineering gas turbines power. Vol.135, 071503.
 Abader.R, (2014). Study on biogasfueled SI engines: effects of fuel composition on emissions and catalyst performance. University of Toronto.
 Toppo. I.O, (2013). CFD analysis of combustion characteristics of jathropha in compression ignition engine. International journal of engineering research & technology, Vol. 2 (10), issn: 22780181.
 Noor, M.M., Wandel, A.P. and Yusaf, T., (2014). The simulation of biogas combustion in a mild burner. Journal of Mechanical Engineering and Sciences, 6(1), pp.9951013.
 VitÃ¡zek, J. KlÃºik, D. UhrinovÃ¡, Z. MikulovÃ¡, M & Moji, (2016).Thermodynamics of combustion gases from biogas, Res. Agr. Eng., 62 (Special Issue): S8S13.
 Mameri A, et al., (2016). Numerical investigation of counterflow diffusion flame of biogas hydrogen blends: effects of biogas composition, hydrogen enrichment and scalar dissipation rate on flame structure and emissions. International journal of hydrogen energy, Vol.41, 201122.
 KuÂ´znia. M, Jerzaka. W, Lykob. P & Sikora. J, (2015). analysis of the combustion products of biogas produced from organic municipal waste. Journal of power technologies, Vol. 95 (2), 158165.
 Pablo.G, (2015) Spark ignition engine performance and emissions in a high compression engine using biogas and methane mixtures without knock occurrence. Thermal science, Vol. 19, 191930.