Available online at www.sciencedirect.com
Fuel 87 (2008) 1941–1948 www.fuelfirst.com
The effect of temperature and pressure on the physicochemical properties of petroleum diesel oil and biodiesel fuel Marzena Dzida a
a,*
, Piotr Prusakiewicz
b
University of Silesia, Institute of Chemistry, Szkolna 9, 40-006 Katowice, Poland b Rafineria Trzebinia S.A., Fabryczna 22, 32-540 Trzebinia, Poland
Received 10 July 2007; received in revised form 13 October 2007; accepted 16 October 2007 Available online 3 December 2007
Abstract Three commercial fuels were studied: biodiesel (based mainly of the fatty acids methyl esters of rapeseed oil), diesel oil Ekodiesel Ultra (standard petroleum diesel oil with sulphur content less than 10 mg/kg), and ON BIO 10 (blend of 20 vol.% of biodiesel with 80 vol.% of standard petroleum diesel oil with sulphur content less than 10 mg/kg). The speeds of sound were measured within the temperatures from 293 to 318 K and at pressures from 0.1 to 101 MPa. The densities and heat capacities were measured under atmospheric pressure in the temperature range from 273 to 363 K and 283 to 359 K, respectively. Using the experimental results, the physicochemical properties such as: density, isentropic bulk modulus, heat capacity, and isobaric thermal expansion were calculated in the same temperature and pressure range as the speed of sound was measured. The results obtained show that although the bulk modulus of ON BIO 10 is higher than that of diesel oil Ekodiesel Ultra over the whole pressure range, the difference is rather small and can be compensated by temperature. Isobaric thermal expansivity of biodiesel decreases with pressure slightly less than that of the diesel oil Ekodiesel Ultra. It is approximately independent of temperature and composition of the fuel at pressures 40 ± 5 MPa. Ó 2007 Elsevier Ltd. All rights reserved. Keywords: Physicochemical properties; Diesel; Biodiesel
1. Introduction In the petrochemical and automotive industries, an important point of interest is the influence of ‘‘bio-component’’ additives, such as alcohols or esters, on the properties of fuels and greases. Biodiesel is an alternative diesel fuel consisting of alkyl monoesters of fatty acids obtained from vegetable oils or animal fats. Mixtures of biodiesel and petroleum-based diesel fuel can be used in diesel engines. Physicochemical properties of bio- and petroleum-based diesels, such as density, bulk modulus, and heat capacity, are different due to dissimilar chemical composition of the fuels. That may influence combustion and
*
Corresponding author. Tel.: +48 032 359 1612; fax: +48 032 259 9978. E-mail address:
[email protected] (M. Dzida).
0016-2361/$ - see front matter Ó 2007 Elsevier Ltd. All rights reserved. doi:10.1016/j.fuel.2007.10.010
exhaust emission. Density is one of the most important properties of fuels, because injection systems, pumps and injectors must deliver the amount of fuel precisely adjusted to provide proper combustion. Isobaric thermal expansion and bulk modulus, or its reciprocal, i.e. compressibility, characterize how temperature and pressure affect density. The bulk modulus of biodiesels is higher than that of petroleum diesels. Since the modulus determines the spray characteristics upon injection, it has been closely related to NOx content in the exhausts [1–5]. That content depends also on the temperature of exhausts, thereby on the heat capacity of fuel [6]. As the fuel injection in an engine is approximately an adiabatic process, the adiabatic bulk modulus seems to be more useful than the isothermal one in estimation of the fuel injection timing. That is particularly important for modern common rail systems where the pressure can reach 250 MPa
1942
M. Dzida, P. Prusakiewicz / Fuel 87 (2008) 1941–1948
instead of ca. 35 MPa. The only experimental method that leads directly to adiabatic modulus is the acoustic one, based on the measurement of the speed of sound. The method has found wide acceptance as a relatively simple tool for determination of thermodynamic properties, especially at high-pressures. Apart from the adiabatic modulus and compressibility, the isothermal compressibility, density, isobaric thermal expansion, and heat capacities of compressed liquids obtained by that method have been reported [7–12]. Sun et al. [7] even claimed that at elevated pressures, density determined by the acoustic method is more nearly accurate than that measured in a direct way. In the present paper, we report experimental speeds of sound in a standard petroleum diesel oil, as well as in a biodiesel and in ON BIO 10, being a blend of 20 vol.% of biodiesel with 80 vol.% of standard petroleum diesel oil. The data cover the pressure range from 0.1 to 101 MPa and the temperature interval from 293 to 318 K. The speeds are completed with densities and heat capacities measured at atmospheric pressure and temperatures 273–363 K and 283–359 K, respectively. The densities and isobaric heat capacities at pressures up to 100 MPa were calculated from experimental data according to the suggestion of Davis and Gordon [8]. To this end, a slightly modified procedure of Sun et al. was applied [7,9]. From the relationships between pressure, temperature and density, the material constants were determined: adiabatic bulk moduli and isobaric thermal expansions of the three fuels studied. We aimed this work mainly at comparison of effects of pressure and temperature on density, heat capacity, adiabatic bulk modulus, and isobaric thermal expansion of three commercial fuels. Furthermore, we tried to find out whether the differences in the physicochemical properties
of the studied fuels are significant from the practical point of view. Data similar to those reported in this work have been used in calculations of other thermodynamic parameters, in simulations, and as a reference material for tests of performance characteristics of diesel engines [5,13,14]. 2. Experimental section 2.1. Samples Three commercial fuels were studied: (i) Ekodiesel Ultra, a standard petroleum diesel oil with sulphur content less than 10 mg/kg, which fulfilled norm EN 590, manufactured by PKN Orlen S.A., (ii) Biodiesel, based mainly of fatty acids methyl esters of rapeseed oil, which fulfilled norm EN 14214, manufactured by Rafineria Trzebinia S.A., and (iii) ON BIO 10, blend of 20 vol.% of biodiesel with 80 vol.% of standard petroleum diesel oil with sulphur content less than 10 mg/kg, also made by Rafineria Trzebinia S.A. Properties of diesel and biodiesel fuels are given in Table 1. 2.2. Ultrasonic speed measurements The speed of sound in liquids under test has been measured at atmospheric and higher pressures using two apparatus designed and constructed in our laboratory [15,16]. Two measuring vessels of the same acoustic path and construction have been used, one of them designed for measurements under atmospheric pressure, the other one for measurements under elevated pressures. A single transmitting–receiving ceramic transducer operating at 2 MHz and an acoustic mirror have been applied. The measuring sets
Table 1 Physical and chemical properties of diesel oil Ekodiesel Ultra, ON BIO10, and biodiesel Property Cetane number Cloud point Flash point Kinematic viscosity at 40 °C Water content Sulfur content Carbon residue on 10% distillation residue FAME content Oxidation stability at 110 °C Acid value Iodine value Linoleic acid methyl ester Methanol content Distillation at 101.3 kPa: 95% recovered up to Up to 250 °C recovered Up to 350 °C recovered
Units
Ekodiesel Ultra from PKN Orlen SA
ON BIO10 from RT SA
European standard for diesel EN 590
Biodiesel from RT SA
European standard for biodiesel-FAME EN 14214
°C °C mm2/s
54.1 14 65 2.7
51 10 68 2.88
>51 <0 >56 2.0–4.5
51.7 – >130 4.5
>51 – 120 3.5–5.0
mg/kg mg/kg % m/m
64 8.8 0.02
110 9.9 0.03
<200 <10 <0.3
88 9.7 0.04
<500 <10 <0.3
% m/m
–
17.7
5.0
97.5
>96.5
h mg of KOH/g g of I2/ 100g % m/m % m/m
– –
– –
<25 –
8.1 0.33
>6.0 <0.50
–
–
–
115
<120
– –
– –
– –
8.7 <0.05
<12 <0.20
°C vol.% vol.%
330 40 98
339 30 98
<360 <65 >85
– – –
– – –
M. Dzida, P. Prusakiewicz / Fuel 87 (2008) 1941–1948
operate on the principle of the pulse-echo-overlap method. The pressure was provided by a hand operated hydraulic press and was measured with a strain gauge measuring system (Hottinger Baldwin System P3MD) with accuracy better than 0.15%. The temperature was measured using an Ertco Hart 850 platinum resistance thermometer (traceable to a NIST standard) with an uncertainty of ±0.05 K and resolution of 0.001 K. The uncertainty of the speed of sound measurements was estimated to be 0.03% at atmospheric pressure, 0.04% under pressures up to 60 MPa and 0.05% under pressures from 60 MPa to 101 MPa.
1943
Other details of the high-pressure device and the method of the speed of sound measurements can be found in the previous papers [15,16]. 2.3. Density measurements The densities at atmospheric pressure were measured using a vibrating tube densimeter Anton Paar DMA 5000. The uncertainty of the density measurements was 0.05 kg m3, whereas the repeatability was estimated to be better than 0.005 kg m3.
Table 2 Speed of sound in diesel oil Ekodiesel Ultra, ON BIO10, and biodiesel measured at pressures up to 101 MPa within the temperature range 293–318 K Ekodiesel Ultra
ON BIO10
Biodiesel
T (K)
p (MPa)
u (m s1)
T (K)
p (MPa)
u (m s1)
T (K)
p (MPa)
u (m s1)
292.88 292.87 292.86 292.86 292.87 292.86 292.86 292.86 298.16 298.15 298.15 298.15 298.15 298.16 298.15 298.15 303.14 303.14 303.14 303.13 303.14 303.13 303.14 303.14 308.10 308.13 308.13 308.13 308.13 308.13 308.13 308.13 313.11 313.11 313.11 313.11 313.11 313.10 313.10 313.10 318.28 318.30 318.29 318.29 318.29 318.29 318.29 318.29
0.10 15.20 30.40 45.60 60.79 75.99 91.18 101.32 0.10 15.20 30.39 45.59 60.79 75.99 91.18 101.32 0.10 15.20 30.39 45.59 60.79 75.99 91.18 101.32 0.10 15.20 30.39 45.59 60.79 75.99 91.18 101.32 0.10 15.20 30.39 45.59 60.79 75.99 91.18 101.32 0.10 15.20 30.39 45.59 60.79 75.99 91.18 101.32
1370.76 1446.68 1513.41 1574.04 1630.60 1682.79 1731.94 1763.07 1350.67 1427.29 1494.57 1555.73 1613.02 1665.71 1715.26 1746.72 1331.91 1408.42 1477.07 1539.25 1597.12 1650.84 1700.62 1732.67 1313.30 1390.52 1460.51 1523.58 1581.90 1636.44 1686.91 1718.90 1294.68 1373.44 1444.26 1508.45 1567.18 1622.22 1673.14 1705.42 1275.61 1356.21 1428.16 1493.14 1552.45 1608.16 1659.76 1692.19
292.91 292.86 292.86 292.86 292.86 292.86 292.86 292.86 298.17 298.15 298.15 298.15 298.15 298.15 298.15 298.15 303.16 303.14 303.14 303.14 303.13 303.13 303.13 303.13 308.13 308.13 308.13 308.13 308.13 308.13 308.13 308.12 313.12 313.10 313.10 313.11 313.10 313.10 313.10 313.10 318.30 318.30 318.30 318.30 318.30 318.30 318.30 318.29
0.10 15.20 30.39 45.59 60.79 75.99 91.18 101.32 0.10 15.20 30.39 45.59 60.79 75.99 91.18 101.32 0.10 15.20 30.39 45.59 60.79 75.99 91.18 101.32 0.10 15.20 30.39 45.59 60.79 75.99 91.18 101.32 0.10 15.20 30.39 45.59 60.79 75.99 91.18 101.32 0.10 15.20 30.39 45.59 60.79 75.99 91.18 101.32
1379.57 1453.58 1518.94 1578.62 1634.08 1685.40 1733.72 1764.43 1359.69 1434.72 1500.96 1561.16 1617.66 1669.47 1717.81 1749.44 1341.03 1416.90 1484.09 1545.04 1602.03 1654.76 1703.72 1735.30 1322.58 1399.36 1467.83 1529.71 1586.94 1640.55 1690.06 1721.65 1304.19 1382.43 1451.79 1514.65 1572.56 1626.52 1676.69 1708.53 1285.28 1365.24 1435.71 1499.39 1557.67 1612.63 1663.24 1695.25
292.92 292.86 292.86 292.86 292.86 292.86 292.86 292.86 298.20 298.15 298.15 298.15 298.15 298.15 298.15 298.14 303.18 303.14 303.14 303.13 303.13 303.13 303.13 303.13 308.16 308.13 308.13 308.13 308.13 308.12 308.12 308.12 313.14 313.10 313.10 313.10 313.10 313.10 313.10 313.10 318.31 318.29 318.28 318.30 318.30 318.29 318.29 318.30
0.10 15.20 30.39 45.59 60.79 75.99 91.18 101.32 0.10 15.20 30.39 45.59 60.79 75.99 91.18 101.32 0.10 15.20 30.39 45.59 60.79 75.99 91.18 101.32 0.10 15.20 30.39 45.59 60.79 75.99 91.18 101.32 0.10 15.20 30.39 45.59 60.79 75.99 91.18 101.32 0.10 15.20 30.39 45.59 60.79 75.99 91.18 101.32
1414.95 1482.08 1541.36 1596.67 1647.96 1695.54 1741.08 1769.94 1395.70 1463.92 1524.19 1580.09 1632.20 1680.47 1726.16 1755.32 1377.71 1446.76 1508.29 1564.48 1617.60 1666.26 1712.43 1741.88 1359.86 1430.19 1492.54 1549.56 1602.96 1652.55 1698.72 1728.68 1342.10 1413.95 1477.24 1534.76 1588.66 1639.09 1685.83 1715.79 1323.81 1397.11 1461.47 1519.87 1574.29 1625.11 1672.44 1702.47
1944
M. Dzida, P. Prusakiewicz / Fuel 87 (2008) 1941–1948
Table 3 Densities and specific heat capacities of diesel oil Ekodiesel Ultra, ON BIO10, and biodiesel measured at atmospheric pressure Ekodiesel Ultra
ON BIO10
Biodiesel
T (K)
q (kg m )
T (K)
q (kg m )
T (K)
q (kg m3)
273.15 278.15 283.16 288.15 293.16 298.16 303.16 308.16 313.15 318.16 323.15 328.16 333.16 338.16 343.15 348.16 353.15 358.15 363.15
845.882 842.363 838.845 835.331 831.816 828.303 824.789 821.277 817.760 814.242 810.719 807.191 803.660 800.127 796.584 793.038 789.482 785.919 782.350
273.15 278.15 283.16 288.16 293.15 298.16 303.16 308.16 313.16 318.16 323.15 328.16 333.16 338.16 343.15 348.16 353.15 358.15 363.15
855.514 851.980 848.449 844.923 841.398 837.874 834.348 830.825 827.301 823.776 820.248 816.715 813.178 809.643 806.099 802.550 798.996 795.433 791.865
273.15 278.15 283.16 288.16 293.15 298.16 303.16 308.16 313.16 318.16 323.15 328.15 333.15 338.16 343.15 348.15 353.16 358.15 363.15
894.390 890.731 887.086 883.445 879.814 876.185 872.563 868.942 865.325 861.709 858.096 854.483 850.871 847.261 843.648 840.034 836.416 832.796 829.181
T (K) 282.60 283.14 288.14 293.15 298.14 303.15 308.14 313.16 318.15 323.14 328.15 333.14 338.15 343.14 348.16 353.15 358.14 358.97
3
Cp (J kg1 K1) 1895.2 1897.8 1916.5 1935.0 1952.7 1971.4 1990.5 2009.9 2028.7 2049.1 2068.6 2089.3 2109.0 2128.7 2147.4 2166.7 2186.0 2188.8
T (K) 282.67 283.14 288.14 293.15 298.14 303.15 308.14 313.14 318.15 323.14 328.15 333.14 338.15 343.14 348.16 353.15 358.14 358.97
3
Cp (J kg1 K1) 1899.1 1901.2 1919.2 1937.2 1954.6 1971.8 1990.2 2008.5 2026.3 2045.4 2064.6 2083.6 2102.5 2121.2 2139.7 2157.6 2175.4 2178.6
T (K) 282.51 283.14 288.14 293.15 298.16 303.15 308.14 313.16 318.15 323.16 328.15 333.14 338.15 343.14 348.16 353.15 358.14 358.97
Cp (J kg1 K1) 1942.7 1945.6 1959.4 1971.9 1985.2 1998.5 2013.1 2028.1 2042.3 2057.6 2073.1 2088.9 2104.6 2119.6 2134.5 2148.5 2163.7 2165.9
Table 4 Coefficients of polynomial Eq. (1) for the speed of sound, density, and specific isobaric heat capacity, respectively and mean deviations from the regression line Fuel
c0 (m s1)
c1 (m s1 K1)
Ekodiesel Ultra ON BIO10 Biodiesel
2727.370 2749.493 2717.464
5.446246 5.564141 5.234999
Ekodiesel Ultra ON BIO10 Biodiesel
q0 (kg m3) 1032.289 1044.171 1095.576
q1 (kg m3 K1) 0.6653943 0.6788419 0.7461680
q2 Æ 105 (kg m3 K2) 6.27435 4.36606 3.48745
Ekodiesel Ultra ON BIO10 Biodiesel
h0 (J kg1 K1) 1035.764 1033.153 1547.327
h1 (J kg1 K2) 2.388477 2.573437 –
h2 Æ 103 (J kg1 K3) 2.31209 1.73489 5.47687
2.4. Specific heat capacity measurements The specific isobaric heat capacity was measured by a high sensitivity differential scanning calorimeter Micro DSC III,
c2 Æ 103 (m s1 K2)
du0 (m s1)
2.78037 3.02869 2.69146
0.009 0.002 0.004 dq (kg m3) 0.016 0.016 0.015 h2 Æ 106 (J kg1 K4) – – 1.82282
dCp (J kg1 K1) 0.8 0.8 1.0
manufactured by Setaram and based on the Tian-Calvet principle. The uncertainty of the measurements was estimated to be ± 0.15%. More details of the measuring set, calibration and measurement procedure can be found in [17].
M. Dzida, P. Prusakiewicz / Fuel 87 (2008) 1941–1948 Table 5 Coefficients of Eq. (2) and mean deviations from the regression line du Fuel
j
a1j (Kj MPa s m1)
Ekodiesel Ultra
0 1 2
0.2620890650
ON BIO10
Biodiesel
– 8.70894267Æ107
a2j (Kj MPa s2 m2) 2.67267574Æ104
du (m s1) 0.51
– 1.00348744Æ109
0 1 2
0.2822989252
2.48484590Æ104
– 1.05453199Æ106
– 7.36412361Æ1010
0 1 2
0.3301165804 – 1.34877055Æ106
2.56840831Æ104 1.95352490Æ107 –
0.26
0.26
1945
3. Measurement results The ultrasonic speeds in diesel oil Ekodiesel Ultra, ON BIO10, and biodiesel were measured at temperatures from 293 to 318 K in about 5 K steps and under pressures 0.1, 15, 30, 45, 60, 75, 90, and 101 MPa. The experimental values are listed in Table 2. The densities of the fuels and biofuels were measured under atmospheric pressure in the temperature range from 273 to 363 K in 5 K steps. The experimental values are collected in Table 3. The specific isobaric heat capacities were measured at atmospheric pressure and at temperatures from 283 to 359 K in about 0.02 K steps. Thus it has given ca. 3400 experimental points
Table 6 The calculated densities of diesel oil Ekodiesel Ultra, ON BIO10, and biodiesel at pressures up to 100 MPa and within the temperature limits 293 and 318 K
Table 7 The calculated isobaric specific heat capacities of diesel oil Ekodiesel Ultra, ON BIO10, and biodiesel at pressures up to 100 MPa and within the temperature limits 293 and 318 K
q (kg m3)
Cp (J kg1 K1)
p (MPa)
T (K) 293.15
p (MPa) 298.15
303.15
308.15
313.15
298.15
303.15
308.15
313.15
318.15
814.24 821.17 827.61 833.59 839.18 844.44 849.41 854.13 858.63 862.93 867.05
Ekodiesel Ultra 0.1a 1935 10 1930 20 1926 30 1923 40 1921 50 1919 60 1918 70 1917 80 1916 90 1916 100 1915
1953 1948 1945 1942 1940 1938 1937 1936 1935 1934 1934
1972 1967 1963 1960 1958 1956 1955 1954 1953 1953 1952
1991 1986 1982 1979 1977 1975 1974 1973 1972 1971 1971
2010 2005 2001 1998 1995 1994 1992 1991 1990 1990 1989
2030 2024 2020 2017 2014 2012 2011 2010 2009 2009 2008
827.31 833.97 840.18 845.97 851.40 856.52 861.38 866.00 870.41 874.63 878.69
823.78 830.60 836.96 842.87 848.41 853.62 858.56 863.26 867.74 872.03 876.14
ON BIO10 0.1a 10 20 30 40 50 60 70 80 90 100
1937 1932 1929 1926 1924 1922 1920 1919 1918 1917 1917
1955 1950 1946 1943 1941 1939 1938 1937 1936 1935 1934
1973 1968 1964 1961 1959 1957 1956 1954 1953 1953 1952
1991 1986 1982 1979 1977 1975 1973 1972 1971 1970 1970
2009 2004 2000 1997 1995 1993 1991 1990 1989 1988 1987
2028 2022 2018 2015 2013 2011 2009 2008 2007 2006 2005
865.33 871.66 877.61 883.20 888.48 893.48 898.25 902.80 907.16 911.35 915.39
861.71 868.19 874.28 879.98 885.35 890.45 895.29 899.91 904.34 908.59 912.68
Biodiesel 0.1a 10 20 30 40 50 60 70 80 90 100
1972 1968 1965 1963 1961 1960 1958 1957 1956 1956 1955
1986 1982 1979 1977 1975 1973 1972 1971 1970 1969 1968
2000 1996 1993 1990 1988 1987 1986 1984 1984 1983 1982
2014 2010 2007 2004 2002 2001 1999 1998 1997 1997 1996
2028 2024 2021 2018 2016 2015 2013 2012 2011 2011 2010
2043 2039 2035 2033 2031 2029 2028 2027 2026 2025 2024
Ekodiesel Ultra 0.1a 831.84 10 837.94 20 843.67 30 849.04 40 854.09 50 858.88 60 863.43 70 867.77 80 871.92 90 875.91 100 879.74
828.32 834.58 840.45 845.93 851.09 855.97 860.60 865.02 869.24 873.29 877.18
824.81 831.23 837.23 842.83 848.10 853.07 857.79 862.28 866.57 870.68 874.63
821.29 827.87 834.02 839.74 845.12 850.18 854.98 859.55 863.91 868.08 872.09
817.77 824.52 830.81 836.66 842.14 847.30 852.19 856.84 861.27 865.50 869.57
ON BIO10 0.1a 10 20 30 40 50 60 70 80 90 100
841.42 847.44 853.12 858.43 863.45 868.21 872.73 877.05 881.19 885.17 888.99
837.89 844.07 849.87 855.31 860.42 865.27 869.88 874.27 878.48 882.52 886.40
834.37 840.70 846.64 852.19 857.41 862.34 867.03 871.51 875.78 879.88 883.82
830.84 837.33 843.41 849.07 854.40 859.43 864.20 868.75 873.09 877.25 881.25
Biodiesel 0.1a 10 20 30 40 50 60 70 80 90 100
879.83 885.59 891.05 896.21 901.10 905.76 910.21 914.47 918.57 922.52 926.33
876.21 882.10 887.68 892.94 897.93 902.67 907.20 911.54 915.70 919.71 923.57
872.58 878.61 884.32 889.69 894.77 899.60 904.21 908.61 912.84 916.91 920.83
868.96 875.13 880.96 886.44 891.62 896.54 901.22 905.70 910.00 914.13 918.10
a
Calculated from Eq. (1).
318.15
T (K) 293.15
a
Calculated from Eq. (1).
1946
M. Dzida, P. Prusakiewicz / Fuel 87 (2008) 1941–1948
for each liquid. Therefore the values of isobaric specific heat capacity every 5 K are collected in Table 3. The dependencies of the speed of sound, density, and isobaric specific heat capacity on temperature at atmospheric pressure were approximated by second-order polynomials of the type: y¼
2 X
bj T j ;
ð1Þ
In order to calculate the density at p2, q, ap and Cp at p1 should be known beside the speed of the ultrasound as a function of pressure. The change of the heat capacity within the pressure limits Dp is given by: " # T 2 oap DC p ap þ Dp: ð5Þ q oT p
j¼0
where y is the speed of sound, u0, density, q, or isobaric specific heat capacity, Cp, at atmospheric pressure p0; bj are the polynomial coefficients bj = cj for the speed of sound, bj = qj for the density, and bj = hj for the isobaric specific heat capacity) calculated by the least squares method. The backward stepwise rejection procedure was used to reduce the number of non-zero coefficients. The coefficients and standard deviations from the regression lines are given in Table 4. The standard deviations are much smaller than the measurements accuracy. Since the sensitivity of the pressure gauge is lower than that of both the ultrasonic apparatus and the thermometer, the equation suggested by Sun et al. [7] was chosen in this work for smoothing out the speed of sound, pressure and temperature: p p0 ¼
m X n X i¼1
A numerical procedure proposed by Sun et al. [7] slightly modified by Marczak et al. [9] was applied. Details of the algorithm were discussed in previous work [9]. The densities and isobaric heat capacities of the fuels under test were calculated for temperatures from 293.15 to 318.15 K and for pressures up to 100 MPa. The calculated density and isobaric specific heat capacity are listed in Tables 6 and 7, respectively. From the densities and Table 8 The isentropic bulk modulus of diesel oil Ekodiesel Ultra, ON BIO10, and biodiesel at pressures up to 100 MPa and within the temperature limits 293 and 318 K bS Æ 106 (Pa) p (MPa)
T (K) 293.15
i
aij ðu u0 Þ T j ;
ð2Þ
j¼0
where aij are the polynomial coefficients calculated by the least squares method, u is the speed of sound at p > 0.1 MPa, u0 is the speed calculated from Eq. (1). The coefficients aij and the mean deviations from the regression lines are given in Table 5. The stepwise rejection procedure was used to reduce the number of the non-zero coefficients. 4. Physicochemical properties at elevated pressures The pressure dependence of the density is given by the following thermodynamic relationship: oq 1 T a2p ¼ 2þ ; ð3Þ op T u Cp where ap is the isobaric thermal expansion calculated from definition: ap ¼ ð1=qÞðoq=oT Þp . The change of the density of the liquid Dq, caused by a pressure increase from p1 to p2 at constant temperature, can be calculated by integration: ! Z p2 Z p2 a2p T 1 a2p T 1 Dq ¼ þ dp þ Dp: ð4Þ dp 2 u2 C p Cp p1 p1 u The approximate relationship Eq. (4) is sufficiently accurate, provided Dp is small enough, because the heat capacity depends rather slightly on pressure. Moreover, the first term on the right hand side of Eq. (4) is significantly larger than the second one since the latter results from the difference between the isentropic and isothermal compressibility that is small in comparison with both the compressibilities.
298.15
303.15
308.15
313.15
318.15
Ekodiesel Ultra 0.1 1561 10 1690 20 1815 30 1937 40 2055 50 2170 60 2283 70 2394 80 2504 90 2611 100 2718
1511 1640 1765 1886 2004 2119 2231 2342 2451 2559 2665
1463 1591 1716 1836 1954 2068 2181 2292 2401 2508 2615
1416 1544 1668 1788 1905 2020 2132 2243 2352 2459 2565
1370 1498 1622 1742 1858 1973 2085 2195 2304 2411 2518
1326 1453 1576 1696 1813 1927 2039 2149 2258 2366 2472
ON BIO10 0.1 10 20 30 40 50 60 70 80 90 100
1599 1728 1853 1974 2091 2206 2319 2430 2539 2646 2752
1549 1678 1802 1923 2040 2155 2267 2378 2487 2594 2700
1501 1629 1753 1874 1991 2105 2217 2328 2436 2543 2649
1453 1581 1706 1826 1943 2057 2169 2279 2387 2494 2600
1407 1535 1659 1779 1896 2010 2122 2232 2340 2447 2553
1362 1490 1614 1734 1851 1965 2076 2186 2295 2402 2507
Biodiesel 0.1 10 20 30 40 50 60 70 80 90 100
1759 1884 2005 2123 2238 2351 2461 2570 2677 2783 2887
1707 1832 1953 2071 2186 2298 2408 2517 2624 2729 2834
1656 1781 1903 2020 2135 2247 2357 2466 2572 2678 2782
1607 1732 1853 1971 2086 2198 2308 2416 2522 2627 2731
1559 1684 1805 1923 2038 2150 2259 2367 2474 2579 2682
1511 1637 1759 1876 1991 2103 2213 2321 2427 2532 2635
M. Dzida, P. Prusakiewicz / Fuel 87 (2008) 1941–1948
5. Discussion The densities of the fuels and biofuels increase monotonically with increasing pressure and decreasing temperature. Temperature dependence of densities of the fuels and biofuels is described by the second-order polynomials (Table 4) while pressure dependence is described by the thirdorder polynomial. Respective slopes of the curves are practically the same. Nevertheless, densities of ON BIO 10 and diesel oil Ekodiesel Ultra are lower than densities of biodiesel, over the whole temperature and pressure range, by about 4.6–4.2% and 5.8–5.3%, respectively. The isobaric specific heat capacity of the fuels decreases with increasing pressure and decreasing temperature. However, the effect of pressure is smaller than that of temperature (Table 7). Difference between isobaric specific heat capacities of diesel oil Ekodiesel Ultra and ON BIO 10 are negligible over the whole temperature and pressure range. The heat capacity of biodiesel is higher than that of diesel oil Ekodiesel Ultra and ON BIO 10 by 2% at the most (Tables 3 and 7). Isobaric thermal expansivity of biodiesel decreases with pressure slightly less than that of the diesel oil Ekodiesel Ultra. It is approximately independent of temperature and composition of the fuel at pressures 40 ± 5 MPa (Fig. 1). The crossing point of isobaric thermal expansion is characteristic for simple liquids. For example, the isotherms of ap of hexane cross each other at 65 ± 2 MPa [18]. However, taking into account the uncertainty of the thermal expansion coefficient equals ± 1%, there is not sig-
nificant difference between ap of all the tested fuels. At temperatures from 293.15 to 318.15 K and pressures from atmospheric one to 100 MPa, the bulk moduli of biodiesel and ON BIO 10 were higher than that of diesel oil Ekodiesel Ultra by 6.2–14.0% and 1.3–2.7%, respectively (Fig. 2 a). Temperature change by 5 K changes the bulk modulus by 1.9–3.1% and by 1.8–3.3% for biodiesel and ON BIO 10, respectively. Thus, the effect of addition of biodiesel to the petroleum fuel on the bulk modulus is of the same order as that resulting from the slight change of temperature (Fig. 2b). With temperature increasing from 293.15 to 318.15 K, at atmospheric pressure, the bulk moduli decrease in a similar way for the diesel oil Ekodiesel Ultra and ON BIO 10: from 1561 to 1325 MPa and from 1599 to 1361 MPa, respectively. These ranges almost overlap each other. Thus, although the bulk modulus of ON BIO 10 is indeed higher than that of diesel oil Ekodiesel Ultra, the difference is rather small and can be compensated by temperature. Similar conclusion is valid for pressures up to100 MPa. These results show that effect of temperature on the bulk modulus should be taken into account over the whole pressure range.
3000
2500
Β S / MPa
speeds of sound, the adiabatic bulk modulus (the reciprocal of adiabatic compressibility jS) were calculated by the 2 1 Laplace formula: B1 S ¼ jS ¼ ðqu Þ . Results of the calculations are given in Table 8. The isothermal compressibility ðjT ¼ jS þ a2p T =qC p Þ and isothermal bulk modulus ðjT ¼ B1 T Þ can be also calculated from the above data.
1947
2000
1500
1000
0
20
40
60
80
100
60
80
100
p / MPa 0.9
3000
293.15 K
2500
Β S / MPa
α p · 10 3/ K -1
0.8
0.7
2000
1500 318.15 K
0.6
0
20
40
60
80
100
p / MPa
Fig. 1. Isobaric thermal expansion of: diesel oil Ekodiesel Ultra: (d) 293.15 K, (s) 318.15 K; ON BIO 10: (j) 293.15 K, (h) 318.15 K; Biodiesel (m) 293.15 K, (4) 318.15 K, lines calculated from empirical P function: ap ¼ 3i¼0 ai pi .
1000 0
20
40
p / MPa
Fig. 2. Isentropic bulk modulus of : a) diesel oil Ekodiesel Ultra: (d) 293.15 K, (s) 318.15 K; ON BIO 10: (j) 293.15 K, (h) 318.15 K; Biodiesel (m) 293.15 K, (4) 318.15 K; b) (s) diesel oil Ekodiesel P Ultra; (h) ON BIO 10; lines calculated from empirical function: Bs ¼ 3i¼0 ai pi .
1948
M. Dzida, P. Prusakiewicz / Fuel 87 (2008) 1941–1948
6. Conclusions 1. Investigated fuels and biofuels differ in density. That facilitates telling one fuel from another. 2. Isobaric specific heat capacity of biodiesel is higher than those of diesel oil Ekodiesel Ultra and ON BIO 10 that are approximately the same in the temperature range from 293.15 to 318.15 K and pressures up to 100 MPa. 3. Isobaric thermal expansions of the three investigated fuels are approximately equal to each other and independent of temperature at pressures 40 ± 5 MPa. 4. Increase of bulk modulus of petrodiesel due to addition of 20 vol.% of biodiesel is small. It is approximately the same as that caused by decrease of temperature by ca. 5 K. That probably holds even for pressures much higher than 100 MPa attained in this experiment. Acknowledgements The authors are profoundly indebted to Dr. W. Marczak for critical reading of the manuscript and Ms. M. Bucek for the participation in the density measurements. References [1] Tat ME, Van Gerpen JH, Soylu S, Canakci M, Monyem A, Wormley SJ. The speed of sound and isentropic bulk modulus of biodiesel at 21 °C from atmospheric to 35 MPa. Am Oil Chem Soc 2000;77:285–9. [2] Tat ME, Van Gerpen JH. Effect of temperature and pressure on the speed of sound and isentropic bulk modulus of mixtures of biodiesel and diesel fuel. J Am Oil Chem Soc 2003;80:1127–30. [3] Szybist JP, Boehman AL, Taylor JD, McCormick RL. Evaluation of formulation strategies to eliminate the biodiesel NOx effect. Fuel Proc Technol 2005;86:1109–26. [4] Boehman AL, Morris D, Szybist JP. The impact of the bulk modulus of diesel fuels on fuel injection timing. Energ Fuels 2003;18:1877–82. [5] Kegl B. Numerical analysis of injection characteristics using biodiesel fuel. Fuel 2006;85:2377–87.
[6] Forson FK, Oduro EK, Hammond-Donkoh E. Performance of jatropha oil blends in a diesel engine. Renew Energ 2004;29:1135–45. [7] Sun TF, Ten Seldam CA, Kortbeek PJ, Trappeniers NJ, Biswas SN. Acoustic and thermodynamic properties of ethanol from 273.15 to 333.15 K and up to 280 MPa. Phys Chem Liq 1988;18:107–16. [8] Davis LA, Gordon RB. Compression of mercury at high pressure. J Chem Phys 1967;46:2650–60. [9] Marczak W, Dzida M, Ernst S. Determination of the thermodynamic properties of 1-propanol and 1-hexanol from speed of sound measurements under elevated pressures. High Temp – High Press 2000;32:283–92. [10] Plantier F, Daridon JL, Lagourette B. Nonlinear parameter (B/A) measurements in methanol, 1-butanol and 1-octanol for different pressures and temperatures. J Phys D: Appl Phys 2002;35:1063–7. _ [11] Dzida M, Zak A, Ernst S. Thermodynamic and acoustic properties of binary mixtures of alcohols and alkanes. I. Speed of sound in (ethanol + n-heptane) under elevated pressures. J Chem Thermodyn 2005;37:405–14. [12] Dzida M, Marczak W. Thermodynamic and acoustic properties of binary mixtures of alcohols and alkanes. II. Density and heat capacity of (ethanol + n-heptane) under elevated pressures. J Chem Thermodyn 2005;37:826–36. [13] Tate RE, Watts KC, Allen CAW, Wilkie KI. The densities of three biodiesel fuels at temperatures up to 300 °C. Fuel 2006;85:1004–9. [14] Yamate K, Ueta A, Shimamoto Y. Influence of physical and chemical properties of biodiesel fuel on injection, combustion and exhaust emission characteristics in a DI-CI Engine. In: The fifth international symposium on diagnostics and modeling of combustion in internal combustion engines (COMODIA 2001), July 1–4, 2001, Nagoya; p. 402–9. _ [15] Zak A, Dzida M, Zore˛bski M, Ernst S. A high pressure system for measurements of the speed of sound in liquids. Rev Sci Instrum 2000;71:1756–65. [16] Dzida M, Chora˛z_ ewski M, Zore˛bski M, Man´ka R. Modification of a high pressure device for speed of sound measurements in liquids. J Phys IV 2006;137:203–7. [17] Go´ralski P, Tkaczyk M, Chora˛z_ ewski M. Heat capacities of a,xdichloroalkanes at temperatures from 28415 K to 353.15 K and a group additivity analysis. J Chem Eng Data 2003;48:492–6. [18] Randzio SL, Grolier J-PE, Quint JR, Eatough DJ, Lewis EA, Hansen LD. n-hexane as a model for compressed simple liquids. Int J Thermophys 1994;15:415–41.