Experimental values of densities (ρ) and viscosities (η) in the binary mixtures of n-octane, n-decane, n-dodecane, and n-tetradecane with octan-2-ol are presented over the whole range of mixture composition at T=298.15 K. From these data, excess molar volume (VmE), deviations in viscosity (Δη), and excess Gibbs free energy of activation ΔG∗E have been calculated. These results were fitted to Redlich-Kister polynomial equations to estimate the binary coefficients and standard errors. Jouyban-Acree model is used to correlate the experimental values of density and viscosity at T=298.15 K. The values of VmE have been analyzed using Prigogine-Flory-Patterson (PFP) theory. The results of the viscosity composition are discussed in the light of various viscosity equations suggested by Grunberg-Nissan, Tamara and Kurata, Hind et al., Katti and Chaudhri, Heric, Heric and Brewer, and McAllister multibody model. The values of Δlnη have also been analyzed using Bloomfield and Dewan model. The experiments on the constituted binaries are analyzed to discuss the nature and strength of intermolecular interactions in these mixtures.
1. Introduction
Density and viscosity data for liquid mixtures are important from practical and theoretical points of view. Experimental measurements of these properties for binary mixtures have gained much importance in many chemical industries and engineering disciplines [1]. Experimental liquid viscosities of pure hydrocarbons and their mixtures are needed for the design of chemical processes where heat and mass transfer and fluid mechanics are important. Prediction of the liquid behavior of hydrocarbon mixture viscosities is not yet possible within the experimental uncertainty. Therefore, experimental measurements are needed to understand the fundamental behavior of this property and then to develop new models [2]. Alkanes are important series of homologous, nonpolar, and organic solvents. They have often been used in the study of solute dynamics because their physicochemical properties as a function of chain length are well-known [3]. They are also employed in a large range of chemical processes [4]. The physicochemical properties play an important role in the understanding of several industrial processes. Properties such as viscosity or surface tension are required in many empirical equations for different operations such as mass and heat transfer processes. For example, it is necessary to know the mass transfer coefficient to design gas-liquid contactors. To determine the equations that modelize the mass transfer process requires knowledge of the density, viscosity, and surface tension of the liquid phase [5].
Viscosities and excess molar volumes of binary mixtures of methylbenzene [6] and methylcyclohexane [7] with octan-2-ol at T=298.15 K are only reported. To the best of our knowledge, the properties of the binary mixtures of n-octane, n-decane, n-dodecane, and n-tetradecane with octan-2-ol have not been reported earlier.
In the present paper, we report density and viscosity data for the binary mixtures of n-octane, n-decane, n-dodecane, and n-tetradecane with octan-2-ol T=298.15 K. This work will also provide a test of various semiempirical equations to correlate viscosity of binary mixtures. The types of used relations are Tamara and Kurata, Heric, Heric and Brewer, Hind et al., Katti and Chaudhri, and McAllister multibody model.
2. Experimental
Chemicals used in the present study were of analytical grade and supplied by S.D. Fine Chemicals Pvt Ltd. Mumbai with quoted mass fraction purities: n-octane (>0.99), n-decane (>99.6), n-dodecane (>99.8), and n-tetradecane (>99.7). Octan-2-ol (purity >99.3) was supplied by E-Merck. Prior to use, all liquids were stored over 0.4 nm molecular sieves to reduce the water content and were degassed. The binary mixtures of varying composition were prepared by mass in special air-tight bottles. The masses were recorded on a Mettler balance to an accuracy of ±1×10-5 g. Care was taken to avoid evaporation and contamination during mixing. The estimated uncertainty in mole fraction was <1×10-4.
The densities of the solutions were measured using a single capillary pycnometer made up of borosil glass with a bulb of 8 cm3, and capillary with internal diameter of 0.1 cm was chosen for the present work. The detailed pertaining to calibration, experimental set up, and operational procedure has been previously described [6–10]. An average of triplicate measurement was taken in to account. The reproducibility of density measurement was ±3×10-5 g/cm3.
The dynamic viscosities were measured using an Ubbelohde suspended level viscometer [6–10] calibrated with conductivity water. An electronic digital stop watch with readability of ±0.01 s was used for the flow time measurements. At least three repetitions of each data reproducible to ±0.05 s were obtained, and the results were averaged. Since all flow times were greater than 300 s and capillary radius (0.1 mm) was far less than its length (50 to 60) mm, the kinetic energy and end corrections, respectively, were found to be negligible. The uncertainties in dynamic viscosities are of the order of ±0.003 mPa·s.
The purity of the samples and accuracy of data were checked by comparing the measured densities and viscosities of the pure compounds with the literature values, which are given in Table 1. Thus, our results are in good agreement with those listed in the literature.
Comparison of experimental densities (ρ), viscosities (η), and speeds of sound (u) of pure liquids with the literature values at T = 298.15 K.
Experimental values of densities (ρ) and viscosities (η) of mixtures at 298.15 K are listed as a function of mole fraction in Table 2. The density values have been used to calculate excess molar volumes (VE) using the following equation:
(1)VmE/m3·mol-1=(x1M1+x2M2)ρ12-(x1M1ρ1)-(x2M2ρ2),
where ρ12 is the densities of the mixture and x1, M1, ρ1, and x2, M2, ρ2 are the mole fractions, the molecular weights, and the densities of pure components 1 and 2, respectively.
Densities (ρ), viscosities (η), excess molar volumes VE, viscosities deviations Δη, and excess Gibbs free energy ΔG*E of binary mixtures at T = 298.15 K.
x1
ρ×10-3
VmE×106
η
Δη
ΔG*E
kg·m−3
m3·mol−1
mPa·s
mPa·s
KJmol−1
n-Octane(1) + octan-2-ol(2)
0
0.81705
0
6.429
0
0
0.0487
0.8114
−0.051
5.956
−0.185
115
0.1002
0.8054
−0.097
5.301
−0.535
149
0.1498
0.79959
−0.133
4.654
−0.888
137
0.2006
0.79361
−0.16
3.998
−1.244
78
0.2516
0.78756
−0.177
3.364
−1.576
−30
0.3003
0.78176
−0.183
2.796
−1.856
−183
0.3538
0.77534
−0.18
2.228
−2.108
−410
0.4003
0.76975
−0.169
1.79
−2.27
−661
0.4466
0.76415
−0.151
1.411
−2.375
−960
0.5001
0.75767
−0.123
1.049
−2.421
−1359
0.5517
0.75142
−0.091
0.778
−2.387
−1776
0.6002
0.74554
−0.057
0.59
−2.287
−2157
0.6539
0.73905
−0.018
0.454
−2.106
−2469
0.6994
0.73357
0.013
0.39
−1.901
−2560
0.7449
0.72811
0.042
0.366
−1.656
−2431
0.8
0.72157
0.068
0.376
−1.32
−2019
0.8465
0.7161
0.081
0.402
−1.018
−1571
0.9003
0.70985
0.081
0.434
−0.668
−1034
0.9521
0.70394
0.061
0.441
−0.355
−670
1
0.69867
0
0.512
0
0
n-Decane(1) + octan-2-ol(2)
0
0.81705
0
6.429
0
0
0.0554
0.81077
0.036
5.222
−0.898
−233
0.0998
0.80583
0.071
4.638
−1.234
−301
0.1554
0.79974
0.119
3.972
−1.589
−403
0.1999
0.79499
0.158
3.491
−1.822
−497
0.2554
0.7892
0.205
2.951
−2.052
−632
0.2999
0.78467
0.239
2.564
−2.19
−755
0.3555
0.77916
0.275
2.136
−2.308
−927
0.3998
0.7749
0.298
1.836
−2.361
−1078
0.4554
0.76969
0.318
1.509
−2.377
−1285
0.4999
0.76564
0.326
1.284
−2.354
−1461
0.5555
0.76074
0.326
1.046
−2.281
−1690
0.5998
0.75694
0.318
0.889
−2.191
−1871
0.6555
0.7523
0.298
0.73
−2.039
−2080
0.7
0.7487
0.274
0.63
−1.89
−2223
0.7554
0.74432
0.236
0.538
−1.673
−2337
0.7998
0.74091
0.2
0.487
−1.476
−2362
0.8554
0.73674
0.149
0.45
−1.203
−2281
0.8999
0.73347
0.106
0.438
−0.966
−2126
0.9554
0.72948
0.052
0.445
−0.649
−1811
1
0.72635
0
0.845
0
0
n-Dodecane(1) + octan-2-ol(2)
0
0.81705
0
6.429
0
0
0.0554
0.81098
0.101
5.513
−0.634
−159
0.0999
0.80637
0.167
4.916
−1.004
−296
0.1554
0.80091
0.24
4.254
−1.383
−404
0.1998
0.79676
0.291
3.786
−1.625
−517
0.2555
0.7918
0.345
3.272
−1.856
−659
0.2998
0.78805
0.38
2.916
−1.986
−770
0.3554
0.78355
0.415
2.532
−2.087
−902
0.3998
0.78012
0.435
2.27
−2.123
−999
0.4554
0.77603
0.451
1.995
−2.115
−1102
0.4999
0.7729
0.456
1.813
−2.071
−1167
0.5555
0.76916
0.453
1.627
−1.973
−1219
0.5998
0.76631
0.444
1.508
−1.866
−1236
0.6555
0.76288
0.422
1.391
−1.699
−1221
0.6998
0.76027
0.398
1.32
−1.545
−1180
0.7554
0.75713
0.359
1.253
−1.328
−1095
0.7998
0.75473
0.32
1.213
−1.142
−1005
0.8555
0.75185
0.263
1.176
−0.896
−869
0.8999
0.74965
0.21
1.152
−0.694
−750
0.9555
0.747
0.134
1.124
−0.439
−599
0.9997
0.74519
0
1.337
0
0
n-Tetradecane(1) + octan-2-ol(2)
0
0.81705
0
6.429
0
0
0.0555
0.81134
0.126
5.619
−0.568
−144
0.0999
0.80721
0.193
5.163
−0.831
−230
0.1555
0.80239
0.268
4.632
−1.121
−330
0.1998
0.7988
0.321
4.238
−1.321
−417
0.2554
0.79459
0.376
3.786
−1.533
−531
0.2998
0.79144
0.413
3.455
−1.67
−628
0.3555
0.78774
0.449
3.081
−1.802
−751
0.3999
0.78496
0.47
2.813
−1.876
−850
0.4554
0.78171
0.487
2.519
−1.93
−968
0.4999
0.77925
0.492
2.314
−1.941
−1054
0.5555
0.77637
0.49
2.098
−1.915
−1144
0.5998
0.7742
0.48
1.957
−1.863
−1196
0.6554
0.77163
0.457
1.82
−1.759
−1226
0.7
0.76969
0.432
1.741
−1.644
−1217
0.7555
0.76741
0.39
1.683
−1.46
−1154
0.8
0.76569
0.348
1.668
−1.282
−1060
0.8555
0.76365
0.287
1.688
−1.02
−887
0.9
0.76211
0.23
1.736
−0.779
−703
0.9553
0.76029
0.15
1.835
−0.439
−424
0.9997
0.75914
0
2.081
0
0
The variation in excess molar volumes, VmE with mole fraction of the binary mixtures of n-octane, n-decane, n-dodecane, and n-tetradecane with octan-2-ol at T = 298.15 K, is displayed in Figure 1. The VmE curve for the mixture of n-octane and octan-2-ol is sigmoidal and tends to change to positive values at higher mole fractions (x1≥0.65) of n-octane while VmE values for n-decane, n-dodecane, and n-tetradecane with octan-2-ol mixtures show positive deviation over the entire composition range.
Plot of excess molar volumes (VmE) against mole fraction of octan-2-ol with (◊) n-octane; (□) n-decane; (▵) n-dodecane, and n-tetradecane (○) at T = 298.15 K. The corresponding dotted (- - -) curves have been derived from PFP theory.
Generally, VmE can be considered as arising from three types of interactions between component molecules of liquid mixtures [11–13]: (1) physical interactions consisting of mainly of dispersion forces or weak dipole-dipole interaction making a +ve contribution, hereby the contraction volume and compressibility of the mixtures, (2) chemical or specific interactions, which include charge transfer, forming of H-bonds and other complex forming interactions, resulting in a −ve contribution, and (3) structural contribution due to differences in size and shape of the component molecules of the mixtures, due to fitting of component molecules into each other’s structure, hereby reducing the volume and compressibility of the mixtures, resulting in a −ve contribution.
The large positive VmE values (Figure 1) for n-decane, n-dodecane, and n-tetradecane with octan-2-ol mixture are attributed to the breaking up of three-dimensional H-bonded network of octan-2-ol due to the addition of solute, which is not compensated by the weak interactions between unlike molecules. The VmE exhibits an inversion in sign in the mixtures f (n-octane + 2-octanol). The values of VmE are negative in the lower region of x1 due to interstitial accommodation is more as compared to the de-clustering of 1-octanol molecules and beyond (x1≥0.65) as the amount of n-octane increases in the mixture, due to dispersion forces, thereby making a positive contribution to VmE.
The viscosity deviations (Δη) were calculated using
(2)Δη/mPa·s=η12-x1η1-x2η2,
where η12 is the viscosities of the mixture and x1, x2 and η1, η2 are the mole fraction and the viscosities of pure components 1 and 2, respectively.
Figure 2 depicts the variation of Δη with mole fraction of the binary mixtures of n-octane, n-decane, n-dodecane, and n-tetradecane with octan-2-ol at T = 298.15 K.
Plot of viscosity deviations (Δη) against mole fraction of octan-2-ol with (◊) n-octane; (□) n-decane; (▵) n-dodecane and n-tetradecane (○) at T = 298.15 K.
The deviations viscosity may be generally explained by considering the following factors [14]. (1) The difference in size and shape of the component molecules and the loss of dipolar association to a decrease in viscosity; (2) specific interactions between unlike molecules such as H-bond formation and charge transfer complexes may cause increase in viscosity in mixtures rather than in pure component. The former effect produces negative in excess viscosity, and latter effect produces positive in excess viscosity. Positive values of Δη are indicative of strong interactions whereas negative values indicate weaker interactions [15].
The negative deviations in viscosity support the main factor of breaking of the self-associated alcohols and weak interactions between unlike molecules. The negative values viscosity deviation decreases in the following sequence: n-octane >n-decane >n-dodecane >n-tetradecane.
Excess Gibbs free energies of activation of viscous flow ΔG*E for binary mixtures can be calculated as
(3)ΔG*E=RT[ln(ηυη2υ2)-x1ln(η1υ1η2υ2)],
where υ is the molar volume of the mixture, υi is the molar volume of the pure component, R is the gas constant, T is the absolute temperature, and η is the dynamic viscosity of the mixture, respectively. ηi is the dynamic viscosity of the pure component i and x1 the mole fraction in component. The ΔG*E values of all binary systems are shown in Table 2. The values of ΔG*E for all binary mixtures are negative over entire mole fraction. According to Meyer et al. [16], negative values of ΔG*E correspond to the existence of solute-solute association.
The excess molar volumes and deviations in viscosity were fitted to Redlich and Kister [17] equation of the type
(4)Y=x1x2∑inai(x1-x2)i,
where Y is either VE or Δη and n is the degree of polynomial. Coefficient ai was obtained by fitting (5) to experimental results using a least-squares regression method. In each case, the optimum number of coefficients is ascertained from an examination of the variation in standard deviation (σ).
σ was calculated using the relation
(5)σ(Y)=[∑(Yexpt-Ycalc)2N-n]1/2,
where N is the number of data points and n is the number of coefficients. The calculated values of the coefficients ai along with the standard deviations (σ) are given in Table 3.
Coefficients ai of (4) and corresponding standard deviation (σ) of (5) at T = 298.15 K.
System
a0
a1
a2
a3
σ
n-Octane + octan-2-ol
VmE×106/(m3·mol−1)
0.125
0.1577
−0.41
0.1523
0.943
Δη/(mPa·s)
−9.6807
0.4921
4.661
−2.7989
0.024
n-Decane + octan-2-ol
VmE×106/(m3·mol−1)
0.129
0.167
−0.462
0.1437
0.0023
Δη/(mPa·s)
−8.4704
2.213
−6.8274
−1.0134
0.1273
n-Dodecane + octan-2-ol
VmE×106/(m3·mol−1)
0.363
−1.276
0.289
0.812
0.529
Δη/(mPa·s)
−7.989
3.365
−2.905
8.812
0.060
n-Tetradecane + octan-2-ol
VmE×106/(m3·mol−1)
0.880
−0.116
0.009
0.708
0.020
Δη/(mPa·s)
−6.322
−4.730
−10.25
15.66
0.377
4. Theoretical Analysis4.1. Semiempirical Models for Analyzing Viscosity of Liquid Mixtures
Several empirical and semiempirical relations have been used to represent the dependence of viscosity on concentration of components in binary liquid mixtures, and these are classified according to the number of adjustable parameters used to account for the deviation from some average [18, 19]. We will consider here some of the most commonly used semiempirical models for analyzing viscosity of liquid mixtures based on one, two, and three parameters. An attempt has been made to check the suitability of equations for experimental data fits by taking into account the number of empirical adjustment coefficients.
The equation of Grunberg-Nissan, Tamara and Kurata Hind et al., and Katti and Chaudhri has one adjustable parameter.
Gruenberg-Nissan provided the following empirical equation containing one adjustable parameter [20]. The equation is
(6)lnη12=x1lnη1+x2lnη2+x1x2G12,
where G12 may be regarded as a parameter proportional to the interchange energy also an approximate measure of the strength of the interaction between the components.
The one-parameter equation due to Tamura and Kurata [21] gave the equation of the form
(7)ηm=x1η1Φ+x2η2Φ2+2(x1x2ΦΦ2)1/2T12,
where Φ and Φ2 are the volume fractions of components 1 and 2, respectively; T12 is Tamura and Kurata constant.
Hind et al. [22] proposed the following equation:
(8)ηm=x12η1+x22η2+2x1x2H12,
where H12 is attributed to unlike pair interactions.
Katti and Chaudhri [23] derived the following equation
(9)lnηV=x1lnηV10+x2lnηV20+x1x2WvisRT,
where Wvis is an interaction term and υi is the molar volume of pure component i.
Heric [24] expression is
(10)lnηm=x1lnη1+x2lnη2+x1lnη1+x2lnη2+ln(x1η1+x2η2)+δ12,
where δ12=α12x1x2 is a term representing departure from a noninteracting system and α12=α21 is the interaction parameter. Either α12 or α21 can be expressed as a linear function of composition:
(11)α12=β12′+β12′′(x1-x2).
From an initial guess of the values of coefficients β12′ and β12′′, the values of α12 are computed.
Heric and Brewer [25] equation is
(12)lnν=x1lnν1+x2lnν2+x1lnM1+x2lnM2-ln[x1M1+x2M2]+x1x2[α12+α21(x1-x2)].M1 and M2 are molecular weights of components 1 and 2, and α12 and α21 are interaction parameters, and other terms involved have their usual meaning. α12 and α21 are parameters, which can be calculated from the least-squares method.
McAllister’s multibody interaction model [26] was widely used to correlate kinematic viscosity (ν) data. The two-parameter McAllister equation based on Eyring’s [27] theory of absolute reaction rates has taken into account interaction of both like and unlike molecules by two-dimensional three-body model. The three-body interaction model is
(13)lnνm=x13lnν1+3x12x2lnZ12+3x1x22lnZ21+x23lnν2-ln[x1+(x2M2)M1]+3x12x2ln[23+M2(3M1)]+3x1x22ln[13+2M2(3M1)]x23ln(M2M1).
And four-body model was given by
(14)lnνm=x14lnν1+4x13x2lnZ1112+6x12x22lnZ1122+4x1x23lnZ2221+x24lnν2-ln[x1+x2(M2M1)]+4x13x2ln[3+(M2/M1)4]+6x12x22ln[1+(M2/M1)2]+4x1x23ln[1+(3M2/M1)4]+x24ln(M2M1),
where Z12, Z21, Z1112, Z1122, and Z2221 are interaction parameters and Mi and νi are the molecular mass and kinematic viscosity of pure component i, respectively.
The correlating ability of each of (6)–(14) was tested as well as their adjustable parameters and standard deviations (σ):
(15)σ(%)=[(1(n-k)∑100(ηexptl-ηcalcd)ηexptl)2]1/2,
where n represents the number of data points and k is the number of numerical coefficients given in Table 4. The interaction parameter G12 is negative for binary systems. Nigam and Mahl [28] concluded from the study of binary mixtures that (1) if Δη>0, G12>0 and magnitude of both are large then strong specific interaction would be present; (2) if Δη<0, G12>0 then weak specific interaction would be present; (3) if Δη<0, G12<0 magnitude of both are large then the dispersion force would be dominant. According to Fort and Moore [29] and Ramamoorty and Varadachari [30], system exhibits strong interaction if the G12 is positive; if it is negative they show weak interaction. On this basis, we can say that there is a weak interaction in the system studied.
Adjustable parameters of (6)–(13) and standard deviations of binary mixture viscosities for x1 n-alkanes + (1-x1) octan-2-ol at T = 298.15 K.
Equation
System including n-alkanes + octan-2-ol
n-Octane
n-Decane
n-Dodecane
n-Tetradecane
Grunberg-Nissan
G12
−2.398
−3.226
−1.997
−2.002
σ
3.759
3.904
0.858
1.075
Tamura and Kurata
T12
0.872
−1.117
0.184
0.675
σ
1.737
1.884
0.812
0.782
Hind et al.
H12
−0.91
−1.409
−0.401
0.079
σ
1.758
2.145
1.172
0.711
Katti and Chaudhri
Wvis
−2.402
−3.198
−1.96
−1.871
σ
3.755
3.905
0.857
1.081
Heric and Brewer
α
12
−2.381
−3.936
−2.109
−2.104
α
21
−4.123
−5.529
−1.506
−0.827
σ
−0.926
2.449
−0.559
0.866
McAllister’s three-body
Z12
0.185
0.094
0.879
1.439
Z21
6.363
7.068
3.929
3.544
σ
0.926
2.449
0.559
0.866
McAllister’s four-body
Z1112
0.269
0.51
0.817
1.140
Z1122
0.936
19.584
3.639
7.334
Z2221
6.926
1.998
3.598
2.871
σ
5.149
5.865
0.622
1.622
Interaction parameter Wvis shows almost the same trend as that of G12. In fact, one could say that the parameters G12 and Wvis exhibit almost similar behaviour, which is not unlikely in view of logarithmic nature of both equations.
Tamara and Kurata and Hind et al. represent the binary mixture satisfactory as compared to Gruenberg-Nissan and Katti and Chaudhri. Use of three parameters equation reduces the σ values significantly below that of two parameters equation. From this study, it can be concluded that the correlating ability significantly improves for these nonideal systems as number of adjustable parameters is increased. From Table 4, it is clear that McAllister’s three-body interaction model is suitable for correlating the kinematic viscosities of the binary mixtures studied.
4.2. Prigogine-Flory-Patterson (PFP) Theory
The Prigogine-Flory-Patterson (PFP) theory [31–34] has been commonly employed to estimate and analyze excess thermodynamic functions theoretically. This theory has been described in details by Patterson and coworkers [35, 36]. According to PFP theory, VmE can be separated into three factors: (1) an interactional contribution, VmE (int.) (2) a free volume contribution, VmE(fv), and (3) an internal pressure contribution, VmE(P*). The expression for these three contributions are given as VmE:
(16)VmE(int)=[(υ1/3-1)υ2/3ψ1θ2(4/3υ-1/3-1)P1*χ12],VmE(fv)=[(υ1-υ2)2(14/9υ-1/3-1)ψ1ψ2][(14/9υ-1/3-1)υ],VmE(P*)=[(υ1-υ2)(P1*+P2*)ψ1ψ2](P1*ψ1+P2*ψ2).
Thus, the excess molar volume VmE is given as
(17)VmE(x1V1+x2V2)=VmE(int)-VmE(fv)+VmE(P*),
where ψ, θ, and P* represent the contact energy fraction, surface site fraction, and characteristic pressure, respectively, and are calculated as
(18)ψ1=(1-ψ2)=Φ1P1*(Φ1P1*+Φ2P2*),(19)θ2=(1-θ1)=Φ2[Φ1(V2*/V1*)],(20)P*=Tυ2ακT.
The details of the notations and terms used in (16)–(19) may be obtained from the literature [31–34, 37, 38]. The other parameters pertaining to pure liquids and the mixtures are obtained from the Flory theory [7, 31, 38] while α and κT values are taken from the literature [39–44]. The interaction parameter χ12 required for the calculation of VmE using PFP theory has been derived by fitting the VmE expression to the experimental equimolar value of VmE for each system under study.
The values of χ12, θ2, three PFP contributions interactional, free volume, P* effect, and experimental and calculated (using PFP theory) VmE values at near equimolar composition are presented in Table 5. Study of the data presented in Table 6 reveals that the interactional and free volume contributions are positive, whereas P* contributions are negative for all the three systems under investigation. For these binary mixtures, it is only the interactional contribution which dominates over the remaining two contributions. The P* contribution, which depends both on the differences of internal pressures and differences of reduced volumes of the components, has little significance for the studied binary mixtures as compared to the other two.
Flory parameters of the pure compounds.
Components
106V*/(m3·mol−1)
106P* (J·m−3)
T*/K
Octan-2-ol
129.9790
535
5563
n-Octane
127.4844
444
4826
n-Decane
155.6091
453
5091
n-Dodecane
183.7700
455
5290
n-Tetradecane
212.1200
460
5479
Calculated values of the three contributions to the excess molar volume from the PFP theory with interaction parameter T = 298.15 K.
Component
χ12×106 (J·m−3)
VE×106 (m3·mol−1) at x=0.5
Calculated contribution
Experimental
PEP
VE×106 (int)
VE×106 (fv)
VE×106 (P*)
n-Octane + octan-2-ol
11.59
0.123
0.124
0.0331
0.1087
−0.4194
n-Decane + octan-2-ol
12.57
0.326
0.328
0.5999
0.0644
−0.2283
n-Dodecane + octan-2-ol
16.57
0.456
0.458
0.5524
0.0175
−0.0666
n-Ttetradecane + octan-2-ol
15.08
0.492
0.494
0.5157
0.0116
−0.0187
Furthermore, in order to check whether χ12, derived from nearly equimolar VmE values, can predict the correct composition dependence, VmE has been calculated theoretically using χ12 over the entire composition range. The theoretically calculated values are plotted in Figure 1 for comparison with the experimental results. Figure 1 show that the PFP theory is quite successful in predicting the trend of the dependence of VmE on composition for the present systems.
In order to perform a numerical comparison of the estimation capability of the PFP theory, we calculated the standard percentage deviations (σ%) using the relation
(21)σ%=[∑{100(expt-theor.)expt}2(n-1)]1/2,
where n represents the number of experimental data points.
4.3. Bloomfield and Dewan Model
There are different expressions available in the literature to calculate η. Here, Bloomfield and Dewan [45] model have been applied to compare calculated Δlnη values using experimental data for each binary mixture by the following relation:
(22)Δlnη=lnη-(x1lnη1+x2lnη2).
Bloomfield and Dewan [45] developed the expression from the combination of the theories of free volumes and absolute reaction rate
(23)Δlnη=f(υ)-ΔGRRT,
where f(υ) is the characteristic function of the free volume defined by
(24)f(υ)=1/(υ-1)x1(υ1-1)-x2υ2-1
and ΔGR is the residual energy of mixing, calculated with the following expression [46]:
(25)ΔGR=ΔGE+RT{x1ln(x1Φ1)+x2ln(x2Φ2)},
where Φ1 and Φ2 are segment fractions defined by
(26)Φ2=1-Φ1=x2[x2+x1(r1/r2)],
where r1 and r2 are in the ratio of respective molar core volumes V1* and V2*.
The excess energy can be obtained from the statistical theory of liquid mixtures proposed by Flory and coworkers [31, 32] and is given by
(27)ΔGE=x1P1*V1*[1(υ1)-1(υ)]+3T~1ln{(υ11/3-1)(υ1/3-1)}+x2P2*V2*[1(υ2)-1(υ)]+[3T~2ln{(υ21/3-1)(υ1/3-1)}]+(x1θ2V1*χ12)υ,
where the various symbols used have their usual meanings.
Using the χ12 values from fitting values of VmE and the values of the parameters for the pure liquid components, we have calculated ΔGR and f(υ) and finally Δlnη, according to the Bloomfield and Dewan relationship, which is compared with the experimental data. Figure 3 shows that the good agreement between the estimated and experimental curves occurs for given binary systems.
Plot of Δlnη against mole fraction of octan-2-ol with (◊) n-octane; (□) n-decane; (▵) n-dodecane, and n-tetradecane (○) at T = 298.15 K. The corresponding dotted (- - -) curves have been derived from Bloomfield and Dewan model.
4.4. Jouyban and Acree Model
Hasan et al. [6, 47, 48] has used this model for various binary system proposed a model for correlating the density and viscosity of liquid mixtures at various temperatures. The proposed equation is
(28)lnym,T=f1lny1,T+f2lny2,T+f1f2∑[Aj(f1-f2)jT],
where ym,T, y1,T, and y2,T is density or viscosity of the mixture and solvents 1 and 2 at temperature T, respectively, f1 and f2 are the volume fractions of solvents in case of density and mole fraction in case of viscosity, and Aj are the model constants.
The correlating ability of the Jouyban-Acree model was tested by calculating the average percentage deviation (APD) between the experimental and calculated density and viscosity as
(29)APD=(100N)∑[(|yexptl-ycalcd|)yexptl],
where N is the number of data points in each set. The optimum numbers of constants Aj, in each case, were determined from the examination of the average percentage deviation value.
The constants Aj calculated from the least square analysis are presented in Table 7 along with the average percentage deviation (APD). The proposed model provides reasonably accurate calculations for the density, viscosity, and ultrasonic velocity of binary liquid mixtures at 298.15 K, and the model could be used in data modeling.
Parameters of Jouyban-Acree model and average percentage deviation for densities, viscosities, and ultrasonic velocities at of binary mixtures at T = 298.15 K.
System
A0
A1
A2
APD
Density
n-Octane + octan-2-ol
3.7446
−1.8027
0.616
0.0293
n-Decane + octan-2-ol
−9.0926
1.0266
0.6448
0.1183
n-Dodecane + octan-2-ol
−11.8504
0.0401
0.4771
0.0816
n-Tetradecane + octan-2-ol
−12.6833
3.5123
−4.3389
0.0559
Viscosity
n-Octane + octan-2-ol
−7.1847
−6.9259
3.3215
0.1423
n-Decane + octan-2-ol
−4.2926
11.2857
−2.7614
0.7668
n-Dodecane + octan-2-ol
−7.5692
−10.9008
4.4397
0.1031
n-Tetradecane + octan-2-ol
−6.4032
−26.9478
−7.7829
0.4296
5. Conclusions
The present paper is a continuing effort towards the understanding of the mixing behavior of binary liquid mixtures comprising of n-alkanes + octan-2-ol. Excess molar volumes and deviations in viscosity were calculated and fitted to the Redlich-Kister equation to test the quality of the experimental values. The negative deviations in viscosity support the main factor of breaking of the self-associated alcohols and weak interactions between unlike molecules. Several empirical and semiempirical relations have been used to represent the dependence of viscosity on concentration of components in binary liquid mixtures. An attempt has been made to check the suitability of empirical and semiempirical relations for experimental viscosities data of n-alkanes + octan-2-ol fits by taking into account the number of empirical adjustment coefficients. Bloomfield and Dewan model and McAllister’s three-body interaction model are suitable for the above binary system. PFP theory is also quite successful in predicting the trend of the dependence of VmE on composition for the present systems.
Acknowledgments
Authors are thankful to Prof B. R. Arbad, Dr. B. A. M. university for their valuable suggestions and discussion. Authors are also thankful to Principal Dr. R. S. Agarwal, J. E. S. College, Jalna for the facilities provided.
LalK.TripathiN.DubeyG. P.Densities, viscosities, and refractive indices of binary liquid mixtures of hexane, decane, hexadecane, and squalane with benzene at 298.15 K20004559619642-s2.0-003425798110.1021/je000103xEstrada-BaltazarJ. A.JuanF. J.GustavoA.Experimental liquid viscosities of decane and octane + decane from 298.15 K to 373.15 K and up to 25 MPa19984334414462-s2.0-0032074418ZhangY.VenableR. M.PastorR. W.Molecular dynamics simulations of neat alkanes: the viscosity dependence of rotational relaxation19961007265226602-s2.0-0000376273AminabhaviT. M.AralaguppiM. I.GopalakrishnaB.KhinnavarR. S.Densities, shear viscosities, refractive indices, and speeds of sound of bis(2-methoxyethyl) ether with hexane, heptane, octane, and 2,2,4-trimethylpentane in the temperature interval 298.15-318.15 K19943935225282-s2.0-0028466730Gómez-DízD.MejutoJ. C.NavazaJ. M.Rodríguez-ÁlvarezA.Viscosities, densities, surface tensions, and refractive indexes of 2,2,4-trimethylpentane + cyclohexane + decane ternary liquid systems at 298.15 K20024748728752-s2.0-003666317010.1021/je010288nHasanM.ShirudeD. F.HirayA. P.SawantA. B.KadamU. B.Densities, viscosities and ultrasonic velocities of binary mixtures of methylbenzene with hexan-2-ol, heptan-2-ol and octan-2-ol at T = 298.15 and 308.15 K20072521-288952-s2.0-3384694032110.1016/j.fluid.2007.01.001IloukhaniH.SamieyB.MoghaddasiM. A.Speeds of sound, isentropic compressibilities, viscosities and excess molar volumes of binary mixtures of methylcyclohexane + 2-alkanols or ethanol at T = 298.15 K20063821902002-s2.0-3014444069610.1016/j.jct.2005.04.019PalA.BhardwajR. K.Excess molar volumes and viscosities of binary mixtures of diethylene glycol dibutyl ether with chloroalkanes at 298.15 K20024147067112-s2.0-0036109103MahajanA. R.MirganeS. R.DeshmukhS. B.Volumetric, Viscometric and Ultrasonic studies of some amino acids in aqueous 0. 2M. LiCIO4. 3H20 solutions at 298.15 K200742345352MahajanA. R.MirganeS. R.DeshmukhS. B.Ultrasonic studies of some amino acids in aqueous 0.02M. LiCI04. 3H20 solutions at 298.15 K200742373378TreszczanowiczA. J.BensonG. C.Excess volumes for n-alkanols + n-alkanes II. Binary mixtures of n-pentanol, n-hexanol, n-octanol, and n-decanol + n-heptane197810109679742-s2.0-0005857565AliA.NainA. K.SharmaV. K.AhmadS.Molecular interactions in binary mixtures of tetrahydrofuran with alkanols (C6,C8,c10): an ultrasonic and volumetric study20044296666732-s2.0-10444229560IloukhaniH.Rezaei-SametiM.Basiri-ParsaJ.Excess molar volumes and dynamic viscosities for binary mixtures of toluene + n-alkanes (C5−C10) at T = 298.15 K—Comparison with Prigogine-Flory-Patterson theory20063889759822-s2.0-3374596643710.1016/j.jct.2005.10.011MehraR.PancholiM.Temperature-dependent studies of thermo-acoustic parameters in hexane + 1-dodecanol and application of various theories of sound speed20068032532632-s2.0-33645919279YangC.XuW.MaP.Thermodynamic properties of binary mixtures of p-xylene with cyclohexane, heptane, octane, and N-methyl-2-pyrrolidone at several temperatures2004496179418012-s2.0-974425880810.1021/je049776wMeyerR.MeyerM.MetzgerJ.PenelouxA.Thermodynamic and physicochemical properties of binary solvent197168406412RedlichO.KisterA. T.Thermodynamics of nonelectrolyte solutions. Algebraic representation of thermodynamic properties and the classification of solutions194840345348IrvingJ. B.Viscosity of binary liquid mixtures, a survey of mixture equations1977630East Kilbride, UKNational Eng LabIrvingJ. B.The effectiveness of mixture equations1977631East Kilbride, UKNational Eng LabGrunbergL.NissanA. H.Vaporisation, viscosity, cohesion and structure of the liquids1949164799800TamuraM.KurataM.Viscosity of binary mixture of liquids1952253237HindR. K.McLaughlinE.UbbelohdeA. R.Structure and viscosity of liquids camphor+pyrene mixtures1960563283302-s2.0-0000779830KattiP. K.ChaudhriM. M.Viscosities of binary mixtures of benzyl acetate with dioxane, aniline, and m-cresol1964934424432-s2.0-33947483159HericE. L.On the viscosity of ternary mixtures196611166682-s2.0-0002590369HericE. L.BrewerJ. G.Viscosity of some binary liquid nonelectrolyte mixtures19671245745832-s2.0-33947333918McAllisterR. A.The viscosities of lquid mixtures19606427431GlasstoneS.LaidlerK. J.EyringH.1941New York, NY, USAMcGraw-HillNigamR. K.MahlB. S.Molecular interaction in binary liquid mixtures of dimethylslphoxide with chlroehanes & chlroehenes197191255FortR. J.MooreW. R.Viscosities of binary liquid mixtures196662111211192-s2.0-11544309057RamamoortyK.VaradachariP. S.Study of some binary liquid mixtures197311238FloryP. J.Statistical thermodynamics of liquid mixtures1965879183318382-s2.0-33745918488AbeA.FloryP. J.The thermodynamic properties of mixtures of small, nonpolar molecules1965879183818462-s2.0-33745928017PrigogineI.1957Amsterdam, The NetherlandsNorth-HollandPattersonD.DelmasG.Corresponding states theories and liquid models197049981052-s2.0-000294603810.1039/DF9704900098TancrèdeP.BothorelP.De St. RomainP.PattersonD.Interactions in alkane systems by depolarized Rayleigh scattering and calorimetry. Part 1.—Orientational order and condensation effects in n-hexadecane + hexane and nonane isomers197773115282-s2.0-3704909800510.1039/F29777300015De St. RomainP.VanH. T.PattersonD.Effects of molecular flexibility and shape on the excess enthalpies and heat capacities of alkane systems197975170017072-s2.0-000072003010.1039/F19797501700RodriguezA. T.PattersonD.Excess thermodynamic functions of n-alkane mixtures. Prediction and interpretation through the corresponding states principle19827835015232-s2.0-001334817210.1039/F29827800501AminabhaviT. M.BanerjeeK.BalundgiR. H.Thermodynamic interactions in binary mixtures of 1-chloronaphthalene and monocyclic aromatics19993887687772-s2.0-0033505258RiddickJ. A.BungerW. B.SakanoT. K.19862New York, NY, USAJohn Wiley & SonsKrishnaiahA.NaiduP. R.Excess thermodynamic properties of binary liquid mixtures of 1,2-dichloroethane with normal alkanes19802521351372-s2.0-0343819579AicartE.TardajosG.Díaz PeñaM.Isothermal compressibility of cyclohexane + n-hexane, cyclohexane + n-heptane, cyclohexane + n-octane, and cyclohexane + n-nonane19802521401452-s2.0-0000981065PandeyJ. D.PantN.Surface tension of ternary polymeric solution198210412329933022-s2.0-33845553859Al-JimazA. S.Al-KandaryJ. A.Abdul-LatifA.-H. M.Acoustical and excess properties of {chlorobenzene + 1-hexanol, or 1-heptanol, or 1-octanol, or 1-nonanol, or 1-decanol} at (298.15, 303.15, 308.15, and 313.15) K20075212062142-s2.0-3384707418710.1021/je060353zAicartE.TardajosG.Diaz PeñaM.Isothermal compressibility of cyclohexane-n-decane, cyclohexane-n-dodecane, and cyclohexane-n-tetradecane198126122262-s2.0-0002566769BloomfieldV. A.DewanR. K.Viscosity of liquid mixtures19717520311331192-s2.0-0344784584JouybanA.KhoubnasabjafariM.Vaez-GharamalekiZ.FekariZ.AcreeW. E.Jr.Calculation of the viscosity of binary liquids at various temperatures using Jouyban-Acree model20055355195232-s2.0-1964439734110.1248/cpb.53.519JouybanA.Fathi-AzarbayjaniA.KhoubnasabjafariM.AcreeW. E.Jr.Mathematical representation of the density of liquid mixtures at various temperatures using Jouyban-Acree model2005448155315602-s2.0-28244490902KrishnaiahA.NaiduP. R.Excess thermodynamic properties of binary liquid mixtures of 1,2-dichloroethane with normal alkanes19802521351372-s2.0-0343819579