1. Introduction
The increasing trend of pollutant emissions reduction has forced engine manufacturers to focus their interest on hybrid and electric vehicles. The advantages of such vehicles are their relatively low fuel consumption and the possibility to recover energy. However, that does not mean that internal combustion engines (ICEs) should be eliminated from commercial use. One of the most important advantages of the ICE over electrically driven vehicles is the range. The limited battery charge of electric vehicles greatly constrains their maximal driving range. There are a lot of places in the ICE that can be improved to recover more energy from the exhaust gases. The hot exhaust gases contain a large portion of the energy that could be recovered with exhaust after-treatment components such as turbochargers which are driven by the hot exhaust gases from the cylinder. The first turbochargers had a single-stage rotor. However, the need for further energy recovery led to the implementation of multi-stage turbochargers. This has led to a greater energy recovery but on the other hand, also made the exhaust systems more complicated. Nowadays, two-stage turbocharging systems are used in ICEs. In such a system, the recovered energy is used to boost the air pressure at the intake of the engine. An analytical approach based on the power balance between the turbine and compressor with minimum energy consumption offers a reliable way to optimize the two-stage turbocharging system [
1]. Key parameters such as the total expansion ratio and the total pressure ratio are often used with an analytical predesigned model [
2,
3]. The analytical approach of a two-stage turbocharged engine with three bypass valves was also studied at different altitudes [
4]. The results showed that the cylinder combustion pressure decreased with the increase of the altitude. The two-stage turbocharged system was also analyzed during the steady-state and transient engine response regimes [
5]. The results showed an improved durability of the two-stage turbocharging system with the reduction of the break-specific fuel consumption (BSFC) at low engine speeds. The analytical design of the Miller-cycle, two-stage turbocharging system leads to the reduction of the in-cylinder combustion temperature [
6]. The 1-D simulations of series two-stage turbocharging systems allowed us to investigate the control of the pressure inside the exhaust system by the opening and closing of the variable guide vanes [
7]. Researchers have shown that with the use of the 1-D simulations the displacement of the spark ignition V8, a two-stage turbocharged engine’s pressure can be lowered by 20% without reducing the brake torque [
8]. Such simulations used a predictive knock model. The minimum BSFC criterion was also used with the 1-D simulations of the two-stage turbocharging system [
9]. The results were very optimistic, indicating an increase in the brake power and a reduction in fuel consumption. The 1-D simulation of the engine performance, with the two-stage turbocharging system at different altitudes, showed a greater increase in the brake power and torque when compared to the single-stage turbocharger [
10,
11]. Such simulations neglect the pressure losses inside ducts caused by the vortices, however, the 1-D calculations are in good agreement with the experimental data and require less computational time. The more complex 3-D calculations require a lot of computational resources. To reduce computational time, the geometry of the problem can be simplified. Such calculations allow investigating a fully developed turbulent flow with the Reynolds-averaged Navier-Stokes solver [
12]. Some computational solvers for multistage turbomachinery issues implement discrete turbulence models with coupled, mixing-plane boundary conditions [
13]. Such an approach favors the simulations under steady-state conditions. A coupled 1-D and 3-D calculations are also performed [
14]. The 1-D calculations are performed in order to establish inlet boundary conditions for the 3-D numerical model. Both 1-D and 3-D solvers can be run simultaneously or independently.
The turbocompounding system is a type of multi-stage turbine system and one of the approaches to recover energy from exhaust gases. In such a system, an auxiliary turbine is installed behind the turbocharger. The exhaust gases drive the turbocharger and the auxiliary turbine. The energy from the exhaust gases can be recovered mechanically or electrically [
15,
16]. The additional electrical energy generated by the electric turbocompounder can be used to supply the vehicle’s electronic units or batteries. Turbocompounding systems are mostly used in heavy-duty engines, where the reduction of fuel consumption plays an important part. Such systems are often assisted by exhaust recirculation valves which reduce the NOx concentration [
17]. The mechanical recovery system consists of the auxiliary turbine connected with the crankshaft via gear train and fluid coupling. The electrical recovery is based on the auxiliary turbine that drives the alternating current (AC) generator. The generated electric energy can be used to supply the battery.
A two-stage turbocharging system further increases the energy recovery of the exhaust gases. However, such systems are complex and require precise control over the position of the variable nozzle vanes. Another interesting fact is that the two-stage turbocharging system increases the backpressure in the exhaust duct thus leading to an increase in the pumping losses. This might affect the in-cylinder combustion process. The turbochargers of the two-stage system operate independently from each other. That is why the rotational speed of each stage is different. Also, very few works have dealt with the 3-D computational approach of such a system.
This paper presents the numerical 3-D modelling results of a six-inlet, two-stage radial inflow turbine under pulse-flow conditions. The novelty of this work is the description of a two-stage model with both rotors placed on the same shaft, thus operating with the same rotational speed. Also, the separation of the exhaust pulses was investigated by implementing a six-inlet housing. Such a two-stage system will be used as a turbocompound unit in the 2-stroke, 6-cylinder opposed piston (OP) engine which is being designed at the Warsaw University of Technology. The two-stage system will be connected mechanically to the engine crankshaft via a gear train and fluid coupling.
The paper is organized as follows.
Section 2 presents the theoretical approach for the estimation of the key turbine operating parameters. A detailed description of the computational geometry, mesh, and simulation setup is provided in
Section 3.
Section 4 presents the results of the model validation against the experimental data.
Section 5 presents the results with an in-depth discussion. The conceptual adaptation scheme of the two-stage turbine system for the OP engine is described in
Section 6. The conclusions are given in
Section 7.
2. Theory
The typical hot section of a turbocharger consists of the volute and turbine wheel. The volute might be vaneless or equipped with the nozzle guide vanes, however, most of the volutes are vaneless as they are cheaper, smaller, less complex, and easier to manufacture. However, for research purposes, the nozzle guide vanes were considered to accelerate hot exhaust gases before the rotor wheel.
Due to the rotational movement of the rotor wheel, the correlation between velocities in the turbine stage is represented by Equation (1):
where
is the relative velocity,
is the absolute velocity and
is the linear velocity at the tip of the rotor blades.
Figure 1 shows the turbine configurations with the following Sections: inlet, behind nozzle vanes, behind first-stage rotor wheel, behind second-stage nozzle vanes, and behind the second-stage rotor wheel which corresponds to the following indexes: 0–4.
The pressure ratio at nozzle vanes of first and second-stage is calculated using Equations (2) and (3) [
18]:
where
and
are the pressures behind the first and second-stage nozzle vanes,
and
are the total pressures before the first and second-stage nozzle vanes,
and
are the total temperatures before the first and second-stage nozzle vanes,
and
are the absolute velocities after nozzle vanes for both stages,
is the vane flow coefficient,
is the exhaust specific heat ratio,
is the exhaust gas constant.
The temperature behind the nozzle vane at both stages can be calculated using Equations (4) and (5):
where
and
are the temperatures behind the first and second-stage nozzle vane.
Turbine isentropic enthalpy drop at the first and second-stage can be calculated using the following equation:
where
is the theoretical velocity during isentropic expansion from pressure
to pressure
for the first-stage turbine and from the
to
for the second-stage turbine. The reaction ratio
for both stages indicates rotor enthalpy drop compared to total enthalpy drop and it is calculated as:
where
is the rotor enthalpy drop. Indexes 1 and 2 represent the first and second-stage. The expansion work in rotor blades, as well as expansion efficiency for both stages, can be calculated with the equations:
Also, the total–static efficiency of both stages can be obtained using equations:
where
and
are total enthalpy at the inlet to the first and second-stage,
,
,
and
are total and isentropic enthalpies at the outlet of the first and second-stage.
The blade tip speed of both stages can be calculated using the formula:
where
outer diameter of the turbine wheel. Index 1 stands for the first-stage rotor and 2 stands for the second-stage rotor.
3. Computational Approach
Figure 2a presents the B&W K44 turbine wheel which was used to create the 3-D model for the computational domain. The presented turbine wheel was scanned using the AKON smartSCAN 3D HE scanner and turntable. A total number of 60 scans per 360 degrees was done to scan as many surfaces as possible. However, the inter-blade zones were difficult to process, especially close to the hub. In those areas, computer-aided design software was used to create accurate surfaces.
Table 1 shows the general parameters of the turbine wheel.
Figure 3a–c shows the set-up of the computational domain which consisted of nine domains: the exhaust pipes, the first-stage nozzle vanes, the first-stage turbine rotor, the tip clearance gap, the diffuser, the inter-stage pipes, the second-stage nozzle vanes (variable), the second-stage turbine rotor and the outlet. The same rotor geometry was used for the first and second-stage. The rotors were installed on a single shaft, thus the rotational velocity of the rotors was the same.
As can be seen from
Figure 3a the six exhaust pipes are equally placed on the circumference of the first-stage rotor. The number of exhaust pipes corresponded to the number of engine cylinders. Each pipe transported the exhaust gases from the cylinder to the rotor. The hydraulic diameter of each exhaust pipe was 80 mm and the length was 1500 mm. The nozzle guide vanes were used to accelerate the exhaust gases before entering the first-stage rotor. The trailing edge (TE) angle of the nozzle vane was set to
. To investigate the leakage between the exhaust pipes during the pulse-flow, a small gap between the rotor blade tip and shroud was modeled and shown in
Figure 3c. The height of the gap was set to 3 mm. To reduce the pressure losses between the two stages, the length and cross-Sectional area of the diffuser were iteratively adjusted. The inter-stage pipes consisted of the upstream straight segment, the elbow, and the downstream straight segment. The second-stage nozzle vanes consisted of 19 vanes and were fixed at the angle of
. The outlet domain was extended up to seven diameters of the second-stage rotor outlet.
A structured hexahedral mesh was generated on every domain except for the turbine rotor of the first and second stages. The ICEM CFD software was used for mesh generation. The rotor domains contained tetrahedral mesh which was later transformed into the polyhedral mesh in FLUENT. The mesh independence study of the rotor domains was performed before the beginning of the calculations.
Figure 4 shows the plot of the volume-averaged total pressure values during the single revolution of the crankshaft for the different number of grid densities.
In order to inspect the mesh independence study for the transient issue, the total pressure values were investigated during the single revolution of the crankshaft. It was found that for the exhaust pipes the volume-average pressure plot remained almost unchanged for the number of the elements of 7200. For the first and second-stage rotor, the number of 1.34 million cells provided suitable pressure changes. The pressure changes inside the second-stage nozzle vanes remained unchanged with the number of 711,000 cells.
Table 2 shows the number of cells for the rest of the domains.
As a result, the complete computational domain contained 2.1 million cells.
Figure 5 shows the polyhedral mesh generated on the rotor domains.
The ANSYS FLUENT software was used for transient simulations. For this reason, the sliding mesh method was used to capture the instantaneous parameters during the single time step. For the inlet of each exhaust pipe, the mass flow inlet boundary condition was set. The mass flow rate profile at the inlet to each exhaust pipe is shown in
Figure 6.
The mass flow rate profile was obtained from the 1-D, transient simulations of 6-cylinder, 1500 kW, 2-stroke engine. The rapid increase of the mass flow rate from 0 to 1 kg/s represented the beginning of the pulse-flow during the opening of the cylinder exhaust valve. The interval between the exhaust pulses was
which was adequate for the firing interval of the 6-cylinder, 2-stroke engine which parameters can be seen in
Table 3.
Table 4 shows the parameters of the boundary conditions.
To accurately capture the transient changes inside the whole exhaust system, a time step of
was chosen. Such a time step was chosen due to the strongly warring cross-section and high velocity of the exhaust gases. Based on the previous simulations, such a small time step leads to stable calculations with the pressure-velocity coupling scheme. Each case was calculated using the single-equational Spalart-Allmaras (SA) turbulence model [
19]. The SA turbulence model is based on the transport equation that is solved for a modified turbulent kinematic viscosity
. The modified turbulent viscosity is assumed to be linear in near-wall region which greatly reduces computational time. Most of the turbomachinery simulations are done with k-ω shear stress transport or k-ε turbulence models. However, the SA model is also capable of capturing near-wall shockwave [
20,
21,
22]. The SA model can also be extended for subsonic and transonic flows which also makes it suitable for capturing wake inside the stator domain [
23]. Due to the linearization of the turbulent kinematic viscosity, the computational time is reduced [
24]. The SA model tends to be more stable during calculations with the unstructured grid. With such a small time step size, the SA turbulence model greatly reduces the computational time.
The calculations were carried out at the three different turbines rotational speed.
Table 5 shows the defined operating points of both rotors and the blade tip speed, which was calculated with Equation (13).
The calculations were carried for three revolutions of the engine crankshaft. The reason behind this was to achieve the same pressure variations inside each exhaust pipe during the single revolution of the crankshaft.
Table 6 shows the parameters of the node device that was used for computing purposes.
With such computational restraints and such a low time step, the time needed for the calculation of the single revolution of the crankshaft was about 1.5 weeks.
Figure 7 shows the mass flow average pressure variations at the first-stage turbine inlet during the 3 revolutions of the crankshaft.
Figure 7 shows that for every turbine speed the pressure variations for each exhaust pipe were almost similar during the 3rd revolution of the crankshaft except for the 40,000 rpm case. The shapes of the pressure plot for each exhaust pipe were different in this case. For the 40,000 rpm case, the maximal pressure values during the 3rd revolution of the crankshaft changed from 485,000 Pa for the exhaust pipe 2 to 502,500 Pa for the exhaust pipe 5. However, such a difference was within the acceptable range. The minor changes can be observed between the 3rd and 2nd revolution of the crankshaft. However, based on the author’s previous experiences, the results for the 4th revolution of the crankshaft were not presented as further changes did not exceed 1%. That is why the results were presented during the 3rd revolution of the crankshaft.
5. Results and Discussion
Figure 9 shows the plots of the area-averaged mass flow rate that was monitored in the Section 1-1 for three different turbine speeds. The plots were captured during the 3rd revolution of the crankshaft. For every turbine speed, the shapes of the pressure plots were almost similar and shifted by
crank angle degree (CAD). From
Figure 8 one can see that the shape of the pressure plot of the exhaust pipe 6 during the exhaust phase differs from the rest of the pressure plots. This pressure plot comes from the 2nd revolution of the crankshaft because of the
CAD shift between the exhaust pulses. The shapes of the pressure plots of 40,000 rpm case were different when compared to the pressure plots of 50,000 rpm and 60,000 rpm cases. The dashed rectangle shows the deteriorated mass flow rate at the end of the pulse for every exhaust pipe. Such a phenomenon is shown in
Figure 9 in much greater detail.
Figure 10 shows the mass flow rate and the total pressure of exhaust pipes 1 and 2 during the 3rd revolution of the crankshaft at three different turbine speeds.
From
Figure 10 it can be seen that for 50,000 rpm and 60,000 rpm case the mass flow rate of the exhaust pipe 2 decreased between the
and
CAD and for the
case, between the
and
CAD. This was associated with the increase of the mass flow rate inside the exhaust pipe 1 due to the opening of the exhaust valve. The lowest value of the pipe 2 mass flow rate was
kg/s at
CAD for
,
kg/s at
CAD for
and
kg/s at
CAD for
The pipe 2 exhaust valve opened at the
CAD for every case, this increased the mass flow rate. The interaction between the pipe 1 and pipe 2 mass flow rate after the opening of the pipe 1 exhaust valve was caused by the leakage between those pipes. As it can be seen from
Figure 10a–c, between the opening of the pipe 1 exhaust valve and opening of the pipe 2 exhaust valve, first the exhaust pipe 2 operated in lower pressure values due to the end of the exhaust pulse, while the pressure inside the exhaust pipe 1 was at its peak values due to the opening of the exhaust valve. Although, the height of the tip clearance gap was only 3 mm, part of the exhaust gases passed from exhaust pipe 1 to exhaust pipe 2. Also
Figure 10d–f show the velocity vectors at the state represented by the points A–C on the mass flow rate plots. The figures clearly show the leakage between the adjacent exhaust pipes with velocity vectors entering exhaust pipe 2 in the upstream direction. The lowest leakage was observed at point C for
while the highest leakage was observed at point A for
. At the
rotor caused the lowest acceleration of the exhaust gases. That is why a larger portion of the exhaust gases passed through the tip clearance gap. However, no negative values of the pipe 2 mass flow rate were observed. This means that no backflow of the exhaust gases from pipe 1 to pipe 2 occurred.
Figure 11 shows the pressure variations in the Section 1-1 during the 3rd revolution of the crankshaft for three different turbine speeds. For each turbine speed, the shapes of the pressure plots of each exhaust pipe were almost similar apart from exhaust pipe 6 which came from the 2nd revolution of the crankshaft. The pressure plots of each exhaust pipe were shifted by
CAD which is associated with the inlet boundary conditions shown in
Figure 6.
Figure 11d–f show the pressure contours and represent the turbine state at lines A–C respectively. By looking at the pressure contours, it can be seen that the highest pressure occurred inside the exhaust pipe 5.
Unlike the ordinary scroll turbine, the rotor was not supplied on the whole circumference during a single pulse. The rotor was supplied with exhaust gases only on the
of its circumference. However, due to the short intervals between pulses, it can be seen that the rotor was constantly supplied by the exhaust gases. The intersection of the pressure plots at the Section 1-1 during the 3rd revolution of the crankshaft for three turbine speeds is shown in
Figure 12.
From
Figure 12 it can be seen that during the 3rd revolution of the crankshaft the pressure was changing in a quasi-constant manner for every turbine speed. The average pressure values were
Pa,
Pa, and
Pa for
,
and
case, respectively. Such differences between the pressure values for 40,000 rpm, 50,000 rpm and 60,000 rpm were caused by the increasing pressure ratio.
The pressure variations behind the first and second-stage rotor are shown in
Figure 13a,b. One can see that the pressure at the outlet of the first-stage rotor was quasi-constant with minor fluctuations during the 3rd revolution of the crankshaft. The pressure at the Section 2-2 increased with the increase of the turbine rotational speed which was caused by the increase of the inlet pressure as shown in
Figure 12. The pressure monitored at the Section 4-4 reached lower values with the increase of the turbine rotational speed. However, comparing the pressure plots from
Figure 13a,b, one can assume that the pressure ratio of the second-stage rotor should further increase with the increase of the turbine rotational speed.
To investigate the pressure losses caused by the inter-stage pipes, the pressure values were monitored at the Sections 2.1-2.1 and 2.2-2.2. The results are shown in
Figure 14 during the 3rd revolution of the crankshaft for three different turbine speeds. The instantaneous pressure losses were marked as percentage loss of the pressure monitored at the Section 2.1-2.1. The highest variations of the pressure losses can be seen in
Figure 14a for 40,000 rpm case. The pressure losses, in this case, reached up to 12%. This might be caused by the lower acceleration of the exhaust gasses behind the first-stage rotor and vortices occurring inside the inter-stage pipes. For the turbine speed of 50,000 rpm and 60,000 rpm case, the pressure losses were the lowest. In those cases, the highest value of the instantaneous pressure losses reached up to 6%.
Figure 15 presents the plots of the pressure ratio for the first and second-stage rotor at three different turbine speeds during the 3rd revolution of the crankshaft. As it was mentioned previously, the pressure ratio of the first-stage rotor increased with the increase of the turbine speed. The same phenomenon was observed for the pressure ratio of the second-stage rotor. Higher changes in the pressure ratio values were observed at the first-stage rotor which was caused by the rapid pressure changes at the rotor inlet. On the other hand, changes in the pressure ratio for the second-stage rotor were lower. This was mainly a consequence of the quasi-constant character of the pressure behind the first and second-stage rotor. For the second-stage rotor, the pressure ratio fluctuations tend to be lower with the increase of the turbine speed. Such a phenomenon allows the second-stage rotor to operate in almost steady-stage conditions.
The temperatures at the outlet from both rotors are shown in
Figure 16. From
Figure 16 one can see that the temperature behind the first and second-stage rotor became quasi-constant during the 3rd revolution of the crankshaft. It was found that, with the increase of the turbine speed, the temperature at the outlet increased. The reason for such a phenomenon is not known. A great difference between the temperatures at the Sections 2-2 and 4-4 can be observed which resulted in the recovery of the energy from the exhaust gases on the first and second-stage rotor.
Figure 16c–e show the contours of the total temperature behind the first-stage rotor at state A for the 40,000 rpm, 50,000 rpm and 60,000 rpm case respectively. On the other hand,
Figure 16f–h show the temperature contours behind the second-stage rotor at state B. It can be seen that the temperature drops behind the second-stage rotor. However, minor temperature changes can be seen behind the first stage rotor.
In
Figure 17 one can see the variation of the total-static efficiency of the first and second-stage rotor during the 3rd revolution of the crankshaft for the three different turbine speeds. The efficiency during the 3rd revolution of the crankshaft was oscillating in a narrow range for every turbine speed. With the increase of the turbine speed, the efficiency of the first-stage rotor decreased. It was mainly caused by the presence of the second-stage rotor which increases the pressure behind the first-stage rotor and thus, lowers its pressure ratio. For the first-stage rotor, the highest efficiency was 0.89 for 40,000 rpm case. A completely different situation was observed for the efficiency values of the second-stage rotor. The efficiency values increased with the increase of the turbine speed. The highest efficiency value was 0.62 for the 60,000 rpm case. The high-efficiency values of the first-stage rotor were caused by the considerably higher inlet pressure and higher enthalpy drop. The efficiency of the complete two-stage system can be calculated using the formula below:
Figure 18 shows the total-static efficiency of the complete two-stage system compared to the efficiency of the single-stage system with the same turbine model which was previously presented by the authors in [
25].
The efficiency values of the two-stage system were higher than those in the single-stage system [
25]. The average efficiency values of the two-stage system were 0.491 for the 40,000 rpm, 0.496 for the 50,000 rpm and 0.504 for the 60,000 rpm case. For the single-stage system, the efficiency values were 0.442, 0.475 and 0.484 for the 40,000 rpm, 50,000 rpm and 60,000 rpm case, respectively. Such differences lead to the conclusion that the two-stage system is more efficient than the single-stage system at the low turbine speed, thus at the low-load operating points of the engine however, the overall efficiency of the two-stage system is relatively low when compared to the modern scroll turbochargers.
7. Conclusions
The paper deals with the unsteady simulation of the two-stage, radial turbine under pulse-flow conditions. The novelty of this research was the development of a transient, 3-D simulation of the two-stage turbine with two rotors placed on a single shaft. The separation of the exhaust gases at the inlet of the first-stage rotor was also investigated. The simulations were carried out for the pulse-flow conditions at the inlet. The turbine wheel was scanned with the use of the 3-D scanner to create a virtual model. The same turbine model was used for the first and second-stage. The computational domain was generated using 3-D CAD software. The simulations were carried out using Ansys Fluent software with a small time step. The simulations lasted until the 3rd revolution of the crankshaft. Three different turbine speeds were investigated. The two-stage system was capable of recovering energy from the exhaust gasses. The temperature between the inlet to the first-stage rotor and behind the second-stage rotor dropped by about . Also, the total pressure values changed from about 450,000 Pa at the inlet of the first-stage rotor to 120,000 Pa behind the second-stage rotor. It was also found that despite the pulsatile operation of the two-stage system, the first-stage rotor was constantly supplied with exhaust gases. The pressure at the inlet to the first-stage rotor and behind the second-stage rotor was changing in a quasi-constant manner. The pressure changes decreased with the increase of the turbine speed. It was also found that the presence of the second-stage rotor limited the pressure behind the first-stage rotor and thus, limited its pressure ratio. That is why the total-static efficiency of the first-stage rotor decreased with the increase of the turbine speed. A great efficiency difference between the two-stage system and single-stage system at the low turbine speed was observed. This showed that such a system might be efficient at low-loads of the ICE. However, at the higher turbine speed, the efficiency difference between the two-stage system and single-stage system was lower and reached up to 2%. On the other hand, the two-stage system with newly designed housing provided proper separation between exhaust pulses. A relatively low leakage between the adjacent exhaust pipes was observed during the simulation. Also, no negative values of the mass flow rate were observed which meant that no backflow occurred during the simulation. The low total-static efficiency values of the two-stage system give room for improvement. However, for multicylinder engines, more than one scroll turbocharger needs to be used. Because of that, the exhaust system is more complex. On the other hand, with such a two-stage turbocharging system, the exhaust system would be easier to design. Such a system satisfies the demands of both the variable turbine geometry turbocharger and scroll turbocharger. The two-stage turbocharging system is a very promising device with its reduced leakage between the exhaust pipes and variable turbine geometry features. That is why further research will focus on the influence of the different positions of the second-stage nozzle vanes to control the backpressure in the exhaust system. The final milestone is the manufacturing of the prototype device with a two-stage expansion system and adjustable nozzle vanes at the second-stage. Such a system will be used in a 2-stroke, 6-cylinder OP engine which is being developed at the Warsaw University of Technology.