Thermodynamic and Transport Properties of Biomass-Derived Furfural, Furfuryl Alcohol and Their Mixtures

The limited reserves and well-known disadvantages of using fossil energy sources have increased the need for appropriate renewable substitutes in the production of various chemicals and materials. Biomass has been shown to be worthy of attention since it can be converted to biofuels and value-added chemicals relatively easily. The design of biomass valorisation process requires knowledge on the thermodynamic behaviour of the biomass-derived compounds, such as furfural and furfuryl alcohol. The thermodynamic and transport properties of the binary system furfural + furfuryl alcohol were studied at various temperatures and pressures. Density, speed of sound and refractive index were measured in the temperature range T = (288.15–345.15) K and viscosity was measured at temperatures up to 373.15 K, all at atmospheric pressure. Further, the density of pure components was obtained in the temperature range (293.15–413.15) K for furfural and (293.15–373.15) K for furfuryl alcohol at pressures up to 60.0 MPa. The obtained density values were correlated using the modified Tammann–Tait equation with an average absolute deviation lower than 0.009% for furfural and furfuryl alcohol. The optimised parameters were used for the calculation of the isothermal compressibility, the isobaric thermal expansivity, the internal pressure and the isobaric and isochoric specific heat capacities. The reported data are a valuable source of information for the further application of the investigated compounds.


Introduction
Due to increasing environmental pollution caused by the use of fossil fuels and the depletion of their reserves, there is an increasing need for renewable energy sources that would meet the basic principles of "green" chemistry, as defined by Anastas and Warner [1].The transition from today's economy based on fossil fuels to a sustainable biobased economy is a great challenge.Biomass has appeared as a promising, carbon-neutral substitute for fossil fuels in terms of energy production and the production of various chemicals and materials.In particular, lignocellulosic biomass has gained significant attention due to its low cost, less waste and being the most abundant biopolymer in nature [2].It mainly consists of three components, i.e., cellulose and hemicellulose that compose of polysaccharides, and lignin, an aromatic polymer [3].Lignocellulose has shown a great potential in production of platform chemicals that can be further converted to added-value chemicals by a hydrodeoxygenation process (HDO).The design of biomass valorisation processes, as well as the collection and handling the obtained products, requires the knowledge on thermodynamic behaviour of all components present in the process in order to maximise efficiency and avoid unnecessary energy consumption.Furthermore, Energies 2021, 14, 7769 2 of 18 HDO consists of complex reactions that overlap and entail accurate microkinetic models involving transfer phenomena and thermodynamic properties of the compounds [4].
Some promising compounds produced from biomass are furfural and furfuryl alcohol.Furfural is one of the furan derivatives formed from the hemicellulose part of lignocellulose and is a very important platform chemical that can be transformed to widely used chemicals and fuels [5].Nowadays, its production is based on hydrothermal conversion of biomass, particularly, pentosan sugars are hydrolysed to monosaccharides which are further transferred to furfural by dehydration [5][6][7][8].The efficiency of the conversion depends strongly on reaction condition, used catalyst and solvent, and the optimisation of these three parameters is subject of many research [7,9].Furfuryl alcohol is often used in the production of resins, pharmaceuticals and fragrances, and it is typically produced by the selective hydrogenation of furfural [5,10].The world annual production of furfural has been estimated to over 300 ktons and about 65% of that is used in the production of furfuryl alcohol [5].
Despite the great range of the application of furfural and furfuryl alcohol, such as solvents, biofuels, an alternative to the production of antacids, fertilizers, plastics and paints [11], the thermodynamic properties of furfural, furfuryl alcohol, and especially their mixtures are not enough studied.Lomba et al. [12] performed thorough study of physicochemical properties of biomass-derived green solvents.They modelled density, refractive index, speed of sound, surface tension, dynamic viscosity, static permittivity and vapour pressure of furfural and furfuryl alcohol at temperatures (279.15-338.15)K and seemingly at atmospheric pressure.Bendiaf et al. [13] studied density and speed of sound of furfural mixtures with different alcohols, and Zaoui-Djelloul-Daouadji et al. [14] reported the same properties for furfuryl alcohol and its mixtures with toluene or ethanol, at temperatures (283.15-313.15)K and 0.1 MPa.Almeida et al. [15] presented density and refractive index data of furfural at (288.15-318.15)K and atmospheric pressure (≈ 95 kPa), whereas Hough et al. [16] reported density of furfuryl alcohol at temperatures (303.15-343.15)K within a study of the heat capacities of organic liquids at atmospheric pressure.Further, Qureshi et al. [17] presented data on density and viscosity of furfuryl alcohol in the temperature range (288.15-318.15)K at 0.1 MPa within their study on biofuel-relevant compounds.Nduli and Deenadayalu [18] reported density, speed of sound and refractive index of furfural and furfuryl alcohol at temperatures (298.15-318.15)K, while Mahi et al. [19] gave the same thermodynamic properties but only for furfuryl alcohol in the temperature interval (293.15-323.15)K, all at atmospheric pressure.More recently, Belhadj et al. [20] published the density, speed of sound and refractive index of furfural at temperatures (293.15-323.15)K and pressure of 0.1 MPa.
In addition to the thermodynamic properties at atmospheric pressure, it is necessary to know the thermodynamic behaviour at high pressures and temperatures, because many processes as well as HDO [21] take place under such conditions.Only two papers reporting thermodynamic properties of furfural or furfuryl alcohol at higher temperatures and pressures were found during a literature review.Guerrero et al. [22] studied density and the derived thermodynamic properties of furfural and furfuryl alcohol at pressures up to 60 MPa, and at temperatures (283.15-338.15)K.In the second mentioned paper, Baird et al. [23] studied vapour pressures and densities of various bio-compounds and reported density of furfural in the wide range of temperature (293.15-448.15)K at pressures up to about 10 MPa.To the best of our knowledge, there are no published data on density, viscosity, speed of sound and refractive index of binary mixtures furfural + furfuryl alcohol.
In this paper, density, speed of sound and refractive index of furfural, furfuryl alcohol and their binary mixtures are reported at temperature range (288.15-343.15)K and 0.1 MPa while their viscosities are obtained under the same pressure over the temperature range (288.15-373.15)K. Further, the densities of the pure components are measured at higher pressures, up to 60 MPa, and in the temperature interval (293.15-413.15)K for furfural and (293.15-373.15)K for furfuryl alcohol.The high-pressure density data were fitted to the modified Tammann-Tait equation [24], which enabled the determination of the derived

Apparatus and Procedure
All experimental data in this paper were obtained using high precision instrumentation from the manufacturer Anton Paar (8054 Graz, Austria).Density, ρ, at atmospheric pressure and speed of sound, u, were determined by means of DSA 5000 M. The DSA 5000 M device (Anton Paar, 8054 Graz, Austria) enables the measurement of density in the range of (0-3) g•cm −3 and speed of sound within the interval (1000-2000) m•s −1 , both at temperatures (273.15-343.15)K.The device contains two measuring cells: one made of stainless steel used for measuring of speed of sound passing through the sample and Ushaped tube made of borosilicate glass where density is measured.Densities are calculated from the oscillation period of U tube filled with sample with automatic viscosity correction.The calibration of the device is performed daily using ambient air and ultrapure water.
Dynamic viscosity, η, at atmospheric pressure was measured using Stabinger viscometer SVM 3000/G2 (Anton Paar, 8054 Graz, Austria).Besides dynamic viscosity in the range (0.2-20,000) mPa•s; this device also measures kinematic viscosity in the range (0.2-20,000) mm 2 •s −1 and density within (0.65-3) g•cm −3 , at temperatures starting from 20 K below room temperature to 378.15 K.The SVM 3000/G2 is a rotational viscometer with a very small measuring cell containing a tube filled with sample that rotates at a constant speed and a measuring rotor, which floats in the sample.The viscosity measurement is based on the measurement of torque and speed of the rotor immersed in the sample.In addition to the cylinder for measuring the dynamic viscosity, the device also has an oscillating U-tube for density measurement.
Refractometer RXA-156 (Anton Paar, 8054 Graz, Austria) was used for the measurements of refractive index, nD, at atmospheric pressure.The refractive index range that can be measured on this device is 1.32-1.56at temperatures (283.15-343.15)K.The work of this refractometer is based on the measurement of the critical angle of the total reflection of light (wavelength of 589.3 nm) after passing through a sample.
The devices and the measurements procedures were described in more detail in our previous publications [25,26].The expanded uncertainties of the performed furfural and (293.15-373.15)K for furfuryl alcohol.The high-pressure density data were fitted to the modified Tammann-Tait equation [24], which enabled the determination of the derived thermodynamic properties, such as the isothermal compressibility, κT, isobaric thermal expansivity, αp, the internal pressure, pint, and the difference in isobaric and isochoric heat capacities, cp-cv.In addition, isobaric and isochoric specific heat capacities for furfural and furfuryl alcohol at 0.1 MPa were calculated.

Apparatus and Procedure
All experimental data in this paper were obtained using high precision instrumentation from the manufacturer Anton Paar (8054 Graz, Austria).Density, ρ, at atmospheric pressure and speed of sound, u, were determined by means of DSA 5000 M. The DSA 5000 M device (Anton Paar, 8054 Graz, Austria) enables the measurement of density in the range of (0-3) g•cm −3 and speed of sound within the interval (1000-2000) m•s −1 , both at temperatures (273.15-343.15)K.The device contains two measuring cells: one made of stainless steel used for measuring of speed of sound passing through the sample and Ushaped tube made of borosilicate glass where density is measured.Densities are calculated from the oscillation period of U tube filled with sample with automatic viscosity correction.The calibration of the device is performed daily using ambient air and ultrapure water.
Dynamic viscosity, η, at atmospheric pressure was measured using Stabinger viscometer SVM 3000/G2 (Anton Paar, 8054 Graz, Austria).Besides dynamic viscosity in the range (0.2-20,000) mPa•s; this device also measures kinematic viscosity in the range (0.2-20,000) mm 2 •s −1 and density within (0.65-3) g•cm −3 , at temperatures starting from 20 K below room temperature to 378.15 K.The SVM 3000/G2 is a rotational viscometer with a very small measuring cell containing a tube filled with sample that rotates at a constant speed and a measuring rotor, which floats in the sample.The viscosity measurement is based on the measurement of torque and speed of the rotor immersed in the sample.In addition to the cylinder for measuring the dynamic viscosity, the device also has an oscillating U-tube for density measurement.
Refractometer RXA-156 (Anton Paar, 8054 Graz, Austria) was used for the measurements of refractive index, nD, at atmospheric pressure.The refractive index range that can be measured on this device is 1.32-1.56at temperatures (283.15-343.15)K.The work of this refractometer is based on the measurement of the critical angle of the total reflection of light (wavelength of 589.3 nm) after passing through a sample.
The devices and the measurements procedures were described in more detail in our previous publications [25,26].The expanded uncertainties of the performed Aldrich Chemicals 0.98

Apparatus and Procedure
All experimental data in this paper were obtained using high precision instrumentation from the manufacturer Anton Paar (8054 Graz, Austria).Density, ρ, at atmospheric pressure and speed of sound, u, were determined by means of DSA 5000 M. The DSA 5000 M device (Anton Paar, 8054 Graz, Austria) enables the measurement of density in the range of (0-3) g•cm −3 and speed of sound within the interval (1000-2000) m•s −1 , both at temperatures (273.15-343.15)K.The device contains two measuring cells: one made of stainless steel used for measuring of speed of sound passing through the sample and Ushaped tube made of borosilicate glass where density is measured.Densities are calculated from the oscillation period of U tube filled with sample with automatic viscosity correction.The calibration of the device is performed daily using ambient air and ultra-pure water.
Dynamic viscosity, η, at atmospheric pressure was measured using Stabinger viscometer SVM 3000/G2 (Anton Paar, 8054 Graz, Austria).Besides dynamic viscosity in the range (0.2-20,000) mPa•s; this device also measures kinematic viscosity in the range (0.2-20,000) mm 2 •s −1 and density within (0.65-3) g•cm −3 , at temperatures starting from 20 K below room temperature to 378.15 K.The SVM 3000/G2 is a rotational viscometer with a very small measuring cell containing a tube filled with sample that rotates at a constant speed and a measuring rotor, which floats in the sample.The viscosity measurement is based on the measurement of torque and speed of the rotor immersed in the sample.In addition to the cylinder for measuring the dynamic viscosity, the device also has an oscillating U-tube for density measurement.
Refractometer RXA-156 (Anton Paar, 8054 Graz, Austria) was used for the measurements of refractive index, n D , at atmospheric pressure.The refractive index range that can be measured on this device is 1.32-1.56at temperatures (283.15-343.15)K.The work of this refractometer is based on the measurement of the critical angle of the total reflection of light (wavelength of 589.3 nm) after passing through a sample.
The devices and the measurements procedures were described in more detail in our previous publications [25,26].The expanded uncertainties of the performed measurements at 0.1 MPa were estimated based on the repeatability of the measurements, including the influence of samples' purities [27].The impurities in the used chemicals were not analysed but, based on the literature [5], it was assumed that for furfural (1%) they mostly are 2-methylfuran and 2-acetyl furan and for furfuryl alcohol (2%), the expected impurities are furfural, 2-methylfuran and resins.The density was determined with an expanded uncertainty, U, of 0.8 kg•m −3 with a 95 % level of confidence (coverage factor, k = 2) and the expanded uncertainties (confidence level 95 %, k = 2) in speed of sound, viscosity and refractive index were estimated to 0.2 m•s −1 , 0.007 mPa•s and 2.8•10 −3 , respectively.
All mixtures were prepared using a Mettler Toledo AG 204 mass balance (8606 Greifensee, Switzerland) with the precision 1•10 −7 kg.The estimated standard uncertainty, u, in mole fraction, x, was less than ±1•10 −4 .Density at higher pressures, p, up to 60 MPa and temperatures, T, in the range (293.15-413.15)K for furfural and (293.15-373.15)K for furfuryl alcohol, were experimentally determined on Anton Paar digital vibrating tube densimeter DMA HP (Anton Paar, 8054 Graz, Austria) by measuring the oscillation period of a U-shaped tube filled with the sample.The temperature of the measuring cell was precisely controlled by an integrated Peltier thermostat and the expanded uncertainty in temperature, U(T) with confidence level of 95% (k = 2) was estimated to 0.02 K.A Pressure Generator model 50-6-15 (High Pressure Equipment Co, Erie, PA 16505, USA) was used to achieve and control the desired pressure by means of acetone as the hydraulic fluid.A Wika S-10 pressure transducer (Alexander Wiegand GmbH & Co, 63911 Klingenberg, Germany) was used for pressure detection and an estimated expanded uncertainty in pressure, 95% confidence level (k = 2) was 0.1 MPa.The classic calibration procedure with one reference fluid, proposed by Comuñas et al. [28], was applied to obtain density values from the measured periods of oscillations of the U-tube.A detailed description of the apparatus, vibrating tube densimeter calibration and high-pressure density measurement procedure is described elsewhere [29].Since the measuring cell DMA HP does not provide automatic correction of density due to effect of sample's viscosity, it had to be estimated after experiment based on the information received from the supplier Anton Paar.For the studied compounds, furfural and furfuryl alcohol, the mentioned correction of the density was less than 0.08 kg•m −3 .The overall combined expanded uncertainty of the reported densities measured at DMA HP measuring cell was estimated taking into account the uncertainties of all properties that affect the density measurements, such as the oscillation period of U tube, temperature and pressure, and the densities of the reference fluids.Besides that, the damping effects on the vibrating tube and the impurities present in the used chemicals were considered, which resulted in the expanded uncertainty in density, U(ρ) of 0.9 kg•m −3 at temperature (293.15-363.15)K and 1.3 kg•m −3 at higher temperatures (373.15-413.15)K (confidence level 95%, k = 2).

Experimental Results
Thermodynamic properties such as density, speed of sound and refractive index of binary mixtures furfural + furfuryl alcohol, were measured at 0.1 MPa over the temperature range T = (293.15-343.15)K.The dynamic viscosity of the mentioned mixtures was studied at somewhat wider temperature range of (288.15-373.15)K and the same pressure, 0.1 MPa.The results of the measurements are listed in Table 2.
[14], Qureshi et al. [17] and Baird et al. [23] are of the same order of magnitude as the uncertainties estimated for the results given in this work.Mahi et al. [19] and Belhadj et al. [20] assessed higher uncertainties while the uncertainties of the other compared data [12,13,15,16,18,22] were lower than those reported here.The agreement of the results presented in this paper with the literature data [12][13][14][15][16][17][18][19][20] is very good and falls mostly within the reported uncertainties.The unsatisfactory results were obtained only when the density and speed of sound of furfuryl alcohol were compared with the data of Zaoui- the data given here and those of Lomba et al. [12] was about 1.11% (0.05 mPa•s).The measured speed of sound of furfuryl alcohol deviates from data of Lomba et al. [12] and Zaoui-Djelloul-Daouadji et al. [14] with AAD of 0.03% and 0.16% (0.4 m•s −1 and 2.3 m•s −1 ), respectively.The comparison of speed of sound data with the values reported by Mahi et al. [19] gave similar results, AAD ≈ 0.11% (0.6 m•s −1 ) but for the results of Nduli and Deenadayalu [18] AAD was about 0.44% (6.2 m•s −1 ).Finally, Lomba et al. [12], Nduli and Deenadayalu [18] and Mahi et al. [19] also reported the refractive indices of furfuryl alcohol, which agree well with the experimental values presented here; the AAD was about 0.09%, 0.02% and 0.04% (0.001, 0.0003 and 0.0006), respectively.The measurement uncertainties reported by Zaoui-Djelloul-Daouadji et al.
The measurement uncertainties reported by Zaoui-Djelloul-Daouadji et al. [14], Qureshi et al. [17] and Baird et al. [23] are of the same order of magnitude as the uncertainties estimated for the results given in this work.Mahi et al. [19] and Belhadj et al. [20] assessed higher uncertainties while the uncertainties of the other compared data [12,13,15,16,18,22] were lower than those reported here.The agreement of the results presented in this paper with the literature data [12][13][14][15][16][17][18][19][20] is very good and falls mostly within the reported uncertainties.The unsatisfactory results were obtained only when the density and speed of sound of furfuryl alcohol were compared with the data of Zaoui-Djelloul-Daouadji et al. [14] and Nduli and Deenadayalu [18].The comparison of experimental results with literature data performed by the mentioned authors in their papers [14,18] gave deviations similar to those obtained comparing the data presented here with their results [14,18].Therefore, the poor agreement of the measured density and speed of sound of furfuryl alcohol with the data published in two mentioned papers [14,18] should not cast doubt on the reliability of the presented results.Density, speed of sound and refractive index decrease linearly as a function of temperature (Figure 2).Density and refractive index increase for increasing concentration of furfural in the mixture because these thermodynamic properties of pure furfural are higher than of furfuryl alcohol.The dependence of the mentioned properties on the furfural fraction in the mixtures is given at Figure S1.The density, speed of sound and refractive index measured for the studied mixtures are higher than the values that would be expected for ideal mixture (based on Kay's rule).Deviation of the mixture's properties from ideal behaviour indicates the specific interactions between components and in this case the increase in ρ, u and n D implies the presence of attractive forces, assumingly hydrogen bonds, between furfural and furfuryl alcohol.The influence of the composition of the mixture on speed of sound is weaker than on the other studied properties.Furthermore, the addition of small amount of furfural (up to 40%) caused the increase in speed of sound while further increase in furfural fraction led to decrease in speed of sound of binary mixture.
Energies 2021, 14, x. https://doi.org/10.3390/xxxxxwww.mdpi.com/journal/energiesbe expected for ideal mixture (based on Kay's rule).Deviation of the mixture's properties from ideal behaviour indicates the specific interactions between components and in this case the increase in ρ, u and nD implies the presence of attractive forces, assumingly hydrogen bonds, between furfural and furfuryl alcohol.The influence of the composition of the mixture on speed of sound is weaker than on the other studied properties.Furthermore, the addition of small amount of furfural (up to 40%) caused the increase in speed of sound while further increase in furfural fraction led to decrease in speed of sound of binary mixture.
As for the dependence of viscosity on temperature, it decreases exponentially with temperature elevation, as expected (Figure 2b).The obtained results showed that the increase in the molar fraction of furfural in mixture led to lower viscosity due to lower viscosity of furfural than of furfuryl alcohol.Additionally, the measured viscosities of the studied mixtures are lower than the values expected for ideal mixture (Grunberg-Nissan rule) indicating the heteromolecular interactions (Figure S1b).Furfural and furfuryl alcohol are both derivatives of furan, where hydrogen at position 2 is substituted with formyl or hydroxymethyl group, respectively (see Table 1).Both compounds have similar densities, higher than water.However, furfuryl alcohol is more viscous comparing to furfural presumably due to hydrogen bonding, especially at lower temperatures.The difference in viscosities diminishes at higher temperatures which might be attributed to the rupture of hydrogen bonding in furfuryl alcohol.The experimentally obtained viscosity data were fitted to the Vogel-Fulcher-Tammann (VFT) Equation (1) [30][31][32] leading to the optimized parameters given in Table 3.
[14], Qureshi et al. [17] and Baird et al. [23] are of the same order of magnitude as the uncertainties estimated for the results given in this work.Mahi et al. [19] and Belhadj et al. [20] assessed higher uncertainties while the uncertainties of the other compared data [12,13,15,16,18,22] were lower than those reported here.The agreement of the results presented in this paper with the literature data [12][13][14][15][16][17][18][19][20] is very good and falls mostly within As for the dependence of viscosity on temperature, it decreases exponentially with temperature elevation, as expected (Figure 2b).The obtained results showed that the increase in the molar fraction of furfural in mixture led to lower viscosity due to lower viscosity of furfural than of furfuryl alcohol.Additionally, the measured viscosities of the studied mixtures are lower than the values expected for ideal mixture (Grunberg-Nissan rule) indicating the heteromolecular interactions (Figure S1b).Furfural and furfuryl alcohol are both derivatives of furan, where hydrogen at position 2 is substituted with formyl or hydroxymethyl group, respectively (see Table 1).Both compounds have similar densities, higher than water.However, furfuryl alcohol is more viscous comparing to furfural presumably due to hydrogen bonding, especially at lower temperatures.The difference in viscosities diminishes at higher temperatures which might be attributed to the rupture of hydrogen bonding in furfuryl alcohol.
The experimentally obtained viscosity data were fitted to the Vogel-Fulcher-Tammann (VFT) Equation ( 1) [30][31][32] leading to the optimized parameters given in Table 3.The obtained good agreement between calculated and measured viscosities (AAD ≈ 0.2%) confirmed the suitability of VFT equation for viscosity correlation (Figure 2b).
Density measurements for pure compounds were also carried out at pressures up to 60 MPa and at temperatures ranging (293.15-413.15)K for furfural and (293.15-373.15)K for furfuryl alcohol (Table 4).Densities of the studied compounds at 0.1 MPa pressure determined using DMA HP device were compared with the values measured by means of DSA 5000 M and the average absolute deviations were 0.01% for furfural and 0.04% for furfuryl alcohol.
Guerrero et al. [22] measured densities of furfural and furfuryl alcohol at temperatures (283.15-338.15)K and pressures up to 60 MPa and their results differ from the data presented in this paper for about 0.02% (0.3 kg•m −3 ) and 0.03% (0.4 kg•m −3 ), respectively.Baird et al. [23] determined the density and vapour pressure of 11 biocompounds in the temperature interval (293.15-423.15)K at pressures up to 10 MPa and one of them was furfural; the agreement between their and the data given here was very good with AAD ≈ 0.03% (0.4 kg•m −3 ).
The experimentally determined densities of the pure compounds as a function of temperature and pressure are presented at Figure 3 showing that the density depends on temperature almost linearly.As expected, an increase in density with pressure rise and decrease in density as a function of increasing temperature was noted for both of the studied compounds.Densities of the studied compounds at 0.1 MPa pressure determined using DMA HP device were compared with the values measured by means of DSA 5000 M and the average absolute deviations were 0.01% for furfural and 0.04% for furfuryl alcohol.Guerrero et al. [22] measured densities of furfural and furfuryl alcohol at temperatures (283.15-338.15)K and pressures up to 60 MPa and their results differ from the data presented in this paper for about 0.02% (0.3 kg•m −3 ) and 0.03% (0.4 kg•m −3 ), respectively.Baird et al. [23] determined the density and vapour pressure of 11 biocompounds in the temperature interval (293.15-423.15)K at pressures up to 10 MPa and one of them was furfural; the agreement between their and the data given here was very good with AAD ≈ 0.03% (0.4 kg•m −3 ).
The experimentally determined densities of the pure compounds as a function of temperature and pressure are presented at Figure 3 showing that the density depends on temperature almost linearly.As expected, an increase in density with pressure rise and decrease in density as a function of increasing temperature was noted for both of the studied compounds.

High Pressure Density Correlation
The correlation of the experimentally determined high-pressure densities was performed applying the modified Tammann-Tait equation [24] (Equation ( 2)).That further enabled the calculation of various derived properties.

High Pressure Density Correlation
The correlation of the experimentally determined high-pressure densities was performed applying the modified Tammann-Tait equation [24] (Equation ( 2)).That further enabled the calculation of various derived properties.
where ρ ref represents the density of sample at the reference pressure, p ref = 0.1 MPa, B(T) is temperature dependent, while C is temperature independent parameter.ρ ref and B(T) are expressed in a form of second-order polynomial: where a i and b i , in addition to C, are adjustable parameters.Firstly, for each individual compound, the density data obtained at reference pressure (0.1 MPa) were fitted to Equation (3) which resulted in the determination of parameters a i .The parameter optimisation was conducted by applying the Levenberg-Marquardt algorithm [33] aiming to minimize the absolute average deviation between the measured and correlated values.The second step was to adjust parameters b i and C of Equations ( 2) and (4) by applying the same optimization procedure to the whole density data set, excluding densities at 0.1 MPa.The optimized parameters of the modified Tammann-Tait equation (Equations ( 2)-( 4)) are given in Table 5.The average absolute percentage deviation (AAD), the maximum percentage deviation (MD), the average percentage deviation (Bias) and the standard deviation (σ) of the experimental data from those calculated using the modified Tammann-Tait equation are also given in Table 5.The low values of AAD (about 0.006%) achieved for both studied pure compounds indicate the good quality of density data modelling.

Derived Thermodynamic Properties
The knowledge of the density at wide ranges of temperature and pressure enables calculation of various derived volumetric properties.These properties are calculated by differentiating density with respect to pressure or temperature.
The change of density as a response to pressure change is described by the isothermal compressibility, κ T : [34] The incorporation of Equation ( 2) into (5) gives: [34] The influence of temperature on density, i.e., the change of density when temperature is changed under constant pressure, is described by the isobaric thermal expansivity, α p : [34] Energies 2021, 14, 7769 13 of 18 The following expression is derived from the modified Tammann-Tait Equations ( 2) and ( 7): [34] where ρ ref '(T, p ref ) and B'(T) are derivatives of the parameters ρ ref (T, p ref ) and B(T) with respect to T, respectively: The thermal pressure coefficient, γ, which represents the ratio between κ T and α p , can be calculated as follows: [34] The internal pressure, p int , which gives insight into intermolecular interaction of the sample can be determined using Equation ( 12): [34] where U stands for an internal energy and V is volume of the sample.Another important thermodynamic property is the difference between the specific heat capacity at constant pressure, c p , and the specific heat capacity at constant volume, c v : [34] The coupling of Equations ( 5) and ( 7) with (13) leads to the following expression of the mentioned property: Knowledge on the isothermal and isentropic compressibility enables the calculation of the isobaric specific heat capacity: [35] which is significant for the determination of the isochoric heat capacity using Equation ( 14).The calculation of isentropic compressibility, κ S , requires the knowledge on density and speed of sound: [35] The calculated isothermal compressibility, the isobaric thermal expansivity, the internal pressure and the difference between the isobaric and isochoric specific heat capacities for both examined pure compounds, in the temperature interval (293.15-413.15)K for furfural and (293.15-437.15)K for furfuryl alcohol at pressures up to 60 MPa, are given in the Supplementary Material to the paper (Tables S1 and S2).The isothermal compressibility and the isobaric thermal expansivity calculated for furfural and furfuryl alcohol are presented in Figures 4 and 5.
The calculated isothermal compressibility, the isobaric thermal expansivity, the internal pressure and the difference between the isobaric and isochoric specific heat capacities for both examined pure compounds, in the temperature interval (293.15-413.15)K for furfural and (293.15-437.15)K for furfuryl alcohol at pressures up to 60 MPa, are given in the Supplementary Material to the paper (Tables S1 and S2).The isothermal compressibility and the isobaric thermal expansivity calculated for furfural and furfuryl alcohol are presented in Figures 4 and 5.The measurement uncertainties reported by Zaoui-Djelloul-Daouadji et al.
[14], Qureshi et al. [17] and Baird et al. [23] are of the same order of magnitude as the uncertainties estimated for the results given in this work.Mahi et al. [19] and Belhadj et al. [20] assessed higher uncertainties while the uncertainties of the other compared data [12,13,15,16,18,22] were lower than those reported here.The agreement of the results presented in this paper with the literature data [12][13][14][15][16][17][18][19][20] is very good and falls mostly within The calculated isothermal compressibility, the isobaric thermal expansivity, the internal pressure and the difference between the isobaric and isochoric specific heat capacities for both examined pure compounds, in the temperature interval (293.15-413.15)K for furfural and (293.15-437.15)K for furfuryl alcohol at pressures up to 60 MPa, are given in the Supplementary Material to the paper (Tables S1 and S2).The isothermal compressibility and the isobaric thermal expansivity calculated for furfural and furfuryl alcohol are presented in Figures 4 and 5.The measurement uncertainties reported by Zaoui-Djelloul-Daouadji et al.
[14], Qureshi et al. [17] and Baird et al. [23] are of the same order of magnitude as the uncertainties estimated for the results given in this work.Mahi et al. [19] and Belhadj et al. [20] assessed higher uncertainties while the uncertainties of the other compared data [12,13,15,16,18,22] were lower than those reported here.The agreement of the results presented in this paper with the literature data [12][13][14][15][16][17][18][19][20] is very good and falls mostly within The measurement uncertainties repo Qureshi et al. [17] and Baird et al. [23] are tainties estimated for the results given in th assessed higher uncertainties while the [12,13,15,16,18,22] were lower than those re sented in this paper with the literature dat The isothermal compressibility and the isobaric thermal expansivity of both studied pure compounds increase as temperature rises at constant pressure and decrease with increasing pressure along the isotherms (Figures 4 and 5).The obtained values of κT are slightly higher for furfural than for furfuryl alcohol so the change of pressure will affect densities of both compounds similarly.The isothermal compressibility is inversely proportional to the ability of molecules to create hydrogen bonds [22], which could explain the lower values of κT obtained for furfuryl alcohol.As for αp, the calculated values are somewhat higher for furfural than for furfuryl alcohol, meaning that the increase in temperature will cause greater expansion, i.e., a density decrease in the case of furfural than of furfuryl alcohol.This indicates that furfuryl alcohol has a better packed structure than furfural, likely because of stronger intermolecular interactions due to hydrogen bonds that are also known to limit the movement of molecules and that way disable the expansion [12,36].The isobaric thermal expansivity shows the typical behaviour where its dependence on temperature becomes weaker with the increase in pressure resulting in the intersection of the isotherms.The intersection point where αp is temperature independent the data given here and those of Lomba et al. [12] was about 1.11% (0.05 mPa•s).The measured speed of sound of furfuryl alcohol deviates from data of Lomba et al. [12] and Zaoui-Djelloul-Daouadji et al. [14] with AAD of 0.03% and 0.16% (0.4 m•s −1 and 2.3 m•s −1 ), respectively.The comparison of speed of sound data with the values reported by Mahi et al. [19] gave similar results, AAD ≈ 0.11% (0.6 m•s −1 ) but for the results of Nduli and Deenadayalu [18] AAD was about 0.44% (6.2 m•s −1 ).Finally, Lomba et al. [12], Nduli and Deenadayalu [18] and Mahi et al. [19] also reported the refractive indices of furfuryl alcohol, which agree well with the experimental values presented here; the AAD was about 0.09%, 0.02% and 0.04% (0.001, 0.0003 and 0.0006), respectively.The measurement uncertainties reported by Zaoui-Djelloul-Daouadji et al.
The measurement uncertainties reported by Zaoui-Djelloul-Daouadji et al.
[14], Qureshi et al. [17] and Baird et al. [23] are of the same order of magnitude as the uncertainties estimated for the results given in this work.Mahi et al. [19] and Belhadj et al. [20] assessed higher uncertainties while the uncertainties of the other compared data [12,13,15,16,18,22] were lower than those reported here.The agreement of the results presented in this paper with the literature data [12][13][14][15][16][17][18][19][20] The coupling of Equations ( 5) and ( 7) with ( 13) leads to the following expression of the mentioned property: Knowledge on the isothermal and isentropic compressibility enables the calculation of the isobaric specific heat capacity: [35] 2 () Which is significant for the determination of the isochoric heat capacity using Equation ( 14).The calculation of isentropic compressibility, κS, requires the knowledge on density and speed of sound: [35] The calculated isothermal compressibility, the isobaric thermal expansivity, the internal pressure and the difference between the isobaric and isochoric specific heat capacities for both examined pure compounds, in the temperature interval (293.15-413.15)K for furfural and (293.15-437.15)K for furfuryl alcohol at pressures up to 60 MPa, are given in the Supplementary Material to the paper (Tables S1 and S2).The isothermal compressibility and the isobaric thermal expansivity calculated for furfural and furfuryl alcohol are presented in Figures 4 and 5   ) but for the results of Nduli and Deenadayalu [18] AAD was about 0.44% (6.2 m•s −1 ).Finally, Lomba et al. [12], Nduli and Deenadayalu [18] and Mahi et al. [19] also reported the refractive indices of furfuryl alcohol, which agree well with the experimental values presented here; the AAD was about 0.09%, 0.02% and 0.04% (0.001, 0.0003 and 0.0006), respectively.The measurement uncertainties reported by Zaoui-Djelloul-Daouadji et al.
[14], Qureshi et al. [17] and Baird et al. [23] are of the same order of magnitude as the uncertainties estimated for the results given in this work.Mahi et al. [19] and Belhadj et al. [20] assessed higher uncertainties while the uncertainties of the other compared data [12,13,15,16,18,22] were lower than those reported here.The agreement of the results presented in this paper with the literature data [12][13][14][15][16][17][18][19][20]  ).Fina [18] and Mahi et al. [19] also reported the agree well with the experimental values pr and 0.04% (0.001, 0.0003 and 0.0006), respe The measurement uncertainties repo Qureshi et al. [17] and Baird et al. [23] are tainties estimated for the results given in th assessed higher uncertainties while the [12,13,15,16,18,22] were lower than those re sented in this paper with the literature dat The isothermal compressibility and the isobaric thermal expansivity of both studied pure compounds increase as temperature rises at constant pressure and decrease with increasing pressure along the isotherms (Figures 4 and 5).The obtained values of κ T are slightly higher for furfural than for furfuryl alcohol so the change of pressure will affect densities of both compounds similarly.The isothermal compressibility is inversely proportional to the ability of molecules to create hydrogen bonds [22], which could explain the lower values of κ T obtained for furfuryl alcohol.As for α p , the calculated values are somewhat higher for furfural than for furfuryl alcohol, meaning that the increase in temperature will cause greater expansion, i.e., a density decrease in the case of furfural than of furfuryl alcohol.This indicates that furfuryl alcohol has a better packed structure than furfural, likely because of stronger intermolecular interactions due to hydrogen bonds that are also known to limit the movement of molecules and that way disable the expansion [12,36].The isobaric thermal expansivity shows the typical behaviour where its dependence on temperature becomes weaker with the increase in pressure resulting in the intersection of the isotherms.The intersection point where α p is temperature independent ∂α p ∂T p = 0 was not observed at the studied range of pressure for furfural, while for furfuryl alcohol, it is expected to occur at pressure slightly above 60 MPa.In addition to this, the isobaric thermal expansivities of both studied compounds were calculated using a pseudo-experimental technique where the measured densities at constant pressure were fitted using polynomial function: The values obtained this way agreed very well with the α p calculated by differentiating the modified Tammann-Tait equation (Tables S1 and S2); the average absolute percentage deviations for α p of furfural was less than 0.2% and for furfuryl alcohol it was about 0.5%.
The dependence of the internal pressure on pressure is given in Figure S2.The p int represents the change in internal energy as a result of a very small change in volume at constant temperature and gives insight mainly into weak intermolecular forces such as dispersive, repulsive and dipolar [36].The internal pressures of the studied compounds are positive implicating the dominant attractive intermolecular forces.The internal pressure decreases with the increase in temperature while pressure does not affect it considerably, especially at lower temperatures, for both studied compounds.The increase in pressure leads to lower values of p int of furfural at temperatures under 343.15K while at temperatures above 343.15K p int increases when pressure rises.In the case of furfuryl alcohol, the internal pressure decreases when pressure increases in the whole studied range of temperature indicating that the higher pressure restricts the change of internal energy as a respond to expansion.The p int values calculated for furfural are higher than those obtained for furfuryl alcohol, which is more noticeable at lower temperatures, possibly due to stronger dipole-dipole interactions within more polar furfural when compared with furfuryl alcohol [12].Although molecules of furfuryl alcohol are linked by hydrogen bonds that are stronger than dipole-dipole interactions between molecules of furfural, hydrogen bonds do not have significant effect on the internal pressure [36].
The isentropic compressibility, κ S , of furfural and furfuryl alcohol (Table 6) increases with the increase in temperature.A comparison of isothermal (Tables S1 and S2) and isentropic compressibility (Table 6) showed that the values determined for κ T are for about 25% higher than those calculated for κ S, for both studied compounds.In general, the ratio κ T /κs is equal to the ratio between isobaric and isochoric heat capacities, C p /C V .While isentropic compressibility of furfural at 0.1 MPa is higher than of furfuryl alcohol (for about 6%), the isentropic compressibilities of the examined samples are almost the same.The κ S of furfuryl alcohol is slightly higher than κ S of furfural at temperatures up to 313.15 K and at higher temperatures, the relation is reversed.
The calculated values of the specific heat capacities at constant pressure and at constant volume are given in Table 6.The obtained values for both, c p and c V , are higher for furfuryl alcohol than furfural.The increase in temperature resulted in the increase in the isobaric and isochoric specific heat capacities and the temperature influence was stronger in the case of furfuryl alcohol.Table 6.The isentropic compressibility, κ S , the isobaric specific heat capacity, c p , and the isochoric specific heat capacity, c V , of furfural and furfuryl alcohol at 0.1 MPa pressure.

Conclusions
In this work, several thermodynamic and transport properties such as density, viscosity, speed of sound, and refractive index of the binary mixture furfural + furfuryl alcohol were reported at various temperatures and 0.1 MPa.In addition, densities of pure furfural and furfuryl alcohol were determined along with the derived thermodynamic properties at pressures up to 60 MPa over the temperature range from (293.15-413.15)K for furfural and (293.15-373.15)K for furfuryl alcohol.
Experimental results at 0.1 MPa show that furfural has higher density and refractive index than furfuryl alcohol, while viscosities are higher for furfuryl alcohol than furfural.Speeds of sound have more or less similar values, which are more noticeable at low temperatures.All of the studied thermodynamic properties decrease as temperature rises, where the density, speed of sound and refractive index decrease linearly while viscosity depends exponentially on temperature.
Concerning (ρ, p, T) data, as expected, the densities of the examined substances decrease with the increase in temperature and rise when pressure increases.The correlation of the measured high-pressure density data was successfully performed applying the modified Tammann-Tait equation.That further led to the determination of the isothermal compressibility, the isobaric thermal expansivity, the internal pressure and the difference between the isobaric and isochoric specific heat capacities for both examined pure compounds.The isothermal compressibility and the isobaric thermal expansivity are higher for furfural than for furfuryl alcohol and they increase as temperature rises or pressure decreases.The higher values of κ T and α p indicate that the change of pressure and temperature will affect density of furfural more than the density of furfuryl alcohol.Furfural also has higher internal pressure in comparison to furfuryl alcohol, which could be the result of the stronger dipole-dipole intermolecular forces within furfural since hydrogen bonds (dominant between molecules of furfuryl alcohol) do not affect the internal pressure.S1: The isothermal compressibility, κ T , the isobaric thermal expansivity, α p , the difference between isobaric and isochoric specific thermal capacity, c p -c v , and internal pressure, p int , of furfural, Table S2: The isothermal compressibility, κ T , the isobaric thermal expansivity, α p , the difference between isobaric and isochoric specific thermal capacity, c p -c v , and internal pressure, p int , of furfuryl alcohol.

Figure 1 .
Figure 1.The percentage deviations of the measured (a) density, ρ, (b) dynamic viscosity, η, (c) speed of sound, u, and (d) refractive index, n D , from the literature data reported by: (

Table 1 .
The chemicals used in the measurements.

Table 1 .
The chemicals used in the measurements.

Table 4 .
Density, ρ, of the pure furfural and furfuryl alcohol at temperatures (293.15-413.15)K and pressures up to 60 MPa.

Table 5 .
Parameters of the modified Tammann-Tait equation and comparison criteria for furfural and furfuryl alcohol at temperatures (293.15-413.15)K and at pressures up to 60 MPa.