1 Introduction

The flow across the gap between the rotor blade and the stationary turbine case, known as blade tip-leakage flows, constitute a major source of inefficiency in the turbine, accounting for up to 1/3 of the overall losses in a turbine stage [1]. The working fluid passes through this gap, thus the potential energy existing in this part of the fluid it is not extracted. Because of the rotational speed and the relative movement between the rotor blade and the casing, in the tip clearance area, a vortex forms and it is directed in the blade passage by the pressure difference between the pressure and the suction side of the rotor blade. Here, the vortex interacts with passage vortices to form a complex flow pattern. In general, the efficiency decrease linear with the increase of the clearance gap size, although the line does not pass through the origin, and the relation is presumably nonlinear at very small gaps (less than 1% chord) [2].

In the designing phase a small clearance is desirable to minimize the flow across this gap, thus minimize the losses, but a minimum gap must be provided to avoid the potential rub between the two elements. Because of the centrifugal forces on the rotor and the different thermal inertia of these elements, the tip clearance gap is not constant across the engine operating envelope, moreover, the variation is not linear. Hany Moustapha et al. [3] described the variation of tip clearance through a typical engine cycle and stated that during acceleration, due to a rapidly growing in rotational speed and centrifugal forces, rub between rotor blade tip and turbine casing is possible. During the deceleration another rub can occur if there is a full power demand at the end of the deceleration phase. Figure 1 describes the variation of the clearance tip during the acceleration and deceleration process.

Fig. 1
figure 1

Variation of tip clearance through a typical engine cycle (constructed using [3])

In the context of developing a 400 N microjet engine, an axial turbine is designed to drive a centrifugal compressor. A particular problem occurs when the fuel mass flow is calculated for the dynamic regimes, especially for the lower engine speeds. Since the turbine theoretical map is calculated with the same clearance gap for all the operational power settings, during the dynamic regimes the clearance gap changes as a result of different elements expansion (or contraction), thus at this operational conditions the turbine map information are not precise.

This fact leads, in case of acceleration, to a high turbine inlet temperature, when the fuel mass flow is too high as a result of underestimating the turbine power, or a slow acceleration when the mass flow is too small as a result of overestimating the turbine power.

The present study numerically investigated the effect of tip clearance and the tip clearance variations to turbine and engine performance. For this study, it was considered that the turbine follows the typical tip clearance variation trend presented in Fig. 1. Thus, no heat transfer analysis was conducted, the blade geometry was considered constant and the tip clearance variation was determined by the variation of the turbine case diameter.

2 Numerical Setup

The small size turbine consists of a single-stage axial flow with both, the vanes and the blades, designed with constant profile across the radius. The turbine geometric information is presented in Table 1 and the 3D models for the vane and the rotor are presented in Fig. 2.

Table 1 Turbine information
Fig. 2
figure 2

3D models of turbine vane and rotor

The grid is structured and has approximately 1,800,000 nodes in total, with the smallest cells close to the walls. To reduce the computational time, the numerical study is conducted on a single channel using periodic boundaries condition. At the interface between the stator part and the rotor part, the frozen rotor model is used. Figure 3 presents the grid for the vane and the blade.

Fig. 3
figure 3

Numerical grid for stator and rotor blade

At the inlet, it was imposed the absolute pressure, 472,600 Pa, and the absolute temperature, 1173 K. The reference pressure is 0 atmosphere and the rotational speed is 80,000 rpm. Also at the outlet, it was imposed the mass flow, namely 0.685 kg/s per machine. In Fig. 4, the channel considered for the numerical study is presented with the respective boundary conditions. The turbulence model used for this numerical study is the k-epsilon model coupled with a wall function.

Fig. 4
figure 4

The computational domain

To verify that a proper mesh was created and to ensure an accurate prediction of parameters, especially in the walls vicinity, the cell distribution near the walls is verifying in terms of y + parameter. Figure 5 presents the distribution of y + across the rotor blade. The y + values are acceptable considering the turbulence model used and the fluid evolution through the stage. As the main objective of this study is to determine the influence of the tip clearance variation on the thermodynamic performances of the turbine and not the precise results for every turbine map, no further analysis of the numerical results precision was conducted.

Fig. 5
figure 5

y + distribution across rotor blade

3 Results and Discussion

For this study, five models were used, corresponding to five values of the clearance gap, from an ideal configuration, without a clearance gap to a large gap, corresponding to approximately 1.1% of blade chord. In each case, the blade geometry is identical, with the exception of the rotor shroud which increases with the value of the respective clearance gap. For each of the five computed cases, major changes in flow next to the tip of the blade and changes in the mass flow that passes through the gap were observed. Figure 6 presents the flow through a section near the rotor shroud and can be deduced, as expected, that the flow through this section and the influence on the flow in the vicinity is proportional with the clearance gap.

Fig. 6
figure 6

The flow near the rotor blade tip for a 0.3 mm, b 0.25 mm, c 0.2 mm, d 0.15 mm, e 0 mm clearance gap

The turbine performances can be compared for different tip clearances by comparing the turbine characteristics. To obtain the complete maps, 20 cases were calculated for each clearance gap. By using specific parameters, calculated from the numerical results, the maps were drawn. The specific parameters are defined in Eq. (1) [4]:

$$\left\{ \begin{gathered} \Delta h_{\text{Tred}}^{*} = f\left( {\dot{M}_{g} \frac{{\sqrt {T_{3}^{*} } }}{{p_{3}^{*} }},\frac{n}{{\sqrt {T_{3}^{*} } }}} \right) \hfill \\ \eta_{T}^{*} = f\left( {\dot{M}_{g} \frac{{\sqrt {T_{3}^{*} } }}{{p_{3}^{*} }},\frac{n}{{\sqrt {T_{3}^{*} } }}} \right) \hfill \\ \end{gathered} \right.,$$
(1)

where: \(\Delta h_{\text{Tred}}^{*}\)—specific enthalpy drop; \(\eta_{T}^{*}\)—turbine efficiency; \(\dot{M}_{g} \frac{{\sqrt {T_{3}^{*} } }}{{p_{3}^{*} }}\)—reduced mass flow; \(\frac{n}{{\sqrt {T_{3}^{*} } }}\)—reduced rotational speed.

Using ANSYS CFX, four rotational speeds were simulated (80, 60, 40 and 20,000 rpm), for each of the five mass flows (representing 100%, 95%, 90%, 80%, respective 70% of nominal mass flow). In Fig. 7, it is shown the turbine map for 0.2 mm clearance gap.

Fig. 7
figure 7

Turbine map for 0.2 mm clearance gap

The influence of the clearance gap in turbine performances can be seen in Fig. 8 for the higher regimes (100% mass flow) and in Fig. 9 (70% mass flow). By increasing the gap, the performances drop.

Fig. 8
figure 8

Turbine performance comparison for five values of the clearance gap (100% mass flow)

Fig. 9
figure 9

Turbine performance comparison for five values of the clearance gap (70% mass flow)

By comparing the effect of the clearance gap on the turbine performances, for higher and lower regimes, can be seen, from Fig. 10, that at higher regimes the clearance gap has a higher impact on the performance. A more than 10% drop in power was obtained for the maximum power regime and approximately 5.3% for a lower regime (70% mass flow). In both cases, the curves presents a linear decrease in performance.

Fig. 10
figure 10

Power variation with the clearance gap for higher and lower regimes

The maximum gap considered in this paper is 0.3 mm which represents approximatively 1.71% of blade height. In Table 2 there are present the percent of blade height and chord.

Table 2 Gap percent of blade height and chord

The present result agrees with the study of Haas and Kofskey [5] and Bindon [6] which determined a linear variation of performance with the clearance gap although they predicted a nonlinear relation for gaps less than 1% of chord.

To determine the influence of the clearance gap on the whole engine the compressor map is needed. Both the compressor and the turbine map are completed by using commercials programs like Smooth C respectively Smooth T and GasTurb. The performance map of the single stage, 5.1 nominal pressure ratio, the radial compressor is presented in Fig. 11.

Fig. 11
figure 11

Compressor map

The turbine map for an ideal configuration without clearance gap is presented in Fig. 12.

Fig. 12
figure 12

Turbine map

Using Transient module in GasTurb 13, the engine acceleration process from idle to maximum power was simulated by the use of engine cycle and component maps. For each clearance gap, a turbine performance characteristic was created and then the acceleration process was simulated. Figure 13 presents the working line on the turbine performance map (determined for a constant 0.2 mm clearance gap) for the process considered.

Fig. 13
figure 13

Working line on the turbine performance map

In Fig. 14, it is presented the turbine inlet temperature evolution for each case during the dynamic process.

Fig. 14
figure 14

Turbine inlet temperature variation during dynamic operation

Because of the drop in turbine performances, with the increase of the clearance gap, to supply the necessary power to the compressor, a higher turbine inlet temperature is needed. Thus, we cannot correctly approximate the working line with only one value of the clearance gap.

By approximating the dynamic of the clearance gap and by using the working lines calculated before we can draw a more precise working line. The variation of clearance gap was approximated by scaling a typical evolution for axial turbines (presented in Fig. 1). On the rotational speed interval, from idle to nominal speed (45,000–80,000 rpm), the tip clearance variation was determined to achieve the designed gap (0.2 mm) at a nominal point. Thus, the turbine maps have been used as presented in Table 3.

Table 3 Turbine maps variation during acceleration

Figure 15 presents the turbine inlet temperature variation during acceleration corrected with the dynamic variation of the clearance gap.

Fig. 15
figure 15

Turbine inlet temperature variation considering multiple clearance gaps

A similar curve describes the fuel flow variation in Fig. 16. Knowing the precise power needed by the compressor, the correct temperature can be determined, thus the fuel flow needed can be predicted. For this simulation, the fuel flow used for the acceleration process, considering also the variation of the clearance gap using multiple turbine maps, was approximately 1.5% lower (globally) than considering only the nominal clearance gap. This value of lower fuel flow needed was calculated across the entire acceleration process, the clearance gap being lower or greater than the nominal one at respective regimes, thus the savings in fuel is limited. As for the maximum difference in the fuel flow, in this simulated process, an approximately 8% higher value was used which would determine a raise in temperature by 20 degrees, compared with the optimal acceleration process.

Fig. 16
figure 16

Fuel flow consumed considering constant and variable clearance gaps

The working line from Fig. 17 was calculated across multiple turbine maps as it is described in Table 2.

Fig. 17
figure 17

Comparison between working lines for constant and variable clearance gap during simulating acceleration process

4 Conclusions

This paper conducted a study of the effect of the clearance gap on the turbine and engine performances during dynamic operations. Using a model of the dynamic response of the clearance gap, a few values were determined and then for each value a complete turbine map was created using CFD calculation and specialized commercial programs like Smooth C respectively Smooth T. By comparing the results, it was determined that the clearance gap effects on the turbine performances are more visible at higher regimes. The turbine inlet temperatures, required for acceleration in this simulation, considering multiple turbine maps, was determined as well as the fuel flow. It was determined that a 1.5% lower fuel flow was needed by considering the variation of the clearance gap during acceleration. The working line was also described in this paper across multiple turbine maps.

Future work should concentrate on determining an empirical coefficient to correct the fuel flow calculated to consider the effect of the clearance gap variation on the turbine performance. In addition, an experimental validation on these results is required.