Large Eddy Simulation of in-cylinder flow in
Internal Combustion Engines
Araz Banaeizadeh and Farhad A. Jaberi Department of Energy & Visteon Corporation
Objective Grid Generation Solution Method
For every manifold one
To have a high quality grid for a block is required, So Navier-
Using Filtered Compressible Navier-Stokes Equations.
Numerical simulation of in- moving piston, complex cylinder three blocks are added
cylinder flow in MSU/Visteon head and moving valves, a unique 12-
to have 12-block grid Fourth order compact differencing for spatial derivatives.
internal combustion engine multi-block grid system is developed
multi- 31-part initial grid system. For solving the
with a new multi-block LES
multi- generated grids two Third order Runga-
Runga-Kutta for temporal derivatives.
for MSU’s 3-valve laboratory engine.
kind of grid overlap
flow solver. This grid was made out of 31 initial are used along block Multi-
Multi-block method to model complex geometries.
parts.The 31 parts were then merged common boundaries.
to a form a nine block grid system. Gradient Models for sub-
sub-grid scale correlations.
Results (2D Slice Contours) Results (3D Iso-Levels) Results (Manifold
U(x) Velocity Temperature 3D iso-levels of vorticity, temperature and pressure velocity profile)
Piston and intake CA : Crank Angle for CA=90’ and CA=115’ CA : Crank Angle ; Vp : Piston Velocity
valves are 0.02
accelerating Vorticity Temperature Pressure CA= 15’, Vp=1.5 m/s
Radial Distance from Manifoid to Valve Shank (m)
CA= 40’, Vp=4.5 m/s
CA= 65’, Vp=8.0 m/s
CA= 90’, Vp=12. m/s
CA=115’, Vp=13. m/s
0.015 CA=136’, Vp=11. m/s
V(piston)=12 m/s CA=165’, Vp=5.0 m/s
Stroke=105 mm Bore=90 mm
Two Intake Valves D=33 mm Piston is
One Exhaust Valve D=37 mm downward and 0.005
intake valves are CA=115’
Direct Injection accelerating
5.4L 3 Valve Engine Lift Profiles of Intake and Exhaust Valves
upward 20 40 60
Piston Location as a function of Engine Crank Angle
Axial Velocity Component (m/s)
Intake Lift (mm)
at TDC V(piston)=13 m/s
Exhaust Lift (mm) 90.00
P is ton L ocation from BD C (m m )
Valve Lift (mm )
Conclusions & Remarks
at BDC at BDC
0 90 180 270 360 450 540 630 720
0 90 180 270 360 450 540 630 720 downward and CA=245’
Crank Angle (degrees)
Due to vortical flow structure that is generated during intake stroke, large eddy simulation is
Crank Angle (degrees)
intake valves are
Valves Motion and Piston Motion accelerating highly recommended for IC engine flow simulation.
versus crank angle scale
LES is able to accurately predict the evolution of the large scale structure in IC engines.