Operaciones Unitarias PDF
Operaciones Unitarias PDF
Operaciones Unitarias PDF
Operaciones
Unitarias I
Grupo BQC-9325 Curso Enero-Junio 2007
OPERACIONES UNITARIAS I
O
OBBJJE
ETTIIV
VOO
C
Coom mpprreennddeerr llooss ffuunnddaam meennttooss ddee llaa ttrraannssffeerreenncciiaa ddee ccaannttiiddaadd ddee
mmoovviim miieennttoo yy ddee llaa m meeccáánniiccaa ddee fflluuiiddooss,, rreellaacciioonnáánnddoollooss yy
aapplliiccáánnddoollooss aa llaass ooppeerraacciioonneess uunniittaarriiaass eenn llaass qquuee ééssttooss ssee vveenn
iinnvvoolluuccrraaddooss
ESTAS NOTAS NO PRETENDEN SER UN REEMPLAZO A LOS LIBROS DE TEXTO RECOMENDADOS,
EN SU LUGAR REPRESENTAN UNA GUIA PARA EL ESTUDIANTE EN SU INTRODUCCIÓN A LA
MATERIA DE OPERACIONES UNITARIAS I, DE LA CARRERA DE INGENIERIA BIOQUÍMICA
IMPARTIDA EN EL INSTITUTO TECNOLÓGICO DE ZACATEPEC.
TEMARIO
CONSIDERACIONES PRELIMINARES...............................................1
UNIDAD I (VISCOSIDAD).......................................................................3
1.1 Procesos de Transferencia. Definición, Tipos, Gradientes y Parámetros................................................................3
EJERCICIOS ...................................................................................................................................................................29
TAREAS ...........................................................................................................................................................................30
DR. J. PORCAYO-CALDERON i
OPERACIONES UNITARIAS I
EJERCICIOS ...................................................................................................................................................................53
TAREAS ...........................................................................................................................................................................55
ANEXOS...........................................................................................................................................................................56
LA ECUACIÓN DE CONTINUIDAD Y LA ECUACIÓN DE MOVIMIENTO EN COORDENADAS CARTESIANAS,
CILÍNDRICAS Y ESFERICAS. ......................................................................................................................................56
COMPONENTES DEL TENSOR ESFUERZO DE CORTE EN COORDENADAS CARTESIANAS, CILÍNDRICAS Y
ESFERICAS. ..................................................................................................................................................................59
EJERCICIOS ...................................................................................................................................................................78
TAREAS ...........................................................................................................................................................................79
EJERCICIOS ...................................................................................................................................................................97
TAREAS ...........................................................................................................................................................................99
ANEXOS.........................................................................................................................................................................100
What is Pumping Head? ..............................................................................................................................................100
DR. J. PORCAYO-CALDERON ii
OPERACIONES UNITARIAS I
BIBLIOGRAFIA
Welty, Wicks y Wilson McNaughtun Kenneth
Fundamentos de Transferencia de Momentum, Bombas
Calor y Masa Ed. McGraw-Hill
Ed. Limusa Greene Richard W.
Foust A. S. and Wenzel L. A. Vávulas
Principios de Operaciones Unitarias Ed. McGraw-Hill
Ed. CECSA Brown
Perry and Chilton Operaciones Unitarias
Manual del Ingeniero Químico Ed. McGraw-Hill
Ed. McGraw-Hill Leva Max
Benett C. D. and Meyrs J. E. Fluidización
Momentum, Heat and Mass Transfer Kunni
Ed. McGraw-Hill Fluidización
Geankoplis Christie J.
Procesos de Transporte y Operaciones Unitarias
Ed. CECSA
CONSIDERACIONES PRELIMINARES
La Ingeniería Química tratan con el cambio de sustancias (gases, líquidos, sólidos) por medio de procesos físicos,
químicos y biológicos para propósitos industriales. 1
Esta definición muestra que la ingeniería química es un campo muy amplio y diverso. La pregunta surge al cuestionarse:
Puede una sola disciplina científica tratar un campo tan amplio?.
La primer respuesta se dio en 1915 por Arthur D. Little cuando acuño el termino "Operación Unitaria":
La introducción y definición del término "Operación Unitaria" por Arthur D. Little, fue el primer paso en el intento de
reducir la variedad de procesos técnicos a las leyes fundamentales de matemáticas, física, química, y mecánica. Al
mismo tiempo, impulsó el desarrollo para reemplazar el conocimiento empírico por modelos matemáticos que den una
descripción aproximada de los procesos reales. La tabla I muestra un resumen de las operaciones unitarias más
importantes. Note que las operaciones unitarias que involucran cambios químicos son referidas como "Procesos
Unitarios", término acuñado por R.N. Shreve.
El posterior desarrollo del concepto de "Operaciones Unitarias" mostró que los procesos de abstracción se pueden
realizar analizando las operaciones unitarias en términos de principios fundamentales como balances de masa y
energía, equilibrio de fases, transporte de momentum, energía y masa. El inicio de esta era se atribuye a la publicación
del libro "Transport Phenomena" de Bird, Stewart y Lightfoot en 1960. La tabla II muestra un resumen de estos
principios fundamentales y su relación con las diferentes clases de operaciones unitarias.
1
Esta sección fue tomada de http://chemengineer.mininco.com/library/weekly/aa033098.htm
DR. J. PORCAYO-CALDERON 1
OPERACIONES UNITARIAS I
(UNIDAD I)
Como podrá observarse de acuerdo al temario de este curso, el contenido de la materia OU-I nos envía al estudio de
los fenómenos de transporte, en particular a la transferencia de momentum.
DR. J. PORCAYO-CALDERON 2
OPERACIONES UNITARIAS I
(UNIDAD I)
UNIDAD I (Viscosidad)
1.1 Procesos de Transferencia. Definición, Tipos, Gradientes y Parámetros 2
Los fenómenos de transporte se refieren a la transferencia de momentum, calor y masa de un punto del universo a otro.
Las ecuaciones de transporte son ecuaciones dimensionales e involucran varios términos. Para que se sostenga la
igualdad en cada término de la ecuación, cada término debe ser dimensionalmente homogéneo.
El desarrollo de las ecuaciones de transporte implican varios tipos de notaciones, las más importantes de ellas son:
¾ Cantidades Escalares. Representadas por letras ordinarias o letras griegas, p.e., s, F, φ, μ, etc.
¾ Vectores. Son representados por letras ordinarias en negrillas, p.e., v, w, A, etc.
¾ Tensores. Se representaran por letras griegas en negrillas, p.e., τ, Δ.
¾ Valores Absolutos. Los valores absolutos de cantidades escalares y las magnitudes de los vectores se
designarán por la notación ⏐⏐, p.e., ⏐-5⏐o⏐v⏐.
¾ Cantidades Especiales. Se designaran por el superíndice * arriba del símbolo, p.e., p*.
¾ Valores Promedio. Se indicará con la notación < >, p.e., <T>, <v>, <θ>, etc.
Desarrollo de Ecuaciones
Los fenómenos de transporte tratan con la transferencia de ciertas cantidades (momentum, masa y energía) de un punto
del sistema a otro. Los tres mecanismos básicos de transporte involucrados en un proceso, son:
Estos tipos de difusión molecular se pueden describir por leyes fenomenológicas (modelos matemáticos). Cada una de
estas leyes se reducen al producto del coeficiente de transporte apropiado y su gradiente.
El momentum, calor y masa, son cantidades que se conservan. Es decir obedecen la ley de la conservación en un
sistema:
La ecuación anterior se puede representar en función del tiempo, con lo que se obtendrían velocidades:
2
Transport Phenomena for Engineers, Author: Louis Theodore, Ed. International Textbook Company (1971), pp 4-9.
DR. J. PORCAYO-CALDERON 3
OPERACIONES UNITARIAS I
(UNIDAD I)
⎧Vel. de entrada ⎫ ⎧ Vel. de salida ⎫ ⎧Vel. de generación⎫ ⎧Vel. de acumulación ⎫
⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪
⎨de la Cantidad ⎬ − ⎨de la Cantidad ⎬ + ⎨ de la Cantidad ⎬ = ⎨ de la Cantidad ⎬
⎪ en el Sistema ⎪ ⎪ en el Sistema ⎪ ⎪ en el Sistema ⎪ ⎪ en el Sistema ⎪
⎩ ⎭ ⎩ ⎭ ⎩ ⎭ ⎩ ⎭
La ley de la conservación a un nivel macroscópico, microscópico o molecular. Se ilustrarán las diferencias de estos
niveles con el siguiente ejemplo. Considere un sistema en el cual fluye un fluido a través de un tubo cilíndrico.
Definamos el sistema como el fluido contenido en el tubo entre los puntos 1 y 2 en cualquier instante (tiempo).
¾ La aproximación macroscópica, se aplica cuando se esta interesado en determinar los cambios que
ocurren a la entrada y salida del sistema. La ecuación resultante describiría los cambios globales que
ocurren en el sistema sin tomar en cuenta las variaciones internas en el sistema. Este tipo de
aproximaciones se aplica generalmente en Operaciones Unitarias.
De lo explicado anteriormente se puede concluir que los análisis que se tratarán en este curso corresponderán a
aproximaciones microscópica.
Condiciones iniciales y frontera. Para poder resolver las ecuaciones diferenciales (derivadas de la aplicación de ley de
la conservación a aproximaciones microscópicas), es necesario especificar las condiciones iniciales y/o frontera para el
sistema. Esta información surge de la descripción del problema o la situación física.
¾ El número de condiciones frontera que se deben especificar es la suma de la derivada de orden más alto por
cada variable independiente que aparezca en la ecuación diferencial. Un valor de la solución en la frontera
del sistema es una condición frontera.
¾ El número de condiciones iniciales que se deben especificar es el orden de la derivada de tiempo más alta
que aparece en la ecuación diferencial. Esta condición es usada solamente si el tiempo es una variable. El
valor de la solución al tiempo igual a cero es una condición inicial.
DR. J. PORCAYO-CALDERON 4
OPERACIONES UNITARIAS I
(UNIDAD I)
Solución de ecuaciones. El procedimiento general que se debe seguir para resolver problemas en fenómenos de
transporte, requiere primeramente la transformación de la palabra ecuación en expresiones matemáticas, y su solución.
El procedimiento sugerido es como sigue:
3
Transport Phenomena for Engineers, Author: Louis Theodore, Ed. International Textbook Company (1971), pp 67-72.
DR. J. PORCAYO-CALDERON 5
OPERACIONES UNITARIAS I
(UNIDAD I)
El signo menos, es debido a que el momentum es transferido en la dirección negativa de z como consecuencia de un
gradiente de velocidad positivo.
Entonces se puede establecer que todo aquel fluido cuyo esfuerzo de corte se describe por la ecuación anterior se define
como un fluido newtoniano.
Unidades 4. Las unidades de los términos de la ley de Newton son las siguientes:
z cm m ft
1 dina = gm-cm/s 2
4
Fenómenos de Transporte, Bird, Steward, Lightfoot.
DR. J. PORCAYO-CALDERON 6
OPERACIONES UNITARIAS I
(UNIDAD I)
μ kgm/m-s lbm/ft-s
νy m/s ft/s
z m ft
2
gc 9.8067 kgm-m/kgf-s 32.2 lbm-ft/lbf-s2
Donde kgm-m/kgf-s2 = newton/kgf
A menudo es conveniente dividir la viscosidad del fluido por su densidad (v = μ/ρ) obteniendo un coeficiente cinemático
de la viscosidad. Esta viscosidad cinemática es una medida de la resistencia del fluido de la transición de flujo laminar a
turbulento. Es decir, mide la resistencia del fluido a efectos de variables de flujo como la velocidad. La viscosidad
cinemática la indica la velocidad de difusión de momentum y tiene unidades de cm2/s (1 stoke = 1 cm2/s).
¾ Los fluidos con baja viscosidad cinemática soportan turbulencia rápidamente que aquellos fluidos con alta
viscosidad cinemática.
¾ Para un fluido newtoniano la viscosidad es una propiedad de estado ya que es una función de la presión,
temperatura, y composición.
EJEMPLO 1. 1
Calcular la densidad de flujo de C.M. en estado estacionario, τyx, en kgf/m2, cuando la velocidad de la lamina
inferior, en la dirección positiva del eje x, es 0.3 m/s, la distancia entre las láminas es 0.0003 m, y la viscosidad del
fluido es de 0.7 cp.
¾ Flujo laminar es un flujo ordenado que ocurre cuando las capas adyacentes de fluido resbalan unas sobre
otras, la mezcla entre capas se presenta solo a nivel molecular. La fuerza dominante para la transferencia
de momentun son las fuerzas viscosas.
5
Fundamentos de Transferencia de Momento, Calor y Masa. Welty, Wicks y Wilson, Ed. Limusa. Pag 203.
DR. J. PORCAYO-CALDERON 7
OPERACIONES UNITARIAS I
(UNIDAD I)
La existencia de los dos tipos de flujo fue descrito
cualitativamente por Reynolds en 1883. En su experimento
clásico hizo que fluyera agua, a una velocidad controlada, a
través de un tubo transparente. En un punto del tubo
introdujo un colorante cuya gravedad específica era la
misma del agua y observó el patrón de flujo a flujos cada
vez mayores.
Para flujo en tuberías se ha encontrado que el régimen de flujo puede definirse por el valor adimensional Número de
Reynolds con base a las siguientes reglas generales:
DR. J. PORCAYO-CALDERON 8
OPERACIONES UNITARIAS I
(UNIDAD I)
1.4 Fluidos Newtonianos, No-Newtonianos y Modelos Reológicos 67
Para propósitos de este curso, los fluidos se pueden dividir en tres grupos de acuerdo a su comportamiento de flujo:
⎛ dv y ⎞
τ zy = μ ⎜⎜ − ⎟ = μγ
⎟
⎝ dz ⎠
Donde:
τzy = esfuerzo de corte
μ = viscosidad.
⎛ dv y ⎞
⎜− ⎟
⎜ dz ⎟ = γ = gradiente de velocidad o velocidad de corte
⎝ ⎠
2. Fluidos No-Newtonianos Independientes del Tiempo. La viscosidad de este tipo de fluidos es dependiente de
la temperatura y de la velocidad de corte, y como su nombre lo indica son ésta no depende del tiempo.
Dependiendo de cómo cambia la viscosidad con la velocidad de corte, el comportamiento de flujo estos se sub-
clasifican en:
i) Peudoplásticos (Esfuerzo-Adelgazantes). Su
viscosidad disminuye al incrementar la velocidad de corte.
(Pinturas, shampoos, lechadas, jugos concentrados de
fruta).
n
⎛ dv y ⎞
τ zy = m⎜⎜ − ⎟ = mγ n
⎟
⎝ dz ⎠
Donde:
n = índice de comportamiento de flujo o grado de desviación del comportamiento newtoniano.
m = índice de consistencia, es equivalente a μ.
Para fluidos pseudoplásticos n < 1.
Para fluidos dilatantes n > 1.
Cuando n = 1, la ecuación anterior se reduce a la ley de Newton de la viscosidad
6
Transfer Operations. R.A. Greenkorn and D.P. Kessler, Ed. McGraw-Hill, pp176.
7
Fenómenos de Transporte. Bird, Steward, Lightfoot, pp 1-10.
DR. J. PORCAYO-CALDERON 9
OPERACIONES UNITARIAS I
(UNIDAD I)
⎛ dv y ⎞
log(τ zy ) = n log⎜⎜ ⎟⎟ + log m
⎝ dz ⎠
⎛ dv y ⎞
τ zy = τ o + m⎜⎜ − ⎟ = τ o + mγ
⎟
⎝ dz ⎠
Donde:
τo = esfuerzo de cedencia.
iv) Fluidos Hershey-Bulkley. El comportamiento de estos
fluidos es del tipo llamado cuasiplástico o mixto. El
modelo reológico que describe este comportamiento esta
descrito por la siguiente ecuación, llamada también
Modelo Generalizado el cual describe todos los
comportamientos de flujo antes citados:
n
⎛ dv y ⎞
τ zy = τ o + m⎜⎜ − ⎟ = τ o + mγ n
⎟
⎝ dz ⎠
Otro modelo propuesto es el Modelo de Casson el cual se adapta a la mayoría de las suspensiones:
1/ 2
⎛ dv y ⎞ 1
τ zy = m⎜⎜ − ⎟
⎟ + τ 10 / 2 = mγ n + τ o 2
⎝ dz ⎠
En general se puede decir que el comportamiento reológico de cualquier fluido independiente del tiempo se
puede expresar por medio de los modelos H-B o Casson.
DR. J. PORCAYO-CALDERON 10
OPERACIONES UNITARIAS I
(UNIDAD I)
En este caso para poder determinar los valores
de m, n y τo, primero debemos conocer el
valor de τo. Para ello tenemos tres casos:
⎛ dv y ⎞
log(τ zy − τ o ) = n log⎜⎜ ⎟⎟ + log m
⎝ dz ⎠
DR. J. PORCAYO-CALDERON 11
OPERACIONES UNITARIAS I
(UNIDAD I)
3. Fluidos No-Newtonianos Dependientes del Tiempo. La viscosidad de este tipo de fluidos depende de la
temperatura, velocidad de corte y tiempo. Dependiendo de cómo cambia la viscosidad éstos se caracterizan
como:
¾ Tixotrópicos (Tiempo-Adelgazantes).
Muestran una disminución gradual en su
viscosidad aparente en función del tiempo al
aplicarles repentinamente un esfuerzo
cortante (Yogurt, pinturas).
¾ Reopécticos (Tiempo-Espesantes).
Muestran un incremento en su viscosidad
aparente respecto al tiempo al aplicarles
repentinamente un esfuerzo cortante.
Una clasificación adicional que podemos tomar en consideración es la de los fluidos viscoelásticos.
Donde λ es el término elástico, y un término de comportamiento dependiente del tiempo, algunos señalan que
este tipo de fluidos tienen memoria ya que "recuerdan" que les sucedió (es decir qué esfuerzo se le aplicó y
cuanto se elongó para volver a su estado inicial. Un ejemplo de este tipo de fluidos es el cemento pegamento).
Algunas características particulares de este tipo de fluidos son 8:
¾ Rebobinado elástico
8
http://wwwce1.kat.lth.se/staff/ulf_b/b5_rheo.htm
DR. J. PORCAYO-CALDERON 12
OPERACIONES UNITARIAS I
(UNIDAD I)
¾ Aumento del jet
¾ Auto vaciado
¾ Efecto Weissenberg
Determinación de los Parámetros Reológicos. Los ejemplos siguientes ilustrarán los procedimientos de obtención de
los parámetros reológicos a partir de datos de esfuerzo de corte y velocidad de corte o gradiente de velocidad.
EJEMPLO 1. 2
DR. J. PORCAYO-CALDERON 13
OPERACIONES UNITARIAS I
(UNIDAD I)
SOLUCION. Al graficar los datos se observa que 0.0025
el esfuerzo de corte es prácticamente cero a
gradiente de velocidad (velocidad de corte) cero.
0.002
0.001
0.0005
0
0 0.001 0.002 0.003 0.004 0.005 0.006 0.007 0.008
Velocidad de corte o gradiente de velocidad [1/sec]
donde:
y = log τ yx a=n
⎛ dv ⎞ b = log m
x = log⎜⎜ − x ⎟⎟
⎝ dy ⎠
De lo anterior se puede observar que el tratamiento del modelo reológico de esta manera, nos puede proporcionar los
valores del coeficiente de consistencia (m) y del índice de comportamiento de flujo (n). Por tanto tomando logaritmos a
los valores experimentales y graficando éstos de acuerdo a la ecuación de la recta, obtenemos el siguiente gráfico:
DR. J. PORCAYO-CALDERON 14
OPERACIONES UNITARIAS I
(UNIDAD I)
-2.6
-3.5 -3 -2.5 -2 -1.5 -1 -0.5 0 • De la pendiente de la recta se obtiene
-2.65 que el valor del n es 0.387.
• El valor de la ordenada al origen es –
LOG (Esfuerzo de corte [dinas/cm ])
2
-3
LOG (Velocidad de corte o gradiente de velocidad [1/sec])
EJEMPLO 1. 3
(Reoeje1). Considere el conjunto de datos de velocidad de corte y esfuerzo cortante dados en la tabla siguiente:
0.00016
0.00014
De acuerdo al gráfico se puede asumir que el
Esfuerzo de Corte (Pa)
0.00012
comportamiento del fluido corresponde a un fluido
dilatante de acuerdo a la siguiente expresión: 0.0001
0.00008
n
⎛ dx ⎞
0.00006
τ yx = m⎜⎜ − ⎟⎟ 0.00004
⎝ dy ⎠ 0.00002
0
0 0.002 0.004 0.006 0.008 0.01 0.012 0.014
Velocidad de Corte (s-1)
DR. J. PORCAYO-CALDERON 15
OPERACIONES UNITARIAS I
(UNIDAD I)
donde:
⎛ dv ⎞
y = logτ yx x = log⎜⎜ − x ⎟⎟ a=n b = log(m)
⎝ dy ⎠
De lo anterior se puede observar que el tratamiento de esta manera del modelo reológico, nos puede proporcionar los
valores del coeficiente de consistencia (m) y del índice de comportamiento de flujo (n). Por tanto tomando logaritmos a
los valores experimentales y graficando éstos de acuerdo a la ecuación de la recta, obtenemos el siguiente gráfico:
-3
-3.1 -2.9 -2.7 -2.5 -2.3 -2.1 -1.9 -1.7 -1.5
y = 1.7002x - 0.4553 -4
R2 = 1 tanto el valor de m es 0.3505 Pa-sn
-4.5
Por tanto el modelo reológico es:
-5
1.7002
⎛ dv x ⎞
τ yx = 0.3505⎜⎜ − ⎟⎟
⎝ dy ⎠
-5.5
-6
Log Velocidad de Corte
Viscosidad Aparente.
La relación τ/γ se considera la viscosidad aparente del fluido. De acuerdo a un análisis de la tendencia de los modelos
reológicos se observa que ésta varía con la velocidad de corte.
DR. J. PORCAYO-CALDERON 16
OPERACIONES UNITARIAS I
(UNIDAD I)
Efecto de la Temperatura.
Un parámetro ambiental que influye considerablemente en la viscosidad o cualquier parámetro reológico, es la
temperatura. El modelo más utilizado para considerar este efecto es el uso de la ecuación de Arrhenius:
Donde:
E k = constante de velocidad de reacción.
− E = energía de activación.
k= Ae RT A = constante de integración referida como factor de
frecuencia.
T = temperatura absoluta.
EJEMPLO 1. 4
(Heldman2.2)El coeficiente de consistencia para puré de durazno varía con la temperatura de acuerdo a los
siguientes datos.
Temperatura m
[°F] [dinas-secn/cm2]
40 130
77 90
130 45
140 38
Desarrolle una relación entre la temperatura y el parámetro reológico usando la ecuación tipo Arrhenius.
• De donde: A = 0.08386
3.5
3669.7
• Por tanto la relación es: m = 0.08386e T 3
0.00165 0.0017 0.00175 0.0018 0.00185 0.0019 0.00195 0.002 0.00205
1/T
Donde T está definida en °R.
DR. J. PORCAYO-CALDERON 17
OPERACIONES UNITARIAS I
(UNIDAD I)
En el gráfico siguiente se muestran los valores 140
100 Datos
Experimentales
Indice de Consistencia
80
Datos Modelo
60
40
20
0
0 20 40 60 80 100 120 140 160
Temperatura [°F]
¾ Modelo de Eyring. Modelo de dos parámetros derivado de la teoría cinética de los líquidos, y predice el
comportamiento pseudoplástico para valores finitos de τ, y tiende asintóticamente a la ley de la viscosidad
de Newton cuando τ tiende a cero ( en este caso μ=A/B):
⎛ 1 dv y ⎞
τ zy = A arcsen h⎜⎜ − ⎟
⎝ B dz ⎟⎠
¾ Modelo de Ellis. Modelo de tres parámetros positivos ajustables. Si alfa mayor que 1, el modelo tiende a la
ley de Newton para valores bajos de τzy. Si alfa menor que 1, la ley de Newton se establece para valores
elevados de τzy. El modelo presenta gran flexibilidad y en el están comprendidas, casos particulares, la ley
de Newton (para ϕ1 = 0), y la ley de la potencia (para ϕ0 = 0).
dv y α −1 ⎞
− = ⎛⎜ ϕ 0 + ϕ 1 τ zy α
⎟ τ zy = ϕ oτ zy + ϕ 1τ zy
dz ⎝ ⎠
¾ Modelo de Reiner-Philippoff. Modelo de tres parámetros positivos ajustables. Este modelo incluye la
observación de que el comportamiento Newtoniano se presenta tanto a valores muy bajos como muy
elevados del gradiente de velocidad. La curva resultante muestra dos puntos de inflexión.
⎛ ⎞
⎜ ⎟
dv y ⎜ 1 ⎟
− =⎜ ⎟ τ zy
dz ⎜μ + μ 0 − μ∞ ⎟
⎜ ∞ 1+ τ τ
⎝ zy s ( )2 ⎟
⎠
DR. J. PORCAYO-CALDERON 18
OPERACIONES UNITARIAS I
(UNIDAD I)
1.5 Estimación de Viscosidad en Gases 910
A continuación se describen brevemente algunos métodos empleados para estimar la viscosidad en gases:
μ c = 61.6(MTC )1 / 2 (V ´c )−2 / 3 ó
μ c = 7.7 M 1 / 2 Pc2 / 3Tc−1 / 6
EJEMPLO 1. 5
Calcular la viscosidad del N2 a 50°C y 854 atm, siendo M = 28 g/g-mol, Pc = 33.5 atm y Tc = 126.2K.
9
Fenómenos de Transporte. Bird, Steward, Lightfoot, Capítulo I.
10
Perry´s Chemical Engineers´ Handbook, Seventh Edition.
DR. J. PORCAYO-CALDERON 19
OPERACIONES UNITARIAS I
(UNIDAD I)
El valor experimental reportado es 0.000455 g/cm-s.
Para calcular la viscosidad de una mezcla de n componentes, mediante la figura 1.3-1 del BSL, se utilizan las
propiedades pseudocríticas definidas empíricamente por las ecuaciones:
n n n
Pc´ = ∑ xi Pci Tc´ = ∑ xi Tci μ c´ = ∑ xi μ ci
i =1 i =1 i =1
Este método no es muy exacto si las sustancias que forman la mezcla tienen distinta constitución química o sus
propiedades críticas difieren notablemente.
EJEMPLO 1. 6
La viscosidad del CO2 a 45.3 atm y 40.3°C es
0.00018 poise. Estimar el valor de la viscosidad a
114.6 atm y 40.3°C usando la segunda gráfica de
esta sección. Considere Pc = 72.9 atm y Tc = 304.2K
DR. J. PORCAYO-CALDERON 20
OPERACIONES UNITARIAS I
(UNIDAD I)
• Gases Puros a baja presión. Método de Stiel y NM 0.5 Pc2 / 3
Thodos μ v = 0. 00046
Para predecir la viscosidad de vapores de Tc1 / 6
hidrocarburos puros a baja presión (debajo de Tr =
0.6), el método de Stiel y Thodos es el más exacto. N = 0.00034Tr0.94 para Tr ≤ 1.5
Sólo se requiere el peso molecular, la temperatura y
presión crítica. Las unidades de la temperatura están
en K, la presión en Pa, y la viscosidad en cp. N = 0.0001778(4.58Tr − 1.67 )0.65 para Tr > 1.5
EJEMPLO 1. 7
Estime la viscosidad del vapor de propano a 101.3 kPa y 80°C. Si la temperatura y presión crítica es de 96.7°C y
4.246 Mpa, respectivamente y el peso molecular del propano es de 44.1.
y:
NM 0.5 Pc2 / 3 (0.0003255)(44.1) 0.5 (4,246,000)2 / 3
μ v = 0.00046 = 0.00046 = 0.0097cp
Tc1 / 6 369.81 / 6
⎛
μ = ⎜ 5 * 10 −8
⎜
M 0.5 Pc2 / 3
Tc1 / 6
[exp(1.439 ρ r ) − exp(− 1.11ρ 1r .858 )]⎟⎟ + μ o
⎞
⎝ ⎠
La presión y temperatura crítica están en Pa y K, la viscosidad en cp. La viscosidad con subíndice cero es la
viscosidad del gas o mezcla de gases determinada de acuerdo a los métodos anteriores considerando baja presión.
ATr
μ=
[1 + 0.36Tr (Tr − 1)]1 / 6
DR. J. PORCAYO-CALDERON 21
OPERACIONES UNITARIAS I
(UNIDAD I)
Para compuestos orgánicos:
M 1 / 2Tc
A=
∑ ni C i
Para gases inorgánicos:
⎡ M 1 / 2 Pc2 / 3 ⎤
A = 1.6104 * 10 −10 ⎢ 1/ 6 ⎥
⎣⎢ Tc ⎦⎥
Las viscosidades se dan en Pa-s (1000 cp = 1 Pa-s) con la Tc = K y Pc = Pa. Los valores de C (contribución de
grupos) para compuestos orgánicos se dan en la tabla 1.III.
ρ r ≤ 0.1 ⇒ (μ − μ )B = 1.656 * 10 ρ
o −7 1.111
r
Una extensión a este método para gases no-polares es el propuesto por Jossi:
DR. J. PORCAYO-CALDERON 22
OPERACIONES UNITARIAS I
(UNIDAD I)
0.267
-0.153
0.768
0.553
0.178
0.741
0.524
0.691
0.116
0.023
0.590
0.359
0.446
1.006
1.283
0.796
0.359
1.202
1.402
1.865
1.341
0.971
0.368
0.497
1.813
0.886
11
Perry´s Chemical Engineers´ Handbook, Seventh Edition
DR. J. PORCAYO-CALDERON 23
OPERACIONES UNITARIAS I
(UNIDAD I)
1.6 Estimación de Viscosidad en Líquidos
• Teoría de Eyring 12
En un líquido puro en reposo sus moléculas están en
movimiento constante, pero debido a su compacto
“empacamiento” este movimiento se debe prácticamente a la
vibración de cada molécula dentro de una “jaula” o retícula
formada por las moléculas más próximas.
KT − ΔGo+ + / RT
k= e
h
Donde K y h son las constantes de Boltzmann y Plank, R es la constante molar de los gases, y ΔGo++ es la energía
libre de activación molar del fluido en reposo.
Donde V´ es el volumen molar del líquido, ±(a/δ)τxyV´/2 es el trabajo transmitido a las moléculas que se mueven
hacia la cima de las barreras de energía (“+” si se mueven a favor del esfuerzo cortante aplicado o “-“ si es contra).
12
Fenómenos de Transporte. Bird, Steward, Lightfoot, pp 1-10.
DR. J. PORCAYO-CALDERON 24
OPERACIONES UNITARIAS I
(UNIDAD I)
Donde Ñ representa el número de Avogadro. Generalmente el
término (δ/a) es la unidad. La expresión anterior también se ha Ñh 3.8Tb / T
correlacionado con la temperatura normal de ebullición, μ= e
V´
obteniéndose la siguiente expresión:
Ambas expresiones indican que la viscosidad disminuye exponencialmente con la temperatura, lo cual está de
acuerdo con el comportamiento observado con la mayoría de los líquidos. Estas ecuaciones no son muy exactas
(hasta 30% de error) pero son útiles para una estimación aproximada. No es recomendable su aplicación para
moléculas lineales muy largas, p.e. n-C20H42.
No se sugiere el uso de este método para el primer miembro de una serie homóloga o para temperaturas arriba del
punto normal de ebullición (Tr ≈ 0.75)
EJEMPLO 1.1
13
Perry´s Chemical Engineers´Handbook, Seventh Edition.
DR. J. PORCAYO-CALDERON 25
OPERACIONES UNITARIAS I
(UNIDAD I)
Por tanto en valor de B, es: B = Ba + ∑ ΔBi = 537.26 + (1)[− 272.85 + 25.041(9.48)] = 501.8
i
El valor estimado de la viscosidad es: ⎛1 1 ⎞ ⎛ 1 1 ⎞
log μ = B⎜⎜ − ⎟⎟ − 3.0 = 501.8⎜ − ⎟ − 3.0 = −2.9569
⎝ T To ⎠ ⎝ 273.15 279.72 ⎠
μ = 0.001104 Pa − s = 1.104cp
El valor experimental a 0°C es 1.224 cp.
14
Perry´s Chemical Engineers´Handbook, Seventh Edition.
15
Perry´s Chemical Engineers´Handbook, Seventh Edition.
DR. J. PORCAYO-CALDERON 26
OPERACIONES UNITARIAS I
(UNIDAD I)
16
Perry´s Chemical Engineers´Handbook, Seventh Edition.
DR. J. PORCAYO-CALDERON 27
OPERACIONES UNITARIAS I
(UNIDAD I)
DR. J. PORCAYO-CALDERON 28
OPERACIONES UNITARIAS I
(UNIDAD I)
EJERCICIOS
1. En la tabla siguiente se muestran los valores de velocidad de corte (dvy/dx) y esfuerzo cortante (τxy) para un
producto lácteo fermentado con un contenido de grasa de 0.5% determinados a 15°C (Datos extraídos de:
http://wwwce1.kat.lth.se/staff/ulf_b/b5_rheo.htm).
a) Grafique estos datos e indique a que modelo reológico se ajustan los datos experimentales.
b) Asumiendo que los datos experimentales de ajustan a un modelo Hershey-Bulkley determine los
parámetros y modelo reológicos considerando un esfuerzo de cedencia de 2 Pa.
c) Asumiendo que los datos experimentales de ajustan a un plástico de Bingham determine los parámetros y
modelo reológicos considerando un esfuerzo de cedencia de 2 Pa.
d) Asumiendo que los datos experimentales se ajustan al modelo de Ostwald-de Wale, determine los
parámetros y modelo reológicos necesarios para predecir el comportamiento del fluido.
e) De acuerdo a la correlación de los datos experimentales, indique que parámetros o parámetro reológico es
más sensible de error.
3. Predecir la viscosidad del oxígeno, nitrógeno y metano moleculares, a la presión atmosférica y 20°C. Expresar todos
los resultados en cp. Compare sus resultados con los valores experimentales que se indican en el capítulo 1 del texto
“Fenómenos de Transporte” de Bird, Steward and Lighfoot.
4. Se conocen los siguientes datos de la viscosidad de las mezclas de hidrógeno y freón-12 (diclorofluorometano) a
25°C y 1 atm. Calcular y comparar los resultados obtenidos mediante las ecuaciones siguientes, para las tres
presiones intermedias, utilizando los datos de la viscosidad de los componentes puros.
X1 μ [g/cm-s]
Fracción molar de H2
0.00 0.000124
0.25 0.000128
0.50 0.0001319
0.75 0.0001351
1.00 0.0000884
2
−1 / 2 ⎡ 1/ 2 1/ 4 ⎤
n xi μi 1 ⎛⎜ M ⎞ ⎛
⎢1 + ⎜ μ i
⎞ ⎛M j ⎞ ⎥
μ mezcla = ∑ φ ij = 1+ i ⎟ ⎟ ⎜
⎜M
⎟
⎟
n ⎜ M j ⎟⎠ ⎢ ⎜μj ⎟ ⎥
i =1
∑ x jφij
8⎝
⎢⎣ ⎝ ⎠ ⎝ i ⎠ ⎥⎦
j =1
5. Estimar la viscosidad del N2 a 20°C y 67 atm, utilizando (a) el gráfico de viscosidades reducidas contra temperatura
reducida y μc = 0.00018 poise, (b) el gráfico de viscosidades reducidas contra presión reducida y μo = 0.000175
poise. Expresar los resultados en kgm/m-s. (Faltan datos de Tc y Pc )
6. Estime la viscosidad de los vapores del alcohol isopropílico a 251°C y presión atmosférica.
DR. J. PORCAYO-CALDERON 29
OPERACIONES UNITARIAS I
(UNIDAD I)
TAREAS
Se establecerán en clases.
DR. J. PORCAYO-CALDERON 30