1.5.

3 Balance de cantidad de movimiento en régimen estacionario en un sistema de

coordenadas cilíndricas.
Se aplica el balance de cantidad de movimiento al flujo que se desarrolla dentro de un tubo en
régimen laminar ya sea con fines de transportar un fluido de proceso o que se lleve a cabo una
reacción química, de tal forma que el tubo se comporta como un reactor químico. Los esquemas
explicativos se muestran a continuación.
Figura 1.10 Esquema para la aplicación del balance de cantidad de movimiento en coordenadas
cilíndricas

Fuente: Adaptado de: (Gaskell, D., 2013, pag.54)

Figura 1.11: Corte longitudinal del flujo vertical descendente en un tubo. (a) Esquema del
elemento diferencial de volumen, (b) Perfil de velocidad y esfuerzo de corte

Fuente: Adaptado de (Bird, et al. 2006)

En a) se observa una porción de tubo de longitud L, donde se practica el balance de cantidad de
movimiento; en b) se observa el elemento de volumen de control, mostrando las principales
contribuciones del balance de cantidad de movimiento; en c) se muestra un corte longitudinal que
permite observar el volumen de control de espesor r y d) permite observar los perfiles de
densidad de flujo viscoso de transporte de cantidad de movimiento y velocidad.
El balance de cantidad de movimiento en régimen estacionario viene dado por la ecuación (1.28);
los términos que intervienen como entrada y salida al y del volumen de control son los términos
de transporte advectivo y difusivo, las fuerzas superficiales o de presión y la fuerza debido al
campo gravitacional se considera como termino de entrada, pues actúa en l dirección del flujo.

A continuación, se muestran los términos que intervienen en el balance de cantidad de

movimiento aplicado al flujo vertical descendente
Tabla 1.4: Términos que componen el balance de cantidad de movimiento en coordenadas
cilindricas
Entrada de cantidad de
movimiento por transporte  rz r  2 rL
Salida de cantidad de
movimiento por  rz r +r  2 (r + r ) L
molecular transporte molecular
Entrada de cantidad de
movimiento por transporte ( 2 r rVz ) ( Vz ) z =0 Salida de cantidad de
movimiento por transporte ( 2 r rVz )( Vz ) z = L
advectivo advectivo

2 r rp
Fuerza de presión a la Fuerza de presión a la
entrada en la superficie 2 r rp z =0 salida de la superficie z=L
anular anular
Fuerza gravitacional que actúa en la envoltura anular: 2 r rL g
Aplicando la ecuación (1.28), se tiene:

 rz r  2 rL −  rz r +r  2 ( r + r ) L + ( 2 r rVz ) ( VZ ) z =0 − - ( 2 r rVz )( Vz ) z= L +

2 r rP0 z =0 − 2 r rP z=L + 2 r rL g = 0 (1.40)

En la ecuación (1.28) todos los términos están expresadas en unidades de fuerza; si se utiliza el
Sistema Internacional de Unidades, estarán expresadas en Newton (N).

2 rL y 2 ( r + r ) L es el área lateral, es decir el área de rozamiento del fluido; 2 r r es

el área de flujo del elemento diferencial; 2 rrL es el volumen del elemento diferencial. Cabe
mencionar la fuerza debido al flujo advectivo: ( 2 r rVz ) ( VZ ) z =0 y ( 2 r rVz )( Vz ) z= L

numéricamente son iguales, puesto que no hay acumulación ni cambio de diámetro, por tanto,
su diferencia es nula; en consecuencia, si se aplica el limite correspondiente a la ecuación (1.49)
reducida, se obtiene:

2 ( r + r ) L rz –2 rL rz 2 r r 2 r rL g

. ( p0 – pL ) –
r +r
=0
r
lim – (1.41)
r →0 2 Lr 2 Lr 2 Lr
La ecuación (1.41) conduce a:

 r rz – r rz r   ( p0 – pL )
r +r 
lim  = + g r (1.42)
r →0
 r   L 
Por otro lado:
P0 = p0 −  g (0) PL = pL −  gL , en consecuencia, la ecuación (1.42) se transforma en:

1d (P – P )
( r rz ) = 0 L (1.43)
r dr L
La ecuación (1.36) representa la ecuación de conservación de cantidad de movimiento (Glasgow,
L., 2010; pp.24-25; )
Soluciones analíticas de la ecuación de balance de cantidad de movimiento
Muchos de los problemas de balance de cantidad de movimiento aplicado a fluidos newtonianos
en sistema de coordenadas cilíndricas pueden ser resueltos por las técnicas clásicas de resolución
de ecuaciones diferenciales ordinarias de variables separables. En este item aplicaremos estas
técnicas a los diversos ejercicios que se dan a continuación.
Problema 1.22
Perfil de velocidad de un fluido newtoniano en un tubo
Considerando que el fluido que se desplaza por el tubo es un fluido newtoniano, resuelva la
ecuación diferencial dada por la ecuación (1.43) por medio de un procedimiento analítico y
encuentre los siguientes parámetros de flujo:
a) Ecuación de distribución de velocidades
b) Ecuación de distribución de la densidad de flujo viscoso de transporte de cantidad de
movimiento
c) Flujo volumétrico
d) Velocidad media
e) Velocidad máxima
f) Relación entre velocidad media y velocidad máxima
g) Fuerza de rozamiento en la pared.
h) Punto donde la velocidad media es igual a la velocidad puntual
Solución
a) Integrando la ecuación (1.36) se obtiene:

 rz =
( P0 – PL ) r + C1'
(1)
2L r
En la ecuación (1) las presiones están escritas en mayúsculas pues incluye tanto la presión
estática como la presión dinámica.
Dado que el fluido es newtoniano se tiene:  rz = –  dVz / dr ; esta ecuación se iguala con la
ecuación (1) y se tiene:

dVz ( P0 – PL ) C1'
– = r+ (2)
dr 2L r
Cuya integración conduce a:

Vz = –
( P0 – PL ) r 2 + C Ln r + C
(3)
4 L
1 2
En la ecuación (1.55) debemos encontrar las constantes de integración, usando las siguientes
condiciones de frontera.

dVz
r =0 o  rz = 0 (a)
dr

r = R Vz = 0 (b)

Si aplicamos la condición (a) en (2), se tendría una indeterminación; sin embargo, este es un
problema real y su solución tiene que estar acotada en el campo real; para levantar dicha
indeterminación, la condición necesaria es que C1 sea igual a cero. Ahora usamos la condición (b)
en (3) y se tiene:

C2 =
( P0 – PL ) R2
(4)
4 L
Luego reemplazando esta constante en (3) se obtiene:

Vz = –
( P0 – PL ) r 2 + ( P0 – PL ) R2 = ( P0 – PL ) R2 1 − (r / R)2 
4 L 4 L 4 L   (5)

La ecuación (5) es la ecuación de distribución de velocidades.

b) La ecuación de distribución de esfuerzos o de la densidad de flujo viscoso de transporte
de cantidad de movimiento se obtiene a partir de (1), es decir

 rz =
( P0 – PL ) r
(6)
2L

Puesto que si C1 = 0 entonces C1 = 0 ; el perfil de esfuerzos resulta ser lineal.

c) El caudal o flujo volumétrico se obtiene integrando el perfil de velocidades

Q=
 A
Vz (r )dA (7)

2 R
( P0 – PL ) R 2 1–  r 2  rdrd
Q=
0 0
4 L


  
 R  
(8)

La integración de (1.60) permite obtener:

 R 4 ( P0 – PL )
Q= (9)
8 L
La ecuación (9) se conoce como la ecuación de Hagen-Poiseville y se utiliza siempre que se
conozca la diferencia de presiones entre dos puntos.
d) La velocidad media se obtiene dividiendo el caudal entre el área de flujo
  Vz  =
Q
=
Q
=
( P0 – PL ) R 2 (10)
2
8 L
 
R
A rdrd
0 0

e) La velocidad máxima se obtiene haciendo r = 0 en (5), es decir:

Vz max =
( P0 – PL ) R2
(11)
4 L
f) La relación entre la velocidad media y la velocidad máxima se obtiene dividiendo la
ecuación (10) entre la ecuación (11)

 Vz  1
= (12)
Vz max 2
g) El perfil de velocidad en términos de la velocidad media

VZ = 2  VZ  R2 (1 − (r / R)2 ) (13)

h) La fuerza de rozamiento en la pared

( P0 – PL ) R =  R2
Fz = 2 RL rz r = R = 2 RL ( P0 – PL ) (14)
2L
i) Finalmente podemos obtener el punto donde la velocidad puntual es igual a la velocidad
media:

( P0 – PL ) R2 (1 − (r / R)2 ) = ( P0 – PL ) R2
(15)
4 L 8 L

2
Resolviendo se obtiene: r = R
2
Figura P_1.22 Parámetros de flujo expresados en forma grafica (a) Variacion de la velocidad en
función al radio del tubo; (b) Variacion del esfuerzo de corte en función al radio del tubo; (c)
Velocidad media del fluido en el tubo.
 Fuente: Adaptado de (Gaskell, D., 2013, pag. 56)

Problema 1.23
Flujo a través de un capilar inclinado
Determine el flujo másico en el tubo de flujo capilar de longitud L y de radio R de la figura. El
fluido que circula por el tubo capilar es agua y como fluido nanométrico se utiliza tetracloruro de
carbono. La Figura P_1.21 representa el sistema de flujo descrito esquematizada para la
realización del balance de cantidad de movimiento.
Figura P_1.23A: Esquema del flujo en un tubo capilar provisto de un manómetro diferencial

Fuente:

Adaptado de (Gaskell, D., 2013, pag. 65; Bird et al., 1993)

Solución
Desarrollamos el balance de presiones para encontrar su diferencia entre los puntos de medición

pA = p1

p1 = p2− a gh2
p2 = p3− a gH

p3 = p4

p4 = p5 + m gH

p5 = p6 + a gh2

p6 = p7 + a gh1 = p7 + a gLsen

p7 = pB

PA − PB = gH ( m − a ) + a gLsen (1)

A continuación, se hace un balance de cantidad de movimiento a partir de la ecuación (1.43)

1 d P –P 
(r rz ) =  A B +  gsen  (2)
r dr  L 
Luego reemplazando en (2) en (1)

1 d  gH ( m − a ) + a gLsen 
(r rz ) =  −  gsen  (3)
r dr  L 
Reemplazando la ecuación (1) en (3)

1 d  gH (  m −  a ) 
(r rz ) =   (4)
r dr  L 
Integrando la ecuación (4)

gH C
 rz = ( m – a )r + 1 (5)
2L r

Una condición de frontera aplicable a este problema, r = 0  rz = 0 ; por tanto para que el
problema tenga una solución acortada en el campo real, se requiere que C1 sea igual a cero. Luego
usando la ley de viscosidad de Newton

dVz gH
– = ( m – a ) r (6)
dr 2L
Integrando la ecuación (7)

gH ( m – a ) 2
Vz = – r + C2 (7)
4 L
La condición de frontera r = R Vz = 0 , permite obtener la constante C2 , es decir
gH ( m – a ) 2
C2 = R (8)
4 L
Reemplazando (9) en (8) se obtiene:

gH ( m – a ) 2   r  
2

Vz = R 1–    (9)
4   R  
El caudal se obtiene mediante:
2 R
Q=  rVz (r )drd (10)
0 0

Reemplazando (9) y (10) e integrando se obtiene:

 R 4 Hg
Q= ( m – a ) (11)
8 L
Luego, el flujo de masa es:

 R 4 Hg a

m = a  Q = ( m – a ) (12)
8 L