https://ntrs.nasa.gov/search.jsp?

R=20180001298 2020-05-14T14:50:00+00:00Z

NASA/TM—2018-219762

Orifice Mass Flow Calculation in NASA’s

W-8 Single Stage Axial Compressor Facility
Richard F. Bozak
Glenn Research Center, Cleveland, Ohio

February 2018
 NASA STI Program . . . in Profile

Since its founding, NASA has been dedicated • CONTRACTOR REPORT. Scientific and
to the advancement of aeronautics and space science. technical findings by NASA-sponsored
The NASA Scientific and Technical Information (STI) contractors and grantees.
Program plays a key part in helping NASA maintain
this important role. • CONFERENCE PUBLICATION. Collected
papers from scientific and technical
conferences, symposia, seminars, or other
The NASA STI Program operates under the auspices
meetings sponsored or co-sponsored by NASA.
of the Agency Chief Information Officer. It collects,
organizes, provides for archiving, and disseminates • SPECIAL PUBLICATION. Scientific,
NASA’s STI. The NASA STI Program provides access technical, or historical information from
to the NASA Technical Report Server—Registered NASA programs, projects, and missions, often
(NTRS Reg) and NASA Technical Report Server— concerned with subjects having substantial
Public (NTRS) thus providing one of the largest public interest.
collections of aeronautical and space science STI in
the world. Results are published in both non-NASA • TECHNICAL TRANSLATION. English-
channels and by NASA in the NASA STI Report language translations of foreign scientific and
Series, which includes the following report types: technical material pertinent to NASA’s mission.

• TECHNICAL PUBLICATION. Reports of For more information about the NASA STI
completed research or a major significant phase program, see the following:
of research that present the results of NASA
programs and include extensive data or theoretical • Access the NASA STI program home page at
analysis. Includes compilations of significant http://www.sti.nasa.gov
scientific and technical data and information
deemed to be of continuing reference value. • E-mail your question to [email protected]
NASA counter-part of peer-reviewed formal
• Fax your question to the NASA STI
professional papers, but has less stringent
Information Desk at 757-864-6500
limitations on manuscript length and extent of
graphic presentations.
• Telephone the NASA STI Information Desk at
757-864-9658
• TECHNICAL MEMORANDUM. Scientific
and technical findings that are preliminary or of • Write to:
specialized interest, e.g., “quick-release” reports, NASA STI Program
working papers, and bibliographies that contain Mail Stop 148
minimal annotation. Does not contain extensive NASA Langley Research Center
analysis. Hampton, VA 23681-2199
NASA/TM—2018-219762

Orifice Mass Flow Calculation in NASA’s

W-8 Single Stage Axial Compressor Facility
Richard F. Bozak
Glenn Research Center, Cleveland, Ohio

National Aeronautics and

Space Administration

Glenn Research Center

Cleveland, Ohio 44135

February 2018
 Acknowledgments

The author would like to thank Jonathan Borman, TFOME/NASA, and Jeffrey Hayes, Honeywell Aerospace, for their assistance
and review of this work. This work was supported by the Advanced Air Transport Technology Project of NASA’s Advanced Air
Vehicles Program.

Trade names and trademarks are used in this report for identification
only. Their usage does not constitute an official endorsement,
either expressed or implied, by the National Aeronautics and
Space Administration.

Level of Review: This material has been technically reviewed by technical management.

Available from
NASA STI Program National Technical Information Service
Mail Stop 148 5285 Port Royal Road
NASA Langley Research Center Springfield, VA 22161
Hampton, VA 23681-2199 703-605-6000

This report is available in electronic form at http://www.sti.nasa.gov/ and http://ntrs.nasa.gov/

 Orifice Mass Flow Calculation in NASA’s W-8 Single Stage
Axial Compressor Facility
Richard F. Bozak
National Aeronautics and Space Administration
Glenn Research Center
Cleveland, Ohio 44135

Abstract
Updates to the orifice mass flow calculation for the W-8 Single Stage Axial Compressor Facility at
NASA Glenn Research Center are provided to include the effect of humidity and incorporate ISO 5167. A
methodology for including the effect of humidity into the inlet orifice mass flow calculation is provided.
Orifice mass flow calculations provided by ASME PTC-19.5-2004, ASME MFC-3M-2004, ASME Fluid
Meters, and ISO 5167 are compared for W-8’s atmospheric inlet orifice plate. Differences in expansion
factor and discharge coefficient given by these standards give a variation of about ±0.75 percent mass
flow except for a few cases. A comparison of the calculations with an inlet static pressure mass flow
correlation and a fan exit mass flow integration using test data from a 2017 turbofan rotor test in W-8
show good agreement between the inlet static pressure mass flow correlation, ISO 5167, and ASME Fluid
Meters. While W-8’s atmospheric inlet orifice plate violates the pipe diameter limit defined by each of the
standards, the ISO 5167 is chosen to be the primary orifice mass flow calculation to use in the W-8 facility.

Introduction
When evaluating the aerodynamic performance of fans and compressors, the accurate measurement of
mass flow is critical. Like many other internal flow facilities, the W-8 Single Stage Axial Compressor
Facility at NASA Glenn Research Center utilizes an orifice plate for the measurement of mass flow
through the facility. The W-8 facility’s historical orifice mass flow calculation is based on equations from
ASME Fluid Meters 6th Ed. (Refs. 1 and 2). Recently, the effect of humidity has been shown to have a
significant impact on the evaluation of the performance of a three-stage axial compressor by Berdanier,
et al. (Ref. 3). Since W-8 primarily uses an atmospheric inlet, it was decided to update W-8’s orifice flow
calculations to include the effect of humidity. Along with the orifice mass flow calculation updates to
include the effect of humidity, updated equations from ISO (Ref. 4) and ASME (Ref. 5) standards were
also included.
This paper will provide an overview of the W-8 Single Stage Axial Compressor Facility’s
atmospheric inlet orifice plate. The calculation methodology for the inclusion of the effect of humidity
and update to ISO 5167 (Ref. 4)/ASME MFC-3M-2004 (Ref. 5) are reviewed. The orifice calculations
provided in ISO 5167 and ASME MFC-3M-2004 were found to be the same and are referred to as the
ISO 5167 calculations from this point forward. In order to incorporate the effect of humidity REFPROP
Version 9.0 (Ref. 6) was utilized for the calculation of thermodynamic properties. Then, the orifice mass
flow calculation is compared against other orifice mass flow calculations from ASME PTC-19.5-2004
(Ref. 7), and ASME Fluid Meters (Ref. 2). Finally, the orifice calculations are compared against two
other flow calculation methods with an inlet static pressure mass flow correlation (Ref. 1) and an
integration fan exit flow from test data obtained during a 2017 turbofan rotor test. All calculations and
data shown are in English units as is standard for W-8 testing.

NASA/TM—2018-219762 1
 Nomenclature
𝑃𝑃 Pressure, psia
∆𝑃𝑃 Differential Pressure, psid
𝑇𝑇 Temperature, Rankine
𝐻𝐻𝑟𝑟𝑟𝑟𝑟𝑟 Relative Humidity, %
𝜌𝜌 Density, lbm/ft3
𝑑𝑑 Orifice Diameter, inches
𝐷𝐷 Pipe Diameter, inches
𝛽𝛽 Orifice-to-Pipe Diameter Ratio, 𝑑𝑑⁄𝐷𝐷
𝑙𝑙1 Distance from Orifice Upstream Face to Upstream Pressure Measurement, inches
𝑙𝑙2′ Distance from Orifice Downstream Face to Downstream Pressure Measurement, inches
𝐿𝐿1 Normalized Upstream Pressure Measurement Location, 𝑙𝑙1 ⁄𝐷𝐷
𝐿𝐿′2 Normalized Upstream Pressure Measurement Location, 𝑙𝑙2′ ⁄𝐷𝐷
𝑅𝑅𝐷𝐷 Reynold’s Number based on Pipe Diameter (𝐷𝐷)
𝜖𝜖 Expansion Factor
𝐴𝐴 Factor in Calculation of Discharge Coefficient
𝑀𝑀2′ Factor in Calculation of Discharge Coefficient
𝐶𝐶 Discharge Coefficient
𝜇𝜇 Viscosity, lbm/ft-s
𝛾𝛾 Ratio of Specific Heats
𝑔𝑔𝑐𝑐 Standard Acceleration Due to Gravity, 32.174 ft/s2
𝑅𝑅 Universal Gas Constant, 1545.34896 ft-lbf/lbm-mol-°R
𝑀𝑀 Molar Mass, g/mol
𝑥𝑥 Mole Fraction
𝑤𝑤 Mass Fraction
𝑞𝑞𝑖𝑖 Assumed Initial Mass Flow, lbm/s
𝑞𝑞𝑚𝑚 Mass Flow, lbm/s

Subscripts:
1 Upstream Orifice Condition
2 Downstream Orifice Condition
𝑣𝑣𝑣𝑣 Vapor Saturation
𝑎𝑎 Dry Air Component
ℎ Water Vapor Component

NASA/TM—2018-219762 2
 W-8 Single Stage Axial Compressor Facility
The W-8 Single-Stage Axial Compressor Facility at the NASA Glenn Research Center is an internal
flow fan/compressor facility. A schematic of the facility is shown in Figure 1. The facility is capable of
delivering up to 7,000 hp at speeds up to 21,240 rpm to a 22-in. diameter fan or compressor. Typically, up
to 100 lbm/s of air is provided from an atmospheric inlet. Flow conditioning screens in the inlet plenum
reduce the turbulence intensity at the fan to less than 1 percent. The flow through the facility is controlled
by a sleeve throttle valve. The air can be exhausted through an atmospheric exhaust system or an altitude
exhaust system. Facility instrumentation can include up to 400 channels of steady pressure and
thermocouple measurement, tip clearance sensing, and up to 96 channels of high speed rotating fan data.
Orifice plates are utilized throughout the facility to measure the mass flow through either an
atmospheric or combustion air inlet, as well as the flow through two boundary layer bleed pipes which
can be utilized to reduce the inlet boundary layer. While the same calculations can be used on each of the
orifice plates, the focus of this paper is on the atmospheric inlet orifice plates.

Figure 1.—Schematic of the W-8 Single Stage Axial Compressor Facility at the NASA Glenn Research Center.

NASA/TM—2018-219762 3
 Orifice Mass Flow Calculation
The methodology for including humidity into the ISO 5167 (Ref. 4) orifice flow calculation is shown
below. These equations utilize REFPROP (Ref. 6) for the calculation of gas properties. The humidity
compensated orifice mass flow calculation is written as a subroutine since it is used to determine the mass
flow through either the atmospheric or combustion inlets, as well as the two boundary layer bleed lines.
Pressure, temperature, and relative humidity are measured in the facility’s inlet plenum. From these
measurements, the partial pressures of air and water vapor are calculated. The partial pressure of water
vapor is the relative humidity multiplied by the pressure of vapor saturation for water at the plenum
temperature, as shown in Equation (1).

𝑃𝑃ℎ = 𝐻𝐻𝑟𝑟𝑟𝑟𝑟𝑟 ∗ 𝑃𝑃𝑣𝑣𝑣𝑣 (1)

The partial pressure of dry air is the remaining portion of pressure measured in the plenum.

𝑃𝑃𝑎𝑎 = 𝑃𝑃 − 𝑃𝑃ℎ (2)

The mole fraction of water vapor can be calculated from the partial pressure of water vapor and the
combined pressure. This fraction is assumed to be constant throughout the fluid stream since no other
fluids are added or removed.

𝑥𝑥ℎ = 𝑃𝑃ℎ ⁄𝑃𝑃 (3)

The orifice mass flow subroutine (𝑤𝑤𝑤𝑤𝑤𝑤𝑤𝑤) requires inputs that define the fluid (pressure, temperature,
mole fraction of water vapor) and define the orifice geometry (orifice diameter, pipe diameter),

𝑤𝑤𝑤𝑤𝑤𝑤𝑤𝑤(𝑃𝑃1 , ∆𝑃𝑃, 𝑇𝑇1 , 𝑥𝑥ℎ , 𝑑𝑑, 𝐷𝐷) (4)

where,

𝑃𝑃1 Upstream Orifice Pressure, psia

∆𝑃𝑃 Orifice Differential Pressure, psid
𝑇𝑇1 Upstream Orifice Temperature, °R
𝑥𝑥ℎ Mole Fraction Water Vapor
𝑑𝑑 Orifice Diameter, inches
𝐷𝐷 Pipe Diameter, inches

The partial pressures of dry air and water vapor are found from the mole fraction of water vapor (𝑥𝑥ℎ )
and the upstream orifice pressure (𝑃𝑃1 ),

𝑃𝑃1ℎ = 𝑥𝑥ℎ 𝑃𝑃1 (5)

𝑃𝑃1𝑎𝑎 = (1 − 𝑥𝑥ℎ )𝑃𝑃1 (6)

These partial pressures and the upstream orifice temperature (𝑇𝑇1 ) are then used to calculate the partial
density of each dry air and water vapor. This can be done with REFPROP or the ideal gas law (Eqs. (7)
and (8)). The partial densities are added together to obtain a mixed density of the fluid.
𝑃𝑃1𝑎𝑎
𝜌𝜌1𝑎𝑎 = 𝑅𝑅 (7)
𝑇𝑇
𝑀𝑀𝑎𝑎 1

𝑃𝑃1ℎ
𝜌𝜌1ℎ = 𝑅𝑅 (8)
𝑇𝑇
𝑀𝑀ℎ 1

NASA/TM—2018-219762 4
 𝜌𝜌1 = 𝜌𝜌1𝑎𝑎 + 𝜌𝜌1ℎ (9)

The partial densities can then be used to find the mass fractions of dry air or water vapor.

𝑤𝑤1𝑎𝑎 = 𝜌𝜌1𝑎𝑎 ⁄𝜌𝜌1 (10)

𝑤𝑤1ℎ = 𝜌𝜌1ℎ ⁄𝜌𝜌1 (11)

The ratios of specific heats for dry air and water vapor are calculated separately with REFPROP and
then combined using a mass weighting.

𝛾𝛾1 = 𝑤𝑤1𝑎𝑎 ∗ 𝛾𝛾1𝑎𝑎 + 𝑤𝑤1ℎ ∗ 𝛾𝛾1ℎ (12)

The viscosity of the dry air and water vapor mixture is found from a subroutine (𝑣𝑣𝑣𝑣𝑣𝑣𝑣𝑣𝑣𝑣𝑣𝑣𝑣𝑣(𝑃𝑃1 , 𝑇𝑇1 , 𝑥𝑥ℎ )).
This subroutine uses REFPROP to find the viscosities of water vapor and dry air at the given pressure and
temperature and mixes them based on the formulation by Tsilingiris (Ref. 7). This subroutine is provided
for Fortran and MATLAB (The Mathworks, Inc., Natick, MA) in Appendix B.
The orifice diameter is normalized by the pipe diameter for use in the expansion factor and discharge
coefficient equations.

𝛽𝛽 = 𝑑𝑑⁄𝐷𝐷 (13)

The expansion factor (𝜖𝜖), is defined in ISO 5167-2 Equation (5).

𝑃𝑃 1⁄𝛾𝛾
𝜖𝜖 = 1 − (0.351 + 0.256𝛽𝛽4 + 0.93𝛽𝛽8 ) �1 − �𝑃𝑃2 � � (14)
1

Since the calculation of Reynolds number requires a mass flow assumption, an iterative calculation is
set up to update the assumed mass flow (𝑞𝑞𝑖𝑖 ) with the calculated mass flow (𝑞𝑞𝑚𝑚 ) until the calculation error
has been minimized.
4𝑞𝑞
𝑅𝑅𝐷𝐷 = 𝜋𝜋𝜋𝜋 𝑖𝑖 (15)
𝜇𝜇
12

From Reynolds number and geometric parameters, the discharge coefficient is calculated from the
empirical equation provided in ISO 5167-2.
0.7 0.3
106 𝛽𝛽 106
𝐶𝐶 = 0.5961 + 0.0261𝛽𝛽2 − 0.216𝛽𝛽 8 + 0.000521 � � + (0.0188 + 0.0063𝐴𝐴)𝛽𝛽3.5 � 𝑅𝑅 � +
𝑅𝑅𝐷𝐷 𝐷𝐷
𝛽𝛽 4 ′ ′ 1.1
(0.043 + 0.080𝑒𝑒 −10𝐿𝐿1 − 0.123𝑒𝑒 −7𝐿𝐿1 )(1 − 0.11𝐴𝐴) 1−𝛽𝛽4 − 0.031�𝑀𝑀2 − 0.8𝑀𝑀2 �𝛽𝛽 1.3
(16)

where,
2𝐿𝐿′
𝑀𝑀2′ = 1−𝛽𝛽2 (17)

19000𝛽𝛽 0.8
𝐴𝐴 = � � (18)
𝑅𝑅𝐷𝐷

For 𝐷𝐷 and 𝐷𝐷⁄2 pressure taps, 𝐿𝐿1 = 𝑙𝑙1 ⁄𝐷𝐷 = 1 and 𝐿𝐿′2 = 𝑙𝑙2′ ⁄𝐷𝐷 = 0.47.

NASA/TM—2018-219762 5
 The mass flow is then calculated with the units specified using Equation (1) from ISO 5167-1.

1 𝜋𝜋 2𝜌𝜌1 (∆𝑃𝑃)𝑔𝑔𝑐𝑐
𝑞𝑞𝑚𝑚 = � � 𝑑𝑑2 𝐶𝐶𝐶𝐶� 1−𝛽𝛽 4
(19)
12 4

The residual error is defined as the difference between the calculated mass flow and the assumed
mass flow. Each iteration, the assumed mass flow is set equal to the calculated mass flow from the
previous iteration.

𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒 = 𝑞𝑞𝑚𝑚 − 𝑞𝑞𝑖𝑖 (20)

The calculation is considered converged when the residual calculation error is less 0.0001 lbm/s.
This calculation has been written into a Fortran subroutine for use in W-8’s Escort Data System
(Appendix A.1.) and into a MATLAB® function (Appendix A.2.).

Orifice Flow Calculation for the W-8 Atmospheric Inlet

The atmospheric inlet for the W-8 Single Stage Axial Compressor Facility is capable of using orifice
plates to measure the mass flow through the facility. The orifice plates are located between the inlet filter
housing inlet valves as shown in Figure 1. Upstream and downstream orifice pressures are measured with
static pressure taps in four circumferential locations, at axial locations 48-1/8 in. upstream and 23-3/4 in.
downstream from the center of the orifice plate. Upstream and downstream orifice fluid temperature is
measured in two circumferential locations; 246-7/8 in. from the upstream orifice face and 96-1/8 in. from
the downstream orifice face, respectively. The orifice plates are 1/2 in. thick aluminum with the diameters
given in Table 1 and the pipe is steel with an inside diameter of 47.5 in.
W-8’s atmospheric inlet orifice plates violate the pipe diameter limit of use (~39 in.) defined by all
standards evaluated. Also, without flow straighteners, the distance between the orifice plate and
downstream elbow is only about 5 pipe diameters which could result in additional uncertainty
(as specified in ISO 5167).
The constant inputs to the 𝑤𝑤𝑤𝑤𝑤𝑤𝑤𝑤 subroutine shown in Equation (4) for W-8’s atmospheric inlet are
shown below:

𝑃𝑃1 Upstream Orifice Pressure, psia

∆𝑃𝑃 Orifice Differential Pressure, psid
𝑇𝑇1 Upstream Orifice Temperature, °R
𝑥𝑥ℎ Mole Fraction Water Vapor
𝑑𝑑 Orifice Diameter from Table 1, inches
𝐷𝐷 47.5 in.

If we assume standard day inlet pressure and temperature and dry air, the orifice flow measurement
capability can be evaluated for the orifice plates available for use in W-8. Figure 2 shows the correlation
between orifice differential pressure and mass flow for the W-8 orifice plates shown in Table 1. To avoid
Mach number effects, it is specified in ISO 5167-2 (5.3.2.2) to limit measurements to 𝑃𝑃2 ⁄𝑃𝑃1 ≥ 0.75,
which for an atmospheric inlet gives a maximum recommended Δ𝑃𝑃 of about 3.5 psid. The five different
orifice plate sizes available cover a wide operating range from about 5 lbm/s to over 100 lbm/s.

NASA/TM—2018-219762 6
 TABLE 1.—W-8 ATMOSPHERIC INLET
AVAILABLE ORIFICE PLATE SIZES
Available orifice plates
Orifice plate diameter, Orifice-to-pipe
d, diameter ratio,
in. β
Designated Actual
35 34.9950 0.737
30 30.0055 0.632
25 25.0000 0.526
23 23.0038 0.484
19 18.9995 0.400
10 10.0050 0.211

120
Orifice Weight Flow, qm (lbm/sec)

100

80

60

35"
40
30"
25"
20 23"
19"
10"
0
0 0.5 1 1.5 2 2.5 3
∆P (psi)
Figure 2.—Flow capability of the W-8 atmospheric inlet orifices with standard day, dry air inlet conditions.

NASA/TM—2018-219762 7
 Comparison With Other Orifice Flow Calculations
The ISO 5167 orifice mass flow calculation is compared with the calculations provided in ASME
Fluid Meters 6th Ed., which has historically been used to calculated orifice mass flow in W-8, and ASME
PTC 19.5. The majority of the calculation is the same, but different equations for the expansion factor (𝜖𝜖)
and discharge coefficient (𝐶𝐶) are provided in each. These difference are evaluated for W-8’s atmospheric
inlet orifice plates.
The expansion factor (𝜖𝜖), used in ASME Fluid Meters 6th Ed. and ASME PTC 19.5 Equation (3-8.2) is:
∆𝑃𝑃
𝜖𝜖 = 1 − (0.41 + 0.35𝛽𝛽 4 ) (21)
𝛾𝛾𝛾𝛾1

The difference between expansion factors is shown in Figure 3 for the six W-8 atmospheric inlet
orifice plates at standard day inlet conditions. The expansion factor shown in Equation (21) differs from
that in Equation (14) (ISO 5167) by up to 0.8 percent for high ∆𝑃𝑃 conditions for all orifice plates except
the 35 in. orifice plate.
While the discharge coefficient equation given in ISO 5167 and ASME MFC-3M are identical, two
other equations are given in ASME Fluid Meters 6th Ed. and ASME PTC 19.5. The discharge coefficient
is calculated in ASME Fluid Meters 6th Ed. as follows:
K
𝐶𝐶 = (22)
�1−𝛽𝛽4

where,

𝐾𝐾 = 𝐾𝐾0 + 𝑏𝑏𝑏𝑏 (23)

0.00025
𝐾𝐾0 = (0.6014 − 0.1352𝐷𝐷−0.25 ) + (0.3760 + 0.07257𝐷𝐷−0.25 ) � + 𝛽𝛽 4 + 1.5𝛽𝛽16 � (24)
𝐷𝐷 2 𝛽𝛽 2 +0.0025𝐷𝐷

0.0011 0.0004
𝑏𝑏 = �0.0002 + �+ �0.0038 + � (𝛽𝛽 2 + (16.5 + 5𝐷𝐷)𝛽𝛽16 ) (25)
𝐷𝐷 𝐷𝐷

𝜆𝜆 = 1000⁄𝑅𝑅𝐷𝐷 (26)

The discharge coefficient given in ASME PTC 19.5 Equation (4-8.1) is:
0.0900𝐿𝐿1 𝛽𝛽 4 91.71𝛽𝛽 2.5
𝐶𝐶 = 0.5959 + 0.0312𝛽𝛽2.1 − 0.1840𝛽𝛽 8 + 1−𝛽𝛽 4
− 0.0337𝐿𝐿′2 𝛽𝛽3 + 0.75
𝑅𝑅𝐷𝐷
(27)

The discharge coefficients are compared in Figure 4 for the W-8 atmospheric inlet orifice plates with
standard day inlet conditions. The ASME PTC 19.5 discharge coefficient varies by up to 3.5 percent from
the other calculations for the larger orifice plate diameters (35 and 30 in.). The other discharge
coefficients are within about 0.8 percent of each other.
Figure 5 shows the combined effect of the differences in expansion factor and discharge coefficient
for standard day inlet conditions for W-8’s atmospheric inlet. For the two larger orifice diameters (35 and
30 in.), the ASME PTC 19.5 discharge coefficient difference resulted in a difference between 1 and
3.5 percent. The remaining mass flows for each of the standards are within about ±0.75 percent of the
ISO 5167 calculation.

NASA/TM—2018-219762 8
 1

0.99

0.98
Expansion Factor

0.97 35"
0.96 30"
25"
0.95 23"
19"
0.94
10"
0.93 ISO 5167
ASME Fluid Meters
0.92
0 0.5 1 1.5 2 2.5 3
∆ P (psi)

0.2
Expansion Factor Difference (%)

0
relative to ISO 5167

-0.2

-0.4

-0.6

-0.8
0 0.5 1 1.5 2 2.5 3
∆ P (psi)
Figure 3.—Comparison of expansion factor between ISO 5167/ASME MFC-3M and ASME Fluid Meters
6th Ed./ASME PTC 19.5 for W-8’s atmospheric inlet orifice plates.

NASA/TM—2018-219762 9
 0.63

0.625 35"
30"
Discharge Coefficient

0.62
25"
0.615 23"
19"
0.61 10"
0.605 ISO 5167
ASME PTC 19.5
0.6 ASME Fluid Meters
0.595

0.59
0 0.5 1 1.5 2 2.5 3
∆ P (psi)

3.5
Discharge Coefficient Difference (%)

2.5
relative to ISO 5167

1.5

0.5

-0.5
0 0.5 1 1.5 2 2.5 3
∆ P (psi)

Figure 4.—Comparison of orifice discharge coefficients between ISO 5167/ASME MFC-3M, ASME Fluid Meters
6th Ed., and ASME PTC 19.5 for W-8’s atmospheric inlet orifice plates.

NASA/TM—2018-219762 10
 100

80
Mass Flow (lb /s)

35"
m

30"
60 25"
23"
40 19"
10"
20 ISO 5167
ASME PTC 19.5
ASME Fluid Meters
0
0 0.5 1 1.5 2 2.5 3
∆ P (psi)

4
Mass Flow Difference (%)

3
relative to ISO 5167

-1
0 0.5 1 1.5 2 2.5 3
∆ P (psi)

Figure 5.—Comparison of calculated mass flow between ISO 5167/ASME MFC-3M, ASME Fluid Meters 6th Ed.,
and ASME PTC 19.5 for W-8’s atmospheric inlet orifice plates.

NASA/TM—2018-219762 11
 Test Data Comparison
The orifice flow is compared against two other calculation methods. The first being an inlet static
pressure correlation described in Appendix A of Reference 1. In this correlation, a ring of static pressures
in the inlet are correlated with mass flow using a curve-fit provided by CFD simulations of the W-8 inlet
bellmouth. This correlation was shown to be within 0.4 percent of the ASME Fluid Meters 6th Ed. orifice
flow calculation at the highest flow rate (95 lbm/s). The second calculation of mass flow utilizes three
seven-element total pressure and temperature rakes behind a rotor and static pressures measured on the
end-walls in the same plane. From these measurements, Mach number and density are calculated.
Assuming an average flow angle allows for the calculation of axial velocity and mass flow. Due to the
constant flow angle assumption and nonuniform flow, this calculation is expected to have the most error.
In February 2017 data was obtained with a 1.5 fan pressure ratio high bypass turbofan rotor installed
in the W-8 facility. Data was obtained with the 35 in. orifice plate installed and the mass flow calculations
are compared at nine conditions across the operating range in Figure 6. At the time this data was obtained,
the downstream pressure taps were installed 30-1/8 in. from the downstream orifice face (D/2), which
violates the location limits provided in ISO 5167 by about 6 in. This measurement location has since been
moved to 23-3/4 in. from the orifice plate, which is within the limits. This change was found to have an
impact of about 0.15 percent on the measured mass flow. The orifice calculations are compared relative to
the ISO 5167 calculation in terms of percent error. The fan exit mass flow calculation over-predicts the
mass flow by up to 5 percent. The ASME PTC 19.5 calculation also over-predicts the mass flow by about
3 percent. The static pressure mass flow correlation and orifice calculation from ASME Fluid Meters have
between 0.5 and 1 percent error relative to the ISO 5167 calculation.

4
Weight Flow Difference (%)
relative to ISO 5167

ASME PTC 19.5

2 ASME Fluid Meters
Static Pressure Correlation
Fan Exit Calculation
1

0
40 50 60 70 80 90 100
Orifice Mass Flow [ISO 5167] (lb /s)
m
Figure 6.—Comparison of orifice calculations with inlet static pressure correlation and
fan exit weight flow calculation relative to the ISO 5167 orifice calculation.

NASA/TM—2018-219762 12
 Conclusions
The W-8 facility’s orifice mass flow calculations have been updated to include humidity and follow a
standard (ISO 5167). A methodology for including the effect of humidity into the inlet orifice mass flow
calculation was provided. Orifice mass flow calculations provided by ASME standards ASME PTC-19.5-
2004 and ASME MFC-3M-2004, ASME Fluid Meters, and ISO 5167 were compared for W-8’s
atmospheric inlet orifice plate. Differences in expansion factor and discharge coefficient give a variation
of about ±0.75 percent mass flow except for the ASME PTC 19.5 discharge coefficient for the larger
orifice diameters (35 and 30 in.) which gave between 1 and 3.5 percent difference. The calculations were
also compared with an inlet static pressure mass flow correlation and a fan exit mass flow integration
using test data from a 2017 turbofan rotor test in W-8. For this case, the ISO 5167 calculation is within
about 0.5 percent of the historic W-8 mass flow calculation (ASME Fluid Meters) and within about
1 percent of the inlet static pressure mass flow correlation. While W-8’s atmospheric inlet orifice plate
violates the pipe diameter limit of use defined by the standards each of the standards, the ISO 5167 was
chosen to be the primary orifice mass flow calculation to use in the W-8 facility.

NASA/TM—2018-219762 13
 Appendix A.—Orifice Mass Flow Subroutines
A.1. ISO 5167 Orifice Mass Flow FORTRAN/ESCORT Subroutine
C-----------------------------------------------------------------------
C WFOR.FOR
C Subroutine to calculate Orifice Weight Flow (WFO) in lbm/sec.
C Called from FACILITY_CALCS.TAG.
C
C Modifications:
C 08/16/2017 – Updated routine to follow ISO 5167. Removed
C correction for static tap locations, thermal expansion.
C Assumed Reynold’s # diameter changed to D from ORPLATE.
C 08/02/2016 - Updated routine to account for actual upstream and
C downstream tap locations, for thermal expansion, and
C for discharge coefficient based on
C ASME PTC-19.5-2004. (J. Borman and R. Bozak)
C 10/03/2016 - Add convergence error, IERR, to calling sequence.
C (Jerri Vokac for Rick Bozak)
C 09/14/2016 - Initial Version submitted by Rick Bozak and modified
C for Escort system. (Jerri Vokac)
C
C Call: WFOR( UPO, DPOAV, TOAV, XH2O, ORFD, PIPED, WFO, WFORERR )
C
C Inputs:
C UPO: Upstream Orifice Pressure (psia)
C DPOAV: Orifice Delta Pressure (psid)
C TOAV: Orifice Temperature (degR)
C XH2O: Mole Fraction of Water Vapor
C ORFD: Orifice Plate Diameter (inch)
C PIPED: Pipe Diameter (inch)
C
C Output:
C WFO: Orifice Weight Flow (lbm/sec)
C WFORERR: Convergence Error:
C - 0 for successful
C - 1 solution not converged
C-----------------------------------------------------------------------
C-----------------------------------------------------------------------

SUBROUTINE WFOR( UPO, DPOAV, TOAV, XH2O, ORPLATE, D,

* WFO, WFORERR )

C-----------------------------------------------------------------------
C Data Recording Parameters
C-----------------------------------------------------------------------
COMMON/QAREC/DONT_USE_RDGN,BAROM,IRDTY,IRECB,IABO,IREC,IESP(2),
* IERLOG,IRCS,IDONT_USE,ISTORCYC,ICYCTY,ICYC,ICCS,
* IBATCH,NETACT,IAPREC,IAPCLR,IAPRCS,IPRCS,IAPSCAN,
* NAVGSCAN,ISTDAVG,ITRANSFL,NOREC_CAL,ISGHEAD(5),
* IDATAVIS,IVIDEOREC,LRNT,IREADNO_MASS,
* IDRTRANS,ICOMFILE,IREADNO_DUMPS,ISSA
INTEGER*2 IRDTY,IRECB,IABO,IREC,IESP,IERLOG,IRCS,IDONT_USE,
* ISTORCYC,ICYCTY,ICYC,ICCS,IBATCH,NETACT,IAPREC,
* IAPCLR,IAPRCS,IPRCS,IAPSCAN,NAVGSCAN,ISTDAVG,
* ITRANSFL,NOREC_CAL,ISGHEAD,IDATAVIS,IVIDEOREC,LRNT,
* IDRTRANS,ICOMFILE
C-----------------------------------------------------------------------
C Include Standard Constants
C-----------------------------------------------------------------------
INCLUDE 'CONSTANTS.FOR'

C-----------------------------------------------------------------------

NASA/TM—2018-219762 15
C Declare Local Airflow Variables
C-----------------------------------------------------------------------
real Md,ko,l,k
real PHP,PAP,Mu
real RHOH, RHOH2, RHOA, RHOA2, RHO2, RHOM, WAIR, WH2O
real GAMH2O, GAMAIR, GAM, XGAM, Y, CD, RED
real UPO, DPOAV, TOAV, XH2O
real WFO, MDOT

C-----------------------------------------------------------------------
C Declare Local Dimensioning Variables
C-----------------------------------------------------------------------
real L1, L2, L1ND, L2ND, BETA, D, ORPLATE, PIPED, ORFD
real ALPHAP, ALPHAO

C-----------------------------------------------------------------------
C Declare Local REFPROP Variables
C-----------------------------------------------------------------------
real*8 wmm, ttrp, tnbpt, tc, pc, Dc, Zc, acf, dip, Rgas
real*8 p, rhol, rhov, xliq, xvap
real*8 rhopa, rhoph, rho, eta, tcx, e, h, s, cv, cp, w, hjt
real*8 x(2)
real*8 TPROP, PPROP, PPROPA, PPROPH, PPROPA2, PPROPH2, RHOPROP

parameter (ncmax=20) ! max number of components in mixture

character hrf*3, herr*255
character*255 hf(ncmax), hfmix

C REFPROP SETUP initialization (cannot run 1/second)

c IF ( IBATCH .NE. 0 ) THEN
c hf(1) = 'AIR.PPF'
c hf(2) = 'WATER.FLD'
c hrf = 'DEF'
c hfmix = 'hmx.bnc'
c call SETUP( 2, hf, hfmix, hrf, ierr, herr )
c if (ierr.ne.0) write(*,*) herr
c ENDIF

x(1) = 0.5000d0 !! Air composition (molar fractions)

x(2) = 0.5000d0

C Initialize error term

ierr = 0

C-----------------------------------------------------------------------
C REFPROP Inputs
C-----------------------------------------------------------------------
PPROP = UPO / KPA2PSI
TPROP = TOAV / K2R

C Partial Pressure of Water Vapor and Air

PHP = XH2O * UPO
PAP = (1 - XH2O) * UPO

PPROPH = PHP/KPA2PSI
PPROPA = PAP/KPA2PSI
PPROPH2 = (XH2O*(UPO-DPOAV))/KPA2PSI
PPROPA2 = ((1-XH2O)*(UPO-DPOAV))/KPA2PSI

C Air Properties
CALL PUREFLD(1)
call info(1, wmm, ttrp, tnbpt, tc, pc, Dc, Zc, acf, dip, Rgas)

NASA/TM—2018-219762 16
 Mair = wmm !! Molecular Weight, g/mol
Rair = R / Mair !! Gas Constant for dry air, lbf*ft/lbm*R*mol

call tprho ( TPROP, PPROPA, x, 2, 0, rhopa, ierr, herr)

rhoa = rhopa * Mair * KGM2LBFT !! Density lb/ft^3

C ADDED 08/02/2017-----------------------------------------------------
CALL TPRHO ( TPROP, PPROPA2,X, 2, 0, RHOPA, IERR, HERR)
RHOA2= RHOPA * MAIR * KGM2LBFT
C----------------------------------------------------------------------

RHOPROP = RHOA/(Mair*KGM2LBFT)
call therm ( TPROP, RHOPROP, x, PPROPA,
* e, h, s, cv, cp, w, hjt)
gamair = cp / cv !! Ratio of Specific Heats

C Water Vapor Properties

CALL PUREFLD(2)
call info ( 2, wmm, ttrp, tnbpt, tc, pc, Dc, Zc, acf,
* dip, Rgas)

Mh2o = wmm !! Molecular Weight, g/mol

Rh2o = R / Mh2o !! Gas Constant for water vapor,
!! lbf*ft/lbm*R*mol

call tprho ( TPROP, PPROPH, x, 2, 0, rhoph, ierr, herr)

rhoh = rhoph * Mh2o * KGM2LBFT !! Density lb/ft^3

C ADDED 08/02/2017----------------------------------------------------
CALL TPRHO (TPROP, PPROPH2, X, 2, 0, RHOPH, IERR, HERR)
RHOH2= RHOPH * MH2O * KGM2LBFT
C---------------------------------------------------------------------

RHOPROP = RHOH/(MH2O*KGM2LBFT)
call therm ( TPROP, RHOPROP, x, PPROPH,
* e, h, s, cv, cp, w, hjt)
gamh2o = cp / cv !! Ratio of Specific Heats

C Calculation of Mixture Properties

rhom = rhoa + rhoh
RHO2 = RHOA2+ RHOH2
wair = rhoa / rhom
wh2o = rhoh / rhom

gam = (gamair*wair) + (gamh2o*wh2o)

M = 1 / ((wair/Mair)+(wh2o/Mh2o))
Rm = R/M

C Subroutine to calculate the viscosity of the mixture

call viscmix(UPO,TOAV,XH2O,Mu)

C Calculated Properties
beta = ORPLATE / D !! Orifice Diameter/ Pipe ID
L1ND = 1 !! Non-dim length from plate to upstream tap
L2ND = 0.47 !! Non-dim length from plate to dwnstm tap

C Assumed Initial Mass Flow

Md = 50.
C Expansion Factor Air (Updated 08/02/2017) [5]
Y=1-((0.351+0.256*(BETA**4)+0.93*(BETA**8))*
* (1-(((UPO-DPOAV)/UPO)**(1/GAM))))

NASA/TM—2018-219762 17
C Initial Error
e = 1.
icount = 0

do while ( e .gt. 0.0001 )

icount = icount + 1

C Reynolds Number Air (based on D)

Red = ( 4. * Md ) / ( pi * ( D / 12. ) * Mu )

C xgam = ( DPOAV / UPO ) / gam

C Discharge Coefficient (Updated 08/02/2017) [4]

A =(19000*beta/Red)**0.8
M2 =(2*L2ND)/(1-beta)

C=0.5961+0.0261*(beta**2)-0.216*(beta**8)
* +0.000521*(((10**6)*beta/Red)**0.7)
* +((0.0188+0.0063*A)*(beta**3.5)*(((10**6)/Red)**0.3))
* +((0.043+0.080*exp(-10*L1ND)-0.123*exp(-7*L1ND))
* *(1-0.11*A)*((beta**4)/(1-(beta**4))))
* -(0.031*(M2-0.8*(M2**1.1))*(beta**1.3))

C Mass Flow Calculations

mdot = pi/48.0*(ORPLATE**2)*CD*Y
* *(2*GC*rhom*DPOAV/(1-beta**4))**0.5

C Calculate Error
e = abs( mdot - Md )

C Set Mass flow estimate for next iteration

Md = mdot

if ( icount .gt. 1000 ) then

WFORERR = 1
C type *,'***** Orifice Mass Flow Not Converged *****'
exit
endif

enddo

C Set Orifice Mass Flow Output

wfo = mdot

C------------------------------------------------------------------------------
return
end

NASA/TM—2018-219762 18
 A.2. ISO 5167 Orifice Mass Flow MATLAB® Subroutine

%% Orifice Weight Flow Calculation

% 8/15/2017: Rick Bozak
%
% Iteratively solves for orifice coefficients and then calculates
% a weight flow for the mixture of dry air and water vapor following
% ISO 5167 and ASME MFC-3M-2004. Uses viscmix subroutine.
% Uses REFPROP for fluid properties.
%
% Requires SETUP of two fluids; (1) 'air.ppf' and (2) 'water.fld'
% For dry air, XH2O = 0.
%
% Usage:
% [WFO IERR] = wfor_iso(P1,DP,T,XH2O,d,D)
%
% Inputs:
% P1: Orifice Upstream Pressure (psia)
% DP: Orifice Delta Pressure (psi)
% T: Orifice Temperature (deg R)
% XH2O: Mole Fraction of Water Vapor
% d: Orifice Diameter (inches)
% D: Pipe Inside Diameter (inches)
%
% Output:
% WFO: Mass Flow of Fluid (lbm/s)
% IERR: Convergence Error, 0 for solved, 1 not converged

%% Test Case
% Inputs
%P1=14.5; % psia
%DP=0.5; % psi
%T=534.39; % R
%XH2O=0.01936; % Mole Fraction of water vapor
%d=35; % Orifice Diameter, inches
%D=47.5; % Inside Pipe Diameter, inches

% Output
% WFO = 87.6443; % lbm/sec
%
% 87.64423 = wfor_iso(14.5,0.5,534.39,0.01936,35,47.5)

%% Function
function [wfo C epsilon IERR] = wfor_iso(P1,DP,T1,XH2O,d,D)

IERR=0;

% Unit Conversion Factors

K2R=1.8;
kPa2psi=0.145038;
M2FT=3.280084;
kgm2lbft=0.062428;
pas2lfts=0.671969;

% Air wand Water Vapor Properties

R=1545.34896; % Universal Gas Constant, lbm-ft/lbf*R
GC=32.1740486; % Standard Acceleration Due to Gravity

% Pressure and Temperature in SI

PPROP=P1/kPa2psi;
TPROP=T1/K2R;

NASA/TM—2018-219762 19
% Partial Pressure of Water Vapor
Php=XH2O*P1;

% Partial Pressure of Air

Pap=(1-XH2O)*P1;

PPROPH=Php/kPa2psi;
PPROPA=Pap/kPa2psi;

% Air
Ma=refpropm('M','T',TPROP,'P',PPROPA,'air.ppf'); % Molecular Weight, g/mol
rhoa=refpropm('D','T',TPROP,'P',PPROPA,'air.ppf')*kgm2lbft; % lbm/ft^3
gama=refpropm('K','T',TPROP,'P',PPROPA,'air.ppf'); % Ratio of specific heats

% Water Vapor
Mh=refpropm('M','T',TPROP,'P',PPROPH,'water'); % Molecular Weight, g/mol
rhoh=refpropm('D','T',TPROP,'P',PPROPH,'water')*kgm2lbft; % lbm/ft^3
gamh=refpropm('K','T',TPROP,'P',PPROPH,'water'); % Ratio of specific heats

% Calculate Mass Fraction of Water Vapor (Specific Humidity)

wa=rhoa./(rhoa+rhoh);
wh=rhoh./(rhoa+rhoh);

% Calculation of Mixture Properties

rho=rhoa+rhoh;
Mu=viscmix(P1,T1,XH2O);
gam=(gama*wa)+(gamh*wh);
M=1/((wa/Ma)+(wh/Mh));
Rm=R/M;

% Orifice Plate Properties

% Calculated Properties
beta=d/D; % Orifice Diameter/ Pipe ID
L1=1;
L2=0.47;

% Assumed Initial Mass Flow

qi=50;

% Initial Error
err=1;
count=0;

% Iterative Calculation of Orifice Mass Flow

while (err > 0.0001)
count=count+1;
Red=(4*qi)/(pi*(D/12)*Mu); % Reynolds Number (based on d)

% Expansion Factor [ISO 5]

epsilon=1-((0.351+0.256*(beta^4)+0.93*(beta^8))*(1-(((P1-DP)/P1)^(1/gam))));

% Discharge Coefficient [ISO 4]

A=(19000*beta/Red)^0.8;
M2=(2*L2)/(1-beta);
C=0.5961+0.0261*(beta^2)-0.216*(beta^8)...
+0.000521*(((10^6)*beta/Red)^0.7)...
+((0.0188+0.0063*A)*(beta^3.5)*(((10^6)/Red)^0.3))...
+((0.043+0.080*exp(-10*L1)-0.123*exp(-7*L1))*(1-0.11*A)*((beta^4)/(1-
(beta^4))))...
-(0.031*(M2-0.8*(M2^1.1))*(beta^1.3));

NASA/TM—2018-219762 20
 % Mass Flow [2-1]
qm=(1/12)*(pi/4)*(d^2)*C*epsilon*...
sqrt((2*rho*DP*GC)/(1-(beta^4)));

% Calculate Error
err=abs(qm-Md);

% Break Loop if Mass Flow Does Not Converge

if count>1000
disp('Orifice Mass Flow Not Converged')
IERR=1;
break
end

% Set Mass flow estimate for next iteration

qi=qm;

end

% Set Mass Flow for Output

wfo=qm;

NASA/TM—2018-219762 21
 Appendix B.—Viscosity Mixture Subroutine
B.1. Viscosity Mixture FORTRAN/ESCORT Subroutine

C-----------------------------------------------------------------------
C VISCMIX.FOR
C Subroutine to calculate the viscosity of a mixture of dry air
C and water vapor.
C Called from WFOR.FOR.
C
C Modifications:
C 11/01/2016 - Initial Version
C ( Rick Bozak)
C
C Call: VISCMIX( P, T, XH2O , Mu)
C
C Inputs:
C P: Pressure (PSIA)
C T: Temperature (R)
C XH2O: Mole Fration of Water Vapor
C
C Output:
C Mu: Viscosity (lb-ft/s)
C
C-----------------------------------------------------------------------
C-----------------------------------------------------------------------

SUBROUTINE VISCMIX( Pdum, T, XH2O, Mu )

real Mair,Mh2o,PHP,PAP,Mu,PHIAV,PHIVA,Muair,Muh2o
real*8 wmm, ttrp, tnbpt, tc, pc, Dc, Zc, acf, dip, Rgas
real*8 p, rhol, rhov, xliq, xvap
real*8 rho, eta, tcx, e, h, s, cv, cp, w, hjt
real*8 x(2)

real*8 TPROP, PPROP, PPROPA, PPROPH, RHOPROP

real*8 k2r, kpa2psi, m2ft, kgm2lbft, pas2lfts

character*255 hrf*3, herr*255

character*255 hfmix

data pi/3.14159/

C-----------------------------------------------------------------------
c For REFPROP

x(1) = 0.5000d0 !! Air composition (molar fractions)

x(2) = 0.5000d0

c-----------------------------------------------------------------------

ierr = 0

C Temperature Conversion Kelvin to Rankine

K2R = 1.8

C Pressure Conversion kPa to psi

KPA2PSI = 0.145038

C Length Conversion m to ft

NASA/TM—2018-219762 23
 M2FT = 3.28084

C Density Conversion kg/m^3 to lbm/ft^3

KGM2LBFT = 0.062428

C Viscosity Conversion Pa-s to lbm/ft-s

PAS2LFTS = 0.671969

C Air Wand Water Vapor Properties

R = 1545.33 !! Universal Gas Constant, lbm-ft/lbf*R
GC = 32.174 !! Standard Acceleration Due to Gravity

C REFPROP Inputs
p = pdum !! jav
PPROP = P / KPA2PSI
TPROP = T / K2R

C Partial Pressure of Water Vapor and Air

PHP = XH2O * P
PAP = (1 - XH2O) * P

PPROPH = PHP/KPA2PSI
PPROPA = PAP/KPA2PSI

C Air Properties
CALL PUREFLD(1)
call info( 1, wmm, ttrp, tnbpt, tc, pc, Dc, Zc, acf, dip, Rgas)
Mair = wmm !! Molecular Weight, g/mol
Rair = R / Mair !! Gas Constant for dry air, lbf*ft/lbm*R*mol

call tprho ( TPROP, PPROPA, x, 2, 0, rho, ierr, herr)

rhoa = rho * Mair * KGM2LBFT !! Density lb/ft^3

RHOPROP = rhoa /(Mair*KGM2LBFT)

call trnprp (TPROP,RHOPROP,x,eta,tcx,ierr,herr)
Muair = eta * 1e-6 * PAS2LFTS

C Water Vapor Properties

CALL PUREFLD(2)
call info ( 2, wmm, ttrp, tnbpt, tc, pc, Dc, Zc, acf,
* dip, Rgas)

Mh2o = wmm !! Molecular Weight, g/mol

Rh2o = R / Mh2o !! Gas Constant for water vapor,
!! lbf*ft/lbm*R*mol

call tprho ( TPROP, PPROPH, x, 2, 0, rho, ierr, herr)

rhoh = rho * Mh2o * KGM2LBFT !! Density lb/ft^3

RHOPROP = RHOH/(MH2O*KGM2LBFT)
call trnprp (TPROP,RHOPROP,x,eta,tcx,ierr,herr)
Muh2o = eta * 1e-6 * PAS2LFTS

C Mixture Properties
rho = rhoa + rhoh
wair = rhoa / rho
wh2o = rhoh / rho

C Equation 22
PHIAV = ((2**0.5)/4)*((1+(Mair/Mh2o))**-0.5)*
& ((1+((Muair/Muh2o)**0.5)*((Mh2o/Mair)**0.25))**2)

C Equation 23

NASA/TM—2018-219762 24
 PHIVA = ((2**0.5)/4)*((1+(Mh2o/Mair))**-0.5)*
& ((1+((Muair/Muh2o)**0.5)*((Mair/Mh2o)**0.25))**2)

C Equation 20 to Combine
Mu = ((wair*Muair)/(wair+(wh2o*PHIAV)))+
& ((wh2o*Muh2o)/(wh2o+(wair*PHIVA)))

C------------------------------------------------------------------------------
return
end

NASA/TM—2018-219762 25
 B.2. Viscosity Mixture MATLAB® Subroutine
% Function for Calculating the Viscosity of an Air/Water Vapor Mixture
% Rick Bozak: 9/16/2016
%
% Follows Tsilingiris, P.T "Thermophysical and transport properties
% of humid air at temperature range between 0 and 100 C", Energy Conversion
% and Management, Vol. 45, Issue 5, p. 1098-1110, May 2008.
%
%
%% Test Case
% Inputs
%P=14.5; % psia
%T=534.39; % R
%XH=0.01936; % Mole Fraction
%
% Output
% Mum = 1.2274e-05
% 1.2279e-05 = viscmix(14.5,534.39,0.01936)

function [Mum]=viscmix(P,T,XH)

% Unit Conversion Factors

K2R=1.8;
kPa2psi=0.145038;
M2FT=3.280084;
kgm2lbft=0.062428;
pas2lfts=0.671969;

% Pressure and Temperature in SI

PPROP=P/kPa2psi;
TPROP=T/K2R;

% Partial Pressure of Water Vapor

Php=XH*P;

% Partial Pressure of Air

Pap=(1-XH)*P;

PPROPA=Pap/kPa2psi;
PPROPH=Php/kPa2psi;

% Air
Ma=refpropm('M','T',TPROP,'P',PPROPA,'air.ppf'); % Molecular Weight, g/mol
Mua= refpropm('V','T',TPROP,'P',PPROPA,'air.ppf')*pas2lfts; % Viscosity, lbm/ft-s
rhoa=refpropm('D','T',TPROP,'P',PPROPA,'air.ppf')*kgm2lbft; % lbm/ft^3

% Water Vapor
Mh=refpropm('M','T',TPROP,'P',PPROPH,'water'); % Molecular Weight, g/mol
Muh=refpropm('V','T',TPROP,'P',PPROPH,'water')*pas2lfts; % Viscosity, lbm/ft-s
rhoh=refpropm('D','T',TPROP,'P',PPROPH,'water')*kgm2lbft; % lbm/ft^3

rho=rhoa+rhoh;

% Calculate Mass Fraction of Water Vapor (Specific Humidity)

wa=rhoa./rho;
wh=rhoh./rho;

% Equation 22
phiav=(sqrt(2)/4)*((1+(Ma/Mh))^(-1/2))*...
((1+((Mua/Muh)^(1/2))*((Mh/Ma)^(1/4)))^2);

NASA/TM—2018-219762 26
% Equation 23
phiva=(sqrt(2)/4)*((1+(Mh/Ma))^(-1/2))*...
((1+((Muh/Mua)^(1/2))*((Ma/Mh)^(1/4)))^2);

% Equation 20, Viscosity of the Mixture

Mumix=((wa*Mua)/(wa+(wh*phiav)))+((wh*Muh)/(wh+(wa*phiva)));

Mum=Mumix;

NASA/TM—2018-219762 27
 References
1. Van Zante, Dale E., Podboy, Gary G., Miller, Christopher J., Thorp, Scott A., “Testing and
Performance Verification of a High Bypass Ratio Turbofan Rotor in an Internal Flow Component
Test Facility,” GT2007-27246 and NASA/TM—2009-215661.
2. Beam, Howard, editor, “Fluid Meters 6th Edition,” ASME, New York, 1971.
3. Berdanier, Reid A., Smith, Natalie R., Fabian, John C., and Key, Nicole L., “Humidity Effects on
Experimental Compressor Performance – Corrected Conditions for Real Gases,” Journal of
Turbomachinery, March 2015, Vol. 137, 031011-1.
4. International Organization for Standardization “Measurement of Fluid Flow by Means of Pressure
Differential Devices Inserted in Circular-Cross Section Conduits Running Full,” ISO 5167 (2003).
5. ASME Measurement of Fluid Flow in Closed Conduits “Measurement of Fluid Flow in Pipes Using
Orifice, Nozzle, and Venturi,” ASME MFC-3M-2004.
6. Lemmon, E.W., Huber, M.L., McLinden, M.O. NIST Standard Reference Database 23: Reference
Fluid Thermodynamic and Transport Properties-REFPROP, Version 9.1, National Institute of
Standards and Technology, Standard Reference Data Program, Gaithersburg, 2013.
7. ASME Performance Test Code “Flow Measurement,” ASME PTC 19.5-2004. Reaffirmed 2013.
8. Tsilingiris, P.T “Thermophysical and transport properties of humid air at temperature range between
0 and 100 C,” Energy Conversion and Management, Vol. 45, Issue 5, p. 1098–1110, May 2008.

NASA/TM—2018-219762 28