The aim of the present study is to evaluate the efficiency of a boundary layer suction technique in case of a centrifugal compressor stage in order to extend its stable operating range. First, an analysis of the flow pattern within the radial vaned diffuser is presented. It highlights the stall of the diffuser vanes when reaching a low massflow. A boundary layer separation in the hub-suction side corner grows when decreasing the massflow from the nominal operating point to the surge and finally leads to a massive stall. An aspiration strategy is investigated in order to control the stall. The suction slot is put in the vicinity of the saddle that originates the main separating skin-friction line, identified thanks to the analysis of the skin-friction pattern. Several aspiration massflow rates are tested, and two different modelings of the aspiration are evaluated. Finally, an efficient control is reached with a removal of only 0,1% of the global massflow and leads—from a steady-state calculations point of view—to an increase by 40% of the compressor operating range extent.
The stable operation of compressors is limited towards low massflow by the occurrence of flow instabilities. To prevent the release of these instabilities, that can be destructive, compressors never operate closer to the stability limit than a chosen massflow. That is called the “surge-margin.” Essential for the stability of the compressor working, it is however in this forbidden part of the operating map that the compressor delivers its maximum pressure ratio, and the surge margin is thus detrimental to the performance of the stage.
The flow phenomena that cause compressor surge are not yet well known. They need to be understood in order to drive design changes that can increase the operating range of the compressor. For the centrifugal compressors, several studies have identified the diffuser as a key zone for the inception of instabilities. In particular, large boundary layer separations may occur in this element as the operating point moves towards surge, leading to rotating stall and/or surge.
A number of investigations have then been done, trying to improve the surge margin using control techniques. Raw [
In the same series of experiments, Skoch also evaluated the effects of hub-side injection into the diffuser throat, since PIV measurements had revealed a hub-corner stall on the pressure side of the vanes of the diffuser, that could suggest the use of air injection into the diffuser passage to energize the flow and control the separation. But these attempts did not demonstrate any significant surge-margin improvement [
The present work also focuses on the development of a control technique in order to extend the stable operating range of a centrifugal compressor. But instead of using injection, and decreasing the flow angle, it aims at controlling the separations that may occur, by using a boundary layer suction technique. That control action has indeed already been used in axial compressor stages or cascades and has provided attractive results [
An analysis of the flow is first presented, based on the solutions of both steady-state and unsteady simulations. Skin-friction patterns are drawn and reveal the topology of the flow.
Considering the localizations of the singular elements and separation skin-friction lines, an aspiration strategy is discussed and then implemented in the steady-state simulations. Several massflow rate removals are tested, and numerical performances of the aspirated compressor stage are presented.
The test case is a centrifugal compressor stage manufactured by TURBOMECA, Groupe Safran, used for propelling helicopters, and that is representative of the present state of Turbomeca’s technical know-how. It is composed of a backswept splitted unshrouded impeller, a radial vaned diffuser, and an axial diffuser. The relative Mach number at impeller leading edge at shroud is about 1, whereas the flow in the diffuser is transonic. The axial diffuser is not taken into account in this study.
Experimental data used for computing the reference performance curves of the test case were provided by Turbomeca and were obtained during the characterization of the compressor stage. In particular, the static pressure was probed at the entry of the impeller, in the middle of the vaneless diffuser, and at the diffuser outlet. As an example, at the diffuser outlet, the experimental value of static pressure is the result of an average over 31 holes equally distributed along the circumference and all connected to the same annular chamber where the pressure is measured.
Computations were performed with the elsA software developed at Onera. This code solves the compressible RANS equations by a cell-centered finite volume method. The turbulence model chosen for this study is the two-equation model of Smith (k-l), that was shown by Rochuon [
For steady-state simulations, the mixing-plane approach is used for modeling the impeller-diffuser interaction, and the computation domain is reduced to one single passage per row thanks to a hypothesis of uniformity of flow within all passages.
For the unsteady simulations, a phased-lagged hypothesis is used. The main assumptions of the phase-lagged approach are the geometrical uniformity of the passages and the scrolling of the rotor blades in front of the stator blades as the unique source of flow temporal fluctuations. Under such conditions, there is only a phase difference between the states of passages, and computational domain may be reduced to a single passage per row.
The exit static pressure is set using a prescribed ratio between the pressure and the massflow values through the mesh exit section. This technique, in comparison with a classic simple prescription of the outlet static pressure, allows modeling the flow in the compressor stage beyond the peak of the compressor characteristic, that is necessary in order to study the limit of a compressor operating range, as pointed out by Hill IV [
In the elsA software, the prescribed ratio is respected thanks to an evolution of the outlet static pressure between two successive iterations according to the following relation:
H, C, and O topologies are employed to create a structured multiblock mesh with the commercial meshing software Numeca Autogrid. The complete mesh is composed of 25 blocks, 18 of which describing the impeller geometry and 7 the diffuser. The tip clearance is meshed with an O-H topology, in order to ensure a matching connection with the C-mesh block around the impeller blade and the H behind its trailing edge. As for the vane of the diffuser, it is rounded by an O-block.
The number of points in the impeller tip clearance height is 21.
Fillets are not included. That could be the point for some future extension of study.
The cell width at walls is set to 1
Computations with a more refined mesh in the diffuser, including 4e6 points in the diffuser passage and a coarser one with only 4e5 in the diffuser passage, have also been performed in order to evaluate the mesh dependency of the numerical flow field and in particular of boundary layer separations that occur in the diffuser. Cell width at walls and mesh in the impeller passage have been kept the same.
Table
Main mesh dimensions without upstream and downstream H-blocks.
Impeller passage | Diffuser passage coarser/normal/refined | |
---|---|---|
Meridional direction | 217 | 73/141/277 |
Tangential direction | 137 | 37/81/153 |
Height direction | 81 | 65/77/89 |
In order to ensure that the flow fields resulting from the numerical simulations are valid, performance curves computed from numerical and experimental data are compared. The numbering of the planes used to express the results is shown in Figure
Numbering of planes used for performance computation.
When comparing experimental and numerical data, the static pressure numerical values are probed at exactly the same locations than the static pressure sensors. When comparing numerical results only, with and without aspiration, for example, the values are averaged in the entire plane, surface-weighted for static quantities and mass-weighted for total ones.
All performance coefficients are calculated assuming the following hypothesis. The area-averaged static pressure values in planes are supposed to be equal to the mean value between the hub static pressure and the shroud one. That allows to calculate total quantities from static pressure measurements at endwalls. The Euler theorem is applied, assuming that the direction of the flow at the impeller inlet is perfectly axial. That is used in order to calculate the flow conditions at diffuser inlet.
One important parameter for characterizing the diffusion efficiency is the static pressure recovery coefficient of the diffuser stage. It is defined as the pressure rise through the diffuser divided by the inlet dynamic pressure:
Performance data relating to the radial vaned diffuser are plotted as functions of the ratio (
Figures
Total-to-static compressor stage pressure ratio.
Total-to-static compressor stage isentropic efficiency.
As already observed in a previous work [
Figure
Diffuser static pressure recovery coefficient.
Regarding the steady-state simulations, a sudden decrease of the diffuser static pressure recovery occurs when decreasing the massflow and suggests a stall of the diffuser vanes. With the coarser mesh, the decrease is less steep and is predicted at a lower massflow than with the two more refined meshes. However, the study of the flow topology does not reveal any significant difference in stall inception, which will be further presented, and the stall of the diffuser can be considered as mesh independent.
Even if the steady-state model is unable to predict correctly the surge massflow, flow structures that limit its stable numerical operating range can be analyzed, since convergence criteria are respected. Moreover, the existence of the boundary layer separation is also predicted by the phase-lagged model.
A common method to analyze flow separations is the topological study of the skin-friction pattern [
Since Prandtl [
Critical points in skin-friction pattern (Délery [
Attachment node
Isotropic node
Saddle point
Focus
Center
Later, Lighthill [
More recently, Surana et al. [
The four types of separation lines (Surana et al. [
In the following, the skin-friction patterns on wall surfaces in the diffuser are discussed for the operating points identified in Figure
Figure
At the nominal operating point A, a particular flow pattern is located on the pressure side surface near shroud. It is not a separation: the node
On the vane suction side, two saddles are located in the hub and shroud corner and correspond to the location of the maximal adverse pressure gradient, that is, the diffuser throat. These saddles separate the core flow from the reverse flow in the corners, which may be strengthened by the modeling of corners as sharp edges instead of smooth fillets. A further comparison should allow to specify their influence.
At maximal pressure ratio operating point B, the reattachment line
When decreasing more the massflow rate, the focus
Hub skin-friction pattern.
Skin-friction patterns on pressure and suction sides of diffuser vanes.
One should notice that the stalled topology differs from the one widely called as “corner stall,” and described in the work of Lei et al., for example [
Hub-corner stall-limiting streamlines and separation lines (Lei et al. [
In the present case, the growth of the corner separation in the hub corner is explained by the cross flow: the overdeviation of the flow near the walls at impeller outlet—see Figure
Flow conditions at diffuser inlet operating point B.
Hub-suction side corner separation 3D streamlines.
Helicity in a section across diffuser passage.
In order to confirm the existence of the boundary layer separation in the diffuser, the time-averaged flow field coming from the unsteady calculation labeled as D in Figure
Figure
Evidence of a boundary-layer separation time-averaged flow field; operating point D.
In addition, a skin-friction pattern including two focuses on both suction and hub surfaces, like the steady stalled one, is periodically observed when looking at the temporal evolution of the skin-friction pattern predicted by the unsteady simulations (Figure
Snapshot of the time-varying skin-friction pattern. Existence of a two-focus topology; operating point D.
That agreement in flow topology coming from steady and unsteady calculations gives confidence in the boundary layer separation prediction and in the fact that the compressor stable operating range is actually limited by a stall of the radial diffuser vanes.
Then, a control technique aiming at avoiding the diffuser stall and then extending the compressor operating range is developed.
In two dimensions, Miller and Chapman [
Location of suction—2D approach.
In three dimensions, main and reverse flows are separated by a three-dimensional surface, whose intersection with the solid wall is a separation line of one of the four types in Figure
Two ways for modeling the effect of the suction device on the flow have been evaluated. The first uses a boundary condition in order to set a surface massflow through the faces of the mesh cells at wall included in the suction slot. This technique is quite easy to use. The precision of the suction slot geometrical description depends on the density of the mesh, that is, the size of the faces describing the hub solid surface. The second is known as “chimera technique” and is available in the elsA software. This one allows to stack up two or more meshes, in order to model complex geometries using a basic mesh and others more refined describing the solid boundary, back and forth interpolations being done between the stacked meshes. It is now used in a wide range of applications for modeling the technological effects in turbomachinery [
In this paper, the study of the most efficient shape for suction slot is not a priority. The main aim is to explore the effect of a fluid removal on the performance of the compressor stage and on the flow topology within the diffuser. Then, the first approach for modeling the suction slots is appropriate.
The second modeling was also tested, using a common shape for slot, without fillets, in order to validate the results obtained with the boundary condition technique.
Different implementations of suction have been tested: a removal of 1 percent of the experimental surge massflow, through the hub surface, along the vane, using boundary limit condition, the same removal of 1 percent of the experimental surge massflow, at the same location, using the chimera technique, a removal of 0,3 percent of the experimental surge massflow rate through the hub surface, with a reduced slot surface compared to the previous case, in order to keep the slot surface massflow constant. The upstream edge of the suction slot is also kept at the same location, a removal of 0,3 percent of the experimental surge massflow, through the vane suction side, keeping again the slot surface massflow constant, a removal of 0,1 percent of the experimental surge massflow, through the same suction slot shape as the 0,3% hub-suction case. The slot surface massflow is then reduced, a removal of 0,1 percent of the experimental surge massflow through a smaller suction slot, so that the slot surface massflow is the same as the case of 0,3% hub suction.
Figure
Locations of suction slots and percentage of removed massflow; operating point C.
Figure
Static pressure recovery coefficient of aspirated diffuser.
Figure
Isentropic efficiency of the aspirated compressor stage.
The
The influences of the aspiration devices on the performance of the stage are predominant for operating points at lowest massflow. The main effect is the withdrawal of the drop of the diffuser static pressure recovery coefficient, the impeller working under the same conditions as without aspiration. Efficiency decrease is moderate and is mainly due to
Chimera and boundary limit conditions techniques give both similar results, ensuring the validity of exploring the effect of suction by using the boundary limit condition technique, prior to designing an optimal suction slot.
The vane suction-side suction slot does not improve the performance more than the other configurations and should then be avoided for technical reasons, due to the vane thickness.
In conclusion, the suction technique seems to be a promising way in order to control the stall of the diffuser static pressure recovery. An analysis of the skin-friction pattern within the aspirated diffuser allows to study more precisely the effect of aspiration on flow field.
Figure
Skin-friction patterns for the different aspirated diffusers.
For the 1% suction cases, the boundary layer separation in the hub-vane suction side corner region is completely removed. Even for the operating point at lowest massflow on Figure
Chimera technique gives similar results, only excepting a small vortex created at the slot downstream edge, which highlights the necessity of an optimization process of the suction slot shape; that is however not a matter of concern for now.
The 0,3% aspirated cases do not achieve a complete control of the hub-corner boundary layer separation and only diminish its extent. Hub and suction-side aspirations have the same effect on the flow field; in particular, the extent of the separation is similarly reduced.
Concerning the 0,1% aspirated cases, the hub-corner boundary layer separation is even less reduced. But the interesting thing is that, at this operating point, the two suction slots lead to a similar effect on the flow, despite the different surface massflows. The total amount of removed fluid is then thought to be the determinant parameter in this case, and not the slot surface massflow.
However, when decreasing more the stage massflow, the degeneration of the separation into a massive stall is not avoided using the smallest suction slot. Figure
Stall of the 0,1% aspirated diffuser.
A skin-friction pattern investigation within a radial vaned diffuser, aiming to highlight the flow phenomena limiting the stable operating range, has revealed the growth of a boundary layer separation in the diffuser hub-suction side corner. When reaching the numerical stability limit, that boundary layer finally led to a fully stalled flow pattern, causing a dropdown of the diffuser static pressure recovery coefficient and thus of the global compressor stage performances.
The existence of that boundary layer separation was confirmed by phase-lagged models results.
Taking into account the location of the skin-friction pattern critical points, a boundary layer suction technique was developed. The main separation saddle was located, and the suction slot was placed in its vicinity.
Several cases were tested, and an effective control of the diffuser stall has been reached with a removal of only 0,1% of the total stage massflow. From a steady-state point of view, that effective control has led to an increase by 40% of the stall-free diffuser operating range.
Comparison between the various suction configurations has not revealed any significant differences between the results obtained with a hub or a suction-side suction slot. And the amount of fluid removal was thought to be the determining factor, not the slot surface massflow.
In spite of these encouraging results, it is well known that the unsteadiness plays a leading role in case of transonic centrifugal compressors operating near surge. The aspiration technique—whose development has been based on the results of steady-state model, that stabilizes operating points for massflows lower than the experimental surge massflow—should then be implemented in unsteady calculations in order to ensure its efficiency and the increase of the operating range.
The validity of single-passage calculations should also be examined, since it may have a stabilizing effect on the numerical simulations, maybe as important as the time-dependent effects due to impeller-diffuser interaction.
The influence of the fillets on the flow topology should also be evaluated in future work, since it may modify the adequate location of the suction slot.
The major role played by the vortex created by the overincidence on diffuser vane near the shroud was also pointed out by the present study. It may then be interesting to evaluate other design techniques, that would take that observation into account. In particular, the vane of the studied diffuser is of 2D shape, with a leading edge that is just a line parallel to axis of revolution, and the hub wall is of planar shape. Studies have already shown that potential improvements could be obtained thanks to a redesign of the diffuser leading edge [
Blade pressure side, suction side
Massflow
Corrected massflow for impeller inlet at standard conditions
Impeller static pressure ratio
Static pressure recovery coefficient of the diffuser
Static pressure
Total temperature
Total pressure
Impeller rotational speed
Absolute flow angle
Compressor specific speed
Value at the
The authors would like to thank Turbomeca, Groupe Safran, which supported this research. Thanks also go to the French Aerospace laboratories, Onéra, that provides the numerical code, elsA, and to the