Enkog Castillo1990
Enkog Castillo1990
Enkog Castillo1990
North-Holland
R. C A S T I L L O and J.V. O R O Z C O
Instituto de Ffsica, Universidad Nacional Autonoma de M~xico, P.O. Box 20-364,
01000 Mexico, DF, Mexico
A method for calculating the thermal conductivity of dense fluids and their mixtures is
described. This method relies on the results of the revised Enskog theory for hard-sphere fluid
mixtures, and the use of the temperature and density dependent diameters given by Mansoori
and Canfield, and Rasaiah and Stell. In addition, a correction in order to take into account
the internal degrees of freedom contribution was introduced. A comparison with the
conformal solution model, and with accurate experimental data was made. A very good
agreement was found.
1. Introduction
In recent years, there has been an increasing demand from industry for
transport properties information of fluid mixtures, and in particular for thermal
conductivity (A). However, it is clear that the acquisition of reliable data for
the great variety of mixtures and thermodynamic states can never be achieved
in a satisfactory way by direct measurements only. There are two ways out,
mainly in engineering calculations, to fulfill this gap of information. The first
one uses empirical correlation schemes. These usually are limited to narrow
ranges of temperatures and pressures and often mainly designed for pure
fluids. A review of these methods are given in Reid et al. [1].
The second way involves the use of predictive methods more or less related
to kinetic theory. In this case, the prediction of transport properties is more
reliable if the m e t h o d adopted is founded on a realistic physical model of the
fluid mixture, and if it is combined with a suitable theory. Even when neither
models nor the theories are exact, it is possible to give reasonably accurate
results depending little on measurements or empirical prescriptions.
In predictive methods there are two groups of procedures for obtaining A's in
the dense regime. The first one is mainly based on the conformal solution
model (CSM), which is related to the corresponding-states principle. The best
suited method for the prediction of A on non-polar fluid mixtures over the
entire range of PVT states was given by Ely and Hanley [2]. They realized this
method using a state equation, and a A surface for a reference fluid (methane).
Although some additional corrections are required in order to include the
internal degrees of freedom contribution, and deviations from the correspond-
ing states model, this method has good predictive features (for details see refs.
[2] and [3]).
The other group of procedures mainly used for predicting A's for dense fluids
is related to the kinetic theory of hard spheres developed by Enskog [4]. Some
examples can be found elsewhere [5-8]. Here, one procedure follows the
so-called modified Enskog theory [5] (MET), an ad hoc modification of the
Enskog theory to include attractive intermolecular forces. Until quite recently
[9] MET was generalized to the binary case, although its ability to predict A's in
binary mixtures has not been reported. Another procedure relies on the
revised Enskog theory [10] (RET) combined with a prescription to obtain
effective state-dependent hard-sphere diameters in terms of the parameters
associated with the potential chosen to model actual systems. We call this
approach the effective-diameter revised Enskog theory (EDRET). We will
show in this paper that this procedure can be generalized to the binary case, in
the same way as we did for other transport properties [11-14]. Here, the basic
assumptions are three: (1) particles in the fluid mixture can be modeled by the
Lennard-Jones potential. (2) The behavior of A in this regime is mainly
determined by the repulsive part of the interaction potential. (3) The hard-
sphere expression for A coming from the RET can give confident results if
some prescription to obtain state-dependent effective diameters is used in
order to reflect the somewhat soft repulsive part of the model potential.
The main purpose of this paper is to show the usefulness of E D R E T to
predict A's of one-component fluids and their mixtures, and to present some
comparison between our approach and the CSM. The paper is organized as
follows: in section 2, we present the details of our method. In section 3, the
numerical results are presented and compared with the CSM, and experimental
data.
2. Theory
RET, first derived by van Beijeren and Ernst [10]. Here, the two hard-sphere
radial distribution functions are the same functionals of the number densities as
the radial distribution functions of a binary mixture in nonuniform equilibrium.
The R E T equations can be solved by the use of the C h a p m a n - E n s k o g solution
method. The molecular fluxes and the transport coefficients for dense hard-
sphere binary mixtures, up to the Navier-Stokes level, can be directly obtained
on the basis of the procedure used in ref. [15].
Thermal conductivity in a binary mixture, 3`, is defined through the following
expression [16]:
a T _ p T D T olx a (1)
Jq = - 3 , -~r Or
where Jq is the heat flux, T is the temperature, p is the mass density, D T is the
thermal diffusion coefficient a n d / z a is the chemical potential of component 1.
Here, we must have a comment. Diffusion and heat conduction in a fluid
mixture become inextricably coupled. In a uniform mixture, as time progresses
in the course of an experimental determination, concentration gradients appear
in the mixture and they themselves provide driving forces for opposing fluxes.
Eventually, in a steady state thermal diffusive flux is exactly balanced by this
back diffusive flux, so that there is no net molecular flow. Hence, to stress this
fact experimentalists define two A coefficients. The first one is 3`0, this is the
thermal conductivity of a uniform mixture (no concentration gradients), and
the second one is 3`=, the thermal conductivity corresponding to a fluid in the
steady state (no mass fluxes). Since many accurate experiments actually
measure 3`=, as in the transient hot-wire technique, we decided to calculate this
quantity. The formal definition of 3`= for a binary mixture is given by [17]
OT
Jq = - A = 0---7 (Ji = O), (2)
+ kk n'n'
i4x:'
508 R. Castillo and J.V. Orozco / Thermal conductivity of dense fluids
and
th th
dj = d j ~rrl°gT" (5)
crPLJ(r)=4eiI(~2)12--(-~)6] , (6)
where Ei is the well depth and o-0i is the minimal separation of two molecules of
species i such that ~LJ(tr0i) = 0, and separately the effective diameters in terms
of o'0i, e i, and of the thermodynamic state of the system are determined. For
the cross interaction, we assumed it to be that of a hard-sphere mixture with
effective diameters, i.e.,
(BH) [19], by Weeks, Chandler and Andersen (WCA) [20], and by Verlet and
Weiss (VW) [21]. From the schemes coming from a variational theory, we have
those given by Mansoori-Canfield [22] and Rasaiah-Stell [23] (MC/RS). All
these schemes give a hard-sphere diameter depending on both temperature and
density, except that given by BH that is temperature dependent only.
In order to obtain numerical estimates of A's for actual fluids using EDRET,
a correction to include the internal degrees of freedom contribution is needed.
But, until now, there is no formal theory to handle transport properties of
polyatomic dense fluids, although several attempts have appeared [24, 25] with
a very limited success. In this paper we will follow an heuristic approach
followed by some authors [5, 26], and suggested by a previous work of Mason
and Monchick [27] on the basis of the Wang Chang-Uhlenbeck-De Boer
theory [28]. This theory was developed by these authors to deal with the
transport properties of polyatomic gases in the dilute regime. Mason and
Monchick showed that the A can effectively be separated into two parts: a part
dealing with the transfer of thermal energy by the translational motion of the
molecules, and a part dealing with the transfer of energy to the internal
degrees of freedom of the molecules. Hence, we will assume that A of a dense
fluid can be split into a part due to the energy transfer by molecular motion
and collisional transfer (A'), given by EDRET, and a part due to the energy
transfer associated with the internal degrees of freedom of the molecules (A").
Thus,
A = (8)
In addition, we assume that A" for the one-component fluid can be represented
by the first order approximation formula given by Mason-Monchick for
quasi-elastic collisions [5]:
tP
for the dilute hard-sphere gas, C v IS the molar heat capacity at constant volume
for the internal degrees of freedom, and M is the molecular weight.
For the case of the binary mixture the evaluation of h" has not been explored
from the point of view of kinetic theory. Hence, a mixing rule is the only way
out. This rule should keep a close analogy with eq. (9) in order to have a fair
density dependence. Hence, we modified the empirical mixing rule given by Li
[29]:.
)tmix(P, T, x,) =
(v xixjAi/. X12(P, T, xi), (10)
- /'=1
510 R. Castillo and J.V. Orozco / Thermal conductivity of dense fluids
where
0
where B0i is the dilute gas viscosity, Cp, i is the ideal gas molar heat capacity, R
is the gas constant, M i is the molecular weight of component i, and f~nt has a
constant value of 1.32.
0
Cp,i values for different systems were obtained through a temperature
expansion up to 6th order, except for the noble gases where the internal
contribution was not considered (actually we used the expansions that appear
in the TRAPP computer program developed by Ely and Hanley [3]).
6O
T r = 2.5 °o Tr=l 1
4O °'~° .~;;;
elll. _llit~l e~
÷
20~Oot ° ". .".' ." . . . . . . 11
e~lol 0 .D°ooo. •
C<3
o....:;~,.,,s o;°
a o
-20 °
~ o
o
-40 o
60
Tr=2 T r = 09
4O
°;~,~"~.
oo ,1o 4
",~a,aR.*..~aoo"OO°°
6<b 0 ooo o-
o
r20 o
-40 o
a
60 a+ +
40
Tr= 1.5 Tr = 0 7
4;
2 0 ~, * e ° °
oo. ~o
0 - I : i i l l ~ l l t .l . .i. .J J ~oooo no •= a
d.<b
-20 .
o
-40
-60 J~, , ~, ,~L ,, . . . . . r ........ , t l l , l l l l l l l l l ~ l l I i t l l l , t l ~
o 1 2 pr3,O 1 2
Fig. 1. Percent deviation for predicted A's from experimental data, BA, versus the reduced density,
Pr, for six different reduced temperatures, and several prescriptions to obtain the state-dependent
hard-sphere diameters. (E3) MC/RS; (+) BH; (O) WCA; (A) VW.
119.8 K), along the phase diagram. As we can see, our procedure increasingly
overestimates A as we progressively go to lower reduced densities. At lower
densities (below Pr ~ 0.7) the other schemes to obtain the hard-sphere diam-
eter ( B H , W C A and V W ) are better to fit the experimental data. At higher
temperatures they become more accurate. This is not surprising because these
schemes were built to fit the low density equilibrium properties, whereas the
M C / R S scheme was directed toward the denser regimes. The important point
to emphasize is that since the crossover occurs at a reduced density of about
0.7, our recommendation is that one uses the other schemes from reduced
densities below 0.7 and the M C / R S scheme from 0.7 onwards. The general
conclusion is: the M C / R S prescription gives better predictions in a wide range
of fluid reduced densities ( P r = 0 . 7 - - 2 . 7 at T r > 1), and at liquid densities
512 R. Castillo and J.V. Orozco / Thermal conductivity of dense fluids
(T, < 1). This is in close agreement with our previous work in other transport
properties [8, 11-14].
When a comparison is made between our A calculations, and accurate
experimental data for nitrogen, methane, ethane, and ethylene, we found a
better agreement in the denser or liquid region of the phase diagram than in
the moderate or low density one, confirming the conclusions given for the case
of argon.
Figs. 2 and 3 show A for liquid nitrogen and ethylene respectively, along the
saturation line. Here, we compare our procedure, the CSM, and experimental
data [32, 33]. In a general way, we can say that our calculations are closer to
the experimental data than CSM. On the other hand, in these figures we can
see that the internal degress contribution as implemented above is less im-
portant as density increases. This agrees with previous results where the role of
this contribution has been discussed [34].
Figs. 4 and 5 show A for mixtures. Fig. 4 shows a nitrogen/oxygen mixture
(X~2 = 0.78569, Xo2 = 0.21431, p = 0.6 g/cm3). Here, A calculations are com-
pared with CSM, and with experimental data for air [35], at several tempera-
tures. As we can see in this figure, EDRET without the internal degrees
contribution gives a better fitting to the experimental data than CSM. If we
include the internal degrees contribution to EDRET, we have a better fitting to
the experimental data at lower temperatures, but at high temperatures our
procedure overestimates A around 15%.
O
180 I
"~ 170 '
E 160 o
E
15o a o
140 ~.,. o
>- 150
I-- 120
> 110
~- 100
9O
g3
z 8O
0
(D TO
6O
_J
<~ 50
oz 4O
uJ 30 I I J I I I J I I I I
-r 70 75 80 85 90 95 100 105 HO 115 120 125
TEMPERATURE (OK)
Fig. 2. Thermal conductivity, A, of liquid nitrogen versus temperature, at the saturation line. (El)
experimental; dashed line, CSM; solid line, EDRET without the internal degrees correction; (O)
EDRET with the internal degrees correction.
R. Castillo and J.V. Orozco / Thermal conductivity of dense fluids 513
,,--
o 280
-%
E 260 ~ ^ T
24O
E 22O - ~.~
OA
+ •
>.- zoo +
>
160 + •
o
D
C]
Z
0
0 - z; :, :,
_J
<[
80 ++
60 + +
n-
LIJ 40 i L L i l
I 105 130 160 190 220 250
I--
TEMPERATURE (OK)
Fig. 3. Thermal conductivity, A, of liquid ethylene versus temperature, at the saturation line. (A)
experimental; solid line, CSM; (+) EDRET without the internal degrees correction; (O) E D R E T
with the internal degrees correction.
92
90
88-
86
84
82
80
78
76
74.
72
70
68
66
4- 4- 4- 4-
4- 4- 4- 4-
62 +
60 4- 4- 4- 4- + 4-
58 I [ 1 I I I 1 I [ I l I I I I
120 125 130 1,35 140 145 150 15,5 160 165 170 17,5 180 185 190 195 200
TEMPERATURE
Fig. 4. Thermal conductivity of air at 0.6 g cm -3 versus temperature. ([~) experimental; (+) CSM;
lower solid line, E D R E T without the internal degrees correction; upper solid line, E D R E T with
the internal degrees correction.
514 R. Castillo and J.V. Orozco / Thermal conductivity o f dense fluids
It EXP
X -- e~et
--- cS
0,12
0,02 I
~0 25 Pmolar
/
T=203.697K
eo
0.12
O.OZ i
0 Z5 Pmolar
Acknowledgements J,
References
[1] R.C. Reid, J.M. Prausnitz and T.K. Sherwood, The Properties of Gases and Liquids, 3rd ed.
(McGraw-Hill, New York, 1977).
[2] J.F. Ely and H.J.M. Hanley, Ind. Eng. Chem. Fundam. 22 (1983) 90.
[3] F. Ely and H.J.M. Hanley, A Computer Program for the Prediction of Viscosity and Thermal
Conductivity in Hydrocarbon Mixtures, Nat. Bur. Stand (U.S.) Tech. Note 1039 (April 1981).
[4] S. Chapman and T.G. Cowling, The Mathematical Theory of Nonuniform Gases, 3rd ed.
(Cambridge Univ. Press, Cambridge, 1970).
[5] H.J.M. Hanley, R.D. McCarty and E.G.D. Cohen, Physica 60 (1972) 322.
[6] J.C.G. Calado, J.M.N.A. Fareleira, U.V. Mardolcar and C.A. Nieto de Castro, Int. J.
Thermophys. 9 (1988) 351.
[7] J.H. Dymond, Q. Rev. Chem. Soc. 3 (1985) 317.
[8] R. Castillo and J.W. Orozco, Int. J. Thermophys., to be published.
[9] J.M. Kincaid, S. Perez and E.G.D. Cohen, Phys. Rev. A 38 (1988) 3628.
[10] H. van Beijeren and M.H. Ernst, Physica 68 (1973) 437.
[11] R. Castillo, M. Lopez de Haro and E. Martina, Int. J. Thermophys. 7 (1986) 851.
[12] R. Castillo, E. Martina and M. Lopez de Haro, KINAM 7 (1986) 61.
[13] R. Castillo and S. Castafieda, Int. J. Thermophys. 9 (1988) 383.
[14] R. Castillo and S. Castafieda, J. Non-Equil. Thermodyn. 14 (1989) 69.
[15] M. Lopez de Haro, E.G.D. Cohen and J.M. Kincaid, J. Chem. Phys. 78 (1983) 2746.
[16] S.R. De Groot and P. Mazur, Non-Equilibrium Thermodynamics (Dover Publications, New
York, 1984).
[17] H.E. Khalifa, J. Kestin and W.A. Wakeham, Physica A 97 (1979) 273.
[18] G.A. Mansoori, N.F. Carnahan, K.E. Starling and T.W. Leland, J. Chem. Phys. 54 (1971)
1523.
[19] J.A. Barker and D. Henderson, J. Chem. Phys. 47 (1967) 4714.
[20] J.D. Weeks, D. Chandler and H.C. Andersen, J. Chem. Phys. 54 (1971) 5237.
[21] L. Verlet and J.J. Weiss, Phys. Rev. A 5 (1972) 939.
[22] G.A. Mansoori and F.B. Canfield, J. Chem. Phys. 51 (1969) 4958.
[23] J. Rasaiah and G. Stell, Mol. Phys. 18 (1970) 249.
[24] S.A. Rice and A.R. Allnatt, J. Chem. Phys. 34 (1961) 2144.
[25] S. Jagannathan, J.S. Dahler and W. Sung, J. Chem. Phys. 83 (1985) 1808.
[26] H.J.M. Hanley and E.G.D. Cohen, Physica A 83 (1976) 215.
[27] E.A. Mason and L. Monchick, J. Chem. Phys. 36 (1962) 1622.
[28] C.S. Wang Chang, G.E. Uhlenbeck and J. de Boer, in: Studies in Statistical Mechanics, vol.
2, part C, J. De Boer and G.E. Uhlenbeck, eds. (North-Holland, Amsterdam, 1964).
[29] C.C. Li, AIChE J. 22 (1976) 927.
[30] J. Karkheck and G. Stell, J. Chem. Phys. 75 (1981) 1475.
[31] B.A. Younglove and H.J.M. Hanley, J. Phys. Chem. Ref. Data 15 (1986) 1323.
[32] K. Stephen, R. Krauss and A. Laesecke, J. Phys. Chem. Ref. Data 16 (1987) 993.
[33] P.M. Holland, B.E. Eaton and H.J.M. Hanley, J. Phys. Chem. Ref. Data 12 (1983) 917.
516 R. Castillo and J.V. Orozco / Thermal conductivity of dense fluids
[34] H.J.M. Hanley, Proc. 7th Symp. on Thermophysical Properties, A. Cezairliyan, ed. (Am.
Soc. Mech. Eng., New York, 1977).
[35] K. Stephan and A. Laesecke, J. Phys. Chem. Ref. Data 14 (1985) 227.
[36] H.M. Roder and D.G. Friend, Int. J. Thermophys. 6 (1985) 607.
[37] L. Mistura, J. Chem. Phys. 62 (1975) 4571.