Centrifugal Compressor Polytropic Performance-Improved Rapid Calculation Results-Cubic Polynomial Methods

International Journal of
Turbomachinery
Propulsion and Power
Article
Centrifugal Compressor Polytropic Performance—Improved
Rapid Calculation Results—Cubic Polynomial Methods
Matt Taher 1, * and Fred Evans 2
1 Principal Engineer, LNG Technology Center, Bechtel Energy, Houston, TX 77056, USA
2 Independent Researcher, Canyon Lake, TX 78133, USA; bfred614@gvtc.com
* Correspondence: mtaherim@bechtel.com or mtaher@ASME.org
Abstract: This paper presents a new improved approach to calculation of polytropic performance
of centrifugal compressors. This rapid solution technique is based upon a constant efficiency,
temperature-entropy polytropic path represented by cubic polynomials. New thermodynamic
path slope constraints have been developed that yield highly accurate results while requiring fewer
computing resources and reducing computing elapsed time. Applying this thermodynamically
sound cubic polynomial model would improve accuracy and shorten compressor performance test
duration at a vendor’s shop. A broad range of example case results verify the accuracy and ease
of use of the method. The example cases confirm the cubic polynomial methods result in lower
calculation uncertainty than other methods.
Keywords: polytropic process; compressor performance test; ASME PTC-10; piecewise cubic polyno-
mial approximation; real gas; temperature-entropy; polytropic path

Citation: Taher, M.; Evans, F.

Centrifugal Compressor Polytropic
Performance—Improved Rapid 1. Introduction
Calculation Results—Cubic A highly accurate centrifugal compressor polytropic performance approximation
Polynomial Methods. Int. J.
method has been developed that is easy to employ. This real gas method is based upon
Turbomach. Propuls. Power 2021, 6, 15.
a constant efficiency, temperature-entropy polytropic path. The elegance of this method
https://doi.org/10.3390/ijtpp6020015
is its exceedingly simple way of calculating polytropic efficiency with sufficiently high
precision as required for compressor performance testing, while providing the polytropic
Academic Editors: Francesco Martelli
compression path on T-s and h-s diagrams to a high degree of accuracy. A constant effi-
and Marcello Manna
ciency polytropic path can be modeled as either a single or several sequential piecewise
Received: 31 December 2020
cubic polynomial segments affording solutions that allow for determining thermodynamic
Accepted: 21 May 2021
state variables along a continuous path [1,2]. New analytic terms have been developed for
Published: 28 May 2021 slope and curvature of temperature versus entropy along the constant efficiency polytropic
path [1]. Furthermore, a new screening method has been developed to assist in determining
Publisher’s Note: MDPI stays neutral how many cubic polynomial segments are required to provide sufficient accuracy for a
with regard to jurisdictional claims in given application. Example cases are reported that demonstrate accuracy and compar-
published maps and institutional affil- isons to other polytropic performance calculation methods as well as a review of results
iations. documented by Evans [3]. Cubic polynomial endpoint path methods are demonstrated to
achieve better accuracy than any other endpoint polytropic efficiency calculation method.
Cubic polynomial sequential segment path methods are shown to be superior to other
multi-point numerical methods described in available literature.
Copyright: © 2021 by the authors.
The Taher–Evans Cubic Polynomial method (TE-CP) provides not only overall poly-
Licensee MDPI, Basel, Switzerland.
tropic compression performance results but has the ability to predict fluid state parameters
This article is an open access article at any arbitrary point along the polytropic compression path. The analysis is based upon
distributed under the terms and an inlet flange to discharge flange constant efficiency polytropic path for a single, uncooled,
conditions of the Creative Commons compressor section that may contain multiple impellers. Thermodynamic state parameters
Attribution (CC BY-NC-ND) license used in the calculations are based upon total conditions at inlet and discharge measurement
(https://creativecommons.org/ locations. Typically, these are pressure and temperature measurements. An Equation of
licenses/by-nc-nd/4.0/). State (EOS) provides all other necessary thermodynamic state parameters based upon
Int. J. Turbomach. Propuls. Power 2021, 6, 15. https://doi.org/10.3390/ijtpp6020015 https://www.mdpi.com/journal/ijtpp

Int. J. Turbomach. Propuls. Power 2021, 6, 15 2 of 43
known fluid composition, total pressures and total temperatures. These guiding principles
align the TE-CP calculation method with requirements of the ASME PTC-10 [4].
1.1. Polytropic History

A compressor polytropic performance calculation method was first documented in
the 18600 s by Zeuner as discussed in the 1906 English version of his thermodynamics
book [5]. In the 1960’s, Schultz [6] documented an expanded version of Zeuner’s work
and added a correction factor to acknowledge that fluids were neither perfect nor ideal.
These polytropic analysis methods were included in the ASME PTC 10 version published
in 1965 [7], resulting in the industry accepted vernacular label of “Schultz Methods”. The
1997 version of ASME PTC 10 [4], retained the basic Schultz methods but changed the path
definition model from constant efficiency dictated by [6,8], to constant polytropic exponent,
which is an incorrect definition for the polytropic process of real gas compression. While
useful for some applications, Schultz’s methods have been shown to provide results that
are less accurate than required, especially near a fluid’s critical point and in the dense phase
region [9].
Since Schultz’s methods were first codified by ASME, analysts have published
suggested refinements, improvements and alternatives. Some notable references are,
Kent [10]; Mallen and Saville [9]; Nathoo and Gottenberg [11,12]; Huntington [13,14];
Hunseid, et al. [15]; Oldrich [16]; Sandberg and Colby [17]; Taher [18]; Wettstein [19];
Plano [20]; Evans and Huble [21,22]; Sandberg [23]; Taher [1]. The variations included
the use of several different polytropic path equations as well as a plethora of numerical
integration techniques. Evans and Huble [21] provided concise reviews of several of these
methods while a tutorial by Evans and Huble [22] provided implementation details for
some of them.
It is important for the reader to keep in mind that any calculated polytropic path is an
approximation based upon an assumed model of the actual thermodynamic process rather
than an absolute knowledge of the exact path. However, in much of the currently available
technical literature on the subject, Schultz’s methods have been described as an almost
universal default definition for “polytropic performance”. Evans and Huble [22], pointed out
the error of assuming that Schultz’s formulas represent “a fundamental truth instead of a model
used for convenience”. George Box, a statistician, said “Remember that all models are wrong; the
practical question is how wrong do they have to be to not be useful” [24]. While Schultz’s methods
are acceptable for introducing polytropic concepts in a classroom, their limitations must also
be taught. When using any mathematical model of a constant efficiency polytropic path, an
analyst must understand the limitations of the model and take care not to apply it outside
of its useful boundaries. This paper thoroughly documents the thermodynamic soundness
and superior features of the Taher–Evans Cubic Polynomial methods for calculation of
polytropic efficiency to satisfy the requirements of ASME PTC-10 [4].
In the quest to provide improved calculation methods for the PTC-10 code, the re-
ality of how OEMs operate compressor test stands was taken into account. Real time
performance analyses and results can be provided by existing, dedicated, data acquisitions
systems, which can significantly shorten testing duration and thus, lower costs. This
imposes two competing requirements on polytropic efficiency calculations within the
computing systems; namely: (1) high accuracy and (2) rapid solutions.
Many publications have documented that numerical methods satisfy the first require-
ment. However, numerical solutions typically involve significant resources and computing
elapsed time. The Taher–Evans Cubic Polynomial methods described in this paper have
been shown to uniquely satisfy both requirements stated above as well as provide a
thermodynamically sound documentation of the theory behind the solutions.
1.2. Polytropic Process

Unlike an isentropic process, the polytropic process accounts for degradation of
energy. Converting mechanical work to gas energy in a compressor results in the useful
compression
Int. J. Turbomach. Propuls. Power 2021, REVIEWvdp, which is always less than the gas enthalpy change, dh, by3 of
6, x FOR PEER work, the45
amount of Tds as shown in Equation (1).
work, 𝑣𝑑𝑝, which is always less thandhthe

=gas + Tds change, 𝑑ℎ, by the amount of 𝑇𝑑𝑠 (1)
vdpenthalpy as
shown in Equation (1).
The polytropic compression process is defined as a reversible and non-adiabatic
𝑑ℎ = 𝑣𝑑𝑝 + 𝑇𝑑𝑠 (1)
process, in which the ratio of useful compression work, vdp, to the change of enthalpy, dh,
remainsTheconstant
polytropic compression
along process
the polytropic is defined as path
compression a reversible and non-adiabatic
as defined by Equation process,
(2) and
in which the ratio of
illustrated in Figure 1. useful compression work, 𝑣𝑑𝑝, to the change of enthalpy, 𝑑ℎ, remains
constant along the polytropic compression path as defined by Equation (2) and illustrated in
vdp
Figure 1. ηp = (2)
dh
Figure1.1.Polytropic
Figure Polytropiccompression
compression path
path on
on the
the h-s
h-sDiagram.
Diagram.
Since the polytropic process is defined as𝑣𝑑𝑝 a reversible process, the degraded energy
𝜂 = R2
can be represented by heat, which is equal to𝑑ℎ1 Tds. The polytropic work, w p , as the
(2)
measure of useful
Since compression
the polytropic processwork, is theas
is defined difference between
a reversible the
process, total
the enthalpy
degraded change
energy can
and the area under the polytropic path on the T-s diagram between inlet and discharge, as
be represented by heat, which is equal to 𝑇𝑑𝑠. The polytropic work, 𝑤 , as the measure of
illustrated in Figure 2. Substituting Equation (1) into Equation (2) yields a convenient form
useful compression work, is the difference between the total enthalpy change and the area
for compressor polytropic efficiency shown in Equation (3).
under the polytropic path on the T-s diagram between inlet and discharge, as illustrated in
Figure 2. Substituting Equation (1) into EquationTds (2) yields a convenient form for compressor
polytropic efficiency shown in Equation = 1−
η p (3). (3)
dh
The compressor polytropic efficiency 𝑇𝑑𝑠 deviates from unity (η p < 1) for an
𝜂 = always
1− (3)
uncooled compressor section (Tds > 0). 𝑑ℎ
The shape of the polytropic path on the T-s diagram, which illustrates the gas behavior
The compressor polytropic efficiency always deviates from unity (𝜂 < 1) for an un-
as temperature and entropy increase during the compression process, can significantly
cooled compressor section (𝑇𝑑𝑠 > 0). R
influence the value of degraded energy, Tds. Additionally, the use of different equations
The shape of the polytropic path on the T-s diagram, which illustrates the gas behavior
of state to predict thermodynamic conditions along the polytropic compression path can
as temperature and entropy increase during the compression process, can significantly influ-
greatly influence
ence the value of the results,energy,
degraded as illustrated by example cases
𝑇𝑑𝑠. Additionally, inofAppendix
the use D.
different equations of state
to predict thermodynamic conditions along the polytropic compression path can greatly in-
fluence the results, as illustrated by example cases in Appendix D.
Int. J. Turbomach. Propuls. Power 2021, 6, x FOR PEER REVIEW 4 of 45
Figure2.2.Area
Figure Areaunder
underthe
thepolytropic
polytropicpath
pathon
onthe
theT-s
T-sdiagram
diagramrepresents
representsdegradation
degradationof ofenergy
energyin
inaa
polytropic
polytropiccompression
compressionprocess. Additionally,
process. seesee
Additionally, Figure 10 10
Figure forfor
different possible
different curve
possible shapes
curve of the
shapes of the
polytropicpath
polytropic pathbetween
betweenendpoints.
endpoints.
2.2.Polytropic
PolytropicPath:
Path:Polynomial
PolynomialApproximation
ApproximationMethodsMethods
The constant efficiency polytropic path
The constant efficiency polytropic path on the on the T-sT-s
diagram
diagram is a is
continuous
a continuousreal-valued
real-valued
function
functionof ofentropy
entropydefined
definedon onaabounded interval [s𝑠1 , ,s𝑠2 ]. . This
bounded interval This actual
actual polytropic
polytropic pathpath can
can be approximated to a high degree of accuracy using a single polynomial
be approximated to a high degree of accuracy using a single polynomial or a series of sequen- or a series
oftial
sequential
piecewise piecewise
segmentsegment polynomials.
polynomials. The concept
The concept is to accurately
is to accurately estimate estimate
the actualthe un-
actual unknown temperature function, T(s), by an approximating
known temperature function, T(s), by an approximating function that is simple but function that is simple
R 2 accurate
but accurate enough (i.e., within acceptably small error
enough (i.e., within acceptably small error tolerances) to calculate tolerances) to 𝑇𝑑𝑠, which 1isTds,
calculate the de-
which is the degraded part of energy transfer in the polytropic compression process.
graded part of energy transfer in the polytropic compression process. Taher [1], documented
Taher [1], documented the mathematics of applying sequential piecewise cubic polynomials
the mathematics of applying sequential piecewise cubic polynomials to approximate a con-
to approximate a constant efficiency temperature-entropy polytropic path.
stant efficiency temperature-entropy polytropic path.
The temperature-entropy relationship along the polytropic path can be approximated
The temperature-entropy relationship along the polytropic path can be approximated
using a single polynomial for the entire range from inlet to discharge. This is referred to as
using a single polynomial for the entire range from inlet to discharge. This is referred to as
an “Endpoint” method since the thermodynamic state variables at both compressor inlet
an “Endpoint” method since the thermodynamic state variables at both compressor inlet and
and discharge are required inputs and are known from testing. For some applications,
discharge are required inputs and are known from testing. For some applications, when a
when a single polynomial cannot provide sufficient accuracy, a series of sequential steps
single
or polynomial
segments cannot provide
of polynomials can besufficient
employed. accuracy, a series ofsections
The following sequential willsteps or segments
discuss and
of polynomials can be employed. The following sections will discuss and compare R2 first- and
compare first- and third-degree polynomial approximants as used to calculate 1 Tds. First
third-degree
degree polynomial
polynomials approximants
are described as used
as using “steps” while third𝑇𝑑𝑠
to calculate . Firstpolynomials
degree degree polynomials
are
described as using “segments”. This distinction is due to the discontinuity of path slope“seg-
are described as using “steps” while third degree polynomials are described as using
ments.”for
change This distinction
linear steps whereasis duesequential
to the discontinuity
piecewise cubic of pathsegmentsslopepossess
changea for linear steps
continuous
whereas sequential piecewise cubic segments possess a continuous
slope change. It will be shown that this continuity of slope at the knots gives rise slope change. It will
to be
shown that this continuity of slope at the knots gives rise
increased accuracy while using fewer segments to achieve the desired accuracy. to increased accuracy while using
fewer segments to achieve the desired accuracy.
2.1. Temperature-Entropy Polytropic Path Approximation: Linear Polynomial Endpoint Method
2.1. ATemperature-Entropy
linear straight-line Polytropic
approximationPath Approximation:
between inlet and Linear Polynomial
discharge Endpoint
conditions hasMethod
been
used by Kent [10]
A linear and Sandberg
straight-line [23] to approximate
approximation the actual
between inlet polytopicconditions
and discharge path on thehasT-sbeen
diagram.
used by Neither
Kent [10] Kentandnor Sandberg
Sandberg explicitly
[23] mentioned
to approximate theemploying a linear path
actual polytopic approxima-
on the T-s
tion for theNeither
diagram. polytropicKentpath
noron the T-s diagram,
Sandberg explicitlybut in fact a employing
mentioned straight-linea approximation of
linear approximation
temperature, T, as apath
for the polytropic function
on theofT-s
entropy, s, is but
diagram, usedinin their
fact linear approximation
a straight-line method.
approximation of tem-
Since pressure
perature, T, asand temperature
a function are known
of entropy, at both
s, is used endpoints,
in their all other thermodynamic
linear approximation method. Since
properties at inlet and discharge can be determined by an equation
pressure and temperature are known at both endpoints, all other thermodynamic of state. properties
Equation
at inlet (4) represents
and discharge can betemperature
determinedas bya an
function of entropy
equation of state. connecting compressor
section inlet and(4)
Equation discharge points
represents on a T-s diagram.
temperature The form
as a function of the function
of entropy connectingis acompressor
classic
straight
section line.
inlet and discharge points on a T-s diagram. The form of the function is a classic
(1)
straight line. 1 T (s) = ms + b (4)
Int. J. Turbomach. Propuls. Power 2021, 6, x FOR

ThePEER REVIEW
constraints applied to Equation (4) are shown in Equations (5) and (6).5 of 45
T (s1 ) = T1 (5)
( )
𝑇(𝑠)
T (=s 𝑚𝑠
) =+T𝑏 (4)
(6)
2 2
The constraints applied to Equation (4) are shown in Equations (5) and (6).
Equations (7) and (8) show the analytically derived values for the slope and intercept
) = 𝑇 are constant and only dependent
of the straight-line approximant. These𝑇(𝑠values (5)upon the
thermodynamic parameters at the compression endpoints.
𝑇(𝑠 ) = 𝑇 (6)
T2 − T1 ∆T
Equations (7) and (8) show the analytically
m= derived
= values for the slope and intercept of (7)
s − s
2 constant
the straight-line approximant. These values are 1 ∆s
and only dependent upon the ther-
modynamic parameters at the compression endpoints.
s2 T1 − s1 T2
b=𝑇 − 𝑇s −∆𝑇 (8)
𝑚= 2 = s1 (7)
𝑠 −𝑠 ∆𝑠
(1)
Substituting Equations (7) and (8) into Equation (4) yields Equation (9), 1 T (s), which
is a first-degree polynomial that approximates 𝑠 𝑇 − 𝑠 the
𝑇 actual polytropic path T (s).
𝑏= (8)
𝑠 −𝑠
− T1 2 T1 − Equation

Substituting Equations ((7)
1) and (8) T2Equation
∼ into (4)syields s1 T2 ( )
(9), 𝑇(𝑠), which (9)
1 T (s) = s+
s2 − the
is a first-degree polynomial that approximates s2 − s1 path 𝑇(𝑠).
s1 actual polytropic
( ) 𝑇 −𝑇
The deviation of the approximated 𝑠 𝑇 − 𝑠from
temperature 𝑇 the actual temperature value
𝑇(𝑠) ≅ 𝑠+ (9)
𝑠 −𝑠 𝑠 −𝑠 (1)
using a linear approximant is represented by the error function 1 ErT ( S ) in Equation (10).
The deviation of the approximated temperature from the actual temperature value using
a linear approximant is represented (1) (1) ( )𝐸𝑟𝑇(𝑠) in Equation (10).
T (by
s) the
= error function
1 T ( s ) + 1 ErT ( s ) (10)
( ) ( )
𝑇(𝑠) = 𝑇(𝑠) + 𝐸𝑟𝑇(𝑠) (10)
As shown in Figure 3, temperature deviation varies along the polytropic path in the
As shown
interval [s1 , s2 ]infrom
Figure
the3,compressor
temperaturesection
deviation
inletvaries along the The
to discharge. polytropic
actual path in the path
polytropic
interval 𝑠 , 𝑠 from the compressor section inlet to discharge. The actual polytropic path
may largely deviate from a straight line connecting the endpoints.
may largely deviate from a straight line connecting the endpoints.
Figure 3.
Figure 3. Temperature
Temperaturedeviation
deviationof of
thethe
linear approximant
linear fromfrom
approximant the actual polytropic
the actual path aspath
polytropic approx-
as approx-
imated with 10-segment piecewise cubic polynomials is shown. The highest temperature deviation is
imated with 10-segment piecewise cubic polynomials is shown. The highest temperature deviation is
4.276 K (1.0734% relative deviation), in the middle of the interval range (s = 3.409 kJ/kg-K). The error
4.276 K (1.0734%
in estimating relative deviation),
the degraded energy is in𝐸𝑟the
( )
𝑇𝐷𝑆middle of the
= 0.1054 interval
kJ/kg (0.698% range (s =deviation),
relative 3.409 kJ/kg-K).
which isThe error
(1) ( ) ( )
in estimatingbythe
represented thedegraded energy
area bounded ErTDS
is 1the
between = 0.1054
two curves 𝑇(𝑠) and
of kJ/kg (0.698%𝑇(𝑠).
relative deviation), which is
(1) (3)
represented by the area bounded between the two curves of 1 T (s) and 10 T ( s ).
(1) R2
The term 1 ErTDS in Equation (11) shows the error in calculating 1 Tds using a
linear approximation for the temperature-entropy relationship between endpoints of the
polytropic compression path.
Z 2
( T1 + T2 ) (1)
Tds = (S2 − S1 ) + 1 ErTDS (11)
1 2
(1)
The 1 ErTDS error can only be estimated by comparing the result against a more
accurate approximation for the polytropic path. In this paper, a 10-segment cubic polyno-
(3)
mial method (i.e., using 10 T (s) approximants) is used as the basis of evaluation. (See
Appendix B for confirmation of the accuracy of the multi-segment cubic polynomial
method as applied with 10 segments).
Employing a straight line as the approximant to model the actual polytropic path on
the T-s diagram overly simplifies the actual compression path resulting in reduced accuracy
R2
of the calculation for the 1 Tds and thus the resulting polytopic efficiency rendering the
polytropic work values less accurate. For example, as illustrated in Figure 3, a linear
approximation for the polytropic path on T-s diagram overestimates the degraded energy
(1)
(i.e., 1 ErTDS = 0.1054 kJ/kg) by the bounded area between the linear approximant and
the actual path as estimated with the cubic polynomial method. A linear approximant is
inherently limited to a fixed constant slope along the entire polytropic path. This is a major
R2
factor in estimating the 1 Tds with large errors using a linear endpoint approximant.
2.2. Temperature-Entropy Polytropic Path Approximation: Cubic Polynomial Endpoint Method

The desire to improve the accuracy of approximating the actual polytropic path on
the T-s diagram has led to employing the approximant function as a cubic polynomial.
Polynomials of the first degree and second degree are limited since they cannot account for
path slope changes and concavity of the actual polytropic path respectively. A polynomial
of the third degree is the simplest form of a polynomial approximant, which allows
the model to account for these two important features. As explained later in this paper
(see Section 4), the slope and concavity of the polytropic path on the T-s diagram reveal
significant thermodynamic insight about the behavior of fluids during the compression
process and must be carefully studied.
Using Equation (12), the actual temperature that increases with entropy, T (s), along
the polytropic path of a compressor section can be approximated using the third-degree
(3)
polynomial 1 T ( s ).
(3)
1 T (s) = As3 + Bs2 + Cs + D (12)
Since thermodynamic conditions at endpoints are known, the constraints applied to
Equation (12) are shown in Equations (13)–(16).
T (s1 ) = T1 (13)
T (s2 ) = T2 (14)
T1 1 + η p X1

dT
= E1 = (15)
ds η,1 C p1 1 − ηp
1 + η p X2

dT T
= E2 = 2 (16)
ds η,2 C p2 1 − ηp
Relationships (15) and (16) determine values of the slope of the temperature-entropy
curve at endpoints of the polytropic path. For details on how these equations are derived
from thermodynamic relationships, see Appendix A.
Int. J. Turbomach. Propuls. Power 2021, 6, x FOR PEER REVIEW (3) 7 of 45

Coefficients A, B, C and D of the third order polynomial 1 T (s) in Equation (12)
are “analytically” derived from the known conditions above. The results are shown in
Equations (17)–(20). polynomial ( )𝑇(𝑠)
Coefficients A, B, C and D of the third order in Equation (12) are
2 E1 + E2 T2 − T1
“analytically” derived from theAknown
= conditions above. The − results are shown in Equations (17)
(17)–(20). ( s2 − s1 )2 2 s2 − s1
2
𝐸 + 𝐸 𝑇 − 𝑇
𝐴 =B = 1 E2 − E1 −− 3(s1 + s2 ) A (17) (18)
(𝑠 − 𝑠 )
2 s2 − 2s1 𝑠 −𝑠
2
E 1+𝐸E2− 𝐸 3− 3(𝑠 + 𝑠 )
C =𝐵 =12 𝑠 − − 𝑠 2 A s22 2+ s21 𝐴 − B(s1 + s2 ) (18) (19)
2
( T +𝐸T2+) 𝐸 A3 3 B C
D = 𝐶1 = − − 𝐴(𝑠 s1 ++s𝑠32 ) −− 𝐵(𝑠 s+
2
+
𝑠 )s 2
2 − (s + s2 ) (19) (20)
2 2 22 2 1 2 1
Coefficient A accounts
(𝑇 + 𝑇 )for𝐴the deviation 𝐵
of the average
𝐶
value of the slope at the end-
𝐷= − (𝑠 + 𝑠 ) − (𝑠 + 𝑠 ) − (𝑠 + 𝑠 )

(20)
points, E1 +2 E2 , from that
2 of the 2 straight line,
2 which connects
2 the endpoints, ( ∆T T2 −T1
∆s = S2 −S1 ).
Coefficient B accounts
Coefficient for for
A accounts thethe
change of slope
deviation of thealong thevalue
average polytropic path.atWhen
of the slope the polytopic
the end-
∆
path on the T-s, from
points, diagram
that ofapproaches
the straight to a straight
line, line that
which connects theconnects
endpoints,the
( endpoints,
= ). Co-(e.g., in the
∆
case of compressing
efficient B accounts forantheideal gas)ofcoefficients
change slope alongA and
the B approach
polytropic path.zero
When andthecoefficients
polytopic C and D
represent
path on the T-s diagram approaches to a straight line that connects the endpoints, (e.g., in the(7) and (8).
the slope and the intercept of the straight line as shown in relationships
case of
Thecompressing
deviationan of ideal gas) coefficients
the temperature A and B approach
approximated zero
using and coefficients
a cubic polynomial C and
approximant
D represent the slope and the intercept of the straight line as shown in relationships
(3) (7) and
from the actual temperature value is shown by the error function 1 ErT (s) in Equation (21).
(8).
The temperature
The deviation deviation varies along
of the temperature the polytropic
approximated pathpolynomial
using a cubic in the interval [s1 , s2 ] from the
approximant
compressor inlet
from the actual to discharge.
temperature
( )
value is shown by the error function 𝐸𝑟𝑇(𝑠) in Equation (21).
The temperature deviation varies along the polytropic path in the interval 𝑠 , 𝑠 from the
(3) (3)
compressor inlet to discharge. T (s) = 1 T (s) + 1 ErT (s) (21)
( ) ( )
𝑇(𝑠) = 𝑇(𝑠) + 𝐸𝑟𝑇(𝑠) (21)
As illustrated in Figure 4, the deviation of approximated temperature from the actual
As illustrated
temperature showsin Figure 4, the improvement
significant deviation of approximated temperature
for the cubic fromapproximation
polynomial the actual over
temperature shows significant improvement for (the
3) cubic polynomial
(1) approximation over
the linear version, which is expected since
( ) 1 ErT ) s)
( ( ErT (s).
the linear version, which is expected since 𝐸𝑟𝑇(𝑠) ≪ 𝐸𝑟𝑇(𝑠). 1
Figure 4. Temperature deviation of the cubic polynomial endpoint method from the actual polytropic
path as approximated with 10-segment piecewise cubic polynomial approximants is shown. The
highest temperature deviation is 1.185 K (0.289% relative deviation), at s = 3.414 kJ/kg-K. The error
(1)
in estimating the degraded energy is 1 ErTDS = 0.0237 kJ/kg (0.157% relative deviation), which is
(3) (3)
represented by the area bounded between the two curves of 1 T (s) and 10 T ( s ).
Figure 4. Temperature deviation of the cubic polynomial endpoint method from the actual polytropic
path as approximated with 10-segment piecewise cubic polynomial approximants is shown. The high-
Int. J. Turbomach. Propuls. Power 2021, est temperature deviation is 1.185 K (0.289% relative deviation), at s = 3.414 kJ/kg·K. The error8 in
6, 15 esti-
of 43
( )
mating the degraded energy is 𝐸𝑟𝑇𝐷𝑆 = 0.0237 kJ/kg (0.157% relative deviation), which is repre-
( ) ( )
sented by the area bounded between the two curves of 𝑇(𝑠) 𝑎𝑛𝑑 𝑇(𝑠).
(1)
( ) s ), and cubic poly-
Figure
Figure 55 compares
compares the
the temperature
temperature deviation of the
deviation of linear, 1 T𝑇((𝑠),
the linear, and cubic polyno-
( )(3)
mial, 𝑇1(𝑠),
nomial, T (sendpoint
), endpoint approximants
approximants from
from thethe highly
highly accurate
accurate polytropic
polytropic path
path as as ap-
approxi-
( ) (3)
mated with with
proximated a 10-segment piecewise
a 10-segment cubic
piecewise cubic polynomials, 𝑇10(𝑠),
polynomials, T (sfor thethe
), for high
highpressure
pressurepro-
pane casecase
propane withwith
the endpoint conditions
the endpoint shown
conditions on the
shown onfigure. The cubic
the figure. The polynomial endpoint
cubic polynomial
approximant
endpoint provides provides
approximant a maximum temperature
a maximum deviation
temperature that is less
deviation thatthan a third
is less than aofthird
that of
the linear endpoint approximant.
of that of the linear endpoint approximant.
Figure 5. Distribution
Distributionofoftemperature
temperaturedeviation
deviationshows
showsthethe
maximum error
maximum of cubic
error polynomial
of cubic end-
polynomial
point approximant is less than one third of the linear endpoint approximation.
endpoint approximant is less than one third of the linear endpoint approximation.
The term ( 1)ErTDS

(3) R2
𝐸𝑟𝑇𝐷𝑆 in Equation (22) shows
shows the
theerror
errorin
incalculating
calculating 1 𝑇𝑑𝑠
Tds using
using the
the cu-
cubic polynomial
bic polynomial endpoint
endpoint approximation.
approximation.
Z 2 (𝑇 ++T𝑇2 ))
(T ((𝐸
E2 −−E𝐸1 )) ( )
𝑇𝑑𝑠== 1
Tds ((𝑠
s2 −−s𝑠1 ))−− ((𝑠
s2 − + (31) ErTDS
− s𝑠1 ))2 + 𝐸𝑟𝑇𝐷𝑆 (22)
(22)
1 22 12
12
The improved
improved accuracy
accuracy of calculating 1 𝑇𝑑𝑠
of calculating using the
the cubic
cubic polynomial
polynomial endpoint
endpointap-
R2
The Tds using
( ) ( )
proximant, (𝑇3(𝑠),
) as compared to the linear approximant, 𝑇 (1(𝑠),
) can be calculated using
approximant, 1 T ( s ), as compared to the linear approximant, 1 T ( s ), can be calculated
Equation (23).
using Equation (23).
( ) ( ) (𝐸 − 𝐸 )
(1) 𝐸𝑟𝑇𝐷𝑆 −(3) 𝐸𝑟 𝑇𝐷𝑆 = − ( E2 − E1 ) (𝑠 −𝑠 ) (23)
1 ErTDS − 1 ErTDS = −
12 − s1 )2
( s2 (23)
12
Equation (23) shows that the accuracy of the cubic polynomial approximant as com-
Equation
pared to a linear(23) shows that the
approximant accuracy
increases for of the cubic polynomial
compression applicationsapproximant
with a largeas com-
difference
pared to a linear approximant increases for compression applications with a large difference
between the slopes at the endpoints. As explained in the next section of this paper, this ana-
between the slopesconveniently
lytical relationship at the endpoints. As explained
provides a means toin the next
evaluate section of
accuracy of the
thislinear
paper, this
endpoint
analytical
method. relationship conveniently provides a means to evaluate accuracy of the linear
endpoint method.
3. Polytropic Efficiency Calculation
3. Polytropic Efficiency Calculation
Since the late 19th century, various calculation methods have been proposed to esti-
mate polytropic efficiency. The polytropic efficiency can only be estimated. The accuracy of
the estimation depends upon the calculation method and the accuracy of thermodynamic
properties used in the calculation. As shown in this paper, different equations of state for a
fixed method can result in variation of polytropic efficiency values, which may be even
larger than the effect of measurement uncertainties. It is important to ensure that the same
equation of state as used for calculating expected performance is used when calculating
polytropic efficiency based on equipment test results.
Regardless of which calculation method is used, usually five significant digits can
well serve the purpose for any compressor performance evaluation. Accuracy of measured
inlet and discharge conditions as well as of the equation of state and test data used to
obtain thermodynamic properties for the calculation method may impose to reduce the
number of significant digits. The identification of significant digits is only possible through
knowledge of the circumstances [25]. However, in this paper, the intended precision
of polytropic efficiency is chosen as five significant digits in order to evaluate different
polytropic efficiency calculation methods.
A novel approach is used in this paper to differentiate between the actual polytropic
efficiency (for a fixed calculation method and EOS) from the calculated polytropic efficiency
by using an error function, Erη. This is believed to help clarify that regardless of the
calculation method, the actual polytropic efficiency, η p , can only be estimated with some
acceptable error, Erη, as compared to a more accurate reference method. In this paper, the
Taher–Evans Cubic Polynomial 10-segment method is used as the basis to compare the
accuracy of other methods.
3.1. Polytropic Efficiency Calculation: Linear Polynomial Endpoint Method

In the simplest case, the polytropic path can be approximated using a linear approxi-
mant, which connects the endpoints (i.e., one step) as follows:
(1)
1 T (s) = ms + b (24)
R2
The degraded energy in the polytropic compression process, 1 Tds, is approximated
(1)
using the linear endpoint approximant 1 T ( s ):
Z 2 Z 2
(1) (1)
Tds = 1 T (s) ds + 1 ErTDS (25)
1 1
(1) R 2 (1)
By replacing 1 T (s) from (24) and analytically performing the integral 1 1 T ( s ) ds using
the relationship (9), the relationship (26) emerges:
Z 2
( T1 + T2 ) (1)
Tds = (s2 − s1 ) + 1 ErTDS (26)
1 2
(1)
Using the relationships (2) and (26), the polytopic efficiency, 1 η, as approximated
with the linear endpoint approximation method is developed as follows:
R 2 (1)
1 T ( s ) ds s2 − s1

(1) 1 T1 + T2
1η = 1− R2 = 1− (27)
h2 − h1 2
1 dh
(1)
The deviation of the approximated endpoint polytropic efficiency, 1 η, from the actual
(1)
polytropic efficiency, η p , using the linear endpoint approximant, 1 T ( s ), is represented by
(1)
the error function 1 Erη in Equation (28).
(1) (1)
ηp = 1η + 1 Erη (28)
A similar relationship to (27) for calculating the polytropic efficiency was suggested
by Stepanoff [26] in 1955.
3.2. Polytropic Efficiency Calculation: Taher–Evans Cubic Polynomial Endpoint Method

In a more accurate case, the polytropic path can be approximated using a cubic
polynomial, which connects the endpoints (i.e., one segment) as follows:
(3)
1 T (s) = As3 + Bs2 + Cs + D (29)
R2
The degraded energy in the polytropic compression process, 1 Tds, can be approxi-
(3)
mated using the cubic polynomial endpoint approximant, 1 T ( s ):
Z 2 Z 2
(3) (3)
Tds = 1 T (s) ds + 1 ErTDS (30)
1 1
(3) R 2 (3)
By replacing 1 T (s) from (29) and analytically performing the integral, 1 1 T (s) ds,
using relationships (17)–(20), the relationship (31) emerges:
Z 2
( T1 + T2 ) ( E − E1 ) (3)
Tds = ( s2 − s1 ) − 2 (s2 − s1 )2 + 1 ErTDS (31)
1 2 12
(3)
Using the relationships (2) and (31), the polytopic efficiency, 1 η, as approximated
with the cubic polynomial endpoint approximation method is developed as follows:
R 2 (3)
1 T (s) ds s2 − s1 E2 − E1

(3) 1 T1 + T2
1η = 1− R2 = 1− − ( s2 − s1 ) (32)
h2 − h1 2 12
1 dh
Equation (32) can be solved recursively by assuming an initial value for the efficiency
to calculate endpoint slopes E1 and E2 , calculating the resulting efficiency from (32), and
minimizing the difference between assumed and calculated efficiency.
(3)
The deviation of the approximated endpoint polytropic efficiency, 1 η, from the actual
polytropic efficiency, η p , using Taher–Evans cubic polynomial endpoint approximant
(3) (3)
1 T ( s ), is represented by the error function 1 Erη in Equation (33).
(3) (3)
ηp = 1η + 1 Erη (33)
By comparing relationships (27) and (32) the efficiency deviation between the linear
and cubic polynomial endpoint methods can be developed as follows:
1 E2 − E1

(3) (1)
1 Erη − 1 Erη = − 12 h − h ( s2 − s1 )2 (34)
2 1
As shown by the relationship (34), the deviation linearly increases with the change of
slopes, ( E2 − E1 ), at endpoints. This analytical relationship conveniently provides a means
to evaluate accuracy of the linear methods as applied to endpoints.
(3)
Once the polytropic efficiency, 1 η is determined, coefficients A, B, C, and D are
calculated, and the polytopic path model on T-s and h-s diagrams can be graphically
represented using relationships (29) and (51).
Figure 6 compares the efficiency deviation of the linear endpoint and cubic polynomial
(3)
endpoint methods from the highly accurate 10-segment cubic polytropic efficiency, 10 η,
using different equations of state for the high pressure propane case with the condi-
tions shown on the figure. Evidently, the linear endpoint method involves significantly
larger deviations for the polytropic efficiency as compared with the cubic polynomial
endpoint method.
Int. J. Turbomach. Propuls. Power 2021, 6, x FOR PEER

Using REVIEW examples in Appendix B, it is shown that Taher–Evans cubic polyno-
illustrative 11 of 45
mial endpoint method provides the highest accuracy among all other endpoint methods.
( )
(1) 𝜂 , and cubic polynomial, ( )
Figure6.6.Percentage
Percentagedeviation
deviationofofpolytropic
polytropicefficiency
efficiencyofof the
the linear,
linear,
(3) 𝜂,
Figure 1 η, and cubic polynomial, 1 η, ( )
endpoint methods from the highly accurate 10-segment cubic polynomials polytropic efficiency,
endpoint methods from the highly accurate 10-segment cubic polynomials polytropic efficiency, 10 η, ,
(3) 𝜂
as calculated using different equation of states.
as calculated using different equation of states.
3.3.Polytropic
3.3. PolytropicEfficiency
EfficiencyCalculation:
Calculation:Linear
LinearMultistep
MultistepMethod
Method
( )
approximant, 1𝑇T(𝑠),
Asillustrated
illustrated inin Figures
Figures 33 and
and 6, (1)
As 6, the
the linear
linear approximant, (s),when
whenapplied
appliedtotoend-
points may introduce large deviations from the actual T-s polytropic
endpoints may introduce large deviations from the actual T-s polytropic path resulting in path resulting in re-
( )(1)
duced accuracy
reduced accuracyfor forthe calculated
the calculated polytropic
polytropic efficiency, 𝜂1. η.
efficiency, ByBy
applying thethe
applying composite
compos- trap-
ezoidal rule, the integral 𝑇𝑑𝑠 Ris2approximated by employing piecewise
ite trapezoidal rule, the integral 1 Tds is approximated by employing piecewise linear linear polynomial
approximants
polynomial that cover the
approximants thatentire
coverpath in multiple
the entire path insteps with smaller
multiple subintervals
steps with 𝑠 ,𝑠
smaller subin- .
In the linear
tervals [si , si+1 ]. multi-step method, temperature deviation along the polytropic path re-
ducesIn as
thethe number
linear multi-step j, whichtemperature
of steps,method, approximatedeviation
the entropy interval
along 𝑠 , 𝑠 increases.
the polytropic path
reduces as the number of steps, j, which approximate the entropy interval [s1 , s2 ] increases.
( ) ( )
𝑇(𝑠) =j 𝑇 (𝑠) + 𝐸𝑟𝑇 (𝑠)
(35)
T (s) = ∑ j T i (s) + j ErTi (s)
(1) (1)
(35)
( )
i =1
where for the 𝑖th step, 𝑇 (𝑠) is defined in the sub-interval 𝑠 , 𝑠 .
(1)
where for the ith step, j T i (s) is defined in the sub-interval [si , si+1 ].
As shown in Figure 7, maximum temperature deviation along the polytropic path re-
ducesAsas
shown in Figure
the number 7, maximum
of steps temperature
increases. deviation
As expected, along the polytropic
linear approximation path
involves larger
reduces as the number of steps increases. As expected, linear approximation
deviation as compared with cubic polynomial approximation for the same number of steps involves
larger deviation
or segments. Asas compared
reviewed with cubicC,
in Appendix polynomial
a maximum approximation for the samesegments
of five cubic polynomial number has
of steps or segments. As reviewed in Appendix C, a maximum of five cubic polynomial
been documented to provide sufficient accuracy for a broad range of example cases by Evans
segments has been documented to provide sufficient accuracy for a broad range of example
[3].
cases by Evans [3].
The degraded energy in the polytropic compression process, 𝑇𝑑𝑠 , is approximated
by dividing the path into j-steps and using a linear approximant for each step.
( ) ( )
𝑇𝑑𝑠 = 𝑇 (𝑠) 𝑑𝑠 + 𝐸𝑟 𝑇𝐷𝑆 (36)
( )
where 𝐸𝑟 𝑇𝐷𝑆 shows the error in approximating the integral, 𝑇(𝑠) 𝑑𝑠, by using the
( )
linear approximant 𝑇 (𝑠) in the subinterval 𝑠𝑖 , 𝑠𝑖+1 . By summing up errors of all j-steps,
the total error for approximating the integral 𝑇𝑑𝑠, with linear approximants using j-steps
is developed:
R2
The degraded energy in the polytropic compression process, 1 Tds , is approximated
by dividing the path into j-steps and using a linear approximant for each step.
Z 2 j Z i +1
∑
(1) (1)
Tds = j T i (s) ds + j Eri TDS (36)
1 i =1 i
(1) R i +1
Int. J. Turbomach. Propuls. Power 2021,where j Eri TDS
6, x FOR PEER shows the
REVIEW error in approximating the integral, T (s) ds, by using
12 of 45
i
(1)
the linear approximant j T i (s) in the subinterval [si , si+1 ]. By summing up errors of all
R2
j-steps, the total error for approximating the integral 1 Tds, with linear approximants
using j-steps is developed:
( ) ( )
𝐸𝑟𝑇𝐷𝑆 = j 𝐸𝑟 𝑇𝐷𝑆 (37)
∑ j i
(1) (1)
j ErTDS = Er TDS (37)
i =1
Figure 7. Maximum
Maximum temperature
temperature deviation
deviationalong
alongthe
thepolytropic
polytropicpath
pathfor
forthe
theexample
examplecase
case14.
14.See
See
Table A1
Table A1 in
in Appendix
Appendix BB for
for gas
gas compositions
compositionsand
andcompressor
compressorinlet
inletand
anddischarge
dischargeconditions.
conditions.
By applying the
the relationship
relationship(26)
(26)for
foreach
eachstep,
step,and
andusing
usingthe
thecomposite
compositetrapezoidal rule,
trapezoidal
𝑇𝑑𝑠12, is
R
the integral,
rule, the integral, approximated
Tds, as follows:
is approximated as follows:
𝑗
( Ti(+ ))
j
!
Z 2 𝑇𝑖 + 𝑇𝑖+1
Ti+ ( )
= =∑
Tds𝑇𝑑𝑠 1
(s(𝑠 + (1) 𝐸𝑟
1 −−s𝑠i )𝑖 ) +
i +𝑖+1 𝑇𝐷𝑆
j ErTDS
(38)
1 2 2
i =1𝑖=1
( )
Using the
Using therelationships
relationships (2)(2)
andand (38),
(38), the the polytopic
polytopic
(1)
efficiency,
efficiency, 𝜂 , is developed as fol-
j η, is developed as follows:
lows:
j
si+𝑠1 − s−i 𝑠 Ti𝑇+ 𝑖 +Ti𝑇+𝑖+1

(1) ( )
j η 𝜂
= =
1 −1 ∑ hi+ℎ1 − h−i ℎ
−
22
1
(39)
(39)
i =1
( )
The deviation
deviation ofofthe
theapproximated
approximatedpolytropic efficiency, 𝜂j η,
polytropicefficiency,
(1) , from the actual polytropic
The from the actual poly-
( )
efficiency, 𝜂 using linear approximants, 𝑇 ((𝑠), 1) at each step is represented by the error
tropic efficiency,
( )
η p using linear approximants, j T i (s), at each step is represented by the
function, 𝐸𝑟𝜂(1, )in Equation (40).
error function, j Erη, in Equation (40).
( ) ( )
𝜂 = 𝜂 + 𝐸𝑟𝜂 (40)
(1) (1)
η p = j η + j Erη (40)
Because the linear method ignores the change in slope and concavity along the compres-
sion path, a large number of steps are needed to achieve the desired accuracy. Sandberg [23]
provided no criteria to determine the required number of linear steps to approximate 𝑇𝑑𝑠.
The uncertainty about the required number of steps in the linear multi-step method could
possibly be overcome by applying a large number of linear steps (50 to 100+ steps) to ensure
sufficient accuracy. However, pre-selecting 50 or 100+ linear polynomial steps has been
Because the linear method ignores the change in slope and concavity along the com-
pression path, a large number of steps are needed to achieve the desired accuracy. Sand-
berg [23] provided no criteria to determine the required number of linear steps to approxi-
R2
mate 1 Tds. The uncertainty about the required number of steps in the linear multi-step
method could possibly be overcome by applying a large number of linear steps (50 to 100+
steps) to ensure sufficient accuracy. However, pre-selecting 50 or 100+ linear polynomial
steps has been shown to be inefficient, (see Figure 11). Utilizing a large number of linear
steps significantly increases the calculation time due to the required nested iteration loops.
Figure A3b documents a comparison of elapsed time measurements for the supercritical
propane example case.
3.4. Polytropic Efficiency Calculation: Taher–Evans Cubic Polynomial Multi-Segment Method

The T-s polytropic path can be accurately approximated with a set of piecewise
thermodynamically coupled cubic polynomials. The T-s polytropic path is divided into “j”
segments over the entropy range, [s1 , s2 ], from compressor section inlet to discharge.
Each subinterval [si , si+1 ] of the actual path, T (s), is approximated by a cubic polynomial
(3)
approximant, j T i (s), which allows the temperature as well as T-s polytropic path slope,
(dT/ds), to be matched at endpoints of adjacent segments, thus increasing accuracy while
smoothing the overall path. The key factor in the accuracy of cubic polynomial multi-
segment method as compared with the linear multistep method is taking advantage of
applying the thermodynamic constraint, dTds , at each intermediate point (knot) to enable
smooth transition from one to the following segment.
The relationship (41) shows the actual temperature-entropy path, T (s), as approx-
(3)
imated with a set of “j” cubic polynomial approximants, j T i ( s ), and the temperature
(3)
deviation at each segment, j ErTi ( s ) :
j
∑
(3) (3)
T (s) = j T i ( s ) + j ErTi ( s ) (41)
i =1
(3)
where j T i (s) at each segment is determined by:
T (s) = Ai s3 + Bi s2 + Ci s + Di


 s i ≤ s ≤ s i +1
T ( s ) = T

i i



T (s 1 ) = Ti+1


 i+

 !
(3)
dT

Ti 1 + j η Xi (42)

 ds i = Ei = C pi (3)
1− j η

 !
(3)

1 + j η Xi + 1

Ti+1

dT
= Ei+1 = Cpi+1


 ds (3)
i +1 1− η

j
The coefficients ( Ai , Bi , Ci , Di ), of each cubic polynomial approximant are calculated

using the relationships (17) to (20), where subscripts 1 and 2 are replaced with endpoint
conditions of the ith segment.
R i +1 (3)
Similar to relationships (30) and (31), the degraded energy is calculated, i j T i ( s ) ds
(3)
and its related error, j Eri TDS, at each segment are defined. The total error for approximat-
R2
ing the integral 1 Tds, using cubic polynomial approximants is determined by summing
up errors of each segment.
j
∑
(3) (3)
j ErTDS = j Eri TDS (43)
i =1
R2
By applying the relationship (32) for each segment, the integral, 1 Tds, is approxi-
mated as follows:
j
Z 2
!
( Ti + Ti+1 ) ( Ei+1 − Ei )
Tds = ∑ 2 (3)
( s i +1 − s i ) − (si+1 − si ) + j ErTDS (44)
1 i =1
2 12
(3)
Using the relationships (2) and (44), the polytopic efficiency jη is developed as follows:
j
s i +1 − s i ( E − Ei )

Ti + Ti+1
= 1− ∑
(3)
jη − i +1 ( s i +1 − s i ) (45)
i =1
h i +1 − h i 2 12
(3)
The deviation of the approximated polytropic efficiency jη from the actual polytropic
(3)
efficiency, η p using cubic polynomial approximants j T i (s) at each step is represented by
(3)
the error function j Erη in Equation (46).
(3) (3)
ηp = jη + j Erη (46)
3.4.1. Polytropic Efficiency Calculation Procedure: Taher–Evans Cubic Polynomial

Multi-Segment Method
The solution requires an initial estimate of efficiency and then recursive solutions
for segment discharge temperature. Equation (45) is solved recursively by assuming an
(3)
initial value for the efficiency, j η to calculate endpoint slopes E1 and E2 , calculating
the resulting efficiency from (45), and minimizing the difference between assumed and
calculated efficiency. The step-by-step procedure is explained as follows:
Note: inlet and discharge pressures and temperatures discussed below are assumed
to be for total conditions at a given compressor section.
(1) Determine an initial estimate of polytropic efficiency by assuming ( E2 = E1 ) in the
relationship (45) for the cubic polynomial endpoint method.
Note: for the purpose of compressor performance testing, the initial value of the
efficiency is the expected efficiency provided by the compressor manufacturer.
(2) Select the number of segments (j).
There are different ways to select intermediate knots. In this paper, the overall pressure
ratio is divided across the number of segments (j) to determine equal pressure ratio
steps as shown in the next step. q
p2
(3) Use “j” number of equal pressure ratio segments. pr = j
p1 .
(4) Use Equation (29) to estimate the initial value of segment outlet temperature. (This
requires a solution to the cubic polynomial endpoint method).
Note: steps 1 through 4 are only required for the initialization of the solution algorithm.
(5) Calculate segment temperature rise, (∆Ti = Ti+1 − Ti )
(6) Calculate the segment polytropic efficiency using the relationship (45).
(7) Calculate efficiency deviation from the initially estimated polytropic efficiency from
step 1.
(8) Iterate segment outlet temperature to match assumed efficiency.
Note: different root finding algorithms may be used to optimize iterations for segment
outlet temperature.
Note: the acceptable efficiency convergence should be set less than 1E−6 .
(9) Move to next segment using outlet conditions of previous segment as inlet conditions.
(10) Repeat steps 5 through 9 for each segment.

(11) Compare final discharge temperature with the given value, Td
(12) If step 11 result is not within a small tolerance, iterate efficiency and determine new
values for coefficients A, B, C and D.
(13) Repeat steps 5 through 11 until agreement is reached within the tolerance.
Note: the acceptable convergence from the given discharge temperature should be set
less than 1E−8 .
In the Taher–Evans Cubic Polynomial method, intermediate points are determined by
setting one thermodynamic variable (such as pressure or entropy) at each end of a segment
and calculating other thermodynamic variables to force the segment endpoint to lie on a
constant efficiency path via nested recursive algorithms (within a small error tolerance).
Using cubic polynomial approximants, the polytropic path slope can be accurately calcu-
lated since the first derivative of each segmental polynomial can be determined at each
intermediate point (see Ei and Ei+1 in Equation (42)).
The difference between final segment calculated discharge temperature and the known
discharge temperature, Td , can be minimized by iterating efficiency. These nested itera-
tions of efficiency (inner loop) and the final segment discharge temperatures (outer loop)
converge quickly.
When a final efficiency is determined, summations of enthalpy changes and Tds losses
can be performed for the entire set of segments. The cubic polynomial multi-segment
method provides intermediate points along the constant efficiency temperature-entropy
path that have been forced to coincide with the actual path. Additionally, the path slope at
each intermediate point is equal at the end of the ith segment and the beginning of the i+1
segment. See Table 1.
3.4.2. Distinction between Cubic Interpolation and Approximation

As the final point in this section, it should be noted that the approximation of functions
is different from interpolation. In the case of interpolation, at certain intermediate points,
the value of an unknown function is given (such as measured values at intermediate points),
but the function itself is not known. Interpolating cubic splines are used to interpolate
between known intermediate points, which is irrelevant in the problem of approximating
an unknown function without any given intermediate point. Accordingly, the descriptor
“spline” is not applied in this paper as it usually is thought of as being used to approximate
a curve passing through previously known knots.
Table 1. Cubic polynomial 5-segment method as applied for Case 10 using AspenTech Hysys V.10 [27] and PR-LK equation of state. See Table A1 in the Appendix B for gas compositions
and compressor inlet and discharge conditions for Case 10.
Rs
T p s h i+1 (3) T (s)ds
Segment Number 5 Ai 5 Bi 5 Ci 5 Di 5 Ei 5 Ei+1 5 ηp,i 5 wp,i si 5 i
i [K] [bar] [kj/kg-K] [kj/kg] - - - - 2
[kg-K /kj] 2
[kg-K /kj] - [kj/kg] [kj/kg]
Inlet 366.483 3.390 20.684 −2276.711 −184,369 1,886,535 −6,432,843 7,310,148 1499 1650 0.81147 13.28549 3.08654
1 Outlet 379.567 3.398 26.316 −2260.339
Inlet 379.567 3.398 26.316 −2260.339 151,457 −1,534,519 5,183,886 −5,838,652 1650 1831 0.81147 13.19405 3.06530
2 Outlet 393.341 3.406 33.480 −2244.080
Inlet 393.341 3.406 33.480 −2244.080 86,118 −868,526 2,921,119 −3,276,023 1831 2020 0.81147 13.03013 3.02721
3 Outlet 407.879 3.414 42.596 −2228.022
Inlet 407.879 3.414 42.596 −2228.022 −1,407,357 14,438,426 −49,373,634 56,277,303 2020 2158 0.81147 12.82697 2.98001
4 Outlet 423.124 3.421 54.193 −2212.215
Inlet 423.124 3.421 54.193 −2212.215 −4,606,378 47,325,393 −162,069,532 185,004,389 2158 2187 0.81147 12.67478 2.94466
5 Outlet 438.706 3.428 68.948 −2196.596
Summary of all 5 Segments 0.81147 65.01143 15.10372
is different from interpolation. In the case of interpolation, at certain intermediate poin
value of an unknown function is given (such as measured values at intermediate point
the function itself is not known. Interpolating cubic splines are used to interpolate be
known intermediate points, which is irrelevant in the problem of approximating an unk
function without any given intermediate point. Accordingly, the descriptor “spline”
applied in this paper as it usually is thought of as being used to approximate a curve p
through previously known knots.
4. Polytropic Path on Temperature-Entropy Diagram
4. Polytropic Path on Temperature-Entropy Diagram
For an uncooled compressor section, entropy increases along the polytropic path.
The rate of changeFor an uncooled with
of temperature compressor
entropysection,
at each entropy increases
point along along thepath
the polytropic polytropic
is pat
rate of change of temperature with entropy
determined with the relationship (A10) from Appendix A. at each point along the polytropic path is
mined with the relationship (A10) from Appendix A.
T 1 + ηp X

dT
= E =𝑑𝑇 𝑇 1+𝜂 𝑋 (A10)
ds η p c p 1=−𝐸η=p (A
𝑑𝑠 𝑐 1−𝜂
The change ofThe temperature
change ofwith entropy from
temperature with the compressor
entropy from the section inlet to section
compressor dischargeinlet to disc
varies as temperature, specific heat, c , and compressibility function, X, change
varies as temperature, specific heat, 𝑐 , and compressibility function, 𝑋,
p at change
each at each
point along the compression
along path. Atpath.
the compression conditions near to the
At conditions critical
near to thepoint, thepoint,
critical rate ofthe
change
rate of change
of c p and X is usually large. This has significant impact on the change of temperature with
and 𝑋 is usually large. This has significant impact on the change of temperature wi
entropy as the relationship (A10) suggests. The supercritical propane example case has
tropy as the relationship (A10) suggests. The supercritical propane example case ha
inlet conditions extremely close to the critical point and the thermodynamic properties
conditions extremely close to the critical point and the thermodynamic properties are v
are varying rapidly as compression begins. Figures 8 and 9 illustrate the variation of
rapidly as compression begins. Figures 8 and 9 illustrate the variation of specific he
specific heat, c p , and compressibility function, X, respectively, in the vicinity of the critical
and compressibility function, 𝑋, respectively, in the vicinity of the critical point. Each
point. Each isobar has a distinct peak. The inlet conditions lie in an area on the flank of a
has a distinct peak. The inlet conditions lie in an area on the flank of a peak while disc
peak while discharge conditions (not shown) are far removed from the critical area peaks.
conditions (not shown) are far removed from the critical area peaks. Equations (15) an
Equations (15) and (16) would show a large difference in slope of the T-s path between
would show a large difference in slope of the T-s path between inlet, 𝐸 , and dischar
inlet, E1 , and discharge, E2 , as documented in Table A2. These large differences are one
as documented in Table A2. These large differences are one indicator of the difficult
indicator of the difficulty of a compression case.
compression case.
Figure 8. Change of specific heat, c p , along the polytropic path for supercritical propane example.
Moreover, the change of specific heat, c p , with temperature at different inlet pressures is shown.
Int. J. Turbomach. Propuls. Power 2021, 6, x FOR PEER REVIEW
Figure 8. Change of specific heat, 𝑐 , along the polytropic path for supercritical propane example. Moreover, the change o
specific heat, 𝑐 , with temperature at different inlet pressures is shown.
Figure 9. Change of compressibility function, 𝑋, along the polytropic path for supercritical propane example. Furthermor
Figure 9. Change of compressibility function, X, along the polytropic path for supercritical propane
the change of compressibility function, 𝑋, with temperature at different inlet pressures is shown.
example. Furthermore, the change of compressibility function, X, with temperature at different inlet
pressures is shown.
While the change of temperature with entropy along the compression path of a
cooled
While the change compressor
of temperaturesectionwith
is always
entropy positive
along the( compression
> 0), the curvature
path ofofanthe polytropi
un-
on thesection
T-s diagram maypositive dT
change(sign>(from
cooled compressor is always ds 0), theconcave
curvatureupward
of theto downward
polytropic or vice ver
path
on the T-s diagram
changemay change
in sign sign (from
happens whenconcave upward
equals to downward
to zero at any point or vice
alongversa). A
the compression
d2 T
change in sign happens when
which indicates a change
ds2
equals to zero at any point along the compression path
of curvature (as seen in the case of HP Ethylene in Appendix
which indicates a Thechange
cubicof polynomial
curvature (as seen in method
endpoint the casecanof HP Ethylene
be used in Appendix
to approximate theB).
concavity
The cubic polynomial
actual endpoint
T-s polytropic method
path. No othercanendpoint
be used to approximate
methods the concavity
can provide of
such a mathematic
the actual T-s thermodynamic
polytropic path. No otherabout
insight endpoint methods can
the behavior provide
of the actualsuch
T-s apolytropic
mathematicalpath from k
and thermodynamic insight conditions
thermodynamic about the behavior
at the inletofand
thedischarge
actual T-sofpolytropic path section.
the compressor from
known thermodynamic conditions
The concavity at the
of the inlet andapproximant
polynomial discharge of canthe compressor
be examinedsection.
using the second
The concavity of of
ative test thethepolynomial
path equation approximant can be examined
with the following criteria: using the second
derivative test• of the path
If (−𝐵 ⁄3𝐴equation
) 𝑖𝑠 out with
of thetherange
following
of 𝑠 criteria:
, 𝑠 and 3𝐴𝑠 + 𝐵 > 0 the upward concav
• ( )
If (− B/3A) isthe
out of the
cubic polynomial 1 , s2 ] and 3As1𝑇+
range of [sapproximant, B>
(𝑠), is 0unchanged
the upward concavity
along of See Fig
the path.
(3)
curve I. approximant, 1 T (s), is unchanged along the path. See Figure 10
the cubic polynomial
curve I. • If (−𝐵⁄3𝐴) is out of the range of 𝑠 , 𝑠 and 3𝐴𝑠 + 𝐵 < 0 the downward con
• If (− B/3A) isofout
theofcubic
the range
polytropic s2 ] and 3As1 (+)𝑇B(𝑠),
of [s1 ,approximant, < 0isthe downward
unchanged concavity
along the path.of See Fig
( 3 )
curve II.approximant, T (s), is unchanged along the path. See Figure 10
the cubic polytropic 1 (
curve II. • If 𝑠 < (−𝐵⁄3𝐴) < 𝑠 the concavity of the cubic polytropic approximant,
( ) (3)
• changes
If s1 < (− B/3A ) < at
s2 the
theinflection
concavitypoint of the cubic (−𝐵⁄3𝐴). See
𝑠 = polytropic Figure 10 curve
approximant, III.
1 T ( s ),
(3)
Downward
changes at the concavity
inflection point j s inisf an
= indication
(− B/3A). of increased
See Figure 10difficulty for the compression c
curve III.
shown and confirmed by example cases in Appendix B. See Appendix
Downward concavity is an indication of increased difficulty for the compression B for several com
case
sion cases evaluated.
as shown and confirmed by example cases in Appendix B. See Appendix B for several
compression cases evaluated.
The actual location of the inflection point on the T-s polytropic path can be more
precisely approximated by applying the cubic polynomial multi-segment method.
The rate of change of dTds with entropy can be considered as a screening criterion to
determine the difficulty of a compression case.
The curve shape of the polytropic path as predicted by the cubic polynomial endpoint
method provides an easy to apply criteria to determine how many multi-segments to
employ to assure an acceptable relative deviation. Based on 115 cases covering 22 fluids
(7 pure fluids and 15 mixtures reviewed in Appendix C) studied by Evans [3], the number
Int. J. Turbomach. Propuls. Power 2021, 6, xof cubic
FOR PEERpolynomial
REVIEW segments required to achieve ≤0.001% relative deviation for the 19 of 45
polytropic efficiency is three segments for Category I and five segments for Categories II
and III.
Figure
Figure 10. Three
10. Three different
different categories
categories areare identifiedbased
identified based upon
upon polytropic
compressionpath shapes
path on on
shapes the the
T-s diagram. The The
T-s diagram.
number
number of cubic
of cubic polynomial
polynomial segments
segments requiredtotoachieve
required ≤0.001% relative
achieve ≤0.001% relative deviation
deviationfor
forthe
thepolytropic
polytropicefficiency is shown.
efficiency is shown.
5. The actualChange
Enthalpy location of the
with inflection
Entropy alongpoint on the T-sPath
the Polytropic polytropic path can be more pre-
cisely approximated
The relationshipby (2)
applying the cubic polynomial
can be rearranged to multi-segment method.
The rate of change of with entropy can be considered as a screening criterion to de-
dh T
termine the difficulty of a compression case. = (47)
ds 1 − ηp
The curve shape of the polytropic path as predicted by the cubic polynomial endpoint
methodAs provides an easyoftothe
the efficiency apply criteriaprocess
polytropic to determine howthe
is constant, manyslope multi-segments
of the tangent to employ
the
to assure an acceptable relative deviation. Based dh
path increases as the temperature rises along theon 115 cases covering
compression path ( ds 22
> fluids (7 pure fluids
0). Rearranging
andthe
15relationship
mixtures reviewed
(47) and in Appendix
taking C) studied
the derivative withbyrespect
Evansto[3],
ds the number of cubic polyno-
gives:
mial segments required to achieve"≤0.001% relative deviation for the polytropic efficiency is
dT dh d2 h
#
three segments for Category I and T segments
d five − T
for ds
Categories
2 II and III.
dh
= ds ds 2 =0 (48)
ds dh
ds
ds
5. Enthalpy Change with Entropy along the Polytropic Path
TheByrelationship
replacing (2)
T can be
with dhrearranged
, the secondto
derivative of the function h(s) along the poly-
1− η p ds
tropic compression path appears: 𝑑ℎ 𝑇
= (47)
d2 h𝑑𝑠 11− 𝜂dT
= (49)
ds2 1 − η p ds
As the efficiency of the polytropic process is constant, the slope of the tangent to the
path increases as the temperature
compressorrises along
dT the compression d2 hpath ( > 0). Rearranging the
For an uncooled section, ds > 0. Therefore, ds2 > 0. As shown in Figure 1,
relationship (47) compression
the polytropic and taking the
pathderivative with
defined with respect
function to is𝑑𝑠
h(s) gives:
always concave upward along
the compression path since the second derivative is positive in an uncooled compressor
𝑑𝑇 𝑑ℎ 𝑑 ℎ
− the
𝑇 thermodynamic value of d 2h at any
2
𝑑 𝑇 A provides
section. Equation (A11) in Appendix 𝑑𝑠 𝑑𝑠 𝑑𝑠 ds
= =0 (48)
𝑑𝑠 𝑑ℎ path. 𝑑ℎ
point along the polytropic compression
𝑑𝑠 𝑑𝑠
By replacing with , the second derivative of the function h(s) along the poly-
tropic compression path appears:
Int. J. Turbomach. Propuls. Power 2021, 6, x FOR PEER REVIEW 2

the compression path since the second derivative is positive in an uncooled compresso
tion. Equation (A11) in Appendix A provides the thermodynamic value of at any p
along the polytropic compression path.
The change ofThe
enthalpy
changewith entropywith
of enthalpy at any pointatalong
entropy the polytropic
any point path can be
along the polytropic path can b
approximated proximated
by substituting the relationship
by substituting of temperature
the relationship from (12) from
of temperature to (47):
(12) to (47):
dh As3 + Bs𝑑ℎ 𝐴𝑠 + 𝐵𝑠 + 𝐶𝑠 + 𝐷
2 + Cs
= +D (5
= 𝑑𝑠 1−𝜂 (50)
ds 1 − ηp
By integratingBy integrating
(50), (50), the of
the relationship relationship of enthalpy
enthalpy and entropyand
canentropy can be approximated
be approximated at a
point along the polytropic
any point along the polytropic path: path:
s4 − s41 + B3 𝐴s3(𝑠−4 −
A
C𝐵 2 𝐶
s31𝑠41+
) +2 s(𝑠3−−s𝑠2131 ) +

+ D(𝑠s2−
− s𝑠121) + 𝐷(𝑠 − 𝑠1 )
h ( s ) = h1 + 4
ℎ(𝑠) = ℎ1 + 4 3 2 (51) (5
1 − ηp 1−𝜂
where, s represents any entropy value between compressor section inlet, s1 and section
where, 𝑠 represents any entropy value between compressor section inlet, 𝑠 and section
discharge, s2 .
charge, 𝑠 .
6. Comparison of Cubic and Linear Polynomial Results
6. Comparison of Cubic and Linear Polynomial Results
The Taher–Evans Cubic Polynomial 10 segment method is used as a comparison
The Taher–Evans Cubic Polynomial 10 segment method is used as a compa
standard to determine relative deviations in polytropic efficiency for all other methods. See
standard to determine relative deviations in polytropic efficiency for all other methods
Equation (A12) for the relative deviation definition. The number of recommended cubic
Equation (A12) for the relative deviation definition. The number of recommended cubic
segments basedments
uponbased
T-s path
uponcurve shapes
T-s path listed
curve in Figure
shapes listed10inwas developed
Figure from results
10 was developed from resul
reported by Evans [3]. Figure 11 shows that determining a similar set of recommendations
ported by Evans [3]. Figure 11 shows that determining a similar set of recommendation
for the numbertheof number
linear polynomial steps required
of linear polynomial stepswas not possible.
required was not The wide The
possible. scatter
wideof scatter o
the linear methods’ results that would yield similar accuracy to the cubic methods
linear methods’ results that would yield similar accuracy to the cubic methods was was
very random. random.
FigureFigure 11. Cubic

11. Cubic and linear
and linear polynomial
polynomial solution
solution results
results comparison
comparison forfor
115115 example
example cases.
cases.
Appendices B and C discuss example cases and compare various polytropic calculation
methods. The 19 example cases in Appendix B were selected for illustration of a range
of relatively easy and difficult polytropic efficiency calculation applications. Several of
these are used to show various characteristics of the calculation results. Appendix C is a
review of 115 example cases studied by Evans [3], and includes the 19 example cases from
Appendix B. This larger population further confirms and documents the trends discovered
when applying the multi-segment cubic and multi-step linear polynomial methods.
7. Conclusions
• The Taher–Evans Cubic Polynomial method (TE-CP), defined, described and tested
herein, illustrates that a highly accurate calculation method for real gas centrifugal
compressor polytropic performance efficiency has been developed and implemented
that employs a temperature—entropy cubic polynomial path function.
# Both endpoint and sequential segment versions are described and tested.
# Previously published polytropic efficiency calculation methods including first
degree linear polynomial methods are shown to be less accurate and/or slower
to achieve a solution than Taher–Evans Cubic Polynomial methods.
# The TE-CP methods’ superior results are due to the cubic polynomial path
having additional thermodynamic constraints applied to more accurately ap-
proximate the actual path slope at any point along the path. Determination
of path slope at the compression endpoints is independent of performance
calculation method.
# T-s polytropic path curve shapes can be determined from derivatives of the
cubic polynomial path equation and are indicative of calculation relative diffi-
culty. The curve shape is used as a criterion to select the required number of
cubic segments. No other polytropic method can provide such a mathematical
and thermodynamic insight into the fluid compression process.
# The number of T-s path cubic polynomial segments required to achieve an ac-
ceptable relative deviation ≤0.001% for polytropic efficiency is three segments
for category I and five segments for categories II and III curve shapes.
# As shown in Figure A10, the cubic polynomial methods provide low uncer-
tainty with only a few cubic segments. Uncertainty is based upon polytropic
efficiency relative deviation, and its magnitude is related to the size of the
statistical population.
# Cubic polynomial methods provide continuous equations to plot the polytropic
path on T-s and h-s diagrams. This is a very unique feature of the cubic
polynomial method that enables visualizing the polytropic path.
# The cubic path coefficients (A, B, C, and D) in Equation (12) provide mean-
ingful insights about the behavior of the polytropic path on the temperature-
entropy diagram.
# The Taher–Evans Cubic Polynomial methods are highly suitable for application
to compressor performance testing according to ASME PTC 10 [4].
• Cubic and linear polynomial calculation methods have been extensively documented
and compared.
# Polytropic efficiency calculations have been reviewed for 17 methods across
115 example cases yielding a total of 1955 independent calculations that validate
superior results are achieved when using both endpoint and sequential segment
Taher–Evans Cubic Polynomial (TE-CP) polytropic path methods. A total of
22 fluids were employed including seven pure fluids and 15 fluid mixtures to
cover wide ranges of mole weight and critical pressures and temperatures.
# Cubic polynomial endpoint path methods were demonstrated to achieve better

accuracy than any other endpoint polytropic efficiency calculation method.
# Cubic polynomial sequential segment path methods were demonstrated to
achieve high accuracy acceptable polytropic efficiency results significantly
faster with fewer segments required than linear polynomial path methods
based upon 115 compression example cases.
# For all 115 example cases, T-s polytropic paths agreed between 10 cubic seg-
ments and 100 linear steps which confirmed the validity of using the result of
10-segment cubic polynomial approximants as the standard for comparison for
polytropic efficiency relative deviation and uncertainty.
# Cubic polynomial path methods show that the required number of segments
to achieve an acceptable result remains within a maximum of five segments
or less for all the cases studied. For the linear polynomial path method, the
required number of steps did not follow a regular pattern as shown in Figure
11. No rule or algorithm can be determined for the required number of linear
polynomial steps to achieve an acceptable result.
# While both cubic and linear polynomial methods can reach acceptable results,
the cubic methods require fewer repetitive iterations and are completed much
faster. Sequential segment cubic polynomial method calculations are typically
at least an order of magnitude faster than sequential step linear polynomial
methods. This is a clear advantage for cubic methods when processing com-
pressor performance test data in real time.
# The cubic polynomial methods provide an order of magnitude lower uncer-
tainty with significantly fewer cubic segments than linear polynomial methods.
(See Figure A10).
# The cubic polynomial methods’ uncertainty rapidly decreases as the number
of segments increases while the linear polynomial methods illustrate a much
more gradual decrease in uncertainty as they approach an asymptote. The
linear polynomial methods uncertainty asymptote appears to be more than an
order of magnitude greater than the cubic methods.
• Identical calculations for selected example cases using several different equations of
state to obtain thermodynamic state point data for application of cubic and linear
polynomial methods revealed significant differences in results for efficiency can occur
based upon the chosen EOS. See Appendix D for selective examples.
Recommendations
• When rapid, highly accurate, low uncertainty, centrifugal compressor polytropic
efficiency calculations are required, the Taher–Evans Cubic Polynomial methods (TE-
CP) should be used.
• Due to significant differences found between results when employing various equa-
tions of state, it is recommended that a project use the same equation of state through-
out the project timeline to avoid potential discrepancies.
Author Contributions: Conceptualization and methodology, M.T.; software, validation, formal

analysis, writing—original draft preparation, writing—review and editing, M.T. and F.E.; All authors
have read and agreed to the published version of the manuscript.
Funding: This research received no external funding.
Institutional Review Board Statement: Not applicable.
Informed Consent Statement: Not applicable.
Data Availability Statement: Not applicable.
Conflicts of Interest: The authors declare no conflict of interest.
Nomenclature
A Polynomial approximant coefficient

B Polynomial approximant coefficient
b Linear approximant coefficient

1 ∂v
β Isobaric expansivity, = υ ∂T p
BWRS Benedict–Webb–Rubin EOS as modified by K. E. Starling
C Polynomial approximant coefficient
CH4 Methane
CO2 Carbon Dioxide
C3 H8 Propane
CP Cubic Polynomial
Cp heat capacity at constant pressure
D Polynomial approximant coefficient
dT
E Slope of the polytropic path on the T-s diagram, ds ηp
EOS Equation of State
Erη Deviation of polytropic efficiency
ErT Deviation of approximated temperature
R
ErTDS Deviation of approximated Tds
EOS Equation of State
h Specific enthalpy
HP High pressure
OEM Original Equipment Manufacturer
PR Peng-Robinson EOS
Peng-Robinson EOS and the Lee-Kesler EOS for the calculation
PR-LKP
of enthalpies
LP Linear Polynomial, or Low Pressure
LKP Lee-Kesler-Plöcker EOS
m Linear approximant coefficient
MP Medium Pressure
n Schultz [6] polytropic volume exponent as applied in pvn = c
T Temperature
TE-CP Taher–Evans Cubic Polynomial
REFPROP Reference properties EOS software
s Specific entropy
SC Supercritical
v Specific volume
X Compressibility function, = βT − 1
XY Schultz method using compressibility functions X and Y
ηp Polytropic efficiency
Subscripts and Superscripts
(q)
j T i (s)
i indicates the ith segment
j total number of segments
q degree of the polynomial approximant
(3)
Example: 10 T 6 (s) indicates the cubic polynomial approximant, which is
used to approximate the segment number 6 of a polytropic path that is
divided into 10 segments.
This symbology provides a clear way to compare different methods
discussed in this paper.
inf inflection
p Polytropic, or pressure
r ratio
d compressor section discharge (measured)
1 compressor section inlet
2 compressor section discharge
Appendix A. Temperature Change with Entropy along the Polytropic Path

2
Derivation of the following two important relationships, dT d h
ds and ds2 along the poly-
tropic compression path is explained in this appendix.
Choose the entropy to be a function of temperature and pressure. Writing the differen-
tial of s = s( T, P) gives:
∂s ∂s
ds = dT + dP (A1)
∂T P ∂P T
The first partial derivative in (A1) can be written as follows:

∂h
cp

∂s ∂T P
= = (A2)
∂T P ∂h T
∂s P

∂h
where = c P and ∂h
∂T ∂s P = T from the first partial derivate of the Gibbs relationship
P
dh = Tds + vdP.
The second partial derivative in (A1) is one of the Maxwell relations:

∂s ∂v
=− (A3)
∂P T ∂T P
Substituting the first and second partial derivatives in (A1) from (A2) and (A3) yields:
cp

∂v
ds = dT − dP (A4)
T ∂T P

Using the definition of isobaric expansivity β = υ1 ∂T
∂v
for the second partial deriva-
p
tive gives:
cp
dT − βvdP
ds = (A5)
T
Employing the Gibbs relationship dh = Tds + vdP for a polytropic process and
substituting dh with vdP
η p gives:
Tds
vdP = 1 (A6)
ηp − 1
Substituting vdP from (A5) to (A4) yields:
cp βTds
ds = dT − 1 (A7)
T ηp − 1
This can be further simplified to

 
βT c
1 + ds = p dT (A8)
1 T
ηp − 1
From (A8) the slope of the polytropic path on a T-s diagram is determined:
1 + η p ( βT − 1)

dT T
= (A9)
ds cp 1 − ηp
Substituting ( βT − 1) with the Schultz compressibility function X gives:
T 1 + ηp X

dT
= (A10)
ds cp 1 − ηp
dT
By replacing ds from (A10) to (47) gives:

d2 h T 1 + ηp X
= (A11)
ds2 cp 1 − ηp 2

The relationships (A10) and (A11) were developed and presented for the first time by
Taher [1].
The relationship (A11) shows that the polytropic path on the enthalpy-entropy dia-
gram is concave upward since the second derivative is positive in an uncooled compres-
sor section.
Appendix B. Comparative Numeric Examples of Application of the Taher–Evans

Cubic Polynomial Methods
The Taher–Evans Cubic Polynomial method for calculation of centrifugal compressor
polytropic performance has been applied to 19 compression example cases which are com-
prised of five pure fluids and four hydrocarbon fluid mixtures. Eight cases are considered
low pressure while eleven cases are considered medium or high pressure. Input data for
ten cases came from existing literature references while nine cases are new. All cases were
chosen based upon the desire to explore a wide range of process conditions for fluids
common to the hydrocarbon industry.
The fluids and process conditions selected represent a combination of relatively easy
and difficult calculations for centrifugal compressor polytropic performance. Table A1
lists input data for the 19 cases. The “easy cases” are found by the results listed below to
be low pressure versions for CO2 , C3 H8 and CH4 /C3 H8 /CO2 . The most difficult case is
supercritical propane due to inlet conditions.
Using results of the above referenced analyses, this appendix will review polytropic
efficiency calculation method comparisons for following three categories:
• Endpoint methods
# Schultz (two variations) [6];
# Mallen Saville [9];
# Cubic Polynomial;
# Linear polynomial.
• Three point methods
# Cubic Polynomial two segment;
# Huntington’s three point [13].
• Multi-step and multi-segment methods
# Cubic Polynomial—sequential segments (2 through 10 segments);
# Linear Polynomial—sequential steps (10, 20, 50, 90 and 100 steps).
The Taher–Evans Cubic Polynomial 10 segment method is used as a comparison
standard to determine relative deviations in polytropic efficiency for all other methods.
The validity of using this standard is confirmed below. Equation (A12) documents the
relative deviation calculation.

ηTE-CP 10 Segment − η p
Relative Deviation, % = 100 × (A12)
ηTE-CP 10 Segment
Industry standardized compressor data sheets document operating conditions and

performance characteristics with polytropic efficiency shown as a percentage, typically
with two decimal places of accuracy, such as 68.25%. One criterion applied in the example
cases below is that the deviation should be known to be less than 0.001% of an acceptable
standard. Since all polytropic efficiency calculation methods are approximations, a standard
determines the magnitude of the difference between two calculations.
As mentioned in the paper (see also Appendix D), polytropic performance results
can be impacted by the choice of an equation of state (EOS). The REFPROP (version 10)
software package [28] produced by the National Institute of Science and Technology (NIST)
was employed for calculations in this appendix. REFPROP provides reference quality,
multi parameter, real gas equations of state for thermodynamic and transport properties of
individual fluids. REFPROP applies optimized binary interaction parameters via mixing
rules to determine mixture properties.
When using the REFPROP software package, all example case inlet and discharge state
point phases are single phase and either superheated or supercritical. For some example
cases, the close proximity of inlet conditions to a fluid’s critical point yielded two phase
or liquid phases for some other equations of state used in Appendix D. The polytropic
performance methods discussed in this paper are neither designed nor recommended for
two phase flows. The impact of equation of state choices will be reviewed in Appendix D.
Table A1. Example case input data.

Case Number 1 2 3 4 5 6 7 8 9 10 11
HP Ethylene
LP Ethylene
HP Propane
Schultz R12
SC Propane
LP Propane
PTC10 CO2
SC Ethane
MP CO2
HP CO2
LP CO2
Case Name
Ethane 1.0000
Composition, mol %
Propane 1.0000 1.0000 1.0000
Carbon
1.0000 1.0000 1.0000 1.0000
Dioxide
Ethylene 1.0000 1.0000
R12 1.0000
p1 , psia 10 360 362.5 750 300.01 400 500 1100.1 17.99 300 650
p1 , bara 0.689 24.821 24.993 51.710 20.685 27.579 34.474 75.849 1.240 20.684 44.816
Compression Conditions
T1 , ◦ F −10 50 98.3 110 100 100 75 98.3 −31.36 200 210
T1 , ◦ C −23.33 10 36.833 43.333 37.778 37.778 23.88 36.833 −35.2 93.333 98.889
p2 , psia 130 1000 7250 3500 487.76 1200 3000 6000.3 60 1000 3500
p2 , bara 8.963 68.948 499.870 241.316 33.630 82.737 206.843 413.706 4.137 68.948 241.316
T2 , ◦ F 210 195 566.3 285 201.59 325 400 368.3 60 330 300
T2 , ◦ C 98.889 90.556 296.833 140.556 94.217 162.778 204.444 186.833 15.556 165.556 148.889
Case Number 12 13 14 15 16 17 18 19
MP CH4 /C3 H8 /CO2
HP CH4 /C3 H8 /CO2

LP CH4 /C3 H8 /CO2
HP Heavy NG
PTC10 HPNG
LP Heavy NG
HP CH4 /CO2
LP CH4 /CO2
Case Name
Methane 0.30294 0.30294 0.30294 0.256274 0.256274 0.688671 0.688671 0.8600
Ethane 0.03748 0.03748 0.03748 0.022871 0.022871 0.119957 0.119957 0.1125
Propane 0.43533 0.43533 0.43533 0.005691 0.005691 0.102964 0.102964 0.0075

Composition, mol %
isoButane 0.00222 0.00222 0.00222 0.000001 0.000001 0.031432 0.031432
Butane 0.00218 0.00218 0.00218 0.000001 0.000001 0.027541 0.027541
isoPentane 0.000001 0.000001 0.008927 0.008927
Pentane 0.000001 0.000001 0.005563 0.005563
Hexane 0.000175 0.000175
Nitrogen 0.00399 0.00399 0.00399 0.002629 0.002629 0.002477 0.002477 0.0040
Carbon
0.21586 0.21586 0.21586 0.712356 0.712356 0.012468 0.012468 0.0160
Dioxide
p1 , psia 650 2071 2071 500 1649.7 500 1750 2520.6

Compression Conditions
p1 , bara 44.816 142.790 142.790 34.474 113.743 34.474 120.658 173.789
T1 , ◦ F 115 160 160 50 101.2 100 50 100
T1 , ◦ C 46.111 71.111 71.111 10 38.444 37.778 10 37.778
p2 , psia 2200 10,025.70 11,010.80 1700 8189.70 1800 4375 6500
p2 , bara 151.685 691.248 759.168 117.211 564.660 124.106 301.646 448.159
T2 , ◦ F 270 291.9 301.5 250 316.2 295 125 280

Appendix B.1. Endpoint Methods

Schultz [6] described several alternative versions of his methodology in his original
ASME paper all of which applied a path equation of pvn = c. The variations involved
calculating his “polytropic head factor”, f , in various ways as well as multiple methods of
averaging his “polytropic volume exponent”, n. For the comparisons below, the following
two variations of Schultz’s method are included for examination since they were included
in the 1997 version of ASME PTC10 [4].
• Schulzn—with n based upon the endpoints of compression
ln( p2 /p1 )
# n= ln(v2 /v1 )
• Schultz XY average n—with n averaged based upon individual calculations at inlet
and discharge
n1 + n2
# n= 2
Mallen and Saville [9] employed a path equation of Tds dT = c. Path equations are
described above in this paper for the liner and cubic polynomial methods.
Figure A1a–e illustrate polytropic efficiency relative deviation results for endpoint
methods for the 19 example cases. The Taher–Evans endpoint method provides superior
results compared to the other endpoint methods for all example cases. Example cases
PTC10 CO2 , LP CO2 , LP propane and LP CH4 /C3 H8 /CO2 are considered easy cases
since the Taher–Evans Cubic Polynomial endpoint method achieves deviations from the
10-segment method of less than 0.001%. Thus, using multi-segment analyses would not
be required for these cases. The two Schultz based methods show significant errors for
some example cases and do not exhibit a consistent pattern. Figure A1f directly compares
two related cases in that the PTC10 CO2 case represents a Type 2 test designed for the
HPNG example case. The PTC10 CO2 case is a low-pressure case and deviations are seen
to be smaller than for the HPNG case as expected. However, the Taher–Evans endpoint
method provides better accuracy for the PTC 10 HPNG case. The HPNG case requires
more segments to achieve less than a 0.001% deviation level.
Notable comparisons for endpoint methods:
• Figure A1a
# HP ethylene; Mallen and Saville [9] highlighted their improvements over
Schultz’s method for high pressure, but still have an approximate 1% deviation;
# LP and HP ethylene; typical and expected patterns for cubic and lin-
ear polynomials.
• Figure A1b,c
# CO2 and propane; as pressure increases accuracy decreases is a typical pattern
• Figure A1d
# CH4 /C3 H8 /CO2 ; medium and high pressure cases are very similar since
discharge conditions are close
• Figure A1f
# PTC10; the HPNG case will require more cubic segments to achieve sufficient
accuracy than the Type 2 test design using CO2
Figure A1. (a–f) Summary comparison of relative deviation results for endpoint methods.
Figure A1. (a–f) Summary comparison of relative deviation results for endpoint methods.
Table A2 lists the T-s path inlet and discharge slope and curve category for the 19 exam-
TableFive
ple cases. A2 cases
lists show
the T-s path
that inlet and
an inflection discharge
point slope
in curvature andalong
exists curvethecategory
path. For for the 19
these
example
five cases, endpoint analyses do not provide sufficient polytropic efficiency accuracy and path.
cases. Five cases show that an inflection point in curvature exists along the
For these five analyses
multi-segment cases, endpoint analysesChanges
are recommended. do not inprovide sufficient
path slope polytropic
and curvature can be efficiency
used
accuracy
as criteriaand multi-segment
for determining analyses
how many cubicare recommended.
segments are needed toChanges in path slope
provide sufficient accu- and
curvature can be 4,
racy (See Section used as criteria
Figure for determining
10 and Appendix C). how many cubic segments are needed to
provide sufficient accuracy (See Section 4, Figure 10 and Appendix C).
Table A2. Taher–Evans Cubic Polynomial endpoint method T-s path slope and curvature.
Case Number Case 1 Case 2 Case 3 Case 4 Case 5

T-s Slope at Inlet, E1 , lbm.R2 /BTU 14,168 10,442 9383 14,537 7026
T-s Slope at Discharge, E2 , lbm.R2 /BTU 19,822 12,253 8302 7946 7881
T-s Slope Change, % 39.90 17.34 −11.52 −45.34 12.17
T-s Curve Shape Category I I III II I
Entropy at Inflection Point, BTU/lbm.R - - 0.65070 - -
Temperature at Inflection Point, ◦ F - - 304.71 - -
T-s Slope at Inlet, E1 , lbm.R2 /BTU 8551 15,181 10,553 7178 10,330
T-s Slope at Discharge, E2 , lbm.R2 /BTU 10,744 19,708 8887 8151 15,523
T-s Slope Change, % 25.66 29.82 −15.79 13.57 50.28
T-s Curve Shape Category I I II I I
Entropy at Inflection Point, BTU/lbm.R - - - - -
Temperature at Inflection Point, ◦ F - - - - -
T-s Slope at Inlet, E1 , lbm.R2 /BTU 14,053 10,664 11,513 11,509 13,135
T-s Slope at Discharge, E2 , lbm.R2 /BTU 3773 9986 5526 5352 16,183
T-s Slope Change, % −73.15 −6.36 −52.01 −53.50 23.20
T-s Curve Shape Category III III III III I
Entropy at Inflection Point, BTU/lbm.R 0.48820 0.64410 0.55540 0.55670 -
Temperature at Inflection Point, ◦ F 288.95 176.68 286.75 294.07 -
Case Number Case 16 Case 17 Case 18 Case 19
T-s Slope at Inlet, E1 , lbm.R2 /BTU 18,783 5147 4601 2786
T-s Slope at Discharge, E2 , lbm.R2 /BTU 12,049 5762 3362 2632
T-s Slope Change, % −35.85 11.94 −26.94 −5.53
T-s Curve Shape Category II I II II
Entropy at Inflection Point, BTU/lbm.R - - - -
Temperature at Inflection Point, ◦ F - - - -
Appendix B.2. Three Point Methods

Figure A2a–e provide a summary of the polytropic efficiency deviation results compar-
ing the Taher–Evans Cubic Polynomial 2 segment method with Huntington’s three-point
method [13]. The two-segment cubic polynomial method shows sufficient accuracy for 13
of 19 example cases while Huntington’s three-point method shows sufficient accuracy for
only 7 of 19 example cases based upon the 0.001% requirement. In all example cases, the
cubic polynomial method is more accurate than Huntington’s method. Figure A2f com-
pares PTC10 CO2 and HPNG cases. Both cubic and the Huntington’s three-point method
are sufficient for the low-pressure test case but only the cubic polynomial is sufficient for
the HPNG case.
Appendix B.3. Multi-Step and Multi-Segment Methods

Table A3 lists the polytropic efficiency calculated for Taher–Evans Cubic Polynomial 2
through 10 multi-segment methods, and linear polynomial 10, 20, 50, 90 and 100 multi-step
methods. The general trend is that using more segments or steps yields higher accuracy. In
many example cases, as more segments or steps are applied, the efficiency only changes
beyond the fourth decimal place as highlighted in the table by the outlined values. The
results are shown to have “stabilized” to four decimal places for additional segments or
steps. In all example cases, the highest accuracy shown agrees between cubic and linear
methods, but the linear method requires many more steps before “stabilizing”. The general
trend of higher-pressure cases requiring more segments or steps to stabilize is confirmed.
For the supercritical ethane example case, the cubic method demonstrates a polytropic
efficiency of 79.9519% for three segments and remains at that value for the four decimal
places shown as the number of segments is increased. The linear method needs between 50
and 90 steps to reach the same accuracy for efficiency. In other words, the linear method
approaches an asymptote much more slowly than the cubic method at the cost of increased
computing resources and elapsed time to reach a solution. This illustrates the usefulness
of the cubic method in being able to take advantage of matching the slope of the T-s path
between sequential segments. Both the cubic and linear methods are approximants that
can reach essentially the same level of approximation, but the cubic method employs more
thermodynamic constraints to achieve the resulting polytropic efficiency much quicker.
It is not necessary to reach a “stabilized” polytropic efficiency to satisfy the maximum
allowable deviation requirement of 0.001% as compared to the Taher–Evans Cubic Polyno-
mial 10 segment results. Deviations are less than 0.001% for 13 of 19 example cases
Int. J. Turbomach. Propuls. Power 2021, 6, x FOR PEER REVIEW
using
31 of 45
2 cubic segments while only 4 of 19 are less than 0.001% for 10 linear steps.
Figure A2. (a–f) Summary comparison of relative deviation results for three point methods.
Figure A2. (a–f) Summary comparison of relative deviation results for three point methods.
Appendix B.3. Multi-Step and Multi-Segment Methods
Table A3 lists the polytropic efficiency calculated for Taher–Evans Cubic Polynomial 2
through 10 multi-segment methods, and linear polynomial 10, 20, 50, 90 and 100 multi-step
methods. The general trend is that using more segments or steps yields higher accuracy. In
many example cases, as more segments or steps are applied, the efficiency only changes be-
yond the fourth decimal place as highlighted in the table by the outlined values. The results
are shown to have “stabilized” to four decimal places for additional segments or steps. In all
example cases, the highest accuracy shown agrees between cubic and linear methods, but the
Table A3. Summary comparison of polytropic efficiency results for linear multi-step and cubic polynomial multi-
segment methods.
Case Number Case 1 Case 2 Case 3 Case 4 Case 5 Case 6 Case 7

Cubic Polynomial 2 Segment 75.0442 81.9199 80.6174 79.9531 59.4118 65.2211 78.4933
Linear Polynomial 10 step 75.0408 81.9194 80.6153 79.9544 59.4111 65.2189 78.4916
Case Number Case 8 Case 9 Case 10 Case 11 Case 12 Case 13 Case 14
Cubic Polynomial 2 Segment 80.3316 82.1055 72.1215 71.3151 59.2985
Linear Polynomial 10 step 80.3306 82.1073 72.1209 71.3164 59.2990
Figure A3a illustrates polytropic efficiency relative deviation results for all cubic multi-
segment methods and linear multi-step methods for the supercritical propane example
case. The more difficult an example case is, the more segments or steps it takes to achieve
a relative deviation of 0.001% or less. This most difficult case achieves this target with
5 segments for the cubic polynomial method while it would take 44 steps for the linear
polynomial method to reach the same deviation value as the 5 cubic segment method. The
cubic polynomial method does not suffer from the requirement of a small step size for each
segment as does the linear polynomial method.
Figure A3b shows the elapsed calculation times for 5 cubic segments and 44 linear
steps are approximately 3.5 s and 55.8 s, respectively. Thus, the linear method takes 16 times
longer to reach the same level of accuracy for this very difficult example case. For the six
example cases that required more than 2 cubic segments to reach a deviation less of than
0.001%, the elapsed time ratio varied from 5 to 16 times longer for the linear method.
deviation of 0.001% or less. This most difficult case achieves this target with 5 segments for
the cubic polynomial method while it would take 44 steps for the linear polynomial method
to reach the same deviation value as the 5 cubic segment method. The cubic polynomial
method does not suffer from the requirement of a small step size for each segment as does
the linear polynomial method.
Int. J. Turbomach. Propuls. Power 2021, 6, 15 Figure A3b shows the elapsed calculation times for 5 cubic segments and 44 linear steps 32 of 43
are approximately 3.5 s and 55.8 s, respectively. Thus, the linear method takes 16 times longer
to reach the same level of accuracy for this very difficult example case. For the six example
cases that required more than 2 cubic segments to reach a deviation less of than 0.001%, the
The time
elapsed elapsed
ratiotime
variedtofrom
reach a solution
5 to is dependent
16 times longer upon
for the linear the algorithm employed to
method.
resolve theelapsed
The inner and
timeouter nested,
to reach recursive
a solution loops forupon
is dependent discharge temperature
the algorithm employedand to
efficiency.
re-
solve
The the innerused
algorithm and outer nested,
for both cubicrecursive
and linearloopsmethod
for discharge
resultstemperature and efficiency.
in this appendix is identical.
The algorithm
Thus, used for
while another both cubic
algorithm and linear
might reachmethod
solutionsresults in or
faster thisslower,
appendixtheisratio
identical.
between
Thus,and
cubic while another
linear algorithm
method might
elapsed reach
times solutions
should fasterapproximately
remain or slower, the ratio
thebetween
same for cubic
a given
and linear
example case. method elapsed times should remain approximately the same for a given example
case.
(a)
(b)
Figure A3. (a) Polytropic efficiency relative deviations for supercritical propane example. (b) elapsed
Figure A3. (a) Polytropic efficiency relative deviations for supercritical propane example. (b) elapsed
time for supercritical propane example.
time for supercritical propane example.
Figures A4a–e illustrate polytropic efficiency relative deviation results for all cubic
multi-segment methods and linear multi-step methods for the 19 example cases. The 10 cubic
segment results are used as the standard for comparison. For each case, fewer cubic segments
are needed than linear steps to achieve the 0.001% required accuracy. Figure A4f directly
compares the HPNG case and the PTC10 CO2 Type 2 performance test designed to prove the
compressor should meet the design requirements within the tolerances allowed by the ASME
Figure A4a–e illustrate polytropic efficiency relative deviation results for all cubic
multi-segment methods and linear multi-step methods for the 19 example cases. The
10 cubic segment results are used as the standard for comparison. For each case, fewer
cubic segments are needed than linear steps to achieve the 0.001% required accuracy.
Figure A4f directly compares the HPNG case and the PTC10 CO2 Type 2 performance
Int. J. Turbomach. Propuls. Power 2021, 6,test designed
x FOR to prove the compressor should meet the design requirements within35the
PEER REVIEW of 45
tolerances allowed by the ASME PTC10 code [4].
Figure A4. (a–f) Summary comparison of results for linear multi-step and cubic polynomial multi-segment methods.
Figure A4. (a–f) Summary comparison of results for linear multi-step and cubic polynomial multi-segment methods.
Appendix
AppendixB.4. High
B.4. HighPressure
Pressure Ethylene
EthyleneExample
Example
This example was first created
This example was first created by by Mallen
Mallen and
andSaville [9][9]
Saville to to
demonstrate
demonstrate differences
differences be-
tween theirtheir
between polytropic work work
polytropic calculation method
calculation and Schultz’s
method method method
and Schultz’s [6] at higher pressures.
[6] at higher
Since then, itSince
pressures. has been
then,used bybeen
it has manyusedauthors as a strong
by many test
authors asof
a other
strongmethods due to
test of other its large
methods
pressure and temperature ranges. The high compression ratio and extreme
due to its large pressure and temperature ranges. The high compression ratio and extreme discharge tem-
perature for this case illustrate it would likely not be feasible to build a single section
discharge temperature for this case illustrate it would likely not be feasible to build a single com-
pressor to accomplish this application. However, its use as an example to test calculation
methods has been documented.
Ethylene has been described as an unruly fluid by Brown, [29]. Evans and Huble, [22],
postulated Brown’s comment could be due to the proximity of ethylene’s critical temperature
to usual inlet temperature ranges for common applications causing rapid, large changes in
section compressor to accomplish this application. However, its use as an example to test
calculation methods has been documented.
Ethylene has been described as an unruly fluid by Brown, [29]. Evans and Huble, [22],
postulated Brown’s comment could be due to the proximity of ethylene’s critical tem-
perature to usual inlet temperature ranges for common applications causing rapid, large
changes in specific heat, cp , and especially density with small inlet temperature changes. In
the vicinity of the critical point, large changes in inlet density can occur due to small changes
in6,inlet
Int. J. Turbomach. Propuls. Power 2021, x FORtemperature
PEER REVIEW causing a centrifugal compressor’s operating point to move36from of 45 one
side of its performance map to the other without any process operator control changes.
The Taher–Evans Cubic Polynomial method provides the ability to calculate and
graphically display the
The Taher–Evans T-s Polynomial
Cubic values at any methodpoint alongthe
provides theability
constant polytropic
to calculate efficiency
and graph-
path
ically display the T-s values at any point along the constant polytropic efficiency path demon-
between inlet and discharge. The HP ethylene example case serves to be-
strate
tweenthe method’s
inlet usefulness
and discharge. The HPtoethylene
visualize the path,
example its curvature
case serves and anthe
to demonstrate inflection
method’spoint.
Figure A5a to
usefulness shows paths
visualize thefor 1 and
path, 10 segment
its curvature andcubic polynomial
an inflection point.methods
Figure A5a and 1 and
shows 100 step
paths
linear polynomial
for 1 and 10 segment methods. The 100-step
cubic polynomial linear
methods method’s
and 1 and 100orange markers
step linear are seen
polynomial to lie on
meth-
ods.
top ofThe
the 100-step
solid bluelinear method’s cubic
10-segment orangepath.
markersTheare
oneseen to lie on
segment top of
cubic the solid
dotted blueblue
line10-
shows
ansegment cubicpoint
inflection path.near
The onethe segment
center ofcubicthe dotted
graph.blue Theline shows an inflection
10-segment point nearpath
cubic polynomial
the center
begins withofan theupward
graph. The 10-segment
concavity that cubic polynomial
transitions path begins
to downward at with an upward
the inflection con- The
point.
cavity that transitions to downward at the inflection point. The 100 step linear
100 step linear markers follow this pattern illustrating the accuracy of the cubic polynomial markers follow
this pattern illustrating the accuracy of the cubic polynomial method’s path.
method’s path.
Figure A5b shows an expanded view of the area outlined by a black box in Figure A5a.
Figure A5b shows an expanded view of the area outlined by a black box in Figure A5a.
The inflection point of the 10-segment cubic polynomial path is shown to occur at a slightly
The inflection
higher point of
temperature forthe
the10-segment
more refinedcubic polynomial
calculation. path is shown
The 10-segment cubictopolynomial
occur at a be-slightly
higher temperature for the more refined calculation. The 10-segment
gins with 5 upward concave segments, and ends with 4 downward concave segments. The cubic polynomial
begins with 5isupward
6th segment “s” shaped concave
since it segments, and ends
starts out concave withand
upward 4 downward concave segments.
shifts at the inflection point
The 6th segment
to concave is “s” shaped
downward. since itofstarts
The curvature out concave
the cubic segmentsupward and shifts
is determined at theofinflection
by values the
point to concave
derivatives of thedownward. The curvature
T-s path equation. The point of the cubic
where segments
the blue 10-segmentis determined
curve crosses bythe
values
oforange dotted linear
the derivatives ofendpoint
the T-s pathline isequation.
not an inflection point,where
The point but rather
the just
bluea 10-segment
crossing of the curve
upward concave segments of the 10-segment cubic method path.
crosses the orange dotted linear endpoint line is not an inflection point, but rather just a
crossing of the upward concave segments of the 10-segment cubic method path.
(a)
Figure A5. Cont.

(b)
FigureA5.
Figure A5.(a)
(a) T-s
T-s compression
compression paths
pathsfor
forHP
HPethylene example
ethylene examplecase. (b) the
case. (b) location of inflection
the location point point
of inflection
on T-s compression polytropic path is compared for different methods—HP ethylene example case.
on T-s compression polytropic path is compared for different methods—HP ethylene example case.
Appendix B.5. Supercritical Propane Example

The results
Appendix shown in
B.5. Supercritical TableExample
Propane A4 illustrated that the supercritical propane example
case isThetheresults
mostshown
difficult. TheA3following
in Table illustrateddiscussion uses thispropane
that the supercritical case toexample
shed light
case isupon
reasons
the mostperformance
difficult. Thecalculations are more
following discussion difficult
uses this casefortocertain applications.
shed light upon reasons Figure A6a,b
perfor-
are pressure-temperature
mance calculations are more and pressure-enthalpy
difficult diagramsFigures
for certain applications. for thisA6a,b
example. Figure A6a
are pressure-
temperature
shows and pressure-enthalpy
inlet conditions are very neardiagrams for this example.
the critical point such Figure
thatA6atheshows inlet pressure
reduced con-
andditions are very
reduced near the critical
temperature are point
very such
nearthat
1.0.the reduced
While thesepressure and reducedmight
inlet conditions tempera-
not be
ture are very near
recommended 1.0. While
in normal these inletindustry
compression conditions might not
practice, the be recommended
example serves asinannormal
“acid test”
forcompression
polytropicindustry practice,
calculation the example
methods. Figureserves as an “acid
A6b shows test” forin
isotherms polytropic
the dense calculation
phase region
methods. Figure A6b shows isotherms in the dense phase region are nearly horizontal at inlet
are nearly horizontal at inlet but are nearly vertical at discharge. For the cubic polynomial
but are nearly vertical at discharge. For the cubic polynomial compression curve shown in
compression curve shown in the figure, near the inlet, pressure changes are more rapid
the figure, near the inlet, pressure changes are more rapid than enthalpy changes while the
than enthalpy
opposite changes
is true while theconditions
near discharge opposite giving
is true the
nearpath
discharge
its shape.conditions giving the
Near discharge, onepath
itsvalue
shape.
of enthalpy could be associated with two values of pressure at the same temperature of
Near discharge, one value of enthalpy could be associated with two values
pressure
making an at EOS
the same temperature
root finding algorithm making an EOStoroot
more difficult finding algorithm more difficult
converge.
to converge.
Figures 8 and 9 in the main body of the paper illustrate the rapid variation of thermo-
Figures
dynamic 8 and 9 𝑐in the
properties 𝑋 in the
and main body of theofpaper
vicinity illustrate
the critical pointtheforrapid variation ofThis
this application. thermo-
dynamic properties
causes large c p and
differences X in
of path thebetween
slope vicinitytheof the
inletcritical point for
and discharge (𝐸 this 𝐸 ).
and application. This
causes large differences of path slope between the inlet and discharge (E1 and E2 ).
(a)
(b)
Figure
FigureA6. (a) (a)
A6. Pressure-temperature diagram
Pressure-temperature for for
diagram supercritical propane
supercritical example.
propane (b) pressure-enthalpy
example. diagram
(b) pressure-enthalpy for super-
diagram for
critical propane example.
supercritical propane example.
FigureA7
Figure A7illustrates
illustratesthe
theT-s
T-spath
pathdetermined
determined byby four
four methods
methods forfor
thethe supercritical
supercritical pro-
propane example case. The cubic and linear polynomial endpoint
pane example case. The cubic and linear polynomial endpoint methods are shownmethods are shown with
with dot-
dotted
ted lines.lines.
The The 10 segment
10 segment Taher–Evans
Taher–Evans Cubic
Cubic Polynomial
Polynomial method
method is shown
is shown with
with a solid
a solid blue
blueThe
line. line.100-step
The 100-step
linearlinear
methodmethod is shown
is shown with with markers
markers which
which laytop
lay on on of
topthe
of cubic
the cubic
poly-
polynomial line. The endpoint methods clearly do not provide a sufficiently
nomial line. The endpoint methods clearly do not provide a sufficiently accurate path for accurate paththis
for this example case. The 10-segment cubic polynomial and the 100-step linear
example case. The 10-segment cubic polynomial and the 100-step linear polynomial agree on polynomial
agree
the path. on the path.
Figure A7.Linear
FigureA7. Linearand
andcubic
cubicpolytropic
compressionmethods
methodscomparisons
comparisonsfor
forsupercritical
supercriticalpropane.
propane.
Thesupercritical
The supercritical propane
propane example
example case
case serves
serves to
to illustrate
illustratethe
theaccuracy
accuracyofofthe
the1010
seg-
seg-
ment Taher–Evans Cubic Polynomial method even when inlet conditions are near
ment Taher–Evans Cubic Polynomial method even when inlet conditions are near the critical the
critical
point of point of a fluid.
a fluid.
AppendixC.
Appendix C.Overview
Overviewof of115
115Example
ExampleCaseCaseResults
Results
Evans[3]
Evans [3]provided
provideda broad
a broad setset
of of
115115 example
example cases
cases thatthat compares
compares endpoint
endpoint and and
multi-
multi-step linear polynomial solutions to endpoint and multi-segment cubic polynomial
step linear polynomial solutions to endpoint and multi-segment cubic polynomial solutions.
solutions.
This This22
set covers setfluids
covers
(7 22 fluids
pure (7 pure
fluids and 15fluids and 15with
mixtures) mixtures)
a widewith a wide
range range
of mole of
weight
mole weight and critical point properties. The 19 example cases discussed in Appendix B
and critical point properties. The 19 example cases discussed in Appendix B above were in-
above were included in the 115 example cases. Seventeen polytropic efficiency calculation
cluded in the 115 example cases. Seventeen polytropic efficiency calculation methods were
methods were included which yielded a total of 1955 independent calculations. This effort
included which yielded a total of 1955 independent calculations. This effort is sufficient to
is sufficient to document the validity of the results for the polytropic performance methods
document the validity of the results for the polytropic performance methods employed and
employed and gives a robust comparison basis for the conclusions and recommendations
gives a robust comparison basis for the conclusions and recommendations included in this
included in this paper.
paper.
One consistent result for all 115 example cases matches the trend shown in Figure A7
One consistent result for all 115 example cases matches the trend shown in Figure A7
above in that the 10 segment cubic polynomial temperature-entropy based polytropic path
above in that the 10 segment cubic polynomial temperature-entropy based polytropic path
and the 100 step linear polynomial path are matched. If fact, these two methods agreed to
and the 100 step linear polynomial path are matched. If fact, these two methods agreed to
four decimal places for polytropic efficiency. The figures and tables below summarize the
four decimal places for polytropic efficiency. The figures and tables below summarize the
superior results provided by the cubic polynomial methods for all of the 115 example cases
superior results provided by the cubic polynomial methods for all of the 115 example cases
due to much more rapid solutions requiring fewer computing resources.
due to much
Figuremore rapid solutions
A8 illustrates requiring
that the fewer computing
cubic multi-segment and resources.
linear multi-step methods
Figure A8 illustrates that the cubic multi-segment
both approach an asymptote as the number of segments and and linear multi-step
steps increase. methods
The cubicboth
approach an asymptote
polynomial methods as the number
approach of segments
the asymptote and
much steps
more increase.
rapidly The cubic
yielding polynomial
an acceptable
methods approach the asymptote much more rapidly yielding an acceptable solution
solution in less than 10 segments while the linear polynomial methods would require more in less
thanthan
10 segments
50 steps towhile thean
provide linear polynomial
equivalent methods
solution. would
The cubic andrequire more than
linear methods 50 steps to
approach
provide an equivalent
the asymptote fromsolution.
opposite The cubic and
directions linear
for this methods
methane, approach
propane, andthe asymptote
carbon dioxidefrom
mixture
opposite examplefor
directions case.
this methane, propane, and carbon dioxide mixture example case.
Figure
FigureA8.
A8.Asymptotic
Asymptoticapproach
approach to
to aa solution.
solution.
ForFor
thethe
115115 examplecases,
example cases,114
114 reached
reached ananacceptable
acceptablerelative deviation
relative of ≤of
deviation 0.001%
≤0.001%in in
three
three or less
or less cubic
cubic polynomial
polynomial segments
segments as as documented
documented inin TableA3.
Table A4.Only
Onlyone
oneexample
examplecase
case required 5 cubic segments due to inlet conditions being extremely close to
required 5 cubic segments due to inlet conditions being extremely close to the critical the critical
point,
point, which is unlikely to occur in reality. Based upon an analysis of these results,
which is unlikely to occur in reality. Based upon an analysis of these results, the recom- the
recommended number of segments required for a particular application can be related to
mended number of segments required for a particular application can be related to the path
the path categories defined in Figure 10 above.
categories defined in Figure 10 above.
• •Category
Category I: 3 cubic polynomial segments required;
I: 3 cubic polynomial segments required;
• Category II: 5 cubic polynomial segments required;
• Category II: 5 cubic polynomial segments required;
• Category III: 5 cubic polynomial segments required.
• Category III: 5 cubic polynomial segments required.
Recall that Figure 11 shows that determining a similar set of recommendations for
the number of linear polynomial
Table A3. Cubic polynomial steps required
solution segments was
required not possible.
versus polytropicThe
pathwide scatter
curve of the
category for
linear methods
115 example cases. results that would yield similar accuracy results to the cubic methods was
very random.
Number of Cubic Polynomial Segments Required
T-s Path Category Totals
1 Segment 2 Segments 3 Segments 4 Segments 5 Segments
I 28 16 44
II 10 28 8 1 47
III 4 13 7 24
Totals 42 57 15 0 1 115
Table A4. Cubic polynomial solution segments required versus polytropic path curve category for
115 example cases.
Int. J. Turbomach. Propuls. Power 2021, 6,T-s PathPEER REVIEW Number of Cubic Polynomial Segments Required
x FOR 41 of 45
Category Totals
1 Segment 2 Segments 3 Segments 4 Segments 5 Segments
I 28 16 44
II 10 28 8 1 47
IIIRecall
that Figure
4 11 shows
that
13 determining
7 a similar set of recommendations
for the
24
number
Totals of linear polynomial
42 steps
57 required
15 was not possible.
0 The wide
1 scatter of the
115 linear
methods results that would yield similar accuracy results to the cubic methods was very ran-
dom.
Figure
FigureA9A9shows
showsdissimilar
dissimilar patterns
patterns for average elapsed
for average elapsedcalculation
calculationtime
timebetween
between cubic
cubic segments and linear steps. The cubic polynomial solutions average
segments and linear steps. The cubic polynomial solutions average time per segmenttime per seg- de-
ment decreases
creases with anwith an increasing
increasing number number of segments.
of segments. However,
However, the linear
the linear polynomial
polynomial solutions
solutions pattern shows an increasing time per step as the number of steps increases.
pattern shows an increasing time per step as the number of steps increases. The averages The are
averages are over all 115 example cases.
over all 115 example cases.
Figure
Figure A9.Average
A9. Average elapsed
elapsed calculation
calculation time
time perper cubic
cubic segment
segment andand linear
linear stepstep
for for
115115 example
example cases.
cases.
Figures
Figures A8A8andand A9A9 combined
combined show
show that
that thethe Taher–Evans
Taher–Evans Cubic
Cubic Polynomial
Polynomial polytropic
polytropic
pathmethods
path methodsare aresuperior.
superior.While
Whileboth
bothcubic
cubicandandlinear
linearpolynomial
polynomialmethods
methods can can reach
reach ac-
ceptable results, the
acceptable the cubic
cubicmethods
methodsrequire
requirefewer
fewer repetitive
repetitive iterations
iterations andand are are completed
completed
muchfaster.
much faster.This
Thisisisa aclear
clearadvantage
advantagewhen whenprocessing
processingtest
testdata
dataininreal
realtime.
time.
Figure A10 compares polytropic efficiency calculation
Figure A10 compares polytropic efficiency calculation uncertainty for uncertainty forcubic
cubicandandlin-
linear
polynomial methods. Both exhibit decreasing uncertainty with an
ear polynomial methods. Both exhibit decreasing uncertainty with an increase in cubicincrease in cubic segments
or linear or
segments steps. However,
linear the cubic polynomial
steps. However, methods uncertainty
the cubic polynomial rapidly decreases
methods uncertainty rapidly as the
number as
decreases of the
segments
number increases whileincreases
of segments the linear polynomial
while methods
the linear illustrate
polynomial a much
methods more
illus-
gradual
trate a muchdecrease
more as they approach
gradual decrease an as asymptote.
they approachThean linear methodsThe
asymptote. uncertainty asymptote
linear methods
appears to asymptote
uncertainty be more than an order
appears to of
bemagnitude
more than greater
an order than the uncertainty
of magnitude of the
greater 9-segment
than the
cubic method.
uncertainty A 9-segment
of the comparison of uncertainty
cubic method. Abetween
comparison115 and 19 examplebetween
of uncertainty case statistical
115 andpop-
19ulations
example revealed that uncertainty
case statistical populationsdecreased
revealedwith a larger
that population
uncertainty for both
decreased withcubic and lin-
a larger
ear polynomial
population methods.
for both cubic andHowever, cubic methods
linear polynomial had lower
methods. uncertainty
However, cubic than linear
methods hadmeth-
ods for
lower both population
uncertainty than linearsizes.
methods for both population sizes.
Int. J. Turbomach. Propuls. Power 2021, 6, x FOR PEER REVIEW 42
Figure A10.
Figure A10. Calculation Calculation
uncertainty uncertainty
comparison comparison
for for cases.
115 example 115 example cases.
Appendix
Appendix D. Equation ofD. Equation
State of on
Impact State Impact on
Polytropic Polytropic Efficiency
Efficiency
When aapplication
When a compression compressionisapplication is first
first designed bydesigned
a user, itbytypically
a user, it involves
typically involves
a a
cess simulator that has many equations of state available.
process simulator that has many equations of state available. One is selected from thoseOne is selected from those avai
available withinwithin a specific
a specific simulator
simulator software
software package
package based based
uponupon
the the overall
overall project
project Type. When
Type.
When various vendors receive an inquiry for a compressor containing the necessary inletand disch
ious vendors receive an inquiry for a compressor containing the necessary inlet
state point operating conditions and other pertinent details, each vendor may select a d
and discharge state point operating conditions and other pertinent details, each vendor
ent equation of state for their detailed compressor design according to individual experi
may select a different equation of state for their detailed compressor design according to
Using different equations of state can yield variations in predicted performance. The s
individual experience. Using different equations of state can yield variations in predicted
holds true for analyzing compressor test data to verify conformance to testing codes an
performance. The same holds true for analyzing compressor test data to verify conformance
contractual guarantees.
to testing codes and/or contractual guarantees.
A user is encouraged to require the use of a consistent equation of state througho
A user is encouraged
project by all to parties
requireinvolved
the use of a consistent
from inception equation of state
to final testing andthroughout
operation a to minimi
project by all parties involved from inception to final testing and operation
avoid EOS related discrepancies. It should be pointed out that variations to minimizeexist in the im
or avoid EOS related
mentationdiscrepancies.
of similarly It should
named be pointed
equations outacross
of state that variations exist inplatforms.
various software the Sh
implementationtheseof similarly
variations named
cause equations
issues, theyofshould
state across various software
be investigated platforms.
and documented.
Should these variations cause issues,
All polytropic they should
efficiency be investigated
calculation results shown andindocumented.
this Appendix D were perfor
All polytropic
withefficiency calculation
the Taher–Evans results
Cubic shown in
Polynomial this Appendix
endpoint method.DThese
were results
performedare influence
with the Taher–Evans Cubic Polynomial
thermodynamic state pointendpoint
parameters method. These
calculated by results are influenced
the various equations ofbystate. From
thermodynamicmain statetext
point parameters
of this paper, the calculated by the various
required equations equations
are (15), (16) andof state. From
(32).
the main text of thisForpaper,
giventheinputs
required equations
of fluid are (15),
composition, (16)
total and (32).
pressure and total temperature at inlet
For given inputs
discharge of fluid composition,
the parameters totaltopressure
required and total
be calculated by antemperature at inlet
EOS to determine polytropic
ciency are:
and discharge the parameters required to be calculated by an EOS to determine polytropic
efficiency are: • Enthalpy, ℎ;
• Enthalpy, h;• Entropy, 𝑠;
• Entropy, s; • Specific heat at constant pressure, 𝐶 ;
• • Compressibility C p ; 𝑋.
function,
Specific heat at constant pressure,
• Compressibility function, X.
Even though various EOS may use different reference states for enthalpy and entropy,
it is the change in these parameters that are used in Equation (32), therefore differences in
individual reference states for different EOS is not an issue for the purpose of this study.

Even though various EOS may use different reference states for enthalpy and entropy,
it is the change in these parameters that are used in Equation (32), therefore differences in
individual reference states for different EOS is not an issue for the purpose of this study.
Theequations
The equationsofofstate
state investigated
investigated are:
are:
• • PR, PR,Peng–Robinson
Peng–Robinson[30];[30];
• • LKP, Lee, Kesler, Plöcker [31,32];
LKP, Lee, Kesler, Plöcker [31,32];
• • PR-LK, Peng–Robinson with Lee–Kesler
PR-LK, Peng–Robinson with Lee–Kesler Enthalpy;
Enthalpy;
• • BWRS, Benedict–Webb–Rubin–Starling
BWRS, Benedict–Webb–Rubin–Starling [33–35] [33–35]
• • REFPROP,
REFPROP,NISTNISTreference
referencefluid
fluidthermodynamic
thermodynamic andand transport
transport properties
properties data-
database
base [28,36,37]
[28,36,37]
o# Employs Employs individual,
individual, reference
reference quality
quality EOSEOS for each
for each purepure
fluid;fluid;
o # Applies mixing rules to the Helmholtz energy
Applies mixing rules to the Helmholtz energy of the mixture of the mixture components;
components;
o# REFPROP REFPROP is a is a program,
program, not anot a database
database of measurements.
of measurements.
Results for the first four EOS above were
Results for the first four EOS above were derived from derived from Aspen
Aspen HYSYS
HYSYS V.10V.10 [27]
[27] while
while
the REFPROP results were calculated via a spreadsheet
the REFPROP results were calculated via a spreadsheet with the EOS embedded. with the EOS embedded.
FigureA11
Figure A11illustrates
illustratesvariations
variationsininpolytropic
polytropicefficiency
efficiency for
for seven
seven example
example cases
cases from
from
Appendix B using the five equations of state listed above. These
Appendix B using the five equations of state listed above. These EOS induced polytropic EOS induced polytropic
efficiencyvariations
efficiency variationscan canexceed
exceedvariation
variationdue
due toto measurement
measurement uncertainty
uncertainty values
values inin some
some
cases. The supercritical propane example case results are misleading
cases. The supercritical propane example case results are misleading in that only REFPROP in that only REFPROP
consideredthe
considered theinlet
inlet conditions
conditions toto
bebe superheated
superheated (see(see Figure
Figure A6a).
A6a). TheThe EOS EOS
usedused within
within the
the process simulator considered the fluid to be either two phase or liquid.
process simulator considered the fluid to be either two phase or liquid. This serves to illus- This serves to
illustrate
trate thethat
the care caremust
that must be exercised
be exercised in selecting
in selecting an EOSan EOS
when when
inletinlet conditions
conditions are near
are near the
the critical point
critical point of a fluid.of a fluid.
Otherthan
Other than thethe supercritical
supercritical propane
propane example
example case,case,
the the largest
largest polytropic
polytropic efficiency
efficiency var-
variation of 12.75% is seen for the MP CH 4 /C 3 H 8 /CO 2 example
iation of 12.75% is seen for the MP CH4/C3H8/CO2 example case. The figure clearly demon- case. The figure clearly
demonstrates
strates the needthe need
to use to use a consistent
a consistent approachapproach
to selection to selection of an of
of an equation equation
state forofcompres-
state for
compression
sion calculations. calculations.
Figure A11. Variations of the calculated polytropic efficiency using different equations of state for seven example cases.
Taher–Evans Cubic Polynomial endpoint method is used for the calculations.
To explore the efficiency variations shown in Figure A11, the EOS supplied parameters
that are used by the Taher–Evans Cubic Polynomial method are compared.
References
1. Taher, M. Mathematical Modeling of the Polytropic Process Using the Sequential Cubic Polynomial Approximation. In Proceed-
ings of the ASME GT2021, ASME Turbo Expo 2021, Virtual, Online, 7–11 June 2021. Paper Number GT2021-59715.
2. Taher, M.; Evans, B.F. Centrifugal Compressor Polytropic Performance Evaluation Using Cubic Polynomial Approximation for
the Temperature-Entropy Polytropic Path. In Proceedings of the ASME GT2021, ASME Turbo Expo 2021, Virtual, Online, 7–11
June 2021. Paper Number GT2021-59678.
3. Evans, B.F. Centrifugal Compressor Polytropic Efficiency—Cubic and Linear Polynomial Solution Methods Comparing 115 Example Case
Calculations; Technology Report; ResearchGate: Berlin, Germany, 2020. [CrossRef]
4. ASME PTC 10. Performance Test Code on Compressors and Exhausters; American Society of Mechanical Engineers: New York, NY,
USA, 1997.
5. Zeuner, G.A. Fundamental Laws of Thermodynamics Theory of Gases. In Technical Thermodynamics; D. Van Nostrand Company:
New York, NY, USA, 1906; Volume 1.
6. Schultz, J.M. The Polytropic Analysis of Centrifugal Compressors. ASME J. Eng. Power 1962, 84, 69–82. [CrossRef]
7. ASME PTC 10. Compressors and Exhausters Power Test Code 10; American Society of Mechanical Engineers: New York, NY,
USA, 1965.
8. ASME PTC 10. Compressors and Exhausters Power Test Code 10; reaffirmed of the 1965 edition; American Society of Mechanical
Engineers: New York, NY, USA, 1992.
9. Mallen, M.; Saville, G. Polytropic Processes in the Performance Prediction of Centrifugal Compressors; Institution of Mechanical
Engineers Conference Publications C183/77; Institution of Mechanical Engineers: London, UK, 1977; pp. 89–96.
10. Kent, R.G. Application of Basic Thermodynamics to Compressor Cycle Analysis. In Proceedings of the International Compressor
Engineering Conference, West Lafayette, IN, USA, 10–12 July 1974. Paper 135.
11. Nathoo, N.S.; Gottenberg, W.G. Measuring the Thermodynamic Performance of Multistage Compressors Operating on Mixed
Hydrocarbon Gases. In Proceedings of the Tenth Turbomachinery Symposium; Turbomachinery Laboratory, Texas A&M University:
College Station, TX, USA, 1981; pp. 15–23.
12. Nathoo, N.S.; Gottenberg, W.G. A New Look at Performance Analysis of Centrifugal Compressors Operating with Mixed
Hydrocarbon Gases. ASME J. Eng. Power 1983, 105, 920–926. [CrossRef]
13. Huntington, R.A. Evaluation of Polytropic Calculation Methods for Turbomachinery Performance. ASME Trans. J. Eng. Gas
Turbines Power 1985, 107, 872–879. [CrossRef]
14. Huntington, R. Limitations of the Schultz Calculation for Polytropic Head—A Proposal for Revision of PTC-10; Technology Report;
ResearchGate: Berlin, Germany, 1997. [CrossRef]
15. Hundseid, O.; Bakken, L.E.; Helde, T. A Revised Compressor Polytropic Performance Analysis. In Proceedings of the ASME
GT2006, ASME Turbo Expo 2006, Barcelona, Spain, 8–11 May 2006. Paper Number 91033.
16. Oldrich, J. Advanced Polytropic Calculation Method of Centrifugal Compressor. In Proceedings of the ASME IMECE2010, the
ASME International Mechanical Engineering Congress and Exposition 2010, Vancouver, BC, Canada, 12–18 November 2010.
Paper Number 40931.
17. Sandberg, M.R.; Colby, G.M. Limitations of ASME PTC 10 in Accurately Evaluating Centrifugal Compressor Thermodynamic
Performance. In Proceedings of the 42nd Turbomachinery Symposium, Houston, TX, USA, 1–3 October 2013; Turbomachinery
Laboratory, Texas A&M University: College Station, TX, USA, 2013.
18. Taher, M. ASME PTC-10 Performance Testing of Centrifugal Compressors—The Real Gas Calculation Method. In Proceedings of
the ASME GT2014, ASME Turbo Expo 2014, Düsseldorf, Germany, 16–20 June 2014. Paper Number GT2014-26411.
19. Wettstein, H.E. Polytropic Change of State Calculations. In Proceedings of the ASME IMECE2014, the ASME International
Mechanical Engineering Congress and Exposition 2014, Montreal, QC, Canada, 14–20 November 2014. Paper Number 36202.
20. Plano, M. Evaluation of Thermodynamic Models used for Wet Gas Compressor Design. Master’s Thesis, Innovative Sustainable
Energy Engineering, Norwegian University of Science and Technology, Trondheim, Norway, June 2014.
21. Evans, B.F.; Huble, S.R. Centrifugal Compressor Polytropic Performance: Consistently Accurate Results from Improved Real
Gas Calculations. In Proceedings of the ASME GT2017, ASME Turbo Expo 2017, Charlotte, NC, USA, 26–30 June 2017. Paper
Number GT2017-65235.
22. Evans, B.F.; Huble, S.R. Tutorial, ‘Centrifugal Compressor Performance: Making Enlightened Analysis Decisions. In Proceedings
of the 46th Turbomachinery Symposium, Houston, TX, USA, 11–14 December 2017; Turbomachinery Laboratory, Texas A&M
University: College Station, TX, USA, 2017.
23. Sandberg, M.R. A More Detailed Explanation of the Sandberg-Colby Method for the Evaluation of Centrifugal Compressor Thermodynamic
Performance; ResearchGate: Berlin, Germany, 2020. [CrossRef]
24. Box, G.E.; Draper, N.R. Empirical Model-Building and Response Surfaces; Wiley: Hoboken, NJ, USA, 1987.
25. IEEE/ASTM SI 10-2016—American National Standard for Metric Practice; Institute of Electrical and Electronics Engineers:
Washington, DC, USA, 2016.
26. Stepanoff, A.J. Turboblowers: Theory, Design, and Application of Centrifugal and Axial Flow Compressors and Fans; John Wiley and Sons:
New York, NY, USA, 1955.
27. HYSYS; Aspentech: Bedford, MA, USA, 2020.
28. Lemmon, E.W.; Bell, I.H.; Huber, M.L.; McLinden, M.O. NIST Standard Reference Database 23: Reference Fluid Thermodynamic and
Transport Properties-REFPROP, Version 10.0; Standard Reference Data Program; National Institute of Standards and Technology:
Gaithersburg, MD, USA, 2018.
29. Brown, R.N. Compressors Selection and Sizing, 2nd ed.; Butterworth-Heinemann: Woburn, MA, USA, 1997.
30. Peng, D.Y.; Robinson, D.B. A New Two-Constant Equation of State. Ind. Eng. Chem. Fundam. 1976, 15, 59–64. [CrossRef]
31. Lee, B.I.; Kesler, M.G. A Generalized Thermodynamic Correlation Based on Three-Parameter Corresponding States. AIChE J.
1975, 21, 510–527. [CrossRef]
32. Plöcker, U.; Knapp, H.; Prausnitz, J. Calculation of High-Pressure Vapor-Liquid Equilibria from a Corresponding States Correlation
with Emphasis on Asymmetric Mixtures. Ind. Chem. Process Des. Dev. 1978, 17, 324–332. [CrossRef]
33. Benedict, M.; Webb, G.B.; Rubin, L.C. An Empirical Equation for Thermodynamic Properties of Light Hydrocarbons and their
Mixtures. J. Chem. Phys. 1940, 8, 334–345. [CrossRef]
34. Benedict, M.; Webb, G.B.; Rubin, L.C. An Empirical Equation for Thermodynamic Properties of Light Hydrocarbons and their
Mixtures II. Mixtures of Methane, Ethane, Propane, and n-Butane. J. Chem. Phys. 1942, 10, 757–758. [CrossRef]
35. Starling, K.E. Fluid Thermodynamic Properties for Light Petroleum Systems; Gulf Publishing: Houston, TX, USA, 1973.
36. Span, R. Multiparameter Equations of State: An Accurate Source of Thermodynamic Property Data; Springer: New York, NY, USA, 2000.
37. Span, R.; Wagner, W. A New Equation of State for Carbon Dioxide Covering the Fluid Region from the Triple-Point Temperature
to 1100 K at Pressures up to 800 MPa. J. Phys. Chem. Ref. Data 1996, 25, 1509–1596. [CrossRef]

Centrifugal Compressor Polytropic Performance-Improved Rapid Calculation Results-Cubic Polynomial Methods

Uploaded by

Copyright:

Available Formats

Centrifugal Compressor Polytropic Performance-Improved Rapid Calculation Results-Cubic Polynomial Methods

Uploaded by

Document Information

Original Title

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

Centrifugal Compressor Polytropic Performance-Improved Rapid Calculation Results-Cubic Polynomial Methods

Uploaded by

Copyright:

Available Formats

International Journal of

Citation: Taher, M.; Evans, F.

Int. J. Turbomach. Propuls. Power 2021, 6, 15. https://doi.org/10.3390/ijtpp6020015 https://www.mdpi.com/journal/ijtpp

1.1. Polytropic History

1.2. Polytropic Process

work, 𝑣𝑑𝑝, which is always less thandhthe

Int. J. Turbomach. Propuls. Power 2021, 6, x FOR

2.2. Temperature-Entropy Polytropic Path Approximation: Cubic Polynomial Endpoint Method

Int. J. Turbomach. Propuls. Power 2021, 6, x FOR PEER REVIEW (3) 7 of 45

The term ( 1)ErTDS

3.1. Polytropic Efficiency Calculation: Linear Polynomial Endpoint Method

3.2. Polytropic Efficiency Calculation: Taher–Evans Cubic Polynomial Endpoint Method

Int. J. Turbomach. Propuls. Power 2021, 6, x FOR PEER

3.4. Polytropic Efficiency Calculation: Taher–Evans Cubic Polynomial Multi-Segment Method

The coefficients ( Ai , Bi , Ci , Di ), of each cubic polynomial approximant are calculated

3.4.1. Polytropic Efficiency Calculation Procedure: Taher–Evans Cubic Polynomial

(10) Repeat steps 5 through 9 for each segment.

3.4.2. Distinction between Cubic Interpolation and Approximation

Int. J. Turbomach. Propuls. Power 2021, 6, 15 18 of 43

Int. J. Turbomach. Propuls. Power 2021, 6, 15 20 of 43

FigureFigure 11. Cubic

# Cubic polynomial endpoint path methods were demonstrated to achieve better

Author Contributions: Conceptualization and methodology, M.T.; software, validation, formal

A Polynomial approximant coefficient

Appendix A. Temperature Change with Entropy along the Polytropic Path

Substituting vdP from (A5) to (A4) yields:

This can be further simplified to

Substituting ( βT − 1) with the Schultz compressibility function X gives:

Appendix B. Comparative Numeric Examples of Application of the Taher–Evans

Industry standardized compressor data sheets document operating conditions and

Table A1. Example case input data.

Propane 1.0000 1.0000 1.0000

Ethylene 1.0000 1.0000

T1 , ◦ F −10 50 98.3 110 100 100 75 98.3 −31.36 200 210

HP CH4 /C3 H8 /CO2

Methane 0.30294 0.30294 0.30294 0.256274 0.256274 0.688671 0.688671 0.8600

Ethane 0.03748 0.03748 0.03748 0.022871 0.022871 0.119957 0.119957 0.1125

Propane 0.43533 0.43533 0.43533 0.005691 0.005691 0.102964 0.102964 0.0075

isoButane 0.00222 0.00222 0.00222 0.000001 0.000001 0.031432 0.031432

Butane 0.00218 0.00218 0.00218 0.000001 0.000001 0.027541 0.027541

isoPentane 0.000001 0.000001 0.008927 0.008927

Pentane 0.000001 0.000001 0.005563 0.005563

Hexane 0.000175 0.000175

Nitrogen 0.00399 0.00399 0.00399 0.002629 0.002629 0.002477 0.002477 0.0040

p1 , psia 650 2071 2071 500 1649.7 500 1750 2520.6

p1 , bara 44.816 142.790 142.790 34.474 113.743 34.474 120.658 173.789

T1 , ◦ F 115 160 160 50 101.2 100 50 100

T1 , ◦ C 46.111 71.111 71.111 10 38.444 37.778 10 37.778

p2 , psia 2200 10,025.70 11,010.80 1700 8189.70 1800 4375 6500

p2 , bara 151.685 691.248 759.168 117.211 564.660 124.106 301.646 448.159

T2 , ◦ F 270 291.9 301.5 250 316.2 295 125 280

Appendix B.1. Endpoint Methods

Case Number Case 1 Case 2 Case 3 Case 4 Case 5

Appendix B.2. Three Point Methods

Appendix B.3. Multi-Step and Multi-Segment Methods

Case Number Case 1 Case 2 Case 3 Case 4 Case 5 Case 6 Case 7

Int. J. Turbomach. Propuls. Power 2021, 6, x FOR PEER REVIEW 34 of 45

Figure A5. Cont.

Appendix B.5. Supercritical Propane Example