1 Introduction

Liquefied biomethane (LBM), also known as liquefied biogas (LBG), bio liquefied natural gas (bio-LNG) or renewable LNG, is a climate-neutral fuel with a comparatively large energy density. LBM essentially combines the benefits of LNG and carbon neutrality. In a recent study, Gustafsson and Svensson report that using LNG as a transport fuel involves less greenhouse gas emissions only under certain circumstances, while “there are clear environmental incentives to use LBM” [1]. However, for LBM, it currently seems problematic “to compete economically with the price of LNG, due to the higher specific production costs”. Nevertheless, in the future, there appears to be a realistic chance that LBG gets competitively viable as the supply chain is evidently simpler than the conventional one for LNG [2, 3]. Since the utilization of LNG and LBG as transport fuel is part of the European Union clean fuel strategy, the Joint Research Project “Metrological support for LNG and LBG as transport fuel” was funded as part of the European Metrology Programme for Innovation and Research (EMPIR) with the overall aim to enable the large-scale roll-out of LNG and LBG as a transport fuel. This includes to achieve smaller uncertainties in the determination of the LNG’s and LBG’s energy quantity during custody transfer, and here, one of the key properties is the density.

In our last paper [4], we presented density measurements of a (0.99 methane + 0.01 butane) mixture and a (0.98 methane + 0.02 isopentane) mixture over the temperature range from (100 to 160) K at pressures up to 10.8 MPa. This was done to test the performance of the EOS-LNG of Thol et al. [5], which is a new fundamental equation of state particularly developed for LNG mixtures that includes binary specific departure functions for butanes and pentanes; for this development, mainly the data of Lentner et al. [6] was used. Within the work of Eckmann et al. [4], we demonstrated that the EOS-LNG could improve the density prediction for LNG containing butane and isopentane compared to the GERG-2008 equation of state of Kunz and Wagner [7, 8], and we delivered accurate data for further improvement of the EOS-LNG.

Here, we report the results of accurate density measurements of two liquefied biomethane-like mixtures over the temperature range from (100 to 180) K at pressures up to 9.0 MPa, utilizing a special single-sinker densimeter for cryogenic liquid mixtures [4, 6, 9,10,11,12]. For both mixtures, the methane mole fraction is larger than 0.95, and the further components are nitrogen, hydrogen, and oxygen; one of the mixtures also contains carbon dioxide. Such measurements of biomethane-like mixtures, in particular in the homogeneous liquid phase and the supercritical region, at low temperatures and over a wide pressure range are unique so far. To the best of our knowledge, measurements like reported here have not been carried out elsewhere before as already discussed in our earlier papers [4, 6, 9,10,11,12]; other authors only measured data for saturated-liquid densities of binary mixtures and LNG-like mixtures [13,14,15,16,17,18,19,20,21] but not in the homogeneous liquid region below the cricondenbar pressure. For the two mixtures studied in the present work, however, no measurements of saturated-liquid densities could be found in literature. Hence, the new (p, ρ, T, x) data that we measured can be used to test the performance of the GERG-2008 equation of state [7, 8] regarding these two methane-rich multicomponent mixtures with compositions typical for biomethane. (We note that the EOS-LNG was solely improved towards the methane + butanes and pentanes binaries, thus, there is no difference in the performance of EOS-LNG and the GERG-2008 equation of state for the mixtures investigated in this work.)

2 Experimental Section

2.1 Apparatus Description

A single-sinker magnetic-suspension densimeter, particularly developed for cryogenic liquid mixtures, was utilized; it covers a temperature range from (90 to 300) K at pressures up to 12 MPa. The densimeter’s design, temperature and pressure measurement as well as the implementation of a special “VLE-cell” were described in detail previously [9,10,11]. A schematic diagram of the densimeter, presented by Richter et al. [11], is given in the Online Resource 1. In one of our latest papers [6], we presented a few improvements of the core apparatus to reduce diffusion effects and the force-transmission error (FTE) [22, 23]. Here, we summarize the description of the apparatus presented by Richter et al. [11]. Overviews of this general type of densimeter were provided by Wagner and Kleinrahm [24] as well as by McLinden [25].

The single-sinker method basically allows an absolute determination of the fluid density. This method is applied in conjunction with a magnetic-suspension coupling and a load compensation mechanism (differential method). A sinker of known volume \({V}_{\text{S}}\left(T,p\right)\) and known mass \({m}_{\text{S}}\) (in the present case: a single-crystal silicon, \({m}_{\text{S}}\) ≈ 60.95 g, \({V}_{\text{S}}\) ≈ 26.17 cm3, \({\rho }_{\text{S}}\) ≈ 2.329 g ⋅ cm–3) is weighed while immersed in fluid of interest inside a pressure-tight measuring cell. Thus, the result of weighing the sinker located in the fluid, \({m}_{\text{S},\text{fluid}}^{*}\), is the difference between the mass of the sinker and the buoyancy of the fluid:

$${m}_{\text{S},\text{fluid}}^{*}={m}_{\text{S}}-{\rho }_{\text{fluid}}\cdot {V}_{\text{S}}\left(T,p\right),$$
(1)

where \({\rho }_{\text{fluid}}\) denotes the fluid density. When weighing the sinker inside the evacuated measuring cell via the magnetic-suspension coupling, the weighing result is not the mass of the sinker, \({m}_{\text{S}}\), but a slightly different result, \({m}_{\text{S},\text{vac}}^{*}\), which is the result of a small FTE of the magnetic-suspension coupling [22, 23]. Rearranging Eq. (1) yields the fluid density:

$${\rho }_{\text{fluid}}=\frac{{m}_{\text{S},\text{vac}}^{*}-{m}_{\text{S},\text{fluid}}^{*}}{{V}_{\text{S}}(T,p)}.$$
(2)

Equation (2) actually requires additional terms since essential details of the measurement procedure (e.g., the correction of the FTE) as discussed by Richter et al. (see reference [11], Appendix A) have to be taken in account. Moreover, the problem of the FTE due to the magnetic-suspension coupling is described in detail by Kleinrahm et al. [23].

For the measurement of fluid density, the sinker is connected to an analytical balance (resolution: 0.01 mg) employing an appropriate coupling/decoupling device. Gravity and buoyancy forces acting on the sinker are transmitted to the balance via the magnetic-suspension coupling, thus, isolating the fluid sample (which may be at high pressure and very low temperature) from the balance, which is at ambient conditions. For compensation of the balance’s zero-point drift, the weighing procedure considers the small drift of the balance reading in the tare position by subtracting it from the balance reading in the measuring position. The balance is operated near its tare point using a load compensation mechanism for reduction of possible errors of the balance due to changes in the slope of the characteristic curve over the weighing range.

2.2 Experimental Material

The two synthetic biomethane mixtures were gravimetrically prepared and delivered in steel cylinders with an internal volume of 50 dm3 from a commercial supplier (Praxair Inc.). The desired compositions of the mixtures to be measured were previously specified by the collaborators within the EMPIR project. The gas mixtures have been produced from high-purity starting materials, to minimize effects from impurities on the density measurements. The mole fractions of the main components have been determined by comparison with the Dutch primary measurement standards in accordance with ISO 6143 [26] using gas chromatography with a thermal conductivity detector; a description of the instrument configuration and measurement procedure can be found elsewhere [27]. This procedure establishes the metrological traceability of the measured mole fractions as well as the smallest uncertainties, both required to obtain credible reference data. Results of the analyses, which are metrologically traceable to primary and internationally accepted measurement standards, are listed in Table 1. The stated uncertainties in composition were reported by the Dutch Van Swinden Laboratorium (VSL) in the calibration certificates as expanded uncertainties (k = 2), whereby the underlying standard uncertainties were determined according to the Guide to the expression of uncertainty in measurement (GUM) [28]. Further information on the mixtures is summarized in Tables 1 and 2. These are relevant for the experimental procedures described in Sect. 2.3.

Table 1 Composition, standard uncertainties (k = 1) for each component, and molar mass M of the studied mixtures
Table 2 Additional information on the studied mixtures relevant for the supercritical liquefaction procedure

To prevent any kind of change of the mixture composition (e.g., due to phase separation), the sample was handled very carefully. Filling the sample into the measuring cell of the densimeter with the correct composition was essential. Hence, the sample cylinder was prepared according to the following procedure: (1) Rolling the sample cylinder for at least two hours to re-homogenize the gas mixture. (2) Heating the cylinder at the bottom part for at least 3 h using a heating jacket to obtain vortices inside the sample cylinder for homogenizing the gas. (3) Filling the gas mixture into the well-evacuated system to a pressure of about 0.2 MPa through the evacuated filling line and leaving the sample with a residence time of about two minutes before evacuating the apparatus again. Step (3) was repeated three times before the final filling. Thereby, residual gas from previously studied samples is removed to prevent an unwanted compositional change of the new sample filled into the densimeter.

2.3 Experimental Procedures

The details of filling the densimeter and the basic procedure of operating the apparatus are reported by Richter et al. [11]. We already used this procedure for our previous density measurements on synthetic LNG mixtures and methane-rich binary mixtures [4, 6, 11, 12]. In the present work, this established procedure was applied for density measurements of two synthetic biomethane-like mixtures. Here, it is important to mention that our densimeter involves the application of a special VLE-cell, which serves as a buffer for the unavoidable phase transition (vapor–liquid) between the very cold measuring cell and the pressure measurement circuit that is kept at a constant temperature of about 313.15 K. The temperatures of the measuring cell and the VLE-cell can be controlled independently of each other at different set points. Since all areas of the measurement system are interconnected, the pressure is everywhere the same (apart from a pressure head correction).

At first, we filled the densimeter at ambient temperature to a pressure higher than the cricondenbar pressure (e.g., pfill = 9.0 MPa for present measurements); for example, see Fig. 1. Then, the measuring cell and the VLE-cell were cooled simultaneously at constant pressure by continuously adding sample to the system until the VLE-cell had reached a slightly subcritical temperature or at least a temperature considerably below the temperature at the cricondenbar pressure pccp. We maintained this temperature of the VLE-cell at a constant value, while the measuring cell was cooled further to the desired temperature and finally controlled to maintain a constant value. Thereby, the filling procedure was completed, and the first density value at p > pccp was measured. The composition determined for the gas mixture in the sample cylinder before filling the densimeter remained unchanged for the following density measurements in the homogeneous liquid phase.

Fig. 1
figure 1

Left: Principle of the filling and liquefaction procedure shown in a p, T-diagram for a methane-rich gas mixture containing hydrogen. As an example, the points to be measured at T = 100 K are plotted. For illustration purposes, a composition was chosen that differs from the mixtures investigated in this work (compare with Fig. S3 in the Online Resource 1). The phase envelope was calculated using the GERG-2008 equation of state of Kunz and Wagner [7, 8]. Right: Schematic presentation of the measurement system consisting of the measuring cell (M), VLE-cell (V) and pressure measurement system (P)

After the first density measurement above pccp (point 2 in Fig. 1), the pressure was reduced by venting gaseous sample from the system. This pressure reduction, at approximately constant temperature of the VLE-cell, results in a vapor–liquid equilibrium forming in the VLE-cell. The venting was stopped, when a liquid volume fraction of about 30 % was established in the VLE-cell (point 3 in Fig. 1). The respective pressure for the set-point temperature of the VLE-cell can be approximately calculated with the GERG-2008 equation of state [7, 8] in advance. At this point, the next density measurement was carried out.

To measure densities at further (p, T) state points in the homogeneous liquid phase, the pressure in the system was no longer adjusted by venting sample from the system but via reducing the temperature of the VLE-cell. Hence, from this point on, the total mass of fluid inside the system remained constant. By reducing the temperature of the VLE-cell, the pressure in the entire system was decreased accordingly to measure the points 4 to 6 (see Fig. 1). The last point 6 was always measured close to the calculated saturated-liquid line. The respective VLE-cell temperature for a desired pressure was determined via mass balance calculations using the GERG-2008 equation of state [7, 8]. The determined temperature was always just a little above the saturation temperature for the desired pressure and mixture composition, resulting in a liquid volume fraction in the VLE-cell between (30 and 70)  %. For each investigated isotherm, we filled the densimeter separately.

2.4 Uncertainty in Density Measurement

The uncertainty in density measurement was determined in accordance with the GUM [28]. Assuming that there is no correlation between the input quantities, the expanded combined uncertainty U for the determination of cryogenic liquid densities using the above described densimeter can be calculated from:

$$U\left( \rho \right) = k \cdot u_{c} \left( {\rho \left( {T,p,\user2{x}} \right)} \right) = k \cdot \left[ {\begin{array}{*{20}c} {u\left( \rho \right)^{2} + \left( {\left( {\frac{{\partial \rho }}{{\partial p}}} \right)_{T} \cdot u\left( p \right)} \right)^{2} + \left( {\left( {\frac{{\partial \rho }}{{\partial T}}} \right)_{p} \cdot u\left( T \right)} \right)^{2} } \\ { + {\mkern 1mu} u\left( {\rho \left( \user2{x} \right)} \right)^{2} + u\left( {\rho _{{{\text{repro}}}} } \right)^{2} + u\left( {\rho _{{{\text{corr}}}} } \right)^{2} } \\ \end{array} } \right]^{{{\raise0.7ex\hbox{$1$} \!\mathord{\left/ {\vphantom {1 2}}\right.\kern-\nulldelimiterspace} \!\lower0.7ex\hbox{$2$}}}} ,$$
(3)

where \(u(\rho )\), \(u(T)\) and \(u(p)\) denote the standard uncertainties in density, temperature and pressure, respectively; \(u(\rho ({\varvec{x}}))\) corresponds to the standard uncertainty in density resulting from the uncertainty of the gas composition; \(u({\rho }_{\text{repro}})\) accounts for an additional uncertainty from the reproducibility of our measurements; and \(u({\rho }_{\text{corr}})\) takes into account the uncertainty of the correction of the force-transmission error. A detailed description of the uncertainty evaluation has been reported in previous works [6, 11, 12]. The expanded combined uncertainty in measurement was estimated for each measured state point. Table 3 shows the contributions to the relative combined expanded uncertainty in density of an exemplary density measurement of the LBG I mixture at T = 140 K and p = 3.02424 MPa.

Table 3 Uncertainty budget for the density measurements. As an example, the uncertainty was calculated for the LBG I mixture at T = 140.000 K, p = 3.02424 MPa, and ρexp = 386.559 kg·m–3 (see Table 4)
Table 4 Experimental densitiesa for two synthetic biomethane-like mixtures (compositions see Table 1) and relative deviations of experimental densities ρexp from densities ρGERG calculated with the GERG-2008 equation of state of Kunz and Wagner [7, 8] (as implemented in the TREND 4.0 software package [33]), where p is the pressure, T is the temperature (ITS-90), and 100 (U(ρ)/ρ) is the relative expanded combined uncertainty (k = 2). The measurements were carried out in the homogenous liquid phase and the supercritical region

One of the main contributions (up to 68 % and 82 % for LBG I and LBG II, respectively) to the combined standard uncertainty in density originates from the uncertainty in the gas composition \(u(\rho ({\varvec{x}}))\). Here, it has to be considered that variations in composition of the individual components are correlated, which requires to consider the covariances between the mole fractions in the uncertainty evaluation. As described in ISO 6976 [31], the combined standard uncertainty in density due to the uncertainty in composition can be calculated by

$$u{\left(\rho \left(x\right)\right)}^{2}={\sum }_{i=1}^{N}\sum_{j=1}^{N}\frac{\partial \rho }{\partial {x}_{i}}\cdot \left({x}_{i}\right)\cdot r\left({x}_{i}, {x}_{j}\right)\cdot u\left({x}_{j}\right)\cdot \frac{\partial \rho }{\partial {x}_{j}},$$
(4)

where \(N\) is the total number of components in the respective mixture (see Table 1) and \(\partial \rho /\partial {x}_{i}\) are the sensitivity coefficients for each respective component. The values for \(\partial \rho /\partial {x}_{i}\) were calculated with the GERG-2008 equation of state [7, 8] for each specific state point measured. The correlation coefficients \(r({x}_{i},{x}_{j})\) for the components \(i\) and \(j\) were calculated according to ISO 14912 [32] by

$$r\left({x}_{i}, {x}_{j}\right)=\frac{u\left({x}_{i}, {x}_{j}\right)}{u\left({x}_{i}\right)\cdot u\left({x}_{j}\right)}$$
(5)

with

$$u\left({x}_{i}, {x}_{j}\right)=-{x}_{i}{u}^{2}\left({x}_{j}\right)-{x}_{j}{u}^{2}\left({x}_{i}\right)+{x}_{i}{x}_{j}\sum_{k=1}^{N}{u}^{2}({x}_{k})/{\left(\sum_{k=1}^{N}{x}_{k}\right)}^{2},$$
(6)

where \(u\left({x}_{i}, {x}_{j}\right)\) is the covariance matrix and \(u\left({x}_{i}\right)\) are the normalized uncertainties of the components in mole fraction (see Table 1).

3 Results and Discussion

The densities of two selected synthetic biomethane-like mixtures were measured in the homogeneous liquid phase and the supercritical region. Furthermore, saturated-liquid densities were determined.

3.1 Results for Homogeneous Liquid and Supercritical Densities

For LBG I, the densities were measured along isotherms at T = (100, 120, 140 and 160) K and for LBG II at T = (140, 160 and 180) K; six state points were investigated along each isotherm. Starting always with the supercritical filling pressure and ending close to the saturated-liquid pressure, a pressure range from about (1.5 to 9.0) MPa was covered. To avoid a change in composition for the liquefied samples, the filling of the apparatus was carried out at pressures at least 0.5 MPa above the cricondenbar pressure pccp (see Table 2). Furthermore, to avoid vaporization of the mixture in the measuring cell, all measurements in the liquid phase were carried out at pressures at least 0.1 MPa above the saturation pressure psat,GERG, calculated with the GERG-2008 equation of state [7, 8]. Two or three replicates were measured at each (p, T) state point, and in total 24 (p, ρ, T, x) data are reported for LBG I and 6 for LBG II.

The density measurements of the LBG II mixture at T = (140 and 160) K revealed that parts of the carbon dioxide contained in the mixture solidified. This assumption could be verified by a p, T-phase diagram of the LBG II mixture, which is attached in the Online Resource 1. The solidification resulted in a change in composition of the studied liquid, and the measured densities show large deviations from densities calculated with the GERG-2008 equation of state [7, 8], namely about − 1.9 % at T = 160 K and − 4.4 % at T = 140 K. Moreover, during these density measurements, the analytical balance readings showed a strong scatter of the values and a temporal drift. Based on the measured densities, we approximately estimated a change of the carbon dioxide mole fraction in the liquid phase from originally about 2.5 % to about (1.6 and 0.3) % at T = (160 and 140) K, respectively. Hence, only the measurements along the isotherm at T = 180 K yielded reliable results.

The experimental results of the (p, ρ, T, x) measurements of the LBG I and the LBG II mixture are listed in Table 4 together with the corresponding state point uncertainties as discussed in Sect. 2.4. The relative expanded combined uncertainty (k = 2) in density was estimated to be equal to or less than 0.027 % for LBG I and between (0.046 and 0.054)  % for LBG II. The difference between these two uncertainties is mainly caused by the higher contribution of the uncertainty in composition of LBG II; i.e., the consideration of five instead of four mixture components.

We would like to mention that a new equation of state for liquefied natural gases (LNG) was recently published by Thol et al. [5]; it is called EOS-LNG. This model is based on the GERG-2008 equation of state [7, 8] and enhances this equation with respect to the prediction of mixtures containing butanes and pentanes. Since these components were not contained in the mixtures investigated within this work, the GERG-2008 equation of state [7, 8] and the EOS-LNG [5] yield the same densities. Therefore, the GERG-2008 equation of state [7, 8] as an established model in the LNG community is used for the comparisons within this section.

The relative deviations of the experimental densities from values calculated with the GERG-2008 equation of state [4, 5] (as implemented in the TREND 4.0 software package [33]) are also reported in Table 4. In Fig. 2, the relative deviations of the experimental densities from values calculated with the GERG-2008 equation of state are plotted versus pressure. The relative deviations for the LBG I mixture are less than ± 0.01 % for T = (100, 120 and 140) K but increase up to − 0.072 % at T = 160 K. The experimental densities of the LBG II mixture at T = 180 K deviate by up to − 0.1627 % from values calculated with the GERG-2008 equation of state, with increasing deviations towards lower pressures.

Fig. 2
figure 2

Relative deviations of experimental densities \({\rho }_{\text{exp}}\) for from densities \({\rho }_{\text{GERG}}\) calculated with the GERG-2008 equation of state [7, 8] (zero line). (pρTx) data measured in the present work. LBG I: ◯, T = 100 K; ✕, T = 120 K; ☐, T = 140 K; + , T = 160 K. LBG II: △, T = 180 K. The compositions of the LBG mixtures are given in Table 1

The authors of the GERG-2008 report uncertainties for binary systems only [7, 8]. Based on these uncertainties, we estimated an uncertainty of the GERG-2008 equation of state of 0.3 % for the calculated densities of the LBG I and II mixtures in the liquid phase. Hence, the GERG-2008 equation of state describes the experimental data clearly within this uncertainty.

3.2 Determination of Saturated-Liquid Densities

For both LBG mixtures, saturated-liquid densities in the temperature range from (100 to 180) K were determined by extrapolation of the experimental densities in the homogenous liquid region to the mixture’s vapor pressure. The results are listed in Table 5. The influence of the extrapolation on the uncertainty of the saturated-liquid densities is relatively small in most cases, except for the higher temperatures T ≥ 160 K.

Table 5 Saturated-liquid densitiesa ρsat,exp for the two synthetic biomethane-like mixtures (compositions see Table 1) and their relative deviations from densities ρsat,GERG calculated with the GERG-2008 equation of state of Kunz and Wagner [7, 8] (as implemented in the TREND 4.0 software package [33]), where psat is the vapor pressurea, T is the temperature (ITS-90), and 100 (U(ρ)/ρ) is the relative expanded combined uncertainty (k = 2)
Table 6 Experimentally determined saturated-liquid pressures psat,exp for LBG I (composition see Table 1) and their absolute deviations from saturated-liquid pressures psat,GERG calculated with the GERG-2008 equation of state of Kunz and Wagner [7, 8], where U(psat,exp) is the absolute expanded uncertainty (k = 2) of the experimentally determined saturate-liquid pressures

The relative deviations of the experimental saturated-liquid densities of the two LBG mixtures from values calculated with the GERG-2008 equation of state [7, 8] are plotted in Fig. 3. The numerical values of these relative deviations are listed in Table 5. The uncertainty of the saturated-liquid densities of the two LBG mixtures calculated with the GERG-2008 equation of state was estimated by us to be 0.3 % (see Sect. 3.1). Hence, the GERG-2008 equation of state describes the new experimental saturated-liquid densities clearly within this uncertainty.

Fig. 3
figure 3

Relative deviations of experimental saturated-liquid densities ρsat,exp for two biomethane-like mixtures from densities ρsat,GERG calculated with the GERG-2008 equation of state [7, 8] (zero line). (psat, ρsat, T, x) data of the present work: LBG I: ◯; BG II: △. The uncertainties of our new saturated-liquid density data are illustrated with error bars (also see Table 5). For the LBG II mixture, no error bars are plotted since they are outside the scaling of the ordinate. The compositions of the LBG mixtures are given in Table 1

For the determination of the saturated-liquid densities, a cubic equation was fitted to the six experimental densities on each isotherm in the homogenous liquid phase. Then, these equations were extrapolated to the mixture’s vapor pressure. The vapor pressures needed for this extrapolation were usually calculated with the GERG-2008 equation of state [7, 8]. For the LBG II mixture at T = 180 K it was psat = 3.763 MPa; its relative expanded uncertainty (k = 2) was estimated by us to be 3 %, based on ref. [7]. For the LBG I mixture, however, the saturated-liquid pressures calculated with the GERG-2008 equation of state for T < 160 K are obviously incorrect. In this region, the saturated-liquid line in a ρ, T-diagram shows a curvature from lower to higher pressures; see Fig. 4. This curvature is also depicted in a p, T- diagram in Fig. S4 in the Online Resource 1. The reason for this unusual behavior is presumably the component hydrogen in the LBG mixture.

Fig. 4
figure 4

ρ, p-diagram of LBG I (see Table 1 for composition). The phase boundary and the isotherms were calculated with the GERG-2008 equation of state of Kunz and Wagner [7, 8]; , (p, ρ, T, x) state points measured in the present work; , saturated-liquid densities determined by extrapolation to the experimentally determined saturated-liquid pressure; ●, critical point

Therefore, we determined the saturated-liquid pressures of LBG I ourselves with our densimeter. The measurement procedure is explained in the Online Resource 1. The results are listed in Table 6. Since our densimeter is not suitable for the measurement of saturated-liquid pressures of fluid mixtures, the uncertainty of the results is relatively large; see Table 6. Moreover, saturation pressures calculated with the GERG-2008 equation of state [7, 8] and the difference of the experimental pressures from the calculated pressures are also given in the table.

The new (p, ρ, T, x) data measured in the homogeneous liquid region are depicted in a p, T-diagram in Fig. 4. Furthermore, the saturated-liquid densities, which were determined by extrapolation of the six experimental densities on each isotherm in the homogenous liquid region to the vapor pressure, are also shown.

It should be mentioned here that the GERG-2008 equation of state [7, 8] also yielded reasonable values for saturated-liquid densities ρsat,GERG(T), even though the respective saturated-liquid pressures are higher than the experimentally determined values (see Fig. 4 and Table 6). Therefore, such calculated values were used in Table 5 and Fig. 3 for comparison with the new experimental saturated-liquid densities.

4 Conclusion

Accurate density measurements of two biomethane-like mixtures, LBG I and LBG II, were carried out along isotherms from T = (100 to 180) K at pressures up to 9.0 MPa, utilizing a special single-sinker densimeter for cryogenic liquid mixtures. Both mixtures consist of more than 95 mol- % methane and the three further components nitrogen, hydrogen, and oxygen; the LBG II mixture additionally contains carbon dioxide. The densities of the LBG I mixture were measured at T = (100, 120, 140, and 160) K at six state points on each isotherm between the vapor pressure and p = 6.6 MPa. The LBG II mixture could only be measured at T = 180 K at six state points between the vapor pressure and p = 9.0 MPa. The reason was presumably the solidification of about (88 and 36)  % of carbon dioxide at T = (140 and 160) K, respectively, of the original content of about 2.5 mol % in this mixture. A problem arose when determining the saturated-liquid densities of the LBG I mixture. These values were determined by extrapolation of the six density state points on each isotherm in the homogeneous liquid region to the mixture’s vapor pressure. Since the vapor pressures calculated with the GERG-2008 equation of state [7, 8] were obviously incorrect for T ≤ 160 K, we measured the vapor pressures by ourselves with our densimeter.

The experimental results were compared with the GERG-2008 equation of state of Kunz and Wagner [7, 8]. The relative deviations of the new experimental values from densities calculated with the GERG-2008 equation of state are less than ± 0.01 % in the homogeneous liquid region of the LBG I mixture at T = (100, 120 and 140) K, and less than 0.08 % at T = 160 K. For the LBG II mixture, the relative deviations in the homogeneous liquid region are between – (0.02 and 0.17)  % at T = 180 K. Thus, the GERG-2008 equation of state is able to describe the new experimental densities of both LBG mixtures within the uncertainty of 0.3 %, where this uncertainty was estimated by us based on the reported uncertainty by the authors of the GERG-2008 equation of state [7, 8] for binary mixtures.