Unclassified Ad Number: Authority
Unclassified Ad Number: Authority
Unclassified Ad Number: Authority
AD011178
CLASSIFICATION CHANGES
TO: unclassified
FROM:
confidential
LIMITATION CHANGES
TO: Approved for public release, distribution unlimited
FROM: Distribution authorized to U.S. Gov't. agencies and their contractors; Administrative/Operational Use; 08 JUN 1953. Other requests shall be referred to National Aeronautics and Space Administration, Washington, DC.
AUTHORITY
NACA notice no. 121 dtd 4 Nov 1957; NASA TR Server website
Reproduced
by
SECLJRiTY
INFCIRrMArIO=N
CONFIDENTIAL
RMoE3Dy
Copy
73
LL.
zRESEARCH
MEMORANDUM
DIFFUSION FACTOR FOR ESTIMATING LOSSES AND LIMITING BLADE LOADINGS IN AXIAL-FLOW-COMPRESSOR BLADE ELEMENTS By Seymour Lieblein, Francis C. Schwenk and Robert L. Broderick Lewis Flight Propulsion Laboratory Cleveland, Ohio
CLASSr=E
DOCU)EnT
meaftal This matertal contab.n WnotnUon offec~ng liD Nk&UO Defense of 00 United States wltta then 25. U.S.C., Secs. 7N3ad 794, the Utnamieton or revelation of wbich In sNW of the septogage bws, -I) Oy "a. uiAnnr to an unnwrithled person to proltibited
WASHINGTON
June 8, 1953
CON FIDENTIAL
NACA RM E53DO1
DIFFUSION FACTOR FOR ESTIMATING LOSSES AND LIMITING BLADE LOADINGS IN AXIAL-FLOW-COMPRESSOR BLADE ELEMENTS By Seymour Lieblein, Francis C. Schwenk, and Robert L. Broderick
A simplified limiting-blade-loading parameter for axial-flowcompressor blade elements was derived from the application of a separation criterion used in two-dimensional boundary-layer theory to a typical suction-surface velocity distribution of a compressor blade element at design angle of attack. The derived limiting diffusion factor comprised a term involving the relative velocity ratio across the blade and a term proportional to the circulation about the element, both of which can be evaluated readily from the design velocity diagram. A correlation was made for the data for NACA 65-series compressor blade sections in low-speed two-dimensional cascade between bladeelement loss at design angle of attack and several commonly used and proposed blade-loading parameters. The derived diffusion factor constituted a satisfactory limiting-loading criterion for these data. Total-pressure loss coefficient at design angle of attack increased slightlyrise in loss was factor up to a value of about 0.6, after which a sharp with diffusion observed. Loss coefficient and diffusion correlations were also made for a large number of conventional experimental axial-flow single-stage compressor rotors and stators. For the hub- and mean-radius regions of the rotors and for the hub, mean, and tip regions of the stators, no consistent variation of loss coefficient with diffusion factor was observed in the low-loss range up to values of diffusion factor of about 0.55 to 0.65, the extent of the available data. For the tip region of the rotors, however, there was a very marked and practically linear variation of loss coefficient with diffusion factor, with resulting blade-element efficiencies less than 0.90 for values of diffusion factor greater than about 0.45. Further verification of the limiting relative inlet Mach number of about 0.70 to 0.75 for subcritical values of diffusion factor was indicated for rotor and stator blades of NACA 65-series profile shapes with about 10 percent maximum thickness ratios.
, i-
.CONFIDENTIAL
I I
CONFIDEN~TIAL
NACA BM E3O
tion
INTRODUCTION
One of the principal limitations and losses in the flow across axial-flow-compressor on the rows is due to the separation of the friction boundary layers blade suction surfaces of the blades. In airfoil theory, the measure of this limiting blade loading or diffusfbn was satisfactorily established by the variation of the lift coefficient of the section. (Lift coefficient is defined as lift force per unit blade area divided by inlet or vector-mean dynamic head.) The use of lift coefficient as a measure of loss variation and limiting loading not, however, generally successful in compressor design application. Furthermore, the use of lift coefficient for the general compressor case of compressible flow with unequal axial velocities at blade-element inlet and outlet presents a difficulty in establishing a significant vector-mean condition. Consequently, various other loading parameters were tried in compressor design with varying degrees of success. These parameters were lift coefficients based on inlet or outlet conditions, solidity multiplied by lift coefficient, ratio of change in tangential velocity to axial velocity, ratio of axial velocity to tip speed, and so forth. In reference 1, Wislicenus proposes an approximate stalling coefficient for axial-flow blades derived from considerations of the stalling suction-surface pressure distribution of isolated airfoils in incompressible flow. The coefficient is developed in terms of vector-mean velocity and indicates that the stalling pressure rise is a function of the sum of the conventional lift coefficient and the over-all change in free-stream velocity across the blade row. Howell, in reference 2, develops an approximate loading parameter for compressor blade elements from the application of a separation criterion used in boundary-layer theory. This parameter involves the product of the mean-velocity lift doefficient and a function of the over-all velocity ratio. No extensive experimental investigations of these various loading parameters have as yet appeared.
*was
2$
Jseparation
-
A similar approach to the development of an improved separation criterion for axial-flow blade elements was made at the NACA Lewis laboratory in order to develop a simplified blade-loading parameter based on boundary-layer theory and independent of lift-coefficient terms. A simplified parameter, called the diffusion factor for design angle of attack, was obtained by deriving an approximate relation between a criterion used in two-dimensional incompressible turbulent boundary-layer theory (as in ref. 2) and the external velocity diagram and solidity of the blade element. The approximation was facilitated by the use of a model blade-surface velocity distribution and of the local pressure characteristics of blades in two-dimensional cascade at design angle of attack.
CONFIDENTI.AL
NACA RM E53DOI
CONFIDENTIAL
An experimental correlation between total-pressure loss coefficient at design angle of attack and a large number of proposed separation parameters is presented for the two-dimensional cascade data of reference 3. In addition, a further correlation between the loss coefficient and the derived diffusion factor is presented for the hub-, tip-, and mean-radius regions of rotors and stators of a large number of experimental singlestage axial-flow compressors. tt
SYMBOLS
The following symbols are used in this report: a ratio of distance from leading edge to point of maximum suctionsurface velocity to chord length constant in diffusion-factor equation drag coefficient lift coefficient stalling coefficient of reference 1 chord length diffusion factor constant in diffusion-factor equation force drag force boundary-layer form factor (ratio of displacement thickness to momentum thickness) angle of incidence, angle between inlet-air direction and blade mean-line direction at leading edge lift force
b CD CL
Cs
c D
*.-
d F
!D
FD H i
lW{W~fTA
p1
4
p q R
Re r s T t V V Vav Vma x x CL
CONFIDENTIAL
static pressure constant reaction or recovery ratio
momentum-thickness Reynolds number radius blade spacing total temperature static temperature relocity circumferentially averaged velocity average velocity between Vma x and V2
IACA RM E53DO1
CO
maximum velocity on suction surface distance along flow path angle of attack, angle between inlet-air direction and blade chord air angle
r' T
-adiabatic
velocity-gradient parameter of reference 2 ratio of specific heats temperature-rise efficiency boundary-layer momentum thickness static density solidity shear stress relative total-pressure loss coefficient, P2,i " P P P1
e
p a T Wi
CONF IDENTIAL
NACA RK E53DO1
Subscripts:
CONFIDENTIAL
a i
m R
S-ref
axial
!I
ideal
stator tangential
inlet
0
1
IThe
2 outlet
Exclusive of three-dimensional effects in the end regions of blade rows, the development of friction boundary layers along the surfaces of compressor blade elements is determined primarily by the surface pressure distributions and therefore by the basic geometry of the element (profile shape, canber, thickness, solidity, and angle of incidence). Therefore, a first approach to the problem of boundary-layer growth and separation can be made by considering the two-dimensional flow about blade elements. Although several techniques are available for computing boundary-layer growth in two-dimensional cascades (refs. 4 to 6), they depend on a prior knowledge of the pressure or velocity distribution about the blade section. However, accurate practical methods of computing compressible pressure distributions about arbitrary compressor cascade sections are not as yet available. In view of the general complexity and incompleteness of current boundary-layer and pressure-distribution theory, the accurate calculation of the boundary-layer development and total-pressure loss of a given blade element was not believed to be currently feasible. It was decided, therefore, as a first approach, to investigate the possibility of deriving a simplified and approximate blade-loading or separation criterion. The approximate separation criterion might then, in conjunction with an extensive experimental correlation, provide a simple and useful means of determining regions of safe or unsafe limiting design conditions for the actual compressor.
CONFIDENTIAL
CONFIDETIAL
NACA RM E53DO1
For conventional compressors, in general, differences in blade shape, boundary-layer Reynolds number, turbulence level, transition from laminar to turbulent flow, and Mach number are not very great; so that, comparatively, the most significant factors influencing the pressure distribution and the loss of a blade element are the angle of incidence and the circulation about the element. Furthermore, if attention is restricted to the design point (design angle of incidence), it would be expected that the basic velocity vector diagram of the element at design
would be the largest factor affecting the loss. The problem is then reduced to deriving a satisfactory, simplified, and approximate relation between the limiting diffusion gradient on the blade suction surface and the over-all velocity change across the .blade row as determined by the
design velocity diagram. Diffusion Factor A parameter frequently used in establishing a separation criterion in incompressible, two-dimensional, turbulent boundary-layer theory (refs. 2 and 7, for example) is given by D
2 V Re
()
where e is the boundary-layer momentum thickness, V is the freestream velocity outside the boundary layer, x is distance along the "1rctionof flow, Re is the momentum-thickness Reynolds number, and n is a constant (negative value). The significance of equation (1) can be seen from inspection of the general momentum boundary-layer equation
de " -dx
pV
_(H+2) V d e dV 2 V
where T is the viscous shear stress at the surface, H is boundarylayer form factor (ratio of displacement thickness to momentum thickness), and p is free-stream density. For configurations with fairly large adverse velocity gradients - as in the case for compressor blade suction surfaces - rapid rates of boundary-layer growth are associated primarily with the last term of equation (2) involving the velocity gradient dV/dx. In fact, for incipient separation, the shear stress T approaches zero, and the form factor H approaches an approximately constant value of about 2.0 to 2.5 (refs. 6 and 7). Equation (2) therefore becomes
dO
d,
constant
(3)
j"
CONIDEINTIAL
NACA RM E53DO1
CONFIDENTIAL
Equation (3) shows that the first factor on the right side of equation (1) is proportional to the gradient, or growth rate, of the iomentum thickness e at the point of incipient separation. In order to evaluate the importance of the second factor, Ren, in equation (1), it is assumed that for typical velocity distributions, the local value of Re is related to the initial value of Re at the
of deceleration.
~start
D
previously accelerating, the initial value of Re defines the condition of ke initial boundary layer (ref. 6). Thus the Ren term in equation (1) reflects the additional effect of the state of the initial boundary ,layer. In general, the greater the magnitude of Re, the smaller -. allowable diffusion at separation (refs. 6 and 8). he For conventional compressors - and certainly for a given cascade investigation - the range of variation of initial boundary-layer Reynolds number, as influenced by such factors as blade size, inlet Mach number, surface roughness, and turbulence level, is generally not large.
The effect of initial boundary-layer Reynolds number can thus be represented by a constant average value. A simplified separation criterion is therefore adopted as the diffusion factor: D V dV (4)
The variation of velocity over the surfaces of a typical compressor blade section set at its design angle of attack is shown in figure 1, and the element velocity diagram is shown in figure 2. By design angle of attack is meant the angle of attack at minimum loss. Diffusion of the suctionsurface velocity must occur from the peak velocity near the leading edge Vmax to the outlet velocity V2 . (It is assumed that V2 at the blade trailing edge is equal to the mean V2 considered by the velocity diagram at the downstream measuring station of the blade row.) For most conventional blade designs, the velocity generally varies approximately linearly from Vma to V2 . The diffusion gradient can then be approximated on the average by idV Vmax - V2
"V -
(1-a)cVav
where Vav is some average surface velocity between Vm and V2 . The use of equation (5) is expected to be valid for moderately nonlinear velocity variations as well, since recent unpublished results of boundary-layer analyses conducted at the NACA Langley laboratory indicate that the over-all diffusion Vmax - V2 is the significant factor affecting limiting loading.
'_5
_-L
CONFIDENTIAL
CONFIDEIAL
NACA RM E5IDOI
From examination of velocity distributions about typical compressor blade sections, it was found that the average velocity Vav in the denominator of equation (5) can be approximated by VI, so that equation (5) can be given by
1ldV
V2 \
(6)
In general, the magnitude of the maximum suction-surface velocity can be considered to be composed of several components so that, as shown in figure 1, Vmax
V + AV
+ AV + AV
(7)
where AV represents the velocity increment due to chordwise variations of the circumferentially averaged velocity (function of ratio of blade thickness to blade spacing and of inlet Mach number), and AV representi the velocity increment above the value of V that arises from the blade loading or circulation (function of profile shape and angle of attack). In terms of inlet velocity, equation (7) can be expressed as
Vmax v1
6 1l
V1
V + AV+(8)
V1
ratio and circulation as indicated in equation (8), use was made of the
characteristics of blade sections in low-speed, two-dimensional cascade. The similarity of surface pressure distributions for a given blade section in cascade and in the compressor configuration in the region of design angle of attack is demonstrated in reference 9. From examination the pressure distributions about the NACA 65-series blade sections in cascade (ref. 3), it was deduced that a functional relation can be established between the velocity ratio and a circulation parameter AVe/oV1 (proportional to lift coefficient) in the form
*of
vT
b-
+ d
(9)
where d represents the surface-velocity rise due to the effect of finite blade thickness indicated in equation (8). In general, both b and d may vary with blade shape, angle of incidence, and inlet Mach
number. If it is then assumed that equation (9) is valid for conventional compressor blades, the diffusion parameter D can be expressed as
D D (l-c +b+ dd] (10)
CONFIDEThIAL
.. ..
iP
.. .. . ....
mlpw l
.-' i !i
? !!
i
NACA RM E53DOl
CONFIDENTIAL
If it can further be assumed, for simplicity, that the average ratio of boundary-layer momentum thickness to chord length, e/c, and the location ratio, a, vary very little for conventional types of blading, then, on the average
e
(1-a)c constant and the diffusion parameter can then be reduced to the form
D=
+d
(ii)
Average values of b and d obtained from examination of the data of reference 3 at design angle of attack were about 0.4 to 0.5 and about 0.1, respectively. The design angle of attack for the cascade data was taken as the value of incidence at the midpoint of the low-loss range. In checking the trend of the effect of compressibility on the value of b and d obtained from the low-speed data, it was found that, if the local pressure coefficient and the circulation parameter both varied in the same manner with Mach number (say, the Prandtl-Glauert relation), the increases in magnitude of b and d at the design angle of attack up to Mach numbers of about 0.75 would be small. Values of b and d were consequently taken as 0.5 and 0.1, respectively, to give
+ 0.l
(12)
Since there are generally small differences in thickness and thickness distribution for commonly used blades, the last term in equation (12) will not vary much, and the design diffusion factor D for a compressor blade element is finally established as
$
*
2 U-T)
V \
Ave
+
(13)
where all velocities are taken relative to the blade. The suction-surface diffusion is thus seen to be a function of two principal factors: the over-all change in relative velocity across the element, and a term proportional to the conventional lift coefficient of the section based on inlet velocity
CONFIDENTIAL
10
CONFIDENTIAL
2 AVe
NACA RM E53DOI
CL,1
V1
cos 01 cos m
For two-dimensional flow with equal axial velocities at inlet and outlet, the diffusion factor of equation (13) becomes (see fig. 2) D '1 cos l + cos 01 cos 02o
(tn-2"ta n 0 )0 2 t n 0 -t
(14)
The general similarity between the derived design diffusion factor of equation (13) and the stalling coefficient of reference 1 is illustrated by the expansion of equation (A8) in appendix A with Vm = (V1 + V2 ) to give
2
l+v
[2
v.)
2t~ve I(
where q is a constant varying between 1 and 2. The velocity gradient parameter r of reference 2, in contrast, involves a product of the lift coefficient and the over-all velocity-ratio terms (eq. (A9), appendix A). RESULTS AND DISCUSSION Two-Dimensional Cascade The general validity of the various blade-loading parameters and separation criteria commonly used in axial-flow-compressor design was first investigated for the extensive low-speed cascade data for the NACA 65-series compressor blade sections as given in reference 3. The analysis was made by examining the variation of blade-element total-pressure loss with various loading parameters at design angle of attack. The total-pressure loss was expressed in terms of a loss coefficient 1, as defined in appendix B. For the cascade data, the total-pressure loss coefficient was computed from the drag coefficients (based on inlet conditions) given in reference 3 according to equation (B17) in appendix B. Design angle of attack for the loss evaluation was taken as the midpoint of the minimum drag range. Loading parameters considered were: (1) theoretical lift coefficient based on vector-mean velocity, (2) solidity multiplied by theoretical lift coefficient based on vector-mean velocity, (3) theoretical lift coefficient based on inlet velocity, (4) theoretical lift coefficient based on outlet velocity, (5) over-all reaction or recovery ratio, (6) relative velocity ratio, (7) velocity-gradient
S..
.CONFIDENTIAL
NACA RM E53DOI
CONFIDENTIAL
11
Aeter
4ing
fdescriptive
A second significant factor in blade-element performance is the general compressibility effect and the formation of surface shocks and boundary-layer separation. An illustrative variation of total-pressure loss coefficient with inlet Mach number for a 65-series blade element is
CONFIDENTIAL
jo
12
CONFIDENTIAL
NACA RM E55DOI
report. The value of inlet Mach number at which the losses tend to rise sharply is called the limiting inlet Mach numberk In gaueral, the limiting Mach number will vary with diffusion factor., inlet angle., and solidity. The effect of blade-chord Reynolds number on loss coefficient is
shown in figure 4(b). The curve was obtained from the data of reference 1T and recomputed in terms of the loss coefficient used in this
shown in figure 4(c) for the data of reference 3. In general, the loss is seen to increase with decreasing Reynolds number, but no sharp rise in loss is observed for the lowest values of Reynolds number obtained in the tests. Thus for the cascade blade, as long as the limiting diffusion factor and limiting inlet Mach number are not exceeded, design blade-element total-pressure loss coefficient will not be excessive for the usual range of solidity, Reynolds number, and blade shape. The exact selection of the magnitudes of these factors of diffusion, solidity, and Reynolds and Mach numbers, of course, will depend on the particular design application and such considerations as range of operation, offdesign characteristics, size, and three-dimensional effects. Compressor Rotor In view of the success of the diffusion factor in establishing a consistent loading limitation for the 65-series cascade blade, a further analysis was made to ascertain the general validity of the diffusion limit for blades in the actual compressor configuration. A correlation was obtained between measured relative total-pressure loss coefficient and calculated diffusion factor across blade elements in the hub, tip, and mean regions of a large number of experimental axial-flow single stages over a range of tip speeds. A transonic rotor, as well as conventional subsonic designs, was included in the comparison. Details of the various rotor designs and blade shapes from which the data were obtained are given in table I.
CONFIDNTAL
NACA RM E53DO
CONFIDENTIAL
13
, W
radius of the streamline flow of the stages investigated. Data for the blade-end regions were taken at points from 10 to 12 percent of the passage height from the outer wall for the tip region, and from 10 to 16 percent of the passage height from the inner wall for the hub region, depending on the particular survey program of each compressor. In all cases these tip and hub locations were outside the wall boundary regions. In order to obtain data representing the design angle of attack, plots of blade-element relative loss coefficient (and efficiency) against diffusion factor were made for each tip speed run as illustrated in figure 5, and the point of minimum loss was determined from the faired curves. Hub- and mean-radius regions. - The variation of minimum relative total-pressure loss coefficient with diffusion factor at the mean- and hub-radius regions of ten single-stage rotor configurations is presented in figure 6 for sub-limiting relative inlet Mach numbers. In the figure are also shown (dotted lines) the expected variation of loss coefficient based on cascade characteristics for average conditions of solidity and - 500 for mean; and a - 1.5, j- 400 for inlet angle (a m 1.2, hub) at inlet Mach numbers of 0.1 and 0.75. The values of loss at the low Mach number were obtained from the data of reference 3, and the values of loss at the high Mach number were obtained from the use of ratios of compressible loss to incompressible loss as determined from figure 4(b). For the hub- and mean-radius regions, no marked variation loss coefficient with diffusion factor is found, with a large accumulation of points between a loss of zero and about 0.1 up to a diffusion factor of about 0.6. It is seen that quite a few points are obtained both below and above the anticipated loss levels for blade-element flow as deduced from the cascade data. Aside from possible differences in solidity, blade shape, and Reynolds number, the data points above the upper cascade limit can be attributed to three-dimensional and unsteady-flow effects in the rotating blade row or to experimental error. In view of the many points at zero or close to zero loss, a good deal of the scatter is suspected to be the result of experimental error. Small inaccuracies in measurements and in determining the correct streamline locations at inlet and outlet can result in significant errors in the loss determination. For example, an error of 10 F in outlet temperature at a pressure ratio of 1.20 will result in computed loss values of from 0.011 to 0.050 for a true efficiency of 0.95 and from 0.042 to 0.082 for a true efficiency of 0.90. At an efficiency of 0.90 and a pressure ratio of 1.4, the spread of the loss error is somewhat smaller (0.033). In view of these considerations, the rotor data of figure 6 for hub and mean regions are interpreted as an adequate confirmation of the loss trend with diffusion factor as established from the cascade data up to
jof
CONFIDENTIAL
14
CONFIDENTIAL
NACA RM E53DOI
tor greater than about 0.55 were not obtained in the tests, so that a comparison of limiting values of diffusion factor could not be made. It appears, however, that no difficulty with limiting loading or separation should be experienced in the hub- and mean-radius regions of the rotor for design diffusion factors up to about at least 0.55, provided the limiting relative inlet Mach number for the particular blade profile is not exceeded. Further confirmation of an inlet Mach number limitation in the compressor configuration was observed for the NACA 65-series rotor blades of about 10 percent maximum thickness at the mean radius, as shown by the loss-coefficient plots of figure 7. The dashed line indicates the expected loss variation as obtained from cascade tests (ref. 11). A marked rise in loss is observed for relative inlet Mach numbers greater than about 0.7 to 0.75. The diffusion factors of the data points at the high Mach number values were less than 0.6. (These high Mach number points were not included in the plot of fig. 6(a).) A similar compressibility effect for the hub region could not be definitely established because of the generally lower levels of inlet Mach numbers in this region. It can reasonably be assumed, therefore, that, for the conventional NACA 65-series blades of 10 percent thickness ratio, Mach number difficulties will be avoided (for subcritical diffusion factors) at huband mean-radius regions of the rotor if the relative inlet Mach number does not exceed a value of about 0.7 to 0.75. For other types of blade and thickness variations, of course-, the limiting Mach number range will be different. Tip region. - The variation of minimum (design) relative totalpressure loss coefficient with diffusion factor in the tip region of the rotor is shown in figure 8. As before, the dashed lines represent the In m 600). range of loss indicated by the cascade data (a - 0.9, marked contrast to the hub- and mean-radius correlations, the tip-region plot shows a greater range of loss values and diffusion factor with a very definite and essentially linear increase in loss coefficient with diffusion factor. The trend of the loss variation below values of D of 0.3 could not be established because of insufficient data. It is expected that the loss coefficient will remain approximately constant as
maximum blade-element temperature-rise efficiency with diffusion factor is shown in figure 9. For a tip-region efficiency of 0.90, diffusion factors no greater than approximately 0.45 are indicated. The early deterioration of the flow in the rotor tip region is believed to be the result of tip losses and three-dimensional effects such as clearance leakage, casing boundary-layer growth and scrubbing, secondary flows, and radial displacement of low-energy fluid out toward
Li
tributes most to the premature rise in loss is not presently known and
C0NFIDMIIAL
NACA RM E53D01
CONFIDENTIAL
15
should constitute an important aspect of future research. In view of the fact that the flow in the tip region is not a "blade-element" or profile flow) the strength of the tip-region correlation with diffusion
factor suggests that the local (or over-all) static-pressure gradients that exist in the blade-end regions (imposed by the blade element) play an important role in the complex tip-loss phenomena. It was also observed that a correlation between loss coefficient and diffusion factor can be obtained for all loss values in the positive range of angle of incidence (increasing diffusion factor) from the region of minimum loss coefficient to values corresponding to a decrease in efficiency of about 0.10 to 0.15 (see fig. 5(a)). The resulting plot is shown in figure 10. The correlation of figure 10 results from the fact that the average slope of the variation of loss against diffusion factor in the positive incidence range (see fig. 5) is roughly equal to the slope of the original plot of minimum loss coefficient (design incidence angle) against diffusion factor shown in figure 8. The data of figure 10 may find use in estimating tip-region losses at positive incidence (low weight flow) off-design conditions. In figure 8, for the rotors with 65-series blade shapes, only the two highest loss points of rotor 1 were at relative inlet Mach numbers substantially greater than the limiting values indicated in figure 7. The effect of inlet Mach numbers above the limiting values is markedly shown in the efficiency plot of figure 9 and in the loss plot of figure 8. Theoretically, the general plots of loss coefficient against diffusion factor as derived herein would be expected to be a family of curves with relative inlet Mach number as parameter. Within the limits of experimental accuracy, a correlation of experimental loss data in terms of diffusion factor may provide a means for separating the individual loss contributions arising from blade loading (conventional surface-velocity gradients) and shock effects (boundary-layer interaction and wave formation). The scatter of the data in figures 8 and 9 is attributed primarily to differences in tip losses, in the locations of the tip-region survey stations (efficiency gradients are large in this region), and to experimental inaccuracies. For example, for a pressure ratio of 1.2 at a true efficiency of 0.85, an error in outlet temperature of 10 F results in total-pressure loss values of 0.079 to 0.138. As an alternate method which circumvents the direct use of temperature-rise efficiency) the relative total-pressure loss can be determined from the computed relative total-pressure ratio as indicated by equations (B2) and (B6). However, experimental inaccuracies are likely to be encountered in this procedure, also. It was impossible to evaluate the relative significance of the two contributing terms (the velocity ratio and the circulation parameter)
SCONFIDENTIAL
16
CONFIDENTIAL
NACA RM E53DOI
in the diffusion factor, as there was no consistent variation between the terms or dominance of one term over the other. At least for the cascade (as indicated by figs. 3(f) and (i)), both terms appear to be about equal]Z significant. The value of 0.5 taken for the constant b in the circulation term of eqjtion-(ll) appears to be quite satisfactory in view of the resulting correlation. Compressor Stator
co
jangle
e
Loss and *iffusion characteristics at design angle of attack for compressor stator blades were also investigated by examining the data from six inlet-stage stator rows, four with NACA 65-series profiles and two with double circular-arc profiles. Details of the blade designs are given in table II. Stator total-pressure loss coefficient was determined from area averages of the circumferential variation of total pressure across a blade spacing as obtained from total-pressure rakes located a short distance downstream of the stator trailing edge. The loss presented is the total-pressure defect due to the blade wake and does not include any possible drop in free-stream total pressure across the stator row. This total-pressure drop was difficult to establish with any degree of acceptable accuracy. Data for the tip region were obtained from locations at 9 to 15 percent of the passage height from the outer wall and at 13 to 16 percent away from the inner wall for the hub. Velocities were obtained from inlet and outlet surveys of pressure and angle. In most cases, only a single midspacing survey was made behind the stators. Plots of total-pressure loss coefficient and diffusion factor against of attack (or incidence) were made as illustrated in figure 11, and values of w and D were taken at the point of minimum loss (design angle of attack). The variation of minimum loss coefficient with diffusion factor for compressor stators over a range of inlet Mach numbers (rotor tip speeds) is shown in figure 12 for hub, mean, and tip regions. Minimum loss values were not attained in the hub and mean regions at some tip speeds for stator 4, because the data were insufficient to define the loss variation over a wide enough range of incidence angle. All points shown are at values of inlet Mach number below the limiting value. The dashed lines, as before, represent the estimated limits bared on cascade performance. (There is very little difference between the minimum bladeelement loss for 65-series and circular-arc blade profiles.) In all three regions, no trend of variation of loss coefficient with diffusion factor is observed, with all points falling between 0.01 and 0.05 up to the highest diffusion factors obtained (about 0.6). A cross plot of minimum loss coefficient against inlet Mach number for all three positions is shown in figure 13 for the stators with NACA 65-series profiles of 10 percent maximum thickness ratio. As in the case
CONFIDENTIAL
'I
NACA RM E53DO1 CONFIDENTIAL 17
of the rotor mean-radius region (fig. 7), a rise in loss for inlet Mach numbers greater than about 0.70 is indicated. The high Mach number points in figure 13 were at values of diffusion factor less than 0.65. On the basis of the limited information available, it is expected that very low stator losses will occur at design angle of attack, if the design diffusion factors -do not exceed about 0M and if the limiting inlet Mach number for the particular blade profile is not exceeded (about 0.70 to 0.75 for the 65-series blade of about 10 percent maximum thickness ratio).
CO
Inasmuch as the data for the preceding analyses have been obtained from a two-dimensional cascade and compressor single stages, the derived limiting values and variations will be most valid in the design of the inlet or early stages of a multistage compressor. It is not known to what extent these values can be applied to the later stages where hub-tip ratios are large and where three-dlmensional effects become more pronounced. It is also to be noted that the correlations were based on actual measured values rather than on original design values of inlet and outlet velocity. However, if the general design control and methods of flow prediction are good, the derived limiting values can be directly applied as design limits. in this respect, the inlet stages again are in a favorable position, inasmuch as design control is usually best for the early stages.
SUM4ARY OF RESULTS
A simplified limiting-blade-loading criterion called the diffusion factor for axial-flow-compressor blade elements was derived from the application of a separation criterion used in two-dimensional boundarylayer theory to the blade suction-surface velocity distribution of a compressor blade element at design angle of attack: 1. From a correlation between minimum blade-element total-pressure loss coefficient and a large number of blade-loading criteria commonly used or proposed for compressor design, it was found that, on the basis
of the data for the NACA 65-series compressor blade sections in lowspeed., two-dimensional cascade:
a. The derived diffusion factor provided a satisfactory limiting-loading parameter for the NACA 65-series compressor blade sections, in that the limiting value was clearly defined and was independent of solidity or inlet-air angle. At given solidity and inlet-air angle, the total-pressure loss coefficient at design angle of attack increased slightly with increasing diffusion factor up to a value of diffusion factor of about 0.6, after which a rapid rise in loss was observed.
-I
A
r-
18
4
CONFIDENTIAL
NACA RM E53DO1
b. The stalling coefficient of Wislicenus and the velocitygradient parameter of Howell were also adequate limiting-loading parameters. 2. From a correlation between minimum blade-element relative totalpressure loss coefficient and derived diffusion factor for data from several experimental single-stage axial-flow compressors at sublimiting inlet Mach numbers, it was found that:
(0
a. For the hub- and mean-radius regions of ten rotor blade rows, the plot of minimum relative loss coefficient against diffusion factor showed essentially no significant variation of loss coefficient with diffusion factor up to values of diffusion factor of about 0.55, which was the extent of the available data. b. In contrast to the rotor hub and mean regions, the loss correlation for the tip region of the rotors revealed a very marked and practically linear variation of loss coefficient with diffusion factor for values of diffusion factor greater than about 0.30. For a tip-region efficiency of 0.90, diffusion factors no greater than approximately 0.45 were indicated. c. For the hub-, tip-, and mean-radius regions of six stator blade rows, the plot of minimum total-pressure loss coefficient egainst diffusion factor revealed no variation of loss coefficient with diffusion factor up to values of about 0.60, which was the extent of the available stator data. d. Further verification of a limiting relative inlet Mach number of about 0.70 to 0.75 at subcritical values of diffusion factor for NACA 65-series profile shapes of about 10 percent maximum thickness ratio was indicated.
Lewis Flight Propulsion Laboratory National Advisory Committee for Aeronautics Cleveland, Ohio
CFEI
Ij
CONFIDENTIAL
I"
NACA R( E5 .DOl CONFIDENTIAL APPENDIX A LOADING PARAMETERS FOR TWO-DIMENSIONAL INCOMRMSSIBLE FLOW Lift coefficient. - The lift coefficient is defined as the ratio of L L lift force per unit area to reference dynamic head, CL 1 2)
CD
19
2(PV )ref
L
and CL
a cos
e m -F D tang m
CD
tan
C tan
2 2Va (tan Al
CL
O cos m( )ref
tan
2)
-%tan 0 m
(A2)
For the purpose of establishing design criteria., the drag term is negcoefficients are then obtained as: lected. The theoretical lift Vecta mean velocity:
2 Cos A
I
Inlet velocity: :
'. CL'
Z cos 2
i
tan i"tan p
)
C"
cos8m o
(A4)
a Cos AMA4
!
C ,Z
2con
.2 2P Cos
(t
a cos Am
anl
(A5)
COif.~TA
20 where
CONFIDENTIAL
o= tf~
v+
'
VO o2
tan
Over-all reaction or recovery ratio. - The over-all reaction or recovery ratio is defined as ratio of static-pressure rise across the blade element to inlet dynamic head, or 1 1 p1 2 (2 2p
R -1
(V2 /Vl 2
(A7)
(A8)
where
velocity-gradient
L,m
~(s
CONFIDENTrIAL
NACA RM E53D01
CONFIDENTIAL
APPENDIX B RELATIVE TOTAL-PRESSURE LOSS COEFFICIENT Compressor Cascade
21
In order to investigate and generalize the basic phenomena and variations involved in the blade-element flow, it is necessary to consider the basic loss in total pressure relative to the blade element for rotating as well as stationary blades. The relative total-pressure loss coefficient c across a blade element is defined as
--
(Bl)
where P' is relative total pressure, p is static pressure, subscripts 1 and 2 refer to blade inlet and outlet, respectively, and all values represent circumferential mass or area averages. When equation (Bl) is divided by P{, the loss coefficient becomes
P J/P1
(P2 P l
P1
1--
or P /Pi
wK
*
r-1
jJ
r)(B3)
+-[2
where M{ is the relative inlet Mach number. The ratio in the numerator of equation (B3) is referred to as the relative total-pressure recovery factor. The ideal relative total-pressure ratio is a function of the element wheel speed and the change in radius across the element, according to the equation
P1
gCONFIDENTIA
CIL
'BI,
22
where M R is
CONFIDENTIAL
NACA RM E53DOI
For stationary blade rows, the total-pressure loss coefficient is reduced in terms of absolute values to
P1 -
Ws
(B5)
C4
-1 + (l
121
For rotating blade rows, since in most cases only stationary (absolute) instrumentation can be employed downstream of rotor rows, the relative loss must be obtained from absolute measurements. The actual relative total-pressure ratio of equation (B2) can be expressed as
1~ +
2 + r 2(B6)
and, since p' = p, the ratio can further be expressed in terms of absolute total pressure as
SP,[l 1 P[
+r_.
+._r_
(M ) 2
+ r-!(Mi)
(T is total temperature
(t1
=2
Tir-, 2
2 ;I 2I/
T'
12
a.
(B7)
Then, with relative and absolute static temperature identical and with
.',:
i #;o-
.,
CONFIDENTIAL
NACA RM E53DOI
CONFIDENTIAL
23
H)
(BB)
The relative total-pressure loss coefficient for rotating blade row is thus obtained from equations (Bl), (K2), and (BB) in terms of measured Ti iF absolute pressure and temperature and relative inlet Mach number as ii'2 P P1 (T/Tr2
(B9)
r-__!l
where the ideal relative total-pressure ratio is obtained from equation (B4). The relative loss coefficient can then further be related to the blade-element efficiency through the conventional adiabatic relation
-2
ST (BlO)
Two-Dimensional Cascade The relation between blade-element total-pressure loss coefficient wand drag coefficient CD can be obtained from examination of the force and velocity diagrams for two-dimensional incompressible (lowspeed) flow. From figure 2, the drag force is given by
"FD
Fe sin Pm - Fa cos Pm
(Bll)
CONFIDENTIAL
24
CONFIDENTIAL
NACA RM E53DOI
With axial velocity equal at inlet and outlet, F a results only from the change in static pressure across the blade) which in the presence of losses is given by a= s(P2 pl
)
s[P
_-
(P1 -
2
1
or
=a
s[-
(B12)
~F
Equation (Bil)
D
tan Om-
(VI2
- V2 2 )
-P]
(B14)
tan Pme= v, =
and
9,1
-2 Va
V
0 =
a 2a v0,2
+V
V1 2
V 222
V ,1 z
V,
~F
D Cos
= 6P V with
(B15)
c/s = a,
cos Pm *-PcV
By definition,
FCD
and
6P
AP
cV1
PI
4CONFIDENTIAL
U
!NAOA RM E53DO1 that, finally CONFIDENTIAL 25
Sso
IfD
I
I
I
_CI
(B17)
W, Cos OM
~In ~dynamic
SVl
equation (BI7) the drag coefficient is related to the inlet head. For drag coefficient based on mean velocity,az V1 = 2 Vm 2
c
0s2 om 2
~so
W = C~
; , Cos Om
(B18)
I
I
REFERENCES
I. Wislicenus, George F. : Fluid Mechanics of Turoomachinery. Hill Book Co-, Inc., 1947. McGraw-
Howell, A. R. : The Present Basis of Axial Flaw Compressor Design. I. Cascade Theory and Performance. R. & M. No. 2095., British
June 1942.
3. Herrig, L. Joseph, Emery, James C., and Erwin, John R. : Systematic Two-Dimensional Cascade Tests of NACA 65-Series Compressor Blades
at Low Speeds.
NACA RM LI
IA, 1951.
Jour. Aero. Sci., vol. 19, no. 6, June 1952, pp. 404-408.
5. Loitsiansklii, Velocities. NACA TM 1r03o 1951. Airfoils in Gas Stream &ynat Subsonic L. G. : Resistance of Cascade of
6. Maskell, E. C. : Approximate Calculation of the Turbulent Boundary Layer in Two-Dimensional Incompressible Flow. Rep. No. Aero.
~Airfoils
"
26
CONFIDENTIAL
NACA RM E53D01
8. Donaldson, Coleman du P., and Lange, Roy H.: Study of the Pressure Rise Across Shock Waves Required to Separate Laminar and Turbulent Boundary Layers. NACA TN 2770, 1952. (Supersedes NACA RM L52C21.) 9. Westphal, Willard R., and Godwin, William R.: Comparison of NACA 65-Series Compressor-Blade Pressure Distributions and Performance in a Rotor and in Cascade. NACA RM L51H20, 1951. 10. Johnsen, Irvi-ng A.: Investigation of a 10-Stage Subsonic AxialI - Aerodynamic Design. NACA RM
11. Briggs, William B.: Effect of Mach Number on the Flow and Application of Compressibility Corrections in a Two-Dimensional SubsonicTransonic Compressor Cascade Having Varied Porous-Wall Suction at the Blade Tips. NACA TN 2649, 1952. 12. Moses, J. J., and Serovy, G. K.: Some Effects of Blade TrailingEdge Thickness on Performance of a Single-Stage Axial-Flow Compressor. NACA RM E51F28, 1951. 13. Standahar, Raymond M., and Serovy, George K.: Some Effects of Changing Solidity by Varying the Number of Blades on Performance of an Axial-Flow Compressor Stage. NACA RM E52A31, 1952. 14. Mahoney, John J., Dugan, Paul D., Budinger, Raymond E., and Goelzer, H. Fred: Investigation of Blade-Row Flow Distributions in AxialFlow-Compressor Stage Consisting of Guide Vanes and Rotor-Blade Row. NACA RM E50G121 1950.
15. Lieblein, Seymour, Lewis, George W., Jr., and Sandercock, Donald M.: Experimental Investigation of an Axial-Flow Compressor Inlet Stage
Operating at Transonic Relative Inlet Mach Numbers. I - Over-All
Performance of Stage with Transonic Rotor and Subsonic Stators up to Rotor Relative Inlet Mach Number of 1.1. NACA RM E52A24, 1952. 16. Jackson, Robert J.: Effects on the Weight-Flow Range and Efficiency of a Typical Axial-Flow Compressor Inlet Stage that Result from the Use of a Decreased Blade Camber or Decreased Guide-Vane Turning. NACA RM E52G02, 1952.
-I~-
CONFIDENTIAL
28
CONFIDENTIAL
NACA PM E53D01
.. .. '..
m0,.
.; ..
.m
. o
. .. + .. +
*+ , 11
.. . ... . .
s.
0 ,
. ...
-,4 +
4
0
,
+
S.. ',.
. .
,
4,
. 4+ 4 .. .
.. + . .. . . .
0. ._40 4( ~444.
'k,
'),+a
0 4+ 0 4 ) . .,, .
. .0a
+ .. 4 ... ... )B 4 ,
ad
(:
. ..... .
40 4. .
..
..
6 40
4~~~ 4. +41 4 4
60
4) 5)
5
-
5
.'4
6
3. 0 .4A
S.
, .-
-'o -1
.++
"
, 4, -.
.-
'-4
-0 f)
-48
d5 + .
....
00
.V
+t
j
+
~
+
(441
.a
.
' '
..
t +,
C NF
,
i .i
, N, A
+
.
+ P. . + .
.~
. .
. I . ' .. + ++ +
ill
;SI'
+
d
__
i++,
;+
:
_+
,.+ ,.,.,',<+',.+,,+
+' m++
"+ ,.N
FI-D . ..
IAL
t~ ..l.! +-++--- T + , + i 7 rt i+
. ..... ..
NACA RM E53DO1
-
CONFIDENTIAL
29
I
0 0 Si i 4
o
3
0-C)
C6
;-1
4'4
C.
04
',04
434
-04
*o
4--
i :
CONFIDENTIAL
30
CONFIDENTIL
NACA RM E53D01
H
_____
H~H (v)--surface
__
Upper (suction)
_____ _____ __________
___________
_____
V-K_
I
Lover (pressure)
surface
Figure 1.- Typical distribution of velocity along surfaces of compressor-cascade blade section in range of design angle of attack.
CONFIDENTIAL
NACA RM E53D)OI
CONFIDENMiAL
31
II
N*N
7i~
~~(d
C0 to
4Y-)
0
>-4
CDH
I-F-3
>-C
rZ4
H4
ca)
CONYIMT~EIAL
I32
CONFIDENTIAL
NACA RM E53D03.
..
0
.04~1
d niha&eg
7 t131,ddeg
I~
inlet-air anle
.04
.....
00
-4--II
.4
.d1.2
1.6
2.0
.4
.8
1.2
1.6
2.0
Figure 35.- Variation of total-pressure loss coefficient at design angle of attack with various blade-loading parameters for NACA 65-series blade sections In two-dimensional low-sveed cascade.
CONFIDENTIAL
NACA PM E53DO1
COM~IDED1TIAL
33
.04
Inlet-air angle,
01, deg
-
.02;
-51
.0
- - -------
Solidity,
I-'-U'1
~.04--*1.0
I 0
0.5
.02
l iI-v
.04-
.4 CLil
(c)
.8
1.2
.4
.8
1.2 CL,
1.6
2.0
2.4
Figure 3. - Continued. Variation or total-preaure loss coefficient at design angle of attack with various blade-loading parameters for NACA 85-series blade sections In two-
1dimensional
low-speed cascade.
CONFIDENIA
34
CONFIDENTIAL
NACA RM E53D01
Y W
45
3 -
s1Iiity, 1
4_ '-" ------
4 . . -
I"
. ..
..
jl0
10
fi
00
f.04
'
oG... t .... .
.2
.4
.L
.0
L
.8 0
LLJ___
.2
__
_____
.4
.6
R -
(V2/Vl)
Figure 3. - Continued. Variation or total-pressure loss coefficient at design angle of attack with various blade-loading parameters for NACA 65-series blade sections In two-dimensional low-speed cascade.
CONFIDENTIAL
NACA RM E53fl01
CONFIDEN~TIAL
35
IQ
.4Inlet-r ge,
~Inltir
0 C 1.2b
Y<zgle,
-----
LO
45
4.5
0.
40
-
I
N CA
o
P/ 4
frfr h tligcefceto
(e.
ar
(g0eoiygain
ene2(e.rfeec
aaee
Bd
FiueS
of toa-.~e blde-ladin
p-aameers
oscefiin or
bld
etosI
w-ieninllwsedcsae
CONFIDENTIA
36
CONFIDENTIAL
NACA RM E53D02.
---~
.02
_
angle~~
_____
6,a
_t
0 _T
60Solidity,
7
.2
j *j-1.-7727777
.0
A0 -
.2.4
6.8
0 Diffusion factor, D
.2
.6
cascade.
CONFIDENTIA
NACA RM E53DOJ.
CONFIZENTIAL
37
.04 X)
.02
deg
60
..
o1.6 Solidity,
2.0 o
(a) Solidity for d Iffu Ion factor of 0.5 (fig. 3(t)). .1?
.04
.0
.4
.6-
TC
--
---
T_
()Ilt
Inletanle ahnubrfrAA
60
l.017urbulence
150
200
250
S00
350
400
450
3 500X10'
4
-,
CONFIDE WUI
38 1>
CONFIDENTIAL
NACA BM E53fl01
'-4
.8
E-
4)
110
____o__
ca
1214_______
_________
S.
4552810 83
Difuio fatr
Fiur
ris
5.
efiiec
Exa04l
an1214e
of vaiaio
ntpreino
oa-rssr
oo
of bld-lmn
oscofiin
rtr9.Br
temeraure
wihdfuinfco
niaemnmmo
aiu
COFDETA
aus
NACA RM E53D01
CONFIDENTIAL
39
Rotor
I 0
ber,
....
-- -
Typical casca.e
Mean-raiu, region.
i
4)
Fi~MathQ
nu... . .. 6 .
:H - :eln enr a rtotal-pressure 10s8 coefficient with diffuI_ sinfactor for axial-flow single-stage rotors. ____ _____
RotorF
ReaiveusInle .2 A 4
-
Macr nubrM I c
.0
.2
.4
.6
.8
1.0
Figure 7. - Variation of minimum (design) blade-element relative total-pressure lose coefficient with relative Inlet Mach number at ean radius for 65-series rotor blades.
CONFIDENTIAL
40 .3
00
CONFIDMn IAL
NACA RM E53D01
Rotor
o
Is
.2
,-' Or > o
04'
.f-b. ...
0
G-4)o
i...al.
inlet-Mach number, _
'
variation
V
/
. .6
.l0
.8
.2
-
Diffusion factor, D Figure 8. Variation of minimum (design) blaelement relative total-pressure loss coefficient
1.0
0 S 1
-1=
Rotor!..-. 4
0
Y.4X
2
b
C
}
4)
*4
6 71
8
0
---
~~
W
04
$4
.9
.8
3ifs~
9
atr"
f
A:
'10
Fiue9.-Vr taiumbae-lmn
j
F
~~.7
0
A~i
.6
.8
Figure 9.
of maximum blade-element temperature-rise efficiency with diffusion factor in tip region of axial-flow singlestage rotors. C01IDMIfAL
NACA RM E53D01
CONFIDENTIAL
41
I.5
.5 JT F
1 2 3 4 5
1
__
Rotor
0
0
19 . 4
'.4,
/V
0 A
/
V
-_ _ _
/
_ _
v
4 V
7 8
10
-- - 0
f0 V
__
o
m.3 -
.... 0
,/
/1
I
o
0
0
0 .8
00
4)W
-_efiiec
of a
ot
i0._
..
4.6
.8 Diffusion factor, Dl
1.0
1.2
of blade-element relative total-pressure Figure 10. loss coefficient with diffusion factor in tip region of axialflow single-stage rotors. All data points in positive incidence range from region of minimum loss to region of drop in efficiency of about 0.10.
-Variation
CONFIDENTIAL
42
CONFIDEN~TIAL
NACA RM E53D01
Rotor tiprspeed,__
___
.4-
100
0.)
00
.2
.127
....
0 9 .08
PO -
V-4 W
ca 4 .4
4J 4)
.04
-12
-8
Angle of incidence, i, deg Figure 11. - Example of variation of wake totalpressure loss coefficient and diffusion factor with angle of incidence in hub region of stator (stator 1). Bare indicate points of minimum loss.
CONFIDENTPIAL
NACA RM E53D01
CONFIDENTIAL
43
Stator 0
---
-- --
K> 4 2
0 3 5 6
-- l
- -
Inlet1
Mach .1 number, M1
1 0.75.-i0 .1o
v
ITypical _variation
cascade
o
.I0-
-K
0
Z
3
4 Typical cascade
oo
variation
.2
.8
1.0
Figure 13. - Variation of minimum wake total-pressure loss coefficient with inlet Mach number at hub, mean, and tip regions of stators with 65-series profiles. CONFIDENTIAL
I[
rr
Z 7C
r4 1 z
I f4
94 V2' *a .8 . :to 64
4 1
0~~
16.
402--D
41
.A
0 I e
-4
04
;a0
F~~
JR4I
ufl
z2
IN
Ic 0
to~Z
0.
JJ.41
esS 0
-0
19
01.2
zI
cl
69
o4
-z
U1 0
41 B v 0a0 .0 p,
02jcs
~2
.2.