Viscosity Measurements on Gaseous Methane: Re-evaluation

Previous experimental viscosity data for gaseous methane, published by Schley et al. (Int J Thermophys 25:1623, 2004) and originally obtained using a vibrating-wire viscometer in the temperature range between 260 K and 360 K, were re-evaluated after an improved re-calibration. For this purpose, a new reference value for argon at 298.15 K and at zero density, proposed by Vogel et al. (Mol Phys 108:3335, 2010) and further updated by Hellmann (Private Communication, 2020), was applied.

Previous measurements of the viscosity for methane, carried out by Schley et al. [1] using a vibrating-wire viscometer with freely suspended weight and additional measurements of temperature T and pressure p for calculating the required density with an equation of state by Setzmann, Wagner, and de Reuck [2,3], have been reevaluated. The re-evaluation concerns the specification of the wire radius by means of an improved calibration. The re-evaluated data, given as pT values, are to be used together with other new accurate pT data to generate a prospective viscosity correlation for methane.
The re-calibration of the applied vibrating-wire viscometer was performed in that way that the radius of the wire was newly specified using previous measurements on argon [4]. The original calibration employed an experimentally based reference value of Kestin and Leidenfrost [5] nowadays considered as obsolete. For the re-calibration, we used a today accepted value for the zero-density viscosity coefficient of argon, derived by Vogel et al. [6] from an ab initio potential on the basis of the kinetic theory of dilute gases and upgraded by Hellmann [7], to 67 Page 2 of 15 be 0,Ar,298.15K = 22.5534 μPa⋅s with a standard uncertainty of ± 0.07 %. The wire radius amounts to 12.7548 μm using the new reference value for argon.
The results reported in Table II of the previous paper of Schley et al. [1] were restricted to p triples along the measured isotherms. In this new report, we include more details in order to make the information comparable to that given for the new viscosity measurements by Humberg et al. [8]. They recommended further a re-evaluation of the measurements by Schley et al. The individual points were not exactly measured at the nominal temperature of an isotherm T nom , but could be kept within small deviations from the nominal temperature each. The experimental viscosity data were adjusted to T nom values at the nominal temperature using a Taylor series expansion restricted to the first power in temperature.
For this, the experimentally determined value of the initial density dependence ( ∕ T) = (0.030 to 0.034) μPa⋅s⋅K −1 advised by Vogel [9] for methane was used. Further, it is supposed that the density values eos(T,p) , computed from the measured data for T and p using the equation of state by Setzmann, Wagner, and de Reuck [2,3], and those for the isotherms are the same. As a result of this, the pressures p T nom , eos at the nominal temperature changed and were recalculated from the densities. The improved experimental pT data of this work for the previous measurements of Schley et al. [1] on methane (six isotherms at 260 K, 280 K, 300 K, 320 K, 340 K, and 360 K) are summarized in Tables 1, 2, 3, 4, 5, and 6, in which the data are given in the sequential arrangement of the original measurements.
The experimental data of each nominal isotherm for methane were correlated as a function of the reduced density by means of a power-series representation restricted to the fourth power: Here, is the reduced density, whereas is the reduced temperature. The values of the critical density c,CH 4 and of the critical temperature T c,CH 4 are those given by Setzmann, Wagner, and de Reuck [2,3]. Weighting factors w i = 100 −2 exp,i were used in the multiple linear least-squares regression to minimize the weighted sum of squares = ∑ i w i ( cor,i − exp,i ) 2 as criterion for the quality of the representation of the considered isotherm. The coefficients i ( ) of Eq. 1 including their standard deviations s.d.
i and the weighted sum of squares for each isotherm are given in Table 7. (1) with c,CH 4 = 162.66 kg⋅m −3 , T c,CH 4 = 190.564 K.