Notas Ing. Spinel
Notas Ing. Spinel
Notas Ing. Spinel
El material aquí presentado hace parte del material de las notas de clase de la asignatura Fluidos
elaborado por el ingeniero Jorge Spinel, son temas que no son fáciles de ubicar en los textos de mecánica
de fluidos, cuyo tratamiento considero es bastante ilustrativo para los estudiantes de la asignatura.
Conceptos previos y complementarios pueden ser consultados en libros de texto o la bibliografía
especializada. Las notas constituyen una guía para trabajar en otro tipo de sistemas, se destacan por las
demostraciones realizadas, no incluyen ejemplos numéricos.
En esta sección el desarrollo se restringe a fluidos newtonianos y presupone que el comportamiento del
sistema es isotérmico.
Un cilindro giratorio se encuentra de manera concéntrica dentro de otro cilindro, en el espacio entre los
dos se ubica un fluido newtoniano a una temperatura constante, la fuerza requerida para la rotación
indica la viscosidad del fluido.
El cilindro interno, articulado a un eje movido por un motor, se hace girar a una velocidad angular
constante (𝜔).
El uso del tacómetro y el potenciómetro permite la medición directa de la velocidad angular constante
(𝜔) y la potencia consumida por el cilindro interno (𝑃𝑟 ).
Si el espaciamiento entre cilindros es muy pequeño (𝑎≪𝑅) el perfil de velocidad en la dirección
radial entre 𝑅 y (𝑅+𝑎) se representa así:
𝑑𝑉 ∆𝑉 𝜔𝑅
= =
𝑑𝑟 ∆𝑟 𝑎
El esfuerzo cortante que se genera entre el fluido que viaja adherido a la pared del cilindro que gira (𝜏𝑅 )
y el fluido adyacente se puede considerar como el resultado relacionar una fuerza diferencial (𝑑𝐹)
actuando sobre un área de contacto diferencial (𝑑𝐴), representada por 𝑅𝑑𝜃𝑑𝑦, siendo 𝑅𝑑𝜃 el lado
(medida del arco en la localización 𝑅) y 𝑑𝑦 la altura de dicha área diferencial.
Si se aplica la Ley de Newton de viscosidad a la acción sobre el elemento de área descrita resulta:
𝑑𝐹 𝑑𝐹 𝜇 𝑑𝑉 𝜇 𝜔𝑅
𝜏𝑅 = = = =
𝑑𝐴 𝑅𝑑𝜃𝑑𝑦 𝑔𝐶 𝑑𝑟 𝑔𝐶 𝑎
El proceso de integración sobre toda el área lateral del cilindro interno debe conducir a la fuerza 𝐹(par
motor) necesaria para mantener su rotación; en consecuencia:
𝐻 2𝜋
𝜇𝜔𝑅 2 𝑑𝜃𝑑𝑦 2𝜋𝜇𝜔𝑅 2 𝐻
𝐹 =∫∫ =
𝑔𝐶 𝑎 𝑔𝐶 𝑎
0 0
De acuerdo con las leyes de los cuerpos en rotación, el producto del par motor (𝐹) por su brazo
(𝑅) provee el momento de torsión (𝑇) asociado, el cual se expresa así:
2𝜋𝜇𝜔𝑅 3 𝐻
𝑇=𝐹𝑅=
𝑔𝐶 𝑎
Finalmente, la asociación del momento de torsión, mediante la velocidad angular de rotación (𝜔), con la
potencia requerida (𝑃𝑟 ) para mantener el movimiento del cilindro interno brinda:
2𝜋𝜇𝜔2 𝑅 3 𝐻
𝑃𝑟 = 𝑇 𝜔 =
𝑔𝐶 𝑎
𝑃𝑟 𝑔𝐶 𝑎
𝜇=
2𝜋𝜔 2 𝑅 3 𝐻
La viscosidad queda dependiente de tres parámetros geométricos (𝑅, 𝐻 y 𝑎), propios del instrumento, de
un parámetro dimensional (𝑔𝐶 ) que tiene que ver con el sistema de unidades adoptado y de dos variables
de operación de muy sencilla medición, mediante tacómetro (𝜔) y el potenciómetro (𝑃𝑟 ).
Para obtener la viscosidad de un fluido newtoniano, se realiza la deformación por acción rotatoria de un
disco, de radio 𝑅, que se encuentra separado una distancia ℎ de otro estacionario y entre los cuales se
coloca el fluido, a una temperatura dada.
Considerando un elemento diferencial de área (𝑑𝐴) sobre la cara inferior del disco que gira:
𝑑𝐴 = 𝑟𝑑𝜃𝑑𝑟
El fluido que viaja con el disco superior en el entorno del área 𝑑𝐴 por causa de la condición de no
deslizamiento, posee una velocidad igual a 𝜔𝑟.
En consecuencia, si se asume que el espacio entre discos (ℎ) es suficientemente pequeño como para
considerarlo flujo de Coutte, la velocidad de deformación, en la localización 𝑟, vendrá dada por:
𝑑𝑉 ∆𝑉 𝜔𝑅
= =
𝑑𝑦 ∆𝑦 ℎ
Si se aplica la Ley de Newton de viscosidad en la localización (ℎ, 𝑟) resulta:
𝜇 𝑑𝑉 𝜇 𝜔𝑟
𝜏|𝑦=ℎ = =
𝑔𝐶 𝑑𝑦 𝑔𝐶 ℎ
Ahora, al relacionar el esfuerzo cortante en (ℎ, 𝑟) con el área de acción de dicho esfuerzo, se puede
estimar la magnitud de la fuerza diferencial que ejerce el fluido en contacto con el disco sobre fluido
adyacente y debajo de él, así:
𝜇 𝜔𝑟 2 𝑑𝜃𝑑𝑟
𝑑𝐹 |𝑦=ℎ = 𝜏|𝑦=ℎ 𝑑𝐴 = 𝜏|𝑦=ℎ 𝑟𝑑𝜃𝑑𝑟 =
𝑔𝐶 ℎ
A diferencia del viscosímetro de discos concéntricos, en donde cualquier área de acción diferencial se
encuentra siempre a la misma distancia radial 𝑅, en el de discos paralelos varía entre 0 y 𝑅.
Por la variación radial se opta, para relacionar la viscosidad con la potencia requerida en el movimiento
del disco superior, por evaluar el momento de torsión diferencial (𝑑𝑇) que genera la fuerza diferencial
𝑑𝐹 en la localización radial 𝑟, así:
𝜇 𝜔𝑟 3 𝑑𝜃𝑑𝑟
𝑑𝑇|𝑦=ℎ = 𝑑𝐹 |𝑦=ℎ 𝑟 =
𝑔𝐶 ℎ
La estimación del momento torsión (𝑇) que se requiere para rotar del disco superior exige la integración
variando el ángulo 𝜃 entre 0 y 2𝜋, y el radio entre 0 y 𝑅. Por lo tanto:
𝑅 2𝜋 𝑅 2𝜋
𝜇 𝜔𝑟 3 𝑑𝜃𝑑𝑟 𝜋𝜇𝜔𝑅 4
𝑇 = ∫ ∫ 𝑑𝑇|𝑦=ℎ = ∫ ∫ =
𝑔𝐶 ℎ 2𝑔𝐶 ℎ
0 0 0 0
Considerando ahora que la potencia requerida para mantener el giro del disco superior (𝑃𝑟 ) se expresa,
en términos de torsión, como 𝑇 𝜔 se llega a:
𝜋𝜇𝜔2 𝑅 4
𝑃𝑟 =
2𝑔𝐶 ℎ
La viscosidad del fluido, a la temperatura que se encuentra, queda determinada, entonces por:
2𝑔𝐶 ℎ𝑃𝑟
𝜇=
𝜋𝜔 2 𝑅 4
Un cojinete cónico gira alrededor de su eje central, en el espacio entre la pared móvil y la estacionaria se
ubica el fluido con viscosidad por determinar.
Parámetros geométricos
• Ángulo de la sección cónica (𝜃)
• Altura (𝐻)
• Espaciamiento entre la pared móvil y la estacionaria (𝑎)
• Medida en dirección perpendicular a la dirección (𝑠)
Cuando el espaciamiento entre las paredes móvil y estacionaria es suficientemente pequeño, el flujo de
Coutte tiene lugar, de modo que en cualquier dirección perpendicular a la coordenada 𝑠 (llámese 𝑠′), la
aplicación de la Ley de Newton de viscosidad, al radio 𝑟, brinda:
𝜇 𝑑𝑉 𝜇 ∆𝑉 𝜇 𝜔𝑟
𝜏|𝑝𝑎𝑟𝑒𝑑 = = =
móvil 𝑔𝐶 𝑑𝑠′ 𝑔𝐶 ∆𝑠′ 𝑔𝐶 𝑎
La coordenada vertical (𝑦), se puede relacionar con la coordenada radial (𝑟), mediante el ángulo de la
sección cónica así:
𝑟 = 𝑦 tan 𝜃
Sustituyendo se produce:
𝜇 𝜔𝑦 tan 𝜃
𝜏|𝑝𝑎𝑟𝑒𝑑 =
móvil 𝑔𝐶 𝑎
El cálculo del momento de torsión que se debe aplicar al sistema para mantener su giro pasa por la
determinación del momento de torsión diferencial asociado con la posición vertical, es función del
esfuerzo cortante asociado con un elemento diferencial de área 𝑑𝐴 en la envolvente cónica. Se cumple:
𝑑𝐹 𝑑𝐹
𝜏|𝑝𝑎𝑟𝑒𝑑 = =
móvil 𝑑𝐴 2𝜋𝑟𝑑𝑠
𝜇 2𝜋𝜔 tan 𝜃 2 𝑦 2 𝑑𝑦
𝑑𝐹 =
𝑔𝐶 𝑎 cos 𝜃
𝜇 2𝜋𝜔 tan 𝜃 3 𝑦 3 𝑑𝑦
𝑑𝑇 = 𝑟𝑑𝐹 = 𝑦 tan 𝜃 𝑑𝐹 =
𝑔𝐶 𝑎 cos 𝜃
𝜇 𝜋𝜔2 tan 𝜃 3 𝐻 4
𝑃𝑟 = 𝑇 ∗ 𝜔 =
𝑔𝐶 2𝑎 cos 𝜃
Al despejar la viscosidad se obtiene:
2𝑔𝐶 𝑎 cos 𝜃 𝑃𝑟
𝜇=
𝜋𝜔 2 tan 𝜃 3 𝐻 4
Micro–manómetros
Se recuerda al lector que un manómetro de tubo en U, está formado por un tubo de vidrio u otro material
transparente doblado en forma de U lleno parcialmente con un líquido de densidad conocida. Uno de sus
extremos se conecta a la zona donde quiere medirse la presión, y el otro se deja libre a la atmósfera.
Cuando este evaluando, inicie desde la presión conocida y trabaje hacia la presión desconocida.
A igual profundidad, la presión es constante en el mismo fluido.
Cuando se mueva hacia abajo del manómetro, las presiones suman.
Cuando se mueva hacia arriba del manómetro, las presiones restan.
Solo incluya la presión atmosférica sobre terminales abiertas
La figura representa un micro–manómetro con ensanchamientos notables en las ramas del tubo en U, el
cual utiliza, dos fluidos líquidos manométricos, inmiscibles entre sí y un fluido de sello (inmiscible con el
fluido con el que entra en contacto), que puede ser líquido o gaseoso.
Los niveles planos 0 y 1 son indicativos de las posiciones de equilibrio de los fluidos manométricos 2 y 3
cuando no existe presión diferencial.
Este manómetro se muestra en posición de desequilibrio cuando una pequeña presión diferencial se
trasmite a sus ramas desde los puntos 𝑎 y 𝑏 (en el caso presente ubicados al mismo nivel).
𝐷2 𝑑 2 ℎ3
𝜋 ∆ℎ = 𝜋
4 4 2
De donde:
𝑑 2 ℎ3
∆ℎ = ( )
𝐷 2
Al sustituir:
𝑑 2
𝑝𝑎 − 𝑝𝑏 = ℎ3 (𝛾3 − 𝛾2 ) + ( ) ℎ3 (𝛾2 − 𝛾1 )
𝐷
Al factorizar:
𝑑 2
𝑝𝑎 − 𝑝𝑏 = ℎ3 [(𝛾3 − 𝛾2 ) + ( ) (𝛾2 − 𝛾1 )]
𝐷
𝑑 2
ℎ𝑎 − ℎ𝑏 = ℎ3 [(𝐺3 − 𝐺2 ) + ( ) (𝐺2 − 𝐺1 )]
𝐷
Se puede observar que los términos entre corchetes son todos dependientes del sistema
manométrico, o de la relación geométrica de diámetros. Esto indica que para un manómetro y
unos fluidos manométricos dados, la presión diferencial solo depende de la lectura ℎ3 .
Cuando el fluido de sello es un gas 𝛾2 ≫ 𝛾1, entonces (𝛾2 − 𝛾1 ) ≈ 𝛾2 y por tanto la ecuación se
reduce a:
𝑑 2
𝑝𝑎 − 𝑝𝑏 = ℎ3 {𝛾3 − 𝛾2 [1 − ( ) ]}
𝐷
La cota cero (0), sobre la rama inclinada, indica la posición de presión diferencia nula y la cota 𝐿
(para el caso mostrado) señala la posición de la interfase fluido de sello – fluido manométrico,
cuando la presión diferencial no es nula.
Igualando los volúmenes de fluido manométrico desplazado en cada una de las ramas del
manómetro resulta:
𝑑 2
∆ℎ = ( ) 𝐿
𝐷
En la figura se aprecia que ℎ2 = 𝐿 sen 𝜃
Sustituyendo:
𝑑 2
𝑝𝑎 − 𝑝𝑏 = (𝛾2 − 𝛾1 ) [sen 𝜃 + ( ) ] 𝐿
𝐷
𝑑 2
𝑝𝑎 − 𝑝𝑏 = 𝛾2 [sen 𝜃 + ( ) ] 𝐿
𝐷
Fluidos en rotación
Ecuación Básica del Estado de Equilibrio Relativo:
𝜕𝑝 𝜕𝑝 𝜕𝑝
𝛻𝑝 = 𝑒𝑥 +
⃗⃗⃗ 𝑒𝑦 +
⃗⃗⃗⃗ 𝑒
⃗⃗⃗
𝜕𝑥 𝜕𝑦 𝜕𝑧 𝑧
𝛻𝑝 = 𝜌(𝑔 − 𝑎)
La naturaleza radial
La identidad, con respecto al eje vertical.
Ecuación básica del estado de equilibrio relativo, en coordenadas cilíndricas se expresa como:
𝜕𝑝 1 𝜕𝑝 𝜕𝑝
𝛻𝑝 = 𝑒𝑟 +
⃗⃗⃗ 𝑒𝜃 +
⃗⃗⃗⃗ 𝑒
⃗⃗⃗⃗
𝜕𝑟 𝑟 𝜕𝜃 𝜕𝑦 𝑦
𝜔 es la velocidad angular
Dirección 𝑟:
𝜕𝑝
= 𝜌𝜔2 𝑟
𝜕𝑟
Dirección 𝜃:
𝜕𝑝
=0
𝜕𝜃
Dirección 𝑦:
𝜕𝑝
= −𝜌𝑔 = −𝛾
𝜕𝑦
𝛻𝑝 = −𝜌(𝑔𝑒⃗⃗⃗⃗𝑦 − 𝜔2 𝑟𝑒⃗⃗⃗𝑟 )
𝜕𝑝 1 𝜕𝑝 𝜕𝑝
𝑑𝑝 = 𝑑𝑟 + 𝑑𝜃 + 𝑑𝑦
𝜕𝑟 𝑟 𝜕𝜃 𝜕𝑦
Si para el campo de presión se remplazan términos, separan variables e integra se obtiene la siguiente
ecuación:
𝜌𝜔2 2 2
𝑝 = 𝑝𝑟 + (𝑟 − 𝑟𝑟 ) − 𝛾(𝑦 − 𝑦𝑟 )
2
Ecuación útil para calcular presiones absolutas, manométricas y diferencias de presión en sistemas
conformados por tanques sometidos a rotación uniforme, especialmente en operaciones de
centrifugación.
El nivel alcanzado por el líquido, sobre la superficie libre en la pared del tanque es 𝑦1 y en el eje de
rotación es 𝑦2 . La superficie libre es una isobara porque soporta la presión atmosférica, por tanto la
diferencia de presión entre dos puntos ubicados sobre la superficie es nula.
𝜌𝜔2 2 2
∆𝑝1−2 = 0 = (𝑟1 − 𝑟2 ) − 𝛾(𝑦1 − 𝑦2 )
2
𝜔2 𝑟 2
𝑦= + 𝑦2
2𝑔
Ecuación que describe la superficie libre, es una parábola que tiene su vértice en 𝑦2 .
Como el volumen diferencial es 𝑟𝑑𝑟𝑑𝜃𝑑𝑦el volumen de líquido contenido se puede determinar por:
𝑅
𝑦 2𝜋
𝑉 líquido = ∫ ∫ ∫ 𝑟𝑑𝜃𝑑𝑦𝑑𝑟
0 0
en rotación 0
y la segunda
𝑅
𝑉 líquido = 2𝜋 ∫ 𝑦𝑟𝑑𝑟
en rotación 0
𝜔2 𝑅 2
𝑉 líquido = 𝜋𝑅 2 ( + 𝑦2 )
4𝑔
en rotación
𝑉 líquido = 𝑉 líquido
en rotación en reposo
Como
𝑉 líquido = 𝜋𝑅 2 𝑦0
en reposo
En consecuencia
𝜔2 𝑅 2
𝑦0 = + 𝑦2
4𝑔
𝜔2 𝑅 2
𝑦0 − 𝑦2 =
4𝑔
𝜔2 𝑅 2
𝑦1 − 𝑦2 =
2𝑔
La disminución de nivel que experimenta la superficie libre, con respecto al nivel de reposo (𝑦0 ) en 𝑟 = 0
es igual al incremento que sucede en 𝑟 = 𝑅.
Los resultados son válidos cuando se conserva el volumen de líquido, lo que implica que la altura del
tanque es mayor que 𝑦1 .
Si para un tanque con el mismo radio, la velocidad angular se incrementa, la localización del vértice de la
parábola que describe la superficie libre desciende proporcional con 𝜔2 .
Ejemplo
Determinar la presión (psia) en el punto A del tanque cilíndrico de la figura
La línea continua la superficie libre en rotación uniforme alrededor del eje vertical.
El punto A se ubica en el fondo del tanque en la pared interior, a esta misma profundidad se ubica el
punto B pero sobre la superficie libre en rotación.
Solución
Evaluar mediante la ecuación de la superficie libre, la localización del vértice de la parábola que está
describe.
𝜔2 𝑅 2
𝑦2 = 𝑦0 −
4𝑔
𝜔2 𝑟 2
𝑦= + 𝑦2
2𝑔
√−2𝑔𝑦2
𝑟𝐵 = = 0.388𝑚
𝜔
Aplicando la ecuación general para el campo de presión a los puntos A y B del sistema, se tiene:
𝜌𝜔2 2
𝑝𝐴 = 𝑝𝐵 + (𝑟𝐴 − 𝑟𝐵2 ) − 𝛾(𝑦𝐴 − 𝑦𝐵 )
2
𝑦𝐴 = 𝑦𝐵 = 0
𝑷𝑨 = 𝟏𝟗. 𝟔 𝒑𝒔𝒊𝒂