TF

Y. Martnez Ochoa, A. Bentez Hernndez facultad de ingenieras Qumica y Mecnica. universidad de Matanzas Camilo Cienfuegos (Cuba)

transporte de fluidos

Estimacin de la viscosidad de lquidos puros y mezclas de lquidos

Mtodos generalizados para bajas y altas temperaturas
Los mtodos de estimacin de la viscosidad de lquidos y mezclas de lquidos, al igual que la estimacin de otras propiedades de transporte, son muy diversos. En este trabajo se realiza un compendio de las ecuaciones recogidas en la bibliografa especializada para la estimacin de la viscosidad de lquidos y de mezclas de lquidos a bajas y altas temperaturas.
el diseo o anlisis de todo proceso qumico se rige por modelos matemticos, y estos a su vez se constituyen haciendo uso de diferentes leyes, como son las leyes de conservacin de la materia, energa y el movimiento, o las leyes cinticas que describen las reacciones qumicas, bioqumicas y las ecuaciones de velocidad de transporte; adems, dependen de las propiedades que tenga el fluido que estar sujeto al transporte, las que pueden ser fsicas, termodinmicas y de transporte. la viscosidad es la propiedad de los fluidos que indica la resistencia que ofrecen a su deformacin ante la accin de una tensin aplicada. depende de la naturaleza del fluido y de su estado termodinmico. segn Perry (1999), la viscosidad absoluta () se define como el esfuerzo absoluto en un punto dividido por el gradiente de velocidad en el punto, y la viscosidad cinemtica () es la relacin entre la viscosidad absoluta y la densidad a la misma temperatura y presin.

1.

Mtodos de estiMaCin
los mtodos de estimacin de la viscosidad de lquidos y mezclas de lquidos, al igual que la estimacin de otras propiedades de transporte, son muy diversos.
abril 2010 ingeniera QuMiCa 95

TRANSPORTE DE FLUIDOS

existen tres grupos de ecuaciones o mtodos que han sido comunicados por varios autores, que son los siguientes: - las ecuaciones basadas en la teora cintica molecular de los lquidos. - las ecuaciones basadas en el teorema de los estados correspondientes. - las ecuaciones basadas en procedimientos estructurales.

2.

donde: n: nmero real de tomos de carbono en la molcula ni: Parmetro de contribucin de cada grupo estructural. los valores de este parmetro para diferentes grupos funcionales se muestran en Tabla 1. el trmino ni se determina sumando las contrii buciones de cada grupo. si un grupo estructural aparece n veces, hay que considerar (n * n) correcciones. Una vez calculado n*, se determinan los valores de B y T0 de la siguiente manera:
4

lQuidos a bajas teMperaturas

la viscosidad de los lquidos puros a bajas temperaturas puede ser estimada por varios mtodos, que se describen a continuacin.

(4) Por otro lado:

5 6

2.1. Ecuaciones basadas en la teora cintica de los lquidos

la teora cintica molecular de los lquidos est mucho menos desarrollada que la de los gases, por lo que el conocimiento que se tiene referente a las propiedades de transporte de lquidos es fundamentalmente emprico. ecuacin de eyring [2, 3] la ecuacin de eyring es vlida para lquidos newtonianos y no se comporta bien para molculas largas o cercanas al punto crtico. Generalmente predice valores dentro de un margen de error de un 40%:
1

donde: n: nmero de avogadro (6,023. 1023 (g/mol)-1) h: Constante de Planck (6,624. 10-27 g cm2/ s) V: Volumen molar (cm3/g.mol-1) Tb: Temperatura normal de ebullicin (K) T: Temperatura (K)

donde: Bi: Parmetro de contribucin. los valores de este parmetro para cada grupo funcional se muestran en la Tabla 1. el trmino Bi se determina sumando las contrii buciones de cada grupo. aunque el grupo funcional aparezca ms de una vez en el compuesto, la contribucin Bi nada ms se tiene en cuenta una sola vez. en la Tabla 2 se muestran los valores para diferentes compuestos lquidos de los parmetros de n*, Bi y T0 que se utilizan en el clculo de la viscosidad por medio de la ecuacin de Van Velzen, Cardozo y langenkamp.

2.3. Efecto de la temperatura

la viscosidad de los lquidos es muy sensible a los cambios de temperatura. la bibliografa recomienda un gran nmero de procedimientos para estimar la viscosidad de lquidos a bajas temperaturas (Tr < 0,8), aunque todos presentan errores considerables. a continuacin se muestran algunas de las ecuaciones empricas que permiten predecir el efecto de la temperatura. ecuacin de Guzman andrade [1, 5]:
7

2.2. Ecuaciones basadas en procedimientos estructurales

Para estimar la viscosidad de lquidos orgnicos, se utiliza la ecuacin de Van Velzen, Cardozo y langerkamp. ecuacin de Van Velzen, Cardozo y langerkamp [1, 7]:
2

donde: : Viscosidad (Pa.s) T: Temperatura (K) B, T0: Parmetros relacionados con la estructura molecular del compuesto. Para determinar los valores de B y T0 se debe calcular la longitud de cadena equivalente (n*)
3

ecuacin de Vogel:

(8) donde: T: Temperatura a, B y C: Constantes, cuyos valores dependen de las unidades de y T y se calculan a partir de datos experimentales.
N 481

96 ingeniera QuMiCa

estiMaCin de la visCosidad de lQuidos puros y MezClas de lQuidos

Tabla 1

ParMetros de coNtribuciN de gruPos fuNcioNales Para el clculo de la viscosidad de lquidos utilizaNdo la ecuaciN de vaN velzeN, cardozo y laNgeNkaMP
Estructura o grupo funcional
n-alcanos isoalcanos Hidrocarburos saturados con dos grupos metil en posicin iso n-alquenos n-alcadienos isoalquenos isoalcadienos Hidrocarburos con un doble enlace y dos grupos metil en posicin iso Hidrocarburos con dos doble enlace y dos grupos metil en posicin iso Ciclopentanos Ciclohexanos alquilbencenos polifenoles Alcoholes primario secundario terciario dioles (correcciones) fenoles (correcciones) -oH en cadena lateral de un anillo aromtico (correcciones) Acidos

Ni
0 1,386-0,238 n 2,319-0,238 n -0,152-0,042 n -0,304-0,084 n 1,237-0,280 n 1,085-0,322 n 2,626-0,518n 2,474-0,560n 0,205+0,069n 3,971-0,172 n 1,48 6,517-0,311 n 0,60 3,055-0,161 n -5,340+0,815 n

Bi
0 15,51 15,51 -44,94+5,410 n* -44,94+5,410 n* -36,01+5,410 n* -36,01+5,410 n* -36,01+5,410 n* -36,01+5,410 n* -45,96+2,224 n* -339,67+23,135 n* -272,85+25,041 n* -272,85+25,041 n* -140,04+13,869 n* -140,04+13,869 n* -188,40+9,558 n*

Notas

Cada grupo CH3, adicional en posicin iso incrementa n en 1,389-0,238 n Cada grupo CH3, adicional en posicin iso incrementa n en 1,389-0,238 n n<16; no recomendado para n= 5 6 n16 n<17; no recomendado para n= 6 7 n17 n<16; no recomendado para n= 6 7 (a, e, f) n16 (a, e, f) (a)

10,606-0,276 n 11,200-0,605 n 3,055-0,161 n ver notas 16,17- n -0,16

-589,44+70,519 n* 497,58 928,83 557,77 213,68 213,68

(b) (b) (b) para n, hacer uso de las contribuciones del alcohol aadindole n 2,50 (a, c, d)

6,795+0,365 n 10,71 ver notas

-249,12+22,449 n* -249,12+22,449 n* -249,12+22,449 n*

isocidos

n<11; no recomendado para n=12 n11 b se calcula como para el cido de cadena lineal. n se calcula como para el cido de cadena lineal, pero se reduce en 0,24 por cada grupo metilo en posicin iso

acidos con anillos aromticos (correccin) Esteres esteres con anillos aromticos (correccin)

4,81 4,337-0,230 n -1,174+0,376 n

-188,40+9,558 n* -149,13+18,695 n* -140,04+13,869 n* si los radicales hidrocarburos tienen configuracin iso, ver (e) adicionarlo a los valores de n y b calculados como steres si los radicales hidrocarburos tienen configuracin iso, ver (e) adicionarlo a los valores de n y b calculados como cetonas si los radicales hidrocarburos tienen configuracin iso, ver (e) el valor de n no es una correccin al valor calculado de manera normal, mientras que b si que lo es, ver (c) si los radicales hidrocarburos tienen configuracin iso, ver (e) la correccin hay que aadirla al calculo como amina (e) si los radicales hidrocarburos tienen configuracin iso, ver (e) si los radicales hidrocarburos tienen configuracin iso, ver (e)

Cetonas Cetonas con anillos aromticos (correccin) Eteres eteres aromticos

3,265-0,122 n 2,70 0,298+0,209 n 11,5 - n

-117,21+15,781 n* -760,65+50,478 n* -9,39+2,848 n* -140,04+13,869 n*

Aminas primarias amina primaria en cadena lateral de compuesto aromtico (correccin) secundaria terciaria

3,581+0,325 n -0,16 1,390+0,461 n 3,27

25,39+8,744 n* 0 25,39+8,744 n* 25,39+8,744 n*

abril 2010

ingeniera QuMiCa 97

TRANSPORTE DE FLUIDOS

Tabla 1 (conT.)
Estructura o grupo funcional
primaria con el grupo nH2 en anillo aromtico

Ni
15,04 - n

Bi
0

Notas
el valor de n no es una correccin al valor calculado para la amina. para el clculo de b se usa el valor de la amina primaria

secundaria o terciaria con, al menos, un anillo aromtico unido al nitrgeno Compuestos nitro 1-nitro 2-nitro 3- nitro 4- nitro; 5- nitro nitroaromticos nitrilo isometil nitrilo nitrilo aromtico dinitrilo Compuestos halogenados fluorados Clorados bromados iodados Correcciones para configuraciones especiales C(Cl)x -CCl-CCl-C(br)x-Cbr-CbrCf3 en alcoholes Cf3 en otros compuestos Aldehdos aldehdos con anillos aromticos (correccin) Anhdridos anhdridos con anillos aromticos (correccin) amidas amidas con anillos aromticos (correccin)

(f)

(f)

7,812-0,236 n 5,84 5,84 5,84 7,812-0,236 n 4,039-0,0103 n -0,7228+0,1755 n 2,321-0,2357 n 10,452-1,1276 n

-213,14+18,330 n* -213,14+18,330 n* -338,01+25,086 n* -338,01+25,086 n* -213,14+18,330 n* -241,66+27,937 n* 286,26-31,009 n* -26,063-11,516 n* 3599,9-199,96 n*

Hay que tener en cuenta la contribucin del grupo alqueno

para la correccin aromtica, ver (f)

1,43 3,21 4,39 5,76

5,75 -17,03 -101,97+5,954 n -85,32

(e,f) (e,f) (e,f)

1,91-1,459 x 0,96 0,50 1,60 -3,93 -3,93 3,38 2,70 7,97-0,50 n 2,70 13,12+1,49 n 2,70

-26,38 0 81,34-86,850 x -57,73 341,68 25,55 146,45-25,11 n* -760,65+50,478 n* -33,50 -760,65+50,478 n* 524,63-20,72 n* -760,65+50,478 n*

(a): Para ms de una sustitucin en anillos aromticos se requieren correcciones adicionales N = 0,51 b = - 571,94 (con -oH) orto b = 54,84 (sin -oH) N = 0,11 b = 27,25 meta N = - 0,04 b = - 17,57 para (b): Para alcoholes, si hay un grupo metilo en posicin iso respecto a un alcohol, hay que aumentar N en 0,24 y b en 94,23 (c): si el compuesto tiene un grupo oH o NH2 en anillo aromtico, o se trata de un ter aromtico, utilizar el valor de N obtenido, pero deben despreciarse los otros sustituyentes del anillo (halgenos, cH3, No2, etc.). en cambio, estos grupos s que se tendrn en cuenta al calcular b. (d): Para los alcoholes aromticos y compuestos con oH en cadena lateral, se deben incluir la contribucin del grupo alcohol. Por ejemplo: Para el o-clorofenol b = b (alcohol primario) + b (cloruro) + b (fenol) + b (correccin orto (a)) con N* = 16,17 (c) b = [-589,44 + 70,519 (16,17)] + (-17,03) + 213,68 + (-571,94) = 175,56 Para el 2 fenil etanol N = N (alcohol primario) + N (correccin) N=8 N* = N + N = 8+ 8,24 = 16,24 b = b (alcohol primario) + b (correccin) b = [-589,44 + 70,519 (16,24)] + 213,68 = 769,47 (e): Para steres, alquilbencenos, hidrocarburos halogenados y cetonas, hay que disminuir N en 0,24 y aumentar b en 8,93 por cada grupo metilo en posicin iso. Para teres y aminas hay que disminuir N en 0,50 y aumentar b en 8,93 por cada uno de estos grupos. (f): Para alquilbencenos, nitrobencenos, bencenos halogenados y aminas secundarias o terciarias con uno o varios anillos aromticos unidos al nitrgeno, hay que aadir la siguiente correccin por cada anillo aromtico. si N<16 aumentar N en 0,60 aumentar N en 3,055 0,161 N si N 16 Para cualquier valor de N aumentar b en 140,04 +13,869 N*

98 ingeniera QuMiCa

N 481

estiMaCin de la visCosidad de lQuidos puros y MezClas de lQuidos

3.

lQuidos a altas teMperaturas

Para estimar la viscosidad de lquidos a altas temperaturas se utilizan las siguientes ecuaciones descritas a continuacin.

la viscosidad de los lquidos es muy sensible a los cambios de TemperaTura

Tabla 2

3.1. Ecuaciones basadas en el teorema de los estados correspondientes

ecuacin de letsou y stiel esta ecuacin es vlida para estimar la viscosidad de lquidos a temperaturas reducidas entre 0,76 y 0, 98:
9

valores Para difereNtes coMPuestos lquidos de los ParMetros de N*, bi y t0 que se utilizaN eN el clculo de la viscosidad Por Medio de la ecuaciN de vaN velzeN, cardozo y laNgeNkaMP .
Compuesto
n-hexano 2-metilbutano 2,3-dimetilbutano

los trminos (.)0 y (.)1 slo dependen de la temperatura reducida segn las expresiones siguientes:
10

N*
6,00 5,20 6,89 7,51 3,36 4,84 4,48 5,52 9,83 18,87 9,48 18,92 9,11 27,44 14,23 12,38 13,42 22,66 12,25 12,01 22,52 9,73 19,21 8,53 12,99 9,97 11,50 7,56 12,70 8,69 15,04 21,58 8,57 10,48 13,00 5,21 8,15 12,36 6,38 13,08 10,79 18,10

B
377,82 351,95 437,37 446,89 211,21 307,40 285,89 347,07 527,30 889,40 501,80 994,10 563,09 1008,70 1113,00 1141,35 1699,10 1399,71 665,40 652,02 1123,29 580,16 1133,17 514,53 645,92 575,96 656,83 545,01 790,47 605,44 904,08 1041,18 442,82 567,43 728,79 319,94 435,85 589,18 383,16 391,19 554,87 931,07

To
209,21 189,83 229,29 242,41 140,15 180,68 171,26 197,74 285,70 392,45 279,72 392,87 273,20 462,58 347,12 324,12 337,49 423,55 322,36 319,06 422,41 284,01 395,31 262,53 322,10 288,04 311,82 243,44 328,36 265,53 356,13 414,74 263,28 296,33 332,23 190,08 255,24 323,86 217,97 333,25 301,19 385,79

11

1-Octeno 1,3-butadieno 2-metil-2-buteno 2-metil-1,3-butadieno

donde: e: Parmetro de viscosidad w: Factor acntrico el parmetro de viscosidad (e) viene definido por la teora de los estados correspondientes propuesta por Yoon y Todos y se determina como:
12

2,3-dimetil-1-buteno n-butilciclo-pentano Tridecilciclo-pentano Etilciclo-hexano Dodecil-ciclohexano o-xileno m-terfenilo 1-pentanol Alcohol isoproplico 2-metil 2-butanol

donde: : Viscosidad (Pa.s) Tc: Temperatura crtica (K) Pc: Presin crtica (Pa) M: Masa molar (kg/kmol) Tr: Temperatura reducida esta ecuacin debe usarse con mucha precaucin, ya que puede llegar a dar errores de ms de un 30%. el error aumenta al acercarse a las regiones crticas.

Propilenglicol Acido n-butrico Acido isobutrico Acido fenil-actico Valerato de etilo Benzoato de benzilo Metil n-butil cetona Acetofenona Etil hexil ter Propil fenil ter Propilamina Bencilamina Etilpropilamina m- Toluidina Benzilfenilamina Nitrometano 2-nitro-2-2penteno Nitrobenceno Cloruro de etilo 1-bromo-2-metilpropano Iodobenceno Aldehdo proplico Benzaldehdo Anhdrido propinico Acetamida

4.

MezClas de lQuidos a bajas teMperaturas

las viscosidades de mezclas lquidas suele presentar un comportamiento fuertemente no lineal y no existe una teora general que aporte buenas predicciones.

4.1. Ecuaciones basadas en procedimientos estructurales

la ecuacin ms recomendada para la determinacin de la viscosidad de mezclas lquidas a bajas temperaturas es la propuesta por Grunberg y nissan [1]. ecuacin de Grunberg y nissan:
13

abril 2010

ingeniera QuMiCa 99

TRANSPORTE DE FLUIDOS

donde: xi: Fraccin molar. si todos los componentes de la mezcla son no polares, se recomienda utilizar fracciones msicas. Gij: Parmetro de interaccin. la ecuacin anterior aplicada a mezclas binarias queda como:
14

5.

MezClas de lQuidos a altas teMperaturas

Para la prediccin de la viscosidad de mezclas de hidrocarburos a altas temperaturas, Perry (1999) recomienda la correlacin de letsou y stiel (ecuacin 9), que se utiliza para la determinacin de la viscosidad de lquidos puros a elevadas temperaturas, sustituyendo en esta ecuacin el trmino de viscosidad del componente puro () por la viscosidad de la mezcla (m).

si no se dispone de datos experimentales para determinar el parmetro de interaccin (Gij), se puede utilizar un valor medio vlido para grandes familias segn se recomienda en la Tabla 3 ecuacin de Kendall Monroe. esta ecuacin es aplicable para la determinacin de la viscosidad de mezclas de lquidos miscibles no electrolticos y no asociados, cuya diferencia de viscosidades (1 - 2) < 15 mPas [5]. adems, Perry (1999) recomienda esta ecuacin para la determinacin de la viscosidad de mezclas de lquidos hidrocarburos:
15

6.

suspensiones lQuido - slido

Para la estimacin de la viscosidad de suspensiones lquido slido, se proponen las ecuaciones de Mori y ototake y de Hatschek [5]. ecuacin de Mori y ototake:
17

donde: l: Viscosidad de la fase lquida s: Fraccin volumtrica de slidos esta ecuacin es vlida para s < 0,1. ecuacin de Hatschek:
18

donde: xi: Fracciones molares de los lquidos que conforman la mezcla i: Viscosidades de los lquidos puros a la temperatura de la mezcla Para la determinacin de la viscosidad de mezclas de lquidos no hidrocarburados, [6] recomienda la siguiente ecuacin: ecuacin recomendada por Perry (1999):
15

esta ecuacin es vlida para 0,5 < s < 0,9.

7.

ConClusiones
se muestra un compendio de los mtodos de estimacin de la viscosidad de lquidos puros, mezclas de lquidos y suspensiones lquido-slido a baja y altas temperaturas comunicados por diferentes autores a partir de la teora cintica de los lquidos, del teorema de los estados correspondientes y de ecuaciones basadas en procedimientos estructurales.
bibliografa [1] aucejo, a. (1999). introducci a ienginyena qumica. biblioteca universitaria. barcelona. 1999. isbN: 84-7306-556-5. [2] betancourt, r. (2003). transferencia Molecular de calor, Masa y/o cantidad de Movimiento. cap 3 estimacin de propiedades de transporte. v.1. p.400. http://www. virtual.unal.edu.co/cursos/sedes/manizales/4090013/links/Pdf/cap_3.pdf [3] bird, r.b.; stewart, W.e.; ligthfood, e.N. (1992) fenmenos de transporte. barcelona: editorial revert, (espaa), isbN 84-291-7060-2. [4] brodkey, r.s.; Hershey, H.c. (1988) transport phenomena - a unified approach. New york (usa): Mcgraw-Hill book company. Mcgraw-Hill: chemical engineering series, isbN 0-07-007963-3. [5] garcell, l. (1992). transferencia de cantidad de movimiento, calor y masa. ciudad de la Habana: editorial Pueblo y educacin [6] Perry r.H. (1997). Perry's chemical engineers Handbook New york (usa): McgrawHill, inc. isbN 0-07-049841-5. (7th ed.). [7] Perry r.H. (1984). Perry's chemical engineers Handbook New york (usa): McgrawHill, inc. isbN 0-07-049841-5. (6th ed.).

donde: xi: Fracciones molares de los lquidos que conforman la mezcla i: Viscosidades de los lquidos puros a la temperatura de la mezcla

Tabla 3

valor Medio del ParMetro de iNteracciN (gij) Para graNdes faMilias

Mezcla
no polar + no polar no polar + polar polar + polar acuosa

Gij
-0,089 -0,218 0.251 1,66

100 ingeniera QuMiCa

N 481