Abstract

The transient response of a radial turbine stage with a variable geometry system is evaluated. Mainly, the consequences of the variations of the aerodynamic performance of the stage on the response time are checked. A simple quasi-steady model is derived in order to formalize the expected dependences. Then an experimental campaign is conducted: a brutal step in the feeding conditions of the stage is imposed, and the response time in terms of rotational speed is measured. This has been reproduced on different declinations of the same stage, through the variation of the stator geometry, and correlated to the steady-state performance of the initial and final operating points of the transient phase. The matching between theoretical expectation and results is surprisingly good for some configurations, less for others. The most important factor identified is the mass-flow level during the transient phase. It increases the reactivity, even far above the theoretical expectation for some configurations. For those cases, it is demonstrated that the quasi-steady approach may not be sufficient to explain how the transient response is set.

1. Introduction

Many applications make use of the interesting properties of radial inflow turbines, namely, the compactness and the ability of reaching high pressure ratios for acceptable efficiencies. The description of the radial turbines operating in steady flows has been widely presented in the literature [1, 2], together with an extensive description of the flow pattern, even for off-design conditions (e.g., in [3]). Some design processes have also been presented to ensure the good performance of the energy-recovering process of the stage (see [4, 5]). The favourable pressure gradient due to the expansion of the flow through the turbine stage stabilizes the boundary layer behaviour and makes separations unlikely if the guidance of the flow is appropriate. As a consequence, a one-dimensional description of the flow is adequate to reach an accurate description of the radial turbine functioning [4]. This is transposable to operational and geometric off-design operation [6]. But the problem gets more complex when the inlet conditions are time dependent. Most of the contributions analysing radial turbines operating in unsteady flows are found in the literature dedicated to possible turbocharging applications. The turbine stage experiences some severe pulsating inlet conditions due to its position in the exhaust manifold. The research activity in that field is currently very active. A review of the state-of-the-art and of the main contributors of the field can be found in [7, 8]. The fact that the rotor part of the stage operates in a quasi-steady manner is one of the most significant conclusions of the different studies. Even if this assessment needs to be generalized (far from the design point, the rotor is likely to be more sensitive to some unsteady processes), it points out that the performances of the static part of the stage (volute and stator) are strongly affected by the time-variation of the inlet conditions.

Another kind of time-dependent operation is defined when dealing with transient phases of the turbine stage. There again, most of the research is motivated by automotive propulsion, since the turbocharger often experiences such a kind of request on a typical run. The literature depicts how the response time of the turbocharger impacts the transient response of the engine. This is presented through the modelling of the engine dynamic response when coupled to a turbocharger, compared with experimentations (in [9, 10]). The turbocharger is then modelled as a system having a first-order transfer function, which is a good representation according to [11]. The general assessment is that the mass-moment of inertia of the rotating assembly is the most dominant factor influencing the response time of the turbocharger (e.g., in [12]). A direct consequence of this variation of the mass-moment of inertia of the turbocharger shaft on the engine overall response is presented in [10]. The issue of reducing the mass-moment of inertia is that the turbine stage must sustain high temperatures and thus requires high-density materials of construction. Some solutions, such as scalloping, are recommended to decrease as much as possible this moment of inertia [12], even if it deteriorates the aerodynamic efficiency. Anyway, the presumed importance of the moment of inertia of the turbocharger shaft is still in discussion. A great part of this discussion is focused on the overall reactivity of the engine, for which the mechanical, chemical, and thermodynamic inertia of the engine itself should be considered (according to [9, 13, 14]). A marginal part of the discussion deals with the isolated turbocharger to understand how its own transient response is affected by other factors. One solution to increase the reactivity is the use of variable geometry devices either on the compressor or on the turbine stage. This is demonstrated in different studies [13, 15, 16] and highlights once again the importance of the static part of the stage, here on the dynamic response of the rotor. The control strategy of the variable geometry device is then crucial as demonstrated by [13]. In [14] is mentioned the fact that the efficiency of the stage is of the same order of importance as the mass-moment of inertia in the establishment of the response time. It is indeed surprising that such an extensive knowledge of flow pattern in radial turbines, as described above, is that poorly correlated to the dynamic response in transient operation. No clear correlation between the dynamic performance and the efficiency of the turbine is presented. Since the definition of efficiency in unsteady flows is a complex matter [17], the only indicator available is the efficiency measured in steady flows. The comparison of time scales of the stage (rotor rotation, propagation of the pressure, or mass wave in the stator) to that of a typical transient request tends to validate this “steady” indicator as representative of the performance. But this is questionable since a pressure wave will not be instantaneously distributed all around the rotor and will definitely affect the boundary layers in the stage. The present work is a first contribution to understand how the efficiency of a radial turbine stage interferes with its dynamic response for different geometric configurations. The first part of the paper focuses on the theoretical background in order to identify the most relevant factors of the problem. The second part of the paper presents the experimental approach and the facility. The last part of the paper presents the results which are discussed in the conclusion.

2. Theoretical Background

The time evolution of the rotational speed during transient operation is based on a rather simple mechanical equilibrium. If the turbine stage is connected to another device through a solid crank arm (compressor, electric generator, etc.) for which the torque is quoted , it can be written as follows: Since if there is no gearbox plugged between the two systems, in steady operation both the torque and the rotational speed of the turbine and of the braking device are equal (the contribution of the bearings is taken into account as part of the “braking” device). For the turbine, in steady-state operation, the power extracted from the flow is given by the first law of thermodynamics: This statement forbids the direct integration of (1), because the dependence of the quantities (mass flow, pressure ratio, and efficiency) to the rotational speed is not explicit. In steady-state operation, a complete map is required to identify this dependence; it does not have a simple analytical expression. And in transient operation, the validity of this map is not demonstrated. Equally, the braking work () depends on the rotational speed and also on the characteristics of the device to which the turbine stage is connected.

The exact modelling of the transient behaviour is not the purpose of this work. But a theoretical reference is needed to meet the two main objectives listed here:(1)the identification of the most influential factors controlling the setting of the response time (in addition to the obvious contribution of the mass-moment of inertia); (2)the verification of the validity of a quasi-steady approach in transient processes. In that purpose, the most simple quasi-steady theoretical expression of the response time for transient operation is built. The transient phase is defined as follows: from an initial steady-state operating point (denoted with “i” indices), a sudden rise (characteristic time of the step < 10 ms) of pressure and mass-flow at inlet of the stage leads to a second steady-state operating point (denoted by “f” indices). One has to drastically simplify the problem to integrate (1). Those simplifications are as follows:(i)an intermediate point is defined (denoted by “i′” indices). The transient phase is then decomposed in two subphases [i] and [f], as illustrated in Figure 1;(ii)the first subphase ([i]) considers the instantaneous variation of the inlet pressure and the mass flow to their final values, at a fixed rotational speed. This is supported by the characteristic velocity of propagation of the pressure perturbation inside the turbine stage which is near the speed of sound (as presented in [7] and discussed in [8]); the time scale is thus shorter than the one of the shaft velocity variation, which keeps its initial value. Then , , and . The torque is discussed later;(iii)the second subphase ([i′–f]) consists in the variation of the rotational-speed due to the difference of torque between turbine and braking device;(iv)a linear evolution of the torque-to-rotational-speed characteristic for both the turbine and the braking device is assumed since it is difficult to evaluate it a priori (see Figure 1).In that quasi-steady situation, a simple expression for both the turbine and braking torques can be obtained. On the braking device side, the functioning characteristic is linear and must be compatible with the final operating point of the turbine (to reach the power equilibrium in steady state operation). It means that The same treatment is applied on the turbine side, namely, linear variation and compatibility with the final operating point: Equation (1) then gives this equation reduces to a first-order differential equation as follows: where is a constant. For such an equation it is possible to identify a time parameter which is representative of the response time of the shaft. Its expression is An estimation of the torque generated by the turbine at point i′ is difficult to give. Such a transient operating point is hardly accessible by measurement. The corresponding steady-state operating point, if it exists, is generally out of the range of conventional test benches, since a large overlapping zone between the speed-lines of the map is required to capture it. This is hardly possible when the braking device is a compressor, limited by blockage and surge/stall regions. A last estimation is then proposed for the torque at this intermediate step of the transient, in order to close the model as simply as possible.(v) The variations of mass-flow, efficiency, and pressure ratio are supposed to be small between the points and f. Then , which leads to the following approximation: The time parameter is written as follows: where is a parameter only depending on the initial and final rotational speeds. It is written as follows: Finally, according to (2) and (3), and to the definition of isentropic speed, one gets The response time to a brutal modification of the operating point from the steady conditions “i” to the steady conditions “f” obviously depends on moment of inertia of the shaft. But it also depends on the initial rotational speed (through ) and on the final operating point performance. It means that a given jump of rotational speed will not have the same duration, depending on the region of the map considered (illustrated in Figure 2). According to the simplified model proposed in (12), decreasing the moment of inertia will improve the reactivity of the stage at every position on the map, but a general improvement in the stage efficiency will have the same effect. Otherwise, if a general improvement cannot be achieved, this simple relation points out the fact that the low mass-flow configuration should be favoured in efficiency, since the high mass-flow regions are naturally reactive.

This result has been derived with the restrictive hypotheses above mentioned. The validity of these hypotheses is also dependent on the position of the operating in the map, and on the geometric configuration of the stator. For example, considering the mass flow constant during the subphase [f] in a transient which is near the blockage region seems more valid than in other regions of the map, or the modification of the stator geometry may shift the reduced tip speed from its optimum value of 0.7. The real dependence of the transient response on the factors identified still needs to be checked.

3. Experimental Approach

The approach is based on a single geometric configuration of the rotating assembly: radial inflow turbine and centrifugal compressor (see Figure 3(a)). The stator of the turbine stage is modified to create different flow configurations and check the dependence of the factors identified in (12). Three modifications are imposed and presented in Figure 3(b): (i) opening angle of the stator blades, (ii) height of the stator blades, and (iii) presence of the spacers and clearances in the stator. The first two modifications aim at the alteration of the mass-flow level inside the stage without changing the final pressure ratio. The last modification interferes with the efficiency of the stage, for the presence of the spacers creates some clearances at the tip and at the hub of the stator blades and generates some losses. It is an artificial method to change the efficiency of a given operating point. When dealing with variable geometry stators, the importance of the reduced section () is demonstrated in [6]. This parameter is a nondimensional expression of the throat section of the stator; its expression is given in (13). The setting of this reduced section defines the operating point of the stage. The stator opening angle and height combinations are chosen so that some different geometric configurations generate similar values of reduced section. For most of the configurations, some tests with and without spacer/clearances are conducted. The different geometric configurations are labelled according to the opening position of the stator (1: closed; 2: intermediate closed/reference; 3: reference5: opened) and subscripted with the height of the stator (s: small; r: reference; h: high). The configurations tested and the associated values of the parameters are presented in Table 1:

3.1. Steady-State Measurements

The “steady-state” identification of the initial and final operating points is obtained on our PTM (Petite TurboMachine) test rig dedicated to turbocharger applications. This steady flow test rig allows either global or local instrumentation of small stages, for a wide range of inlet temperatures (from 20°C to 620°C), up to 250 g/s. The air is supplied by a pressure source (6 bars) stabilized by two regulation stages. The compressor and turbine flows are independent. They are thermally isolated from the ambient air and from each other. A specific circuit supplies hot and pressurized oil to the bearings, with a proper regulation. Inlet/outlet pressure and temperature are measured in order to get the performance of both stages. Mass-flow measurement is operated through a Coriolis flowmeter. Rotational speed is measured using a proximity sensor on the compressor stage. For more details about the test facility, see [18]. The repeatability quality of the experimental device is the good indicator for the experimental uncertainties since the experimental study deals with comparative results. The confidence ranges are given in Table 2 for the nondimensional quantities presented in the results.

For machining convenience, the prototype stators were built with aluminum alloy (AU4G) which does not sustain high temperatures. The air-flow temperature supplied to the turbine stage is thus left at ambient value. The different pressure ratio lines are then described ( from 1.35 to 1.6) by modifying the loading of the compressor, which is controlled by a discharge valve.

3.2. Unsteady Measurements

An improvement of the PTM test rig is developed, to meet unsteady measurement requirements. A rapid solenoid valve generates a sharp pressure variation at the inlet of the turbine stage. A flow derivation supplies air for the stage when the valve is closed, so that the rotor never stops. This flow supply defines the initial operating point. Then the valve is opened. A rapid modification of the inlet conditions (the typical time of variation of the inlet pressure is 10 ms) forces the turbine stage to reach the final operating point. A typical transient phase imposed to the stage by the valve opening, in terms of initial and final pressure ratios, is from to . Some unsteady measurement devices were implemented: two fast-response pressure probes at the inlet of the stage and at the outlet of the compressor stage; a fast-response proximity sensor near the compressor impeller to detect the blade passage.

The pressure sensors are conditioned by in-house amplifiers designed to set the cut-off frequency 10 times greater than the maximum frequency expected in the pressure excitation signal. The sampling frequency of the data acquisition system is set at 200 KHz to ensure an accurate resolution in the observation of the blade passage by the proximity sensor. The raw signal of the blade passage recorded is postprocessed to extract the temporal evolution of the rotational speed.

The initial and final operating points are identified by the use of the steady-state metrology. The extraction of the response time of the shaft is performed through the postprocess of the unsteady measurements. For safety reasons, some regulation of the pressure upstream of the solenoid valve has been implemented. This regulation process is less rapid than the solenoid valve and suffers some lag. This makes the pressure step at the inlet of the stage unperfect (see Figure 4). Thus the response time is not accessible through direct measurements of the raw signal. The shaft should behave like a first-order system according to (7) since the initial and final operating points are fixed. It has thus been decided to use an identification function to express the impulse response of our system. The entry of the system is the pressure rise at inlet, and the response is the rotational speed evolution. The mean square method is applied to match the coefficients to the experimental evolution of the two quantities. The production of the final model is then compared with the data. The matching between the prediction and the experimental data is in the range of 98% for all of the configurations. The characteristic response time of the first-order model is taken as the estimation of the actual response time of the stage. Ten repetitions of each experiment give access to the standard deviation of the response time and to the associated confidence range at 95%. This range is expressed in the following results as error bars around the response time data (typical value: ±2.5 ms).

3.3. Numerical Simulations

The intermediate torque () is only estimated in the model and never measured. As stated above, the validity of the estimation depends on the location of the transient phase in the map and cannot be experimentally confirmed. Some numerical simulations were punctually conducted in order to reach estimations of the torque for some operating points out of reach of the test bench. Those simulations were performed with the Euranus solver, of the Fine/Turbo software suite of Numeca Int. The stator-rotor configuration is meshed with Autogrid 5 (1.5 million cells, ). The steady simulations (mixing plane condition at the rotor/stator interface) use the Spalart-Allmaras turbulence model. The inlet total conditions, outlet static pressure, and rotational speed are imposed as boundary conditions.

4. Results

In this section, the analysis of the results is divided in two parts. First, one focusses on the factors highlighted by the theoretical development for their supposed importance in the establishment of the response time (mainly the final mass-flow and the final efficiency). Second, the discrepancy of the theoretical model is analysed to feed the discussion about the relevance of the quasi-static approach for transient operation.

4.1. Identification of the Influential Factors

Figure 5 presents the general influence of both the mass-flow and the efficiencies of the final operating point on the transient response of the stage. The mass-flow and efficiency plotted values are relative to that of the final point for the reference configuration (, with spacers). The steady-state part of the results (Figure 5(a)) shows the performance of the stage for three opening configurations (closed: , nominal: , and fully opened: ). The value of the final pressure ratio is the same for the three configurations (). The results confirm that it is possible to downgrade the efficiency by the use of the spacers and the clearances, without drastic modifications of the operating point. The penalty induced by the spacers and the clearances is far more important in the closed configuration than in the fully opened one. This has already been explained long ago by [19] since the relative importance of the clearances increases as the throat section of the nozzle decreases. In response to the alteration of both the mass-flow and the efficiency levels, the reactivity of the stage is modified as presented in Figure 5(b). The general trend is in agreement with the theoretical expectation: when the mass-flow level in the stage increases at the end of the transient phase, the response time decreases. It also decreases when the efficiency is increased (for the same mass-flow level across the stage), even if the tendency is less clear. For the closed configurations, the theoretical prediction based on the steady-state measurements of the initial and final operating points is surprisingly accurate. This accuracy is lost for the two other configurations, for which the prediction strongly overestimates the response time: the shaft is almost twice more reactive than expected by the model. The strong dependence of the response time on the final operating point definition is examined in Figure 6. The steady-state results show the position of the initial point and that of the three final operating points in the turbine map, for the stage configuration . The gap in rotational speed between initial and final operating points is successively increased (from = 3000  for transient 1 to = 6000  for transient 3), together with the final mass-flow rate and final pressure ratio, as shown in Figure 6(a). Figure 6(b) presents the response time of the three imposed transient phases. It decreases despite the apparent increase of the distance in the map between the initial and final points. The model also predicts this trend. Equation (12) shows how the “geographic” distance between the initial and final operating points in the map is irrelevant in the prediction of the response time; the proximity obviously does not ensure a short transient phase. The good reactivity is better ensured by a large mass-flow through the stage at the end of the transient and to a lesser extent by a good efficiency.

The results examined in Figure 6 also reinforce the conclusion that in certain conditions the experimental value of the response time is strongly overestimated by the quasi-static approach. The accuracy of the theoretical prediction is thus investigated.

4.2. Validity of the Quasi-Static Approach

A general evaluation of the theoretical model is proposed in Figure 7. All the configurations are classified in terms of reduced section, as presented in Figure 7(a), where the corrected mass-flow of the final operating point is plotted. This diagram shows the importance of the value of the reduced section () for variable geometry stages: for a given pressure ratio, the mass-flow values of the operating point collapse on a single trend. This occurs despite the fact that the different stators differ geometrically (see [6] for more details). The reduced section actually defines the permeability of the stage, as long as the minimum section of the stage is in the stator. This classification in terms of reduced section is important, because the validity of the theoretical prediction of is partially related to it. A reversal formulation of (12) is

This specific formulation is convenient because there is a separation between the response time and the aerodynamic performance on the left-hand side (which are results), and the request in terms of final rotational speed with the characteristics of the shaft on the right-hand side (which are imposed by the transient operation). Both terms of the equation have the dimension of a mass. On the left-hand side, this “transient mass” is roughly the mass of fluid crossing the stage during the transient phase and is supposed to evolve linearly with the square of the final reduced tip speed of the rotor. The slope of the line () is a characteristic mass of the shaft, which is common to all configurations. This supposed linearity is checked in Figure 7(b) where the “transient mass” is plotted as a function of the square of the reduced tip speed. What was observed in Figure 5(b) is here generalized to every configuration tested. The experimental results and the theoretical predictions are in a good agreement for some configurations, those having a small value of the reduced section. A strong discrepancy appears for the others, when the reduced section of the geometry is increased. The actual response time is then much shorter than the prediction of the quasi-static approach, almost halved. For the theoretical model it means that(1)either the quasi-static approach is correct, and all the other approximations (such as the decomposition of the transient phase in two subphases, the linearity of the torques, and the estimation of the intermediate torque ()) lead to the deviation of the model; (2)or those approximations are acceptable but with the quasi-static approach some phenomena in the transient phase are out of reach.

Whatever its origin is, the deviation is inactive for the configurations having a small value of the reduced section and it is very intense for the others. This duality is surprising since the different stages tested are not extremely different from one another in terms of geometry, and the hypotheses of the model are always the same.

Some additional tests have thus been carried out, to incriminate one of the two previous propositions. For the first one, the possible decomposition in two subphases is supported by the literature, and by recent unsteady measurements in the intermediate locations of the stage; it seems reasonable enough. The linearity of the torque during the phase [f] has been checked for the braking torque, by reproducing the expected transient phase through a succession of steady points. In that quasi-static approach, the linearity has been observed. For the turbine, the same approach has not been possible on the complete segment [i′–f] because of the limitations imposed by the blockage and surge regions of the compressor used as “braking device,” but the portion observed was also linear.

Only two possible explanations remain: a strong underestimation of for some configurations and the occurrence of an unsteady phenomenon activating the response time. No measurement of can be reached experimentally (the steady representation of point corresponds to a high mass-flow () and pressure-ratio () and a small rotational speed (); it is far over the possible absorption of energy of our compressor at this rotational speed, for steady conditions). But though (8) it is possible to back-deduce a value of from the experimental results of the response time. This back-deduction quoted is thus the value that the instantaneous torque should take at point to match the measured response time. It is compared in Figure 8(a) to the initial estimation of given by (9) and implemented in the quasi-static model. For the configuration presenting a small reduced-section value (white symbols in Figure 8(a)), as another expression of the good matching between the prediction and the results, we find that . For intermediate and high values of the reduced section (grey and black symbols), the torque back deduced from the response time measurement is far more important, up to 30% more than the approximation. Indeed, the approximation proposed in (9) is quite simplistic and does not have an equal validity for the different stator configurations or the different regions of the map. This validity is checked through the simulation of the operating point with CFD calculations, since it is out of the possible range of the test bench. Some simulations of the operating points and were performed for the stators and (see Section 3 for the details of methodology). The simulation of the point aims at validating the torque predicted by CFD, compared with the steady results. Then the same methodology and postprocess are applied for the prediction of the torque for . This prediction is hopefully more accurate than (9) for the point and compared with the results of the quasi-steady approach ( simulated compared with back-deduced). The results are presented in Figure 8(b) and show that the approximation of (9) is acceptable. More important, even with the proper estimation of the quasi-static prediction still fails to match the measurement, even if it gets closer to it. For the nominal configuration a gap is still present. For the full-opened configuration, the torque at point is still more than 20% lower than the value compatible with the measured response time.

Summarizing, the quasi-steady inaccuracy increases as the characteristic time of the transient decreases. Since the short response-time cases are also the ones for which the mass-flow is important, this inaccuracy is observed for the configurations presenting the less restrictive stator geometries (full-opened geometries). The fact that the response time is significantly shorter than the one predicted is noteworthy. Thus the possibility of a phenomenon out of the scope of the quasi-steady approach rises but still needs to be verified. The creation of an unsteady “overtorque” due to the passage of the pressure front across the stage during the brutal transient could be considered. Since the pressure front is probably damaged by the very restrictive closed nozzles, it is not active for the small reduced-section configurations. When the front keeps some coherence through the stator, the transient phenomenon appears, and the torque is reinforced to an instantaneous value . This has not been demonstrated and is still speculative. It has not been observed but it is not surprising since most of the unsteady results found in the literature concern periodic pulsating flows, for nominal geometries of the stator (for those configurations the transient phenomenon is not strongly active).

5. Conclusion

Some measurements of the transient response of a variable geometry turbine have been conducted and compared to some theoretical predictions. The matching is good for some opening configurations of the stator, less for others. More specifically when the reduced section of the stator is important, some very strong reactivity has been observed. The different conclusions can be summarized as follows.(1)It is possible to change the transient response of a shaft without changing its polar moment of inertia. A substantial gain of efficiency will produce a gain of reactivity. (2)The main factor is the mass flow rate. To increase the reactivity of a given shaft through a transient operation, the most efficient way is to increase the mass flow level. The geometry of the stator should be adapted to the increase of mass-flow rate. (3)For configurations presenting an important value of the reduced section, the actual response time is significantly smaller than the quasi-steady prediction. (4)An unsteady creation of torque is suspected, if not demonstrated. The origin of this “overtorque” is conjecturally related to the passage of the pressure front. (5)If the geometry of the stator is too restrictive, the “overtorque” is ineffective and the quasi-steady approach is accurate enough to predict correctly the response time.

This process needs to be verified in detail. But should it be demonstrated, this could explain some of the unsteady measurements of efficiency presented in the literature which are above 1, even if it is also agreed that the inappropriate definition of the efficiency in an unsteady context is also part of such amazing results.

Nomenclature

Quantities
: Specific heat at constant pressure ()
: Isentropic speed ()
: Rotor inlet diameter ()
H: Stator-blade opening height ()
J: Mass-moment of inertia ()
: Stator throat section
: Stator reduced section
T: Temperature ()
t: Time ()
: Rotor inlet tip speed
: Power
: Mass-flow rate
: Stator-blade opening angle
: Rotational speed
: Efficiency
: Specific heat ratio
: Response time
: Turbine pressure ratio
: Torque
: Transient parameter.
Subscripts
B: Relative to the braking device
f: Final point of the transient
i: Initial point of the transient
: Intermediate point of the transient
rt: Reduced quantity
t: Total state
T: Relative to the turbine.