2015 26 0095 PDF
2015 26 0095 PDF
2015 26 0095 PDF
2015-26-0095
Published 01/14/2015
Bhagyada Dhingra
Sparsh Sharma
Kamalkishore Vora
ARAI
B Ashok
Abstract
Introduction
Downloaded from SAE International by Sparsh Sharma, Friday, January 16, 2015
improve the properties of biodiesel by using additives [17], preheating the fuel to level its viscosity to that of diesel [18] and
adjusting the injection timing [19].
The current research is focused on the measures to be done on engine
geometry so as to not only deliver advantage to any specific fuel. The
most popular techniques employed to improve the combustion
efficiency using HVF (Higher Viscous Fuel due to their higher
viscosity compared to petroleum diesel fuel) are preheating HVF
before injection [8, 20] and changing the injection timing [21].
According to the CI combustion theory, better engine performance
may be achieved with HVF by improving the in-cylinder air-fuel
mixing by increasing the swirl (rotation of air view from top of the
cylinder) and tumble (rotation of air view from front of the cylinder)
of in-cylinder air inside the fuel-injected region. Various methods
have been used to improve the in-cylinder airflow in diesel engines
such as: closing part of the intake port using a butterfly valve
[22,23,24,25] modifying intake port/manifold geometry inserting
shroud around the intake valve and modifying the piston bowl.
Alternatively, in-cylinder airflow may be improved through the use of
a guide vane.
One of the primary approaches to get the fuel entirely oxidized is to
make it mix properly with the air is swirling or commonly known as
turbulence. The more is the swirling of air-fuel mixture, the more
breakdown of viscous fuel particle is. The breakdown of these heavy
fuel particles occurs during the compression/combustion stroke. The
concept of organized air flow is not new and was first found for
carbureted engine using gasoline using gasoline by using rotating
blades placed between carburetor and intake manifold to provide more
air swirl for efficient air-fuel mixing [26]. Afterwards, atomizers and
fuel molecule breaker were patented with different identities [14].
Instead of duplicating what has been done, this research explores the
potential of swirling in CI Engine and designing of a new mechanism
to deliver it to improve the in-cylinder air characteristics to
eventually improve mixing and combustion process to improve the
engine performance. The figure shown below shows the basic shape
of airfoils which has been tested in present work.
(1)
(2)
In this equation p, and SM are the fluid pressure, strain rate and
momentum source, respectively. Since this problem is quite a
complex thus the involvement of Energy equation in the calculations
is inevitable. The rate of change of energy inside the fluid element
which is equal to the combined net flux of heat into the element and
the net rate of work done on the element due to the body and surface
forces comprises the Energy equation which is further mathematically
denoted as follows:
(3)
Where htot and represent the total enthalpy and thermal conductivity
respectively of the fluid.
Downloaded from SAE International by Sparsh Sharma, Friday, January 16, 2015
essential drawings of the intake port, intake runner, intake valve, and
cylinder volume. The dynamic analysis is carried out in ANSYS
Fluent. The CFD simulation is carried out for only cold flow without
combustion. At the boundaries i.e. at the intake as well as exhaust
manifold, constant pressure conditions are given. Logarithmic law of
wall unified with No-slip wall boundary condition is used. Moreover
the simulation is dynamic and time dependent thus turbulence
modeling comes into play.
(4)
(6)
The first term on the RHS is the production of k, the second term ()
is the specific dissipation per unit mass. The last terms describe the
transport of k by molecular and turbulent diffusion. This model is
best suited for flows with high Reynolds number where turbulence is
almost isotropic. The scale elimination procedure in RNG theory
results in a differential equation for turbulent viscosity:
Spring-Based Smoothing
In the spring-based smoothing method, the edges between any two
mesh nodes are idealized as a network of interconnected springs. The
initial spacing's of the edges before any boundary motion constitute
the equilibrium state of the mesh. A displacement at a given boundary
node will generate a force proportional to the displacement along all
the springs connected to the node. Using Hooke's Law, the force on a
mesh node can be written as [45]:
(8)
(7)
where
and
are the displacements of node i and its neighbor j, ni
is the number of neighboring nodes connected to node i, and kij is the
spring constant (or stiffness) between node i and its neighbor j. The
spring constant for the edge connecting nodes i and j is defined as
(9)
Downloaded from SAE International by Sparsh Sharma, Friday, January 16, 2015
The simulation goes projected when the geometry gets the required
set of physical parameters associated with the operational conditions.
In preceding sections, it is decided to go with k- Turbulence
modelling with dynamic meshing supported by spring based
smoothing and dynamic layering.
Boundary Conditions
This section describes the associated boundary conditions with the
suction stroke in CI engine. There is no turbocharging linked with the
in-take stroke. The following figure shows all the attached boundary
conditions to every important component.
The simulation goes projected when the geometry gets the required
set of physical parameters associated with the operational conditions.
In preceding sections, it is decided to go with k- Turbulence
modeling with dynamic meshing supported by spring based
smoothing and dynamic layering.
Figure 3. Dynamic layering at the moving boundary layer (piston-cylinder
wall)
If the cells in layer j are expanding, the cell heights are allowed to
increase until
(10)
Here hmin is the minimum cell height of cell layer j, hideal is the ideal
cell height, and as is the layer split factor.
If the cells in layer j are being compressed, they can be compressed
until
(11)
where ac is the layer collapse factor. When this condition is met, the
compressed layer of cells is merged into the layer of cells above the
compressed layer; i.e., the cells in layer j are merged with those in
layer i.
(12)
Downloaded from SAE International by Sparsh Sharma, Friday, January 16, 2015
Figure 7. Cut plot showing velocity distribution in S1-S4 (clockwise from top
left) configurations
Adding up the results from all the cut plots, it is inferred that S2 and
S4 are the best possible designs for ASD, but to move on for further
iteration process it is suggested to select the best so that advance
design changes can be introduced into the selected design. A thorough
analysis of the plots suggest a draw between S2 and S4 whereas
figure 7 and 8 characterize the drawbacks in S2 thus marking S4 as
the best suited ASD design.
Figure 6. Cut plot showing velocity (left) and vorticity (right) distribution in
baseline configuration
Figure 8. Cut plot showing vorticity distribution in S1-S4 (clockwise from top
left) configurations
Downloaded from SAE International by Sparsh Sharma, Friday, January 16, 2015
Figure.11 shows the vorticity plot for both the configurations. The
swirl intensity is high in L2 configuration at every place. On the walls
where velocity seems to less in previous graph, vorticity plot
evidences the presence of high swirl at those positions. In L1 the
swirl is critically less the left side of the cylinder whereas in L2 the
red and yellow color confirms the high swirl.
From above interpretation, ASD-S4-L2 comes out to be the best and
can be moved on the next the iteration phase.
Till now S4 shape with L2 positioning is confirmed by the robust
computational analysis. This chapter deals with the optimization of
the component for number of blades. A parametric study is performed
with the component maintaining the constant engine geometry
throughout the analysis. Three different ASDs with different number
of blades have been modeled, in SoildWorks 2014, namely B1, B2,
and B3.
Figure 12. Cut plot showing velocity distribution in B1, B2, and B3 (from left
to right)
Figure 10. Cut plot showing velocity distribution in L1 (left) and L2 (right)
configurations
Figure 13. Cut plot showing vorticity distribution in B1, B2, and B3 (from left
to right)
With the most of red color inside the cylinder, configuration with four
blades can be considered as the best in terms of velocity distribution.
The little yellow spot near the left wall shows some static region but
overall it has better performance capabilities.
Figure 11. Cut plot showing vorticity distribution in L1 (left) and L2 (right)
configurations
Now comes to the vorticity plot. This plot also favors the coin to B1.
The high pressure swirl is present around the wall of B1. B2 also has
good turbulence relative to B1. Upon looking B3, there is no sign of
high swirl.
Downloaded from SAE International by Sparsh Sharma, Friday, January 16, 2015
Figure 14. Cut plot showing vorticity distribution in H1, H2, and H3 (from left
to right)
The final plot in the series of this long iteration process is the
vorticity plot accommodating H1, H2, and H3 as shown in figure 15.
Looking at the extreme right, H3 configuration, the vorticity around
the wall is observed which diminishes when moving in to the central
axis. Two dead zones are clearly visible. H2 delivers a surprisingly
high vortical flow with yellow portion motivating the presence of
swirls. Finally coming to H1, i.e. 2.5mm long blades, the vortical
flow is strongest in this configuration with least concentration of dead
fluid zones. The red portion on the right side of the cylinder wall is a
great indication of high velocity swirl.
The final ASD as shown in figure 16 has been manufactured using CNC
The engine simulation model is modelled conferring to the Kirloskar
TAF-1 diesel engine that operated at a constant speed of 1500 rpm is
coupled with a eddy current dynamometer. The technical
specifications of the engine are listed in Table 2. SolidWorks 2014 is
used to model the essential drawings of the intake port, intake runner,
intake valve, and cylinder volume.
The analysis is done for 5 load variations namely 20%, 40%, 60%,
80%, and 100%; with and without ASD.
Figure 15. Cut plot showing vorticity distribution in H1, H2, and H3 (from left
to right)
Downloaded from SAE International by Sparsh Sharma, Friday, January 16, 2015
The engine test is made possible not only because of engine but also
the auxiliaries like electronic display board which shows engine
speed in rpm, data acquisition system recording every bit, fuel tank
holding enough diesel fuel for the test, gas analyzer for analyzing the
emission constituents, and dynamometer. Schematic diagram of the
experimental set up is depicted in figure 18.
Both the graphs show no deviation to each other. At 20% load the
baseline engine run produces 233ppm CO2 whereas run with ASD
produces only 193ppm which is 40ppm less CO2 emissions at minimal
load. At 40% load, ASD engine run emits 50ppm less CO2. This trend
is dominant for all the load values and on an average engine with ASD
emits 45ppm less CO2 than an engine without ASD.
Since the NOx formation is more complex than any other pollutant so
are its measures. At nominal load of 20%, instead of a drop in NOx
emissions, an unexpected rise of 35ppm is observed which could be
estimated due to more mixing inside the cylinder. At 40% load the
emissions get lower in ASD engine whereas the baseline engine keep
following the same trend. It can be simply inferred from the figure 20
that NOx emission increases in engine with ASD which could a
possible outcome of more mixing. So it is advisable for future work if
catalytic converter can be incorporated in the research work for
further development of this new idea.
Downloaded from SAE International by Sparsh Sharma, Friday, January 16, 2015
Figure 25. Cut plot showing velocity comparison between baseline (left) and
ASD-S4-L2-B1-H1
Downloaded from SAE International by Sparsh Sharma, Friday, January 16, 2015
Figure 26. Cut plot showing vorticity comparison between baseline (left) and
ASD-S4-L2-B1-H1
References
1.
2.
3.
4.
Russo, D., et al., State of the art of biofuels from pure plant oil.
Renewable and Sustainable Energy Reviews, 2012. 16(6): p.
4056-4070.
5.
6.
7.
8.
9.
Downloaded from SAE International by Sparsh Sharma, Friday, January 16, 2015
Contact
Bhagyada Dhingra
Fiat-Chrysler Automotive pvt Ltd.
bhagyada.dhingra@fcagroup.com
Sparsh Sharma
Moscow Institute of Physics & Technology
sssparsh14@gmail.com
The Engineering Meetings Board has approved this paper for publication. It has successfully completed SAEs peer review process under the supervision of the session organizer. The process
requires a minimum of three (3) reviews by industry experts.
All rights reserved. No part of this publication may be reproduced, stored in a retrieval system, or transmitted, in any form or by any means, electronic, mechanical, photocopying, recording, or
otherwise, without the prior written permission of SAE International.
Positions and opinions advanced in this paper are those of the author(s) and not necessarily those of SAE International. The author is solely responsible for the content of the paper.
ISSN 0148-7191
http://papers.sae.org/2015-26-0095