<<
Home , Gas

NATL INST. OF STAND & TECH

AlllDS T7M3flS gg I^EFEREIMCE Publi- cations

I/) Q NBS TECHNICAL NOTE 1186 *

U.S. DEPARTMENT OF COMMERCE/National Bureau of Standards

Interpolation Formulas for of Six Gases: Air, , , , , and

103

. LI5753 Mo.llSC 1984 NATIONAL BUREAU OF STANDARDS

The National Bureau of Standards' was established by an act ot Congress on March 3, 1901.

The Bureau's overall goal is to strengthen and advance the Nation's science and technology and facilitate their effective application for public benefit. To this end, the Bureau conducts research and provides: (1) a basis for the Nation's physical measurement system, (2) scientific and technological services for industry and government, (3) a technical basis tor equity in trade, and (4) technical services to promote public safety. The Bureau's technical work is per- formed by the National Measurement Laboratory, the National Engineering Laboratory, and the Institute for Computer Sciences and Technology.

THE NATIONAL MEASUREMENT LABORATORY provides the national system of physical and chemical and materials measurement; coordinates the system with measurement systems of other nations and furnishes essentia! services leading to accurate and uniform physical and chemical measurement throughout the Nation's scientific community, industry, and commerce; conducts materials research leading to improved methods of measurement, standards, and data on the properties of materials needed by industry, commerce, educational institutions, and Government; provides advisory and research services to other Government agencies; develops, produces, and distributes Standard Reference Materials; and provides calibration services. The Laboratory consists of the following centers:

Absolute Physical Quantities^ — Radiation Research — Chemical — Analytical — Materials Science

THE NATIONAL ENGINEERING LABORATORY provides technology and technical ser- vices to the public and private sectors to address national needs and to solve national problems; conducts research in engineering and applied science in support of these efforts; builds and maintains competence in the necessary disciplines required to carry out this research and technical service; develops engineering data and measurement capabilities; provides engineering measurement traceability services; develops test methods and proposes engineering standards and code changes; develops and proposes new engineering practices;

and develops and improves mechanisms to transfer results of its research to the ultimate user. The Laboratory consists of the following centers:

Applied Mathematics — Electronics and Electrical Engineering-^ — Manufacturing Engineering — Building Technology — Fire Research — Chemical Engineering^

THE INSTITUTE FOR COMPUTER SCIENCES AND TECHNOLOGY conducts research and provides scientific and technical services to aid Federal agencies in the selection, acquisiuon, application, and use of computer technology to improve effecfiveness and economy in Government operations in accordance with Public Law 89-306 (40 U.S.C. 759), relevant Executive Orders, and other directives; carries out this mission by managing the Federal Information Processing Standards Program, developing Federal ADP standards guidelines, and managing Federal participation in ADP voluntary standardization activities; provides scientific and technological advisory services and assistance to Federal agencies; and provides the technical foundation for computer-related policies of the Federal Government. The Institute consists of the following centers:

Programming Science and Technology — Computer Systems Engineering.

'Headquarters and Laboratories at Gaithersburg, MD, unless otherwise noted; mailing address Washington, DC 20234. ^Some divisions within the center are located at Boulder, CO 80303. KATIONAL BUREAU OF STANDARDS UBHAPY

Interpolation Formulas for ^^/

Viscosity of Six Gases: lOo ue)-?53 Air, Nitrogen, Carbon Dioxide, ^>o i^^ Helium, Argon, and Oxygen "i^f

Frank E. Jones

Center for Chemical Engineering National Engineering Laboratory National Bureau of Standards Washington, D C 20234

Q ft- \ ^^ftEAU O^

U.S. DEPARTMENT OF COMMERCE, Malcolm Baldrlge, Secretary

NATIONAL BUREAU OF STANDARDS, Ernest Ambler, Director

Issued February 1984 National Bureau of Standards Technical Note 1 1 86

Natl. Bur. Stand. (U.S.), Tech. Note 11 86, 26 pages (Feb. 1984) CODEN: NBTNAE

U.S. GOVERNMENT PRINTING OFRCE WASHINGTON: 1984

For sale by the Superintendent of Documents, U.S. Government Printing Office, Washington, DC 20402 INTERPOLATION FORMULAS FOR VISCOSITY OF SIX GASES: AIR, NITROGEN, CARBON DIOXIDE, HELIUM, ARGON, AND OXYGEN

Frank E. Jones National Bureau of Standards Washington, DC 20234

Equations for the calculation of viscosity for dry air, nitrogen, carbon dioxide, helium, argon, and oxygen have been developed as interpolation formulas fitted to experimental data. The approximate ranges of strict application of the equations are the ranges of (20°C £ t _< SO^C) and pressure (0.04 £ p £ 4 MPa; 0.4 <_ p < 40 atm) for the experimental data. The estimates of relative resTdual standard deviation for the fits (0.05% for air, 0.03% for nitrogen, 0.02% for carbon dioxide, 0.02% for helium, 0.03% for argon, and 0.03% for oxygen) are in close agreement with estimates of precision for the experimental data.

Key words: Air; argon; calculation; carbon dioxide; helium; nitrogen; oxygen; viscosity.

1. INTRODUCTION

In the calibration and use of flow metering devices the precise, accurate value of the viscosity of the gas being metered is often required. The availability of simple equations relating viscosity, y, to temperature, t, and , p, would be a convenience to the engineer in calculating y. In the present work, simple interpolation formulas have been developed which enable the engineer to conveniently make precise, accurate calculations of y for dry air, nitrogen, carbon dioxide, helium, argon, and oxygen using small readily available hand-held calculators. The formulas are fitted to experimental data.

The data were selected from the literature on the basis of claimed accuracy and precision, and of internal consistency. The last of these criteria is particularly important for the development of empirical equations for flow with a dominant term of the Poiseuille form, the subject of a later paper. The sets of data published by Kestin and his collaborators [1-10] meet these criteria and have been used to develop the formulas.

2. EXPERIMENTAL DATA

The ranges of temperature and pressure of Interest in the present

work are 0°C <_ t _< 50°C and < p <_ 4 MPa (0 < p <_ 40 atm). The experimental data used in developing the equations for u are listed in

Tables 1-6. In Table 1 for dry air, the data in the first and third groups are taken from Table III of [1]; the second and fourth groups are taken from Table II of [2]; and the 50.00°C point is inferred from the data in [2], In Table 2 for nitrogen, the first and third groups are taken from Table XI of [1]; the second, fifth, and sixth groups are taken from Table VIII of [2]; the fourth group is taken from Table III of [3]; the 20.58°C point is taken from Table XI of [1]; and the

24.920°C point is taken from Table III of [3]. In Table 3 for carbon dioxide, the 20.00°C data are taken from Table V of [1]; the second group is taken from Table IV of [2]; and the 25.00°C point is taken from Table III of [4]. In Table 4 for helium, the first group is taken from Table VII of [1]; the 25.00°C point is taken from Table IV of [5]; and the third group is taken from Table V of [2]. In Table 5 for argon, the first and second groups are taken from Table IV of [1]; the third group is taken from Table IV of [3]; and the fourth group is taken from Table III of [2]. In Table 6 for oxygen, the first and second groups are adjusted values from Table XII of [1]; the third group is taken from Table V of [10]. 3. DEVELOPrCNT OF EQUATIONS

The equations for y (t,p) were developed using a procedure [6] in which u(t,p) is expressed as

y(t,p) = Pq (t,0) + Ay(p). (1)

where y (t,0) is the viscosity at "zero pressure" and Ay(p) is the

experimental value of viscosity, Vr^^p.^ > minus y (t,0). The procedure will be outlined below.

Air

The data pairs, (20.00°C, 181.94 yg/cm^s), (23.44°C, 183.75 yg/cm-s),

(23.90°C, 184.21 yg/cm-s), (25.00°C, 184.62 yg/cm-s), and (50.00°C, 197.31 yg/cm-s), taken from Table 1, were fitted by least squares to an equation quadratic in t to enable calculation of y at 0.101325 MPa, y-,. The small departures of p from this value of pressure for these pairs are not significant. The resulting equation is

= 10"^ y-j 170.368 + 0.605434 t - 1.33200 x t^, (2)

where y-, is in yg/cm-s and t is in °C. To reduce y-, to "zero pressure,"

y , the value of the increase in y per 0.10325 MPa, 0.11 yg/cm-s, estimated from the tabulated data was subtracted from y-. resulting in

y^ = 170.258 + 0.605434 t - 1.33200 x lO'"^ t^. (3)

Calculated values of y are listed in Table 1 . . . , . .

The difference, Ay, between the experimental values of y and calculated values of y is also listed in Table 1. The values of Ay were fitted by least squares to an equation quadratic in p. The resulting equation is

10"^ Ay = -2.44358 x + 1.17237 p + 0.125541 p^. (4)

Equations (3) and (4) were added together to synthesize the final equation for calculating y:

-3 .2 = + Ay = 170.256 + 0.605434 t - 1.33200 x 10 " t y ca 1^ y^

+ 1.17237 p + 0.125639 /, (5)

where y ,, is in ug/cin«s, t is in °C, and p is in MPa. Values of y^^i^ and caic. — Co I c differences between y^^^^ and experimental values, y^^^. , are listed ca I c me as . in Table 1

The estimate of residual standard deviation (RSD), that is, the

0-^^ estimate of the standard deviation of y^^i^ -y^^oc » ''s yg/cm-s for ca I c . me a 5 the 19 differences (n=19). The estimate of the relative residual standard

deviation (RRSD), that is, the ratio of RSD to the mean , is 0.05^. y^^^^me as

Kestin and Leidenfrost [1] estimated that for their measurements "a precision ranging from +0.01% to +^0.07%, depending on the gas, has been achieved. The final accuracy of the measurements is estimated to be of the order of +0.05%." DiPippo and Kestin [2] estimated the precision of their data to be +0.05%; they reached the conclusion that "no meaningful assessment of the accuracy of the present data can be given if by accuracy we mean the irreducible discrepancy between the best measurements available at any particular time." For y in yPa«s , t in ^, and p in MPa , equation (5) becomes

y.--,_ = 17.0256 + 6.05434 x 10 " t -1.33200 x 10 " t calc.

10'^ + 0.117237 p + 1.25639 x p^. (6)

10" For y in lb ft" sec" , t in ^, and p in PS I , equation (5) becomes

p = 11.4407 + 4.06833 x 10"^ t -8.95063 x 10"^ t^ ca I c.

+ 5.43165 X 10"^ p + 4.01338 x lO"'^ p^. (7)

Nitrogen

The data pairs, (20.58°C, 175.93 yg/cm-s), (21.98°C, 176.56 yg/cm-s),

(23.15°C, 177.05 yg/cm-s), (24.920°C, 177.94 yg/cm-s), (29.75°C, 180.08 yg/cm-s), and (49.37°C, 189.71 yg/cm-s), taken from Table 2, were fitted by least squares to an equation quadratic in t to enable calculation of y-,. The resulting equation is

y, = 167.335 + 0.392728 t + 1.22474 x 10"^ t^. (8)

To reduce y-, to y , 0.11 yg/cm-s estimated from the tabulated data was

subtracted from y, resulting in

10"^ y^ = 167.225 + 0.392728 t + 1.22474 x t^. (9)

The values of Ay were fitted by least squares to an equation quadratic

in p. The resulting equation is .

10'^ 10'^ Ay = -1.12860 X + 1.24165 p + 9.87206 x /. (10)

Equations (9) and (10) were added together to synthesize the final

equation for calculating y:

10"^ Vr..^. = 167.214 + 0.392728 t + 1.22474 x t^ ca I c

+ 1.24165 p + 9.87206 x 10'^ p^, (11)

where y-,-]- is in yg/cni'S , t is in ^, and p is in MPa,

The RSD is 0.05 yg/cm-s for n = 50; the RRSD is 0.03%. Kestin et. al.

[3] estimated the accuracy of their experimental measurements to be of the

order of +0.2% and the relative precision to be 0.03%. The accuracy and

precision estimates quoted above for air from [1] and [2] apply also to

nitrogen.

For y in yPa'S , t in ^, and p in MPa , equation (11) becomes

-2 . , . oo...M M It I I .. .n-4 .2*. y„-,. =^ 16.7214II-* /xi/i +-^ 3.92728< \Jk y I y X xV 10III " t + 1.22474 x 10 " t calc.

10""^ p^' + 0.124165 p + 9.87206 x (12)

10" For y in lb ft" sec" , t in °C, and p in PSI_, equation (11) becomes

-2 r, ^or.r.-7JfllJ / »* Tn-5I f 1 x2-f- y.,n. =^ 11.2363II • <.v\ < +-*- 2.63901y tr^

+ 5.75263 X 10"^ p + 3.15351 x 10"^ p^. (13) . .

Carbon Dioxide

The data pairs, (20.00°C, 146.63 yg/cm-s), (25.00°C, 149.09 yg/cm-s),

(30.40°C, 151.81 yg/cm-s), and (49.12°C, 161.75 yg/cm-s), taken from Table 3, were fitted by least squares to an equation quadratic in t to enable calculation of y-,. The resulting equation is

y, = 137.335 + 0.441133 t + 1.12987 x 10"^ t^. (14)

To reduce y-, to y , 0.13 yg/cm«s estimated from the tabulated data was subtracted from y-, resulting in

u = 137.205 + 0.441133 t + 1.12987 x 10"^ t^. (15)

The values of Ay were fitted by least squares to an equation quadratic in p. The resulting equation is

10"^ Ay = 0.133827 + 6.28105 x p + 0.562974 p^. (16)

Equations (15) and (16) were added together to synthesize the final equation for calculating y:

10"^ v^r.i. = 137.339 + 0.441133 t + 1.12987 x t^ Cd I C

10"^ + 6.28105 X p + 0.562974 p^, (17)

where y-^-i- is in yg/cm'S , t is in %^, and p is in MPa.

The RSD is 0.04 yg/cm-sfor n = 11 ; the RRSD is 0.02%. Kestin et. al

[4] estimated the accuracy of the 25.00°C point, "expressed as the maximum .

estimated uncertainty in the quoted values," to be +0.1%. The accuracy and precision estimates quoted above for air from [1] and [2] apply also to carbon dioxide.

For y in yPa'S , t in ^, and p in MPa_, equation (17) becomes

-2 -..-4 . . . ,onn. .2 y_,-,_ = 13.7339 + 4.41133 x 10 " t + 1.12987 x 10 " t calc.

+ 6.28105 X 10"^ p + 5.62974 x 10"^ p^. (18)

10" For y in lb ft" sec" , t in °C, and p in PS I , equation (17) becomes

10"^ 10"^ Uo.T. = 9.22876 + 2.96428 x t + 7.59238 x t^ ca I c

+ 2.91005 X 10"^ p + 1.79836 x 10"^ p^. (19)

Helium

The data pairs, (20.00°C, 196.14 yg/cm-s), (24.28°C, 198.19 yg/cm-s),

(25.00°C, 198.59 yg/cm-s), (30.56°C, 201.18 yg/cm-s), and (52.79°C, 211.04 yg/cm-s), were fitted by least squares to an equation quadratic in t to enable calculation of y-,. The effect of pressure on y for helium is so

small that no correction is made to reduce y-, to y . The resulting equation is

= = - 10"^ Mt Vr. 185.946 + 0.530773 t 1.04982 x t^. (20)

The values of Ay were fitted by least squares to an equation linear in p. The resulting equation is

Ay = 2.91664 X 10"^ - 6.99813 x 10"^ p. (21) Equations (20) and (21) were added together to synthesize the final equation for calculating y:

= 185.975 + 0.530773 t - 1.04982 x 10"^ t^ y^calc.T

10"^ 6.99813 X p, (22)

where y , is in yg/crri'S, t is in °C, and p is in MPa. calc. ^^^^ —

The RSD is 0.04 yg/cm-s for n = 15; the RRSD is 0.02%. The accuracy and precision estimates quoted above for air from [1] and [2] apply also to helium. The 25.00°C °C value is estimated to be accurate within + 0.1% [5],

For y in yPa^s , t in ^, and p in MPa , equation (22) becomes

-2 . . r.nr.nnfi/iiii.# .. -,^-4 .2.*. y„,-,_ =^ 18.5975IX *-*\-l/»-» +H- 5.30773t-.

- 6.99813 X 10"^ p. (23)

10" For y in lb ft" sec" , t in °C, and p in PS I , equation (22) becomes

10"^ 10"^ y^.i. = 12.4969 + 3.56663 x t - 7.05446 x t^ ca I c.

3.24228 X 10"^ p. (24)

Argon

The data pairs, (20.00°C, 222.86 yg/cm-s), (24.39°C, 225.79 yg/cm-s),

(25.00°C, 226.36 yg/cm.s), and (50.31°C, 243.43 yg/cm-s), taken from Table 5 were fitted by least squares to an equation quadratic in t to enable

calculation of y-. . The resulting equation is .

-^ 10"^ y-, = 208.940 0.702190 t - 3.30712 x t^. (25)

To reduce y-, to y , 0.21 yg/cm-s estimated from the tabulated data was subtracted from y-, resulting in

^ y^ = 208.730 + 0.702190 t - 3.30712 x 10 t^. (26)

The values of Ay were fitted by least squares to an equation quadratic in p. The resulting equation is

10"^ Ay = 3.17420 X + 1.73987 p + 0.152358 p^. (27)

Equations (26) and (27) were added together to synthesize the final equation for calculating y:

10"^ y.=i. = 208.762 + 0.702190 t - 3.30712 x \} ca 1 c

+ 1.73987 p + 0.152358 p^, (28)

where u^^-i- is in yg/cm^s , t is in ^, and p is in MPa .

The RSD is 0.07 yg/cm-s for n = 38; the RRSD is 0.03%. The estimates of accuracy and precision of the experimental data for argon are those quoted above from [1], [2], and [3].

For y in yPajs, t in °C, and p in MPa, equation (28) becomes

-2 . . .n^no in-5 .2 y__-,- = 20.8762 + 7.02190 x 10 ^ t - 3.30712 x 10 " t calc.

+ 0.173987 p + 1.52358 x 10'^ p^. (29)

10 ,

10"^ 10"^ y , = 14.0282 + 4.71850 x t - 2.22228 x t^ calc.

+ 8.06091 X 10"^ p + 4.86689 x lO"'' p^. (30)

Oxygen

The data pairs, (20.00°C, 203.17 ug/cm-s), (25.00°C, 206.25 yg/cm-s), (35.00°C, 212.18 yg/cm-s), and (50. 00^*0, 220.80 yg/cm-s) taken from Table 6, were fitted by least squares to an equation quadratic in t to enable calculation of yi- The resulting equation is

10"^ y^ = 190.539 + 0.650043 t - 8.97542 x t^. (31)

To reduce y] to \Sq, 0.12 yg/cm-s estimated from the tabulated data

was subtracted from y-] resulting in

y^ = 190.419 + 0.650043 t - 8.97542 x 10"^ t^. (32)

The values of Ay were fitted by least squares to an equation quadratic in p. The resulting equation is

Ay = -0.0239595 + 0.130897 p + 1.323418 x 10"^ p^. (33)

Equations (32) and (33) were added together to synthesize the final equation for calculating y:

= 190.395 + 0.650043 t - 8.97542 x 10'^ t^ ^..1.calc.

10""^ + 0.130897 p + 1.32418 x p^ (34)

11 where locale. ^^ ^" Mg/cm-s , t is in ^, and p is in atm.

The RSD is 0.06 g/cm s for n = 16; the RRSD is 0.03%. The estimates of accuracy and precision of the experimental data for oxygen are those quoted above from [1], [2], and [3].

For y in uPa-s , t in ^, and p in MPa , equation (34) becomes

" '•9.0395 + 6.50043 x 10"^ t - 8.97542 x 10"^ ^calc t^

10"^ + 0.129185 p + 1.28975 x p^. (35)

— fi -1 1 10" For y in lb ft" sec" , t in **C, and p in PSI, equation

(34) becomes

,^ = 12.7940 + 4.36809 x 10"^ t - 6.03121 x 10"^ t^ calc.^

10"^ 10"^ + 5.98525 X p + 4.11587 x p^. (36)

4. RANGES OF APPLICATION OF THE EQUATIONS

The equations developed in this paper are interpolation formulas fitted to experimental data. The range of strict application is indicated by the range of t and p of the experimental data listed in

Tables 1-5. In the absence of measurements of y for below

20°C of comparable quality to that of the data in the tables, there are two options in extending calculations below 20" C (0 < t < 20*0): 1) apply the extended law of corresponding states [3,7-10]; 2) extrapolate using the equations developed here, with probable loss in accuracy.

Either option might be followed until suitable experimental data at temperatures below 20''C became available.

12 .

Hanley et al . [11] have developed a functional form to represent critically evaluated viscosity and thermal conductivity coefficient data, and have generated tables. The gases treated by Hanley et al include nitrogen and argon. Values of y calculated using the formulas in the present work have been compared with interpolated values from the tables in [11] at 0°C, S^C, 10"*C, IS^'C, 20''C, and 25°C, and 0.1

MPa, for nitrogen and argon. The deviation of the values in the present work from the values from reference [11], expressed in percent are: for nitrogen, +0.26% at O^C, +0.02% at S^C, -0.12% at 10°C,

-0.28% at IS^'C, -0.32% at 20*0, and -0.38% at 25''C; for argon, -1.16% at 0°C, -1.10% at S^C, -0.99% at 10*^0, -0.86% at 15°C, -0.78% at 20''C, and -0.66% at 25°C. These deviations are all well within the uncertainty, j^ 2%, estimated by Hanley et al . for their tables.

5. CONCLUSIONS

Equations (interpolation formulas fitted to experimental data) for the calculation of y for dry air, nitrogen, carbon dioxide, helium, argon, and oxygen have been developed. The estimates of relative residual standard deviation for the fits are in close agreement with the estimates of precision for the experimental data in the above stated ranges.

6. ACKNOWLEDGMENTS

The author is deeply grateful to J. Kestin for his helpful suggestions and to him and his collaborators for the experimental data on which the present work is based. The typing of the manuscript by

Susan Johnson is gratefully acknowledged.

13 .

Nomenclature p = pressure

RRSD = estimate of relative residual standard deviation

RSD = estimate of residual standard deviation t = temperature in °C

Ay = u -u meas. o y = viscosity

= y ^i calculated value of viscosity ca I c

= y_ experimental value of viscosity-^ meas. ^ y = viscosity at "zero pressure" y^ = viscosity at 0.101325 MPa

14 References

1. Kestin, J., and Leidenfrost, W., "An Absolute Determination of the Viscosity of Eleven Gases," Physica, Vol. 25, 1959, pp. 1033-1062.

2. DiPippo, R., and Kestin, J., "The Viscosity of Seven Gases up to 500°C and Its Statistical Interpretation," Proceedings Fourth

Symposium on Thermophysical Properties , American Society of Mechanical Engineers, College Park, Md., 1958, pp. 304-313.

3. Kestin, J., Paykoc, E., and Sengers, J.V., "On the Expansion for Viscosity in Gases," Physica, Vol. 54, 1971, pp. 1-19.

4. Kestin, J., Ro, S.T. , and Wakeham, W.A., "Viscosity of Carbon Dioxide in the Temperature Range 25-700°C," J. Chem. Phys., Vol. 58, No. 8, 1972, pp. 4114-4118.

5. Kestin, J., Khalifa, H.E., Ro, S.T., and Wakeham, W. A., "The Viscosity and Coefficients of Eighteen Binary Gaseous Systems," Physica, Vol. 88A, 1977, pp. 242-260.

6. Kestin, J., and Whitelaw, J.H., "The Viscosity of Dry and Humid Air," Int. J. Heat Mass Transfer, Vol. 7, 1964, pp. 1245-1255.

7. Kestin, J., Ro, S.T., and Wakeham, W.A., "An Extended Law of Corres- ponding States for the Equilibrium and Transport Properties of Noble Gases," Physica, Vol. 58, 1972, pp. 165-211.

8. Kestin, J., Ro, S.T., and Wakeham, W.A., "Viscosity of the Noble Gases in the Temperature Range 25-700°C," J. Chem. Phys., Vol. 56, No. 8, 1970, pp. 4119-4124.

9. Kestin, J., and Mason, E.A. , "Transport Properties in Gases (Comparison Between Theory and Experiment)," AIP Conf. Proc. No. 11, 1973, Edited by J. Kestin, pp. 137-192.

10. Hellemans, J.M., Kestin, J., and Ro, S.T., "The Viscosity of Oxygen and Some of Its With Other Gases," Physica 65_, 1973, pp. 362-375,

11. Hanley, H.J.M., McCarty, R.D., and Haynes, W.M., "The Viscosity and Thermal Conductivity Coefficients for Dense Gaseous and Argon,

Krypton, , Nitrogen, and Oxygen," J. Phys. Chem. Ref. Data 3^, 1974, pp. 979-1018.

15 ^ —

w fO <— OJ CO E • i—t— I— O'^i— Ovj"^ o CM r— csj r~^ CO UD "^ ^ 3. b: OOOOOOOi— r— O r- I— o o O r— 1— O 1 O oooooooo o o o o o o o o o o o c_> cn + 11 I I I I + + I I I 1 + + + + + n— ;i

fO '— o 3-

LD r-. OO OJ '^ CTi (^ un r-. CM cri LO CO CM CM CTi CO CO CO CO CO CT. CO "=d- CM o o o .— CO CO o

03 CM CM CO CO "^ '=:3- CO -^ ';:r ^ r^ o C71 CO CO OD CO CO CO CO CO CO CO CO CO CO CO CO CO CTi cn CTl 7X

w

E '— LDi— cT^Lnr^Lnr^ COVJDr— ^r^OO IsOCMOJ CTiCO ;1CJ r-LnO>^0(~Dr— CO OOCOCMi— O 0"^iX) o-^ OO cm ooi— 1— cMCMcoLn o o o o o o o CM Ln o o ;3.

QS-

OJ t/1 o COCOCOCOCOCOCOCO CMCnCMI^UDCD 1JD <^Di OE COCOCOCOCOCOCOCO r^roocTicno LOLOun — cn 3. O . .

+-> 1 1 1 [ 1 -^ — — — I— I— — — r- CO'^'=:d-COCO'^ ^^d'"^ r^vo cn cococococx)cococo cocooocococo cococo cncn oI/) o

CO "^ OO 00 O CO o CO 1— CO CO CM CO CO LD cn CO u cn CO CO OO CO 1^ '=^ OO CM r- O (n CM cn CO OJ CM CM CO CO o C7) 5 CO CO CO CO co CO CO OD CO OO CO CO OO CO CO CTi cn cn cn r— r—

O

oi

ocncncMi— oTor^ 1^ cn r— ^ CO CM O I— CM O CO o c ocor^cMcdO':d-cn o uD r^ o CO "^ o o .— r^ o o I/) to =d- r-^ .— -Jd- CO r- <^ O c— I— o o r— CO un cu 1— CM OO .— CO s_ CI.

n3

o oooooooo ^ CO o o cn r^ o cn o^ I— cn o oooooooo «^ >,D o cri CO cn o CTi cn cn CO o -a Ci.C_) cu oooooooo OO <=j- -^j- CO CO OO LO >^ "^d- cn cn o to OviCMCMOJCMCMCMCM CM CM CM CM OJ CM CM CM CM LD «3

16 —

00 fO ^-^ OJ CO e • r>.ooc\jooooor^oro'=d- cvji— roLOi— "^1— ooi— CO LD ^ ^ I— ^ 3. 1= OOOOi— I— oooo 000000(3000 o o o o o 1 o oooooooooo oooooooooo o o o o o u O) I I + + + + + + + + + I I I + I + I + + + I + I r— 3. > ^ fOo 3.

r^uDUDi— CTiOoocytLni— LD O CO 00 ^ UD LD -^ C\J E 1— LOi— VOOUDi— IX><>DV£> LO Ln 00 CT. a\ o 1— o 00 >— CO CO CTi o 1— cocor^r>.(£><>DLnLDLO r^ CO <=d- (NJ cri CO r^ o CD CO 00 r^ r^ r^

CO £ r^r^r-^i— o^oocy>csaLno r-i— Lncoooojtxjr— ^00 CO I— CO CM LO 3.0 u30Lnoroo^"vi-i— oo I— I— C\J

OJ CO =d-"=d-coor^ o o LO LO LO un un LD un LO LO LO ro«^LO(X)i>Dr--oocooa-i 00 00 OO 5 00 fO >) I— 4J CD LO LO LO LO LO LO LO LO LO LO i^ r-. t^ r^ •r— :3. r^ r^ r^ 1-^ r-^ r^ r-> 1^ r-. r-^ r-. r^ oCO u to

^''d-'^COtOOtXJO^CVI c\j <* LO OJ cn CO LD CVJlX)i— LOCnLOOUD>£> LO LO LO CO 00 00 00 O (T> 1— o r— CO CO CTl

I— cocor^<^<£)LOLO LO ^o vo r». r-« CO r-~ OJ cy> CO r>>. oor-^r^r-^r-^r^r^r^r^ I-, r^ t-^ r^ CO CO r-^ r^ r^ 3-

I^ LO ojcO"^cy>r^cT>coococo r^Lor^nrooj'^LOO'^l- LO I— o o cr> oo«*OLOi— r>.roOLoro LOOLDLO'^COOJr— I— O 00 OJ LD CO CTl to +-> CO fO 'd-O00<::fi— r^-vfr— OO r- CO "s^f r^ O O) cn OJ 1— I— 1— «* CM I— s_

OOOOOOOOOO LOOO-^J-LDCM-^r^LOCMLO o o o o o r~~cn^Lor-~oooi— cvji o o o o o Q.O oooooooooo E^- oooooooooo 1— 1— '^"vd-'^^LOLOOOrO LO LO LO LO LO cu CMOJOJCSJOJOJOJCMCSJCVJ OJCSJOJCMOJOJOJOJOJOJ OJ CM CVJ CM OJ

17 ' — —

C/1 ooooooocoooooooooooo om cCO "=^o oO^ I—o COCOOC'OOOOOOOOOOOOCO o o o o o o iir I1+IIII1II+ + +I1 + + + + I + 3-

ro c\j cTi a~v CO r^ LD ro CO 1— "^ .— o cvj CO OJ O £ cr. o o 1— c\j oo ^ tn ^£5 1— Cvj o cTi r--. 1— u 03 ^-^ i^ CO CO CO CO CO CO CO CO cocoo-icria^criOi— c\j CO ^ o o o en CT. <_) CD r^ r^ r^ 1^ r-- r-- r^ r--. r^ r^r^r-^r^r^-r^cococo 00 CO CO CO cr» CO CO ^ 3.

r^LDcvicOi— orocoooocricovocor^Ln":d-co^CM V£) CO 1 I— Ovj(r)coLnuDor-^C7~iCorvjcor^oco<.DuD^£><.D o o t— CVJ !—

o OJ 00 o E r— r— r— r— r— I— r— r— r-r-i— r— i—r— I—.— r— r— I—I— CTiCO OJCOO o u cocococococococococococooDcocooococococo cno cor-^vo 3- OJ en r^r^r^r^r^r^r^r^r^r^r^r^r-~.r-->.r^i^r^r^r~^r~~. CTiO criC3^CT> 3. r^r^i^r^t~v.r^r^r^r^r^r^r^r^i-^r^r^r^r^r>,r-^ r^oo ODCOOO

CO VO CO oicvji— «:j-cnooi— o-j^r-- CTi CO vo LO CTi un CO CO OJ CO .— CO «:^

.— 1 Oi o — cO'^d-'STLnr-^cocnOi Ln CO ix) "^d- «^ >* -^ O r- o^ o r^ C5^ (T> -~^ OJ t~-- CO CO cococococococococricri a~i CTi o r— OJ CO ^ o O CTi o E CD r^ 1^ r^ r^ r-- CO CO CO CO 00 CO CO CO o-i CO 3. 3.

o

o=^oc3^Lnc— a~)CTicor-^cri"*i— r^o^— «^covocn CD Oi— ooi— =:d-r->.r~-i— oojLncoixicOLn'^coi— <— ix)ix> "^-vfiX) oo or^cocDr^cocriLn*j3ooocricriuncoo'X)cococr> oo uD^^o oo •p- +-> ojcocO':d-'=d-Ln'X)r^cococnc\j<^0'^cri'^CJ^ I— I— OJ OJ OJ CO CO 3 u «ou os- o H- (OOOOOOOOOOOOOOOOOOOO un CO o CO r>. CO OJ OOOOOOOOOOOOOOOOOOOO CO r-- CO -o Q-CJ CD E— LnLnLoi^u^u^i^LnLrsLnLnLnLnLnununLnLnLnLn cTi en CT) o~i en o "^ to OJ CMOJOJOvJOJOJOJOJOJOJOvJOJOJOJOvlOJOvJOsJOJOJ OJ OJ <* "^ <* OJ OJ ^) to

18 —

l/l JC-^ a> to- C2 « n— n- CNI O CO o (r> ^ o o o 3- E O O O O o o o o o o o . • • • • • • • • • • 1 o • ^^^ o o o o o o o o o o o U CD 1 + + 1 + 1 1— 3- Oro 3-

r- «^ r- .— LD CO en r— LO IX) CT. "^ O r-N. CO O r— O ea ^^ 00 r^ r^ r^ "^ vi-

T3 CO «:;f Lf) r— OO O Ln CO o tn LO LO •r- 3. • I— "5a- en tn ro >— o X < E o O CMi—OOOOO o o o o CJ) p. O CO ofO J2 n3 M- O CO CO CO CO CO CO CO CO IX) o O •^ "^d- >^ «^ "vl- ">!- "^ ^O 1X3 1X5 IX) 1X5 IX) IX) 1X3 VO I— I— I— 00 "=1- "^ '^ -vl- "vt- Lf) "* COO >=d- =:J- LD *X) O (/I

VJ c\j ro en 1— CO CO IX) CO I— Lf) en t/) IX) cr> CO o r^ IX) LD oo CO r^ o 03 CO r^ r^ r^ IX) uD 1X3 I— r— r— Cn O) ^ -53- «* «=f "Sd- «* '!^ IX) LO U3 •* ;3.

IX> O **^"'~^ csj oo CO CO 00 r^ "^ 0~i 1X3 CO O • E o vo I— r-. oo cr> IX) IX) O O o CO 4-> CO fO CO -^ .— r-^ ^ o o 1— r— 1— ,— a)-v_^ r^ f"^ r'^ s- D.

o o o o o o o >=d- O CVJ o o o o o o o o CO "^ 1 o Q.C_3 E-— o o o o o o o o o cn LO C\J C\J C\J C\J CM CVJ C\J CO CO * CM

19 — — — —

C/1 rtJ. ooooooooc\jroc\JCsji>DU3Lnc\j C\Jo o'^ oCO COo I—o CVIo oooooooo o o o o o o o I + + I + I I + I I + I + I + CD 3- o 3-

CO o CO r-^ 1— Ln uD CTlOOOOr— I I— 1— (XI I— CM I— o I— o CO CO 00 — 1— I— r—

1 U CD criCT>cricriCT>cTicTi CTi C7> CT> O O — I 72- 3. r— 1— 1— CM (NJ CM CM

OLOCMLDi— councD "vi" CM Ln 1— CO O 3- E CMi— r— Or— OOO O o o o o o o

OJ

O)

J2 >> . s eo -l-> to 1— •1 • ^'*'*"v]-«:d-':t-<=d-'^ CO I— un CT) CO ^ to O E O CM r— I— r— O O 3- O a ~~^ O O r— I •1— 1 3- I— — I— CM CM CM CM

^crioocricor--cT>'^ CM CTi O CO to 'd- E cncriooOi— 1— I — I— CM 1— r- O o LnLOtotoiotototo CO 00 CO r— r— I— I— C7) cncT>cy>c3~icricricricri CTi C3^ Cr> O O I— 1 3. I— r— r— CM CM CM CM

LDOOi— OCrii— r^ vD 1— r-^ 1— r-^ 00"^OU3i— COCOO Ln o r-~ o to O to to "^ :— CO "* I— r~- >*

1 OJ CO CM I— — . s- a.

oooooooo o 1— CO cTi to o en oooooooo o o CM >=d- Lo o r^ Q.t_> oooooooo 5*- <^ O O OO CM CMCMCMCMCMCMCMCM CM CM CM CO CO un LO

20 —%^ — —

CO fO s CO ci; I— E • r-^i—'vt-i— r— s;j-roi— I— C\J 00 Vi3 .— CT) Lnco>^rooocM<;i-r— coi-n>^0"^c\jroc\ji 3- b oooooooooo O I— o o o ooooooooooooooooo ' o oooooooooo o o o o o ooooooooooooooooo o en I + I + I I I I 1 + + + + + I 1+1+111111+ I+II + ^— 3. *.— uto ;i

air^r^crvLn^p^LOO":^- CO O CM 00 CT^ oO(r)LOLnroLO^'vr>5:fr^cooooO'^c\j LOi— roLocOi— '^oocor-^ CVJ (13 O LD OJ cj-vi-i^coocM^r^OLni—coo^cooocjr^

to crir^vl-«*roc\ic\jCM CM o cn r~-. <>D uD<>D<>Dv£ir-^r-^r^r~^oooocnooi— cM«d-un o C^JC^JCMC^JC^JCMC^aCM<^JC^J ro OO CM CM CM CMCMCMCMCMOJCMCMCMCMCMCOrOrOCOrOrO 3- CMCM

CMCMr^^cM*i— LnLnuDoouncoco OLnr^cncMLOcocMr— o r— CO 00 "* CO CMCouDr-^Oi— "^v£)Oi-nocMCOc\jr^i— =3- t^«*COCMCMi— oooo (X) OO I— O OOOOi— 1— I— r— CMCMCO«^<5i-Lnv£)COCri CJ» + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + C 3- o s„^^^ CD s- » *-— fO M 01 1— •1— • 00 CO CO 00 CO COOOCOCOCOOOCOCOOOCOCOOOCOCOCOCOOO CO o E O O O O o OOOOOOOOOOOOOOOOO O 3. o O "- CMCMCMCMCMCMCMCMCMCM 1^ <£i KD <^ KO <^(^VO<^V£)V£)VOVOVOV£>V£)VOUD<^ to CD CMCMOJCMCMCMCMCMCMCM CM CM OJ CM CM CMCMCMCMCMCMCMCMCMCMCMCMCMCMCMCMCM •1— :n. OJCMCMCMCMCMCVJCMCMCM CM CM CvJ CM CM CMCMCMCMCMOsJCMCMCMCMCMCMCMCMCMCMCM

(nUDi— CX)<-DCOOVOi— CO I— CM U3 r-~ 00 ^LnocMCOLncriLncMcncoco-^vDcoi^r— V£)i— "vj-LnOOr— LOCOCOr^ CM "^ CTi en CO co>*r^00i— CMuDvDr^r^i~^r-^oococr>ooi— cm'^lo CD CMCMCMCMCMCMCMCMCMCM CO CO CM CM CM CMCMCMCMCMCMCMCMCMCMCMCOOOCOCOCOCO :3- CMCMCMCMCMCMCMCMCMCM CM CM CM CM CM CMCMCMCMCMCMCMCMCMCMCMCMCMCMCMCMCM

CTl >^ Or- ^oooLnoov£)cncooo<.Dr^cMLnvor— • £ c— CMUDCOCMCTir^OOCO O'vl-co>*i— coLnoocMOOO<^co"^cOLncj^ CO -M 00'=:l-0 ^ li3 CO O ooooi— Oi— unO'j^cMCMcoLncrir— 1 1/1 ra 1— O I— 00 "* I— r--. ^ — o o r^ I— <:a- r^ I— cMco>^Lnu3r^coocMLnocM'^(TiLno 1 1 1 S- CO CM — I— I CM CM r— — — I— CMCMCMCMCO=^

oooooooooo o o o o o ooooooooooooooooo00000000000C300000 0-C_) oooooooooo LO LD LD LO LO LnLOLnLOLnLnLDLOLnLnLnLnLnuoLnLOun CD CMCMCMCMCMCMCMCMCMCM CM CM CM CM CM CMCMCMCMCMCMCMCMCMCMCMCMCMCMCMCMCM

21 > — ^'^ ———

to (O.p— O) trt * s o en un ro c\j I :^ t= I— o C\J f— o o ' u o o o o o o o en + + I I I +

I ;i •^^ rOa p-

n CO oo c\j Ln "^ cri CO CTi CTi 1— "^ r— O Ln Ln CM C\J "^ "^ «:d- "=3- 3- CsJ C\J CM CM CvJ CM

• 1— oo cTi CM r-~ 1 3- E CM 1— O -^ OO CM

4-> c C o o u en S- ^—

V) oo cn CO Ln r-^ oo CO CO r--

• E VO I^ CO kO CM CO CO +-> LD o r^ -^d- CM o CO (O Q)- s- D.

CO cn r-^ CO CO I CO CO CO '::^ 00 CO

E — ^ "^ cn cn cTi o cu CM CM ^ "^ '=^ cn

22 in

(U in S* rO<*'d-fOl— OO r— 1— •SfCOCOOO co r— 3-e I— OOOOO O OOOi— o o o • • • I {^ •••••• ••••• "^ • ooQOoooo cpoooo ooo «— 3. o 3.

'JCVJCMO0CVJCVi CMCMCMCSJCM CVJCVJCVJ

trt

r— O00'?^<«0L0l— r-> VO 00 «?1- IT) «;J- o

'>oy3y3i^<.o^o<^^o 1— r— i— r— 1— f— r~co 1— I— I— .— f— f— OU3 o o 00000000 u 3- cooorororororooo y3v£5U3i«o^o CMO to 00000000 OOOOO Or— CM CMCMCMCMCMCMCVICM CMCMCSJCMCSJ CMCMCM

• • r~-i.OCT>ocMi— r^co T^ <7\ 10 00000000 r— I— O O O O r— CM XL 3- CMCMCMCMCMCMCMCM CM CM CM CM CM CM CM CM

00 OCMC0r-~CT»0000 rO CM r- O CT^ ^OLf)>— r>-oocriv£) 00 CM ^O 00 CT) (/) »-> • ••••••• • • • • • r— OO-vTi— p-~^00 I— 00 «:J- r«- o ^ CM I— a.

00000000 OOOOO 000 Q.O 00000000 Lf) Lf) un LD LO Lf) LD O E — CMCMCMCMCMCMCMCM CM CM CM CM CM CM 00 UO

23 NBS-n4A (REV. 2-8C) U.S. DEPT. OF COMM. 1. PUBLICATION OR 2. Performing Organ. Report No, 3. Publication Date REPORT NO, BIBLIOGRAPHIC DATA TT^ 1186 February 1984 SHEET (See instructions) NBS

4. TITLE AND SUBTITLE

Interpolation Formulas for Viscosity of Six Gases: Air, Nitrogen, Carbon Dioxide: Helium, Argon, and Oxygen

5. AUTHOR(S) Frank E. Jones

6. PERFORMING ORGANIZATION (If joint or other than NBS. see instructions) 7. Contract/Grant No.

NATIONAL BUREAU OF STANDARDS DEPARTMENT OF COMMERCE 8. Type of Report & Covered WASHINGTON, D.C. 20234 Final

9. SPONSORING ORGANIZATION NAME AND COMPLETE ADDRESS (Street. City. State, ZIP)

Same as item 6.

10. SUPPLEMENTARY NOTES

[3J Document describes a computer program; SF-185, FlPS Software Summary, is attached.

11. ABSTRACT (A 200-word or less factual summary of most significant information. If document includes a significant bibliography or literature survey, mention it here) Equations for the calculation of viscosity for dry air, nitrogen, carbon dioxide, helium, argon, and oxygen have been developed as interpolation formulas fitted to experimental data. The approximate ranges of strict application of the equations are tne ranges of temperature (20°C < t < 50°C) and pressure (0.04 < p < 4 MPa; of relative residual 0,4 £ p 5 40 atm) for the experimental data. The estimates standard deviation for the fits (0.05% for air, 0.03% for nitrogen, 0.02% for carbon dioxide, 0.02% for helium, 0.03% for argon, and 0.03% for oxygen) are in close agreement with estimates of precision for the experimental data.

12. KEY WORDS (S/x to twelve entries; alphabetical order; capitalize only proper names; and separate key words by semicolons) air; argon; calculation, carbon dioxide; helium; nitrogen; oxygen; viscosity

13. AVAILABILITY 14. NO. OF PRINTED PAGES \~X\ Unlimited 26 For Official Distribution. Do Not Release to NTIS I I rX\ Order From Superintendent of Documents, U.S. Government Printing Office. Washington, D.C 20402. 15. Price

Q^ Order From National Technical Information Service (NTIS), Springfield, VA. 22161

USCOMM-DC 6043-P80

NBS TECHNICAL PUBLICATIONS

PERIODICALS Building Science Series— Disseminates technical i-formation developed at the Bureau on building materials, components, JOURNAL OF RESEARCH—The Journal of Research of the systems, and whole structures. The series presents research results, develop- National Bureau of Standards reports NBS research and test methods, and performance criteria related to the structural and ment in those disciplines of the physical and engineering sciences in environmental functions and the durability and safety charac- which the Bureau is active. These include physics, chemistry, teristics of building elements and systems. engineering, mathematics, and computer sciences. Papers cover a broad range of subjects, with major emphasis on measurement Technical Notes—Studies or reports which are complete in them- methodology and the basic technology underlying standardization. selves but restrictive in their treatment of a subject. Analogous to Also included from time to lime are survey articles on topics monographs but not so comprehensive in scope or definitive in closely related to the Bureau's technical and scientific programs. treatment of the subject area. Often serve as a vehicle for final As a special service to subscribers each issue contains complete reports of work performed at NBS under the sponsorship of other citations to all recent Bureau publications in both NBS and non- government agencies. NBS media. Issued six times a year. Annual subscription: domestic Voluntary Product Standards Developed under procedures $18; foreign $22.50. Single copy, $5.50 domestic; $6.90 foreign. — published by the Department of Commerce in Part 10, Title 15, of the Code of Federal Regulations. The standards establish NONPERIODICALS nationally recognized requirements for products, and provide all concerned interests with a basis for common understanding of the Monographs Major contributions to the technical literature — on characteristics of the products. NBS administers this program as a various subjects related to the scientific Bureau's and technical ac- supplement to the activities of the private sector standardizing tivities. organizations. Handbooks— Recommended codes of engineering and industrial Consumer Information Series— Practical information, based on practice (including safety codes) developed in cooperation with in- NBS research and experience, covering areas of interest to the con- terested industries, professional organizations, and regulatory sumer. Easily understandable language and illustrations provide bodies. useful background knowledge for shopping in today's tech- Special Publications Include proceedings of conferences spon- — nological marketplace. sored by NBS, NBS annual reports, and other special publications Order ihe above NBS puhlicaiions from: Superiniendeni of Docu- appropriate to this grouping such as wall charts, pocket cards, and nienis. Government Priming Washingion. 20402. bibliographies. Office, DC Order the following NBS publications— Fl PS and NBSIR's—from Applied Mathematics Series— Mathematical tables, manuals, and the National Technical Information Service, Springfield, VA 22161. studies of special interest to physicists, engineers, , biologists, mathematicians, computer programmers, and others Federal Information Processing Standards Publications (FIPS engaged in scientific and technical work. PUB)— Publications in this series collectively constitute the National Standard Reference Data Series— Provides quantitative Federal Information Processing Standards Register. The Register data on the physical and chemical properties of materials, com- serves as the official source of information in the Federal Govern- piled from the world's literature and critically evaluated. ment regarding standards issued by NBS pursuant to the Federal Developed under a worldwide program coordinated by NBS under Property and Administrative Services Act of 1949 as amended. the authority of the National Standard Data Act (Public Law Public Law 89-306 (79 Slat. 1127), and as implemented by Ex- 90-396). ecutive Order 11717(38 FR 12315, dated May II, 1973) and Part 6 of Title 15 CFR (Code of Federal Regulations).

NOTE: The principal publication outlet for the foregoing data is NBS Interagency Reports (NBSIR)—A special series of interim or the Journal of Physical and Chemical Reference Data (JPCRD) final reports on work performed by NBS for outside sponsors published quarterly for NBS by the American Chemical Society (both government and non-government). In general, initial dis-

(ACS) and the American Institute of Physics (AlP). Subscriptions, iribuiion is handled by the sponsor: public distribution is by the reprints, and supplements available from ACS, 1 155 Sixteenth St., National Technical Information Service , Springfield, VA 22I6I, NW, Washington, DC 20056. in paper copy or microfiche form. U.S. Department of Commerce National Bureau of Standards

Washington, D.C. 20234 Official Business

Penalty for Private Use $300