Dependencia de la temperatura de la viscosidad.

La viscosidad depende en gran medida de la temperatura. En los líquidos suele disminuir al aumentar la temperatura, mientras que, en la mayoría de los gases, la viscosidad aumenta al aumentar la temperatura. Este artículo analiza varios modelos de esta dependencia, que van desde cálculos rigurosos de primeros principios para gases monoatómicos hasta correlaciones empíricas para líquidos.

Comprender la dependencia de la viscosidad con la temperatura es importante para muchas aplicaciones, por ejemplo, diseñar lubricantes que funcionen bien en condiciones de temperatura variables (como en el motor de un automóvil), ya que el rendimiento de un lubricante depende en parte de su viscosidad. Los problemas de ingeniería de este tipo caen bajo el ámbito de la tribología .

Aquí la viscosidad dinámica se denota por y la viscosidad cinemática por . Las fórmulas dadas son válidas sólo para una escala de temperatura absoluta ; por lo tanto, a menos que se indique lo contrario, las temperaturas están en kelvins . ${\displaystyle \mu }$ ${\displaystyle \nu }$

Causas fisicas

La viscosidad en los gases surge de las moléculas que atraviesan capas de flujo y transfieren impulso entre capas. Esta transferencia de impulso puede considerarse como una fuerza de fricción entre capas de flujo. Dado que la transferencia de impulso es causada por el movimiento libre de las moléculas de gas entre colisiones, el aumento de la agitación térmica de las moléculas da como resultado una mayor viscosidad. Por tanto, la viscosidad gaseosa aumenta con la temperatura.

En los líquidos, las fuerzas viscosas son causadas por moléculas que ejercen fuerzas de atracción entre sí a través de capas de flujo. El aumento de la temperatura da como resultado una disminución de la viscosidad porque una temperatura más alta significa que las partículas tienen mayor energía térmica y son más fácilmente capaces de superar las fuerzas de atracción que las unen. Un ejemplo cotidiano de esta disminución de la viscosidad es el aceite de cocina que se mueve con mayor fluidez en una sartén caliente que en una fría.

gases

La teoría cinética de los gases permite un cálculo preciso de la variación de temperatura de la viscosidad de los gases. La base teórica de la teoría cinética viene dada por la ecuación de Boltzmann y la teoría de Chapman-Enskog , que permiten un modelado estadístico preciso de las trayectorias moleculares. En particular, dado un modelo para interacciones intermoleculares, se puede calcular con alta precisión la viscosidad de gases monoatómicos y otros gases simples (para gases más complejos, como los compuestos de moléculas polares , se deben introducir suposiciones adicionales que reducen la precisión de la teoría). ). ^[1]

Las predicciones de viscosidad para cuatro modelos moleculares se analizan a continuación. Las predicciones de los primeros tres modelos (esfera dura, ley potencial y Sutherland) pueden expresarse simplemente en términos de funciones elementales. El modelo de Lennard-Jones predice una dependencia -más complicada , pero es más preciso que los otros tres modelos y se utiliza ampliamente en la práctica de la ingeniería. ${\displaystyle T}$

Teoría cinética de esfera dura

Si se modelan las moléculas de gas como esferas duras elásticas (con masa y diámetro ), entonces la teoría cinética elemental predice que la viscosidad aumenta con la raíz cuadrada de la temperatura absoluta : ${\displaystyle m}$ ${\displaystyle \sigma }$ ${\displaystyle T}$

${\displaystyle \mu =1.016\cdot {\frac {5}{16\sigma ^{2}}}\left({\frac {k_{\rm {B}}mT}{\pi }}\right)^{1/2}}$

¿Dónde está la constante de Boltzmann ? Si bien se predice correctamente el aumento de la viscosidad gaseosa con la temperatura, la tendencia no es exacta; la viscosidad de los gases reales aumenta más rápidamente que esto. Captar la dependencia real requiere modelos más realistas de interacciones moleculares, en particular la inclusión de interacciones atractivas que están presentes en todos los gases reales. ^[2] ${\displaystyle k_{\text{B}}}$ ${\displaystyle T^{1/2}}$ ${\displaystyle T}$

Fuerza de ley de potencia

Una mejora modesta con respecto al modelo de esfera dura es una fuerza repulsiva de ley potencia inversa, donde la fuerza entre dos moléculas separadas por la distancia es proporcional a , donde es un parámetro empírico. ^[3] Este no es un modelo realista para gases del mundo real (excepto posiblemente a alta temperatura), pero proporciona una ilustración simple de cómo las interacciones intermoleculares cambiantes afectan la dependencia predicha de la viscosidad con la temperatura. En este caso, la teoría cinética predice un aumento de temperatura como , donde . Más precisamente, si es la viscosidad conocida a temperatura , entonces ${\displaystyle r}$ ${\displaystyle 1/r^{\nu }}$ ${\displaystyle \nu }$ ${\displaystyle T^{s}}$ ${\displaystyle s=(1/2)+2/(\nu -1)}$ ${\displaystyle \mu '}$ ${\displaystyle T'}$

${\displaystyle \mu =\mu '(T/T')^{s}}$

Tomando se recupera el resultado de esfera dura, . Para finito , correspondiente a una repulsión más suave, es mayor que , lo que resulta en un aumento más rápido de la viscosidad en comparación con el modelo de esfera dura. Al ajustarse a los datos experimentales de hidrógeno y helio se obtienen predicciones para y se muestran en la tabla. El modelo es modestamente preciso para estos dos gases, pero inexacto para otros gases. ${\displaystyle \nu \rightarrow \infty }$ ${\displaystyle s=1/2}$ ${\displaystyle \nu }$ ${\displaystyle s}$ ${\displaystyle 1/2}$ ${\displaystyle s}$ ${\displaystyle \nu }$

modelo sutherland

Otro modelo simple para la viscosidad gaseosa es el modelo de Sutherland, que agrega atracciones intermoleculares débiles al modelo de esfera dura. ^[4] Si las atracciones son pequeñas, pueden tratarse perturbativamente , lo que lleva a

${\displaystyle \mu ={\frac {5}{16\sigma ^{2}}}\left({\frac {k_{\text{B}}mT}{\pi }}\right)^{1/2}\cdot \left(1+{\frac {S}{T}}\right)^{-1}}$

donde , llamada constante de Sutherland, se puede expresar en términos de los parámetros de la fuerza de atracción intermolecular. De manera equivalente, si se conoce la viscosidad a la temperatura , entonces ${\displaystyle S}$ ${\displaystyle \mu '}$ ${\displaystyle T'}$

${\displaystyle \mu =\mu '\left({\frac {T}{T'}}\right)^{3/2}{\frac {T'+S}{T+S}}}$

Los valores obtenidos del ajuste de los datos experimentales se muestran en la siguiente tabla para varios gases. El modelo es modestamente preciso para una serie de gases ( nitrógeno , oxígeno , argón , aire y otros), pero inexacto para otros gases como el hidrógeno y el helio . En general, se ha argumentado que el modelo de Sutherland es en realidad un modelo pobre de interacciones intermoleculares y sólo es útil como una fórmula de interpolación simple para un conjunto restringido de gases en un rango restringido de temperaturas. ${\displaystyle S}$

Lennard Jones

En condiciones bastante generales en el modelo molecular, la predicción de la teoría cinética se puede escribir en la forma ${\displaystyle \mu }$

${\displaystyle \mu ={\frac {5}{16\pi ^{1/2}}}{\frac {(mk_{\text{B}}T)^{1/2}}{\sigma ^{2}\Omega (T)}}}$

donde se llama integral de colisión y es función de la temperatura y de los parámetros de la interacción intermolecular. ^[5] Está completamente determinado por la teoría cinética, expresándose en términos de integrales sobre trayectorias de colisión de pares de moléculas. En general, es una función complicada tanto de la temperatura como de los parámetros moleculares; Los modelos de ley potencial y Sutherland son inusuales porque pueden expresarse en términos de funciones elementales. ${\displaystyle \Omega }$ ${\displaystyle \Omega }$ ${\displaystyle \Omega }$

El modelo de Lennard-Jones supone un potencial de par intermolecular de la forma

${\displaystyle V(r)=4\epsilon \left[\left({\frac {\sigma }{r}}\right)^{12}-\left({\frac {\sigma }{r}}\right)^{6}\right]}$

donde y son parámetros y es la distancia que separa los centros de masa de las moléculas. Como tal, el modelo está diseñado para moléculas esféricamente simétricas. Sin embargo, se utiliza con frecuencia para moléculas no esféricamente simétricas siempre que no posean un gran momento dipolar . ^{[5] [6]} ${\displaystyle \epsilon }$ ${\displaystyle \sigma }$ ${\displaystyle r}$

La integral de colisión del modelo de Lennard-Jones no se puede expresar exactamente en términos de funciones elementales. Sin embargo, se puede calcular numéricamente y la concordancia con la experimentación es buena, no sólo para moléculas esféricamente simétricas como los gases nobles , sino también para muchos gases poliatómicos. ^[6] También se ha sugerido una forma aproximada de : ^[7] ${\displaystyle \Omega }$ ${\displaystyle \Omega }$

${\displaystyle \Omega (T)=1.16145\left(T^{*}\right)^{-0.14874}+0.52487e^{-0.77320T^{*}}+2.16178e^{-2.43787T^{*}}}$

dónde . Esta ecuación tiene una desviación promedio de sólo el 0,064 por ciento del rango . ${\displaystyle T^{*}\equiv k_{\text{B}}T/\epsilon }$ ${\displaystyle 0.3<T^{*}<100}$

En la siguiente tabla se muestran los valores y estimaciones a partir de datos experimentales para varios gases comunes. ${\displaystyle \sigma }$ ${\displaystyle \epsilon }$

Líquidos

A diferencia de los gases, no existe una teoría microscópica sistemática para la viscosidad de los líquidos. ^[9] Sin embargo, existen varios modelos empíricos que extrapolan una dependencia de la temperatura basándose en las viscosidades experimentales disponibles.

Exponencial de dos parámetros

Una correlación empírica simple y generalizada para la viscosidad del líquido es una exponencial de dos parámetros:

${\displaystyle \mu =Ae^{B/T}}$

Esta ecuación se propuso por primera vez en 1913 y se conoce comúnmente como ecuación de Andrade (llamada así en honor al físico británico Edward Andrade ). Describe con precisión muchos líquidos en un rango de temperaturas. Su forma puede ser motivada modelando el transporte de impulso a nivel molecular como un proceso de velocidad activada, ^[10] aunque los supuestos físicos subyacentes a tales modelos han sido cuestionados. ^[11]

La siguiente tabla proporciona valores estimados de y para líquidos representativos. En la literatura se pueden encontrar tablas completas de estos parámetros para cientos de líquidos. ^[12] ${\displaystyle A}$ ${\displaystyle B}$

Exponenciales de tres y cuatro parámetros

También se pueden encontrar exponenciales tabulados con parámetros adicionales, por ejemplo

${\displaystyle \mu =A\exp {\left({\frac {B}{T-C}}\right)}}$

y

${\displaystyle \mu =A\exp {\left({\frac {B}{T}}+CT+DT^{2}\right)}}$

Los valores representativos se dan en las tablas siguientes.

Modelos de viscosidad cinemática.

El efecto de la temperatura sobre la viscosidad cinemática también se ha descrito mediante varias ecuaciones empíricas. ^[15] ${\displaystyle \nu }$

La fórmula de Walther generalmente se escribe en la forma

${\displaystyle \log _{10}(\log _{10}(\nu +\lambda ))=A-B\,\log _{10}T}$

donde es una constante de desplazamiento y y son parámetros empíricos. En las especificaciones de los lubricantes normalmente sólo se especifican dos temperaturas, en cuyo caso normalmente se supone un valor estándar de = 0,7. ${\displaystyle \lambda }$ ${\displaystyle A}$ ${\displaystyle B}$ ${\displaystyle \lambda }$

El modelo de Wright tiene la forma

${\displaystyle \log _{10}(\log _{10}(\nu +\lambda +f(\nu )))=A-B\,\log _{10}T}$

donde se ha agregado a la fórmula de Walther una función adicional , a menudo un ajuste polinomial a datos experimentales. ${\displaystyle f(v)}$

El modelo de Seeton se basa en una curva que ajusta la dependencia de la viscosidad de muchos líquidos ( refrigerantes , hidrocarburos y lubricantes) versus la temperatura y se aplica en un amplio rango de temperatura y viscosidad:

${\displaystyle \ln \left({\ln \left({\nu +0.7+e^{-\nu }K_{0}\left({\nu +1.244067}\right)}\right)}\right)=A-B\ln T}$

donde es la temperatura absoluta en kelvins, es la viscosidad cinemática en centistokes , es la función de Bessel modificada de orden cero de segundo tipo y son parámetros empíricos específicos de cada líquido. ${\displaystyle T}$ ${\displaystyle \nu }$ ${\displaystyle K_{0}}$ ${\displaystyle A}$ ${\displaystyle B}$

Para la viscosidad del metal líquido en función de la temperatura, Seeton propuso:

${\displaystyle \ln \left({\ln \left({\nu +0.7+e^{-\nu }K_{0}\left({\nu +1.244067}\right)}\right)}\right)=A-{B \over T}}$

Ver también

• Índice de viscosidad
• tribología
• Fenómenos de transporte
• Modelado molecular
• Fuerza intermolecular
• Campo de fuerza (química)
• método de retroceso

Notas

1. Chapman y Cowling (1970)
2. Chapman y Cowling (1970), págs. 97–98, 226–230
3. Chapman y Cowling (1970), págs.
4. Chapman y Cowling (1970), págs. 232-234
5. Reid, Prausnitz y Poling (1987), págs. 391–392
6. Bird, Stewart y Lightfoot (2007), págs.
7. Neufeld, Jansen y Aziz (1972)
8. Bird, Stewart y Lightfoot (2007), págs. 864–865
9. Reid, Prausnitz y Poling (1987), pág. 433
10. Bird, Stewart y Lightfoot (2007), págs. 29-31
11. Hildebrando (1977)
12. Véase Viswanath y Natarajan (1989); Viswanath y otros (2007); Reid, Prausnitz y Poling (1987); y referencias en el mismo
13. Reid, Prausnitz y Poling (1987)
14. Viswanath y Natarajan (1989)
15. Seeton (2006)

Referencias

• Pájaro, R. Byron; Stewart, Warren E.; Lightfoot, Edwin N. (2007), Fenómenos del transporte (2ª ed.), John Wiley & Sons, Inc., ISBN 978-0-470-11539-8.
• Chapman, Sídney; Cowling, TG (1970), La teoría matemática de los gases no uniformes (3.ª ed.), Cambridge University Press
• Hildebrand, Joel Henry (1977), Viscosidad y difusividad: un tratamiento predictivo , John Wiley & Sons, Inc., ISBN 978-0-471-03072-0
• Neufeld, Philip D.; Janzen, AR; Aziz, RA (1972). "Ecuaciones empíricas para calcular 16 de las integrales de colisión de transporte Ω (l, s) * para el potencial de Lennard-Jones (12-6)". La Revista de Física Química . 57 (3): 1100-1102. Código Bib :1972JChPh..57.1100N. doi :10.1063/1.1678363. ISSN  0021-9606.
• Reid, Robert C.; Prausnitz, John M.; Poling, Bruce E. (1987), Las propiedades de los gases y líquidos , McGraw-Hill Book Company, ISBN 0-07-051799-1
• Seeton, Christopher J. (2006), "Correlación viscosidad-temperatura para líquidos", Tribology Letters , 22 : 67–78, doi :10.1007/s11249-006-9071-2, S2CID  97897211
• Viswanath, DS; Natarajan, G. (1989). Libro de datos sobre la viscosidad de los líquidos . Corporación Editorial del Hemisferio. ISBN 0-89116-778-1.