1.

2 ESTIMACIN DE LA VISCOSIDAD
Es frecuente encontrar en la bibliografa las propiedades fsicas de las sustancias puras o de
sus mezclas en las condiciones deseadas, en cuyo caso deber recurrirse a su estimacin
con la mxima precisin posible.
La estimacin de las propiedades fsicas en general y concretamente de las propiedades de
transporte que nos ocupan, puede basarse en la teora, en determinadas correlaciones de
valores experimentales o en una y otras.
En cualquier caso, debe extremarse la prudencia en no aplicar una teora fuera de los limites
para los que se haya establecido y en no extrapolar las correlaciones deducidas
experimentalmente mas all de los intervalos de las variables para los que fueron
formuladas.
Los mejores mtodos estimativos implican casi siempre ecuaciones tericas y correlaciones
empricas para el clculo de las constantes de las mismas que la teora no ha conseguido
llegar a establecer. Se comprende que a medida que disminuye la aportacin emprica en
cualquier correlacin aumentar el grado de validez de la misma.
A continuacin pasaremos revista a los mtodos ms importantes de estimacin de la
viscosidad de las sustancias puras y de sus mezclas, tanto en estado gaseoso como lquido.
1.2.1 VISCOSIDAD
i.Gases puros a pequeas presiones
a) Mtodo Terico
La teora cintica rigurosa conduce a la ecuacin (1.2.1), que de acuerdo con la bibliografa
puede expresarse tambin as:

1.2.1

mkT
5
2 2 2,2 *
16

MRT
5

16 1/2 N 2

1.2.2

26, 69

MT

1/2

1/2


expresndose la viscosidad
temperatura en T en K.

en micropoises, el dimetro de colisin

en

y la

Los dimetros y las integrales de colisin,

y
, son funciones de las fuerzas
intermoleculares y, por consiguiente, para su evaluacin se requiere establecer hiptesis
sobre la naturaleza de las mismas, que lgicamente variaran segn que las molculas del
gas sean o no polares.
-

Gases no polares

En este caso las fuerzas intermoleculares no dependen de la orientacin relativa de las

molculas. En estos gases se ha comprobado experimentalmente que para grandes
distancias entre las molculas dichas fuerzas son atractivas e inversamente proporcionadas
a la distancia y que cuando las molculas se aproximan a distancias muy pequeas, de
0

algunos
, surgen fuerzas repulsivas entre ellas que varan mucho mas rpidamente con la
distancia.
EP
Esta conducta se explica si la energa potencial ( ) del sistema formado por dos
molculas varia con la distancia entre ellas, r, del modo que se indica en la figura 1.2.

Figura 1.2 Variacin de la energa potencial con la distancia en un sistema formado

por dos molculas.

Se considera nula dicha energa potencial cuando la distancia entre las molculas es
infinita. A medida que las molculas se aproximan, la energa potencial disminuye debido
al trabajo realizado por las fuerzas de atraccin, llegando un momento en que empiezan a
manifestarse las fuerzas de repulsin, por lo que la indicada energa potencial alcanzara un

rm

mnimo
al alcanzarse una cierta distancia
entre las molculas a la que se
equilibran unas y otras. Para seguir aproximado las molculas deber realizarse un trabajo
sobre el sistema para vencer las fuerzas repulsivas, con lo que su energa potencial

( )
aumentar. La distancia entre ambas molculas
, para la que la energa potencial vuelve
a alinearse, recibe el nombre de dimetro de colisin: distancia hasta la que se
aproximaran dos molculas inicialmente en reposo y a una distancia infinita si se pusiesen
en movimiento a causa de las fuerzas de atraccin.
La forma exacta de la funcin representada en la figura 1.2:

1.2.3

Ep f r

no se conoce, pero los resultados experimentales indican que una funcin emprica muy
satisfactoria para las molculas no polares es el llamado potencial de Lennard-Jones:
12 6
E p 4
r
1.2.4
r

EP min
con los dos parmetros: energa potencial mnima

( )
, y dimetro de colisin

( / r )12

,a

( / r ) 6

los que nos acabamos de referir, y con los trminos

y
fuerzas repulsivas y atractivas entre las molculas, respectivamente.

, que representan las

Puede comprobarse que:

EP
-

La funcin
(1.2.4) solo se anula para valores de la distancia intermolecular
r r
y
.

EP
-

La funcin

rm 2

1/6

(1.2.4), derivada e igualada a cero, conduce a un valor

EP min

EP
, que sustituido en la propia funcin

puesto que su segunda derivada es el citado valor

rm 2

conduce a
1/6

, resulta positiva.

Ambas circunstancias demuestran la idoneidad del potencial de Lennard-Jones (1.2.4) para

representar la funcin (1.2.3), aunque se han propuesto muchos otros con tal fin.

Numerosos investigadores han determinado la integral de colisin

en funcin del
potencial de Lennard-Jones. Recientemente, Neufeld y col. propusieron la siguiente
ecuacin emprica apropiada para su aplicacin en ordenador:

1.2.5

representando

T*

1,16145 0,52487 2,16178

*
*
T *0,14874 e0,77327T
e 2,43787T

la temperatura adimensional:
kT
T
*

1.2.6

0,3 T * 100
que segn sus autores, es aplicable a un intervalo
0,064%.

, con error medio del

Mediante datos experimentales, de viscosidad a distintas temperaturas, y las ecuaciones

/k
(1.2.2), (1.2.5) y (1.2.6), algunos investigadores han determinado los parmetros
y
para muchas sustancias. En la tabla 1.1 se resumen los correspondientes a algunas de ellas.
Si se dispone de tales parmetros, mediante las ecuaciones (1.2.2) y (1.2.5) puede
calcularse la viscosidad de los gases puros no polares a pequea presin y cualquier
temperatura con error inferior al 1%. De no disponerse a tales parmetros, aunque los
mismos pueden calcularse aproximadamente a partir de las constantes crticas y los factores

acntricos, resulta ms convenientemente el mtodo de los estados correspondiente que se

trata mas adelante.
Gases polares
EP
Para molculas polares debe utilizarse unafuncin potencial
mas adecuada que la de
Lennard-Jones. En la bibliografa se propone la sugerida por Stockmayer, que modifica la
ltimamente citada mediante un trmino adicional representativo de las interacciones
dipolo-dipolo entre las molculas.

Tabla 1.1 Parmetros y / k

para diversas sustancias

Compuesto
(
Ar
He
Kr
Ne
Xe
Aire
BCl3
BF3
B(OCH3)3
Br2
CCl4
CF4
CHCl3
CH2Cl2
CH3Br
CH3Cl
CH3OH
CH4
CO
COS
CO2
CS2
C2H2

Argn
Helio
Kriptn
Nen
Xenn
Aire
Cloruro de boro
Fluoruro de boro
Borato de metilo
Bromo
Tetracloruro de carbono
Tetrafluoruro de carbono
Cloroformo
Cloruro de metileno
Bromuro de metilo
Cloruro de metilo
Metanol
Metano
Monxido de carbono
Sulfuro de carbonilo
Dixido de carbono
Disulfuro de carbono
Acetileno

3,542
2,551
3,655
2,820
4,048
3,711
5,127
4,198
5,503
4,296
5,947
4,662
5,380
4,898
4,118
4,182
3,626
3,758
3,690
4,130
3,941
4,483
4,033

/ k
)

(K)
93,3
10,22
178,9
32,8
231,0
78,6
337,7
186,3
396,7
507,9
322,7
134,0
340,2
356,3
449,2
350
481,8
148,6
91,7
336,0
195,2
467
231,8

C2H4
C2H6
C2H5Cl
C2H5OH
CH3OCH3
CH2CHCH3
CH3CCH
C3H6
C3H8
n-C3H7OH
CH3COCH3
CH3COOCH3
n-C4H10
iso- C4H10
C2H5OC2H5
CH3COOC2H5
n-C5H12
C(CH3)4
C6H6
C6H12
n-C6H14
Cl2
F2
HBr
HCN
HCl
HF
HI
H2
H2O
H2O2
H2S
Hg
HgBr2
HgCl2
HgI2
I2
NH3

Etileno
Etano
Cloruro de Etilo
Etanol
ter dimetlico
Propileno
Metil acetileno
Ciclopropano
Propano
Alcohol n-proplico
Acetona
Acetato de metilo
n-Butano
IsoButano
ter dietlico
Acetato de etilo
n-Pentano
2;2-Dimetilpropano
Benceno
Ciclohexano
n-Hexano
Cloro
Flor
cido bromhdrico
cido cianhdrico
cido clorhdrico
cido fluorhdrico
cido iodhdrico
Hidrgeno
Agua
Perxido de Hidrogeno
Acido sulfhdrico
Mercurio
Bromuro de mercurio
Cloruro de mercurio
Ioduro de mercurio
Yodo
Amoniaco

4,163
4,443
4,898
4,530
4,307
4,678
4,761
4,807
5,118
4,549
4,600
4,936
4,687
5,278
5,678
5,205
5,784
6,464
5,349
6,182
5,949
4,217
3,357
3,353
3,360
3,339
3,148
4,211
2,827
2,641
4,196
3,623
2,969
5,080
4,550
5,625
5,160
2,900

224,7
215,7
300
362,2
396,0
298,9
251,8
248,9
237,1
576,7
560,2
469,8
531,4
330,1
313,8
521,3
341,1
193,4
412,3
297,1
399,3
316,0
112,6
449
569,1
344,7
330
288,7
59,7
809,1
289,3
301,1
750
686,9
750,0
695,6
472,2
558,3

NO
N2
N2O
O2
SO2
SiF4
SnBr4
UF6

Oxido nitroso
Nitrgeno
Oxido ntrico
Oxigeno
Dixido de azufre
Tetrafluoruro de silicio
Bromuro estnnico
Hexafluoruro de uranio

3,492
3,798
3,828
3,467
4,112
4,880
6,388
5,967

116,7
71,4
232,4
106,7
335,4
171,9
563,7
236,8

Brokaw, tras considerar detenidamente la utilizacin del potencial de Stockmayer,sugiere

para el caso de molculas polares la siguiente expresin de la integral de colisin

0, 2 2

Stockmayer

Lennard-Jones

1.2.7
T*
T*

En la que (Lennard-Jones) y
vienen dadas por las ecuaciones (1.2.5) y (1.2.6) y
un parmetro adimensional con expresin:
p2

1.2.8
2 3

p
Siendo

el momento dipolar y

es

los parmetros indicados anteriormente.

A partir de los datos experimentales viscosidad-temperatura y las ecuaciones (1.2.2) y

/k
(1.2.5) a (1.2.8), se han determinado los parmetros
,
y de las molculas polares
mas representativas, resumindose los mismos en la tabla 1.2, pudiendo encontrarse otros
valores en la bibliografa.
Tabla 1.2 Parmetros de la funcin potencial de Stockmayer

p
(debyes)

( A )

Agua

1,85

Amoniaco

1,47

Compuestos

/ k

2,52

(K)
775

1,0

3,15

358

0,4

cido clorhdrico

1,08

3,36

328

0,34

cido bromhdrico

0,80

3,41

417

0,14

cido iodhdrico

0,42

4,13

313

0,029

Dixido de azufre

1,63

4,04

347

0,42

Acido sulfhdrico

0,92

3,49

343

0,21

Cloruro de nitrosilo

1,83

3,53

690

0,4

Cloroformo

1,013

5,31

355

0,07

Cloruro de metileno

1,57

4,52

483

0,2

Cloruro de metilo

1,37

3,94

414

0,5

Bromuro de metilo

1,80

4,25

382

0,4

Cloruro de etilo

2,03

4,45

423

0,4

Metanol

1,70

3,69

417

0,5

Etanol

1,69

4,31

431

0,3

Alcohol n-proplico

1,69

4,71

495

0,2

Alcohol i-proplico

1,69

4,64

518

0,2

ter dimetlico

1,30

4,21

432

0,19

ter dietlico

1,15

5,49

362

0,08

Acetona

1,20

4,50

549

0,11

Acetato de metilo

1,72

5,04

418

0,2

Acetato de etilo

1,78

5,24

499

0,16

Nitrometano

2,15

4,16

290

2,3

Si se dispone de los parmetros acabados de indicar, mediante las ecuaciones (1.2.2),

(1.2.5) y (1.2.7) puede estimarse la viscosidad de los gases puros polares a pequea presin
y cualquier temperatura, con error inferior al 1,5%. De no ser as, aunque pueden calcularse
dichos parmetros a partir del momento dipolar, de la temperatura de ebullicin normal y
del volumen molar de la sustancia en estado lquido, resulta ms aconsejable recurrir al
mtodo de los estados correspondientes que se expone a continuacin.
b) Mtodo de los estados correspondientes
De la ecuacin (1.2.2) se deduce la siguiente proporcionalidad entre la viscosidad y las
variables que la afectan:

M 1/2 RT

1.2.9
2

Evidentemente, el dimetro de colisin

1/2

1
3

puede considerarse proporcional a la potencia

c
del volumen crtico de las molculas,

vc1/3 RTc

1.2.10
Tc

1/3

/ pc1/3

RTc

2/3

pc2/3

pc

siendo y
la temperatura y presin criticas, sobre las que existe amplia informacin en
la bibliografa.
Por tanto, de las proporcionalidades anteriores se deduce, para la viscosidad en las
condiciones criticas:

1.2.11

M 1/2 pc2/3

RTc

1/6

Tr
Podr definirse ahora una viscosidad reducida, en funcin de la temperatura reducida (
de la siguiente forma:

),

1/2 2/3
f Tr
c M pc / RTc 1/6

1.2.12

M 1/2Tc1/6 pc2/3 f Tr

Se han propuesto diversas correlaciones de este tipo, pero aqu solo indicamos las ms
modernas y tiles.

Gases no polares

Thodos y col. ha propuesto la siguiente correlacin:

1.2.13

M 1/2Tc1/6 pc2/3 4, 61 Tr0,618 2, 04 e 0,449Tr 1, 94 e 4,058Tr 0,1

Tc

pc

expresando la viscosidad en micropoises,

en Kelvins y

en atmsferas.

Reichenberg sugiri la siguiente correlacin:

1.2.14

a*Tr

1 0,36Tr Tr 1

1/6

siendo para la mayora de los compuestos inorgnicos, de acuerdo a la proporcionalidad

(1.2.11):

1.2.15

a*

3,5M 1/2 pc2/3

Tc1/6

Tc
Expresando la viscosidad en micropoises,

pc
en Kelvins y

en atmsferas.

Si se trata de compuestos orgnicos, Reichenberg propone el siguiente valor del parmetro

a*

1.2.16

a*

M 1/2Tc
niCi
i

ni
representando

Ci
el nmero de grupos atmicos del tipo i y

grupo que se resumen en la tabla 1.3.

las contribuciones de cada

De no poder estimarse la viscosidad por el mtodo terico, debe utilizarse cualquiera de los
representados por las relaciones (1.2.13) o (1.2.14), que implican un error entre 1 y 3%

Tabla 1.3 Valores de las contribuciones de los grupos Ci para la estimacin del
parmetro a* [ecuacin (1.2.16)], micropoises.
Grupo
-CH3
>CH2 (no en anillo)
>CH- (no en anillo)
>C< (no en anillo)
=CH2
=CH- (no en anillo)

Contribucin
Ci
9,04
6,47
2,67
-1,53
7,68
5,53

>C=

1,78

CH
C- (no en anillo)
>CH2 (en anillo)
>CH- (en anillo)
>C< (en anillo)
=CH- (en anillo)
>C= (en anillo)

7,41
5,24
6,91
1,16
0,23
5,90
3,59

Grupo
-Br
-OH (alcoholes)
>O (no en anillo)
>C=O (no en anillo)
-CHO (aldehdos)
-COOH (acidos)
-COO- (esteres) o
HCOO- (formiatos)
-NH2
>NH (no en anillo)
=N- (en anillo)
-CN
>S (en anillo)
-F
-Cl

Contribucin
Ci
12,83
7,96
3,59
12,02
14,02
18,65
13,41
9,71
3,68
4,97
18,13
8,86
4,46
10,06

Gases polares

Aunque Thodos y col. tambin proponen correlaciones del tipo de la (1.2.13) para gases
puros polares a pequea presin, de no poder utilizarse el mtodo terico para la estimacin
de la viscosidad de los mismos, se recomienda el mtodo de Reichenberg, ecuacion
(1.2.14), que implica un error inferior al 4%.
Ejemplo 1.1
Mediante los mtodos tericos y de los estados correspondientes, estmese la viscosidad del
SO2

135 p

a 40C y 1 atm. En tales condiciones el valor experimental es

Datos:

Tc 430,8K pc 77,8 atm Vc 122, 0 cm3 / mol

Condiciones criticas:
,
y
Solucin:
a) Mtodoterico.Teora cintica rigurosa. Ecuacin (1.2.2)

26, 69

TM
2 v

De la tabla 1.2 se deduce:

_

p 1, 63 debyes; 4, 04 A ; / k 347 K; 0, 42
La temperatura adimensional T*, ecuacin (1.2.6) ser:
_

k T 313
T

0,902

347
*

que sustituida en la ecuacin de Neufeld y col (1.2.5) conduce a:

r L.J .

1,16145

0,902

0,14874

0,52487
2,16178
2,43787 0,902 1, 680
0,7732 0,902
e
e

Como el SO2 es un gas polar, de la ecuacin (1.2.7) se obtiene:

r S .

0, 2 0, 42
1,680
0,902

1, 719

y sustituyendo estos valores en la ecuacin [1] resulta:

26, 69

313 64
2
4, 04 1, 719

Comparando este valor con el experimental:

134, 6 p

Error

134, 6 135
100 0,3%
135

b) Mtodo de los estados correspondientes

En este caso, aunque se trata de un gas polar, aplicaremos la ecuacin de Reichenberg
(1.2.14):

a*Tr

1 0,36Tr Tr 1

1/6

siendo:

a*

3,5M 1/2 pc2/3

Tc1/6

Sustituyendo valores en las ecuaciones [5] y [4], se obtiene:

a
*

3,5 64

77,8
1/6
430,8
1/2

2/3

185, 7

185.7 0, 727
1 0,36 0, 727 0, 727 1

1/6

136, 7

Comparando este valor con el experimental:

Error

136, 7 135
100 1, 26%
135

ii. Mezclas gaseosas a pequeas presiones

a) Mtodo terico
La aplicacin de la teora cintica rigurosa a la determinacin de las viscosidades de
mezclas gaseosas multicomponentes a pequeas presiones conduce a complejas ecuaciones
que pueden simplificarse si se desprecian los efectos de segundo orden. De esta manera, se
llega a una solucin en serie de la forma:

yi i

m
i 1

1.2.17

y
j 1

i ij

yj

yi
Representando

las fracciones molares de los c componentes,

la viscosidad de

ij
cada componente i puro y

un parmetro calculable, como se indica a continuacin.

ij
Mtodo de Wilke de estimacin de
Wilke, utilizando el mtodo de la teora cintica de Sutherland, propuso la expresin:

1/4
M j

M i
1/2
M i
1
M j

1 i
j

ij

1.2.18

1/2

e intercambiando subndices se deduce que:

1.2.19

ji

j Mi
ij
i M j

Las ecuaciones (1.2.17) y (1.2.18) han sido comprobadas reiteradamente tanto con mezclas
M i M j

i j

de componentes no polares como polares. Salvo en los casos en que

y
, como sucede si uno de los componentes es el hidrogeno o el helio, los errores fueron
siempre inferiores al 2-4%. Para mezclas con los dos elementos ltimamente citados,

ij
Reichenberg propuso una correlacin para

ij
Mtodo de Brokaw de estimacin de

bastante ms compleja.

Tras un amplio estudio sobre las viscosidades de las mezclas gaseosas, Brokaw propuso el

ij
siguiente mtodo de estimacin del parmetro

ji i

1.2.20
j

1/2

Aij Sij

Aij

Donde

es funcin exclusiva de la razn de pesos moleculares:

Aij mij M

1.2.21

M ij M ij0,45

1/2
ij

1 M m

0,45
ij

2 1 M ij

1 mij

1/2
ij

siendo:

1.2.22

mij

0,25

1
1 M ij 1 M ij

1.2.23

M ij

Mi
Mj

expresin que puede simplificarse a:

M i
Aij
M

1.2.24
j

0,37

si

0, 4

Mi
1,33
Mj

con error inferior al 1%, en el indicado intervalo de la razn de pesos moleculares.

Sij

El trmino

de la ecuacin (1.2.20) es la unidad si los compuestos i y j son apolares, es

decir, si
y , ecuacin (1.2.8) o tabla 1.2, son inferiores a 0,1. De no ser as, se calcula
mediante la expresin.

Sij S ji

1.2.25

T *
1 i *
T
j

i2
1 Ti
4

i j
4

j2
1 T
4

1/2

*
j

, as como los de

/k

con los que se puede

T j*

Ti *

calcular

i
La tabla 1.2 proporciona los valores de

1/2

1/2

mediante la ecuacin (1.26).

Los errores de este mtodo son similares a los del mtodo de Wilke, pero es preferible a
este ltimo en el caso de que la mezcla contenga compuestos polares.
b) Mtodo de los estados correspondientes
Pueden utilizarse los mtodos propuestos para gases puros evaluando convenientemente el
peso molecular y las constantes crticas correspondientes a las mezclas. Ninguno de estos
mtodos parece ser tan preciso como los tericos anteriormente explicados.
Ejemplo 1.2
SO2
Estmese la viscosidad de unos gases de tostacin dela siguiente composicin molar:
O2
(8%),

N2
(13%),

(79%), que se encuentran a 1 atm y 20C.

Datos:
Viscosidad de los componentes puros a 1 atm y 20C (p)
SO2 125 O2 203 N 2 175
;
;

Solucin:
Como se trata de una mezcla gaseosa a presin atmosfrica, se aplicara el mtodo terico,
ecuacin (1.2.17):

yi i

m
i 1

y
j 1

i ij

ij
evalundose

mediante la ecuacin de Wilke (1.2.18):

1/4
M j

1/2
M
1 i
M j

1 i

ij

1/2

SO2 O2 N 2
Representado con los subndices 1, 2 y 3 los componentes
,
,
, respectivamente
mediante las Ec. [1] y [2] se construye la siguiente tabla:
c

M i /M j

i / j

ij

yi

ij y i

1
1
1
2
2
2
3
3
3

1
2
3
1
2
3
1
2
3

1
2
2,286
0,5
1
1,143
0,437
0,875
1

1
0,616
0,714
1,623
1
1,16
1,400
0,862
1

1
0,562
0,555
1,827
1
1,007
1,778
0,992
1

0,08
0,13
0,79
0,08
0,13
0,79
0,08
0,13
0,79

0,08
0,0073
0,438
0,146
0,13
0,795
0,142
0,129
0,79

Finalmente, la ecuacin [1], con los valores tabulados, conduce a:

0, 08 125 0,13 203 0, 79 175

0,591

1, 071

1, 061

172 p

ij y i
j=1

0,591

1,071

1,061

iii. Gases puros a presiones elevadas

Slo en determinadas regiones de presin y temperatura, la viscosidad de los gases es
funcin de la presin. Habitualmente, tales variaciones, no son significativas a temperaturas
reducidas muy elevadas o a presiones reducidas muy bajas.
Childs y Hanley han establecido un criterio par saber si debe corregirse o no el efecto de la
presin sobre la viscosidad de los gases. En la figura 1.3 se esquematiza el mismo en
funcin a presin y temperatura reducidas. Segn que las condiciones de los gases
conduzcan a puntos situados por encima o por debajo de la recta de la figura, habr que
efectuar o no dichas correcciones; si aquellas condujesen a puntos sobre la misma recta y
no efectuasen tales correcciones, el error solo seria del 1%.

pr

p
pc

0.5

0
0

Figura 1.3 Intervalos de temperaturas y presiones reducidas para la posible

correccin por los efectos de presin.

La figura 1.4slo aproximada, refleja bien las variaciones de la viscosidad de los gases con
la temperatura y la presin. Mediante ella puede conseguirse una rpida estimacin
aproximada de la viscosidad de los gases a presiones elevadas. Para ellos debe calcularse la
viscosidad crtica mediante la siguiente ecuacin:
7, 7 M 1/2 pc2/3

1.2.26 c
Tc1/6

c
en la que se expresan:

Tc
en micropoises,

pc
en Kelvins y

en atmsferas.

La estimacin ms rigurosa de la viscosidad de los gases a elevadas presiones se basa en

correlaciones de las diferencias de viscosidad de los mismos a elevada y baja presin

0
(viscosidad residual,
), y la densidad reducida
Thodos y col. proponen las siguientes ecuaciones:

r / c Vc / V
. Concretamente,

Figura 1.4
Viscosidades reducidas
generalizadas

Tr

T
Tc

Gases no polares

1.2.27

Tc1/6
1
1/2 2/3
M pc

0,25

1, 023 0, 23364 r 0, 58533 r2 0, 40758 r3 0, 093324 r4

0,1 r 3

Gases polares
1/6
c
1/2 2/3
c

1.2.28
M
0

1, 656 r1,111

r 0,1
,

1/6
c
1/2 2/3
c

1.2.29
M
0

0, 0607 9, 045 r 0, 63

1,739

0,1 r 0,9
;

1.2.30

log 4 log

Tc1/6

M 1/2 pc2/3

T 1/6
log 4 log 0 1/2c 2/3
M pc
1.2.31

0, 6439 0,1005 r

0,9 r 2, 2
;

0, 6439 0,1005 r 4, 75 10 4 r3 10, 65

2, 2 r 2, 6

1.2.32

Tc1/6
90
r 2,8
M 1/2 pc2/3
;

1.2.33

Tc1/6
250
r 3
M 1/2 pc2/3
;

ecuaciones en las que las viscosidades

Tc
se expresan en micropoises.

enkelvins y

pc
en atmosferas, debindose recurrir a la bibliografa para encontrar los valores del

Vc
volumen crtico

siendo preferibles valores experimentales de densidades

Excepto a presiones muy elevadas, el error del mtodo es siempre inferior al 10%.

/ 0
Recientemente, Reichenberg propone una correlacin de la razn
algo compleja,
pero, al parecer, bastante precisa tanto para gases apolares como polares. Aunque todava
no muy probada, parece ser que el error de la misma es similar al de la correlacin anterior.

Ejemplo 1.3
Estmese la viscosidad del n-butano a 171C y las siguientes presiones: 6,8 atm; 34 atm y
68 atm.
Datos:
Coeficientes de comprensibilidad y viscosidad experimental del n-butano a 171C y
diferentes presiones:
Presin (atm)
1
6,8
34
68

Z
1,00
0,93
0,54
0,29

(p)
112
113
134
405

Tc 425, 2 K pc 37,5atm Vc 255cm3 / mol

Condiciones criticas:
;
y
Solucin
a) Mtodo grafico
En la tabla 1 se relacionan las condiciones reducidas para el n- butano a 171 C y diferentes
presiones. En todos los casos, del grfico de Childs y Hanley (Fig.1.3), se deduce que el gas
debe considerarse como denso. En la tabla 1 tambin se relacionan las viscosidades
reducidas generalizadas estimadas grficamente (Fig.1.4):
Tabla 1
Presin (atm)

pr= p / p c

T r =T /T c

(p)

Error (%)

6,8

0,181

1,044

0,45

108

-4,4

34

0,907

1,044

0,7

168

+25,4

68

1,813

1,044

2,2

527

+30,1

La viscosidad crtica se calcula por la ecuacin (1.2.26):

1
y en el caso del n-butano se tendr:

7, 7 M 1/2 pc2/3
Tc1/6

7, 7 58

37,5
1/6
425, 2
1/2

2/3

239, 6 p

y, por tanto:

r c 239, 6 p
obtenindose a cada presin los valores que se indican en la tabla 1. Se observa que los
errores son considerados as como que aumentan al hacerlo la presin.
b) Mtodo de Thodos y col (1.2.27)

Tc1/6
1
1/2 2/3
M pc

0,25

1, 023 0, 23364 r 0, 58533 r2 0, 40758 r3 0, 093324 r4

0,1 r 3
ecuacin aplicable para el siguiente intervalo

Vc
Vc
255 103
p
r

7, 004 10 3
c V ZRT Z 0, 082 444
Z
p
3
p

Sustituyendo valores en las ecuaciones [2] y [3] se obtienen los valores que se indican en la
tabla 2.
Tabla 2
P (atm)

r [3]

(p)

Error (%)

6,8

0,051

No es aplicable

34

0,441

147

+9,8

68

1,642

457

+13

Pudindose observar que los errores disminuyen apreciablemente con este mtodo de
estimacin.

iv. Mezclas gaseosas a presiones elevadas

Dean y Stiel proponen la siguiente correlacin para la estimacin de las viscosidades de
mezclas gaseosas a presiones elevadas:

1.2.34
m
en la que

m m0

1/6
Tcm
1,858

1, 08 exp1, 439 rm exp 1,111 rm

1/2 2/3

M m pcm

m0
y

son las viscosidades de la mezcla (micropoises) a elevada y baja presin

rm m / cm
respectivamente, y
la densidad seudorreducida de la mezcla, cociente de la
densidad de la mezcla con la seudocrtica. Los parmetros seudocrticos de la mezcla son
calculables de la forma:
c

1.2.35

Tcm yiTci Z cm yi Z ci Vcm yiVci

i 1

i 1

i 1

pcm

Z cm RTcm
Vcm

El error de estimacin, si se trata de mezclas de gases apolares de bajo peso molecular, es

menor al 10%. En el caso de mezclas de gases apolares de elevado peso molecular o en las
que existan uno o ms gases polares, se sigue utilizando el mtodo por no haber otro mejor,
pero lgicamente el error es superior.
Para mezclas gaseosas puede seguir utilizndose la figura 1.3,para decidir si ha de tenerse
en cuenta el efecto de la presin elevada sobre la viscosidad.
Ejemplo 1.4
Estmese la viscosidad de una mezcla gaseosa de metano (40%) y propano (60%) a 225C y
100 atm.
Datos:
Condiciones criticas:
T c (K)
Metano
Propano

190,6
368,8

pc (atm
)
45,4
41,9

V c (cm3/
mol)
99,0
203,0

Zc
0,288
0,281

Viscosidad de la mezcla gaseosa a 1 atm y 225C, = 142 p. Factores de

comprensibilidad a 225C y 100 atm: metano, 0,98; propano, 0,76.
Solucin
Segn la ecuacin de Dean y Stiel(1.2.34), se tiene:

m m0

Tc1/6
1,858

1, 08 exp1, 439 rm exp 1,111 rm

1/2 2/3

M m pcm

siendo:

m0 142 p
Tcm 0, 4 190, 6 0, 6 369,8 298,12 K
M m 0, 4 16 0, 6 44 32,80
Vcm 0, 4 99 0, 6 203 161, 40cm3 / mol
Z cm 0, 4 0, 288 0, 6 0, 281 0, 284

pcm

Z cm RTcm 0, 284 0, 082 298,12

43, 01atm
Vcm
0,1614

Vmet

V prop

0, 082 498 0,98

100

0, 082 498 0, 76
100

0, 400 l / mol

0,310 l / mol

Vm 0, 4 400 0, 6 310 346 cm3 / mol

 rm

Vcm
Vm

161, 40
0, 466
346

Sustituyendo valores en la ecuacin[1] se llega a:

m 142

2,585
1, 08 exp 0, 671 exp 0, 269
5, 727 12, 276
m 177p

v. Lquidos puros
As como las viscosidades de los gases diluidos pueden estimarse mediante tcnicas
basadas en una teora segura, como es la cintica rigurosa, en el caso de los lquidos carece
de una base terica comparable para tal estimacin.
Las viscosidades de los lquidos son mucho mayores que las de los gases y descienden
abruptamente con la temperatura. Como se ha indicado al principio de este captulo, de las
dos etapas de que constan los transportes de cantidad de movimiento y energa, la primera
de desplazamiento de las molculas, esta muy limitada en los lquidos. En stos, el
transporte realmente se debe a los campos de fuerzas entre molculas suficientemente
prximas como demuestra su elevada densidad.
Las teoras para la estimacin de las viscosidades de los lquidos los asimilan a los gases o
slidos. En cualquier caso, ninguna de ellas conduce a un mtodo terico de estimacin,
sino a tcnicas empricas basadas en la evaluacin de algunas constantes a partir de la
estructura o propiedades fsicas del lquido de que se trate.Otras veces, sin tomar como base
modelo de liquido alguno, los mtodos se basan enteramente en el estudio y correlacin de
los datos experimentales.
La viscosidad de los lquidos a temperaturas inferiores a la normal de ebullicin no esta
muy afectada por la presin para valores moderados de la misma; a presiones muy elevadas
la viscosidad puede aumentar considerablemente. Al parecer, este ltimo efecto crece con la
complejidad molecular, pero no se dispone todava de ningn mtodo seguro de estimacin
de la viscosidad de los lquidos a baja temperatura y elevadas presiones.
En cuanto a la influencia de la temperatura sobre la viscosidad de los lquidos entre las
normales de ebullicin y congelacin, aun habindose propuesto ltimamente numerosas
relaciones, sigue siendo la de Andrade la ms utilizada:

 1.2.36

L Ae B /T

aunque puede conducir a resultados errneos en la proximidad de los puntos de

congelacin, temperaturas a las que la viscosidad de lquidos crece bruscamente. Con dos
valores de la viscosidad a dos temperaturas puede calcularse las constantes A y B, y con
ellas ya resulta posible determinar la viscosidad del liquido puro de que se trate a cualquier
otra temperatura.
Si slo se dispone de un valor experimental de la viscosidad a una cierta temperatura, la
correlacin grafica de Lewis y Squires (Fig 1.5), permite su extrapolacin aproximada. Se
fija el dato conocido sobre la curva y sobre el eje de abscisas la diferencia de temperaturas
para la que se desea el nuevo valor de la viscosidad, que podr leerse en las ordenadas.
T 100C

100C

100C

100C

100C L

100000

10000

1000

100

10

0.1

L T
Figura 1.5 Correlacin grafica

de Lewis y Squires

a) Estimacin de la viscosidad de los lquidos a bajas temperaturas

De los muchos mtodos empricos propuestos para la estimacin de la viscosidad de los

Tr
lquidos a bajas temperaturas (
inferior a 0,75-0,80), en ausencia total de datos
experimentales, se explican a continuacin los cuatro que se juzgan mejores y ms
generales. Todos ellos para lquidos orgnicos y se basan en la contribucin de los distintos
grupos funcionales, y establecen la misma relacin entre viscosidad y temperatura que la
relacin de Andrade (1.2.36), interviniendo la densidad de lquidos en dos de ellos.
Al comparar los datos experimentales con los calculados mediante estos mtodos, se
aprecian errores variables que pueden llegar al 15%. Por ello, en ausencia de datos
experimentales, es aconsejable utilizar preferentemente el ltimo, ms moderno.
Mtodo de Thomas
Thomas propuso la siguiente correlacin:

1
L
log 8,569 1/2
1

L
1.2.37
r

Tr
3

siendo

la viscosidad del liquido en centipoises,

la densidad del liquido en g/cm ,

la temperatura reducida y un parmetro estructural calculado aditivamente a partir de los

valores de la tabla 1.4, teniendo en cuenta la contribucin de los distintos tomos y grupos
atmicos.
Tabla 1.4 Contribuciones estructurales para el clculo del parmetro
de Thomas (1.2.37).
C
-0,462
H
0,249
O
0,054
Cl
0,340
1
Cetonas y esteres.

Br
I
Doble enlace
C6H5
2
Nitrilos.

0,326
0,335
0,478
0,385

S
CO1
CN2

de la ecuacin
0,043
0,105
0,381
|

El mtodo es aconsejable para compuestos aromticos (excepto benceno), compuestos

monohalogenados, compuestos no saturados y n-parafinas de elevado peso molecular. No
debe utilizarse para alcoholes, cidos, naftenos, heterociclos, animas, aldehdos o
compuestos polihalogenos.
Mtodo de Morris

Morris propuso la correlacin:

1.2.38

1
L
J
1

log

L
siendo

la viscosidad del lquido en centipoises,

una constante funcin del tipo de

Tr
compuesto (tabla 1.4),
la expresin:

la temperatura reducida y J un parmetro estructural definido por

1.2.39

J 0, 0577 bi ni

i 1

bi

1/ 2

ni

donde , es la contribucin de cada grupo (tabla 1.4) y

de cada grupo en la molcula.

el numero de veces de aparicin

La imprecisin del mtodo es mayor en los primeros miembros de las series homlogas, en
presencia de cadenas ramificadas y a bajas temperaturas.

Mtodo de Orrick y Erbar

Estos autores proponen la siguiente correlacin:

1.2.40

ln

L
B
A
L M

siendo
la viscosidad del lquido en centipoises,
la densidad del liquido a 20C
3
(g/cm ), M el peso molecular, T la temperatura en K y A y B dos parmetros estructurales
(tabla 1.4).
El mtodo no es aplicable a compuestos con azufre o nitrgeno.

Tabla 1.5 Constante

bi
y contribuciones

de la ecuacin de Morris (1.2.38)

Hidrocarburos
Derivados
halogenados
Derivados
bencnicos

0,087
5

Derivados
bencnicos
halogenados

0,123

teres, cetonas,
aldehdos, acetatos

0,148

Alcoholes

0,081
9

Fenoles

0,089
5

cidos orgnicos

0,117

Varios

Grupo

bi

0,09
6
0,12
6
0,10
0

Grupo

bi

CH3, CH2, CH

0,0825

CH2: miembro anillo saturado

0,1707

Halgeno en CH3

0,0

CH3, CH2, CH: junto a un anillo

0,0520

Halgeno en CH2

0,0893

NO2: junto a un anillo

0,4170

Halgeno en CH

0,0667

NH3: junto a un anillo

0,7645

Halgeno en C

0,0

F, Cl: junto a un anillo

0,0

Br

0,2058

OH: alcoholes

2,0446

Cl

0,1470

COOH: cidos

0,8896

0,1344

C=O: cetonas

0,3217

0,1908

O=CO: acetatos

0,4369

Doble enlace

-0,0742

OH: fenoles

3,4420

C6H4, anillo bencnico

0,3558

O: teres

0,1090

H adicional en un anillo

0,1446

Tabla 1.6 Parmetros A y B de la ecuacin de Orrick y Erbar (1.2.40)

Grupo
tomos de C1

RCR

A
-(6,95 + 0,21n)

B
275+99n

-0,15

35

RCR

R
Doble enlace
Anillo de cinco miembros
Anillo de seis miembros
Anillo aromtico
Sustitucin orto
Sustitucin meta
Sustitucin para
Cloro
Bromo
Iodo
OH
COOO
C=O
COOH

-1,20

400

0,24
0,10
-0,45
0
-0,12
0,05
-0,01
-0,61
-1,25
-1,75
-3,00
-1,00
-0,38
-0,50
-0,90

-90
32
250
20
100
-34
-5
220
365
400
1600
420
140
350
770

n es el nmero de tomos de carbono, excluyendo los pertenecientes a los distintos grupos que se relacionan
en la tabla.

Mtodo de Van Venzel, Cardozo y Langenkamp

Estos autores estudiaron meticulosamente la influencia de la estructura del lquido sobre su
viscosidad y propusieron la correlacin:

1.2.41

1 1

T T0

log L B

T0

siendo

la viscosidad del liquido en centipoises, T la temperatura en K y B y

dos
N*
parmetros estructurales funcin de la longitud de cadena equivalente
. La complejidad
T0
del clculo de esta ultima, y por ello de los dos parmetros B y

, es tal que remitimos a la

citada bibliografa para el mismo, no sin indicar que en ella se encuentran ya calculados y
T0
tabulados los valores de B y

para numerosas sustancias.

El mtodo resulta bastante impreciso con los primeros trminos de las series homlogas.
Ejemplo 1.5
Estmese la viscosidad y densidad de la acetona a 0C, 30C y 60C. y dedzcase una
ecuacin emprica que permita tal estimacin a cualquier temperatura en el intervalo
indicado.
Datos:
Valores experimentales de la viscosidad y densidad de la acetona a diferentes temperaturas.
(cp)
0,389
0,292
0,226

T (C)
0
30
60

(gr/cm3)
0,80
0,78
0,75

Condiciones crticas:
TC 508,1

pC 464
K;

VC 209 cm3 /mol

atm;

Solucin
En todos los casos la temperatura reducida es menor que 0,7, por lo que la acetona puede
considerarse como liquido a baja temperatura.
a) Ecuacin de Thomas (1.2.37):

1

log 8,569 1/L2
1
L
1
r

La constante

se evala mediante la tabla 1.4:

2 0, 462 6 0, 249 0,105 0,675

y sustituyendo valores en la ecuacin (1) se obtiene:


L
508,1
log 8, 569 1/2
1
0, 675
L
2

expresin que relaciona la viscosidad con la temperatura. En la tabla 1 se relacionan

L
los valores de

obtenidos con la misma.

b) Ecuacin de Morris (1.2.38):

log

1
L

J
1

0, 096
En la tabla 1.5 se encuentra para las cetonas:
parmetro estructural J, mediante la ecuacin (1.2.39):
n

J 0, 0577 bi ni
4

i 1

; obtenindose el valor del

1/2

bi
Con los valores

de la tabla 1.5. As, resulta:

J 0, 0577 2 0, 0825 0,321

1/2

0, 737

Sustituyendo valores en la ecuacin [3] se llega a

log

L
508,1
0, 737
1
0, 096
T

ecuacin que relaciona la viscosidad con la temperatura. En la tabla 1 se relacionan

L
los valores de

obtenidos con esta ecuacin.

c) Ecuacin de Orrick y Erbar (1.2.40):

ln

L
B
A
L M

siendo A y B dos parmetros estructurales que se evalan mediante la tabla 1.6

A 6,95 2 0, 21 0,50 7,87

B 275 2 99 350 823
Sustituyendo valores en la ecuacin [6] se obtiene:

ln

L
823
7,87
L 58

ecuacin que relaciona la viscosidad con la temperatura. En la tabla 1 se presentan

L
los valores de

con ella obtenidos.

d) Ecuacin de Andrade (1.2.36):

L Ae B /T

ecuacin que permite correlacionar la viscosidad con la temperatura, pudindose

estimar A y B a partir de datos experimentales. De ella se deduce:

9 ln L ln A T

Ajustando mediante el mtodo de mnimos cuadrados los datos experimentales de la

viscosidad frente a la temperatura, de acuerdo con la ecuacin [9] se obtiene:

B 821, 64

A 1,9247 102
;
es decir:

10

L 1,9247 102 e821,47/T

L
En la tabla 1 se relacionan los valores de

evaluados con esta ecuacin.

Se puede observar que en este caso los errores ms bajos se obtienen con la
ecuacin de Thomas y los ms elevados con la de Orrick y Erbar. Por otro lado, la
mejor correlacin emprica es la ecuacin [10], basada en la de Andrade.

Tabla 1
T
(K)
273
303
333

L exp
(cp)
0,389
0,292
0,226

L [2]
(cp)
0,398
0,295
0,229

%
+2,3
+1,0
+1,3

L [5]
(cp)
0,414
0,303
0,234

%
+6,4
+3,8
+3,5

L [7]
(cp)
0,361
0,261
0,197

%
- 7,12
- 10,6
-12,8

L [9]
(cp)
0,390
0,290
0,227

%
+0,26
-0,68
+0,44

b) Estimacin de la viscosidad de los lquidos a elevadas temperaturas

Para lquidos saturados a elevadas temperaturas parece plausible resulten tiles las
correlaciones deducidas en el apartado 1.2.1.iii) por el mtodo de los estados
correspondientes para los gases a presiones elevadas, siempre que no se excedan los limites
de densidad para ellos establecidos.
Basndose en ellos, Letsou y Stiel proponen una correlacin que se represente grficamente
en la figura 1.6. Esta correlacin constituye el mejor mtodo disponible para la estimacin
de la viscosidad de los lquidos saturados a temperaturas reducidas superiores a 0,76; con
Tr 0,92
errores del 3% si no se sobrepasa

 L : Viscosidad (cp)
: Tc1/6 / M 1/ 2 pc2/3
Tc : Temperatura crtica (K)
pc : Presin crtica (atm)
M : Peso molecular
pis
1
pc

i : Factor acntrico log

p : Presin de saturacin a Tr 0, 7 i 0,3

s
i

Tr

T
Tc L 104

40

8
0.7

0.8

0.9

L Tr
Figura 1.6 Correlacin
vi. Mezclas lquidas.

de Letsou y Stiel

Tr 0, 75
A temperaturas bajas (
) la viscosidad de los lquidos es muy sensible a su
estructura, de la que dependen las posibles interacciones moleculares.
Se ha propuestos numerosas correlaciones para la estimacin de la viscosidad de las

m L
mezclas liquidas
generales:

, respondiendo casi todas ellas a dos posibles expresiones

c

1.2.42

f m L xi x j f ij
i 1 j 1

1.2.43

f m L xi f i L
i 1

ij
representando por

f L

la viscosidad de la mezcla binaria de los componentes i y j, y en las

ln L 1 / L
xi
que
puede ser: ,
,
, etc., y
puede ser la fraccin volumtrica,
msica o molar del componente i en la mezcla. Las correlaciones del tipo de la (1.2.43),
ms sencillas, no son recomendable si para determinadas composiciones la mezcla presenta
valores mximos o mnimos de la viscosidad.
Lobe revis las numerosas correlaciones de la bibliografa y propuso la expresin:

1.2.44

i 1

m
i i exp

i j

RT
j 1

, ji

i
siendo

la viscosidad cinemtica del componente i (centistokes),

la fraccin

j
volumtrica del componente j,
un parmetro de la viscosidad caracterstico del
componente j en la mezcla (cal/(mol)(K)), R la constante de los gases (1,987cal/(mol)(K)) y
T la temperatura(K).
La utilizacin de la ecuacin (1.2.44) para una mezcla de c componentes, requiere al menos

j
de c valores de la viscosidad de la misma que permitan calcular todos los valores de

necesarios. El autor indica la validez de su ecuacin para cualquier mezcla aunque presente
mximos o mnimos de viscosidad.
En el caso de mezclas binarias, si A es el componente de menos viscosidad y la viscosidad
de la mezcla no presentase mximos ni mnimos a determinadas composiciones, se cumple:

1, 7 ln B
RT
A

1.2.45

B

0, 27 ln B / A 1, 3ln B
A
1.2.46 RT

1/ 2

con lo que seria posible estimar la viscosidad de la mezcla a partir exclusivamente de las
viscosidades de los componentes puros.
La estimacin de las viscosidades de las mezclas liquidas a temperaturas elevadas (
Tr 0, 75
) ha sido poco investigada. Se sugiere el uso de la figura 1.6, utilizando para la
mezcla de que se trate de constantes criticas calculadas mediante las expresiones (1.2.35)
Ejemplo 1.6
Estmese la viscosidad a 20C de una disolucin acuosa de metanol al 80% en peso.
Datos:
Viscosidad y densidad de los componentes puros a 20C:

Agua
Metanol

(cp)
1
0,592

(g/cm3)
0,9922
0,7917

Solucin:
a) Se utilizar la ecuacin (1.2.43), suponiendo funcin logartmica:

1 ln m xA ln A xB ln B
es decir:

ln m 0,8ln 0,592 0, 2 ln1

m 0, 657 cp
b) Ecuacin de Lobe (1.244)
Aplicndola a una mezcla binaria se tiene:

2 m A A e / RT
B B

B B e A A / RT

siendo:
0,8
0, 7917
A
0,834
0, 2
0,8

0,9922 0, 7917

B 1 A 0,166

1.2.45
y de acuerdo con las ecuaciones

1.2.46
y

A
1 / 0,9922
1, 7 ln
0, 507
RT
0, 592 / 0, 7917

B
1 / 0,9922
1 / 0,9922

0, 27 ln
1,3ln
0, 703
RT
0,592 / 0, 7917
0,592 / 0, 7917

Sustituyendo valores en la ecuacin [2] se llega a:

1
0,592 0,703 0,166

0,507 0,834
e
0,166
0,810 cst

e
0,
7917
0,9922

m 0,834

considerando una densidad media para la mezcla liquida:

m 0, 2 0, 9922 0,8 0, 7917 0,832 g/cm3

se tiene:

m m m 0,810 0,832 0, 674 cp

Pudindose observar que los dos valores encontrados difieren en un 2,5%

BIBLIOGRAFA:
-

BIRD R. B., STEWART W. E., LIGHTFOOT E. N. (2006): Fenmenos de

Transporte 2da Edicin, Edit. Limusa Wiley., Mxico.

COSTA NOVELLA, E. (1984): Ingeniera Qumica. Tomo 2: Fenmenos de

Transporte, 1ra Edicin, Edit. Alhambra., Espaa.

FAHIEN R. W. (1983): Fundamentals of Transport Phenomena, 1ra Edicin, Edit.

McGraw-Hill. Inc., USA.