Modeling and CFD Analysis of Air Flow
Modeling and CFD Analysis of Air Flow
Abstract: In this study an effort was made to develop a flow simulation modeling and
performance prediction for a centrifugal compressor stage of a heavy-duty D. I. diesel engine.
The model is implemented in simulation software in MATLAB language. Additionally, a
contribution of this paper to demonstrate that off-design performance of a centrifugal
compressor stage in a turbocharger system can be accurately simulated using commercial CFD
software, with design software, FLOEFD to generate high quality meshes and model solving.
The vector plots, contour plots and stream line plots are generated for better understanding of
fluid flow through centrifugal compressor stage. Correlation coefficients have been introduced
in the calculation program, in order to be closer to the CFD simulation results. The results
obtained from mathematical computation model were validated with the CFD analysis and
experimental results performed using a test bench for the variation of the performance
parameters such as isentropic efficiency, power input, and total pressure ratio with mass flow
rate, the results are also presented in graphical form. The results reveal that reasonable
agreement between mathematical models, the numerical results obtained from the CFD
simulations and the real measurements; the maximum difference never exceeds 5%. The results
indicate that the developed mathematical computation model can yield better predictions of
performance for a centrifugal compressor stage in a turbocharger system.
Keywords: centrifugal compressor stage, computational fluid dynamics (CFD), heavy duty
diesel engine, experimental work, mathematical computation model.
1. Introduction
In the last decades, charging the engine is an almost diffused strategy towards the performance
rise. This may be archived with a compressor before the intake air box increasing the air (or
mixture) pressure and introducing a higher mass in the cylinders; this is a great advantage,
improving engine power, reducing fuel consumption and decreasing emissions. However,
because of complex geometry of the turbocharger, it cost a lot of time and labor power to adopt
experiment to gain relevant data of performance. So, it becomes more economical to simulate
turbocharger internal flow field and analysis data. Mathematical modeling and Computational
Fluid Dynamics (CFD) is a cost affected tool to provide detailed flow information inside the
*
Research Eng., Mechanical Power Engineering Department, Faculty of Engineering,
Mataria, Helwan University, Cairo, Egypt; [email protected]
†
Dr., Mechanical Power Engineering Department, Faculty of Engineering, Mataria, Helwan
University, Cairo, Egypt; [email protected]
‡
Assoc. Prof., Mechanical Power Engineering Department, Faculty of Engineering, Misr
University for Science & Technology, Cairo, Egypt; [email protected]
§
Prof., Mechanical Power Engineering Department, Faculty of Engineering, Mataria, Helwan
University, Cairo, Egypt; [email protected]
1/15
Paper: ASAT-17-129-PP
2. Methodology
This study consists of three main parts and they are done in a specific order. The first part has
the purpose to develop a mathematical computation model for the flow behavior and
performance prediction from inlet to the exit of a centrifugal compressor stage consisting of all
the components in place. The second part, involves an investigation of numerical flow
simulation and performance characteristics through turbocharger (compressor) stage using
FLOEFD flow simulation as computational fluid dynamics (CFD) software. Simulation of the
off-design performance of the centrifugal compressor stage was considered for six off- speeds
(rpm):60%,70%, 80%, 90%, 100% and 110% of design speed and the various performance
parameters concerning: pressure, temperature distribution and velocity profiles on the blades,
isentropic efficiencies, Power input and pressure ratio are presented. In addition, correlation
coefficients have been introduced in the mathematical calculation program, in order to be closer
to the numerical simulation. A discussion for the second part will be done before continuing to
the third part. In the third and final part mathematical computation model were validated with
numerical results of the simulation obtained from CFD analysis and the experimental results
performed using a test bench. Furthermore, discussion and conclusions can be found.
2/15
Paper: ASAT-17-129-PP
𝑝5
𝜋𝐶 = ⁄𝑝𝑎 (1)
𝑖𝑛
𝐾−1
𝜂𝑎𝑑.𝐶 = 𝑇0 (𝜋𝑐 𝑘 − 1) /(𝑇5 − 𝑇0 ) (2)
𝐿 𝑘−1 𝑘/(𝑘−1)
𝑃5 = 𝑃4 (1 − 𝑅𝑟,𝑠𝑐 ⋅ ) (3)
𝑇 𝑎 5 𝑘
(𝑛𝑑 −1)/𝑛𝑑
𝑇4 = 𝑇2 (𝑃4 /𝑃2 ) (4)
As shown from the equations (5 -8), Where Ga (kg/s) and 𝐺𝑎𝑐𝑜𝑟𝑟𝑒𝑐𝑡𝑒𝑑 (kg/s) are the discharge
air mass flow rate and corrected discharge air mass flow rates respectively at the impeller inlet.
𝑇1 Total, inlet (K) is the total (stagnation) temperature at the impeller inlet, and 𝑃𝑎 𝑖𝑛 total, inlet
(Pa) is the total (stagnation) pressure at the impeller inlet. The model described in this section
is implemented in simulation software in MATLAB language to predict the turbocharger
compressor stage performances. The developed simple model offers a relatively accurate
performance prediction. However not accurate enough compared to experimental results, thus
all performance characteristic magnitudes will be mathematically calculated and presented as a
function of mass flow rate and turbo speed in order to be correlated based on the CFD and
mathematical results.
𝛼𝜑𝑠 𝑙0 𝑁𝑒 𝑔𝑒
Ga = (5)
36×105
𝐺𝑎𝑐𝑜𝑟𝑟𝑒𝑐𝑡𝑒𝑑 = 𝐺𝑎 + δ (6)
1
𝑣𝑎𝑖𝑛 −𝑣12
𝑇1 = 𝑇𝑎𝑖𝑛 + (7)
2𝑐𝑝
𝑃𝑎𝑖𝑛 = 𝑃0 − ∆𝑃𝑖𝑛 (8)
3/15
Paper: ASAT-17-129-PP
4/15
Paper: ASAT-17-129-PP
5/15
Paper: ASAT-17-129-PP
Density increases in thedownstream the compressor due to the increase of the pressure in other
words as shown in Figure.7. It happens this way due to the conversion of the kinetic energy
into a static energy which leads to an increase in both of the pressure and density causing
maximum air density 1.35 kg/𝑚3 downstream the compressor wheel However, at the volute
exit the density intend to decrease due to the effect of engine suction in the air intake manifold.
6/15
Paper: ASAT-17-129-PP
P4(Mpa)
P2(Mpa)
0.14
0.08 0.12
0.13
0.075 0.12
0.11
0.11
0.07 0.1 0.1
60 80 100 60 80 100 60 80 100
RPM % RPM RPM
The flow static pressure across the compressor is lower than calculated mathematically; at
100% of load the pressure at the impeller inlet is lower by 10 % of the calculated due to
compressor suction effect and high-speed entry flow, at the impeller outlet is lower by 8 % and
at the volute pressure is 16 % lower than the mathematical calculated because of the engine
suction pressure in the intake manifold.
355
130
135 305
V1 (M/S)
V2 (M/S)
V4 (M/S)
110
115 255
95 90
205
75 155 70
60 80 100 60 80 100 60 80 100
RPM RPM RPM
Due to the boundary conditions and the pressure reduction the flow velocity across the
compressor is higher than the calculated, following the first law of thermodynamics the velocity
head increased with the reduction of the pressure, at the impeller inlet the flow is higher with
9%, 9.6 % at the impeller outlet, and 9 % at the volute higher velocity than the calculated.
7/15
Paper: ASAT-17-129-PP
1.4 1.55
1.01 1.5
ρ1 (kg/𝑚^3 )
ρ2 (kg/𝑚^3 )
ρ4 (kg/𝑚^3 )
1.3 1.45
0.96 1.4
1.2
1.35
0.91 1.1 1.3
1.25
0.86 1 1.2
60 80 100 60 80 100 60 80 100
RPM RPM RPM
As a result of static head reduction and velocity, the density across the compressor flow passage
is less than the calculated mathematically, Making 9.5% difference at the impeller inlet, 3.8 %
difference at the impeller outlet and 9.8 % difference at the volute discharge.
T2 (K)
T4 (K)
285 350
330 340
280 330
320
320
275 310 310
60 80 100 60 80 100 60 80 100
RPM RPM RPM
CFD temperature results are less than the mathematical predictions as an effect of static pressure
reduction across the whole system, at impeller inlet a 4% difference, at impeller outlet less than
1%, and 2.5% difference between numerical and mathematical calculations.
8/15
Paper: ASAT-17-129-PP
45
0.76 40
35
Power (kw)
0.71
30
η
0.66
25
0.61
20
0.56 15
0.3 0.5 0.7 0.3 0.5 0.7
Ga (kg/m^3) Ga (kg/m^3)
2 0.75
1.9
0.65
1.8
1.7
𝑀 ̇ (kg/s)
0.55
∆P
1.6
1.5 0.45
1.4
1.3 0.35
1.2 0.25
0.3 0.5 0.7 60 80 100
Ga (kg/m^3) RPM
𝐹 𝑁
Since 𝐵𝑃 = ƞ𝑏𝑡ℎ × ƞ𝑉𝑜𝑙 × ƞ𝑀 × 𝜌𝑏.𝑐 × 𝐴 × 𝐶𝑉𝐹 × 𝑉𝑠𝑡 × 60𝑛 × 𝑖 (9)
Thus, by neglecting the changes in the efficiencies we may simplify the pervious equation to
compare engine break power expectations as a function of air density before the cylinder and
Fuel to air ratio
𝐹
𝐵𝑃𝑁 𝜌𝑏.𝑐 𝑁 𝐴𝑁
=𝜌 × 𝐹 (10)
𝐵𝑃𝑀 𝑏.𝑐 𝑀
𝐴𝑀
So, if the desired output break power is to be 360 hp it shall be 241 hp instead.
Correcting the inlet pressure will in turn change the inlet density.
𝑃4 = 𝑃𝑎,𝑖𝑛 𝜋𝑐 (12)
𝑃4𝑐 = ϑ𝑃4
Obtaining a correct value of the discharge pressure will in turn correct the values of the density
and temperature
where
∝φs l0 Ne ge
𝐺𝑎 = (13)
36×105
is the discharge mass flow rate.
10/15
Paper: ASAT-17-129-PP
Using these correction factors within the Demidov, V and kolchin. A mathematical model we
will obtain correct values of mass flow rate ــ3 % error ــand pressures, temperatures, densities,
and in turn correct values of efficiency, power required.
7. Experimental Setup
The turbocharger chosen to perform this study with the reference Holset HX 80 was fitted on a
heavy-duty D. I. diesel engine; (360 hp@2000 rpm and 1300 Nm@ 1450 rpm), featuring a
mechanical injection system, more details about the main features of the engine and the
turbocharger are listed in Table 1. The engine was instrumented on a high-dynamic test bed
which was installed at the Faculty of Engineering R&D Center in Helwan University, Cairo,
Egypt. The engine was fully equipped for measurements of all standard operating parameters
(torque, speed, air mass flow), but also of the static pressure and temperature in the most
relevant locations along the intake and exhaust system, as well as of high frequency pressure
data at the same locations and inside each cylinder. Additionally, Cooled pressure transducers
were used in order to obtain an accurate measurement of the gas pressure conditions at the
compressor inlet and outlet. In addition, the compressor rotating speed was also recorded in
order to complete the information related to the turbocharger. Tests were carried out at different
engine speeds and loads corresponding to the whole operating range of the engine. In every
performed test, pressures were recorded in different points and the corresponding air mass flow
and fuel flow were measured. The compressor discharge and turbine inlet headers were fitted
with insulation to minimize heat loss to the surroundings, and maintain accuracy of the high-
temperature measurements. A sketch of the engine layout and locations of the transducers used
during the measurement campaign is shown in Fig.15. All instruments are connected to a Data
Acquisition Card (DAQ) to continually monitor the turbocharger under test. In addition to
monitoring temperatures and pressures of the turbine and compressor conditions, the DAQ
monitors bearing and jacket water conditions for stability and safety. The complete performance
profile of a turbocharger is determined from temperature and pressure instruments, and the
related thermodynamic and fluid mechanic properties. A DAQ is programmed to compute flow
and performance calculations of the operating turbocharger and record the performance data
11/15
Paper: ASAT-17-129-PP
for test reporting, in accordance with the Test Standard for Testing Large-Bore Engine
Turbochargers Chapman et al. 2005[12]. For every test condition, the air mass flow, fuel mass
flow, and compressor inlet temperature and compressor outlet pressure were measured by the
sensors of the installation. As the different tests were performed under steady state operations,
the power consumed by the compressor is equal to the power developed by the turbine including
turbocharger mechanical losses. So, the adapted power was obtained for the compressor power
as:
𝑘 (𝑘−1)⁄𝑘
[𝑘−1 𝑅𝑎 𝑇𝑎𝑖𝑛 𝑄𝐸𝑋𝑃 (𝜋𝑐 − 1)]
𝑃𝑐𝑜𝑚𝑝 = ⁄ (14)
1000𝜂𝑎𝑑.𝑐
Small differences can be observed between measured and interpolated data. On the other hand,
inherent sensor uncertainties exist in both type of installation. As the parameters calculated to
describe the compressor performance are combinations of various measurements, results can
be slightly influenced. On the other hand, some heat transfer effects [10], especially in the low
speed and low load region, and differences between compressor behavior under steady or under
pulsating flow conditions [12], affect the compressor efficiency measurement.
12/15
Paper: ASAT-17-129-PP
The results obtained from the corrected mathematical computation model, were validated with
CFD analysis and the experimental results determined on a test bed for performance parameters
such as polytrophic efficiency, power input, and total pressure ratio versus mas flow rate.These
results will be presented in this section as shown in Figures 15 to 18 for turbocharger design
speed of 39000 rpm(Engine Full load).The first performance measure investigated was the total
to static pressure ratio versus the corrected mass flow rate. This is a ratio of the static pressure
ratio at the exit divided by the total pressure at the inlet respectively. As shown in Figure15
total pressure ratio of a centrifugal compressor stage as estimated by mathematical computation
code almost complies, with negligible variation of 0.3%, 0.9% with CFD tools and experimental
results respectively. Figure16 shows mathematical computation code on the polytrophic
efficiency of a centrifugal compressor predicted both CFD results and the experimental results
closely with a variation of 0.3% and 1.4% respectively. Similarly Figure.17 shows power input
to a centrifugal compressor stage is predicted by mathematical computation code which
compares closely with CFD and experimental results with a variation of 1.7% and 4.4. It is clear
from these figures that and mathematical computation results are closely compared to the CFD
numerical simulation results and the obtained experimental values and the difference doesn't
exceed 5% or less for all the calculation points. This means that the flow in the compressor
stage can be calculated with accuracy using both CFD and mathematical calculations methods.
Once again, a good overall agreement is obtained. As justified before some uncertainties can
be observed but they are always included in a narrow confidence interval. This comparison
confirms that the representation defined is valid to predict the turbine performance in the whole
turbocharger operating conditions.
13/15
Paper: ASAT-17-129-PP
2 0.8
1.9 0.75
1.8
1.7 0.7
∆P
1.6 0.65
η
1.5 0.6
1.4
1.3 0.55
1.2 0.5
0.3 0.4 0.5 0.6 0.3 0.4 0.5 0.6
Ga (kg/s) Ga(kg/s)
45
40
Power (kw)
35
30
25
20
0.3 0.35 0.4 0.45 0.5 0.55 0.6
Ga(kg/s)
9. Conclusion
Mathematical or physical molding for engineering designing is the easiest and fastest method
of designing However, it is not accurate enough for designing and matching, on the other hand
numerical molding is more accurate giving a closer vision to the designer but still it's a time-
consuming process and any small change in the model will requires more time to adopt with it.
Having the accuracy of numerical analysis in a Mathematical Code will guarantee fast and easy
designing and matching method for any manufacture.
This paper presents an effort to model the flow from inlet to the exit of a turbocharger
compressor stage consisting of all the components in place and performance prediction by
providing mathematical computation model and numerical analysis using CFD tools and these
were verified by experimental work. Using of the computational fluid dynamics (CFD) gave a
better understanding of the behavior of flow through turbocharger compressor stage and how it
impacts the turbocharger efficiency and how the turbocharger compressor perform combined
with the diesel engines.
1. A good agreement was achieved between the computational fluid dynamics (CFD),
calculated and the experimental results, on the isentropic efficiency of a centrifugal compressor
stage with maximum variation of 1.4%. Similarly, power input to the centrifugal compressor
stage predicted by CFD results which compares closely with experimental and computed results
with a variation of 1.7% and 4.4% respectively. As well as total pressure ratio of a centrifugal
compressor stage as estimated by CFD tools almost complies, with negligible variation of 0.3%,
and 0.9% with experimental and computed results respectively.
14/15
Paper: ASAT-17-129-PP
2. An agreement between correlated mathematical model, CFD computations and the real
measurements, is obtained, in addition the capability of the correlated mathematical model to
predict satisfactorily the performance of a centrifugal compressor stage in a turbocharger
system with maximum error 5% for most of the performance range.
3. The comparison between experimental and simulated results shows a good correlation. This
result is to a large extent to be expected due to the use of adaptation coefficients, the
establishment of such adaptation coefficients are however a delicate process, firstly because
they must represent the real physical phenomena and secondly due to their independence. The
precision obtained could be considered excessive but is quite necessary in the setting up of the
complete turbocharged diesel engine simulation.
10. References
[1] Prasad, V.V.; Kumar, L.M. and Reddy, B.M., “Centrifugal Compressor Fluid Flow
Analysis Using CFD”, Science Insights, Aug 2011.
[2] Bhardwaj, S. and Buke, Y., “Computational Fluid Dynamics Analysis of a Turbocharger
System”, International Journal of scientific research, May 2014, ISSN No 2277 - 8179
[3] Kolchin A., and Demidov V., “Design of Automotive Engines” MIR, Publishers, Moscow,
1984.
[4] Abdelmadjid, C. Ameur M., Boumeddane B., “CFD analysis of the volute geometry
effect on the turbulent air flow through the turbocharger compressor”, ELSEVIER,
Energy Procedia, 36, (2013), 746 – 755,
[5] Harley, P. Spence, S. Early, J., “Inlet recirculation in automotive turbocharger centrifugal
compressor,”
[6] Stephanie Berqvist, “Prediction of Turbo Compressor Maps using CFD”, Chalmers
University of Technology (2014)
[7] Heywood, J. B. “Internal Combustion Engine Fundamentals” McGraw-Hill, 1988.
[8] A. J. Mohamed, K. Gueorgui, and S. A. Mohamed, “Influence of intake manifold design
on in-8CFD analyses and experimental investigations”, Energy 36, 2701-2715, 2011.
[9] Chapman, K. S., Keshavarz, A., Shultz, J., Sengupta, J., 2005, “Test Standard for Testing
Large- Bore Engine Turbochargers,” Gas Research Institute, Report No. GRI-05/0026,
Des Plaines, IL.
[10] S. Shaaban, J. Seume, “Analysis of turbocharger non-adiabatic performance”, in:
C647/027, Proceedings of 8th International Conference on Turbochargers and
Turbocharging, London, 2006.
[11] Ding M.Y, Groth.C, Kacker.S and Roberts.D, “CFD Analysis of Off-design Centrifugal
Compressor Operation and Performance”. International ansys Conference (2006)
[12] P. Le Sausse, P. Fabrie1, D. Arnou, and F. Clunet, “CFD comparison with centrifugal
compressor measurements on a wide operating range”,EPJ Web Conference 45 , 01059
(2013).
[13] Nakhjiri, M, Pelz, P and Matyschok, B , “Physical Modeling of Automotive Turbocharger
Compressor: Analytical Approach and Validation”, SAE International, 2011-01-
2214(2011).
[14] Galindo.J, Fajardo.P , Navarro.R, Garca-Cuevas.L.M. “Characterization of a radial
turbocharger turbine in pulsating flow by means of CFD and its application to engine
modeling”, Applied Energy 103 (2013) 116–1
15/15