Viscosity Measurements on Natural Gas: Re-evaluation

Previous experimental viscosity data for natural gas, 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 320 K, were re-evaluated after an improved re-calibration. For this purpose, a new reference value for the viscosity of 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, University of the Federal Armed Forces Hamburg, Hamburg, 2020), was applied. In addition, the density computed from the measured temperatures and pressures was determined using the equation of state by Kunz and Wagner (J Chem Eng Data 57:3032, 2012) instead of employing a calculation according to the International Standard ISO 12213 nowadays out of date.


Introduction
We appreciate the collaboration with Professor Span, particularly when we worked together on the viscosity correlation for ethane [1].Professor Span's group performed measurements and modeling of the viscosity of nitrogen-carbon dioxide [2] and methane-ethane mixtures [3] which we used to derive diffusion coefficients.Previous measurements of the viscosity for two samples of natural gas with different compositions, carried out by Schley et al. [4] using a vibratingwire viscometer with freely suspended weight and additional measurements of temperature T and pressure p for the calculation of the required density applying the International Standard ISO 12213 [5], have been re-evaluated.The re-evaluation 1 3 177 Page 2 of 24 concerns the determination of the wire radius by means of an improved calibration and the density calculation.In this work, the density is calculated using the equation of state by Kunz and Wagner [6].
During and since the investigation by Schley et al. [4], viscosity measurements on natural gases or on mixtures with at least three components of natural gases were performed by Nabizadeh and Mayinger [7], Assael et al. [8], Langelandsvik et al. [9], Atilhan et al. [10], Kashefi et al. [11], Jarrahian et al. [12], Nazeri et al. [13], and Al Ghafri et al. [14][15][16].The re-evaluated data of this work are reported as pT values and can be used together with the other data to update the database for natural gases.Moreover, mixture models using the current database for natural gases were discussed in the literature, e.g., by Chichester and Huber [17], by Heidaryan and Jarrahian [18], by Yang et al. [19], by Theran-Becerra et al. [20], or by Xiong et al. [21], and should enable to generate a prospective viscosity correlation for natural gases.

Re-evaluation of the Data
The re-calibration of the applied vibrating-wire viscometer was performed such that the radius of the wire was newly determined using previous measurements on argon [22].The original calibration employed an experimentally based reference value of Kestin and Leidenfrost [23] nowadays considered as obsolete.For the re-calibration, we applied a today accepted value for the zero-density viscosity coefficient of argon, derived by Vogel et al. [24] from an ab initio potential on the basis of the kinetic theory of dilute gases and upgraded by Hellmann [25], to be η 0,Ar,298.15K= 22.5534 μPa•s with a standard uncertainty of 0.07 %.The wire radius amounts to R = 12.7548 μm using the new reference value for argon.
The results reported in Tables III and IV of the previous paper of Schley et al. [4] were restricted to p triples along the measured isotherms.In this work, we include more details in order to make the information comparable to recent viscosity measurements.The individual points were not exactly measured at the nominal temperature of an isotherm T nom , but were kept within small deviations from the nominal temperature each.The experimental viscosity data were adjusted to T nom values at the nominal temperature employing a Taylor series expansion restricted to the first power in temperature.For this, the experimentally determined value of the temperature coefficient ( ∕ T) = (0.020 to 0.038) μPa•s•K −1 calculated from the fit of a double polynomial in temperature and in density to the re-evaluated experimental viscosity data for natural gas was used.In addition, it is supposed that the density values eos (T, p) computed from the measured data for p and T applying the equation of state by Kunz and Wagner [6] via their implementation in REFPROP (Lemmon et al. [26]) 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.Furthermore, it has to be stated that the original composition of the two samples natural gas H and L, given by Schley et al. in Table I of their paper including the uncertainties of the underlying gas analysis, were used in this work, but with some combinations of components to yield 20 components as input parameters for the mixture in REFPROP.Table 1 shows the mole fractions of the components of the two natural gas samples H and L that means the composition stated by Schley et al. and that used for the re-evaluation of the density and viscosity data in this work.Schley et al. [4] summed up the mole fractions of oxygen and argon, whereas Table 1 reveals that we split the mole fractions in two parts with equal values and treated separately these components.To meet the requirement of REFPROP (Lemmon et al. [26]) to limit the components to the maximum number of 20, we added the mole fractions of the xylenes (0.0006 mol% for natural gas H and 0.0010 mol% for natural gas L) to that of toluene (an additional 0.0001 mol% was added to the mole fraction of toluene in the case of natural gas H to obtain unity) and the mole fraction of neopentane to that of isopentane (0.0032 mol% for natural gas H and 0.0056 mol% for natural gas L).Moreover, the composition for the sample of natural gas L given by Schley et al. does not yield unity.Therefore, we normalized the mole fractions of all components of this sample as shown in Table 1.Finally, the molar masses of each sample were recalculated, but the numbers do not change compared with the values given by Schley et al.
The improved experimental pT data of this work for the previous measure- ments of Schley et al. [4] on the two samples of natural gas at four isotherms a In Ref. [4]      where is the reduced density, whereas is the reduced temperature.The values of the critical densities c,H and c,L as well as of the critical temperatures T c,H and T c,L were determined using REFPROP (Lemmon et al. [26]).Weighting factors w i = 100 −2 exp,i were applied 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 10 for natural gas H and in Table 11 for natural gas L.

Discussion of the Results
Humberg et al. [3] recently stated that the experimental values published by Schley et al. [4] for methane are assumed to be about 0.15 % too high due to the use of the former reference value for the viscosity of argon given by Kestin and Leidenfrost [23] for the calibration of the viscometer used by Schley et al.The same calibration value was employed for the measurements of the two natural gas samples discussed herein.Therefore, the re-evaluated values reported in this work are expected to be about 0.15 % lower than the older ones.Figures 1 and 2 show deviations of the experimental viscosity data for the two samples natural gas H and L, respectively, given by Schley et al. from the re-evaluated viscosity data reported in this work.
In connection with Figs. 1 and 2, it has to be noted that the density values of the re-evaluated data were re-calculated from the experimental data for pressure and temperature employing the equation of state by Kunz and Wagner [6] via their implementation in REFPROP (Lemmon et al. [26]).Figure 1   density for the four isotherms of natural gas sample H, does not show constant values near to − 0.15 % as expected due to the re-calibration of the radius of the wire.The isotherm at 260 K is characterized by the largest oscillations in the deviations around − 0.15 %.The oscillations decrease with increasing temperature.The maximum and minimum values occur at ≈ 75 kg ⋅ m −3 and ≈ 180 kg ⋅ m −3 , respectively.They are certainly caused by the change in the calculation of the densities from the ISO 12213 [5] to the equation of state by Kunz and Wagner.The corresponding deviations between the original Schley data and the re-evaluated data for the natural gas sample L are displayed for the four isotherms in Fig. 2. The figure shows, similar to Fig. 1, also oscillations of the deviations around − 0.15 %.These oscillations are less at lower densities (  < 100 kg ⋅ m −3 ) and comparable at higher densities (  > 120 kg ⋅ m −3 ) compared with natural gas H.The oscillations in the deviations decrease again with increasing temperature.
Both figures indicate that the differences of the original Schley data and the re-evaluated data amount to extra ± 0.2 % due to the changed density calculation dependent on the density itself.Figures 1 and 2 show that the re-calibration of the a U c,r ( ) = 0.8 % b Influenced by slip  radius of the vibrating wire results in a reduction in the viscosity of − 0.15 %, particularly at low densities.
In addition, a remark concerns the outliers in Figs. 1 and 2. The reason for these deviations are typing errors in Tables III and IV of Schley et al. [4].In Table III (natural gas H) at the isotherm 260 K, the viscosity value at ρ = 4.3077 mol•L −1 (ρ = 77.4kg•m −3 ) should read 12.315 μPa•s instead of 12.335 μPa•s and in Table IV (natural gas L) at the isotherm 320 K, the viscosity value at ρ = 5.3865 mol•L −1 (ρ = 99.3 kg•m −3 ) should read 15.910 μPa•s instead of 15.900 μPa•s.
Moreover, we compare our re-evaluated data with the extended corresponding states (ECS) model [27] implemented in REFPROP [26].This is done in Figs. 3 and 4, which present deviations between the re-evaluated data for the two samples of natural gas H and L and the calculated values for the ECS model.The zero line corresponds to the ECS model developed by Ely and Hanley.The calculation of the ECS values results from the experimental temperature and pressure data.We note that for natural gas sample H as well as sample L no viscosity values applying the ECS model could be calculated at 260 K in low density ranges, for sample H at 4 < ∕kg•m −3 < 84 and for sample L at 2 < ∕kg•m −3 < 103, respectively.This is due to the retrograde behavior of the samples at low temperatures.Figure 3, in which the deviations of the re-evaluated data for the natural gas sample H from the values calculated for the ECS model [27] are illustrated, makes evident that for the higher isotherms 280 K, 300 K, and 320 K the deviations are between +0.09 % at the lowest densities and ≈ − 0.45 % up to high densities.For the isotherm 260 K, the largest deviation amounts to +0.84 % at = 85 kg ⋅ m −3 .The deviations decrease with increasing density for this isotherm, pass through zero, and become negative with a minimum of − 0.45 % at = 181 kg ⋅ m −3 .The deviations at  > 181 kg ⋅ m −3 decrease to − 0.15 % on the average, corresponding to the re- calibration of the wire.
Figure 4, which displays the deviations of the re-evaluated data for the natural gas sample L from the values calculated for the ECS model [27], shows that larger deviations occur compared to natural gas sample H.For the isotherms 280 K, 300 K, and 320 K, the deviations at the lowest densities amount to + 0.80 %, +1.23 %, and +1.44 %, respectively.They decrease with increasing den- sity down to +0.60 % (280 K) and to ≈ 1.0 % (300 K and 320 K) at ≈ 100 kg ⋅ m −3 , followed by a further decrease to +0.67 % (320 K) and to +0.33 % (300 K) at the highest density each.The deviation of the 260 K isotherm at = 104 kg ⋅ m −3 amounts to +0.63 %, whereas the further deviations and those for the iso- therm 280 K show the similar course with increasing density.They pass through a U c,r ( ) = 0.8 % b Influenced by slip zero at ≈ 180 kg ⋅ m −3 and become negative (− 0.13 % for 280 K and − 0.18 % for 260 K at the highest density each).

Conclusion
A new reference value for argon at 298.15 K and at zero density, which was derived by Vogel et al. [24] and further updated by Hellmann [25], was used to re-evaluate viscosity measurements of Schley et al. [4] for two natural gas samples H (highly caloric) and L (lowly caloric) for the four isotherms each.Schley et al. carried out relative measurements using a vibrating-wire viscometer and calibrated their measuring device employing a nowadays outdated reference value of Kestin and Leidenfrost [23].Furthermore, Schley et al. [4] applied the international standard ISO 12213 [5] for calculating the densities from the experimental pressure and temperature data.In this work, the densities were computed employing the equation of state    [27] using REFPROP [26] for the four isotherms, as a function of density .Data: ○ , 260 K; ▵ , 280 K; ◻ , 300 K; and ◊ , 320 K of Kunz and Wagner [6] from the pressure and temperature data reported by Schley et al.
The re-evaluated experimental data of each nominal isotherm for the two natural gas samples H and L were correlated as a function of the reduced density by means of a power-series representation restricted to the fourth power.Finally, the re-evaluated viscosity data were discussed in comparison with values calculated using the approach of the ECS model by Ely and Hanley [27] via its implementation in REFPROP [26].The re-evaluated data of this work reported as pT data should be employed to update the database for natural gases.
Author Contributions SH performed the calculations regarding the re-evaluation and prepared the tables and figures.All authors wrote and reviewed the manuscript.
Funding Open Access funding enabled and organized by Projekt DEAL.There is no funding to declare.

Declarations
Conflict of interest The authors declare that they have no conflicts of interest.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material.If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission Fig. 4 Comparison of the re-evaluated viscosity data of this work exp for the sample natural gas L with values REFPROP calculated from the approach of the ECS model developed by Ely and Hanley [27] using REFPROP [26] for the four isotherms, as a function of density .Data: ○ , 260 K; ▵ , 280 K; ◻ , 300 K; and ◊ , 320 K separately treated b In Ref. [4] treated together with argon c In this work separately treated d In this work treated together with toluene e In this work treated together with isopentane

a
U c,r ( ) = 1.0 % b Influenced by slip with j = H for natural gas H, j = L for natural gas L, and c,H = 208.45kg⋅m −3 , T c,H = 213.26K for natural gas H, c,L = 186.21kg⋅m −3 , T c,L = 191.88K for natural gas L, a U c,r ( ) = 0.8 % b Influenced by slip , which illustrates the deviations between the original Schley data and the re-evaluated data as function of a U c,r ( ) = 0.8 % b Influenced by slip

a
U c,r ( ) = 1.0 % b Influenced by slip

a
U c,r ( ) = 0.8 % b Influenced by slip

Fig. 1 3 177 Page 22 of 24 Fig. 2 Fig. 3
Fig.1Comparison of the re-evaluated viscosity data of this work exp,re−eval for the sample natural gas H with the data exp,Schley for the four isotherms, as a function of density .Data: ○ , 260 K; ▵ , 280 K; ◻ , 300 K; and ◊ , 320 K

Table 1
[4]ponents and their mole fractions of the natural gas samples H and L given by Schley et al.[4]and used for the re-evaluated viscosity measurements of this work

Table 2
Re-evaluated experimental pT data for natural gas H at 260 K

Table 3
Re-evaluated experimental pT data for natural gas H at 280 K

Table 4
Re-evaluated experimental pT data for natural gas H at 300 K

Table 5
Re-evaluated experimental pT data for natural gas H at 320 K

Table 6
Re-evaluated experimental pT data for natural gas L at 260 K

Table 7
Re-evaluated experimental pT data for natural gas L at 280 K

Table 8
Re-evaluated experimental pT data for natural gas L at 300 K

Table 9
Re-evaluated experimental pT data for natural gas L at 320 K a U c,r ( ) = 0.8 % b Influenced by slip

Table 10
Coefficients of Eq. 1 for the re-evaluated viscosity measurements on the sample of natural gas H

Table 11
Coefficients of Eq. 1 for the re-evaluated viscosity measurements on the sample of natural gas L