CFD Analysis of Petrol Internal Combustion Engine: Mechanical Engineering University of Technology, Baghdad-Iraq
CFD Analysis of Petrol Internal Combustion Engine: Mechanical Engineering University of Technology, Baghdad-Iraq
CFD Analysis of Petrol Internal Combustion Engine: Mechanical Engineering University of Technology, Baghdad-Iraq
Journal of University of Babylon for Engineering Sciences by University of Babylon is licensed under
a Creative Commons Attribution 4.0 International License.
132
Journal of University of Babylon for Engineering Sciences, Vol. (26), No. (9): 2018.
1. Introduction
The internal combustion engine (ICE) is considered more important machine that used to provide
mechanical power by converting the fuel power, this takes place with many complex processes. The
modeling of internal combustion engine (ICE) represents one of the highest level of complexity and a
challenging task. This is due to more parameters that play a key role of operation such as, fuel injection,
flame propagation, ignition process, and variety of fluid properties due to unsteady state. Therefore, the
study and development of internal combustion engines need to characterize deeply what happens from the
complex physical phenomenon within engine, especially with the advancement of technology.
The flow characteristic plays one of vital roles in the engine performance. So, several researchers
have been involved in the study of ICE, Heywood [1] mentioned that the air turbulence is better for air -
fuel mixture. Kuo [2] predicted the gas pressure changes within the cylinder by a simple thermodynamic
principle model. Most significant weakness of this model is that it does not take the changes in fuel
reactivity during combustion into account. Movva [3] simulated the fluid flow in combustion chamber of
ICE by wave action simulation, with FORTRAN code, and obtained the pressure distribution, and velocity
at variant boundary condition. Omaraa [4] studied the cold flow in a dual engine, used fluent 2
diamensional transient simulation, and obtained pressure, velocity, and temperature with crank angle by
segment analysis into specific angles. Morauszki et.al [5] investigated the fluid properties in GDI engine
with numerical computation using ANSYS and shows flow pattern, pressure during compression and power
stroke. Ravichandra et.al [6] converted CI into SI by using dynamic mesh technique with fluent 16.0 for
engine numerically and experimentally, recording the temperature and pressure inside combustion
chamber. Abdul Rahiman et.al [7] studied an intake manifold by CFD and showed the variation of Swirl
Ratio (SR) inside the cylinder with respect to crank angle for different manifolds.
This paper represents a modern method in terms of modeling of geometry and 3D simulation by
using the finite volume method. The objective of this research is to develop simulation to calculate the flow
characteristic in combustion chamber of spark ignition engine by using dynamic mesh technique to
visualize the flow within combustion chamber by using ICE code with engine speed 2500 r.p.m. The
velocity, pressure and temperature contours will be studied with respect to different crank angles.
2. Methodology
2.1 CFD Tool
The design and manufacture of internal combustion engines are under significant pressure for
improvement. The generation of engines requires being light, reliable, robust, flexible, and powerful.
Innovative engine designs will be required for satisfying these requirements. The ability to accurately the
performance of multiple engine designs is too much crucial and critical, because IC engines consist of
complex fluid dynamic interactions between air flow, fuel injection, moving parts, and combustion. Using
CFD results, the flow phenomena can be visualized on a 3D geometry and analyzed numerically, providing
tremendous insight into the complex interactions that occur inside the engine. CFD simulation is used as a
part of the design process in automotive engineering, especially with the rise of modern technology.
133
Journal of University of Babylon for Engineering Sciences, Vol. (26), No. (9): 2018.
134
Journal of University of Babylon for Engineering Sciences, Vol. (26), No. (9): 2018.
The events are specified for one complete engine cycle, according to the engine operation
𝐴 𝐴2
𝑃𝑆 = 𝐿 + (1 − 𝐶𝑂𝑆𝜃𝐶 ) − √𝐿2 − 𝑠𝑖𝑛2 𝜃𝑐 (3)
2 4
135
Journal of University of Babylon for Engineering Sciences, Vol. (26), No. (9): 2018.
Continuity Equation
𝜕𝜌
+ 𝑑𝑖𝑣. (𝜌𝑢) = 0 (3)
𝜕𝑡
Momentum Equation
𝐷𝑢 𝜕(−𝜌+𝜏𝑥𝑥 ) 𝜕𝜏𝑦𝑥 𝜕𝜏𝑧𝑥
𝜌 = + + + 𝑆𝑀𝑥 (4)
𝐷𝑇 𝜕𝑥 𝜕𝑦 𝜕𝑧
Energy Equation
𝜕ℎ 𝐷𝑝
𝜌⌊ + 𝑑𝑖𝑣. (ℎ𝑉)⌋ = − + 𝑑𝑖𝑣(𝑘𝑔𝑟𝑑𝑇) + 𝜑 (5)
𝜕𝑡 𝐷𝑡
Momentum
Equations Solution No
Process PISO Convergence
loop Test
Resolve Intermediate Yes
Velocity
Ye No
No
Divergence
Time
Yes
step
Save file
STOP ERROR
End
136
Journal of University of Babylon for Engineering Sciences, Vol. (26), No. (9): 2018.
137
Journal of University of Babylon for Engineering Sciences, Vol. (26), No. (9): 2018.
Pressure at 390 deg. Crank angle Pressure at 630 deg. Crank angle
Pressure at 720 deg. Crank angle Pressure at 990 deg. Crank angle
138
Journal of University of Babylon for Engineering Sciences, Vol. (26), No. (9): 2018.
139
Journal of University of Babylon for Engineering Sciences, Vol. (26), No. (9): 2018.
Temperature at 390 deg. Crank angle Temperature at 420 deg. Crank angle
The investigation of velocity field in cold flow was showed in figure (11)and (12), that the velocity
at 360 CA about 11m/s in chamber and its maximum at the same angle in the end of intake port , due to
partial opening of intake valve (11-a). It is increases in the chamber until reach to maximum value at 450
CA (11-b), then reduces during 2nd and 3th as shown in (11-c, d). With coincidence with exhaust valve
movement at 870 CA where the valve starts to open partially, according to valve profile the velocity inside
the combustion chamber increases. And then the trend is referring to decreasing value of velocity during 4th
stroke.
(a) Velocity of cold flow at 360 CA (b) Velocity of cold flow at 450 CA
(c) Velocity of cold flow at 630 CA (d) Velocity of cold flow at 780 (e) Velocity of cold flow at 990 CA
CA
141
Journal of University of Babylon for Engineering Sciences, Vol. (26), No. (9): 2018.
142
Journal of University of Babylon for Engineering Sciences, Vol. (26), No. (9): 2018.
Figure (15) shows the temperature distribution at different crank angles, the temperature increases
and the maximum value during the flame propagation as shown in the figure below due to combustion.
Then, it reduces when volume of combustion chamber is increased due to expansion during moving the
piston to BDC.
143
Journal of University of Babylon for Engineering Sciences, Vol. (26), No. (9): 2018.
Temperature at 390atCA
Temperature 390 CA TemperatureTemperature
at 720 CA at 720 CA
144
Journal of University of Babylon for Engineering Sciences, Vol. (26), No. (9): 2018.
145
Journal of University of Babylon for Engineering Sciences, Vol. (26), No. (9): 2018.
146
Journal of University of Babylon for Engineering Sciences, Vol. (26), No. (9): 2018.
5. Conclusion
The flow characteristics within SI engine were investigated utilizing CFD ICE CODE. As a result,
the following conclusions were summarized:
The adopted approach of the modeling of combustion simulation via ANSYS ICE CODE with dynamic
mesh technique can be used to develop of internal combustion engine.
The two equations standard k-ε model showed acceptable model to simulate the motion of flow within
engine.
Through this work, the in-cylinder flow characteristics were obtained as expected trends.
The temperature, pressure, and velocity were plotted with respect to the crank angle.
The value of temperature and pressure for firing simulation were bit higher than researches
experimental due the fact this simplified CFD simulation doesn’t include friction losses induced
interactions of engine components.
CONFLICT OF INTERESTS.
- There are no conflicts of interest.
References
[1] HEYWOOD, J.B., Internal Combustion Engine Fundamentals. New York: McGraw- Hill, 1988.
[2] Paulina. S. Kuo. “Cylinder Pressure in a Spark-Ignition Engine a Computational Model “.Engineering
Sciences, J. Undergrad. Sci. 3: pp141-145 fall 1996.
[3] Venkata Suman Movva. “Simulation Of Fluid Flow In Internal Combustion Engines Using Wave
Action Simulation“. M.Sc. thesis.2004.
[4] Ehsan Omaraa “Analysis of Flow Characteristics in Inlet and Exhaust Manifolds of Internal
Combustion Engine“. M.Sc. Thesis.2009.
[5] T.Morauszki, P.Mandli, Z.Horvath, and M.R.Dreyer. “Simulation of Fluid Flow Combustion and Heat
Transfer in Internal Combustion Engine“. Hungarian Journal of Industrial Chemistry Vol.39 (1) pp. 27-
30-2011.
[6] K. M. Ravichandra, D. Manikanta, M. Kotresh. “CFD Simulation of an Engine by Producer of Gas“.
International Journal of Civil Engineering and Technology (IJCIET), Volume 8, Issue 10, October
2017.
[7] Abdul Rahiman, Abdul Razak R. K., Mohammad Samee A. D., Ramis M. K. “CFD Analysis of Flow
Field Development in a Direct Injection Diesel Engine With Different Manifolds”. American Journal of
Fluid Dynamics, 4(3): 102-113, 2014.
[8] ANSYS / Work Bench / V15.0. ICE Guide.
[9] Kakaee, Gharloghi, Foroughfar, and Khanlari. “Thermo-mechanical analysis of an SI Engine Piston
Using Different Boundary Condition Treatments” J. Cent. South Univ. 22: 3817−3829 DOI:
10.1007/s11771-015-2926-7 Springer, 2015.
[10]S. Krishna Addepalli, Mallikarjuna J.M. “Parametric Analysis of 4-Stroke GDI Using
CFD“.Alexandria Engineering Journal.2016.
147
Journal of University of Babylon for Engineering Sciences, Vol. (26), No. (9): 2018.
abc_logo@yahoo.com MahmoodMashkoor@hotmail.com
ANSYS /ICE
SolidWorks
17
. ICE :
148