Fenómenos de transporte Juan Esteban

UTN – Facultad Regional La Plata Zarza Galván

Trabajo Práctico Nº 2

Balances envolventes de Materia y

Cantidad de Movimiento para
Flujo en estado estacionario

Fecha de presentación: 26 de junio de 2021.

Fecha de entrega: 7 de julio de 2021.

Nombre y Apellido: Juan Esteban Zarza Galván

Número de alumno: 30181

Profesor adjunto: Ing. Gustavo Fabián Labianca.

Página -1-
 Fenómenos de transporte Juan Esteban
UTN – Facultad Regional La Plata Zarza Galván

Ejercicio Nº 1: Un fluido viscoso circula con régimen laminar por una rendija formada por dos
paredes planas separadas una distancia 2B. Realizar los balances de materia y cantidad de
movimiento para obtener las expresiones para la distribución de densidad de flujo de cantidad de
movimiento y de velocidad.

Coordenadas: Uso se coordenadas Rectangulares

Balance de Materia:
𝑉𝑒𝑙𝑜𝑐𝑖𝑑𝑎𝑑 𝑉𝑒𝑙𝑜𝑐𝑖𝑑𝑎𝑑 𝑉𝑒𝑙𝑜𝑐𝑖𝑑𝑎𝑑 𝑉𝑒𝑙𝑜𝑐𝑖𝑑𝑎𝑑 𝑉𝑒𝑙𝑜𝑐𝑖𝑑𝑎𝑑
[𝑑𝑒 𝑎𝑐𝑢𝑚𝑢𝑙𝑎𝑐𝑖ó𝑛] = [𝑑𝑒 𝑒𝑛𝑡𝑟𝑎𝑑𝑎 ] − [ 𝑑𝑒 𝑠𝑎𝑙𝑖𝑑𝑎 ] + [𝑑𝑒 𝑔𝑒𝑛𝑒𝑟𝑎𝑐𝑖ó𝑛] − [𝑑𝑒 𝑐𝑜𝑛𝑠𝑢𝑚𝑜 ]
𝑑𝑒 𝑚𝑎𝑡𝑒𝑟𝑖𝑎 𝑑𝑒 𝑚𝑎𝑡𝑒𝑟𝑖𝑎 𝑑𝑒 𝑚𝑎𝑡𝑒𝑟𝑖𝑎 𝑑𝑒 𝑚𝑎𝑡𝑒𝑟𝑖𝑎 𝑑𝑒 𝑚𝑎𝑡𝑒𝑟𝑖𝑎

Suposiciones:
- Estado estacionario: no hay acumulación de materia
- No hay consumo ni generación de materia
- Flujo laminar
- µ y ρ constantes (Isotérmico)
- Sin perturbaciones de entrada y salida
- Fluido newtoniano
𝑉𝑒𝑙𝑜𝑐𝑖𝑑𝑎𝑑 𝑉𝑒𝑙𝑜𝑐𝑖𝑑𝑎𝑑
0 = [𝑑𝑒 𝑒𝑛𝑡𝑟𝑎𝑑𝑎] − [ 𝑑𝑒 𝑠𝑎𝑙𝑖𝑑𝑎 ] → 𝜌∆𝑥∆𝑦𝑉𝑧 | − 𝜌∆𝑥∆𝑦𝑉𝑧 | = 0
𝑧 𝑧 + ∆𝑧
𝑑𝑒 𝑚𝑎𝑡𝑒𝑟𝑖𝑎 𝑑𝑒 𝑚𝑎𝑡𝑒𝑟𝑖𝑎

Se divide la expresión anterior por el volumen de control, y se toma el limite cuando ΔV=0
𝜌𝑉𝑧 𝜌𝑉𝑧
lim ( | − | )= 0
∆𝑥∆𝑦∆𝑧→0 ∆𝑧 𝑧 ∆𝑧 𝑧 + ∆𝑧
Recordando la definición de derivada:
𝜌𝑉𝑧 |𝑧 − 𝜌𝑉𝑧 |𝑧+∆𝑧 𝜌𝑉𝑧 |𝑧+∆𝑧 − 𝜌𝑉𝑧 |𝑧
lim ( ) = 0 → lim (− )=0
∆𝑥∆𝑦∆𝑧→0 ∆𝑧 ∆𝑥∆𝑦∆𝑧→0 ∆𝑧
𝑑𝑉𝑧
−𝜌 = 0 (𝑒𝑐 1)
𝑑𝑧
La ecuación (1) indica que la Vz no depende de la coordenada z, sino de x e y.

Balance de cantidad de movimiento (CM):

𝑉𝑒𝑙𝑜𝑐𝑖𝑑𝑎𝑑 𝑉𝑒𝑙𝑜𝑐𝑖𝑑𝑎𝑑 𝑉𝑒𝑙𝑜𝑐𝑖𝑑𝑎𝑑 𝑆𝑢𝑚𝑎𝑡𝑜𝑟𝑖𝑎 𝑑𝑒
[𝑑𝑒 𝑎𝑐𝑢𝑚𝑢𝑙𝑎𝑐𝑖ó𝑛] = [𝑑𝑒 𝑒𝑛𝑡𝑟𝑎𝑑𝑎] − [ 𝑑𝑒 𝑠𝑎𝑙𝑖𝑑𝑎 ] + [𝑓𝑢𝑒𝑟𝑧𝑎𝑠 𝑞𝑢𝑒 𝑎𝑐𝑡ú𝑎𝑛]
𝑑𝑒 𝐶𝑀 𝑑𝑒 𝐶𝑀 𝑑𝑒 𝐶𝑀 𝑠𝑜𝑏𝑟𝑒 𝑒𝑙 𝑠𝑖𝑠𝑡𝑒𝑚𝑎

Página -2-
 Fenómenos de transporte Juan Esteban
UTN – Facultad Regional La Plata Zarza Galván

Hay 2 tipos de fuerzas que actúan sobre el sistema: las normales y las tangenciales.
𝑆𝑢𝑚𝑎𝑡𝑜𝑟𝑖𝑎 𝑑𝑒 → 𝑁𝑜𝑟𝑚𝑎𝑙𝑒𝑠 → 𝑃𝑒𝑠𝑜
[𝑓𝑢𝑒𝑟𝑧𝑎𝑠 𝑞𝑢𝑒 𝑎𝑐𝑡ú𝑎𝑛]
→ 𝑇𝑎𝑛𝑔𝑒𝑛𝑐𝑖𝑎𝑙𝑒𝑠 → 𝑃𝑟𝑒𝑠𝑖ó𝑛
𝑠𝑜𝑏𝑟𝑒 𝑒𝑙 𝑠𝑖𝑠𝑡𝑒𝑚𝑎 →𝐸𝑠𝑓𝑢𝑒𝑟𝑧𝑜𝑠 𝑐𝑜𝑟𝑡𝑎𝑛𝑡𝑒𝑠

Velocidad de Cantidad de movimiento: 𝐶𝑚 = 𝑚̇ × 𝑉𝑧 = 𝑉𝑧 × (𝑉𝑧 × 𝜌 × 𝐴) = 𝜌∆𝑥∆𝑦𝑉𝑧 2

 Ingreso de CM= 𝜌∆𝑥∆𝑦𝑉𝑧 2 |𝑧
 Egreso de CM= 𝜌∆𝑥∆𝑦𝑉𝑧 2 |𝑧+∆𝑧
Sumatorias de fuerzas:
 Peso: 𝐹𝑔 = 𝑚 × 𝑔 = (𝜌 × 𝑉) × 𝑔 = 𝜌𝑔∆𝑥∆𝑦∆𝑧
 Diferencia de Presión: 𝐹𝑝 = 𝑝 × 𝐴 = 𝑝∆𝑥∆𝑦|𝑧 − 𝑝∆𝑥∆𝑦|𝑧+∆𝑧
 Esfuerzos cortantes: 𝐹𝑡 = 𝜏 × 𝐴 como solo existe velocidad en sentido de z, entonces los
esfuerzos de corte serán 𝜏𝑥𝑧 y 𝜏𝑦𝑧 (el término 𝜏𝑧𝑧 = 0). Sin embargo, se realiza la
consideración que W>>2B, por lo tanto 𝜏𝑦𝑧 = 0.
Todas estas expresiones se ingresan en la ecuación de Balance de CM, y luego se realiza el mismo
procedimiento que con la ecuación (1).
𝜌∆𝑥∆𝑦𝑉𝑧 2 | − 𝜌∆𝑥∆𝑦𝑉𝑧 2 | + 𝜌𝑔∆𝑥∆𝑦∆𝑧 + 𝑝∆𝑥∆𝑦 | − 𝑝∆𝑥∆𝑦 | + 𝜏𝑥𝑧 ∆𝑦∆𝑧 |
𝑧 𝑧 + ∆𝑧 𝑧 𝑧 + ∆𝑧 𝑥
− 𝜏𝑥𝑧 ∆𝑦∆𝑧 | =0
𝑥 + ∆𝑥
𝜌∆𝑥∆𝑦𝑉𝑧 2 |𝑧 − 𝜌∆𝑥∆𝑦𝑉𝑧 2 |𝑧+∆𝑧 + 𝜌𝑔∆𝑥∆𝑦∆𝑧 + 𝑝∆𝑥∆𝑦|𝑧 − 𝑝∆𝑥∆𝑦|𝑧+∆𝑧
lim (
∆𝑉→0 ∆𝑥∆𝑦∆𝑧
𝜏𝑥𝑧 ∆𝑦∆𝑧|𝑥 − 𝜏𝑥𝑧 ∆𝑦∆𝑧|𝑥+∆𝑥
+ )= 0
∆𝑥∆𝑦∆𝑧

𝜌𝑉𝑧 2 |𝑧+∆𝑧 − 𝜌𝑉𝑧 2 |𝑧 𝑝| − 𝑝|𝑧 𝜏𝑥𝑧 |𝑥+∆𝑥 − 𝜏𝑥𝑧 |𝑥

lim (− ) + 𝜌𝑔 + lim (− 𝑧+∆𝑧 ) + lim (− )= 0
∆𝑉→0 ∆𝑧 ∆𝑉→0 ∆𝑧 ∆𝑉→0 ∆𝑥
𝜌𝑑𝑉𝑧 2 𝑑𝑝 𝑑𝜏𝑥𝑧
− + 𝜌𝑔 − − =0
𝑑𝑧 𝑑𝑧 𝑑𝑥
El primer término se simplifica de acuerdo con la definición de la ecuación (1).
𝑑𝜏𝑥𝑧 𝑑𝑝 (𝑝0 − 𝑝𝐿 ) (𝑝0 − 𝑝𝐿 ) + (𝜌 𝑔)(𝐿 − 0) (𝑝0 − 𝑝𝐿 ) + (𝜌𝑔𝐿 − 𝜌𝑔0)
= 𝜌𝑔 − = 𝜌𝑔 + = =
𝑑𝑥 𝑑𝑧 (𝐿 − 0) (𝐿 − 0) (𝐿 − 0)
(𝑝0 − 𝜌𝑔0) − (𝑝𝐿 − 𝜌𝑔𝐿)
=
𝐿
Recordando las definiciones de presión hidrostática:
 𝒫0 = (𝑝0 − 𝜌𝑔0)
 𝒫𝐿 = (𝑝𝐿 − 𝜌𝑔𝐿)

𝑑 𝜏𝑥𝑧 (𝒫0 − 𝒫𝐿 )
=
𝑑𝑥 𝐿
(𝒫0 − 𝒫𝐿 )
∫ 𝑑 𝜏𝑥𝑧 = ∫ 𝑑𝑥
𝐿
(𝒫0 − 𝒫𝐿 )
𝜏𝑥𝑧 = ∙𝑥+𝐶
𝐿

Página -3-
 Fenómenos de transporte Juan Esteban
UTN – Facultad Regional La Plata Zarza Galván

Condiciones criticas
 CC1: cuando 𝑥 = 0; 𝜏𝑥𝑧 = 0 → 𝐶 = 0
(𝒫0 − 𝒫𝐿 )
𝜏𝑥𝑧 = 𝑥
𝐿
Ahora, de acuerdo con la ley de la viscosidad de Newton, hallo el perfil de velocidades:

𝑑𝑉𝑧 𝜏𝑥𝑧
𝜏𝑥𝑧 = −𝜇 → 𝑑𝑉𝑧 = 𝑑𝑥
𝑑𝑥 −𝜇

(𝒫0 − 𝒫𝐿 )
𝑥
∫ 𝑑𝑉𝑧 = ∫ 𝐿 𝑑𝑥
−𝜇

(𝒫0 − 𝒫𝐿 )𝑥 2
𝑉𝑧 = − +𝐶
2𝜇 𝐿
(𝒫0 −𝒫𝐿 )∙𝐵2
CC2: cuando 𝑥 = 𝐵; 𝑣𝑧 = 0 → 𝐶 = 2𝜇 𝐿

(𝒫0 − 𝒫𝐿 ) ∙ 𝑥 2 (𝒫0 − 𝒫𝐿 ) ∙ 𝐵 2
𝑉𝑧 = − +
2𝜇 𝐿 2𝜇 𝐿

(𝒫0 − 𝒫𝐿 )
𝑉𝑧 = × (𝐵 2 − 𝑥 2 )
2𝜇 𝐿

Página -4-
 Fenómenos de transporte Juan Esteban
UTN – Facultad Regional La Plata Zarza Galván

Ejercicio Nº 2: Flujo laminar en una película que desciende por el exterior de un tubo circular.
En una experiencia de absorción de gases, un fluido viscoso asciende por el interior de un pequeño
tubo circular, para descender después por la parte exterior del mismo. Aplicar un balance de cantidad
de movimiento a una envoltura de película de espesor r, tal como se indica en la figura. Obsérvese
que las flechas de entrada de cantidad de movimiento se toman siempre en la dirección r positiva al
efectuar el balance, aun cuando en este caso ocurre que la cantidad de movimiento fluye en la
dirección r negativa. Demostrar que la
distribución de velocidad en la
película descendiente (despreciando
los efectos finales) es:
Salida de cantidad de movimiento en la envoltura de espesor Δr

Distribución de velocidad en la película exterior.

r
Distribución de velocidad en el interior del tubo

Fg Fuerza de gravedad que actúa sobre el volumen 2rΔrL

R
a.R

Entrada de cantidad de movimiento en la envoltura de espesor Δr

Fuerzas existentes:

Coordenadas: Uso se coordenadas cilíndricas.

Balance de Materia:
𝑉𝑒𝑙𝑜𝑐𝑖𝑑𝑎𝑑 𝑉𝑒𝑙𝑜𝑐𝑖𝑑𝑎𝑑 𝑉𝑒𝑙𝑜𝑐𝑖𝑑𝑎𝑑 𝑉𝑒𝑙𝑜𝑐𝑖𝑑𝑎𝑑 𝑉𝑒𝑙𝑜𝑐𝑖𝑑𝑎𝑑
[𝑑𝑒 𝑎𝑐𝑢𝑚𝑢𝑙𝑎𝑐𝑖ó𝑛] = [𝑑𝑒 𝑒𝑛𝑡𝑟𝑎𝑑𝑎 ] − [ 𝑑𝑒 𝑠𝑎𝑙𝑖𝑑𝑎 ] + [𝑑𝑒 𝑔𝑒𝑛𝑒𝑟𝑎𝑐𝑖ó𝑛] − [𝑑𝑒 𝑐𝑜𝑛𝑠𝑢𝑚𝑜 ]
𝑑𝑒 𝑚𝑎𝑡𝑒𝑟𝑖𝑎 𝑑𝑒 𝑚𝑎𝑡𝑒𝑟𝑖𝑎 𝑑𝑒 𝑚𝑎𝑡𝑒𝑟𝑖𝑎 𝑑𝑒 𝑚𝑎𝑡𝑒𝑟𝑖𝑎 𝑑𝑒 𝑚𝑎𝑡𝑒𝑟𝑖𝑎

Página -5-
 Fenómenos de transporte Juan Esteban
UTN – Facultad Regional La Plata Zarza Galván

Suposiciones:
- Estado estacionario: no hay acumulación de materia
- No hay consumo ni generación de materia
- Flujo laminar
- µ y ρ constantes (Isotérmico)
- Sin perturbaciones de entrada y salida
- Fluido newtoniano
𝑉𝑒𝑙𝑜𝑐𝑖𝑑𝑎𝑑 𝑉𝑒𝑙𝑜𝑐𝑖𝑑𝑎𝑑
0 = [𝑑𝑒 𝑒𝑛𝑡𝑟𝑎𝑑𝑎] − [ 𝑑𝑒 𝑠𝑎𝑙𝑖𝑑𝑎 ] → 𝜌𝑉𝑧 2𝜋𝑟∆r⌋𝑧=𝑧 − 𝜌𝑉𝑧 2𝜋𝑟∆r⌋𝑧=𝑧+∆𝑧 = 0
𝑑𝑒 𝑚𝑎𝑡𝑒𝑟𝑖𝑎 𝑑𝑒 𝑚𝑎𝑡𝑒𝑟𝑖𝑎

Se divide la expresión anterior por el volumen de control (2𝜋𝑟∆r∆z), y se toma el limite cuando
ΔV=0
𝜌𝑉𝑧 𝜌𝑉𝑧
lim ( | − | )= 0
∆𝑉→0 ∆𝑧 𝑧 ∆𝑧 𝑧 + ∆𝑧
Recordando la definición de derivada:
𝜌𝑉𝑧 |𝑧 − 𝜌𝑉𝑧 |𝑧+∆𝑧 𝜌𝑉𝑧 |𝑧+∆𝑧 − 𝜌𝑉𝑧 |𝑧
lim ( ) = 0 → lim (− )=0
∆𝑥∆𝑦∆𝑧→0 ∆𝑧 ∆𝑥∆𝑦∆𝑧→0 ∆𝑧
𝑑𝑉𝑧
−𝜌 = 0 (𝑒𝑐 1)
𝑑𝑧
Balance de cantidad de movimiento (CM):
𝑉𝑒𝑙𝑜𝑐𝑖𝑑𝑎𝑑 𝑉𝑒𝑙𝑜𝑐𝑖𝑑𝑎𝑑 𝑉𝑒𝑙𝑜𝑐𝑖𝑑𝑎𝑑 𝑆𝑢𝑚𝑎𝑡𝑜𝑟𝑖𝑎 𝑑𝑒
[𝑑𝑒 𝑎𝑐𝑢𝑚𝑢𝑙𝑎𝑐𝑖ó𝑛] = [𝑑𝑒 𝑒𝑛𝑡𝑟𝑎𝑑𝑎] − [ 𝑑𝑒 𝑠𝑎𝑙𝑖𝑑𝑎 ] + [𝑓𝑢𝑒𝑟𝑧𝑎𝑠 𝑞𝑢𝑒 𝑎𝑐𝑡ú𝑎𝑛]
𝑑𝑒 𝐶𝑀 𝑑𝑒 𝐶𝑀 𝑑𝑒 𝐶𝑀 𝑠𝑜𝑏𝑟𝑒 𝑒𝑙 𝑠𝑖𝑠𝑡𝑒𝑚𝑎

2𝜋𝑟𝜌𝑉𝑧 2 ∆𝑟⌋𝑧=0 − 2𝜋𝑟𝜌𝑉𝑧 2 ∆𝑟⌋𝑧=𝐿 + 𝜏𝑟𝑧 2𝜋𝐿𝑅⌋𝑟=𝑟 − 𝜏𝑟𝑧 2𝜋𝐿𝑅⌋𝑟=𝑟+∆𝑟 + 2𝜋𝑟𝜌𝑔𝐿∆𝑟 = 0

𝜏𝑟𝑧 2𝜋𝐿𝑟⌋𝑟=𝑟 − 𝜏𝑟𝑧 2𝜋𝐿𝑟⌋𝑟=𝑟+∆𝑟 = −2𝜋𝑟𝜌𝑔𝐿∆𝑟

𝜏𝑟𝑧 𝑟⌋𝑟=𝑟 − 𝜏𝑟𝑧 𝑟⌋𝑟=𝑟+∆𝑟
= −𝑟𝜌𝑔
∆𝑟
Aplico el límite cuando Δr=0 para obtener una expresión diferencial:
𝜏𝑟𝑧 𝑟⌋𝑟=𝑟 − 𝜏𝑟𝑧 𝑟⌋𝑟=𝑟+∆𝑟 𝜏𝑟𝑧 𝑟⌋𝑟=𝑟+∆𝑟 − 𝜏𝑟𝑧 𝑟⌋𝑟=𝑟
lim = −𝑟𝜌𝑔 → lim − == −𝑟𝜌𝑔
∆𝑟→0 ∆𝑟 ∆𝑟→0 ∆𝑟
𝑑(𝑟𝜏𝑟𝑧 )
− = −𝑟𝜌𝑔 → 𝑑(𝑟𝜏𝑟𝑧 ) = 𝑟𝜌𝑔. 𝑑𝑟 → ∫ 𝑑(𝑟𝜏𝑟𝑧 ) = 𝜌𝑔 ∫ 𝑟𝑑𝑟
𝑑𝑟
𝜌𝑔𝑟 2 𝜌𝑔𝑟 𝐶
𝑟𝜏𝑟𝑧 = + 𝐶 → 𝜏𝑟𝑧 = +
2 2 𝑟
Condiciones de borde:
𝜌𝑔(𝑎𝑅) 𝐶 𝜌𝑔(𝑎𝑅)2 + 2𝐶 𝜌𝑔(𝑎𝑅)2
𝑟 = 𝑎𝑅 ; 𝜏𝑥𝑦 = 0 → 0 = + →0= →𝐶=−
2 (𝑎𝑅) 2𝑟 2
𝜌𝑔𝑟 𝜌𝑔(𝑎𝑅)2
𝜏𝑟𝑧 = −
2 2𝑟
Ahora, de acuerdo con la ley de la viscosidad de Newton, hallo el perfil de velocidades:

𝑑𝑉𝑧 𝜏𝑟𝑧
𝜏𝑟𝑧 = −𝜇 → 𝑑𝑉𝑧 = 𝑑𝑟
𝑑𝑟 −𝜇

𝜌𝑔(𝑎𝑅)2 𝜌𝑔𝑟
∫ 𝑑𝑉𝑧 = ∫ ( − ) 𝑑𝑟
2𝜇𝑟 2𝜇

Página -6-
 Fenómenos de transporte Juan Esteban
UTN – Facultad Regional La Plata Zarza Galván

𝜌𝑔(𝑎𝑅)2 𝜌𝑔𝑟 2
𝑉𝑧 = ln 𝑟 − + 𝐶2
2𝜇 4𝜇

Condiciones de borde:

𝜌𝑔(𝑎𝑅)2 𝜌𝑔𝑅 2 𝜌𝑔 𝑅2
𝑉𝑧 = 0 ; 𝑟 = 𝑅 → 0 = ln 𝑅 − + 𝐶2 → 𝐶2 = − × ((𝑎𝑅)2 ln 𝑅 − )
2𝜇 4𝜇 2𝜇 2

𝜌𝑔 𝑟2 𝜌𝑔 𝑅2
𝑉𝑧 = × ((𝑎𝑅)2 ln 𝑟 − ) + × ( − (𝑎𝑅)2 ln 𝑅)
2𝜇 2 2𝜇 2

𝜌𝑔 𝑟2 𝑅2
𝑉𝑧 = × ((𝑎𝑅)2 ln 𝑟 − − (𝑎𝑅)2 ln 𝑅 + )
2𝜇 2 2
2 2
𝜌𝑔 𝑟 𝑅
𝑉𝑧 = × ((𝑎𝑅)2 (ln 𝑟 − ln 𝑅) − + )
2𝜇 2 2
2 2
𝜌𝑔 𝑟 𝑟 𝑅
𝑉𝑧 = × (𝑎2 𝑅 2 ln ( ) − + )
2𝜇 𝑅 2 2
Saco factor común 𝑅 2 /2 para llegar a la expresión deseada:
𝜌𝑔𝑅 2 𝑟 𝑟2
𝑉𝑧 = × (2𝑎2 ln ( ) − 2 + 1)
4𝜇 𝑅 𝑅

𝜌𝑔𝑅2 𝑟 2 𝑟
𝑉𝑧 = × [1 − ( ) + 2𝑎2 ln ( )]
4𝜇 𝑅 𝑅
Ejercicio Nº 3: Un fluido viscoso se mueve entre 2 placas horizontales debido a un gradiente de
presión constante La placa inferior está en reposo y la superior se mueve con velocidad V0 en sentido
opuesto al flujo.

Calcular:

1. Distribución del perfil de velocidad. V0

2. A qué distancia de la placa inferior Vx = 0.
3. Velocidad media. P Y Z H PL
0
4. Punto donde la velocidad es máxima. X

W
5. El caudal de circulación.
6. ¿Cuánto debe valer V0, para que el caudal neto sea 0? L
7. Fuerza que se ejerce en la placa inferior.

Uso de Coordenadas rectangulares

1) Distribución del perfil de velocidad: Ecuación de continuidad
𝜕𝜌 𝜕(𝜌𝑉𝑥 ) 𝜕(𝜌𝑉𝑦 ) 𝜕(𝜌𝑉𝑧 )
+ + + =0
𝜕𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑧
Suposiciones:
 Estado estacionario
 Flujo laminar
 µ y ρ constantes (Isotérmico)
 Fluido newtoniano

Página -7-
 Fenómenos de transporte Juan Esteban
UTN – Facultad Regional La Plata Zarza Galván

Observaciones:
𝜕𝜌
1. El término ⁄𝜕𝑡 se anula ya que la densidad no varía con el tiempo.
2. Al no haber movimiento del fluido en las direcciones Z ni Y, sus respectivas derivadas se
anulan.
𝜕(𝜌𝑉𝑥 )
= 0 𝑒𝑐. (1)
𝜕𝑥
Balance de CM para coordenadas rectangulares sobre el eje x:
𝜕𝑉𝑥 𝜕𝑉𝑥 𝜕𝑉𝑥 𝜕𝑉𝑥 𝜕𝑝 𝜕𝜏𝑥𝑥 𝜕𝜏𝑦𝑥 𝜕𝜏𝑧𝑥
𝜌( + 𝑉𝑥 + 𝑉𝑦 + 𝑉𝑧 )=− −( + + ) + 𝜌𝑔𝑥
𝜕𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑥 𝜕𝑥 𝜕𝑦 𝜕𝑧
Observaciones:
𝜕𝑉
1. El término 𝜕𝑡𝑥 se anula ya que estamos en estado estacionario.
𝜕𝑉𝑥
2. De acuerdo con la ecuación (1), el término 𝑉𝑥 es cero.
𝜕𝑥
𝜕𝑉𝑥 𝜕𝑉𝑥
3. Nuevamente, al no haber movimiento del fluido en las direcciones Z ni Y, 𝑉𝑦 𝑦 𝑉𝑧 se
𝜕𝑦 𝜕𝑧
anulan.
𝜕𝜏𝑦𝑥
4. El único esfuerzo realizado es , por los que los demás son cero.
𝜕𝑦
5. El término 𝜌𝑔𝑥 también se desprecia ya que no hay variaciones significativas de las fuerzas
gravitatorias.

𝜕𝜏𝑦𝑥 𝜕𝑝 (𝑃𝐿 − 𝑃0 ) (𝑃𝐿 − 𝑃0 ) (𝑃𝐿 − 𝑃0 )

=− =− → 𝜕𝜏𝑦𝑥 = − 𝜕𝑦 → ∫ 𝜕𝜏𝑦𝑥 = ∫ − 𝜕𝑦
𝜕𝑦 𝜕𝑥 𝐿 𝐿 𝐿
(𝑃𝐿 − 𝑃0 )
𝜏𝑦𝑥 = − 𝑦 − 𝐶1
𝐿
Por la ley de Newton de la viscosidad:
𝑑𝑉𝑥 𝜏𝑦𝑥
𝜏𝑦𝑥 = −𝜇 → 𝑑𝑉𝑥 = 𝑑𝑦
𝑑𝑦 −𝜇

(𝑃𝐿 − 𝑃0 ) 𝐶1 (𝑃𝐿 − 𝑃0 ) 𝐶1
𝑑𝑉𝑥 = ( 𝑦 + ) 𝑑𝑦 → ∫ 𝑑𝑉𝑥 = ∫ ( 𝑦 + ) 𝑑𝑦
𝜇𝐿 𝜇 𝜇𝐿 𝜇
(𝑃𝐿 − 𝑃0 ) 2 𝐶1
𝑉𝑥 = 𝑦 + 𝑦 + 𝐶2 𝑒𝑐. (2)
2𝜇𝐿 𝜇
Condiciones de borde:
 𝑦 = 0 ; 𝑉𝑥 = 0 → 𝑎𝑙 𝑟𝑒𝑒𝑚𝑝𝑙𝑎𝑧𝑎𝑟𝑙𝑎𝑠 𝑒𝑛 (2), 𝑙𝑎 𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡𝑒 𝐶2 = 0
 𝑦 = 𝐻 ; 𝑉𝑥 = 𝑉𝑜
(𝑃𝐿 − 𝑃0 ) 2 𝐶1 𝜇𝑉0 (𝑃𝐿 − 𝑃0 )𝐻
𝑉0 = 𝐻 + 𝐻 → 𝐶1 = −
2𝜇𝐿 𝜇 𝐻 2𝐿
Reemplazando ambas constantes en 2, obtenemos el perfil de velocidades:

(𝑃𝐿 − 𝑃0 ) 2 𝑉0 (𝑃𝐿 − 𝑃0 )𝐻
𝑉𝑥 = 𝑦 + 𝑦− 𝑦
2𝜇𝐿 𝐻 2𝜇𝐿
2) A qué distancia de la placa inferior Vx = 0.
Como podemos observar, el perfil de velocidad es una ecuación cuadrática, por lo tanto tendrá 2
soluciones. Se toma factor común “y” se procede a resolver:

(𝑃𝐿 − 𝑃0 ) 𝑉0 (𝑃𝐿 − 𝑃0 )𝐻
0 = 𝑦[ 𝑦+ − ]
2𝜇𝐿 𝐻 2𝜇𝐿

Página -8-
 Fenómenos de transporte Juan Esteban
UTN – Facultad Regional La Plata Zarza Galván

𝑦1 = 0 → 𝑐𝑜𝑛𝑑𝑖𝑐𝑖ó𝑛 𝑑𝑒 𝑏𝑜𝑟𝑑𝑒
(𝑃𝐿 − 𝑃0 ) 𝑉0 (𝑃𝐿 − 𝑃0 )𝐻 (𝑃𝐿 − 𝑃0 )𝐻 𝑉0 2𝜇𝐿
𝑦2 + − = 0 → 𝑦2 = [ − ]×
2𝜇𝐿 𝐻 2𝜇𝐿 2𝜇𝐿 𝐻 (𝑃𝐿 − 𝑃0 )
2𝜇𝑉0 𝐿
𝑦2 = 𝐻 −
𝐻(𝑃𝐿 − 𝑃0 )
3) Velocidad media
𝑊 𝐻 𝑊 𝐻 (𝑃𝐿 − 𝑃0 ) 2 𝑉0 (𝑃𝐿 − 𝑃0 )𝐻
2𝜇𝐿 𝑦 + 𝐻 𝑦 −
∫0 ∫0 𝑣𝑥 𝑑𝑦 𝑑𝑧 ∫0 ∫0 𝑦 𝑑𝑦 𝑑𝑧
2𝜇𝐿
〈𝑉𝑥 〉 = 𝑊 𝐻 = 𝑊 𝐻
∫0 ∫0 𝑑𝑦 𝑑𝑧 ∫0 ∫0 𝑑𝑦 𝑑𝑧
(𝑃𝐿 − 𝑃0 ) 3 𝑉0 (𝑃𝐿 − 𝑃0 )𝐻 3
𝐻 + 𝐻 − (𝑃𝐿 − 𝑃0 ) 2 𝑉0 (𝑃𝐿 − 𝑃0 )𝐻 2
6𝜇𝐿 2 4𝜇𝐿
〈𝑉𝑥 〉 = = 𝐻 + −
𝐻 6𝜇𝐿 2 4𝜇𝐿
(𝑃𝐿 − 𝑃0 ) 2 𝑉0
〈𝑉𝑥 〉 = − 𝐻 +
12𝜇𝐿 2
4) Velocidad máxima:
Para encontrar la velocidad máxima, vamos a suponer que el término del esfuerzo de la ecuación de
newton es igual a cero.
(𝑃𝐿 − 𝑃0 ) (𝑃𝐿 − 𝑃0 ) 𝜇𝑉0 (𝑃𝐿 − 𝑃0 )𝐻
0=− 𝑦 − 𝐶1 = − 𝑦− +
𝐿 𝐿 𝐻 2𝐿
𝜇𝑉0 (𝑃𝐿 − 𝑃0 )𝐻 𝐿
𝑦=[ − ] × (− )
𝐻 2𝐿 (𝑃𝐿 − 𝑃0 )
𝐻 𝜇𝐿𝑉0
𝑦= −
2 𝐻(𝑃𝐿 − 𝑃0 )
5) Caudal de circulación:
𝑉̇ =< 𝑉𝑥 >× 𝑊 × 𝐻
(𝑃𝐿 − 𝑃0 ) 𝑉0
𝑉̇ = − 𝑊𝐻 3 + 𝑊𝐻
12𝜇𝐿 2
̇
6) V0 cuando 𝑉 = 0:
(𝑃𝐿 − 𝑃0 ) 𝑉0
0=− 𝑊𝐻 3 + 𝑊𝐻
12𝜇𝐿 2
(𝑃𝐿 − 𝑃0 ) 2
𝑉0 = 𝐻
6𝜇𝐿

7) Fuerza que se ejerce en la placa inferior:

𝑊 𝐿 𝑊 𝐿
𝑑𝑉𝑥
𝐹𝑥 = ∫ ∫ 𝜏𝑦𝑥 | 𝑑𝑥 𝑑𝑧 = ∫ ∫ 𝜇 | 𝑑𝑥 𝑑𝑧
0 0 𝑦=0 0 0 𝑑𝑦 𝑦 = 0
𝑊 𝐿 (𝑃
𝐿 − 𝑃0 )
𝐹𝑥 = ∫ ∫ 𝑦 + 𝐶1 𝑑𝑥 𝑑𝑧
0 0 𝐿

𝑊 𝐿 (𝑃
𝐿 − 𝑃0 ) 𝜇𝑉0 (𝑃𝐿 − 𝑃0 )𝐻
𝐹𝑥 = ∫ ∫ 𝑦+ − 𝑑𝑥 𝑑𝑧
0 0 𝐿 𝐻 2𝐿

𝜇𝑉0 (𝑃𝐿 − 𝑃0 )𝐻
𝐹𝑥 = 𝐿 × 𝑊 ( − )
𝐻 2𝐿

Página -9-
 Fenómenos de transporte Juan Esteban
UTN – Facultad Regional La Plata Zarza Galván

Ejercicio Nº 4: Dos fluidos inmiscibles e incompresibles circulan, debido a un gradiente de presión,

en la dirección z de una estrecha rendija horizontal de longitud L y anchura w. Las velocidades de los
fluidos están ajustadas de tal forma que una mitad de la rendija está llena de fluido I (la fase más
densa), y la otra mitad está ocupada por el fluido II (la fase menos densa). Se pide analizar la
distribución de velocidad y de densidad de flujo de cantidad de movimiento en este sistema.

( p0  pL )b     3  
 xz      
2.L    
x 1       
  
b 2       

( p0  pL )b  3     
 xz      
2.L    
Uso de Coordenadas rectangulares
1) Distribución del perfil de velocidad: Ecuación de continuidad
𝜕𝜌 𝜕(𝜌𝑉𝑥 ) 𝜕(𝜌𝑉𝑦 ) 𝜕(𝜌𝑉𝑧 )
+ + + =0
𝜕𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑧
Suposiciones:
 Estado estacionario
 Flujo laminar
 µ y ρ constantes (Isotérmico)
 Fluido newtoniano
Observaciones:
𝜕𝜌
6. El término ⁄𝜕𝑡 se anula ya que la densidad no varía con el tiempo.
7. Al no haber movimiento del fluido en las direcciones X ni Y, sus respectivas derivadas se
anulan.
𝜕(𝜌𝑉𝑧 )
= 0 𝑒𝑐. (1)
𝜕𝑧
Balance de CM para coordenadas rectangulares sobre el eje x:
𝜕𝑉𝑧 𝜕𝑉𝑧 𝜕𝑉𝑧 𝜕𝑉𝑧 𝜕𝑝 𝜕𝜏𝑥𝑧 𝜕𝜏𝑦𝑧 𝜕𝜏𝑧𝑧
𝜌( + 𝑉𝑥 + 𝑉𝑦 + 𝑉𝑧 )=− −( + + ) + 𝜌𝑔𝑧
𝜕𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑥 𝜕𝑥 𝜕𝑦 𝜕𝑧
Observaciones:
𝜕𝑉
3. El término 𝜕𝑡𝑧 se anula ya que estamos en estado estacionario.
𝜕𝑉𝑧
4. De acuerdo con la ecuación (1), el término 𝑉𝑧 es cero.
𝜕𝑧
𝜕𝑉𝑧 𝜕𝑉𝑧
8. Nuevamente, al no haber movimiento del fluido en las direcciones X ni Y, 𝑉𝑦 𝑦 𝑉𝑥 se
𝜕𝑦 𝜕𝑥
anulan.
𝜕𝜏
9. El único esfuerzo realizado es 𝜕𝑥𝑥𝑧, por los que los demás son cero.
10. El término 𝜌𝑔𝑧 también se desprecia ya que el flujo avanza en dirección perpendicular al
sentido de la gravedad.
𝜕𝜏𝑥𝑧 𝜕𝑝
=−
𝜕𝑥 𝜕𝑥
𝑑 𝜏𝑥𝑧 (𝑝0 − 𝑝𝐿 )
=
𝑑𝑥 𝐿
Al integrar las ecuaciones para las dos regiones se obtiene:

Página -10-
 Fenómenos de transporte Juan Esteban
UTN – Facultad Regional La Plata Zarza Galván

(𝑝0 − 𝑝𝐿 )
𝜏 𝐼 𝑥𝑧 =
𝑥 + 𝐶𝐼
𝐿
(𝑝0 − 𝑝𝐿 )
𝜏 𝐼𝐼 𝑥𝑧 = 𝑥 + 𝐶𝐼𝐼
𝐿
CC1: la densidad de flujo de cantidad de movimiento 𝜏𝑥𝑧 es constante en la interface liquido-liquido
CCL: 𝑥 = 0 ; 𝜏 𝐼 𝑥𝑧 = 𝜏 𝐼𝐼 𝑥𝑧
Lo anterior afirma que 𝐶𝐼 = 𝐶𝐼𝐼 ambas constantes se denominan 𝐶1

Cuando La ley de viscosidad de Newton se sustituye resulta:

𝑑𝑣𝑧 𝐼 (𝑝0 − 𝑝𝐿 )
−𝜇 𝐼 = 𝑥 + 𝐶1
𝑑𝑥 𝐿
𝑑𝑣𝑧 𝐼𝐼 (𝑝0 − 𝑝𝐿 )
−𝜇 𝐼𝐼 = 𝑥 + 𝐶1
𝑑𝑥 𝐿
Al integrarse queda:
(𝑝0 − 𝑝𝐿 ) 2 𝐶1
𝑣𝑧 𝐼 = − 𝑥 − 𝐼 𝑥 + 𝐶2𝐼
2 𝜇𝐼 𝐿 𝜇
(𝑝0 − 𝑝 𝐿 ) 𝐶 1
𝑣𝑧 𝐼𝐼 = − 𝑥 2 − 𝐼𝐼 𝑥 + 𝐶2𝐼𝐼
2 𝜇 𝐼𝐼 𝐿 𝜇
Las tres constantes de integración pueden determinarse a partir de las siguientes conducciones límite:
CL2: 𝑥 = 0 ; 𝑣 𝐼 𝑧 = 𝑣 𝐼𝐼 𝑧

CL3: 𝑥 = −𝑏 ; 𝑣 𝐼 𝑧 = 0

CL4: 𝑥 = 𝑏 ; 𝑣 𝐼𝐼 𝑧 = 0

Se obtienen tres ecuaciones simultáneas para las constantes de integración:

Por CL2: 𝐶2𝐼 = 𝐶2𝐼𝐼

(𝑝0 −𝑝𝐿 ) 𝐶1
Por CL3: 0 = − (−𝑏)2 + 𝑏 + 𝐶2𝐼
2 𝜇𝐼 𝐿 𝜇𝐼

(𝑝0 −𝑝𝐿 ) 𝐶1
Por CL4: 0 = − 𝑏2 − 𝑏 + 𝐶2𝐼𝐼
2 𝜇 𝐼𝐼 𝐿 𝜇 𝐼𝐼
A partir de esas tres ecuaciones se obtiene:
(𝑝0 − 𝑝𝐿 )𝑏 𝜇 𝐼 − 𝜇 𝐼𝐼
𝐶1 = − ( 𝐼 )
𝐿 𝜇 + 𝜇 𝐼𝐼

(𝑝0 − 𝑝𝐿 )𝑏 2 2𝜇 𝐼
𝐶2 𝐼 = + ( 𝐼 ) = 𝐶2𝐼𝐼
𝐿 𝜇 + 𝜇 𝐼𝐼
Los perfiles de velocidad y densidad de flujo de cantidad de movimiento son:

(𝑝0 − 𝑝𝐿 ) (𝑝0 − 𝑝𝐿 )𝑏 𝜇 𝐼 − 𝜇 𝐼𝐼
𝜏𝑥𝑧 = 𝑥− ( 𝐼 )
𝐿 𝐿 𝜇 + 𝜇 𝐼𝐼
(𝑝0 − 𝑝𝐿 )𝑏 𝑥 1 𝜇 𝐼 − 𝜇 𝐼𝐼
𝜏𝑥𝑧 = [( ) − ( 𝐼 )]
𝐿 𝑏 2 𝜇 + 𝜇 𝐼𝐼

Página -11-
 Fenómenos de transporte Juan Esteban
UTN – Facultad Regional La Plata Zarza Galván

(𝑝0 − 𝑝𝐿 )𝑏 𝜇 𝐼 − 𝜇 𝐼𝐼
(𝑝0 − 𝑝𝐿 ) 2 [− ( 𝐼 )] (𝑝0 − 𝑝𝐿 )𝑏 2 2𝜇 𝐼
𝐼 𝐿 𝜇 + 𝜇 𝐼𝐼
𝑣𝑧 = − 𝑥 − 𝑥+ ( 𝐼 )
2 𝜇𝐼 𝐿 𝜇𝐼 𝐿 𝜇 + 𝜇 𝐼𝐼
(𝑝0 − 𝑝𝐿 )𝑏 2 2𝜇 𝐼 𝜇 𝐼 − 𝜇 𝐼𝐼 𝑥 𝑥 2
𝑣𝑧 𝐼 = [( ) + ( ) ( ) − ( ) ]
2 𝜇𝐼 𝐿 𝜇 𝐼 + 𝜇 𝐼𝐼 𝜇 𝐼 + 𝜇 𝐼𝐼 𝑏 𝑏

(𝑝0 − 𝑝𝐿 )𝑏 𝜇 𝐼 − 𝜇 𝐼𝐼
(𝑝0 − 𝑝𝐿 ) 2 [− ( 𝐼 )] (𝑝0 − 𝑝𝐿 )𝑏2 2𝜇 𝐼
𝐿 𝜇 + 𝜇 𝐼𝐼
𝑣𝑧 𝐼𝐼 =− 𝑥 − 𝑥+ ( 𝐼 )
2 𝜇 𝐼𝐼 𝐿 𝜇 𝐼𝐼 𝐿 𝜇 + 𝜇 𝐼𝐼
𝐼𝐼
(𝑝0 − 𝑝𝐿 )𝑏 2 2𝜇 𝐼𝐼 𝜇 𝐼 − 𝜇 𝐼𝐼 𝑥 𝑥 2
𝑣𝑧 = [( ) + ( ) ( ) − ( ) ]
2 𝜇 𝐼𝐼 𝐿 𝜇 𝐼 + 𝜇 𝐼𝐼 𝜇 𝐼 + 𝜇 𝐼𝐼 𝑏 𝑏

Ejercicio Nº 5: Flujo en tubos concéntricos con movimiento axial del cilindro interior (varilla).
Considerar el sistema representado en la figura, en el que la varilla cilíndrica se mueve con velocidad
V. La varilla y el cilindro son coaxiales. Hallar la distribución de velocidades en estado estacionario
y la velocidad volumétrica de flujo. Este tipo de problemas se presenta en el recubrimiento de
alambres con barniz.

Cilindro exterior fijo

Vz
Fluido a presión P0
Varilla en movimiento R
r
Vvarilla k.R

Fluido a presión P0

Resolución: Uso de coordenadas CILINDRICAS

Ecuación de cantidad de movimiento respecto de z en coordenadas cilíndricas:

𝜕𝑣𝑧 𝜕𝑣𝑧 𝑣𝜃 𝜕𝑣𝑧 𝜕𝑣𝑧 𝜕𝑝 1𝜕 𝜕𝑣𝑧 1 𝜕 2 𝑣𝑧 𝜕 2 𝑣𝑧

𝜌( + 𝑣𝑟 + + 𝑣𝑟 )=− +𝜇[ (𝑟 )+ 2 + ] + 𝜌𝑔𝑧
𝜕𝑡 𝜕𝑟 𝑟 𝜕𝜃 𝜕𝑧 𝜕𝑧 𝑟 𝜕𝑟 𝜕𝑟 𝑟 𝜕𝜃 2 𝜕𝑧 2

𝑣𝑟 = 0; 𝑣𝜃 = 0 ; 𝑣𝑧 = 0

1𝜕 𝜕𝑣𝑧
0 = 𝜇[ (𝑟 )]
𝑟 𝜕𝑟 𝜕𝑟
1𝑑 𝑑𝑣𝑧
(𝑟 )=0
𝑟 𝑑𝑟 𝑑𝑟
𝑑 𝑑𝑣𝑧
(𝑟 )=0
𝑑𝑟 𝑑𝑟
𝑑𝑣𝑧
∫ 𝑑 (𝑟 ) = 0 ∫ 𝑑𝑟
𝑑𝑟

Página -12-
 Fenómenos de transporte Juan Esteban
UTN – Facultad Regional La Plata Zarza Galván

𝑑𝑣𝑧
𝑟 = 𝐶1
𝑑𝑟
𝐶1
𝑑𝑣𝑧 = 𝑑𝑟
𝑟
𝑑𝑟
∫ 𝑑𝑣𝑧 = 𝐶1 ∫
𝑟

𝑣𝑧 = 𝐶1 𝑙𝑛 𝑟 + 𝐶2 Ec (1)

Condiciones de borde:
-Sobre la varilla: 𝑉𝑧 = 𝑉𝑚𝑎𝑥 𝑦 𝑟 = 𝐾𝑅
-Sobre la cara interior del cilindro 𝑣𝑧 = 0 𝑦 r=R

𝑉𝑚𝑎𝑥 = 𝐶1 𝑙𝑛 𝐾𝑅 + 𝐶2 (2)
{
𝐶1 𝑙𝑛 𝑅 + 𝐶2 = 0 (3)

De (2): 𝐶2 = −𝐶1 𝑙𝑛 𝑅
Ahora sustituyo en (3)
𝑉𝑚𝑎𝑥 = 𝐶1 𝑙𝑛 𝐾𝑅 − 𝐶1 𝑙𝑛 𝑅

𝐾𝑅
𝑉𝑚𝑎𝑥 = 𝐶1 𝑙𝑛 ( 𝑅 ) 𝑉𝑚𝑎𝑥 = 𝐶1 𝑙𝑛 𝐾
𝑉𝑚𝑎𝑥
𝐶1 = ( )
𝑙𝑛𝐾
𝑉𝑚𝑎𝑥 𝑙𝑛𝑅
𝐶2 = − ( )
𝑙𝑛 𝐾
Con estas constantes hallo las velocidades:
𝑉𝑚𝑎𝑥 𝑉𝑚𝑎𝑥 𝑙𝑛𝑅
𝑉𝑧 = ( ) 𝑙𝑛 𝑟 − ( )
𝑙𝑛𝐾 𝑙𝑛 𝐾

𝑉𝑚𝑎𝑥 𝑉𝑚𝑎𝑥 𝑟
𝑉𝑧 = ( ) (𝑙𝑛 𝑟 − 𝑙𝑛 𝑅) = ( ) 𝑙𝑛 ( )
𝑙𝑛𝐾 𝑙𝑛 𝐾 𝑅
𝑟
𝑉𝑧 𝑙𝑛 ( )
= 𝑅
𝑉𝑚𝑎𝑥 𝑙𝑛 𝐾

Cálculo de velocidad media:

2𝜋 𝑅
∫0 ∫𝐾𝑅 𝑉𝑍 𝑑𝑟 𝑑𝜃
〈𝑉𝑧 〉 = 2𝜋 𝑅
∫0 ∫𝐾𝑅 𝑑𝑟 𝑑𝜃
𝑟
2𝜋 𝑅 𝑙𝑛 (𝑅 )
∫0 ∫𝐾𝑅 (𝑉𝑚 𝑙𝑛 𝐾 ) 𝑑𝑟 𝑑𝜃
〈𝑉𝑧 〉 = 2𝜋 𝑅
∫0 ∫𝐾𝑅 𝑑𝑟 𝑑𝜃

𝑉𝑚
[(𝑅 × ln 𝑅 − 𝑅) − (𝐾𝑅 ln 𝐾𝑅 − 𝐾𝑅)][𝑅 × ln 𝑅 − 𝐾𝑅 ln 𝐾𝑅]
〈𝑉𝑧 〉 = 𝐾
𝑙𝑛
𝑅 − 𝐾𝑅

Página -13-
 Fenómenos de transporte Juan Esteban
UTN – Facultad Regional La Plata Zarza Galván

Cálculo de caudal volumétrico:

𝑉̇ =< 𝑉𝑧 >× 𝐴

𝑉̇ =< 𝑉𝑧 >× 𝜋(𝑅 2 − 𝐾 2 𝑅 2 )

𝑉̇ =< 𝑉𝑧 >× 𝜋 × 𝑅 2 (1 − 𝐾 2 )

Ejercicio Nº 6: Calcular la constante de transformación de velocidad de rotación en rpm

(revoluciones por minuto) a velocidad de deformación en s-1 de un viscosímetro rotacional, sabiendo
que el radio interior del rotor es de 1,8415 cm y el radio exterior del cuerpo sólido (Bob) es de 1,7245
cm.
Cilindro exterior gira
a velocidad 

V= f(r)

r

Resolución: Utilización de coordenadas cilíndricas.

Suposiciones:
 Flujo incompresible
 Flujo estacionario en sentido de 𝜃
 Flujo laminar
 Flujo viscoso
 Flujo desarrollado
 No hay cambios significativos ni por presión ni por gravedad

Balance de materia – Ecuación de continuidad:

𝜕𝜌 1 𝜕 1 𝜕 𝜕
+ (𝜌𝑟𝑉𝑟 ) + (𝜌𝑉𝜃 ) + (𝜌𝑉𝑧 ) = 0
𝜕𝑡 𝑟 𝜕𝑟 𝑟 𝜕𝜃 𝜕𝑧
𝜕𝜌
 El término se desprecia ya que el fluido es incompresible.
𝜕𝑡

Página -14-
 Fenómenos de transporte Juan Esteban
UTN – Facultad Regional La Plata Zarza Galván

1 𝜕
 El término 𝑟 𝜕𝑟 (𝜌𝑟𝑉𝑟 ) se desprecia ya que no hay flujo de materia en el sentido de r.

1 𝜕
 El término 𝑟 𝜕𝜃 (𝜌𝑉𝜃 ) no se desprecia, ya que hay movimiento de flujo de materia en el
sentido de 𝜃.

𝜕
 El término (𝜌𝑉𝑧 ) se desprecia ya que no hay movimiento de flujo de materia en el sentido
𝜕𝑧
de z.

Por lo tanto:
1 𝜕
(𝜌𝑉𝜃 ) = 0
𝑟 𝜕𝜃

Balance de cantidad de movimiento

Sólo se realiza para 𝜃, ya que el resto de los términos fueron despreciados.

𝜕𝑉𝜃 𝜕𝑉𝜃 𝑉𝜃 𝜕𝑉𝜃 𝑉𝑟 𝑉𝜃 𝜕𝑉𝜃

𝜌( + 𝑉𝑟 + + + 𝑉𝑧 )
𝜕𝑡 𝜕𝑟 𝑟 𝜕𝜃 𝑟 𝜕𝑧
1 𝜕𝑝 𝜕 1𝜕 1 𝜕 2 𝑉𝜃 2 𝜕𝑉𝑟 𝜕 2 𝑉𝜃
=− +𝜇[ ( (𝑟𝑉𝜃 )) + 2 + + ] + 𝜌𝑔𝜃
𝑟 𝜕𝜃 𝜕𝑟 𝑟 𝜕𝑟 𝑟 𝜕𝜃 2 𝑟 2 𝜕𝜃 𝜕𝑧 2

 El término 𝜌 sale fuera ya que al considerar al fluido como incompresible, la densidad se hace
constante.

𝜕𝑉
 El término 𝜕𝑡𝜃 se desprecia ya que la variación de V en el tiempo a lo largo de 𝜃 es cero, y la
derivada es cero.

𝜕𝑉𝜃
 El término 𝑉𝑟 se desprecia ya que la derivada se hace cero.
𝜕𝑟

𝑉𝜃 𝜕𝑉𝜃
 El término ya que el balance de materia nos indica que vale cero.
𝑟 𝜕𝜃

𝑉𝑟 𝑉𝜃
 El término se desprecia ya que las componentes 𝑉𝑟 y 𝑉𝑧 de velocidad son cero.
𝑟

𝜕𝑉𝜃
 El término 𝑉𝑧 se desprecia ya que la derivada se hace cero.
𝜕𝑧

1 𝜕𝑝
 El término 𝑟 𝜕𝜃 se desprecia ya que no hay variación significativa.

 El término 𝜇 sale como constante porque así lo considera Navier-Stokes.

𝜕 1 𝜕
 El término 𝜕𝑟 (𝑟 𝜕𝑟 (𝑟𝑉𝜃 )) no se desprecia ya que hay esfuerzo en 𝜃 con respecto a r.

Página -15-
 Fenómenos de transporte Juan Esteban
UTN – Facultad Regional La Plata Zarza Galván

1 𝜕2 𝑉𝜃
 El término 2 se desprecia ya que no hay esfuerzo y, por lo tanto, la derivada se hace
𝑟 𝜕𝜃2
cero.

𝟐 𝝏𝑽𝒓
 El término 𝒓𝟐 se desprecia ya que no hay esfuerzo y, por lo tanto, la derivada se hace cero.
𝝏𝜽

𝝏𝟐 𝑽
 El término 𝝏𝒛𝟐𝜽 se desprecia ya que no hay esfuerzo y, por lo tanto, la derivada se hace cero.
 El término 𝝆𝒈𝜽 se desprecia ya que no hay efecto considerable de la gravedad sobre el fluido.
 Por lo tanto:

𝜕 1𝜕
0= ( (𝑟𝑉𝜃 ))
𝜕𝑟 𝑟 𝜕𝑟

1𝜕
0. 𝜕𝑟 = 𝜕 ( (𝑟𝑉𝜃 ))
𝑟 𝜕𝑟

1𝜕
∫ 0. 𝜕𝑟 = ∫ 𝜕 ( (𝑟𝑉𝜃 ))
𝑟 𝜕𝑟

1𝜕
(𝑟𝑉𝜃 ) = 𝐶1
𝑟 𝜕𝑟
𝜕(𝑟𝑉𝜃 ) = 𝐶1 . 𝑟. 𝜕𝑟

∫ 𝜕(𝑟𝑉𝜃 ) = ∫ 𝐶1 . 𝑟. 𝜕𝑟

𝐶1 𝑟 2
𝑟𝑉𝜃 = + 𝐶2
2
Si 𝑉𝜃 = 0

𝑟=𝐵

𝐶1 𝐵 2
0= + 𝐶2
2
𝐶1 𝐵 2
𝐶2 = −
2
𝐶1 𝑟 2 𝐶1 𝐵 2
𝑟𝑉𝜃 = −
2 2
𝐶1 2
𝑟𝑉𝜃 = (𝑟 − 𝐵 2 )
2
Si 𝑉𝜃 = Ω. 𝑅 y 𝑟 = 𝑅

Página -16-
 Fenómenos de transporte Juan Esteban
UTN – Facultad Regional La Plata Zarza Galván

𝐶1 2
Ω. 𝑅 2 = (𝑟 − 𝐵 2 )
2
2Ω. 𝑅 2
𝐶1 = 2
𝑅 − 𝐵2

2Ω. 𝑅 2
𝑟𝑉𝜃 = ( (𝑟 2 − 𝐵 2 ))
2(𝑅 2 − 𝐵 2 )

2Ω. 𝑅 2 𝑟 2 − 𝐵2
𝑉𝜃 = ( ( ))
2(𝑅 2 − 𝐵 2 ) 𝑟

Derivando con respecto a r:

𝜕𝑉𝜃 𝑅2 𝑟 2 − 𝐵2
= (( 2 ) ( )) Ω
𝜕𝑟 (𝑅 − 𝐵 2 ) 𝑟2

𝑟=𝐵

𝜕𝑉𝜃 2𝑅 2
= 2 Ω = 𝐾1 Ω
𝜕𝑟 (𝑅 − 𝐵 2 )

2𝑅 2
𝐾1 = 2
(𝑅 − 𝐵 2 )

2(1,8415 )2
𝐾1 =
((1,8415 )2 − (1,7245)2 )

2(1,8415 )2
𝐾1 Ω = [ ]
((1,8415 )2 − (1,7245)2 )

𝑟𝑎𝑑 1 𝑣𝑢𝑒𝑙𝑡𝑎 60 𝑠
𝐾1 Ω = 16.255721 × ×
𝑠 2𝜋 𝑟𝑎𝑑 1𝑚𝑖𝑛
𝑲𝟏 𝛀 = 𝟏𝟓𝟓, 𝟐𝟑𝟎𝟕 𝒓𝒑𝒎

Constante de transformación de velocidad en rpm

Página -17-