Informe Laboratorio Viscosidad de Fluidos Newtonianos-MF2 (Grupo F)
Informe Laboratorio Viscosidad de Fluidos Newtonianos-MF2 (Grupo F)
Informe Laboratorio Viscosidad de Fluidos Newtonianos-MF2 (Grupo F)
FACULTAD DE INGENIERÍA
DEPARTAMENTO DE INGENIERÍA QUÍMICA
CONCEPCIÓN
Informe
LABORATORIO DE PROCESOS QUÍMICOS
1
Sumario
Los fluidos en la actualidad tienen gran importancia dentro de la industria, ya que sus propiedades
son aprovechadas para el funcionamiento de maquinarias como combustibles o lubricantes, o bien
en la transferencia de calor dentro de un intercambiador, sin embargo, para su uso y
almacenamiento es necesario conocer sus características, en particular su viscosidad, es decir su
resistencia a fluir. El objetivo general de este trabajo es determinar la viscosidad de fluidos
newtonianos y observar cómo afectan los electrolitos a la viscosidad del solvente puro.
En este informe se presenta la viscosidad cinemática y dinámica de los fluidos con un análisis
cualitativo de los datos, además se comprobó que estos fluidos presentan un comportamiento
newtoniano.
2
Nomenclatura
3
1. Introducción
La viscosidad es una medida de la resistencia de un fluido, originada por la fricción interna entre
capa y capa, cuando una capa de fluido se mueve con respecto a otra. Esta propiedad debe ser
considerada en la operación de cualquier proceso, ya sea en el ámbito de la ingeniería química o la
industria en sí, ya que todos los fluidos que participan en estos, presentan comportamientos
característicos debido a su viscosidad. Por lo tanto para optimizar cualquier proceso es necesario
conocer este factor, entender de qué manera afecta y tomar las medidas necesarias para
optimizar el proceso y conseguir la mejor calidad del producto deseado. Algunos casos de
aplicación son por ejemplo en los equipos industriales, estos poseen ciertas piezas, tales como
engranajes, cilindros y rodamientos, que requieren un lubricante adecuado (con una viscosidad
específica) para cumplir sus funciones de manera correcta. También está el dimensionamiento de
bombas, donde se necesita conocer la viscosidad de un fluido para determinar la carga que se
debe suministrar. Por otro lado, conocer la viscosidad de los combustibles derivados del petróleo
es necesario para su almacenamiento, manejo y condiciones de operación.
Para explicar el concepto de viscosidad, se debe pensar en el fluido como un conjunto de capas
mono moleculares que interaccionan entre sí al momento de aplicar un esfuerzo cortante en la
superficie de este (Streeter, 1999). Esta fuerza provoca movimiento en la primera capa, a su vez
ésta mueve la capa adyacente y así sucesivamente. Cada capa tendrá un movimiento menor a la
de la anterior, lo que nos indica que existe una pérdida de energía entre capa y capa, esta pérdida
se debe al roce interno que existe. A este roce se le denomina como viscosidad dinámica, la cual
dependerá de las propiedades específicas del fluido.
4
2. Objetivos
5
3. Marco Teórico
3.1 Efecto de soluciones electrolíticas en la viscosidad
La estructura del agua se basa en que su molécula presenta un dipolo capaz de crear enlaces
puentes de hidrogeno entre sí, además los polos (negativo y positivo) pueden interactuar con
otros iones que sean añadidos al agua. A nivel molecular las moléculas de agua rodean al ión
creando una esfera alrededor de él, conocida como “hydration shell”. Existen iones que son
capaces de crear una segunda capa de moléculas de agua alrededor de ellos, es decir, favorecen al
orden de las moléculas de agua, estos iones son llamados “structure making”. Otros iones
entorpecen la estructura del agua, al no ser capaces de mantener una esfera de moléculas
alrededor de ellos, debido a que las moléculas de disolvente que rodean al ión tienden a ser
atraídas más fuertemente por otras moléculas de agua, esto crea un desorden en la estructura del
agua, estos iones son llamados “structure breaking” (Tongraar A., 1998).
En 1929, Jones & Dole estudiaron las soluciones acuosas diluidas con electrolitos fuertes,
determinaron que la viscosidad de una solución diluida puede ser descrita de la forma:
𝜈=𝐾𝑡 (𝐸𝑐. 3)
Donde K es la constante propia de cada viscosímetro indicada por el fabricante, esta constante
esta descrita por la ley de Hagen-Poiseuille y t es el tiempo teórico.
Sin embargo este cambio de presión afecta en la precisión de la medida. Como consecuencia de
esto, el tiempo medido (tg) es mayor que el obtenido por la ley de Hagen-Poiseuille. La diferencia
entre el tiempo medido y el tiempo teórico es la corrección de Hagenbach (ϑ) (Wilke J):
6
𝜗 = 𝑡𝑔 − 𝑡 (𝐸𝑐. 4)
𝜈 = 𝐾 (𝑡 − 𝜗) (𝐸𝑐. 5)
106 𝜋𝑔𝐷 4 𝐻 𝐸
𝜈= 𝑡− 2 (𝐸𝑐. 6)
128𝑉𝐿 𝑡
Donde:
g: Aceleración de gravedad (m/s2)
D: Diámetro del capilar (m)
L: Largo del capilar (m)
H: Distancia promedio entre el menisco superior e inferior (m)
V: Volumen que pasa a través del capilar (m3)
t: Tiempo (s)
E: Factor de la energía cinética (mm2/s)
106 𝜋𝑔𝐷 4 𝐻
𝐾= (𝐸𝑐. 7)
128𝑉𝐿
𝐸
𝜗= (𝐸𝑐. 8)
𝐾 𝑡2
7
4. Materiales y métodos
4.1 Equipo Experimental
- Viscosímetro capilar Cannon-Fenske.
- Balanza Analítica.
- Cronómetro.
- Fluidos Newtonianos.
- Sal (NaCl, KCl, MgCl2).
8
4.2 Procedimiento experimental
Nota: Para lograr una limpieza adecuada entre medición y medición, es necesario utilizar un
solvente completamente miscible con la muestra, seguido de un solvente volátil y la corriente de
aire.
Nota: Tal como se mencionó en el punto 4, este es el momento adecuado para aclimatar la
muestra a la temperatura deseada.
7. Antes de quitar el tapón, se debe preparar el cronometro, para medir los tiempos (en
segundos) que tarda el líquido en ir desde E hasta F, y desde F hasta I. Para lograr una
correcta medición se deben tomar los tiempos cuando el líquido afora cada marca. Si los
tiempos medidos son menores al mínimo especificado para el viscosímetro, se recomienda
usar un viscosímetro con un capilar más fino, para los demás casos se puede utilizar la
corrección de Hagenbach.
Nota: Es recomendable grabar este suceso para un análisis más detallado y exacto.
9
Nota: Si las viscosidades obtenidas por los tiempos transcurridos en los bulbos C y J son similares
el fluido es Newtoniano; si difieren de manera significante el fluido es no-Newtoniano (Manual
glass viscometer P10-0103).
Tal como se mencionó en el procedimiento experimental, los datos recabados serán los
tiempos (medidos en segundos) que tarda el fluido en ir desde la marca E a la F, y desde la F a la I.
Es decir, por cada muestra serán obtenidos dos tiempos independientes el uno del otro, con los
que será posible calcular dos viscosidades por pasada, sin embargo, para el estudio en cuestión la
viscosidad informada será el promedio de las viscosidades calculadas para el bulbo C (en medidas
sucesivas). Para que dos mediciones de viscosidad (determinadas por el bulbo C) sean aceptadas
no deben diferir por más de un 1.5% la una de la otra (ASTM D445), en caso de que esto ocurra
deben realizarse más mediciones. La importancia de la viscosidad obtenida en el bulbo J radica
que con ella puede discriminarse si el fluido es Newtoniano o no-Newtoniano (Manual glass
viscometer P10-0103), esto por simple comparación de los valores de viscosidad obtenidos para
los bulbos C y J.
Donde:
10
4.3.2 Cálculo de la viscosidad dinámica:
Para cada tiempo medido existe una incertidumbre (ԑ) que se traspasa en el cálculo de la
viscosidad de la siguiente forma:
𝜕𝜈 2
𝜀𝜈2 = ( ) ∗ 𝜀𝑡2
𝜕𝑡
𝐸
Dado que la viscosidad se expresa como: 𝜈 = 𝐾(𝑡 − 𝜗) = 𝐾𝑡 − 𝑡2
𝜕𝜈 2𝐸
( )=𝐾+ 3
𝜕𝑡 𝑡
A su vez, la masa medida durante la preparación de la solución posee una incertidumbre que se
traspasa al cálculo de densidad:
𝑚𝑠𝑜𝑙𝑢𝑡𝑜 + 𝑚𝑠𝑜𝑙𝑣𝑒𝑛𝑡𝑒
𝜌=
𝑉𝑜𝑙𝑢𝑚𝑒𝑛
𝜕𝜌 2
𝜀𝜌2 = ( ) ∗ 𝜀𝑚
2
𝜕𝑚
𝜕𝜌 1
( )=
𝜕𝑚 𝑉
Finalmente el para la viscosidad dinámica:
𝜂 = 𝜈 ⋅ 𝜌 ⋅ 10−3
Cuyo error relativo estará dado por:
𝜕𝜂 2 𝜕𝜂 2 𝜕𝜂 𝜕𝜂
𝜀𝜂2 = ( ) ∗ 𝜀𝜈2 + ( ) ∗ 𝜀𝜌2 Donde: ( )= 𝜌 ( )= 𝜈
𝜕𝜈 𝜕𝜌 𝜕𝜈 𝜕𝜌
11
5. Resultados y Discusión
5.1 Resultados
Tiempo aforo
Cantidad sal Cantidad
Normalidad ± 0.01 (s)
± 0.01 (g) agua (L)
E-F F-I
3.5 15.34 0.075 35.77 48.08
4.5 15.77 0.06 38.98 54.11
5.5 19.28 0.06 41.51 57.58
Tiempo aforo
Cantidad sal Cantidad
Normalidad ± 0.01 (s)
± 0.01 (g) agua (L)
E-F F-I
2.0 9.11 0.06 24.94 34.38
3.1 13.76 0.06 24.17 33.28
3.8 17.2 0.06 25.41 34.80
12
Las viscosidades cinemáticas y dinámicas resultantes son:
13
Efecto de la viscosidad cinemática
sobre soluciones de NaCl
1,6
Visc cinemática (mm^2/s)
1,5
1,4
1,3
1,2 Viscosidad cinemática
1,1 Agua
1
0,9
0,9 1,9 2,9 3,9 4,9 5,9
Normalidad (mol/L)
0,95
0,8
0,8 1,8 2,8 3,8
Normalidad (mol/L)
Gráfico 2. Medición en laboratorio para KCl a Normalidades cercanas a KClsat; Tamb: 20,5°C
14
Efecto de la viscosidad dinámica
sobre soluciones de NaCl
1,9
Visc dinámica (cP)
1,7
1,5
1,3
Viscosidad dinámica
1,1 Agua
0,9
0,9 1,9 2,9 3,9 4,9
Normalidad (mol/L)
1,05
1
Viscosidad dinámica
0,95 Agua
0,9
0,9 1,9 2,9 3,9 4,9
Normalidad (mol/L)
Gráfico 4. Medición en laboratorio para KCl a Normalidades cercanas a KClsat; Tamb: 20,5°C
15
5.2 Discusión
Los resultados obtenidos en las tablas 4, 5 y 8, indican que el NaCl es un soluto “structure
maker” o “Kosmotropes” porque la viscosidad de la solución aumenta con respecto a la
del agua a las mismas condiciones ambientales. Este tipo de solutos provocan un
ordenamiento molecular del agua, disminuyendo la entropía, y afectando el valor de la
constante B en la ecuación 1 (Lencka et al., 1998), que depende fuertemente de esta. La
constante B maximiza cuando la entropía disminuye, aumentando de esta manera la
viscosidad. En términos físicos, este aumento de la viscosidad se debe a que el mejor
ordenamiento molecular aumenta el roce entre capas adyacentes.
En cambio, el KCl puede ser clasificado como “structure breaker” o “Chaotropes”, porque
la viscosidad cinemática de la solución tiende a disminuir con respecto al agua pura (tablas
6,7 y 8), aumentando la fluidez del líquido, desestructurándola (aumentando la entropía)
y por lo tanto disminuyendo la viscosidad, es decir, genera valores negativos de la
constante B.
En específico, las sales usadas en este laboratorio no muestran grandes cambios, ya que
los iones involucrados son medianamente hidratados, y los efectos de los electrolitos se
intensifican de acuerdo a si son fuertemente hidratados (structure maker) o débilmente
hidratados (structure breaker).
En general la concentración de los iones en la solución afecta proporcionalmente en su
efecto como structure maker o structure breaker (Goldsack and Franchetto, 1976), como
es el caso del NaCl quien aumenta la viscosidad con la concentración, o como el KCl quien
disminuye la viscosidad (para soluciones diluidas). El caso de los structure breaker es
especial, ya que en una gráfica viscosidad contra concentración, la viscosidad dinámica en
soluciones diluidas alcanza un valor mínimo y luego aumenta a medida que la
concentración va aumentando, esto se debe porque inicialmente la presencia de estos
iones aumenta la fluidez, quien actúa opuestamente (en signo) a los valores de B, y a
medida que la concentración va aumentando la fluidez disminuye producto que el fluido
se vuelve más denso, convirtiendo a B en una cantidad positiva (Goldsack and Franchetto,
1976).
El comportamiento del agua como fluido Newtoniano, no se vio afectado por la adición de
las sales estudiadas. Esto se verifica por la similitud que hay en las viscosidades
determinadas en el bulbo C y en el bulbo J para cada concentración, donde los resultados
tan solo difieren por bajo un 4% entre ellos.
16
6 Conclusión
Comparando lo planteado con los resultados obtenidos, se puede observar que la teoría
fue correcta, en la solución de NaCl se aprecia que la viscosidad aumentó y está siguió
incrementando a medida que se aumentaba la concentración.
Por otro lado, para la solución de KCl ocurrió lo contrario, disminuyó la viscosidad
cinemática por efecto del comportamiento del ión potasio. Sin embargo, la concentración
de la solución de cloruro de potasio provoca un leve aumento de la viscosidad dinámica
respecto al solvente puro, producto que la densidad de la solución en esta situación toma
el papel determinante en el valor de la viscosidad.
Las soluciones se comportan como fluidos newtonianos, aun cuando se añade una
cantidad de electrolito a la solución.
17
7 Anexos
Para calcular la viscosidad dinámica es necesario conocer la densidad. Con los datos presentes
es posible calcularla como el cociente de la masa y el volumen.
𝑘𝑔
𝑚𝑎𝑔𝑢𝑎 = (𝑉𝑠𝑜𝑙 − 𝑉𝑁𝑎𝐶𝑙 ) × 998
𝑚3
0.01534 𝑘𝑔 𝑘𝑔
𝑚𝑎𝑔𝑢𝑎 = 0.000075 𝑚3 − ( ) × 998 3 = 0.06776 𝑘𝑔
𝑘𝑔 𝑚
2160 3
( 𝑚 )
𝑚 (0.01534 + 0.06776) 𝑘𝑔
𝜌= = = 1108.0 𝑘𝑔/𝑚3
𝑉 0.000075 𝑚3
Normalidad Masa sal (kg) Masa agua (kg) Volumen (m3) Densidad (kg/m3)
3.5 0.01534 0.06776 0.000075 1108.0
4.5 0.01577 0.05259 0.00006 1139.4
5.5 0.01928 0.05097 0.00006 1170.9
Normalidad Masa (kg) Masa agua (kg) Volumen (m3) Densidad (kg/m3)
2.0 0.00911 0.05567 0.00006 1079.7
3.1 0.01376 0.05352 0.00006 1121.4
3.8 0.0172 0.05193 0.00006 1152.2
18
7.1.2 Cálculo de viscosidad
La constante K del viscosímetro en este caso es 0.03644 mm2/s2. En el bulbo C (aforo E - F), el
tiempo que tarda en desplazarse el fluido es 35.77 s; y la corrección de Hagenbach para este
tiempo es aproximadamente 0.1112. Luego la viscosidad cinemática es:
𝜈 = 𝐾(𝑡 − 𝜗)
𝑚𝑚2 𝑚𝑚2
𝜈 = 0.03644 ( ) × (35.77 − 0.1112)(𝑠) = 1.299 ( )
𝑠2 𝑠
𝜂 = 𝜈 ∙ 𝜌 ∙ 10−3
𝑚𝑚2 𝑘𝑔
𝜂 = 1.299 ( ) × 1202.5 ( 3 ) × 10−3 = 1.563 𝑐𝑃
𝑠 𝑚
En el bulbo J (aforo F - I), la constante es 25% menor que en el bulbo C, entonces K es igual a
0.02733 mm2/s2. El tiempo en este caso es 48.08 s. Luego la viscosidad cinemática es:
𝑚𝑚2 𝑚𝑚2
𝜈 = 0.02733 ( ) × (48.08 − 0.1112)(𝑠) = 1.311 ( )
𝑠2 𝑠
𝑚𝑚2 𝑘𝑔
𝜂 = 1.311 ( ) × 1202.5 ( 3 ) × 10−3 = 1.577 𝑐𝑃
𝑠 𝑚
19
7.2 Anexo 2. Corrección de Hagenbach
𝐸 = 𝐾 𝜗 𝑡2
20
Para los valores de tiempo y corrección de Hagenbach conocidos en el anexo anterior:
Tabla 12. Factor de energía cinética de acuerdo al tiempo medido y a la corrección de Hagenbach
Conocido este valor, es posible extrapolar o interpolar el valor del factor de corrección de
Hagenbach para un capilar número 150 de la forma:
𝐸
𝜗=
𝐾 𝑡2
21
Referencias
Cannon Instrument Company, “Instructions for the use of The Cannon-Fenske Opaque (Reverse-
Flow)” Viscometer. https://www.cannoninstrument.com/Image/GetDocument/440?language=en
[Visitado abril de 2016].
Goldsack D., Franchetto R., “The viscosity of concentrated electrolyte solutions. I. Concentration
dependence at fixed temperature”, Can J. Chem. 55, 1062 (1977).
Jones G., Dole M., “The viscosity of aqueous solutions of strong electrolytes with special reference
to barium chloride”, 51(10), 2950-2964, 1929.
Lencka M., Anderko A., Sanders S. J. ,2 and Young R. D., “Modeling Viscosity of Multicomponent
Electrolyte”, International Journal of Thermophysics, Vol. 19, No. 2, 1998.
Marcus, Y., “Effect of Ions on the Structure of Water: Structure Making and breaking”, Chem. Rev.,
2009, 109 (3), pp 1346–1370.
Mott, R., Mecánica de Fluidos, 6° Ed., Pearson Education Prentice-Hall, México (2006).
Norma (ASTM) No. D445-15a “Standar Test Method for Kinematic Viscosity of Transparent and
Opaque Liquids” Título 71, Sección 1/97.
Norma (ASTM) No. D446-12 “Standar Specifications and Operating Instructions for Glass Capillary
Kinematic Viscometers” Título 71, Sección 2/95.
Sandoval, C., Caram, J., Salinas, J., The misleading simplicity of Stokes method for viscosity
measurement, Argentina (2009).
Streeter, V., Wylie, E.B., Bedford, K., Mecánica de Fluidos, 9° Ed., McGraw-Hill, Colombia (1997).
Tongraar A., Liedl K., Rode B., “Born-Oppenheimer ab Initio QM/MM Dynamics Simulations of Na+
and K+ in Water: From Structure Making to Structure Breaking Effects”, J. Phys. Chem. A., 102,
10340-10347, 1998.
Welty, J., Wicks, C., Wilson, R., Rorrer, G., Fundamentals of Momentum, Heat, and Mass Transfer,
5° Ed., John Wiley & Sons, Inc., Estados Unidos (2007).
Wilke J., Kryk H., Hartmann J., “Theory and Praxis of Capillary Viscometry”, Schott-Geräte,
Alemania.
22