Asme PTC 10 PDF
Asme PTC 10 PDF
Asme PTC 10 PDF
ADOPTION NOTICE
ASME-PTC10, "Compressors and Exhausters,'' was adopted on
(DoD).
October 3 , 1994 for use by the Department of Defense
Proposed changes by DoD activities must
be submitted to the
DoD Adopting Activity: Director, US Army Mobility
Technology Center/Belvoir, ATTN: AMSTA-RBES, Fort Belvoir,
VA 22060-5606. DoD activities may obtain copies of this
standard from the Standardization Document Order
Desk, 700
Robbins Avenue, Building 4D, Philadelphia, PA 19111-5094.
The private sector and other Government agencies
may
purchase copies from the
American Society of Mechanical
Engineers, 345 East 47th Street, New York, NY 10017.
Custodians:
Army - ME
Navy - YD-1
Air Force - 99
Adopting Activity
Army - ME
FSC 4310
ASME
PTC 10-1997
Performance
Test Code on
Lompressors
and Exhausters
This document will be revised when the Society approves the issuance
edition. There will be no addenda issued to ASME PTC 10-1997.
of a n e w
PTC 10-1997 is being issued with an automatic subscription service to the interpretations that will be issued to it up to the publication of the next edition.
FOREWORD
(This Foreword is not a part of ASME PTC 10-1997.)
PTC 10 was last revised in 1965 andit has been reaffirmed many times
in the intervening
period. The PTC 1O Committee has been in various states of activity for approximately
the past20 years. During that time the Codehas been completely rewrittento be far more
explanatory in nature.
The performance testing of compressors is complicated by the need in virtually every
case to consider and make correction for the differences between the test and specified
conditions.Thetechniques used to do so arebased upontherules of fluid-dynamic
similarity. Some familiarity with this fundamental technique will be a significant aid to
the users of PTC IO.
Compressors and exhausters come in all sorts of configurations. A very simple case is
a single section compressor with one impeller, and single inlet and outlet flanges. Many
morecomplex arrangementsexist with multiple inlets,outlets,impellers,sections,
intercoolersandside seams. Typical gases handled areair, its constituents,andvarious
hydrocarbons. Tests are commonly run in the shop or in the field, at speeds equal to or
different from the specified speed, and with the specified or a substitute gas. In order to
handle this vast array ofpossibilities PTC 10 reduces the problemto the simplest element,
the section, and provides the instructionsfor combining multiple sections to compute the
overall results.
Uncertainty analysis can play a very important role in compressortesting, from the
design of the test to interpretation of the test results. In all but the very simplest of cases
the development of an analytic formulation, .e., in simple equation form, for overall
uncertainty computation is formidable. The test uncertainty will always be increasingly
more complex to evaluate with the complexity of the compressor configuration, and by
the very nature of the test will be a function of the performance curves.
The modern personal computer is readily capable of completing the calculations required. The Committee developed software and used it to perform both the basic code
calculations and uncertainty analysis computationsfor
widea range of possible compressor
configurations.
This Code was approved by
the PTC 1O Committee on January 18,1991.
It was approved
and adopted by the Council as a standard practice of the Society by action of the Board
on Performance Test Codes on October 14, 1996. It was also approved as an American
National Standard by the ANSI Board of Standards Review on April 22, 1997.
iii
STD-ASME
NOTICE
All PerformanceTestCodesMUST
adhere to therequirements ofPTC 1, GENERAL
INSTRUCTIONS. The following information
is based on that document and
is included
here for emphasis and for the convenience of the user of this Code. It is expected that the
Code user is fully cognizant of Parts I and III of PTC I and has read them prior to applying
this Code.
ASME Performance Test Codes provide test procedures which yield results of the highest
level of accuracy consistent with the best engineering knowledge and practice currently
available. They were developedby balanced committees representing all concerned interests.
Theyspecifyprocedures,instrumentation,equipmentoperatingrequirements,calculation
methods, and uncertainty analysis.
When tests are run in accordance with thisCode, the test results themselves, without adjustment for uncertainty, yield the best available indication of the actual performance of the
tested equipment. ASME Performance Test Codes do not specify means to compare those
it is recommended that the partiesa commercial
to
results to contractual guarantees. Therefore,
test agree before starting the test and preferably before signing the contract on the
method to be used for comparing the test results to the contractual guarantees. It is beyond
the scope of any code to determine or interpret how such comparisons shall be made.
Approved by Letter Ballot #95-1 and B E C Administrative Meeting of March 13-14, 1995
IV
PERSONNELOFPERFORMANCE
OFFICERS
Gordon J. Gerber, Chair
Richard J. Gross, ViceChair
jack H. Karian, Secretary
COMMITTEEPERSONNEL
Helmut B. Baranek, Public Service Electric & GasCompany
John J. Dwyer, Consultant
Gordon J. Gerber, Praxair
Richard J. Gross, The University of Akron
Jack H. Karian, ASME
Robert E. Lawrence, Consultant
Jack A. Lock, LockEngineering
Vincent J.Polignano, IMO Delaval
Frank H. Rassmann, ElliottCompany
Norman A. Samurin, DresserRandCompany
Joseph A. Silvaggio,Jr., Alternate to Polignano, IMO Delaval
S T D - A S M E P T C L O - E N G L L777 m 0 7 5 7 b 7 0 Ob05427 A T 3 m
D. R. Keyser, Chair
P. M. Cerhart, Vice Chair
W. O. Hays, Secretary
C O M M l l T E E PERSONNEL
R. P. Allen
R. L. Bannister
B. Bornstein
J.M. Burns
J. R. Friedman
G. J. Gerber
P. M. Gerhart
R. S. Hecklinger
R. W. Henry
D. R. Keyser
S. Korellis
J.
J. W. Milton
G. H. Mittendorf, ]r.
S. P. Nuspl
R. P. Perkins
A. L. Plumley
vi
S. B. Scharp
J. Siegmund
J. A. Silvaggio, Jr.
R. E. Sommerlad
W. G. Steele, Jr.
J. C. Westcott
J. G. Yost
S T D - A S M E P T C LO-ENGL
CONTENTS
Foreword . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
CommitteeRoster ................................................
BoardRoster ....................................................
Section
1
2
3
4
5
6
Figures
3.1
3.2
3.3
3.4
3.5
3.6
3.7
4.1
4.2
4.3
4.4
4.5
4.6
4.7
4.8
4.9
4.1O
4.1 1
4.1 2
5.1
Objectandscope ........................................
Definitions
and
Description
of Terms ..........................
Guiding Principles ........................................
Instruments
and
Methods
of Measurement ......................
Computation
of
Results ....................................
ReportofTest ...........................................
.........................
iii
V
vi
1
3
11
23
39
55
14
16
18
19
20
21
22
24
24
25
25
26
26
27
29
32
33
33
35
49
Tables
3.1
Permissible
Deviation
From
Specified
Operating
Conditions
for
Type1 Tests...........................................
12
3.2
3.3
3.4
5.1
5.2
5.3
5.4
b7b W
......................................
............................................
..........................
...........................
........................
.................
Nonmandatory Appendices
AUse
of Total Pressureand TotalTemperature to Define Compressor
Performance
B
Properties of Gas Mixtures ..................................
Sample
C
Calculations
C.l
Type 1 Test for a Centrifugal Compressor Using anIdeal Gas
C.2Type2Test
for a CentrifugalCompressorUsinganIdealGas
........
C.3
Ideal Gas Application to Selection of TestSpeed and TestGas and
Methods of PowerEvaluation
C.4
Treatment of Bracketed TestPoints
C.5
Selection of aTestGas for aType2Test Using Ideal and RealGas
Equations.............................................
C.6
Type 2Test Using RealGasEquations forDataReduction
C.7
Treatment of a Two SectionCompressor With ExternallyPiped
Intercoolers,CondensateRemoval ..........................
C.8
Application of Uncertainty
Analysis ...........................
D
References
E
Rationale
for
Calculation Methods ............................
F
Reynolds Number Correction
G
Refined
Methods
for
Calculating
Total
Conditions
H
SIUnits ................................................
..........................................
......................................
........
..............................
............................
..........
..............................................
................................
................
viii
12
13
14
40
41
48
50
59
61
63
65
85
109
119
123
139
151
159
165
167
183
185
187
~~
S T D - A S M E P T C 10-ENGL
PTC
~~
COMPRESSORS
ASMEAND EXHAUSTERS
SECTION 1
1.1
10-1 997
- OBJECTANDSCOPE
nally pipedintercoolersand
for compressors with
interstage side loadinletsoroutlets.
Internally cooled compressorsare includedprovided that test conditions are held nearly identical
to specifiedconditions.
Compressors, as thenameimplies,areusually
intended to produceconsiderabledensitychange
as aresultof
the compressionprocess.Fansare
normally considered to be air or gas moving devices
and are characterized by minimal densitychange.
A distinction betweenthe
two at timesmaybe
unclear. As a veryroughguide,either
PTC 10 or
PTC 11 maybeused for machines falling into the
approximatepressure ratio rangeof 1.05 to 1.2.
Themethodsof
PTC 10, which provide for the
pronouncedeffectsofdensitychangeduringcompression,have no theoretical lower limit. However,
practical considerations regarding achievable accuracy become important in attempting to apply PTC
10 to devicescommonlyclassified
as fans.
For
example,the low temperatureriseassociated with
fans may lead to large uncertainty in power requirement if theheatbalancemethod
i s chosen.Fans
also may require traversing techniques for flow and
gas state measurementsdue to the inlet and discharge
ducting systems employed. Refer to PTC 11 on Fans
for further information.
OBJECT
Theobjectof
this Code is to providea
test
procedure to determine the thermodynamic performance of an axial or centrifugal compressororexhauster doing work on a gas of known or measurable
propertiesunderspecifiedconditions.
This Code is written to provide explicit test procedures which will yield the highest level of accuracy
consistent with the best engineering knowledge and
practice currently available. Nonetheless, no single
universal value of the uncertainty
is, or should be,
expected to apply to everytest.The
uncertainty
associated with any individual PTC 10 test will
dependupon practical choicesmade in terms of
instrumentationandmethodology.
Rulesare provided to estimate the uncertainty for individual tests.
1.2
SCOPE
1.3
BY THIS
EQUIPMENTNOTCOVERED
CODE
1.4
TYPES OF TESTS
ThisCodecontainsprovisions
for two different
types of tests. A Type 1 test must be conducted on
the specified gas with a limited deviation between
test and specifiedoperatingconditions.
A Type 2
test permitsthe use ofasubstitutetest
gas and
extends the permissible deviations between test and
specifiedoperatingconditions.
1.5
PERFORMANCE
RELATION
TO
GUARANTEE
1.7
1.8
1.6
ALTERNATE PROCEDURES
REFERENCES
Unlessotherwisespecified,references
to other
Codes refer to ASME Performance Test Codes. Literaturereferencesareshown
in Appendix D.
INSTRUCTIONS
~~
SECTION 2
- DEFINITIONSANDDESCRIPTION
OF TERMS
2.1
Description
Units
A
a
ft2
b
C
C
C
CP
CV
e
f
gc
H
HR
h
h,
j
K
k
log
In
MW
Mm
M
m
m
N
n
n
ns
P
P
Pv
Qext
Om
ftlsec
ft
dimensionless
Btu/lbm mole "R
Btu/lbm "R
Btu/lbm "R
Btu/lbm "R
in.
In.
dimensionless
dimensionless
Ibm ft/lbf sec2
Btu/lbm-mole
Ibm H20/lbm dry air
Btu/lbm
Btu/hr f t 2 "R
e
ft Ibf/Btu
dimensionless
dimensionless
dimensionless
dimensionless
Ibmllbmole
dimensionless
dimensionless
dimensionless
Ibm
rPm
dimensionless
lb mole
dimensionless
hP
psia
PSi
Btu/min
Btu/min
~~~
0,
QSl
R
RA, RB, RC
Re
Rem
RH
RP
Rt
'r
'P
'9
rt
r"
S
Sc
S
T
t
U
U
V
V
W
W
X
X
Y
Y
Btu/min
Btulmin
ft3/min
ft IbWlbm . "R
dimensionless
dimensionless
dimensionless
percentage
dimensionless
dimensionless
dimensionless
dimensionless
dimensionless
dimensionless
dimensionless
dimensionless
Btullbmmole
fi2
Btu/lbm "R
"R
"F
Btu/lbm
filsec
ftlsec
ft3/lbm
ft Ibf/lbm
Ibm/m in
dimensionless
dimensionless
dimensionless
ft Ibfllbm
dimensionless
144 pv = ZRT
rl
P
EL in
PP
PS
Y
c
T
E
R
d
dimensionless
dimensionless
dimensionless
dimensionless
Ibm/ft sec
dimensionless
dimensionless
dimensionless
ft2/sec
Ibm/ft3
dimensionless
I bf-ft
in.
dimensionless
dimensionless
- "R
L997 D 0 7 5 9 b 7 0 Ob05434T33
COMPRESSORS A N D EXHAUSTERS
Subscripts
2.2.2
Gage
Pressure.
The
gage
pressure
i s that
pressure which i s measured directly with the existing
barometric pressure as the zero basereference.
Ambient
a,b,c,j Component of gas mixture (Appendix B
only)
Average
av
C
Casing
Correction
corr
crit.
Fluids critical point value
d
Compressor discharge conditions
Dry air
da
db
des
dg
sh
Shaft
g
hb
i
lu
Id
m
P
rotor
Dry-bulb
Design
Dry gas
Gas
Heat balance
Compressor inlet conditions
Leakage upstream
Leakage downstream
Gas mixture
Polytropic
Flow location reference
Isentropic
Specified conditions
su
sidestream upstream
sd
sidestream downstream
SV
Saturated vapor
Test conditions
t
Wb
Wet-bulb
1, l n Upstream of fluid meter
2, 2n Downstream or at throat of fluid meter
(Y
Compressor inlet conditions (static,
Appendix A only)
Compressor
discharge conditions (static,
Y
Appendix A only)
static Static
meas. Measured
SP
Superscripts
(
()
2.3
2.2
is
TEMPERATURES
2.3.1AbsoluteTemperature.
The absolute temperature is thetemperaturemeasuredaboveabsolute
zero. It is stated in degreesRankineorKelvin.The
Rankinetemperature i s the Fahrenheittemperature
plus 459.67 and the Kelvin temperature is the Celsius
temperature plus 273.1 5.
PRESSURES
2.2.1AbsolutePressure.
Theabsolutepressure
the pressuremeasuredaboveaperfectvacuum.
S T D - A S M E P T C LO-ENGL 1977 9
COMPRESSORS AND EXHAUSTERS
2.3.2StaticTemperature.
The statictemperature
is the temperature determined in such a way that no
effect is produced by the velocity of the flowing fluid.
2.4.4AbsoluteViscosity.
Absolute viscosity is that
property of any fluid which tends to resist a shearing
force.
2.4.5KinematicViscosity.
The kinematic viscosity
of a fluid is theabsoluteviscosity divided by the
fluid density.
2.4.6 Specific Heat at Constant Pressure. The specific heat at constant pressure, (c,) = (dh/aT), is the
change in enthalpy with respect to temperature at
a constantpressure.
2.4.7Specific
Heat at Constant Volume. Thespecific heat at constant volume, (c,,) = (au/aT), is the
change in internal energy with respect to temperature
at a constant specific volume.
2.3.5 Inlet Total Temperature. The inlet total temperature is the absolute total temperature that exists
at the inlet measuring station(see para. 4.7.7). Unless
specifically stated otherwise, this i s the compressor
inlet temperatureused in this Code.
AcousticVelocity(SonicVelocity).
A pressure wave or acoustic wave of infinitesimal amplitude is described by an adiabaticandreversible
(isentropic) process. The corresponding acoustic velocity for suchwaves in any medium is given by:
2.4.9
2.3.6 Inlet Static Temperature. The inlet static temperature is the absolute static temperature that exists
at the inlet measuring station.
($)
2.3.7DischargeTotalTemperature.
Thedischarge
total temperature is the absolute total temperature
that exists at the dischargemeasuringstation
(see
para. 4.7.8). Unless specifically statedotherwise,
this i s the compressor discharge temperatureas used
in this Code.
2.4.10Fluid Mach Number. The Fluid Mach number is the ratio of fluid velocity to acoustic velocity.
2.5
2.4
a2
MACHINE CHARACTERISTICS
2.5.1Capacity.
Thecapacity of a compressor is
the rate of flow which is determined by delivered
mass flow rate divided byinlet total density. For
an exhauster it is determined by the inlet mass flow
rate dividedbyinlet total density. Forsidestream
machines, this definition must be applied to individualsections.
2.4.1
2.5.3Pressure
Ratio. Pressure ratio i s the ratio of
the absolute discharge total pressure to the absolute
inlet total pressure.
S T D - A S M E P T C 1 0 - E N G L 1997
0 7 5 9 b 7 0 Ob05V3b AOb
streams enteringandleavingcanbequantitatively
defined as well as the powerinput and heat exchange
by conduction and radiation. Such a region can be
considered to be in equilibrium for both a mass
andenergybalance.
2.5.5Temperature
Rise. Temperaturerise
is the
differencebetweenthedischarge
total temperature
andthe inlet total temperature.
2.5.8Specific
Volume Ratio. Thespecificvolume
ratio is the ratio of inlet specific volume to discharge
specific volume.
2.5.15ChokePoint.
Thechoke point is the point
where the machine is run at a given speed and the
flow is increased until maximum capacity is attained.
2.6
2.6.3PolytropicCompression.
Polytropic compression is a reversiblecompressionprocessbetween
the inlet total pressure andtemperatureandthe
discharge total pressureandtemperature.The total
pressures and temperatures are used to account for
thecompressionofthe
gas and thechange
in
the kinetic energy of the gas.Thechange
in the
gravitational potential energy is assumed negligible.
The polytropic process follows a path such that the
polytropic exponent is constant during the process.
2.6.4
Polytropic
Work (Head). Polytropic work
(head) is thereversible work required to compress
a unit mass of gas by a polytropic process from the
inlet total pressure and temperature to the discharge
total pressureandtemperature.
2.6.8
Polytropic
Efficiency.
The polytropic efficiency is the ratio of the polytropic work to the gas
work.
2.7.3Specified
Operating Conditions. Thespecified operating conditions arethose conditions for
which the compressor performance is to be determined. Refer to paras.6.2.3 and 6.2.4.
2.6.9ShaftPower(BrakePower).
The shaft power
(brake power) is the power deliveredto the compressor shaft. It is the gas power plusthemechanical
losses in thecompressor.
2.6.10Isentropic
Work Coefficient. The isentropic
work coefficient is thedimensionless ratio of the
isentropic workto the sum of the squares ofthe
blade tip speeds of all stages in a givensection.
2.7.5Equivalence.
The specifiedoperating conditions and the test operating conditions, for the purpose of this Code, are said to demonstrate equivalence when, for the sameflow coefficient the ratiosof
the three dimensionless parameters (specific volume
ratio, Machine Machnumber, and Machine Reynolds
number) fall within the limits prescribed in Table 3.2.
2.7.6Raw
Data. Raw data is the recorded observation of an instrument taken during the testrun.
2.6.13 WorkInput
Coefficient. The work input
coefficient is the dimensionless ratio of the enthalpy
rise to the sum of the squares of the tip speeds of
all stages in a givensection.
2.7.7Reading.
A reading i s theaverageofthe
corrected individual observations (raw data) at any
givenmeasurementstation.
2.7.8TestPoint.
Thetest point consists of three
or more readings that have been averaged and
fall
within thepermissiblespecified fluctuation.
2.6.14
Total
WorkInput
Coefficient. The total
work input coefficient is the dimensionless ratio of
the total work input to the gas to the sum of the
squares oftheblade tip speeds of all stages in a
givensection.
2.7
2.7.9
Fluctuation.
The fluctuation of a specific
measurement is defined as the highest reading minus
thelowestreading
divided by the average of all
readingsexpressed as a percent.
MISCELLANEOUS
2.7.1FluidReynolds
Number. The Fluid Reynolds
number is the Reynolds number for the gas flow in
a pipe. It is defined by theequation Re = VD/v,
where the velocity, characteristic length, and static
kinematic viscosity are to be used as follows: velocity
V is the average velocity at the pressure measuring
2.8INTERPRETATION
2.8.1 Certainvaluesforthermodynamicstateand
mass flow rateareused in the computation of the
dimensionless performance parameters M, Re, r,, 4,
P,,, pi, T,, and
Unlessotherwise specifically
s1.
OF SUBSCRIPTS
L997
0 7 5 9 b 7 0 Ob0543B bA9
COMPRESSORS A N D EXHAUSTERS
SECTION 3
3.1
- GUIDING
3.2
TYPES OF TESTS
ThisCodedefines
two typesof test which are
based on the deviations between test and specified
operating conditions.
3.2.1 Type 1 testsare conducted with the specified
gas at or very near the specified operating conditions.
Deviations in the specified gas and operating conditions are subject to the limitations imposed by Table
3.1. These limitations are subject to thefurther
restriction that their individual and combined effects
shall not exceedthe limits ofTable 3.2.
3.1.7 Thenumberoftestpersonnelshouldbe
sufficient to assure a careful and orderly observation
ofall instruments with timebetweenobservations
to check for indications of error in instrumentsor
observations.
3.3
LIMITATIONS
PRINCIPLES
10-1997
TABLE 3.1
PERMISSIBLEDEVIATION FROM SPECIFIEDOPERATINGCONDITIONS
TYPE 1 TESTS
FOR
Permissible
Symbol
Variable
Inlet pressure
Inlet temperature
S P d
Molecular weight
Cooling temperature
difference
Coolant flow rate
Capacity
596
8%
psia
Pi
Ji
OR
rPm
Ibm/lbmole
MW
2%
2%
5%
O R
gal/min
ft3/m in
4i
4%
3yo
GENERALNOTES:
(a) Type 1 tests are to be conductedwith the specifiedgas. Deviations are based on the specified values
where pressures and temperatures are expressedin absolute values.
(b) The combined effect of inlet pressure, temperature and molecular weight shall not produce more
than an 8% deviation in the inlet gas density.
(c) The combined effect of the deviationsshall not exceed the limited of Table 3.2. Cooling temperature
difference is defined as inlet gas temperature minus inlet cooling water temperature.
TABLE 3.2
PERMISSIBLE DEVIATION FROM SPECIFIED OPERATING PARAMETERSFOR
TYPE 1 AND 2 TESTS
limit of Test Values as Percent of
Design Values
Parameter
Specific volume ratio
105
Flow coefficient
Symbol
Min
vhd
95
96
Milx
104
90
[Note (111
10
200
NOTE:
12
PTC
COMPRESSORS
ASME
AND EXHAUSTERS
10-1997
TABLE 3.3
LIMITS OF DEPARTURE FROM IDEAL CAS LAWS
GASES
OF SPECIFIED A N D TEST
~~
1.12
Pressure
Maximum
Ratio
Ratio
k rnaxlk min
1.4
2
8
16
32
1.10
1.o9
.O8
1.O7
1 .O6
~~
~~
Min
Max
Min
-0.344
-0.1 75
-0.073
-0.041
-0.031
-0.025
0.279
0.1 67
0.071
0.050
0.033
0.028
0.925
0.964
0.982
0.988
0.991
0.993
Max
1.O71
1 .O34
1 .o1 7
1 .o1 1
1 .O08
1.006
GENERALNOTES:
(a)Where:
X =
T av
- 1 and
v a-r
ap T
(b) Maximum and minimum values of k shall apply to both the specified andtest gas over the complete
range of conditions.
(c) When these limits are exceeded by either the specified gas or the test gas at any point along the
compression path real gas calculation methods shall be usedfor that gas. Ideal or real gas method may
be used if these limits are not exceeded.
Themethodsof
this Codemay be applied
for conversion of testresults to specified operating
condition results for compressors which maybe
treated as oneormoresections.
A section is that
portion of a compressorwherenointermediate
streamleavesorentersbetweenone
impeller inlet
andthe sameoranother
following impeller discharge. See Table 3.2. Heat exchangers are excluded
from the interior of the section boundaries. Section
boundaries are indicateddiagrammatically in Fig.
3.1. The gas state and flow rate shall be established
for each stream where it crosses the section boundary.Thepowerabsorbedandheatloss
or gain by
naturalambientheattransfermustalso
be determined.
3.3.2
in practice differ
between test and specified operating conditions due
to leakagedifferences.Forexample,
it is common
to test at reduced inlet pressure andthereduced
differential pressureacross
a seal to atmosphere
could result in zero or negative leakage. As a result,
volume flow ratio equalitycannot
be achieved
betweentestandspecifiedconditions.
Therefore, it shall be necessary to estimatethe
leakage ratio; that is, the leakage mass flow divided
bythe inlet mass flow for both test andspecified
conditions. If theleakage ratio differencebetween
test andspecified is significant,theseeffectsshall
be appliedto the calculations of capacity and power.
PTC
ASME
1 O- 1997
TABLE 3.4
PERMISSIBLE FLUCTUATIONSOF TEST READINGS'
Symbol
Measurement
~~~
Inlet pressure
Inlet temperature
Discharge pressure
Nozzle differential
pressure
Nozzle temperature
Speed
Torque
Electric motor input
Molecular weight
Cooling water inlet
temperature
Cooling water flow
rate
Line voltage
MW
T
2%
psia
"R
psia
0.5%
2%
PSi
"R
2%
0.5%
rPm
0.5%
Ibf ft
1 Yo
kW
lbdlbmole
0.25%
"R
gal/min
volts
2%
1Yo
2%
GENERALNOTES:
(a) A fluctuation is the percent difference between theminimum and maximum test reading divided by
the average of all readings.
(b) Permissible fluctuations apply to Type 1 and Type 2 tests.
NOTES
(1) Seepara.5.4.2.3.
(2) See para. 4.1 6 for further restrictions.
Power in
r--
-7
-
-------
/1
lest
section
boundary
Multiple
entry
streams
I
I ""-
"
I
I
I
"
"
"
Heat transfer
Multiple
exit
S T D - A S M E P T C 10-ENGL L997
0 7 5 9 b 7 0 Ob05YY3 T q b
ASMEPTC 10-1997
3.5
3.4
subscript
whelre:
go
h2-2)t
93
= 42
q3-2
Section 1 diiharge
from measurementa
side&earn
1 -
2 =
3 =
(r~-2)sp
95
- (rq7-2lt
(rql-2M.J
q
92
Min.
ql-2
Section
flange
computed
before
flange
,,o
105
Max.
FIG. 3.2
(rpl-5)t
9
44
= 45
Section
TYPICAL SIDELOAD
46
= 45
subscript7
Q8
rq7a = 41
QL%5
SECTIONAL COMPRESSORS
flange
6 -
inlet computed
110
105
110
Max.
Section 2 discharge
computed
from internal
measurements
before ridestream
2 mixed
go
95
9o
Min.
5 -
4 = Section
(r,&sp
(rq8-5)t
(r++p
(rq3-2)t
(rql.2)sp
(13
= -
subscript
re5=$
w-5
rq3-2
---------
h&t
3 mixed
105
110
Max.
from
inlet computed
95
Min.
8 = Section 3 discharge
flange measurements
= Section
&8-5)w
section
S T D - A S M E P T C LO-ENGL L997 m 0 7 5 9 b 7 0 0 b 0 5 4 4 5 8 1 9 m
COMPRESSORS AND EXHAUSTERS
ASMEPTC 10-1997
3.6
SAFETY
3.6.1 Thetest
gas usedshallbe
in compliance
with local regulationsandprudentpractice
with
regard to flammability and/or toxicity.
Testgasesused
in a closed loop shall be
continuously monitored for composition and avoidance of combustible mixtures. Air or other oxidizing
gases shallnotbeused
in aclosedloop.
3.6.2
3.7
PIPING
3.5.6 Extraction Sidestreams. Whentheintermediate flows are removed (.e., bleed-off) fromthe
compressor,they will cross asectionboundary.
17
EXHAUSTERS
AND
COMPRESSORS
O .3
0.2
0.1
BI
0.0
-0.1
-0.2
I
-0.3 I
O
I
I
0.2
I
I
I
I
0.4
I
I
I
I
0.6
0.8
FIG. 3.3
1.o
Mmsp
A L L O W A B L EM A C H I N EM A C HN U M B E R
CENTRIFUGAL COMPRESSORS
I
I
1.2
1.4
1.6
DEPARTURES,
3.8
INSTRUMENTATION
FIG. 3.4
3.10
A L L O W A B L EM A C H I N EM A C HN U M B E R
COMPRESSORS
(i) lubricant temperatures, inlet and outlet of bearings, seals, and speed changing gear, if applicable
(j) coolant and lubricant flows, if applicable
fk/ barometric pressure
(I) gas analysis, if atmospheric air is not the test gas
(m) time
PRETEST R U N
3.10.1 Thecompressorshallbeoperated
for sufficient time at the required conditions to demonstrate
acceptable mechanical operation and stable values
of all measurements to be takenduringthe
test,
Preliminarydatashallbetaken
to familiarize test
personnel, to determine if all instrumentsare functioning properly,and
to ascertain if thereading
fluctuations fall within the limits prescribed in Table 3.4.
3.10.4 Thepretestrunmaybeconsidered
as part
of the test if it meets all requirementsofthetest.
3.11
TEST OPERATION
DEPARTURES, AXIAL
~~
S T D - A S M E P T C 10-ENGL L797
m
EXHAUSTERS
FIG. 3.5
essentialinstruments.Three
sets of readingsshall
be taken during eachtest point.
3.1 1.3 The minimum duration of a test point, after
stabilization, shall be 15 minutesfromthe start of
the first set of readings to the end of the third set
of readings.
3.11.4 When a test is only to verify a single specified condition, the test shall consist of two test points
which bracket the specified capacity within a range
of 96 percent to 104 percent.
DEPARTURES,
PTC 10-1997
COMPRESSORS
ASME
AND EXHAUSTERS
12
10
9
FIG. 3.6
SCHULTZ COMPRESSIBILITYFACTOR
- FUNCTION
3.12
system re-
INCONSISTENCIES
3.12.1 Where
four
independent
instruments
are
used to measure a pressure ortemperaturevalue
andonerecordedobservation
is inconsistentdue
to measurementerror, its valueshall be discarded
and the value determined from the average of the
otherthree.Wherefewerthanfourindependent
measuring devices are used, all values shall be used
and averaged to determine the measurement value.
3.11.7
Thechoke
flow maybedeterminedby
gradually opening the discharge throttle valve while
maintaining speed and inlet pressure until the flow
remainsessentiallyconstant
with decreasing dischargepressure.
If the compressoris to be operated as an exhauster
ortested with an opendischarge,thechoke
flow
may be determined by gradually opening the
inlet
valve whileholding speed anddischargepressure
constant.
If choke flow is to bedetermined,the
facilities
shall be designed so as not to limit maximum flow.
Reduced Pressure, Pr
3.13
- FUNCTION X VERSUSREDUCEDPRESSURE
22
SECTION 4
- INSTRUMENTSANDMETHODSOF
MEASUREMENT
12
taps displaced 45 deg. fromthemandatleast
in. downstream.
In special cases when atmospheric conditions satisfytherequirements,thecompressormayberun
without an inletpipe as shown in Fig.4.2.The
inlet opening shall be protected with a screenand
bellmouth suitably designed to eliminate debris and
minimize entrancelosses(seepara.4.4).The
total
inlet pressure is equal to atmospheric pressure. Temperature measuring devices shall be located
on the
screen to measure the temperature of the air stream
at thecompressor inlet.
Forcompressors with an axial inlet, the impeller
may,undersomeconditions,
produce a vortexat
the pressure station to cause substantial error in the
measurement of inlet pressure.Users of this Code,
by agreement,mayusevanes suitably designed for
low pressure loss to prevent rotation at the pressure
taps.Thestaticpressurestationsshall
notbe less
than four pipe diameters upstream of the compressor
flange as shown in Fig.4.3.
4.1 METHODS
4.1.1 The choice of methods provided in this Code
will depend on the compressor,thespecified
gas,
andthetypeoftestselected.
4.2
INSTRUMENTATION
4.2.1 ThePerformanceTestCodeSupplements
in
the PTC 19 series on InstrumentsandApparatus
provide authoritative information concerning instrumentsand their use andshouldbeconsultedfor
suchinformation.Theselection
of instrumentation
shall be determined by the uncertainty limit requirementsofthetest
as well as suitability forthetest
siteconditions.Theinstrumentselectionshallbe
justified by calculation that the uncertainty in results
meets the statedtestobjectives.
Instrumentation is required to determine the inlet
and discharge gas states, flow rate, and compressor
speed. Depending upon the method selected, additional instrumentation may be required to determine
test power.
4.3
PIPING
6 minimum
12 in. minimum
12 in. minimum
6 in. minimum
8 in. minimum
h
temperature
Inlet
Inlet static
pressure
4 taps spaced
90 deg. apart
Inlet Opening
Preceded By
Straight run
Elbow
Reducer
Valve
Flow device
Discharge static
pressure
4 measuring
from
taps spaced
90 deg. apart
4 measuring stations
spaced 90 deg. apart
(45 deg. from static
pressure)
Discharge Opening
Followed By
Minimum Dimension
A
2D
30
20
30
60
IOD
30
ED
3D
Straight run
Elbow
Reducer
Valve
Flow device
5D
Discharge temperature
4 measuring
taps
spaced
90 deg. apart (45 deg.
pressure)
static
Minimum Dimension
20
20
30
30
80
30
30
50
5D
IOD
Inlet Configuration
FIG. 4.1
INLETANDDISCHARGECONFIGURATION
Protecting screen
Inlet pressure
by barometer
Inlet temperature
4 measuring stations
spaced 90 deg. apart
ASMEPTC
10-1997
minimum dimensions
1D min
Inlet temperature
4 measuring taps spaced
90 deg. apart (45 deg.
from static pressure)
A-
4Dmin.
L Straightener (optional)
See Fig. 4.8, para. 4.5
FIG. 4.3VORTEX
PRODUCING AXIALINLET
Discharge temperature
4 measuring locations
spaced 90 deg. apart
FIG. 4.4
OPEN DISCHARGE
4.3.4 Figures 4.6 and 4.7 show a typical arrangement for testing with ageneralclosed
loop and
closed loop with sidestreams.
4.4
PROTECTIVE
SCREENS
S T D - A S M E PTC L O - E N G L L797
ASME PTC 10-1997
EXHAUSTERS
AND
COMPRESSORS
minimum dimensions
Discharge temperature
4 measuring tapsspaced
90 deg. apart (45 deg.
from static pressure)
Straightener (optional)
See Fig. 4.8, para. 4.5
FIG. 4.5
DIFFUSINGVOLUTEDISCHARGEWITHNONSYMMETRICFLOW
arrangement
straightener
and
III
Drain pot
a
Drain tank
FIG. 4.6
TYPICALCLOSED
26
Gauge
LOOP
S T D . A S M E P T C LO-ENGL
L777
Nozzle
arrangement
(See PTC 19.5)
Y-
Drain pot
-i!&+
Drain
tank
Gauge
II-I Cooling
0757b70 O b 0 5 4 5 5 7bB
ASMEPTC 10-1997
Inlet
Inlet
measuring
stations
(See Fig. 4.1
'
""" jJ
Valves for
charging and
venting gas
Nozzle
measuring
arrangement
(See PTC 19.5)
Throttle
valve
Relief
valve
n
-
I
m
Flow equalizer
and straightener
-"
LJ
\ Flow nozzle
Drain
tank
Gauge
?
i
l
+
FIG. 4.7
TYPICALCLOSEDLOOP
27
WITH SIDESTREAM
S T D - A S M E P T C 10-ENGL 2 9 9 7
0759b70 O b 0 5 4 5 b bTq
COMPRESSORSANDEXHAUSTERS
4.5
4.6
PRESSURE
MEASUREMENTS
4.6.1 Reference
should
be
made
to PTC 19.2,
forgeneral information on instruments to measure
pressure.For
therangeofpressures
likely tobe
measured in compressortest, the manometerand
thedeadweightgage
shall beused as standards.
Pressure transducers and other pressure measurement
devices can be used. These can
be calibrated using
deadweight testers or manometers. Deadweight testers shall be certified by a competentlaboratory.
Where gage lines are filled with liquids, means shall
be provided to measure the liquid level, and a
correction shallbe applied forunbalanced liquid
head.
4.6.2 Bourdontubesorsimilar
gages shouldbe
selected to operate in themid-range of the scale.
Thediameters ofthe scalesand thearrangement
of thegraduations shall permit easy reading.The
temperature of the gage during calibration shall be
within 40F of theambienttemperature prevailing
during the test.
4.5.2
4.6.3 Manometerscanbeeither
U-tube orsingle
leg design. Smallboremanometers are subject to
appreciableerrorresultingfrom
capillary forces,
variable meniscus,andrestrictedseparationofentrained gas bubbles. These errors vary with the type
of fluid, the tube diameter, and the tube cleanliness.
Singlelegmanometers
shall be checkedforzero
position before and after test. Manometer fluid shall
be chemically stable when in contact with the test
gases andmetalpartsof
the instrument.
The specific gravity andthe coefficient of temperatureexpansion
of the fluid shall be determined
before the test. See PTC 19.2 for further guidance.
where
Ah= total area of holes in plate, sq in.
Ap= area of cross section or pipe, sq in.
q = inlet volume flow, cfm
pi= inlet density,Ibmper CU ft
Dp= diameter of pipe, in.
pp= densityof gas in pipe upstreamofplate,
Ibm per CU ft
Ap= pressure drop acrossplate, psi
4.6.4 Deadweight gages andtesters shall be selected to suit the pressure range. Deadweight gages
cannotmeasure rapid pressurechangesand where
necessarytheyshallbe
installed in parallel with a
Bourdon tube gage, transducer, or other instrument.
4.6.5 Transducersshallbeselected
with pressure
ranges appropriate for the expectedtestpressures.
They must be calibrated before and after each test.
When automated data collection equipment is
used with a pressure switching device, and a single
transducer, that transducer shall be selected to cover
theentirerangeofpressure.Whenusingpressure
switching devices, sufficienttime between successive
switch points shall be allowed so that the transducer
pressure will reach equilibrium for the selected pres-
4.5.3 A combined flow equalizer and flow straightener is used with flow nozzleswhererequiredby
PTC 19.5. See Fig. 4.8, sketch (dl. The flow straightener shall be the multitube type as shown in Fig.
4.8, sketch (b), preceded by a flow equalizer onehalf pipediameter upstream. Alternatively three flow
equalizers spaced one pipe diameter apart may be
used as shown in Fig. 4.8, sketch (e).
28
L14
(bl Multi-Tube Flow Straightener
min.
IfD-"-J
(e) Multi-Tube Type Equalizer and Straightener
FIG. 4.8
STRAIGHTENERS AND
EQUALIZERS
29
4.6.6
4.6.7
4.7
TEMPERATURE
MEASUREMENTS
Resistancetemperaturedetectorsorthermistorsshouldbeselected
for the appropriate range.
Cautionshould be takenbecausesome
of these
deviceshave a relatively slow responsetime.
4.7.5
Nozzle temperature
m
I
Nozzle pressure
2 taps spaced 90 dag.
"
4
"
para.
Protecting screen
See
4.4
FIG. 4.9
SD minimum
uncertainty shouldbereflected
in theuncertaintybarometer
pressure. Differential pressure i s measured
analysisand in the value of the uncertainty assigned
from two statictaps located 112Ddownstream of the
to thesestations.The internal temperaturemeasurenozzle flange.Temperature is measuredbysensors
ment i s always required when sidestreamand mainat the screen.
stream flows mix internally.
4.8
CAPACITY
MEASUREMENTS
4.8.1 Flowmaybemeasuredbyusingan
ASME
flow nozzle, concentric square edgeorifice, Herschel
type venturi tube, or alternative devices of equal or
betteraccuracy.Reference
shall be made to PTC
19.5, Flow Measurement, for general instruction and
detailed description of the various primary elements
and their applications. Other references are provided
in Appendix D. The interested parties shall mutually
agree upon the type of metering device to be used
and the choice shall bestated in the test report.
4.8.2 The flow measuring device may be located on
either the inlet or discharge side of the compressor.
It
shall be usedto determine the net capacity delivered,
or in the case of anexhauster,the
netcapacity
exhausted, which excludeslosses by shaftleakage,
balancing pistons,condensation,andothernormal
leakage that maybeinherent
in the compressor
design. Multiple devicesare required for multiple
inlet or discharge flow sections.
The nozzle may be used with an open inlet.
The nozzle arrangement shown in Fig. 4.9 may be
used for the test of compressors as exhausters. The
4.8.3
4.8.5 Formulasfor
calculating mass flow for a
variety offlow measuringdevices as provided in
PTC 19.5 shallbeused.Methodsare
included for
the determination of the discharge coefficient, fluid
expansionfactor,andmeteringelementthermal
expansion coefficient for various flow elements.
4.8.4 The nozzle maybeused with an open discharge: Figs. 4.1 O and 4.1 1 show optional arrangements of the flow nozzle on the outlet end of a
pipe for use where it is convenient to discharge the
gas to atmosphere. For a subcritical flow, the nozzle
differential pressure, A,, will be less than the barometric pressure and it shall be measured fromimpact
tubes, as shown in Fig. 4.10. Where theavailable
gaspressures permit,the nozzle maybesizedfor
operation at critical flow. In this case the differential
pressure will begreaterthanbarometric
pressure,
and it shall be measured fromstatic taps located
1 D upstream of the nozzle as indicated in Fig. 4.1 l .
In both cases the minimum length of straight pipe
precedingthe nozzle shall be 1OD andthe pipe
diameter shall be a minimumof 1.66 timesthe
nozzle throat diameter. Temperature measuring Stations shall be located6 0 upstream. The flow straightenerand/or flow equalizer, as described in para.
4.5, shall be used.Users of these arrangements are
cautioned to observe the distinction between critical
and subcritical flow. It shouldbenotedthatthe
velocity ofapproach is included in measurements
made with impact tubes.
S T D - A S N E P T C 1 0 - E N G L L797
0 7 5 7 b 7 0 ObOSLIbL T b l
COMPRESSORS
ASME
AND EXHAUSTERS
PTC 10-1997
pressure NozzletemperatureNozzle
2 measuring stations
spaced 90 deg.
II
Il
II
II
I I
Il
.-
20
====
1..
.
r1
"
"
"
"
I.
"
"
"
"
Efzz
"
"
0.23d
1OD minimum
SPECIAL NOTE: d not more than 0.6D for any nozzle arrangement
FIG. 4.10
DISCHARGENOZZLE ON ANOPENLOOP,SUBCRITICALFLOW
Nozzle temperature
2 measuring stations
spaced
deg-
7
1ODminimum
, - / , - -
FIG. 4.11
4.9
Nozzle pressure
2 measuring stations
spaced 90 deg.
""-I
DISCHARGENOZZLE ON ANOPENLOOP,
CRITICALFLOW
CAS COMPOSITION
4.9.1 Thetest gas must be defined.At the minimum, sampling will betaken at thestartand end
of each test.
4.9.2 Special precautions shall be taken when testingwith theclosed loopto eliminate all liquids
fromthe
gas
stream
andstaticinstrumentlines.
When dealing with gas mixtures subject to variation,
samplesshallbetaken
at eachtest point andbe
analyzed by spectrographic,chromatographic,or
chemical methods. The sample shall be taken from
the piping such that there i s no condensation before
torque meters shall be of a type suitable for calibration. The torsion member shall be selected for readability and accuracy at thespeed and load prevailing
during test.
4.14SHAFTPOWER
BYELECTRICAL
MEASUREMENTS
4.10.2 Thespeed of a compressor driven by synchronous motors may be determined from the number of poles in themotorandthefrequency
of
the power systems. If gearsareused
betweenthe
measuring point and the compressor shaft, the speed
ratio shall be computed from a count of the number
of teeth.
4.12
4.14.2 Efficiency
determination
by input-output
measurements may not be practical for large motors.
For large motors the loss method may be used. The
segregated losses of an induction motor shall include
friction and windage, core loss, 12R loss of the rotor
and the stator, and a load loss. These measurements
shallbemade
in accordance with current ANSI
standards.
METHODS OF SHAFTPOWER
MEASUREMENT
4.1 2.2
sible applicationsfor
4.13
4.15SHAFTPOWER
BY HEATBALANCE
MEASUREMENTS
SHAFTPOWER BY TORQUE
MEASUREMENTS
S T D - A S M E P T C LO-ENGL
L977
m
ASMEPTC 10-1 997
k-
Location C
Location A
Location B
FIG. 4.12
TYPICALSIDESTREAMINLETAREA
surement of the temperature rise (suchas with differential thermocouples). Evidence of nonuniform temperature distribution morethan 2 percentofthe
temperature rise at either the inlet or the discharge
measurement station, may require one of the following procedures
be
used
at the
offending
measurement
station:
(7) Apply insulation to the piping upstreamof
the temperature measurement station in an effort to
minimize thermal gradient.If successful, the temperature measurement installation need not be changed.
(2) Move the temperature measurement station
away from the compressor and add pipe insulation.
This might beparticularly effective when temperature
4.15.2
''
with
4.15.3 Theheatbalancemethod
shall be used
with the following precautionsand limitations.
(a) The inlet and discharge temperatures shall be
measured with instruments suitably selected and applied toprovide combined accuracywithin 1 percent
of the temperature rise. When the rise
is less than
SOOF, consideration should be given to direct mea35
'STD-ASME
P T C 10-ENGL L 7 7 7 W 0757b70Ob054b4
770
COMPRESSORS AND EXHAUSTERS
4.16
SHAFTPOWER
BY HEATEXCHANGER
METHODS
4.16.1 When it is not possible
or
practical to
measureshaft power directly or by a compressor
heat balance, and a heat exchanger is incorporated
in the test arrangement, the heat transferred to the
cooling watermaybeused
to determinethenet
compressorshaft power.
4.16.2 Methods to accountforthemechanical
losses are discussed in para. 4.18. External heat loss
from the casing, piping, and coolermay be evaluated
in accordance with para. 4.17.
4.17
HEAT LOSS
4.16.3 Theheatexchanger
methodshall beused
with the following precautionsand limitations.
(a) The cooling fluidsupply shall be stablein pressure and temperature so that the fluctuation of flow
rates will not deviate more than 2 percent and the
36
4.18MECHANICAL
LOSSES
4.18.1 Whenpractical,theheatequivalent
of the
mechanical losses (integral gears, bearings, and seals)
shall bedeterminedfromthetemperatureriseof
the cooling fluid. The quantity of fluid flowing shall
be determined by calibrated flow meters.Theheat
equivalentoftheexternal
losses as well as the
frictional loss in the mechanical seals, if used, shall
be determined and included in the total mechanical
losses.
Where the mechanical losses are well known and
documented, the calculated values or those values
determinedfrom
prior testing
may
be
used by
agreement by testparties.
4.1 9
INSTRUMENT CALIBRATION
4.19.5 Torque
meters
shall be calibrated by
applying torque with certified standard weights, load
cells, or other appropriate devices spaced to cover
the working range. For strain gage types, the calibration shall includethe brushes, leadwires,andthe
indicatinginstrument.
37
SECTION 5 - COMPUTATIONOFRESULTS
percent
one
in efficiency
and
two percent in dischargespecificvolume.'The
ideal gas equation of
The Calculation
The process Of
state,
14dPv
= RT, and the correspondingderived
establishing
compressor
performance
from
test
data
equations in Tables
and
may be
involves a number of calculation steps. This Section
Forgases
with variablespecific
heats,average
is presented in the following chronologicalorder:
properties are calculated at the arithmetic mean
FORMAT
5.1
0
0
0
0
0
0
0
Perfect
real
or
gas treatment is selected.
section
temperature.
The appropriate test speed is calculated if a Type
2 test is to be performed.
The raw test data is processed.
Test performance is calculated.
Test performance is expressed in dimensionless
form.
Reynolds number correction is applied.
The corrected dimensionless expressions are used
to predict performance at specified operating conditions.
Theimportantsubjectofuncertainty
is treated
separately. The section format
is intended to guide
the user in basic calculation procedure and
to present
the necessary equations. Appendix E is provided as
a background theory source and further explanation
of the equations.
5.2COMPUTATIONALMETHODS
AND REALGASES
5.2.1.3
Tabulated
Properties
and
Equation
of
StateMethods.
Puregases
and gas mixtures for
which tabulated data properties exist may be treated
as real gases.
There are many gas property correlation equations
of state for purecomponentsand
gas mixtures.
Many of the generalized equations of state provide
sufficiently accurate predictions of gas properties to
be used in conjunction with the calculation methods.
Theuse ofeitherofthesemethods
will require
iterative procedures to satisfy the equations in Tables
5.2 to 5.4.
FORIDEAL
5.2.1Choice
of Methods. Thetest andspecified
gases maybetreated as eitherperfectorrealdependingupon their respectivethermodynamicbehavior. For thepurposesofthisCode
ideal gases
are those which fall within the limits of Table 3.3.
Gases which exhibit deviations beyond these limits
are considered real. Threedistinct treatments of gases
are recognized in the computational procedures. The
appropriatechoice will dependupontheselected
gas, knowledge of itsproperties,andthedesired
accuracy.
39
S T D * A S M E P T C LO-ENGL L777
0757b70 Ob054b7 4 8 T
TABLE 5.1
IDEAL GAS DIMENSIONLESS
Parameter
PARAMETERS
Eq. No.
Assumption
[Pinlsp
[5.1T-21
Flow coefficient
-RT;.
-31
15.1
[PinIr
k-1
I5.1 T-41
where nt =
[1-
(5.1T-51
In
PiTd
Isentropic efficiency
[%If =
Polytropic efficiency
hplf
Total
work
input coefficient
k.1
and,
[rlplr =
40
ASMEPTC 10-1997
TABLE 5.2
REAL GASDIMENSIONLESS PARAMETERS
Parameter
Mathematical
Operating
Description
Test
Conditions
at
!
Eq. No.
Flow coefficient
[5.2T-1j
Assumption
4sp
4f
where [ndr =
[T]
15.2T-41
In -
(5.2T-51
and f, =
Irplr =
where nt =
[T]
-
[5.2T-7j
In
[5.2T-8]
Isentropic efficiency
[ d r =
hd
[vSlsp= I
d t
Remccllr
- hi
[Table continued on next pagel
41
PTC
ASME
1O- 1997
Parameter
Mathematical
Operating
Test
Description
Conditions
at
Assumption
Polytropic efficiency
Total
work
input
coefficient
GENERAL NOTE: Appropriate units mustbe chosen to render the parameters dimensionless. Further explanation
of the equationsis available
in Appendix E.
.e.,
(E)
5.3.2 Test Speed Selection. The volume ratio limitation of Table 3.2 may be met by controlling the test
speed. The appropriate test speed i s calculated from
where
[5.3.61
n
l():
and,
Foran
ideal gas,
n
-T
"
n- 1
k
p
i5.3.71
k"1t
= 1r"lsp
n =
[5.3.4]
42
1
Y-m(l
+x)
I5.3.81
S T D - A S M E P T C LO-ENGL
L777
0757b70 Ob05Y70
T7Y
COMPRESSORS A N D EXHAUSTERS
where
m =
?(L+
x)
CP TP
[5-3.91
AF=
1 O0 (AH - AL)
15.4.1 1
where
TestSpeed Validation.When the actualtest
conditions differ from the estimated values, the most
appropriate test speedwill depart fromthe previously
calculated testspeed.Some deviation is allowable.
Thetestspeed
is acceptablewhenthedeviation
satisfies the limits of Table 3.2.
5.3.3
5.4
AF=
fluctuation expressed in
YO (Table 3.4)
ProcessingRawData.Acceptableraw
5.4.2.5 TotalConditions.
Gas
state
static
test
point data shallbeconverted
to total condition
values for thecomputationalprocedure.This
does
notpreclude final presentation in terms of static
conditions, but total values are usedin the intermediate computations.
The relationship between static and
total properties
is velocitydependent.Average
total properties are
estimatedhereinfromtheaveragevelocity
at the
measurementstation.
Theaverage velocity atthemeasurementstation
is givenby
data
M =
5.4.3
5.4.4.2Simplified
Method. Thedifferencebetween total andstatictemperaturesmaybeevaluated from
Test
Pressure
5.4.3.1 Simplified Method. Formeasurement station Fluid Mach numbers of 0.2 orless theeffects
of compressibility are small, A good approximation
of velocity pressuremaybe obtained by assuming
incompressible flow at the measurement station and
calculating anapproximatedensityfromthemeasuredstaticpressureandmeasuredtemperature.
Thus
Va, = ~ / 6 0 p A
[5.4.61
15.4.31
5.4.4.3Refined
Method. For
cases
wherethe
measurementstationFluid
Mach number exceeds
0.2 for areal gas, thediscussion in Appendix G
gives guidelines for more accurate methodsfor relating total temperature to measured temperature. For
cases involving extremevariationfrom
ideal gas
behavior,such as near the critical point, the total
temperature may differ greatly from the value indicated by para. 5.4.4.1 and the methods outlined in
Appendix G shouldbe used.
5.4.3.2Refined
Method. For
cases
wherethe
measurementstation Fluid Mach numberexceeds
0.2, or when abetteraveragevelocityestimate
is
desirable,therefinedmethodofAppendix
G may
be used. This method is based upon the assumption
of uniform compressible flow at themeasurement
station.
5.4.4
Test
Temperature
5.4.4.1RecoveryFactor.
Thetemperature indicated by asensing element is normallyavalue
somewhere between thestatic and total temperature,
depending upon the ability of the sensor to recover
the converted kinetic energy of the gas stream. This
ability is defined in terms of a recoveryfactor,
5.4.5TestDensityandSpecificVolume.
Thetest
total density is calculated from the test total pressure
and total temperature as
[5.4.81
[5.4.9]
TC
COMPRESSORS
ASME
AND EXHAUSTERS
v, = -
10-1997
[5.4.10]
Pt
where
[Coutwh - Lnwhl
5.4.6.2Volume
Flow Rate. ThisCodeuses
a
flow rate definition in the calculation process which
has the units of volume flow rate. It is
q=-
Psh,
[5.4.14]
Pg, iPparasitic,
[5.4.11]
where
mass flow rate
P= total density
This definition is consistent with the use of total
properties in the calculation procedure. It does not
represent the actual local volume flow rate because
it is based upon total rather than static density. All
references to calculated volume flow rate imply this
definition unless otherwise stated.
W=
Pg,
Qextl
[5.4.15]
33000
W,=
CP',
Qr=
5.4.7.4Casing
Heat Transfer. Theexternalheat
loss or gain fromthesectionmaybecomputed
from measurements of the exposed surface area, the
average temperature of the surface, and the ambient
temperature from
5.4.7.1ShaftPowerMethods.
Whenpower input i s measured by instrumentssuch as a torque
meter, dynamometer, or calibrated motor, the shaft
power is calculated using the appropriate formula.
Gas power is calculated by subtracting the parasitic
losses fromthe shaft power (seepara. 5.4.7.5 for
parasiticlosses).
Qr
= ISA tc -.
ta)
hl60
15.4.1 61
where
Sc= heat transfer surface area of exposed compressorandadjoiningpipeforsection
of
interest
tc= casing surface temperature
ta= ambient temperature
hr= coefficient ofheattransfer for area (combinedconvectionand radiation)
Where
the
casing
surface
temperature
varies
be
widely,theaccuracy
of thiscalculationmay
improved by treating small areas of the surface
= measured value
- tl) + Qr
where
5.4.7TestPower.
The calculationof test power
dependsupon the methodof measurement.Both
shaft power and gas power may be of interest. Shaft
power is the power input to thecompressordrive
shaft.Gas power is the power delivered to the gas
in thesection(s1ofinterest.
Psh,
= [wwcpW(r2
S T D - A S M E P T C LO-ENGL L777
Pg = Psh
- Ppararitic
[5.4.19]
33000
(
h
7
Qr)
yielding
Pparasitic
Pmech
Pother
15.4.1 71
(a) Mechanical Losses. Mechanical lossesare always considered to be parasitic losses. Those losses
due to lubricated gears, bearings, seals, etc., may be
estimated from the lubricating oil temperaturerise.
Othermechanical losses from seals, bearing,etc.,
which do not contribute to the lubricating oil temperature rise shall be determined separately. Thatportion
of the mechanical loss evident in the lubricating oil
temperature rise is given by:
5.5
Prnech
= [wcpAtl 33000
[5.4.18]
5.5.1
where
mass flow rate of the lubricating or sealing fluid
cp= specific heat of the lubricating orsealing
fluid
At= temperature rise of the lubricating or sealing fluid
(b) Other Parasitic Losses. When the shaft power
method is used, power supplied to drive auxiliary
equipment is treated as parasitic. Also, power supplied
to sectionsof a multisection compressor other than
the section being tested i s considered parasitic.
When the heat balance method is used, and total
shaft power i s defined toinclude power to drive
auxiliary equipment, theauxiliary power requirement
is treated as parasitic.
W=
M m = U/ai
i5.5.1 I
di
Yi
=G
(5.5.31
Rem = Ublu
[5.5.41
46
DIMENSIONLESS PARAMETERS
ASMEPTC 10-1997
5.6.1TheSingleSectionCompressor
5.6.1.1 Description. Thesinglesectioncompressor from inlet to outlet measurement stations experiences no gas cooling otherthannaturalradiation
and convection. No gas flow is added or removed
other than that lost through
seal or balance piston
leakage. No condensationoccurs.
5.5.3Specific
Volume Ratio. Thespecificvolume
ratio is the ratio of inlet to discharge total specific
volume.
5.5.4VolumeFlowRatio.
Thevolume flow ratio
betweenany two points x and y in thesection is
given by
5.6
CALCULATIONS FORSPECIFIED
OPERATING CONDITIONS
is calculated
di-
mensionlessparametersarecalculated
for thetest
conditions and extended to specified conditions.
47
10-1997
TABLE 5.3
TOTAL WORK INPUT COEFFICIENT, ALL GASES
Parameter
Mathematical
Operating
Conditions
Test
Description
at
Total
work
input
coefficient
(Note
(heat balance method)
Eq. No.
llspwmption
(1 )I
(Pd
coefficient
input
work
Total
(heat balance method)
[arlt=
- fpdtie)33000
It
[5.3T-21
1nSh1rp = Iadlt
GENERAL NOTE Appropriate units mustbe chosen to render the parameters dimensionless. Further explanation
of the equationsis available
in Appendix E.
NOTE
(1 1 This equation applies to a particular model as presented in Appendix E, para. E.3.12. Some of the terms may not apply in a particular
case. Additional terms may apply. The analysis of para. E.3.12 may be followed to develop appropriate equations.
equationsforthosequantitiesofinterest.
Typical
equationsused to do this areshown in Table 5.4.
For example, to calculate the discharge pressure
at the specified condition flow the following steps
are taken: (1 ) the pressure ratio is calculated from the
now known specified operating condition polytropic
efficiency and polytropic work coefficients, and (2)
the discharge pressure is the product of this pressure
ratio and the specifiedoperating conditioninlet
pressure.
5.6.2The
Multisection Compressor
5.6.2.1Description.
A multisection compressor
is a compressor which may be treated as a number
ofindividual singlesectioncompressorsoperating
in series.The output from each single section provides input to the next section. The section boundaries may be drawn to exclude intermediate components such as externalheatexchangers.
The following conditions shall bemet to treat a
compressor as a multisection compressor.
5.6.2.2
Calculation
Method for Multisection
Compressors. The specified operating condition per48
S T D - A S M E PTC LO-EFJLL
L197
m
ASME PTC 10-1 997
I m p
Test F
FIG. 5.1
ASME
10-1 997
TYPICAL CONVERSION
TABLE 5.4
OF DIMENSIONLESS PARAMETERS
Parameter
Mathematical
Description
Test
Operating
at Conditions
Eq. No.
15.4T-1]
[5.4T-21
(5.4T-31
Capacity
[5.4T-4]
I5.4T-51
[5.4T-6]
n-1
where
n
k
= -( t l d s p
(n - 1
(k - 1)sp
I5.4T-71
[5.4T-81
[5.4T-9]
[5.4T-1 O]
I5.4T-11I
In
[5.4T-12]
[Table continued on next page]
50
ASMEPTC
Parameter
Mathematical
Description
Test
Operating
at Conditions
Pressure ratio (real gas)
(Cont'd)
m5P
[-ZR (-1
Cf
[5.4T-141
x)]sf
Discharge pressure
[5.4T-15]
Pressure rise
I5.4T-161
Discharge enthalpy
15.4T-171
15.4T-191
[5.4T-201
Shaft power
Assumption,
n = &h
or,
Sf
ash
Sf
51
10-1997
~~
S T D ASME P T C 1 0 - E N G L L777
5.6.3
RA = 0.066
Rem
Rem
RB =
RC
[5.6.2]
13.67\
15.6.31
log ( + -)
Rem
RC =
0.988
~~~0.243
i5.6.41
where
of the flow
passage, in.
The polytropicwork coefficient should be correctedfor Machine Reynoldsnumber in the same
ratio as the efficiency.
(4.8 x lo6 x b)
+ 0.934
0.000125 + -
5.6.3.1
0757b70 Ob05477 1 T 1
5.6.3.3 limits of Application. Since the performance variations increase substantiallyas the Machine
Reynoldsnumberdecreases,tests
of compressors
designed for operation at low MachineReynolds
numbers should be tested
at conditions closeto those
specified.Therefore,themaximumand
minimum
permissibleratiosbetween
Reml and Rem,
are
shown in Fig. 3.4. Also,
see
Appendix F and
Table E.2.
5.6.4 Mechanical losses. Whenthemechanical
losses at specifiedoperatingconditions
are not
known they may be determined from the following
equation:
in this
ASME
10-1997
'
The UADD
and URSS
models are the mathematical models which
are usedto combine bias and precision errors
to a single uncertainty
value. UAODprovides approximately 99 percent coverage while
URSS
provides approximately 95 percent coverage when neither
bias errors nor
precisionerrors are negligiblecompared to the other.
If the bias error
is negligible, both UADD
and U R S S provide95 percent
coverage.
54
ASMEPTC 10-1997
SECTION 6
6.1
- REPORT OF
CONTENTS
6.2TYPICAL
6.2.3Specified
Operating Conditions
(a) Gas composition and source for properties
(b) Inlet gasstate
(7) Total and static pressure'
(2) Total and static temperature'
(3) Total and static density',
(4) Relative humidity if applicable'
(c) Gas flow rate
(7) Inlet and discharge mass flow rate
(2) Inlet and discharge volume flow rate
(3) Capacity
(dl Discharge static and/or total pressure
(e) Coolant type, properties,flow rate, and temperature for cooled compressors
(0 Speed
(g) Others as needed
REPORT INFORMATION
6.2.4
Expected
Performance
Specified
at
eratingConditions
(a) Developed head
(b) Efficiency
(c) Power requirement
(d) Discharge total temperature
(e) Others as needed
6.2.1
General
Information
(a) Date of test
(b) Location of test
(c) Manufacturer
(d) Manufacturer'sserialnumbersandcomplete
identification
(e) Party or parties conducting test
ff) Representatives of interested parties
(g) Detailed written statement of the test
(h) Agreement made by parties to the test
Op-
6.2.2Description
of Test Installation
(al Type of compressor;radial flow, axial flow, etc.
( 1 ) Type of impellers; open, shrouded, cast, fabricated, etc.
(2) Number of stages
(3) Arrangement of casing and piping
(4)Pipe sizes; inlet and discharge
(51 Arrangement of intercoolers, if used
(6) Impeller diameter and blade tip widths
(b) Description of lubricating system and lubricant
properties
(c) Type of shaft seals
(d) Type and arrangements of driver; turbinedirect
connected,
motor
direct connected,
motor
and
gear, etc.
6.2.6SetupofInstrumentsandMethodsof
Mea-
suring
(a) Description of all allowed departures from this
Code which have been authorized by agreement
(b) Piping arrangementwith sketches and diagrams
(c) Location ofall measuring stationswith diagrams
and sketches
'Pressures, temperatures, and densitiesshould be clearly identified
as static or total conditions.
55
TEST
ASME
10-1997
(m
(e) Instruments usedfor the measurement of pressure, temperature, speed, composition ofgas, density,
and power
(0 Procedures andfacilities used for the calibration
of instruments
(g) Calibration data
(h) Instrument accuracy
(i) Source of test gas property data
0) Method of determining powerlosses, if any, between the power measurement station and the compressor input shaft
(k) Description of sampling and analysis method
for test gas
u)
*Iterativesolutionmay be required.
56
ASMEPTC 10-1997
(2) Enthalpy
(3) Average mixed gas state
(4) Energy loss or gain
(S) Rotor mass flow rate
(t) Mechanical loss
(u) Heat transfer loss
(v) Gas power
(W) Shaft power
(x) Head
6.2.9Computed
Test PerformanceParameters
(a) Isentropic total discharge conditions
(7) Temperature
(2) Density
(3) Enthalpy
(6) Polytropic work coefficient
(7) Overall isentropic volume exponent
(2) Polytropic work factor
(3) Polytropic exponent
(4) Polytropic work
(5) Impeller blade tip velocity
(6) Polytropic work coefficient
(c) Isentropic work coefficient
(7) Isentropic exponent
(2) Isentropic work
(3) Isentropic work coefficient
(d) Polytropic efficiency
(e) Isentropic efficiency
(0 Work input coefficient
(g) Total work input coefficient
( 1 ) Energy lost or gained via leakage
(2) Energy lost or gained via secondary flows
(3) Energy lost via casing heat transfer
(4) Mechanical loss
(h) Flow coefficient
(i) Volume ratio
ci., Machine Mach number
(k) Specific heat ratio, inlet and discharge
57
Analysis
APPENDIX A
USE OF TOTAL PRESSURE A N D TOTAL
TEMPERATURE TO DEFINECOM.PRESSOR
PERFORMANCE
(This Appendixisnot
[A-41
A.2
The difference between static and stagnation conditions is shown graphically on a Mollier Diagram,
Fig. A.l.
A.4 As will be notedfrom Fig. A.1, theprocess
ofcompressiontakesplacebetweenstates
(Y and
y. Some calculations regarding theinternal compression process might requirethe use ofstatic states
intermediate to (Y and y. However, as shown by
Eqs. [A-1I through [A-51, use of the stagnation properties for the external energy balanceof the compressor
is an excellent approximation because:
(a) "Charging" the compressor with receipt of gas
at the stagnation enthalpy hi (at stagnation pressure
pi) i s equivalent to charging it withreceipt of gas at the
static enthalpy h,a (at static pressure pa) plus kinetic
energy
Thisresultinvolvesstaticenthalpiesdetermined
by staticpressures andtemperatures.
A.3 When the stagnation concept is employed, Eq.
[A-21 becomes
59
FIG. A.l
and,
A.6
v;
2gcJ
A.5
The precedinganalysiscanbe
applied only
becausethesystemboundarieswerecarefullydefined so as to preclude any consideration of events,
Theotheruse
of thestagnationpressureand
stagnation temperaturein this Code is for the determination
capacity.
Capacity
of
is a volumetric flow
related
rate
inlet to
conditions.
Capacity
is defined
herein as thedelivered mass flow rate divided by
inlet total density correspondingto total pressure and
temperature. This is convenient because it permits a
clear definition of volume flow rate consistent with
mass flowwithout referring to thedesign of the
compressor.
60
APPENDIX B
PROPERTIES OF GAS MIXTURES
(This Appendix is not a part of ASME PTC 10-1997.)
B.4 With properties of the individual gases determined, the equivalent value of the property
for the
gas mixture maybecalculated
by summingthe
individual propertyvalues
on a total basis, .e.,
quantity of the gas times property value. The equations are summarizedbelow.
Enthalpy:
B.2 Whenthethermodynamic
state is such that
the gas mixtureand its constituent gases must be
treated as real gases, themethodofdefiningthe
thermodynamic state oftheconstituent
gases and
thus arriving at their properties shall be agreed upon
in writing prior to thetest.
Once the state of the gas i s defined, presumably
by pressureandtemperature,theotherproperties
of interestmaybeobtainedfromcharts,tables,or
equations of state.
Entropy:
e + n,Hj
18-31
62
APPENDIX C
SAMPLE CALCULATIONS
(This Appendix is not a part of ASME PTC 10-1 997.)
The sample calculations contained in this Appendix demonstrate the basic calculation principles
of this Code. Each sample highlights one or morefacets of the necessary procedures for application
oftheCode to realmachines.Thedatapresented
is typical and does notrepresentanyactual
operating unit. Additionally this data should not be taken as expected for any actual conducted test.
Sample C.l
Sample
Sample
Sample
Sample
Sample
Sample
Sample
63
Simpliving AssomptionsforThisSample
(al The gas (air) may be treated as an ideal gas with a constant specific heat (evaluated at the
average of the inlet and discharge temperatures).
(b) The Reynolds number correction i s negligible.
W = 600 Ibm/min
(2)inlet and discharge volume flow rates have to be calculated
(HRi)sp
= (0.6220
-(
")
0.6220 Ibm
(0.7825 psia)
Ibm da
(1 4.00
- 0.7826) psia
lbmole W
(28.97 ibm da)
180.02 ibrn W
lbrnole da
= 0.05921
lbmole W
lbrnole da
mole da ( M W d + mole
mole da + mole
(MW,)
Ibm da
) + 0.05921 lbmole
1 .O00lbrnole da
66
+ 0.05921 lbmole
S T D - A S M E P T C LO-ENGL L777
m 0757b70
Ob05492 h35
Ibm
= 28.36 lbmole
(d) Air specific heat at constant pressure is found using dry air and steam properties (Fig. C.l)
(cp)sp
mass da
(Cpda)
mass da
= 0.247
+ 0.03683
Btu
Ibrn da "R
'
(D.20)]
Btu
Ibrn "R
0.247
(0.247
m)Btu
(1.986
0.252
-)
Btu Btu
(0.252 Ibrn "R
)(
Btu
lbmole
Btu
[note: ( H R i ) s p =
Ibrn "R
Btu
= 0.252
(kd)sp
Ibrn
Ibrn "R
(kAp
Wpdsp
+ mass W (cpw)
+ mass W
"R
Btu
Ibrn "R
- (1.986
lbmole "R
lbmole
)=
28.36 Ibrn
)(
lbmole
)=
28.36 Ibrn
(0 Static specific volume at inlet and discharge flanges is found using the ideal gas law
67
(vstatic Jsp
S) $)
(144
fi3
= 15.13 Ibm
S)
(144
$-)
= 7.086 I bm
ft3
fg)
Average fluid velocity at inlet and discharge flanges (see para. 5.4.3.1)
( 6 0 0 2 ) (15.13
(Vilsp
ft3
(600
(vd)sp
= 97.40
g) L)
(7.086
= 129.4
ft
sec
ft
sec
(h) Fluid Mach number at inlet and discharge flanges (see para. 5.4.2.5)
97.40
(M;)sp
ft
sec
lb
68
) (28.36
lbmole ) (560 "R)
Ibm
= 0.0832
129.4
(Md)sp
Ibm2)(1545
Ibf sec
1.385 32.1 74
ft
sec
ft lb
lbrnole "R
) (28.36
lbmole ) (842.8 "R)
Ibm
= 0.0905
(i) Totaltemperaturesatinletanddischargeflangesarefoundusingtheenergyequationand
assuming an adiabatic process (see Eq. [5.4.6])
(97.40
(TJsp
560.0 "R
k)
Ibf
sec2 (0.247
= 842.8 "R
E)(32.1
2 778.1 7 Btu
Btu
k)
(1 29.4
(TdlSp
= 560.8 "R
ft Ibrn
2 7 7 8 . 1 7Btu
E ) (32.174-)
= 844.1 "R
ft Ibrn
747)
(0.252
Ibf sec
Ibm "R
ci,, Since the Fluid Mach number isless than 0.2, the total pressure may be calculated according
to the simplified Eq. i5.4.41
(97.40
(pi)sp
= 14.00 psia
J:,'
2 1 5.1 3 -
(32.174-
(1 29.4
(pd)sp
= 45.00 PSia
i-
2 (7.086
S)
ft Ibrn
in2
ft2)
Ibf sec2) 144 -
= 14.07 psia
A)
-)
2)
g)
= 45.26 psia
(k) Total density at the inlet and discharge flanges is found using the ideal gas law
1
(40.7):
(pi)sp
= ($)*P
Ibm
( 1 4 4in2
~lbf)
) (lbmole)
1
ft Ibf
= 0.06632
(560.8 "R)
fi3
"
28.36Ibm
(45.26
(Pd)sp
( g ) s p =
(1 545
ft Ibf
lbmole "R
S) $)
(144
)(
lbmole)
1
28.36.
Ibm
I bm
= 0.141 7 (844.1 "R)
ft3
"
\ lo4 min)
j= 1
= 2.983
fi2
lo6-
j= 1
sec2
mln
70
0.3667 psia
Ibrn W
= 0.01661
14.1 O psia 0.3667 psia
Ibrn da
= 0.6220
Btu
Ibrn da
(540.0
O R
Btu
- 530.0) R t (0.01661 G
Ibm ") (1092.2 - 38.35) Ibrn
Btu
(1095.5 - 38.35) Ibrn
W
0.01881
=)
(L
lbmole 7 (28.97
Ibrn
lbmole
da
Ibrn18.02
da
W
0.03024
lbmole
lbmole da
(MWdJ + mole
[mole damole
+ mole
da
(MW,)
I,
) + 0.03024 lbmole
Ibm da
lbmole da
1.00 lbmole da
28.65
+ 0.03024 lbmole W
I bm
lbmole
(c) Air specific heat is found using dry air and steam properties (see Fig. C.1)
(c,), =
I,
= 0.244
Ibrn da "R
Btu
Ibrn " R
71
Btu
+ 0.01 881 Ibrn (0.460
da OR
1.O00Ibrn da + 0.01881 Inrn W
1 .O00Ibrn da (0.244
cc,,,, =
Btu
= 0.248 Ibrn "R
+ 0.248
Btu Btu
"
Ibrn "R
(Mt
(0.244
S)-
(Z),
Btu
0.244 Ibrn "R
(0.1 986
0.248
(kd)t
lbmole
Btu O R )
Btu
Ibrn
O R
Btu
(o'248
= 1.397
Btu
(G
1 r
lbmole )
)(
lbmole)
28.65
Ibrn
= 1.388
(e) Static specific volume at inlet and discharge flanges is found using the ideal gas law
(1 545
(vstatic d t
ft Ibf
lbmole O R
) (28.65
1 lbmole) 540.0 "R
Ibm
"
(14.10%)
(144$)
72
= 14.34
ft3
I brn
(1 545
(Vstatic
d)r
ft Ibf
lbmole "R
) (28.65
1 lbmole) 830.0 "R
Ibrn
"
S) $)
(47.00
CO
ft3
= 6.613
I bm
(1 44
?)hr
(1 4.34
4 (T16.87ft)
(3600
(38,000
(VJt =
TT
(38,000%)
hr (6.613
(vd)t
ft
777 (F
10.020ft)'
E)
Ibm
T)
ft
= 97.45sec
E)
Ibrn
(3600
y)
= 127.5 sec
[g) Fluid Mach numbers at inlet and discharge flanges (see para. 5.4.2.5)
ft
97.45 -
(Mi),
=
sec
-)
ft m
J I .397 32 174Ibf!ec2
(I545
ft
lb
) (L28.65-) Ibm
lbmole "R
ft
sec
= 0.0852
(540.0 "R)
127.5 -
(Md)(
=
-)
ft Ibrn
41.388 32.1 74 Ibf sec2 (1 545
lb
lbmole "R
(h) Totaltemperatureatinletanddischargeflanges
adiabatic process (see Eq. [5.4.6])
73
= 0.0902
)(
1 e
(830.0)
"R)
28.65 Ibrn
is foundusingtheenergyequationforan
~~
~~
S T D - A S M E P T C LO-ENGL L997
0759b70 Ob05499 T 9 T
(97.45
= 830.0 "R
(Td)t
Z)
0 . 2 4 4 Btu
~ ) ( 7 7 8 . 1 7 E ) ( 3 2 . 1 7 4 - ft Ibm
Btu
Ibf sec2)
ft ,bm
-) Btu
= 540.8 "R
) = 831.3 "R
E)
7
(32.174 Btu
Ibf secz
(i) Since the Fluid Mach number is less than 0.2, the total pressure may be calculated according
to the simplified method of Eq. (5.4.41
(97.45
( p h = 14.10 psia
&)
(Pd)(
(14.34
2 6.61 3
-)Ibm
ft3
= 14.1 psia
lbfsec~(144%)
A)'
(32.1 74
74
(32.174 ' I b m )
(1 27.5
= 47.00 psia +
&)
Ibf Ihm)
sec2 (144
= 47.27 psia
$)
fj) Total density at the inlet and discharge flanges is found using the ideal gas law
(Pu)sp =
(14.17
(&),=
S) 5)
) (28.65 lbmole)
(540.8"F)
Ibm
Ibf
(1545lbmole "R
(Pd)sp
(&),
S)
)(
ft Ibf
= 0.06997 ft3
"
(47.26
Ibm
(144
$)
I bm
(1 44
lbmole)
= 0.1518 ft3
(831.3"F)
"
(1545lbmole O R
28.65 Ibm
[x
5
[TE
5
Ur],=
[Tz
5
Dr], =
j= 1
j= 1
j=l
(I) The shaft power was measured by the shaft power method to be (Psdt= 1097 hp (shaft power
method)
Theshaftpowercanalsobedeterminedusing
Eq. 15.4.141:
Ihm) (
to bemechanical losses(represented
(38,000
0.2459
) (60 F)
(42.440mln hp
~
(59.5
%)
(0.462
~~
~~
Btu
) (31.O "R)
= 1065 hp + 20.1 hp
= 1085 hp (heat balance method)
75
(m) The gas power can be calculated from the heat balance method using ,Eq. [5.4.13]
shaft powerusingtheshaft
power method.
(38,000
9f =
T) -k)
(60
(0.06997
*)
ft3
= 9051 min
ft3
limits prescribedfor
Thetest inlet total pressure is within the Table 3.1 limit of 5%.
(b) Inlet temperature departure
(T;)sp
( T;)sp
560.8
( N;)sp
10,000
S T D - A S M E P T C L O - E N G L L997
0 7 5 9 b 7 0 D b 0 5 5 0 2 304 9
36,000 38,000
"
0.06632
0.06997
36,000
100 = -0.049%
0.06632
(pdsp
(P ilsp
kr=
0.246
-=
CP
(cp -
Btu
Ibm "R
t Ibf
-)Ibm "R
Btu
)(
lbmole)
28.65 Ibm
1
47.26 psia
In (14.1 7 pia)
In
= 1.555
= 1.392
P)(5
(-) )
60hr
min
2.rrrad
-)
4t =
(0.06997 F
Ibm)(1 0,000 rev
(
18.4y
ft)
= 0.03996
min
(F)
ft Ibf
lbmole
(S
K)
(2.983 x l o 6 sec L
(32.174 *)ft Ibm
- o 555
[ 47 26 -- 1 1
0.555
(540.8 "R)
q)
= 0.4734
(d) Total work input coefficient using the shaft power method (see Eqs. [5.4.18] and [5.3T-21)
(59.5
(om) t
wotp0ATo
e)
min
(0.462
(42.44
Btu
m)
(31.O "R)
-) Btu
= 20.1 hp
minhr
--
T)
2)(h=)
-)mm hp
(38.000
T)
(2.983
lo6
(60
= 0.6052
(e) Total work input coefficient using the heat balance method (see Eq. [5.3T-11)
78
= 0.5996
ft3
0.6997
Ibrn
fi3
= 2.170
Btu
hasbeenused
(Thedesigndischargetemperature
to estimate cpd)
0.250
-)
Btu
Ibrn "R
)(
1 lbmole)
28.36 Ibrn
= 1.389
Polytropic exponent is found assuming equality of the polytropic efficiency at test and specified
conditions (see Eq. [5.4T-71)
79
= 0.7905
1.389
(-)0.389
= 2.823
2.823
nSP
= 0.4734
(2.983
(-)0.549
(1 545
lbmole
1 549
(pd)sp
= 1.823 = 1.549
1 O6
5)(L
sec
ft lb
"R
=)
32.174 ft Ibm
)(
1
28.36
Ibm
lbmole)
2.823
= 3.196
(560.8 "R)
(b) Capacity at specified conditionsis found using the definition of theflow coefficient and equating
the flow coefficients at test and specified conditions (see Eq. (5.4T-11).
0.03996 10,000
2)
(2.")
80
rad
18.4
ft3
(T
ft)3 = 9051 min
WSP
(E)sp
(pi)sp = (9051
min Ibm
ft3
,in)
(0.06631 F)
(60 7)
Ibrn
Ibm
= 600.3 hr
min
= 36,020
(cf) The specific volume ratio based on total conditions is (for information only)
Since this temperature is nearly equal to the design value of 844.1R, the average specific heat
chosen for the calculations is assumed appropriate.
(0 Gas power is found using the equality of the total work input coefficient between the test and
the specified operating condition. Using the shaft power method and Eq. (5.4T-201 gives
Usingtheheatbalancemethod,
Wsp (&b)sp
(PghJsp
Eq. [5.4T-201gives
(y)sp
c. u2 (600.3 x)
Ibm (0.5996)
-
33,000
(2.983
X lo6
-)
ft Ibf
(33,000mm hp
(60
= 1011 hp
$)
(g) Since the specified speed and the test speed are equal, the mechanical losses are assumed
equal. The shaft power is then
( p s h ) ~ h= (Pg,,,
+ Q , A P = 1021
or
(Ps&,
= (Pghb+
method)
(h) Static discharge temperature and pressure may be calculated from the
mass flow rate, flow
area, and total temperature and pressure. Since the flow Mach number is below 0.2,Eqs. [5.4.21,
[5.4.31,f5.4.41,
and [5.4.61 may be used.
81
(Tstatic d)sp
(Td)sp
- 2 ) gVdc c p
= 846.5
Ibm
(pstatic d)sp
(pd)sp
ft'
(130.5)'
- 2 778.1 7
sec
(Pstatic d)sp
( Ibf
-ft) Ibf
sec2
= 845.2
-)
Btu
--
Vd = 44.97 I bf
2gc (1 44)
"R
F)
(130.5)'
(0.1 40
ft Ibm)
2 32.174-
sec2Ibf
ft2
-
Ibf
( sec2
inn) = 44.7 in'
144
Checking
V=Ibstatic
\(0.140-Ibm) - (1;;2)
A/
ft3
(Pstatic d)sp
144
R
I = 130.5
II
(pstatic d)sp
(Trtatic d)sp
--
(1 44
-)
ft Ibf
(1545 ibm "R
82
ft2/
ft
sec
$) $)
(44.7
(Lm)
28.36 Ibm (854.2 "R)
= 0.140-
Ibm
ft3
OR
S T D - A S M E P T C LO-ENGL L777 U 0 7 5 7 b 7 0 0 b 0 5 5 0 8 8 2 2
TABLE C.l.l
CALCULATION SUMMARY
Test
Corrected
to
specified
Operating
Test
Operating
Quantity
Condition
Symbol
ConditionValue
Expected
at
specified
Units
Ibm/hr
38,000
36,000
36,000
rise 2. Pressure
*P
psi
33.1
(total)
3. Head
WP
ft Ibf/lbm 44,100
43,900
43,900
(Psh)sh
hP
1020 1040 1100
(Psh)hh
hP
1020 1030 1080
30.9
4. Shaft power
Shaft
(a)
method
method
(b) Heat
5. Polytropic
efficiency
ratio
0.80
VP
0.790 0.790
6. Flow coefficient
Mm
8. MachineReynoldsno.
Rem
(vi/v&
heat10. Specific
coefficient
11.work
Polytropic
0.0400 0.0400
2.17
1.39
1.39
pp
0.473
0.473
Pmin
0.600
0.600
ash
nhP
0.605
0.600
9 = (w/oJ
ft3/min
2.1 1
0.605
0.600
9050
9050
(bl Staticpressure
P
T
P
560 O R 560
14.0psia 14.0
O R 561
561
14.1psia 14.1
540
14.1
541
14.2
(b) Staticpressure
(c)Totaltemperature
(d) Total pressure
T
P
OR
psia
"R
psia
T
P
845 830
44.7 47.0
831
45.0 47.3
847
843.5
45.0
844.8
45.3
(pg)sh
( pg)hb
hP
hP
1o m
1O20
1010
1 O00
1060
10,000
Not applicable
19. Speed
rPm
10,000
10,000
Qm
hu
20.120.0
20.1
83
1 000
1200
1100 400
lo00SOO
900
700
1300
0.600
0.490
0.450
0.440
600
600
800
700
900
loo01200
1100
FIG. C.l(b)
IDEAL GASSPECIFICHEAT
84
FOR STEAM
1300
(1)
hatic i
-(
0.6220 Ibm
(1.456 psia)
Ibm da
(7.50 - 1.456) psia
= 0.2408
lbmole W
lbmole da
Ibm da
) + 0.2408 lbmole
= 26.04
Ibm
lbmole
86
(d) Air specific heat at constant pressureis found using dry air and steam properties. The specific
heat at constant pressure
for both the dry air
(da)and water vapor( W ) are given in Sample Calculation
C.1. (Fig. C.1)
(cplsp
mass da
Btu
+ 0.1498 Ibrn W O 448
Ibm da "R
*
Ibm W "R
1 .O00 Ibrn da + 0.1 498 Ibrn W
0.267
Btu
Ibrn "R
da
OR
1 .O00 Ibrn da
= 0.281
Btu
Ibrn "R
0.267
k)sp
=
(0.267
BtuBtu
- (1.986
G)
0.281
(kd)sp
(0.281
Btu
Ibrn "R
lbmole "R
26.4 Ibrn
Btu
Ibrn "R
-) Btu - (1.986
Ibrn "R
) ( lbmole ) =
Btu
lbmole "R
)(
lbmole
26.84
)=
Ibrn
(0 The inlet flange kinetic viscosity is found from Ref. (D.20) and is assumed to be that of dry air
at the inlet pressure and temperature
( ~ i ) s p=
4.00
87
1 O4
ft*
sec
(g) Static specific volume at inlet and discharge flanges is found using the ideal gas law
(Vstatic i)sp
(1545 ft lbf)
lbmole "R
5)
(1 44
= 31.98
$)
Ibm
lbmole
(1 103.0 O R )
26.84 Ibm
(48.00
ft3
$)(1 44 $)
= 9.186 Ibm
ft3
(h) Average fluid velocity at inlet and discharge flanges (see para. 5.4.3.1)
Ihm) (
(Vi)sp
(1 7,300 - -h r ) (31.98
hr 3600 sec
=
?r i16.876 , \ 2
88
$)
= 98.94
ft
sec
{i) Fluid Mach number at inlet and discharge flanges (see para.
(-)W V
""
sec
(Mi)sp =
=
(32'174
= 0.0798
lbmole
lb "R
80.61
(Md)sp
T SP
ft
-
98.94.
(1545
S)
2 778.1 7
) (26.84
lbmole
Ibm) (1 103.0 "R)
is foundusingtheenergyequationforan
sec/
ft Ibm
E)
(32.1 74 7)
(0.267 Btu
Btu
Ibf sec
(80.61
(Td)sp = 1103.0
"R +
= 0.0484
lbmole
lb "R
"R +
) (26.84
lbmole
Ibm) (600.0 "R)
ft
sec
ci,, Totaltemperatureatinletanddischargeflanges
adiabatic process (see Eq. [5.4.6])
(Tilsp= 600.0
5.4.2.5)
-)
(778.1 ft Ibf
7 (32.1
Btu
k)2
-)
= 600.7 "R
= 1 1 03.46 "R
"R
(k) Since the Fluid Mach number is less than 0.2, the total pressure may be calculated according
to the simplified Eq. [5.4.4]
89
7.53 psia
2(31.98$)(32.1741bfre;i)(144$)
ft Ibm
48.00 psia +
2 (9.186
$)
(32.174
ft Ibm
m)
(144
G)
= 48.08 psia
(I) Total density at the inlet and discharge flanges is found using the ideal gas law
(7.53
S) $)
Ibm
(1 44
(48.08
(Pd)sp = (&)sp
h Ibf
= 0.03136 ft3
"
)(
z) S)
1
Ibm
= 0.1090 ft3
(1 44
lbmole)
(1 103.5 "R)
"
26.84 Ibm
j=l
= 7.037
lo6 h2
sec2
(E),(3)
c u2 sp
CU2
= 9841 rpm
Nt = 10,000 rpm
as the approximate appropriate test speed. Not that no Reynolds number correction (as used later
in converting the testdata to the specified operating condition) is used in this estimation of the
testspeed.Also, note that the Code speed rule (para. 5.3.2) reduces to the equality of Machine
Mach numbers between the test and the specified operating conditions
for ideal gases with equal
values of the specific heat ratios.
See para. 6.2.7.
fa) Test run number 4
fbl Duration of test = 40 minutes
(c) Compressorspeed = 9,500 rpm
fd) Inlet temperature = TStatic
dbi = 540.0 "R
91
= 0.6220
0.240
Btu
Ibm da "R
0.3667 psia
Ibm W
= 0.01661
14.10 psia - 9.3667 psia
Ibm da
(540.0 - 530.0) R
(1095.5
0.01881
= 0.03024
Btu
- 38.35) Ibm W
lbmole W
lbmole da
92
Btu
(1092.2 - 38.35) Ibm W
lbmole da
I,
) + 0.03024 lbmole
Ibm da
lbmole da
= 28.65
W
W
Ibm
lbmole
(c) Air specific heat is found using dry air and steam properties
= 0.244
W "R
Btu
Ibrn "R
Ibrn da "R
1.O00 Ibrn da
= 0.253
I,
Btu
Ibrn da "R
+ mass W
Btu
Ibrn "R
(C,),
Btu
- = 0.244 +2 0.253 IbrnBtu"R - 0.249 Ibrn "R
('pi
"
93
ratio
Btu
0.244 (kilt
Ibrn "R
Btu
986
lbrnole "R
m)- (0.1
(0.244
Btu
0.253
0.1 986
Btu
Ibrn "R
)(
1 lbmole)
28.65 Ibrn
)(
Btu
lbmole "R
1 lbmole)
28.65 Ibrn
= 1.397
= 1.37
(e) The inlet flange kinematic viscosity is found from Ref. (D.20) and is assumed to be that of dry
air at atmospheric pressure and the existing temperature
= 1 .TOx 10-4ft2
sec
(0 Static specific volume at inlet and discharge flanges is found using the ideal gas law
(vstatic i ) t
(lS4'
ft Ibf
lbmole "R
)(
1 lbmole)
540.0 "R
28.65 Ibrn
"
(14.10
ft Ibf
5)
(144
)(
$)
lbmole)
= 14.34
ft3
Ibrn
1042.3 "R
"
(vstatic d)r
28.65 Ibrn
5)
(99.6
(1 44
$)
[g) Fluid velocity at inlet and discharge flanges (see para. 5.4.3.1)
94
ft3
= 3.919 I bm
(36,500
(Vi),
%)
hr
(1 4.34
E)
Ibrn
4(
16.8767
fi)'(3600 F)
= 93.60
-)hr
=
77
Ibrn
= 72.56
r)
10.020
4 (F
ft
sec
ft
sec
h) (3600
(h) Fluid Mach numbers at inlet and discharge flanges (see para. 5.4.2.5)
93.60
(Milt =
f-t
sec
lb
lbmole "R
72.56
(&)r
ft
sec
) (-S)
(540.0 "R)
28.65 Ibm
) (-
lbmole "R
(i) Totaltemperatureatinletanddischargeflanges
adiabatic process
= 0.0818
28.65
-)
(1 042.2
= 0.0459
"R)
Ibrn
is foundusingtheenergyequationforan
(93.60
0.244
E)Btu
(778.1
-ft)7Ibf
Btu
95
- = 540.7 "R
(32.1 74 ft Ibrn
Ibf sec2)
(72.56
(Td)l
1042.2 'R
+2
m)
(0.253 Btu
A)'
= 1042.6 "R
E)(32.174 ft Ibm
(778.17
Ibf secz)
Btu
(j) Since the Fluid Mach number is less than 0.2, the total pressure may be calculated according
to the simplified method of Eq. [5.4.41
Pl
193.60
\-
2 14.34
(Pd),
(&I
= 99.6 pSa +
f t z
-\secl
-)Ibm (32.1 74 m
Ihm)(144 $)
ft3
14.1 7 psia
= 99.74 psia
)
(
inn)
2 (3.91 9 2)
(32.1 74 lbfsec~ 144 ft Ibm
Total density at the inlet and discharge flanges is found using the ideal gas law
S) $)
(14.17
ft Ibf
(1545 lbmole "R
(Pd)sp
)(
lbmole)
(540.7 )
"
28.65
(99.74
(&),
Ibm
(144
Ibm
E)
(144
= 0.06993
fi3
O R
$)
Ibm
= 0.2555
ft3
z)'
I6 (20 in)'
+ 4 (1 8 in)']
= 6.35
j=1
lo6
ft2
5ec2
mm
S T D e A S M E P T C 30-ENGL 1 7 7 7 II 0 7 5 1 b 7 0 Ob05522 L T 2
Theshaft powercanalsobecalculatedfromthe
and Eq. [5.5.14]
gas powerusingtheheatbalancemethod
Equations L5.4.171 and [5.4.181 show the parasitic losses to be mechanical losses. Also, using
-)
E)
(1042.6 - 540.7) "R
Ibm "R
hr
(psh)t
(42.44
&) F)
x) -
-)
= (1791.4
(60
(143.1 min
Ibrn (0.462 ,bm
(568.5
6740
(42.44
Btu
hr
S) T)
(60
525.0) "R
+ 2.65 + 67.8) hp
(n) The gas power can be calculated from the heat balance methodas done above to get
(P,),
(W C p ) t (Td - Ti), + Qr
= (1791 4 + 2.65)hp
= 1794 hp (heat balance method)
Thegas power can also be calculated from the shaft power using the shaft power method
(P,),=
(Psh)r
- W o Cpo ATO
= (1 862 - 67.8) hp
= 1794 hp (shaft power method)
- (piIr x
(pi)sp
100 =
7.53
- 14.16 X
7.53
in the inlet
100 = -88.0%
which is outside the range of the Table 3.1 limit of 5%; therefore, the test is not a Type 1 test.
Therefore, wemustconduct
a Type 2 test; however, we will verify that this i s aType 2 test,
.e.,satisfies the Table 3.2 requirements.
97
S T D - A S M E P T C 10-ENGL 1997
0 7 5 9 b 7 0 0605523 O34
0.249
fi Ibf
Btu
Ibm "R
Btu
)(
lb mole)
28.65
Ibrn
(99.74 psia)
14.1 7 psia
= 1.507
(99.74 psia) (540.7 "R)
In
(14.17 psia) (1042.6 "R)
2)(L
(-) )
(36,500
rad
(0.06993 F
Ibm )(9,500
S)( 1 2 ft)
98
hr
60 min
20.0
= 0.03148
3
= 1.386
075'1b70
Ob05524 T 7 5
0.507
)(
28.65
L
F
Ibm )
(540.7"R) [(-)99.74 1.507 - 1 1
14.17
ft2 ) ( 1
Ibf sec2)
6.35 lo6--
(-)
1 507
(1 545 lbmole
ft Ibf "R
0.507
(143.1
(Q,,,)t= wocpJT0 =
?)
32.174 ft Ibm
sec2
(0.462 G)Btu
mm
(42.44
= 67.8 hp
-mln
) hr
Btu
Ibf ) (60
mrn hp
F)
(36,500
W-
gc
(6.35 x
F)
= 0.4901
lo6 E)(LF)
sec2
32.174 ft Ibm
(e) Total work input coefficient using the heat balance method (see Eq. [5.3T-l])
(0.249
-)llbmBtu"R
(1042.6- 540.7)"R +
ft'
(6.35 lo6 -()-
sec2
(0
Btu
36,500 hr
= 0.4935
32.174 ft Ibm
b"h
- Ti) I - (0.249E)Btu
CU'
-
cp (Td
gc
99
*)Btu
= 0.4927
c)r
Ibm
-
0.2555
(;)r
ft3
0.6993
I bm = 3.65
ft3
Btu
(Thedesigndischargetemperaturehasbeenused
Average specific heat ratio
to estimate cps)
0.274
Btu
Ibrn "R
Btu
Ibrn "R - (1.986 lbmole "R
-)
Btu
Ibrn "R
Btu
)(
lb mole)
28.36 Ibrn
1
= 1.370
(F)
= 2.73 x 105
SP
NDib
T)
(=E)
20 rad
sec
1O0
(E)
1 .5
ft
= 6.10 x
105
or
Remt < 4.775 Remsp = 4.775 (2.73x
Remt > 0.2094 Rem,
= 0.2094
lo5) =
1.3 x
lo6
(2.73
x 1Os) = 5.72 x lo4
where
= 0.988(2.73 x 105)4.243= 0.04718
0.988
(6.10
x 105)-0.243
= 0.03881
RA,=
0.066 + 0.934
RA,=
0.066 + 0.934
RA,=
0.066 + 0.934
46 X 10'
Rem
1,
RCsp
4 (I.5 in) x 1 os
2.73 X 105
46 X
On4718
= 1.0354
lo5 RCt
Rem
Il
4 (I.5 in) x 1 o5
6.10 X 105
0~13881
= .O99940
0.000125 + Reml
RBI =
13.67
( Remt)
log c +
13.67
6.10 x lo5)
= 0.9961
13.67
log 0.00012 +
6.10 x 105)
1 o1
Then,
-(
~ =) 11~
1 .O354 ) (0.9967)
= 0.1785
-p0.828) (0.99940 0.9961
and
0.822
(vsp)sp
or
3.054
2.054
nsp = -- 1.49
The polytropic work coefficient ratio for the specified operating condition is
0.822
(pplsp
(qp)r
Discharge pressure ratio is found using the definition of the polytropic work coefficient to give
[(E)
(7.037x
r.4044
(1 545
(-Ibf secz)
f
i
l
z
)1
lo6 -
3.042
= 6.400
)(
'
3
)
(b) Capacity at specifiedconditions is found using thedefinition ofthe flow coefficient and equating
the flow coefficients at test and specified conditions (see Eq. [5.4T-11)
z)
= 0.03148 10,000
102
rad
fi)
20
(277~)
fi3
= 9157- min
= 17,230-Ibm
hr
(2) = [
Vd sp
(F)q
1
1
-
= 6.4007.49 = 3.48
*P
(E)7]
n-1
(TdlSp
= [Ti
0.49
Since this temperature is nearlyequal to the designvalue of 1103.5"R, the average specific
heatchosen for the calculations is assumed appropriate.
ffl Gas power is found using the equality of the total work input coefficient between the test and
the specified operating condition. Using theshaft power method, Eq. [5.4T-20], and Table 5.3 gives
Usingtheheatbalancemethod,
33,000
(1 7,230
F)
min
ft2
) (1
Ibf sec2)
-
sec2
32.174 ft Ibm
0 7 5 9 b 7 0 Ob05529 557
(g., The shaft power is found by assuming the mechanical losses are proportional to a power of
the rotational speed (see Eq. 15.6.81)
N
(2)
=
2.5
(Qmlsp
(Qm)r
67.8 hp
10,000
rev
-
2.5
= 77.1
hp
9,500 min
or
(h) Static discharge temperature and pressure may be calculated from the mass flow rate, flow
area, and total temperature and pressure. Since the flow Mach number is below 0.2, Eqs. 15.4.21,
(5.4.33,i5.4.41, and 15.4.61may be used.
With a guessed velocity of 80.2 fthec, obtained by trial anderror,
(802)2
= 1106.
sec2
ft2
-)
ft Ibm
2 778.17- Ibf) (32.1 74
(0.281
Btu
Ibf sec
= 48.2
Ibf
--
(0.1 O9
F)
= 1 105.5 'R
104
ft* e c ~
(80.2)2~
ft Ibm
2 32,174
Ibf sec2
Btu
Ibf
S T D - A S M E P T C LO-ENGL L797
0 7 5 7 b 7 0 Ob05530 2 7 7
Checking
ft
sec
= 80.2 -
(144
(1545
$) $)
Ibrn
(48.2
ft Ibf
-)Ibm
('W)
28.36 Ibm
= 0.109 -
(1105.5 "R)
ft3
O R
,oo -
3.65
3.48
100 = 1.049%
F)
(36,500 hr
-)
(0.06993 Ibm
ft3
(9,500
x 100 =
E)
min
X
F)
(1 7,300 hr
")
(0.03 136 TI bm
fi) (1 0,000 min
105
100 = 99.6%
Thetest capacity
speed ratio is within the Table 3.2range
(c) Test Machine Mach number (see para. 5.5.1)
of 96% to 104%.
= 0.710
Mm, =
(9,500
Mmt
2)($
(-
h) (2
TE)(Az)
= 0.708
1
moe
26.84 Ibm
(540.7 R)
2)
Rem,
of 0.625(0.710
= 2.73
- 0.085)to 0.815
105
sec
?)
(9,500 mm
(hg)
(S (gfi)
h)
= 6.10 x 105
Remt =
Remt
6.10 x 105
x 100 =
2.73 x 105
Remsp
100 = 223.4%
TABLEC.2.1
CALCULATION SUMMARY
Test
QuantityValue
Units
Symbol
Test
Corrected
to
Specified
Operating
Condition
Expected
at
Specified
Operating
Condition
lbmhr
36,500
17,230
17,300
2.Pressurerise
psi
85.6
40.7
40.6
3. Head (total)
ft . Ibf/lbm
80,400
88,450
88,200
4.Shaftpower
(a)Shaft method
(b) Heat method
hP
1851
1862
1010
1016
1025
1025
5. Polytropic efficiency
0.828
0.822
0.82
6. Flow coefficient
0.031 5
0.031 5
0.031 6
0.724
0.703
0.703
610,000
273,000
273,000
3.48
3.48
3.48
1.39
1.37
1.37
0.408
0.44
0.493
0.493
0.490
0.494
0.490
0.494
8700
9160
9190
540
14.1
S41
14.2
660
7.50
601
7.53
600
7.50
601
7.53
psia
1042
99.6
1043
99.7
1106
48.1
1106
48.2
1103
48.0
1103
48. I
hP
hP
1794
1794
933
939
hP
2.65
19. speed
rPm
9,500
10,000
hP
67.8
77.1
14.Capacity
15. Inlet gasstate
(a) Static temperature
(b) Static pressure
(c) Total temperature
(dl Totalpressure
16. Discharge gasstate
(a) Static temperature
(b)Staticpressure
(cl Totaltemperature
(d)Totalpressure
hP
Whin
OR
psia
OR
psia
"R
psia
O R
107
10,000
At an inlet
flow
of 3000 ft3/min
= 90 psia
= 0.76
Shaft power = 690 hp
Thedata in the left hand column above indicate the specified operating conditions. This data
describes the compressor geometry, the operational speed, and inlet gas conditions. The data in
theabove right hand column describestheintendedperformanceofthecompressor
at the
specifiedoperatingconditions.
It is thepurpose of thetest to verify these intendedvaluesor
establishtheactualvalues.
It i s assumed thatcircumstances prohibit testing with methane. Air i s available at 14.7 psia,
520"R, and 50 percent relative humidity. The driver has variable speed capability.
The following assumptions are made to simplify the calculation process so that focus may be
madeondemonstrationpoints.
(a) Both the test gas, air, and the specifiedgas, methane, will be treatedas ideal gases with constant
specific heats. Average values will be used. (The alternative is to use actual gas thermodynamic data
and the Type 2 calculation procedure. This would lead to slightly more accurate results.)
(b) Leakages will be assumed negligible at both test and specifed conditions. The rotor mass flow
rate is then the inlet mass flow rate.
The test speed required to provide equivalence betweentest and specified conditions is obtained
fromthe speed selection rule. For ideal gases,
1o9
TABLE C.3.1
PRETEST C A L C U L A T I O N S U M M A R Y
Methane
CaS
Pd
Pi
90.
30.
psia
psia
Ti
"R
570.
ft-lbf/lbmOR
96.31
1.28
0.078693
0.769 x
16000.
800.
0.532
3.411 X lo5
0.76
k
P
P
Ibm/ft3
Ibmh . sec
ft/sec
rPm
Mm
Rem
TP
1.6
pd/pi
n
qi/qd
51.597
14.7
520.
53.53
1.396
0.0761
1.27x 10-5
12704.
635
0.5675
1.583 X lo5
0.76
(excludes Reynolds
Number correction)
3.51
3.0
1.4
2.19
0.00343
Air
2.1 9
with
1
and
n
n-1
"
'Ip-
k - 1'
and Remco,, = 1
Precise values of pressure ratio, efficiency, and polytropic exponent for both specified and test
conditions are of course unknown before test. However,theappropriatetestspeedmaybe
estimated by making the following assumptions:
(al Thepressure ratio and efficiency at specified operating conditions are equal to the design
values.
(b) The efficiency at test conditions is also equal to the design value.While the Reynolds number
effect might be taken into account here, it is small and the current calculation is only an estimate.
It is ignored simply for computational ease.
The firstassumptionallows
calculation of thespecified conditionpolytropic exponent.The
second allows calculation of the test polytropic exponent. With these a test pressure ratio estimate
and a required testspeedestimatemaybecalculated.This
speed may be used to calculate
Machine Mach andReynoldsnumbers.
Thegasdataused
andresults of the computations indicated abovearesummarized in Table
C.3.1. Thevalues in thistable may be used to determine if it is possible to accomplish the
proposedtest within the allowable deviations in similarity parameters.
Mach NumberCheck:The test Mach number is = 6.6 percentgreaterthanthedesign
Mach
number. This is an unavoidableconsequenceof
gas selection with
110
different k valueswhenspecificvolume
ratio equality i s maintained.
Thedeviation is, however, within the limits of Fig.3.2.
ReynoldsNumberCheck:ThetestReynoldsnumber
is -46 percent of thedesignReynolds
number.This
is within thedeviation
limits of Fig.3.4,
andthe
correction relationship applies. The correction
has not been applied
to the tabulated values,sincethe computations are preliminary.
The compressor is run to obtain a bracketing point. A bracketing point lies within 2 4 percent
of the specified operating condition flow coefficient of interest, which is
3000
= 0.03427
(y)3
12
2~(16000)
The desired test inlet flow may be calculated from test and specified operating
coefficient equality, which yields
9;,= 9i
'P
(9
NS
= 3000
12704
(-)16000
= 2383
condition flow
ft3
rnln
RH;= 50%
50.4 Ibf/in2
Td= 832 "R
20 hp (from lubricating oil temperature rise and flow rate)
5574.5 Btu/hr(calculatedcasingheat
loss)
Psh= 339. hp (shaftpower,perhaps from a torquemeter)
N= 12690 rpm
R = 53.53 ft-lbf/lbm*"R
The next step is to compute the following dimensionless parameters from the test
SpecificVolumeRatio:
pd=
om=
or=
(3;
1
r", =
= 2.14286
Flow Coefficient:
111
data.
(-)n -n 1
ClPt =
RF
xu2
- 1 1 = 0.62702
Work InputCoefficient:
PolytropicEfficiency:
gc
MachineMach Number:
Mmt =
U
-
Remt =
- = 1.583 X lo5
Machine ReynoldsNumber:
Ub
V
whichhavebeenevaluated
using
S T D - A S M E P T C LO-ENGL 1777
W;
0 7 5 7 b 7 0 Ob05537 b 2 3
N = 12690 rev/min
D = 11.459 in.
R = 53.53 ft-lbf/lbm*"R
CU2
k = 1.396
Btu
Ibm "R
cp =
113
ft
lo4 Ibm
14.7
p i = - -Pi
=
R Ti
53.53
Ibf
(5)
144
):(
ti Ibf
(-)Ibm "R
2.9595
Wi
q,=-=
Pi
(20hp) 33000
o
m
.... =
f-\ti
778.1 6
u =" 2lrN
60 24
- (2.-) (12690;)Rad
rev
mm
1.27x
10-5
y="=
ft3
520 "R
(E)
60
sec
0.076047
Ibm
= 0.076047 -
(=)
min
I bm
(F)
(-)minft Ibfhp
Ibf
Btu
min
(--)601 min
sec
( r Ibm
) sec
11.459
(?in.)
= 1.67 X
Ibm
(-121-)fin.t
min
= 848.2 -
0.076047 (F)
b = 0.5 in,
ft3
= 2335 -
(--)121 fin.t
= 634.5
-)secft
ti2
-
sec
= 0.041 7 ft
114
or
6.41 x
lo4 5
1 -
VPsp
= (1
- VPt)
RASPRBsp
Rem, 5 1.81 x
lo6
0.97798
(0.99648)
= 0.2476
184 (0.9971 8)
- (1 - 0.744)1 .O1
where
RAI = 0.066
RB,
(4.8X 1 O6 b) RCsp
Rem,
+ 0.934
= 0.97798
RCr
= 1.01184
= 0.99648
with
RC, = 0.988 Remsf-0.243= 0.044696
RCt= 0.988 Rem,p-0.243 = 0.053862
b= 0.5 in.
Remsf= 3.411 x lo5
Remr= 1.593 x lo5
E= 0.000120 in.
so,
PSP
0.752
Pt
0.744
115
drp
p*p
= 4, =
qpt Remcorr=
0.3363
0.744 (1 .O1
18)
= 0.752
= [0.6344
1.11006 x lo5
((3.4395)
) + 11
96.31 (570)
3.4395
= 2.9750
where
n
1.28
116
0.7524 = 3.4395
S T D O A S M E P T C 20-ENGL 2 7 7 7
0759b70 Ob05542 O 5 4
(3,
Tdsp =
($sp
D 3
gi = [ 4 s p 2 7 ( ~ )
SP
(3"
n-1
flow coefficient
= 0.03363 2
min
wj = pigi =
0.4099
11.459 3
ft 3
16000 rev 1 2 ft) = 2944 rev
min
[96.31
(570)
ft3
Ibm
2944 - = 231.7 mm
mln
F
I
117
rad
7-
The power requirement at the specified operating conditions is determined from the total work
input coefficient.
Ibm
sec
sec
min
ft Ibf
1.11006~lo5Ibm
= 3.862-60-0.849
min
hp
(G=)= 697.3 hp
where,
Om,
-="
33000
Om,
Nsp)2'5
= 20 hp
33000 Nt
In this example both the shaft method and heat balance method give exactly the
same power
requirement. This may not always be true because of differences in the independent measurements
which areused.Thisexamplewas
specifically constructedusingvaluessuchthatthepowers
would match.
118
ft3
3000 min
d =
= 0.03427
in Fig.C.4.1.Thedata
points are valid
which fallsabout midway betweenthedatapoints
bracketing points in that they are well within 4 percent of the flow coefficient of interest (Table
3.2). In the absence of additional data points, the values of the dimensionless coefficients at the
flow coefficient ofinterestaredetermined
by linear interpolation. Linear interpolation gives
These values are used to calculate the compressor performance in dimensional terms as follows:
Flow rate: 3000 ft3/min as above
Dischargepressure:
1.11006 x l o 5
(3.430196.31(570) +
119
'1,
3.4304
= 2.9497
TABLE C.4.1
Point
Data
1st
Units
Test Data
Data Point
12690.Irpml
3.0799 [lbm/sec]
N
Wi
Pi
Ipsial
520.
Ti
PR1
RH
[%I
Pd
Td
[psial
Om
psh
QI
12690.
2.9595
14.7
520.
50.
50.4
832.
20.
339.
5574.5
53.53
1.396
IORI
20.
lhpl
Ihpl
330.39
5495. [Btuhrl
[fi-lbf/lbmoR]
1.396
14.7
50.
49.4
828.
53.53
Calculation Summary:
2nd Data Point
Test
Specified
Opentill#
Operating
Operating
Conditions Conditions
Test
opentill#
conditions Conditions
0.034 4
0.843Ph
0.035
PP
0.838
0.794
0.748
'IP
0.849 nhb
0.84900
0.5320
Mm
ah
Rem
q&d
0.035
0.034
0.843
0.627 0.61
0.739
0.744
0.849
0.84900
0.5674
0.5675
1.583 x Id
2.1429
0.832
0.634
0.752
0.832
5
0.5320
3.411 X 105
2.1351
2.
1668
1 .S83 x 1Os
2.1105
where
(E)
k
tsp
= 0.750
(z
= )
1.28
3.4304
n = 1.41145
and,
Pd
specified
= (pd/pi) pi =
(2.9497)30
120
88.49 lbf/in2
0.623
3.41 1 x 1Os
S T D * A S M E P T C L O - E N G L L977 M 0 7 5 9 b 7 0 Ob05545 7 T T M
FIG. C.4.1
Power requirement:
Heat balance method:
=+
"- wi n h b
PShhb
xu2
33000 g,
Qml
33000
where
121
S T D - A S M E P T C 10-ENGL L997
ASME
10-1997
- (236.08)0.823(1.11006
x 1Os)
33000
t 35.7 hp
= 689.43 hp
Calculated
16000
Methane
16000
Methane
570.
570.
30.
3000.
90.
0.76
690.
30.
2.1 9
2.1 52
3000.
88.49
0.750
705.79
Comparison of the design and actual results indicates that the compressor falls short of meeting
its designpressuregoal at design flow. The implication of this fact to the parties involved are
beyond the scope of this Code, as they would be also had the compressor exceeded its design
goals.However, typical industrial reaction in lieuof mutualacceptance as tested is hardware
modification or specified condition speed adjustment, In the event of hardware modification the
testmust be repeated.Forsmallspeedadjustmentsthetestresultsmayremainvalid.This
is
determined by conversion of the testresults to the new specified condition speedand verifying
that the limits in departurebetweentestandspecified
conditions are not exceeded.
122
STDmASME P T C 1 0 - E N G L
1 9 9 7 D 0 7 5 9 b 7 0 06055i7 5 7 2 D
l It is recognizedthatthere
is apotentialenvironmentalproblem
demonstrate the calculation method.
123
of usingrefrigerant
22. Theusehere
is only to
TABLE C.5.1
SPECIFIEDOPERATINGCONDITIONSANDPREDICTEDCONDITIONS
Inlet
Pressure,psia
Temperature, O R
0.2602
Specific volume ft3/lbm
0.7981
Z
373
0.01
Viscosity centipoise
0.6266 heat Btu/lbm-"R
Specific
Specific heat ratio
820. Sonic velocity hhec
209.8
Enthalpy Btu/lbm
Entropy Btu/lbm-"R
Discharge
200
575
0.7578
0.8768
0.01
0.4894
1.1 28
830.
164.9
1
704.8
021
.S77
650
1 .O98
1 .S92
in.
Impeller
2nd
Diam., in. 36
Width,
Tip
f, in.
st
36
36
36
3rd
5th
36
1.25
4th
-
1 .5
0.000125
for R134a.Thetestspecific
volume ratio comparedtospecified
indicated that the R134agas
conditions are very close (within the accuracy of estimated gas properties) to that of the specified.
Furthercheck on the assumed efficiency also indicated it was within 4 percent of specified.
The polytropic headwas computed along with Reynoldsnumber correction factorand a new
speedwasalsocomputed.Furthercheck
of this speedagainst the critical s p e e d of the unit
indicated a margin of 8.6 percent, which should be within a reasonable range for unit operation;
therefore,R134a could be used.
Table C.5.6 has the basic R22 checkdata.Thetestspecificvolume
ratio is considerably off
from that specified.Thetest speed is 2556 rpm which is too close to the first critcal speed.
This problem demonstrated the extent of calculation necessary to come up with the test speed
for a given unit. The final test power may be increased by changing the inlet pressure and then
re-computing all the values.Effectivelythetemperature
ratio shouldremainconstantand
test
speedmay vary slightly with increase of inlet pressure.
Thetestspeed computed is only an estimate. Once the unit is on test, the 9/N should be set
and the specific volume ratio,
r", checked from test data. If the volume ratio is not correct, the
testspeed should be adjusted and the
q/N reset.
124
~~
TABLE C.5.2
GENERALFLOWCHART
FOR TESTGASSELECTION
Y
I
+
+
P. f
Estimate New
Discharge Pressure
Compute Test
Polytropic Exponent
Compute Corrected
Test Speed
Temperature andPressure
1
Verify Test
Temperature Pressure
Lower Than Max. Design
No
1.
125
Compute
Mach
No.,
Verify Test Volume Ratio
No
TABLEC.5.3
TEST GASINLET CONDITIONS
560
560
pi, psia
Ti, 'R 560
vi, ft3/lbm
0.982
p, centipoise
c, Btu/lbm-'R
k
ai, ft/sec
204.8
721.9
102.
673.8 TO 'R
547.7
590.3 po psia 1069.9
MW
44.01
h&Btu/lbm
si
Xi
Yi
R134a
c02
N2
20
560
10.73
1 .o0
0.975
0.017
0.2499
1.396
1 178.
227.4
493.
28.01
20
20
6.778
0.993
0.015
0.21 03
1.273
894.
121.2
0.02 -
R22
20
2.871 6
3.41
0.01 o9
0.2098
1 .O98
538.8
0.01 1
0.161
1.166
598.5
122.3
0.07
1 .O27
1 .O3
GENERAL NOTE It is recognized that thereis a potential environmental problem of using refrigerant22.
The use here is only to demonstrate the calculation method.
TABLE C.5.4
FIRST PASS FORGASSELECTION
specified
Ca
3.25
2.912
1.226
1 .O273
0.781
0.509/1 .O56
1.1 50b.287
Real
1.1 027
2731 0.0
2.266 x lo7
0.681
565.5
CO2
N2
R22
5.358
2.912
1.840
4.309
2.91 2
1.495
3.296
2.91 2
1.157
3.551
0.781
0.781
0.01
1 .o1
0.781
0.08
1 .O3
0.781
0.02
1 .O3
2.91 2
1.261
Ideal
Real
Real
Real
1 S74
71 422
1.366
34860
1.185
13750
107.2
1030.
[Note (111
86.2
836.6
1.116
1 o499
3.473 x 106
0.65
350
65.9
648
[Note (1 11
5822
4067
[Note (111
[Note (111
3600
R134a
0.983
2232
70.0
706
2554
[Note (311
NOTES
(1 1 Test values exceed the mechanical designlimit for the tested unit.
(2)No Reynolds number correction or verification of volume ratio, efficiency, or real gas correction,
(3)Test speed too close to rotor critical speed.
126
22
S T D - A S M E P T C 10-ENGL L997
0759b70 O b 0 5 5 5 1 T T 3
TABLE C.5.5
Specified
CaS
Polytropic work factor f
RemdRem,
1.004
2.912
0.781
2731O
Remcon
N, rpm
3600
0.681
22734
31790
Mm
9
hP
R134a
1 .o1
0.1 53
0.1
2.92
0.779
10605
1.O03
2247
0.655
14190
2039
Supplement C.5.A
PredictedConditionsSpecifiedGas
PressureRatio rp = pd/pi = 650/200 = 3.25
Volume Ratio r, = Vi/Vd = 0.7578/0.2602 = 2.912
kmaxkmin
= 1.1283 /1 .O975 = 1 .O28
S/N
= 2273413600 = 6.315
CheckSpecifiedGasforTypeofCalculation
FromSchultzCharts
Xmi"
Xmax
UseRealGas
CalculationMethod forSpecified
Gas
= 830.2 ftfsec.
127
Tip Speed
(36)
(3600)
fi
= 565.5 sec
720
MachineMach Number
Mm =
U/&
= 565.5f830.2= 0.681
Rem = Ulbl/pv
,U
021
/(1488.2)
= 6.86 x lov6Ibm/ft-sec
2.5
(565.5)(E)
Rem =
6.86 x
= 2.266 X
(0.7678)
lo7
Ideal Gas
Y = 1.0
= 0.0
f = 1.0
ComputePolytropicExponent
Remco,, = 1 .O
=-
(1.396
- 1)
1.396
0.781
= 0.3632
nt = 1.574
Compute TestGasPressureRatio
rvt = r,
SP
1.574
= 3.25
1.1027
n-1
= 5.358
129
W, =
W, =
=pi
vi
n-1
n - 1) 144
(rp
1 .S704
1.5704
- 1) 144
- 1 (20) 10.725 (5.358*3632
= 71422 ft-lbf/lbm
TestSpeed
Nt = 3600
NOTE Testtemperatureexceedsmechanicaldesignlimit.
Test Gas
= 5822 rpm
Co;!
Initial Estimate:
Assume NonidealGas
Useinletconditionsfor
initial calculations.
ReducedTemperature = R, = Ti /Tc,it = 560/547.7 = 1.022
Reduced
Pressure
= Rp = pi
= 20/1069.9 = 0.0187
x = 0.01
Y = 1.01
Z = 0.993
cp = 0.2103
130
ComputePolytropicExponent
n=
1
Y - m (1
+ X)
1
(
778.1 7 (0.21 03) (44.01 0.781
0.993 (1 545)
+ 0.01)
= 0.2749
n=
1.01
Compute TestPressureRatio
Pdr
z T ~ , 0.993 (-)44.01
= -1545
"'t
P i(
144 (20)
131
(560)
= 6.778
h119
n1
Wp,
Wpt=
= -pi .(rp
" - 1) 144
n-1
(1.3655
1*3655 ) 20 (6.778)(4.304.2677- 1 ) 144
-1
= 34860 ft-IbWlbm
Test Speed
= 3600
NOTE:Testtemperature
-- 4067
rpm
is marginal.Test s p e e d exceedsmechanicaldesign.
Assume NonidealGas
UseInletConditions
ReducedTemperature =
RI
= Tiflcrit
= 5601673.8 = 0.8311
ReducedPressure = Rp =
pdpcrit
= 201590.3 = 0.0339
X = 0.07
132
Y = 1.027
ComputePolytropicExponent
--
1
(
+ 0.07)
778.1 7 (0.2098) (1 02) 0.781
0.975 (1 545)
= 0.1 221 5
1
Y - m (1
n =
+X)
1
1 .O27 - 0.1 221 5 (1
+ 0.07)
= 1.1157
Pd
fp,
fp
1.1157
-
3.25'.'02' = 3.2955
pi
r( =
r6
= (3.2955)0.'22'5= 1.1568
vit
= --
1545
(K)
(560)
= 2.8716
144 (20)
P i,
W, =
0.9753
..n - 1.
pi vi (rp
n-1
- 1) 144
) 20 (2.871 6) 3.29!jS(*)
(1.1157 - 1
57
- 1) 144
fi-l bf
Ibrn
= 10499 -
Compute TestSpeed
Nt = 3600
J"=
= 2232 rprn
Vd
647.8
- 0.951 (1 545)
144 p M W - 144 (65.911 102
ZRT
= 0.9832
r, = 2.871610.9832 = 2.921
ComputePolytropic Work (RealGas)Factor
Isentropic
P= 65.91 psia
T'= 168.7 "F (628.7 OR)
V = 0.9205 ft3/lbm
h'= 135.46 Btu/lbm
n,= Inrplln r,'
r:= VdvJ = 2.8716 /0.9205 = 3.1 196
n,= In 3.2955/1n3.1 196 = 1 .O482
134
lb
ft3
Ws=
ns
ns- 1
(Pd
Vd
- pi
Vi)
144
10241
f=
1.O482
1.O482 - 1
= 1.01
ComputeReynoldsNumberandReynoldsNumberCorrection
Use Preliminary TestSpeed
Rem = Ubfp v
U= rDN/720
= IT 36 (2232)/720
= 350.6 ft/sec.
3 5 0 . 62 (.5~ )
Remt=
= 3.473 x 106
ReynoldsNumberRatio
Remt
"
2.266 x l o 7
Rem,
Allowable Ratio
Remt/Remsp 2 0.1
Therefore,theReynoldsnumber
ratio of 0.153meets
ComputeReynoldsnumberconditions.
RA = 0.66
conditions.
lo6 bIRC
+ 0.934 (4.8 xRem
RA,=
0.66 + 0.934
RB,=
1 .O
RC,=
0.988
= 0.02543
[3.473 X 1O6lo."3
0.66 t 0.934
R
A
=
I
RBI= 1.0
(1
- S&=
(1
- .781),
1 565 1 O
(-)1.548
(-)1.0
= 0.224
qPr= 0.7786
(z
p ()
r p 7n
"P=
1.1157
n-1
1) 144
0.1 157
.O11
(20) 6)(2[3.2955'."57
.871
- 1] (144)
= 10605 ft-lbf/lb
Correct Preliminary Test Speed
NI
= 3600
= 2247 rpm
136
= -0.026
TP,=
?Psp
Remco,,= 1.O
f= 1.0
Usereal gas calculation
Useinletconditionsfor
initial estimate
ComputePolytropicExponent
m=
(-
0.982 (1 545)
2
+ 0.02) = 0.1821
778.17 (0.161) 86.48 0.781
n = 1/[Y - m (1 + XII
t~=1/[1.03 - 0.1821(1
n= 1.1 845
+ 0.0211
Compute TestPressureRatio
'P(
1.1 845
3.251.1027
= 3.546
137
Compute TestTemperatureRatio
rtt
Compute Test
= rpT = 3.546L0~'82''
= 1.259
Head
Z Ti,
1545
0.982 86.48 (560)
(-)
=P it -
144 (20)
= 3.4113
n- 1
1.1845
= 13750 ft-lbf/lbm
Compute TestSpeed
= 3600
= 2554 rpm
Testspeed
is tooclosetothefirst
138
-1
S T D - A S M E P T C L O - E N G L L997
0 7 5 9 b 7 0 Ob85563 715
REDUCTION
A mixed hydrocarbon compressor which wasset up in the Sample Calculation C.5 was tested
on refrigerant134ausingaType2test.
Table C.6.1 outlinestheconditions for which this compressorwasdesigned.
It tabulatesthe
inlet anddischargeconditions, as well as theisentropicdischargeconditions.
Table C.6.2 showsthe gas compositionand critical properties for this hydrocarbonmixture.
The molecular weight and the calculated
gas constantareshown in this table.
Table C.6.3 shows thederiveddesignfunctions,
specifically pressureratio,temperatureratio,
volumeratio, as well as polytropicexponent,volumeflow,
head, efficiency, andpower.The
test is supposed to verify thesevalues.The
calculationof thesevalues is shown in Sample
Calculation CS.
The unit was tested on refrigerant 134a.' The test data are shown
in Table C.6.4. It wasat a
test speed of 2245 rpm. The inlet pressure was held at 20 psia at an inlet temperature of 100OF.
Dischargeconditionsachievedwere
67.5psiaand187.4"F.Thedatashown
is the average of
the actual test data readings. It is assumed that all scatter was within the allowable test requirements
for thesedata point positions.
The derived test functions, pressure ratio, temperature ratio, volume
flow ratio, etc., are shown
in Table C.6.5.This is the reduced data from the test point of Table C.6.4.The calculations are
shown in Supplement C.6.A for obtainingeach of the individual items.
TABLE C.6.1
SPECIFIC DESIGNCONDITIONS
Discharge
Isentropic
~~
650
Pressure,psia
227.7
Temperature
650
244.8
0.2465
Specific volume, ft3/lbm
0.7749
Compressibility factor
P
T
V
0.2602
0.7981
Z
Viscosity, centipoise
P
0.6266
Specific heat, Btu/lbm-OR
CP
Specific heat ratio
k
a
Sonic velocity, ft/sec 820
199.05
209.84h
Enthalpy, Btu/lbm
1.577
1.592S
Entropy,Btu/lbm-"R
200.
115
0.7578
0.8768
373
0.01
021
0.01
0.4894
2831.1
830
164.9
1.577
1 .O975
It is recognized that there is a potential environmental problem of using refrigerants. The use here is only to demonstrate
the calculation method.
139
TABLE C.6.2
GAS COMPOSITION A N D PROPERTIES
Composition:
Methane
Ethane
Propane
N-Butane
20%
25%
50%
5%
Critical properties:
PC
TC
--
v,
Molecular weight
Gas constant R
646.4 psia
577.2 "R
0.7943 fi3/lbm
35.67
43.31 ft-lbf/lbm O R
TABLE C.6.3
DERIVED DESIGN FUNCTIONS
~~~
~~~
Pressure ratio
Temperature ratio
Volume ratio
hax/kmin
9,ICFM
S/N,ICFM/rpm
Reduced temp min/max
Reduced pressure min/max
Schultz factors
XmdXmax
YmdYnux
3.25
1.226
2.91 2
1.O281
22734
6.31 5
0.996/1.221
0.30911 .O06
0.509/1 .O56
1.1 50/1.287
565.5
0.681
2.266 x 10'
1.1027
1 .o495
1 .o04
26570
2731O
0.781
31 790
Table C.6.6 compares the test data and the test data converted to specified operating conditions
with the predicted performance at the specified operating conditions. Supplement C.6.B demonThe calculation ofdischargeconditions,
pressure,
strates the calculationsfortheconversion.
temperature, and volume is shown in SupplementC.6.C,
which alsoillustrates the use of an
iterative procedure.
As canbe seen the inlet capacityfortheconverted
test conditions was within 1 percent of
the original specifieddesign point andthehead was within 2 percent.Theconvertedspecific
volume ratio was within the specified 4 percent allowed.
Supplement C.6.A
S T D - A S M E P T C L O - E N G L L777
0757b713 Ob055b5 5 7 8
TABLE C.6.4
TEST DATA
Mass flow
SPd
4,923 - Ibm/min
2,245 - rprn
Inlet
P
Pressure, psia
T
Temperature
Specific volume,
ft3/lbm0.9639
0.9234
V 6
Compressibility
factor
0.944
0.955
Z
Viscosity, centipoise
CL
Specific heat, Btu/lbm-OR
=P
Specific heat ratio
k
Sonic velocity, ft/sec
a
Enthalpy, Btu/lbrn
140.04
h
Entropy,Btu/lbrn-OR
0.2639
0.2731
S
Isentropic
Discharge
67.520
1 O0 167.49
2.871
0.975
0.01 o9
0.2098
1.O98
538.8
122.3
0.2639
67.5
187.4
135.80
TC
590.3 psia
213.8
OF
TABLE C.6.5
DERIVED TEST FUNCTIONS
14143
Pressure ratio
Temperature ratio
ratio
Volume
q ICFM
qhV ICFWrpm
1st stage tip speed
(ft/sec)
Machine Mach no.
Machine
3.49
Reynolds
Rem
no.
exponent
Isentropic
factor
Polytropic
work
exponent
Polytropic
Polytropic head (ft-lbf/lbm)
efficiency
Polytropic
Unit gas
(hp)
power
Pressure Ratio
3.375
1.1 56
2.980
'P
'1
'V
6.3
352.6
0.654
x 106
1.O718
7 1.o01
1.1139
10735.2
0.778
2059
U
Mm
ns
f
n
WP
VP
p.
rp = p&i
67.5
20
= 3.375
TemperatureRatio
rt = T&;
(460 + 187.4)
(460 + 100)
647.4
=-
560
= 1.156
141
0759b70 O b O 5 5 b b 42q D
S T D - A S M E P T C 10-ENGL L977
TABLE C.6.6
DATA SUMMARY
Test Data Converted
To Specified
Test Data
operating
Conditions
N
9
9m
WP
7lP
PS
Pi
660.8
Pd
ti
246.7
b
0.7578 Vi
4.8
0.7578
Vd
rPm
ICFM
2245
141 37
6.297
10736
0.778
2059
20
67.5
115
1O0
187.4
2.871
0.9635
2.98
122.3
140.04
ft-lbf/lbm
hP
psia
psia
OF
OF
ftVbm
6
V/lbm
Vr
164.9
hi
hd
164.9 Btu/lbrn
Btu/lbrn
InletSpecific
Volume
vi
Predicted Performance
At Specified
Operating Conditions
3600
22670
6.297
27690
0.780
32180
200
3600
22734
6.31 5
2731O
0.781
31790
200
650
115
0.2562
2.958
0.2602
2.91 2
210.5
ZiR Ti
-=
209.84
0.975
1545
(
x
(560)
)
= 2.8716
(20)
144pi
144
r, =
V&
2.8716
0.96347
= 2.9805
InletCapacity
Capacity/Speed Ratio
q/N =
141
37
2245
= 6.297
142
TD N
=
U=
a 3 6 . 0 x 2245
720
= 352.6 ft/sec
Mm = U/a =
MachineMachnumber
352.6
538.8
= 0.654
= Ub/pv
352.6 (2.5/12)(1488.2)
0.01 O9 (2.871 6)
n,
Isentropic Exponent
= In(pdpi)/ln
=
f=
(1 38.50
1 .O7212
0.072 1 2
In (67.5/20)
In (2.871 6/0.9234)
1.07212
- 20 (2.871 6)]144
= In (rp)/ln(rv)=
=
1.1138
143
(vM)
- 122.3) 778.1 7
(67.5 (0.9234)
Polytropic Exponent
3.493 x 106
= 1.002
In (3.3750)
In (2.9805)
S T D - A S M E P T1C0 - E N G L
Polytropic Head
0 7 5 9 b 7 0 Ob055bB 2 T 7 9
1997
n
wp = (7)
(PdVd - pivil
1.1138
[67.5
(m)
(0.96347)
- 20
= 10736 ft-lbf/lbm
= 0.7777
pg=
Gas Power
0.778 (33000)
TIP
= 2059 hp
Supplement C.6.B
Calculation:Conversion From TesttoSpecifiedPerformance
Inlet Capacity
= 141 35
(-)3600
2245
= 22670 ICFM
RC, =
0.988
Remp243
--
0.988
144
= 0.0254
(3.49 x 106~243
(2.871
611
144
~~~~
~~
0 7 5 7 b 7 0 Ob055b9 L33 D
S T D - A S M E P T C L O - E N G L L977
Remf
= 1.0
13.67
RBf =
[4.8X lo6 X b]
Remf
2.5
4.8 X
= 0.66 + 0.934
RC, =
RCt
lo6 X -
0.0254
12]
3.493 x 106
= 1.565
0.988
0.988
= 0.01 61 2
(Rem,p)0.243(2.266 x 1 07)0.245
--
(4.8 X lo6 X
Remsp
b)"sp
("'1
4.8 x 1 0 6 x -
= 0.66 + 0.934
RB,
= 1.548
2.266 x lo7
= 1.0
1 .S48 1 .O
= (1 - 0.7777)
-
1.565
- qp,
- VPsp = 0.2199
145
o.o1612
1.0
VPSP
= 0.780
Remcorr= 1.O03
Polytropic Head
(2)
2
Wpsp= Wpt
= 110736
Remcorr
3600
(-)2245
1.O03 = 27690
ft-lbf
Ibm
Power
3600
(-)
= 14137 2245
Pg =
w~sp
SP
(-)0.7578
- 27690(2991
= 29915
5) = 32, 8o hp
0.780
(33000)
Supplement C.6.C
The conversion from test conditions to computed specified conditions involves an iteration to
obtain the discharge pressure from theknown head and discharge enthalpy. Theiteration procedure
and calculation involves assuming a dischargepressureat
the known dischargeenthalpyand
finding the corresponding temperature and specific volume.
The polytropic exponent and polytropic
head is then calculated for the assumed discharge pressure. This polytropic head is then compared
to the actual and, if not the same, then a new discharge pressure is assumed.The new assumed
pressure is evaluationforproperties
at the known dischargeenthalpy,and
a new discharge
volume is evaluated and polytropic exponent are computed. This iteration procedure is continued
until the conditions matchtherequired head.
146
~~
630
650
660
670
680
Pressure, psia
148
690
qp= 0.780
f= 1 .O04
Step 7 - Calculatedischargeenthalpy.
= 165.9 +
Step 2
- Assume
Btu
27690
= 210.5 0.78 (778.1
7)
Ibm
a discharge pressure.
Pout = 660.8 psia
0.2562 ft3/lbm
head.
1 .O04
[660.8
(0.2562)- 200 (0.7578)]
144
fi-l bf
= 27705 Ibm
149
This sample calculation illustrates the computational procedure, at specified operating conditions,
for a multisection compressor havingexternally piped intercoolers.
Consider a two stageaircompressor
equipped with oneintercoolerand an aftercooler.
Section 2
U 1065 D = 12.204
N = 20,000
Section 1
U = 1200 D = 13.751
N = 20,000
-:
After
After
Cooler
Final discharge
discharge
Final
conditions
"
"
Leakage
Condensate
Leakage
Condensate
Power requirements
* Measurement stations
It is desired to calculatethe compressorperformance at thespecifiedoperating
conditions
shown. The compressor has been tested and the test data reduced
to the following dimensionless
form. Thedatawas
collected with pressureandtemperaturebeingmeasured
at the inlet and
outlet of eachsection.The
flow coefficientswerecalculated
based upon test rotor flow rates.
The selectionof testmethodandthemeansofestablishingleakageandcondensate
flow rates
were subject to prior agreementby parties to the test.
The first step in calculating thespecifiedoperating condition point of interest is to establish
the first section performance, starting
with the flow coefficient. Taking the saturation pressure of
water vapor at 560"R to be approximately 0.949 Ibf/in2, with the remaining specified operating
conditions at the inlet, we obtain
and
STDOASME PTCLO-ENGL
L997
0 7 5 9 b 7 0 Ob05575 437 D
= 0.0250 Ibrn
w/lbm
da
The rotor flow rate is the same as the inlet mass flow rate.The flow coefficient is then
Wrotor
Q=
w~RT~
(F)
13.751
= 0.0295
3
With the flow coefficient established the corresponding polytropic efficiency, polytropicwork
1 dimensionlesscurves
coefficient,and total workinput coefficient areread fromthesection
(see Fig.C.7.1).That
is
= 0.83, ,up = 0.599, and
= 0.722 at C$ = 0.0295.To continue
the calculationstheproperties of air at thespecifiedoperatingconditionsmust
be known. For
the purpose of this example we assume that the air-water vapor
mixture maybetreated asan
r),
Section 2
Section 1
0.50
E
I
U,
0,040
Q
The Mm, Rem, k, and v;/vdfor the data are assumed to match the specified operating conditions
within Table 3.2 Limits. The Rem match is assumed sufficiently close
so as to render the Reynolds
number correction negligible.
FIG. C.7.1
152
ideal gas with constant specific heat. k = 1.395 is used as being representative of the end result
which might be obtainedby considering the properties of the constituent
gases over the compression
range. This value will be used for both sections for this example. The average constant pressure
specificheat is
Btu
0.2457 Ibm "R
vP
= 0.83 gives
k
"
n-1
or,
n =
= 0.83
1.395
= 2.931
0.395
1.5178.
= [l
Thedischargepressure
= 0.599 gives
,up
0.599 (1 20012
(2.931) 54.13 (560) 32.174
2.931
= 2.166
is
are
and
The power absorbed in the compressor section is obtained using the total work input coefficient
S1 = 0.722.
The shaft seal which is located downstream of the rotor leaks 0.03 Ibm/sec for these conditions,
is
Weak
= 6.50
- 0.03
= 6.47 Ibm/sec
The intercooler is known to cool the flow to 560"R atthemass flow rate, gas entry state, and
specifiedoperating condition coolanttemperatureand
flow rate.Theairstreamexperiences
a
total pressure loss of 0.8 psi across the intercooler. It must now be determined if and how much
153
condensation occurs in the cooler. Since the cooler exit velocity is assumed, very low stagnation
valuesareused in the analysis. Thesaturationpressure of the vapor at 560"R is approximately
0.949 psia. If the exit air is at 100 percent relative humidity, the humidity ratio is
HRd =
2 zsv
i::;:
Ibrn water
--)
= 0.0196
Ibrn da
p
(31 .O4 - 0.949
where p = 31.84
0.8 = 31.04 psia
Sincethesaturated humidity ratio is less thanthecoolerentry
humidity ratio,condensation
must occur. The difference between the
two is the ratio of condensate to dry air
COndenSate/Wda = HRi
w/lbrn da
W/
= 0.0054 (6.312)
= Ibm/sec
Wda
- condensate
is
The intercooler exit conditions are the second section inlet conditions.
The previous calculation sequence is repeated for the second section, starting
with calculation
of the flow coefficient. Thegas constant changes slightly due to the water vapor removal.
6.436
--
(53.96)
(560)
144 (31.04)
(3
2000 (
ni
= 0.0197
3
S T D * A S M E P T C LO-ENGL 1 7 7 7
0757b70 Ob05578 L 4 b
53.96 1.395
0.2449
Btu
Ibm
O R
yields
pk E =
= 0.81 1.395
0.395
yielding
n
- 2.861, and n
"
n- 1
= 1.537
n
n-1
rp=
1+
fip
gc
(L)
RTi
n- 1
=L1+
0.56
1 0652
-
12.861
= 1.801
ZU2
-60
Wrotor
Gas Power =
gc
3300
6.436
(0.691
1 10652) 60
32.1 74
= 285.1 hp
3 O00
3
Theshaftsealdownstream
ofthe rotor leaks 0.06 Ibm/sec for theseconditions,
flow rateattheaftercoolerentry
is
Wcooler entry
so the mass
Theaftercooler is known to cool the flowto 580"R at this mass flow rate,gasstate,
and
specifiedoperating condition coolanttemperatureand flow rate.Theaftercoolerpressure
drop
is 1 psi. Assuming a saturation pressure of 1.692 psia and following the intercooler condensation
analysisscheme,
155
HR=
-(R,
Rda
psv
1.692
) =A(
53 34
85.76
54.91
- 1.692
= 0.0198
Ibm
W
Ibm da
156
TABLE C.7.1
SUMMARY OF RESULTS
Conditions:
re
Section:
Specified
Inlet mass flow rate
total
Inlet
pressure
total
Inlet
relative
Inlet
humidity
constant,
air Gas
dry
water
constant,
Gas
1 st Section:
Rotational speed
Tip diameter
2nd
Rotational speed
Tip diameter
lbmhec
6.500
14.7
560.
60.
53.34
85.76
20000
psia
"R
%
ft-lbf/lbm-"R
ft-lbf/lbm-OR
rPm
13.751
20000
rPm
in.
in.
12.204
Intermediate CalculationResults:
1st Section
54.13
0.2457
0.0295
0.83
0.599
0.722
0.722
2nd Section
1.51 78
1200
14.7
2.1 66
31.84
560.
729.
381.9
6.5
0.03
6.47
0.0341
0.8
1.537
1065
3 1 .O4
l. 801
55.91
560.
687.9
285.1
6.436
0.06
6.376
O.
1 .o
Overall results:
Final discharge pressure
Final discharge temperature
Total gas power
Delivered mass flow rate
54.9
580.
667
6.38
psia
OR
horsepower
Ibm/sec
157
53.96
0.2449
0.01 97
0.81
0.56
0.691
0.691
ft-lbf/lbm-"R
Btu/lbm-"R
ft/SeC
psia
psia
"R
"R
horsepower
Ibm/sec
lbmkec
lbm/sec
lbmkec
psia
Assumingasinglesectionand
no leakageorsidestreams
and,
Sincetheshaftpower
wspand
EU2
(F)
159
Theterms
ashsp= ash, =
(/%h
- Ppararitic)
5.3, as
33000
and
where
In general the procedure would now be to break down the individual variables in this equation
if the shaft powerweredetermined from
a torque meter, that power would be expressed as the product of measured torque and measured
speed. Similarly,therotor
massrate offlowmight
beexpressed in terms of nozzle pressure
drop,pressure,temperature,and
gas composition. For brevity, in this example Pdt, Q m l r Ut, and
wf are treated as individually measuredelementalquantities.
Steps 2 rhru 4 - Assume thattheelementalerrorsourceshavebeenlisted,theelemental
errorsestimated, and the corresponding bias and precision errors calculated. Many examples of
this procedure may be found in PTC 19.1.
Thisprocessdepends upon theactualinstrumentation systemand data collection techniques
used.Theresultsmay
be expressed as follows. The bias limits and precision indicesrepresent
the combined effects of the independentmeasurementsforeachparameter.
in terms ofindependentmeasurements.Forexample,
Parameter
Absolute
Bias
wt
'
limit
Bwr = 0.01
S+
Bpfif
Qm,
W m ,
ut
Absolute
Precision
Index
sw,= 0.01 W,
W,
= o.o1 Phr
= 0.01
Om,
BU, = 0.01 Ur
B, = 0.2
S, only in this section, is the Absolute Precision Index
MeanNQ
of the
SU,= 0.01 Ut
of the mean = S / a .
PshsP
[-$($)
#Om, = a
ag,
by
u 2
+(%)']J
33000
A bias error is assumed in the mechanical loss conversion equation dueto an assumed unknown
error in the exponent . It is estimated for this example as 0.2.
The precision index for the result is the square root of the sum of the squares of the product
of sensitivity coefficients and averageindependentparameterprecisionindices.Thus
The bias limit for the result is the square root of the sum of the squares of the product of the
sensitivity coefficients and average independentparameterbias limits. Thus
UA00
psh
BpShsp + t9sSpSh
SP
SP
The value r is called the Students r. It is assigned depending upon the degrees of freedom of
the sample, which is usuallyone less thanthenumber
of points averaged.SeePTC
19.1 for
furtherexplanation.Assuming
a large sample, r = 2 maybe used.
Results
In order to allow expression of the results of this example numerically, assume
dpsh,
4QM1 -
-0.00456
33000
= -1.20 (1 .0512(1
- 0.1) (%)
= -1.1907
w t
- 0.1)
(3)
-
= -2.6638
(2)
I
(2)
S T D - A S M E P T C LO-ENGL L777
= 41.418 X
0757b70 Ob05585 3 8 b M
of the result is
+ 1.75 X
+ 7.10 X
+ 3.74 X
+ 1.20 X
10" = 0.0325
Psh,
... (b)
Theuncertaintiesare
Discussion
The UADDand URSS
uncertaintiesmaybeinterpreted
as follows. For UADD= 0.097 P+ the
measurement Pshr f 0.097 PSht will be expected to contain the true value 99 percent of the time.
the uncertaintiesbecome
It is clear that such analysis is of great value in both planning a test and evaluating test results.
163
APPENDIX D
REFERENCES
(ThisAppendixisnot
(D.9)Daugherty,
R. L., and J. B. Franzini. Fluids
Mechanics with Engineering Application.
McGraw Hill Book Co.; 1977.
M.
Huber, M. L., and M. O. McLinden. Thermodynamic Properties of R134a (1,1,1r2-Tetrafluoroethane). July 14-1 7,1992 proceedings,
International
Refrigeration
Conference,
Purdue University, West Lafayette, IN.
(D.12) J.E. Lay. Thermodynamics. Charles E. Merrill Books, Inc., 2nd printing, 1964.
Simon, H., and A. Bulskamper. On the Evaluation of Reynolds Number and Related Surface
Roughness Effects on Centrifugal Compressor
Performance Basedon Systematic Experimental Investigations. ASME paper no. 83 GT-118:
Transactions of the journal of Engineering for
Power, presented March 27, 1983.
Nathoo, N. S., and W. G. Gottenberg. Measuring the Thermal Dynamic Performance of
Multi-Stage Compressors Operating onMixed
Hydrocarbon Gases. December1981proceedings of the 10th Turbomachinery Symposium, Texas A and M.
(D.13) J. M. Schultz.ThePolytropicAnalysis
of
CentrifugalCompressors.Transactions
of
the ASME, Series A. Vol. 84, lournal of Engineering for Power. January 1962, p. 69-82
and April 1962, p. 22.
(0.14)
(D.15) A. Schaffler. Experimental and Analytical lnvestigation of the Effects of Reynolds Number
and Blade Surface Roughness on Multistage
Axial FlowCompressors. Transactions ofthe
ASME, January 1980,Vol. 102, p. 5-1 3,journa/ of Engineering for Power.
Herd, T. C., and E. J.Hipp. Accuracy Expectations for Gas Turbine and Centrifugal Compressor Performance Testing. Paper
ASME 83GT-128.
(D.16)
R. A. Strub. Influence of the Reynolds Number on the Performance of Centrifugal Compressors. Final Reportof theWorking Group
of the Process Compressor Subcommittee of
the International Compressed Air and Allied
MachineryCommittee(ICAAMC)Zurich,
October 1982.
(D.17)Nathoo,
W. S., andW. G. Gottenberg."A
New Look at Performance Analysis
of Centrifugal Compressors Operating
with Mixed Hydrocarbon Gases.', Transactions of the ASME,
October 1983, Vol. 105,p. 920-926,lournal
of Engineering for Power.
(D.18) Skoch, Gary J., and Royce D.Moore. NASA
Technical Memorandum 7007/5 AJAA-877745. AVSCOM Technical Report 87-C-21
"Performanceof two 10 Ib/sec centrifugal
compressors with different blade and shroud
thickness operating over a range
of Reynolds
Numbers." 23rdJointPropulsion Conference
cosponsored by theAIAA, SA,ASME and
ASEE, San Diego, CA, June 29-July 2, 1987.
(D.19) Moore, M. J.,and H. S. Shapiro. Fundamentals ofEngineeringThermodynamics. John
Wiley & Sons, Inc., 1988.
(D.20)
166
APPENDIX E
RATIONALE FOR CALCULATION METHODS
(This Appendixisnot a part of ASME PTC 10-1997.)
E.l
PURPOSE
The purpose of this Appendix is to describe the Code problem model, background theory,
simplifying assumptions.
E.2
and
PROBLEM MODEL
Dependentparameter
a given set of
= F (many independentparameters)
Examples of dependentparametersaredischargepressure,
head, efficiency,etc.Amongthe
independentparameters aregeometry,speed,
flow rate, inlet gas state, gas properties,etc.The
functional relationship f i s unknown. It is to be determined by the test.
Themanyindependentparametersrepresentthespecifiedoperatingconditions.Practicaltest
situations are often suchthatoneormoreoftheseparameters
is precluded from taking on the
desired specified operating condition value. Means must then be sought to determine the effects
of thedepartures. One method is to employ dimensional analysis.
E.2.1 Dimensional Analysis. The theory of dimensionalanalysisand similarity arediscussed in
PTC 19.23, Guidance Manual for Model Testing,andmany fluid dynamics textbooks. In essence
it provides a means to reduce the numberofparameters
in a problemwhich is expressed in
dimensional terms. This is done by appropriate grouping of dimensional
terms in dimensionless
groups.
.. .
n-
n
-f144piVj(F)
t7p
n-1
hd
167
- hi
" -1
with,
In
Vi
Vd
dJ=
Rem= ULIv
Mm= U/a
gas properties
. . . according to
gastypes
k= C&
Z= pv/RT
It is presumed in performing a dimensionalanalysisthat
all ofthevariablesaffecting
the
thermodynamicand fluid dynamicperformanceofthecompressorhavebeen
included. If so,
different sets ofdimensionalvariables
which combine to form identical sets of independent
dimensionlessgroups (4, Mm, Rem,gas
properties) will haveassociated with them identical
valuesfor r], ,
u, and p i .
This basic model is generally accepted to adequately describe the main features of compressor
performance. It has the immediate advantage of reducing the number of parameters which must
beconsidered in developing a test. But of at leastequalimportance
it provides a means of
accounting for unavoidable departures from desired specified operating conditions.
For example,
it may be used to establish an appropriate test speed to compensate for the effect of a test inlet
temperature which differs from the specified operating
condition temperature.
E.2.3 Allowable Departures. It oftenremains impractical to satisfy all theindependent dimensionless parameter requirements. This situation may be addressed by allowing controlleddepartures
in certain independent dimensionless groups. The assumption is that the limits placed upon these
departures render the effects upon the dependent parameters either
negligible or predictable. The
following approach is taken in this Code.
E.2.3.1 Mach Number. Mach numberdeparturesareassumed
to be of increasing relative
importance as the Machine Mach number increases, This is reflected in the allowable departures
shown in Table 3.2 or Figs. 3.2 and 3.3. It i s assumed that negligible effect is associated with
departure within these limits.
E.2.3.2 Cas Properties. Allowable departuresfromtheideal
gas lawsfor both the test and
the specified gasesare given in Table 3.3. When these limits areexceeded the gas must be
treated as real.
E.2.3.3Reynolds
given
STD-ASME P T C 1 0 - E N G L
L997
Tz=
(hh,)rn
773 = (p/px)rn
where x denotesthesidestreamvalue,and
m denotesthereferencemainstreamvalue.
Theapproachtaken in this Code is to require that the ratio of sidestream to reference flow
rates remain within the limits of Table 3.5 or leakagesperpara.
3.3.6. When these limits are
observed it i s assumed that the effects upon the dependent dimensionless groups are negligible.
No specific restriction is placed upon the density or enthalpy ratios. It is assumed that departures
in theseratios will produce negligible effects upon the dependent dimensionless groups. Where
thorough mixing of inlet streams beforethecompression
is doubtful, this assumptionmay not
be valid. In suchcases the parties to the testmay elect by mutual agreement to further restrict
theseratios as well.
Departures in these secondary dimensionless groups do affect results in the dimensional sense.
This is accounted for in the calculation procedure.
E.2.5
qpsp= q p ,Remcorr=
F1
( 4 , -,4 x M,,,,
as follows:
qm
ppSp= pptRemcorr=
F2
( 4 ,-,
4 x M,,,/
4m
as,,,
169
E.3
Appropriateunitsanddimensionalconstants
computations.
E.3.1 Inlet and Exit Conditions. Thestructureofthe problem model is such that it is necessary
to carefully define the inlet and exit conditions which are used in calculating the dimensionless
groups.The exit conditions are thestagnation condition at thedischargemeasurementstation.
The inlet condition is the stagnation state assigned to the flow stream entering the impeller, and
is denoted by the subscript i on thermodynamic properties.
For a simple single inlet flow stream this is the stagnation state at the inlet flange. For multiple
inlet streams it is the stagnationstate computed from the mixing of the individually determined
streams. A standard calculation scheme is given in subpara. E.5.
E.3.2
Pi 2
N(;)
where
wrotor=mass flow rate entering rotor (mass flow ratecompressed)
pi= inlet total density
N= rotor rotational speed
Dis the blade tip diameter of the 1st impeller for centrifugal compressors
Dis the diameter at the leading edge of the 1ststage rotor blade for axial compressors.
Themass flow rateenteringthe
rotor i s determined giving dueconsideration to all section
inlet and outlet flow streams and leakages.
where
E.3.5 Ratio of Flow Rates. The ratio of flow rates is the ratio of flow ratesat two points in the
flow. It i s given by
170
where
W=
E.3.6
where
The Machine Mach number so defined i s not an actual Fluid Mach number. It is nearly directly
proportional to actualFluid Mach numberswhentheCodedeparture
limits areobserved.The
Codedeparture limits shown in Figs. 3.2 and 3.3 for centrifugal and axial compressorsarealso
given in equationform in Table E.l.
E.3.7
i s given by
Rem = Ub/v
For centrifugal compressors, 6 is the exit width ofthefirst
stage impeller in thesectionof
interest.For axial compressors, b is the chord length at the tip of thefirst stage rotor blade in
the section of interest. The viscosity u is taken for inlet (stagnation) conditions. The Code departure
limits shown in Fig. 3.4 for centrifugal compressorsaregiven in equationform in Table E.2.
E.3.8
171
S T D - A S M E P T C LO-ENGL L777
0759b70 Ob05593 4 5 2
TABLE E.l
MACHINE MACH NO. LIMITS
CENTRIFUGAL COMPRESSORS
Specified Mach No.
Upper
Range Limit
Lower
0-0.214
0.215-0.86
0.86 5 Mm,
Limit
-Mmsp
(0.266 MITI,,- 0.271)
-0.042
AXIAL COMPRESSORS
Specified Mach No.
Limit
Upper
5
Range Limit
Lower
0-0.1
(0.266Mmq
6-0.6
0.1
0.6 5 Mmq
-Mm,
- 0.1 9)
-0.03
+ 0.20)
+ 0.20)
TABLE E.2
REYNOLDS NUMBER APPLICATION LIMITS FOR CENTRIFUGAL
COMPRESSORS'
UPPER LIMIT
Remt/Rem, 5 1.O
x = (Rem,,/107)o.3
Application Range
9x
1x
lo4c Rem,<
lo7c Rem,
Equation
1 x lo7
RemJRem,, =
RemJRemp = 100
LOWER LIMIT
Remt/Remv c 1.O,
x = (Rem,,/107)n.3
Application Range
Equation
9x104cRem,,cl x106
1 x l o 6 c Rem,
RemJRem,,
Remt/Rem,
NOTE:
(1) See Fig. 3.3.
172
= (0.olP
= 0.1
where
W= isentropic work per unit mass
Z U 2 = sum of rotor tip speeds
The isentropic work forthepurposesof
this Code is the work done in an isentropicprocess
between the inlet stagnationstateandthedischargestagnationstate.Theisentropic
work per
pound mass foranideal gas is given by
The isentropic work for a real gas mayalso be calculated from the following:
W, =
ns
f 144 pi vi[(:)
n,- 1
n,- 1
"' - I ]
Equation (a) differs from equation (b) by substituting n, for k and introducing f. For a real gas
the isentropicvolumeexponent is notthe same as k. On test, n, canbe calculated from
In P d
n, =
Pi
In
Vd
Vi
Substitutingthis n, for k in equation (a) would produce a small error unless the isentropic
exponent were constant and equal to
n, along the compression path. The polytropic work factor
f compensates for the difference between nSand the actual isentropic exponent.It is computed from
E.3.9
where
[(?)T
- 11
W, =
f 144 pi Y [(;)y
- 13
n-1
where
or
n =
1
Y-m(1 + X )
and
It is assumedherethat
a variation in n affects W, just as thevarying n, affects W,. The
polytropic work factor f i s assumed to have the same value as computed in para. E.3.8.
E.3.10
Efficiencies. Efficiencies are in general defined as the ratio of ideal to actual work required
in a givencompressionprocess.Thestandardideal
work chosenforthisCode
is the work
required in a polytropic compressionprocess occurring betweenthe total pressure at the inlet
reference stagnation state and the total pressure at the discharge stagnation state. The actual work
is taken as thechange in total enthalpybetween thesestates. It represents theactual work in
the process only in the absence of heattransferandsecondary flow effects.Thedischarge gas
state calculated for specified operating conditions using this efficiency definition therefore assumes
the same relative heattransferandsecondary flow effects as those prevailing attest.
The polytropic efficiency is then
WP
- hi)
(hd
v 5
Ws
(hd - h ) I
174
E.3.11 Work Input Coefficient. The work input coefficient is defined in terms of thestagnation
enthalpy rise. It is a dimensionless representation of the actual gas work not including the effects
ofheattransferandsecondary
flow. The work input coefficient is given by
input
E.3.12 Total Work Input Coefficient. Relative differencesin heat transfer and leakage or sidestream
flow rates often will occurbetween test andspecifiedoperatingconditions.
It is assumedthat
these relative differences are sufficiently small so as to produce negligible changes in the polytropic
work coefficient and efficiency. Theleakageandsidestream
flow ratedifferences,however,can
produce relative differences in actualpowerrequirement.The
following model i s presented to
establish a method to account for these effects as they relate to power consumption. The method
is based on relating the total work input to rotor mass flow rate.
The problem modelandnomenclature
are shown in Fig. E.l. Thedimensionless total work
input coefficient is determined as follows.
The first law of thermodynamicsfor a control volume surrounding the rotor (in Fig. E.l) is
so
175
I
Qrotorl
33000
measurement
P
Figure Nomenclature
W ] = mass flow rateat the inlet measurement station.
hi= enthalpy per unit mass at the inlet measurement station.
Wd= mass flow rate at the discharge measurement station.
hd= enthalpy per unit mass at the discharge measurement station.
W/,,= leakagemass flow rate for gas escaping before the rotor,
.e., upstream of the rotor.
hlu= enthalpy per unit mass for escaping gas. When the gas does not escape, but inlet leakage
occurs, hlu is the enthalpy of the gas outside the seal.
W/d= leakage mass flow rate for gas escaping after the rotor,
.e., downstream of the rotor.
hid= enthalpyper unit mass of theescaping gas.For
inward leakage, it is theenthalpy of
the gas outside the seal.
W,,,= mass flow rateforsidestream
flow enteringafterthemeasurementstation
but before
the rotor.
h,,= corresponding enthalpy per unit mass.
Wsd= mass flow rateforsidestream
flowexiting after the rotor but beforethedischarge
measurementstation.
hd= corresponding enthalpy per unit mass.
Qrotw= net mass flow rate through rotor.
h ~ 1 =enthalpy per unit mass at rotor inlet.
h R Z = enthalpy per unit mass at rotor exit.
Qraor= heat loss rate from rotor.
Q= heat loss rate from section.
P= work input rate excluding mechanical loss.
FIG. E.l
176
Rearrangingandnon-dimensionalizing
with CU2
This equation represents the total work input to the gas in dimensionless form. It is called the
total work input coefficient and is giventhesymbol fl, .e.,
a=
33000 PS
CU2
wrotor
gc
This coefficient bears a closerelationship to thework input coefficient, pin, but accounts
additionally fortheenergylostthrough
heat transferandsecondary flow effects. Like the work
input coefficient it i s assumed to be invariant between test and specified operating conditions at
the same flow coefficient. Its purpose is to aid in properlyaccountingforheattransferand
secondary flow effects in power calculations.
Forheatbalance method tests,
Forshaft
where Pparasitic
represents all power in the shaft power measurement which does not represent
work input to the gas in the compressor section of interest, for example, mechanical losses and
power input to othersections.
177
E.4
SPEED SELECTION
or
1
[ ( y ] ] , = [(3]
SP
or
where
wpt= [(fi)
n - 1.
..
fZi R Ti (;)y
- 11
and
n- 1
178
FIG. E.2
E.5
gases from
mixing.
Introducing the continuity of mass equation yields the mixed inlet staticpressure
The inlet stagnation pressure is obtained by adding the dynamic head deduced from the average
Fluid Mach number at c.
The actual form of the equations to be solved depends upon the choice of gas.The following
set amenable to iterative solution may be written for ideal gases.
Pa
TJT,
"
"
PC
Pc/Pa
180
(1
PC
+) :
(9'
Ac
A,pa
"
[l
(T)
M:]
[1 + (y)
M$]
k
k
k-1
as
where
A=
144 pa
181
0.97
0.2
0.4
0.6
0.8
1.0
W b h
0.98
20
40
80
60
ab
1.07
1.06
1.05
/I///NI//I////////////N///I/N//
2?zzhzl
B
- 0.2
Ab is implied
1.O4
1.O3
1.o1
1.o0
0.99
wblwa
FIG. E.3
APPENDIX F
REYNOLDSNUMBERCORRECTION
(ThisAppendixisnot
TheReynoldsnumbercorrectionfor
centrifugal
compressorsrecommended in thisCode hasbeen
changed significantly from the previous issue of PTC
1 O. The changes resulted from new references not
previously available. The old correction for centrifugalcompressorswasadaptedfromworkon
axial
compressors[Ref.(D.811
since no centrifugal compressor data was available. The
correction for axial
compressorsremainsunchangedfromtheprevious
issue of the Code.
The methodof correction, for centrifugal compressors, recommended in this Code is based on the
work done by Weisner[Ref. (D.211but has been
simplifed for ease of application. The data presented
by Weisner suggests thatthe
Machine Reynolds
number at which acompressoroperates
has an
effect not only on the efficiency, but onthe flow
coefficient and work input coefficient as well. The
correctionsare all based on thedeparturefroma
nominal Machine Reynolds number which may vary
fromonemanufacturer to another.Thecorrection
used in this Code, for centrifugal compressors, is
simplified in that the correction is only applied to
theefficiencyandpolytropicwork
coefficient. No
correction is applied to the flow coefficient or the
workinput coefficient. Additionally,thenominal
condition has been standardized to a Machine Reynoldsnumber
= 4.8 x lo6 x b andthesurface
roughness to 0.000125 in.
Another correction method has been documented
by Simon and Bulskamper [Ref. (D.411. This method
is developed by analogy with the turbulent flow in
rough pipes. Semi-empirical correlations are derived
for efficiency, flow coefficient, head coefficient, and
work coefficient. The equations developed include
a correction to the head, work, and flow coefficients.
Similar correction methods have been proposed by
Casey
[Ref.
(D.1411
andStrub[Ref.(D.1611.Test
datasupportingReynoldsnumbercorrections
has
beenpublished by NASA[Ref. (D.l8)].
183
APPENDIX G
REFINEDMETHODS FOR CALCULATING
TOTALCONDITIONS
(This Appendix is not a part of ASME PTC 10-1997.)
Guidelines are presentedfor calculating total pressure and total temperature with accuracies that exceedthosedeterminedbythesimplermethodsof
paras.5.4.3and5.4.4.
The detailsof thermodynamic property evaluations
are not a part ofthis Code. The calculations outlined
below for both ideal gases and real gasesmust be
supplemented by the data andproceduresneeded
to provide the required properties.
Compressible flow and uniform velocity areassumed for themeasurementstation
in question.
Static absolute pressure, pStatic,
and measured absolute temperature, Tmeas, are the only local measurements. The mass flow rate, W, and pipe flow area,
A, are known andthermodynamicpropertiesare
calculable as needed. The remaining key assumption
is the recovery factor, rf, which is defined in para.
5.4.4 in terms oftemperature.
- Compute velocity
V = (w/60)/pA
Step 4
M = d z m E
Step 6
1+
Pstatic
[1 + ( k - 1) I]&
(k - 1) M
Step 7- Comparefrom
step 6 with thevalue
used in previous steps. If agreement is unacceptable,
use
from step 6 andrepeat steps 2 through 7
until thechange i n
i s acceptable (for example,
Ahtatic < O.OSOF).
Step 8 - Calculate total pressure
p =
c,,
21gccp
IDEAL GASES
- rd V 2
G.l
(1
G.2
k in the above
Step 3
- Compute velocity
Step 4
REAL GASES
Theterm"real
gas" normally suggests thatthe
compressibility factor, Z, is not unity andmust be
included in p-v-T calculationsaccording to pv =
ZRT.
The purpose of an equation of state is to provide
a mathematical expression relating pressure, volume,
andtemperature which corresponds as closely as
possible to known or expected p-v-t behavior.
Most equations of state use temperature and specific volume (or density) as independent variables,
from which pressuremay be calculated. That is,
Step 6
p = pkv)
Step 7 1
t(p, h)
186
APPENDIX H
SI UNITS
(This Appendix is not a part of ASME PTC 10-1997.)
U.S. Customary
Symbol
Description
Units
m2
dimensionless
Btu/lbmole "F
Btu/lbm "F
Btu/lbm "F
0.3048
0.3048
1
1
41 83
41 83
m
dimensionless
N-mkgmole K
N-m/kg K
Btu/lbm "F
41
N.m/kg,K
in.
in.
dimensionless
dimensionless
0.0254
0.0254
1
1
m
m
dimensionless
dimensionless
dimensionless
m/s2
ft/sec2
1
0.3048
dimensionless
0.031 08
1-
2324
1
2324
0.04896
N m/kg mole
kg w/kg da
f
g
gc
Dimensional
constant
32.1 74
H
HR
Molal enthalpy
Humidity ratio
Enthalpy
Coefficient of heat transfer
per unit area (for
combined convection
and radiation)
Btu/l bmole
Ibm w/lbm da
Btu/lbm
Btu/hr ft2 "F
Mechanical equivalent of
heat
ft * Ibf
778.1 7 Btu
C
C
CP
cv
D
d
e
Fa
h
h,
J
K
coefficent
Flow
specific
Ratio
of
heats,
SI Units
0.0929
Conversion
X Factor =
ft2
ft/sec
m/s
ft
83
ft Ibm
Ibf + sec2
N m/kg K
a
m kg
N 'S
N.mkg
N m/s m2 K
not used
dimensionless
dimensionless
1
1
dimensionless
dimensionless
dimensionless
dimensionless
C&"
log
Common
logarithm
1O)
(base
187
S T D - A S M E P T C 10-ENGL 1997
Symbol
Description
U.S. Customary
Conversion
Units
X Factor =
dimensionless
Naperian (natural)
logarithm
Ibm/lbmole
Molecular weight
MW
Machine Mach number
dimensionless
Mm
dimensionless
Fluid Mach number
M
dimensionless
Polytropic exponent for a
m
path on the p - T
diagram
m
Mass (Appendix B only)
Ibm
N
Rotative speed
Tm
dimensionless
Polytropic exponent for a
n
path on the p - v
diagram
lbmole
Number of moles
n
(Appendix B only)
Isentropic exponent foran dimensioness
n,
isentropic path on a
p - v diagram
P
Power
hP
psi
Pressure
P
psi
Velocity pressure
Pv
Total mechanical losses
Btu/min
Om
(equivalent)
Casing heat transfer
Btu/min
Or
Btu/min
External
seal
loss
Qsl
equivalent
ft3/m in
Capacity
9
ft3/min
Volume
flow
rate
9
ft
Ibf/lbm "R
R
Gas constant
RA, RB, RC Machine Reynolds number dimensionless
correction constants
Re
Fluid Reynolds number
dimensionless
Rem
Machine Reynolds number dimensionless
RH
Relative humidity
dimensionless
r
Pressure ratio across fluid
dimensionless
meter
Recovery factor
dimensionless
u
Pressure
ratio
dimensionless
'P
dimensionless
Flow rate ratio
rq
Temperature ratio
dimensionless
rt
Specific volume ratio
dimensionless
rv
S
Btu/lbmole "R
Molal entropy
Heat transfer area of
ft2
SC
exposed compressor
casing and adjoining
in
Pipe
Entropy
Temperature
Temperature
BTU/lbm "R
"R
"F
188
SI Units
dimensionless
kg/kgmole
dimensionless
dimensionless
dimensionless
1
1
0.4536
0.01 667
kg
Hz
dimensionless
0.4536
kgmole
dimensionless
0.746
6895
6895
0.01 757
kW
N/m2(Pa)
N/m2(Pa)
kW
0.01757
0.01 757
kW
kW
0.000471 9
0.0004719
5.381
1
m3/s
m3/s
N.m/kg.K
dimensionless
dimensionless
dimensionless
dimensionless
dimensionless
1
1
1
1
1
1
1
1
41 83
0.09294
dimensionless
dimensionless
dimensionless
dimensionless
dimensionless
N mkgmole K
mz
41 83
N.m/kg.K
K
0.5556
K
0.5556
("F + 459.67)
S T D * A S M E P T C LO-ENGL 1 9 9 7
U 7 5 9 b 7 0 Ob05b09
bliT
U.S. Customary
Description
Conversion
X Factor =
Symbol
U
U
V
V
W
W
X
X
Y
Y
a
77
P
P in
PP
PS
Y
$
Y
e
fl
dl
ftlsec
Btullbm
ftlsec
ft3/l bm
ftllbf/l bm
Ibm/min
dimensionless
dimensionless
dimensionless
SI Units
dimensionless
0.3048
2324
0.3048
0.06243
2.989
0.00756
1
1
1
0.3048
1
dimensionless
dimensionless
dimensionless
m
dimensionless
dimensionless
dimensionless
dimensionless
dimensionless
Ibmlft sec
dimensionless
dimensionless
dimensionless
ft2/sec
Ibm/ft3
dimensionless
ft Ibf
in
dimensionless
1
1.488
1
1
1
0.09294
16.02
1.356
0.0254
1
dimensionless
dimensionless
kg/m S
dimensionless
dimensionless
dimensionless
m2/s
kg/m3
dimensionless
N.m
m
dimensionless
dimensionless
dimensionless
ft
189
m/s
N m/kg
m/s
m3/kg
N m/kg
kgls
ASME Services
ASME is committed to developing and delivering technical information. At ASMEs Information Central, we make
every effort to answer your questions and expedite your orders. Our representatives are ready to assist you in the
following areas:
ASME Press
Member Services & Benefits
Public
Information
Other ASME Programs
Codes & Standards
Self-study Courses
InquiriesPayment
Shipping Information
Credit Card Orders
IMechE Publications
SubscriptionsNournalsMagazines
Development
Professional
Meetings & Conferences
Short Courses
Volumes
Symposia
Member Dues Status
Papers
Technical
Publications
There are four options for making inquiries*or placing orders. Simply mail, phone, fax,
or E-mail us and an Information
Central representative will handle your request.
Mail
ASME
Law 22
Drive, Box 2900
Fairfield, New Jersey
07007-2900
Fax-24 hours
973-882-17 17
973-882-5155
E-Mail-24 hours
Infocentral
@asme.org
Information Central staff are not permitted to answer inquiries about the technical content of this code or standard.
Information as to whether or not technical inquiries are issued to this code or standard is shown on the copyright
page. All technical inquiries must be submitted in writing to the staff secretary. Additional procedures for inquiries
may be listed within.