A Method To Estimate The Performance Map of A Centrifugal Compressor Stage
A Method To Estimate The Performance Map of A Centrifugal Compressor Stage
Introduction than those tested, and does not provide a useful reference based
on the best experience derived from many different stages.
An important requirement during the preliminary design of a
The common alternative to these approaches is to use empirical
centrifugal compressor stage is the calculation of a reliable per-
correlation based methods for predicting performance maps. Such
formance map as a guide to the expected operating flow range and
methods require fairly detailed information of the stage geometry,
the sensitivity to speed variations when the design is completed.
at least on a one-dimensional (1D) basis. Examples of 1D methods
With this information the designer can assess if the design will be
are given by Swain [1], Aungier [2], Oh et al. [3], and Cumpsty
suitable for the application. For example, it is possible to check if
[4]. The authors’ experience with such methods, however, is that
a new design will provide adequate efficiency, pressure ratio, and
they often require fairly tedious tweaking of coefficients in the
surge margin on the low speed characteristics and sufficient choke
empirical models to generate a satisfactory performance map, and
margin at high speeds.
that the characteristic curves often have unreliable shapes.
This paper presents a new engineering approach to estimate the
The approach proposed here is based on the fact that well-
performance map for such preliminary design applications. This is
designed state-of-the-art compressor stages for a particular aero-
a difficult task as the final detailed geometry (such as throat areas,
dynamic duty tend to have fundamentally similar shapes of their
blade angles, blade number, blade thicknesses, etc.) and many
efficiency performance maps. This indicates that the duty itself is
aspects of the aerodynamic design (diffusion levels, flow angles,
an excellent guide to the form of the performance map. The
etc.) have not been finalized at this point. An additional difficulty
method proposed is entirely novel in that it uses four key nondi-
is that the performance map is needed very quickly, typically dur-
mensional parameters at the design point to determine the per-
ing engineering discussions on a proposed new development. The
formance map, rather than the geometry of the stage.
most accurate methods of determining the map for competing
The method has been applied to a wide range of stage types,
design solutions, by first completing the design and then making
which are categorized here as turbocharger or process style impel-
CFD simulations or even measurements, are not suitable in this
lers, with vaned and vaneless diffusers, giving four categories.
time frame.
Turbocharger impellers (which are similar to gas turbine impel-
In some cases, the new stage may be sufficiently similar to a
lers) are relatively long and un-shrouded, have an axial inlet
stage that has previously been tested such that the map of this
designed with an inducer, and are equipped with thin blades and
stage can be quickly scaled to give a good estimate of the
splitter vanes to make them suitable for transonic flow. Process
expected performance. The scaling approach is limited in that it
style impellers in multistage compressors are shrouded and
only works if such data is available, and will not work well if the
shorter, they usually have radial inlets, have thick blades most of-
design point of the new stage is sufficiently different from previ-
ten with no splitter vanes, and are generally designed for lower
ous experience, especially if the back-sweep and work coefficient
Mach numbers. The emphasis of the paper is on the application of
have changed as this changes the slope of the flow versus work
the method to predict performance maps for turbocharger style
characteristic. In addition, the scaling of maps provides no para-
impellers with vaneless diffusers, but examples are also given of
metric description, gives no speed lines at higher or lower speeds
some other categories of stage to demonstrate the validity of the
method over many applications.
The paper is organized as follows. First, some background to the
1
Also at Institute of Thermal Turbomachinery, University of Stuttgart, Pfaffen- key nondimensional parameters and their significance in this
waldring 6, D-70569 Stuttgart, Germany. method is given. This is followed by a description of the algebraic
Contributed by the International Gas Turbine Institute (IGTI) of ASME for publi-
cation in the JOURNAL OF TURBOMACHINERY. Manuscript received September 5, 2011;
equations used to parameterize the efficiency performance maps,
final manuscript received October 17, 2011; published online November 8, 2012. which explains how these equations allow the physical features of
Editor: David Wisler. the maps to be modeled for variations of flow and speed. The
the operating range is largest and the flow at peak efficiency is 2 0 1H 31=H
/ /p
then furthest from the choke conditions. 6
g 6 B/ c / c C 7
The same information is presented in Fig. 2, but plotted relative / > /p ; ¼ ð1 GÞ þ G61 B C 77 (6)
to the flow at peak efficiency (///p). The form shown in Fig. 2 is gp 4 @ /p A 5
1
needed in this method, as in the preliminary design phase the /c
design flow is known and the choke flow is not known as the
throat areas have still to be determined. This equation has a similar ellipselike shape to the equations
Rodgers [5] gave no equations for these curves but in several pub- below the peak efficiency point and this leads to a smooth transi-
lications (Swain [6], Casey [7], and Swain [1]) some analytic equa- tion between the two arcs around the peak efficiency. The curves
tions for such curves have been described. After extensive automatically have the maximum efficiency at the flow coefficient
experimentation with these published equations, and many other for peak efficiency and have the maximum flow with a vertical
forms, a new more physically based equation structure has been characteristic at the flow coefficient for choke. The coefficients H,
developed. Clearly the equations used are not actually a real physical G, and the ratio of the flow coefficient at peak efficiency to that at
model of the losses but a particular structure has been chosen in order choke (/p//c) are not constant but also vary with the tip-speed
to reproduce important physical effects in the stage characteristics. Mach number to give narrower characteristics of different shape
In the current work the variation of stage efficiency with flow as the Mach number increases, see Fig. 2.
along a speed line is a modified form of an elliptic curve: The exponent H has a similar function to the exponent D in Eqs.
y2 1=2 (5). This exponent takes into account the fundamentally different
x2 y2 x shapes of efficiency characteristics for low Mach number and high
þ ¼ 1; ¼ 1 (4)
a2 b2 a b Mach number impellers. Low Mach number impellers tend to have
a smooth drop in efficiency ratio related to incidence losses as the
For flows below the peak efficiency point this equation is modi- flow increases above peak efficiency. A value of H near to 2, giving
fied to have a variable exponent to give the characteristic curves elliptical curves, is needed to match low-speed experimental data.
for efficiency ratio as a function of flow ratio as follows: High Mach number stages tend to have a small plateau of high effi-
" !D #1=D ciency close to the peak efficiency point and then drop much more
g /=/c sharply into choke. To match experimental data on the shape of
/ < /p ; ¼ 1 1 (5) high Mach number characteristics a value of H ¼ 3.5 is needed,
gp /p =/c
which is close to the exponent of 4 used by Swain [1].
Equations (6) includes an additional change to account for the
fact that the location of maximum flow at choke is not necessarily
at zero efficiency, but at an efficiency ratio of (1 G). High Mach
number impellers reach choke at relatively high efficiency (with
0 > G > 1), or even at peak efficiency if all sections operate at
unique incidence conditions, see Lohmberg et al. [8]. The vertical
choked characteristics start immediately below this point and not
when zero efficiency has been reached, as shown in the high
Mach number curves in Fig. 2. In contrast, low Mach number
impellers fall to an efficiency of zero well before the maximum
choke flow is reached (with G > 1), as high losses due to negative
incidence stall occur before choke, see Casey and Schlegel [9].
The flow continues to increase below the point of zero efficiency,
as also shown in the low Mach number curves of Fig. 2. In this
way G is also a function of the tip-speed Mach number.
The analysis of test data for typical turbocharger style stages indi-
cates that at high speeds a value of G ¼ 0.3 is needed—indicating
that choke occurs at an efficiency ratio of 0.7. At low speeds a value
Fig. 2 Variation of efficiency ratio with the ratio of the flow of G ¼ 2 is needed, which indicates that choke occurs well down the
coefficient relative to that at peak efficiency for a range of characteristic when the efficiency ratio reaches 1. The diagrams
tip-speed Mach numbers shown by Casey [7] include this effect but are simply elliptical at all
1 D2 qt1
/2 ¼ / (10)
p b2 q2
speed Mach number, as depicted by the dotted line. This is not an and the Euler work input coefficient can be expressed in terms of
aerodynamic effect but is an error in efficiency measurement due the inlet volume flow coefficient as
to ignoring the heat transfer effects from the turbine in the effi-
ciency determination on a turbocharger test rig; see Casey and cs /D2 tan b02
Fesich [11]. The strength of this effect depends on the amount of kEuler ¼ 1 þ (12)
u2 b2 p ½1 þ ðc 1Þc kM2 nd11
heat transfer, and for adiabatic stages there is only a small fall in d
efficiency as the Mach number is reduced below the design point.
The experimental error in the efficiency determination in these When operating at a flow coefficient near to the peak efficiency
tests is probably not better than 62%. point, backswept impellers have a fairly constant value of the
degree of reaction (around 0.6), so this and the slip factor are
Summary of the Calculation Procedure for the Efficiency. For taken as constant along the characteristic. The polytropic expo-
the calculation of efficiency on a given speed line, the first step is nent can be approximated from the stage efficiency.
to calculate the values of the variable coefficients (D, G, H, and The design value of the work coefficient includes not only the
/p//c) that correspond to this tip-speed Mach number, from Eqs. Euler work but also the disk friction work so that
(7) and (8). The peak efficiency and the flow coefficient at peak
efficiency need to be determined from the equations for the curves kdf
kd ¼ kEuler d 1 þ (13)
shown in Figs. 6 and 7. For a given inlet flow coefficient this then /
allows the efficiency to be determined from Eqs. (5) for flows
below the peak efficiency or Eqs. (6) for higher flows. where the value of kdf is a disk friction coefficient and can typi-
At this point it is useful to summarize what has been achieved cally be taken as 0.003 or calculated from the design value of the
with regard to the prediction of efficiency. The objective was to disk friction power using correlations from Daly and Nece [12].
calculate the value of the polytropic efficiency for specific values The Euler work coefficient varies linearly with flow so that from
of the independent variables (the flow coefficient and tip-speed Eq. (12) at the design point we derive a linear expression at the
Mach number). A completely new set of parametric equations has design tip-speed Mach number and then have
been developed for this purpose, based on the key nondimensional
parameters which determine the performance at the design point. cs
kEuler d ¼ 1 þ /d kkd (14)
The equations describing this functional dependency are physi- u2
cally realistic. They also include many variable coefficients
(which are a function of the tip-speed Mach number) and some The value of the coefficient in this equation, which is the slope
fixed constants, and typical values of these have been given of the Euler work coefficient curve with a variation in the flow,
above. The excellent agreement shown in Fig. 4, which is typical can be determined at the design point as
of other cases, demonstrates that the objective has been achieved
with these equations and coefficients for turbocharger style stages cs
to within an error band of roughly 62% in the efficiency, which kkd ¼ kEuler d 1 þ =/d (15)
u2
probably corresponds to the accuracy of these measurements.
The slip factor is assumed to stay constant and needs to be
Variation of Work Input Coefficient. The second part of the specified, but the coefficient for the slope changes with the Mach
functional relationship of Eq. (1) is the variation of the work input number, as can be seen in Eq. (12). An approximate relationship
coefficient with flow and tip-speed Mach number. This variation for this slope, which is related to the change in density at impeller
is based on the Euler equation following an approach described by outlet with speed follows from Eq. (12), and assumes that the
Casey and Schlegel [9]. Under the assumption that the flow has no impeller efficiency is similar to the stage efficiency so that the
swirl at inlet to the impeller, then the Euler work due to the adia- polytropic exponent can be calculated:
batic work input of the impeller can be approximately estimated
from the velocity triangles. The work is related directly with the 1
flow coefficient at impeller outlet /2 ¼ cm2 =u2 , the assumed slip kki 1 þ ðc 1Þcd kd Md2 nd 1
¼ 1 (16)
velocity cs =u2 , and the impeller outlet blade angle b02 , as follows: kkd 1 þ ðc 1Þcd kd Mi2 nd 1
cu2 cs
kEuler ¼ ¼ 1 þ /2 tan b02 (9) Equations (15) and (16) provide the slope of the work coeffi-
u2 u2
cient curve at each speed and the complete equation for the varia-
Note that the impeller outlet angle is negative for a backswept tion of the work coefficient with flow and tip-speed Mach number
impeller with the notation used here. For a constant value of the is then given as
The coefficients used for the realistic surge line shown in Fig. 9
are as follows: (/s//c)LO ¼ 0.225, (/s//c)HI ¼ 0.835, As ¼ 0.0,
Bs ¼ 1.25, and Cs ¼ 4.75.