1.4-B Integración de las ecuaciones diferenciales de los bces.

microscópicos de continuidad y
cantidad de movimiento.
1.4.B-3 Hallar el perfil de velocidades, la velocidad máxima, el perfil de esfuerzos viscosos y la fuerza ejercida por el fluido
sobre el sólido, para el flujo axial a través de dos tubos concéntricos verticales.

Suposiciones:
Cumplimiento de la hipótesis del continuo.
Estado estacionario.
Fluido newtoniano.
Densidad constante.
Flujo laminar. R
L
kR

r
2.4, Fenómenos de Transporte. R. B. Bird, W. E. Stewart, E. N. Lightfoot

1
 Solución

a) Ecuación de Continuidad en coordenadas cilíndricas (r, , z):

𝜕𝜌 1 𝜕 1 𝜕 𝜕
+ 𝜌𝑟𝑣𝑟 + 𝜌𝑣𝜃 + 𝜌𝑣𝑧 = 0 𝑣𝜃 = 𝑣𝑟 = 0
𝜕𝑡 𝑟 𝜕𝑟 𝑟 𝜕𝜃 𝜕𝑧

𝛻. 𝜌𝑣

𝜕 𝜕𝑣𝑧
𝜌𝑣𝑧 = 0 𝜌 = 𝑐𝑡𝑒 𝜌 =0 𝑣𝑧 = 𝑐𝑡𝑒
𝜕𝑧 𝜕𝑧

El flujo “es desarrollado”; los elementos de volumen fluyen sin cambiar su velocidad en la dirección del
movimiento, ni se aceleran, ni se frenan.

2
a) Ecuación de Balance de Cantidad de Movimiento en coordenadas cilíndricas (r, , z)
Componente z:

𝜕𝑣𝑧 𝜕𝑣𝑧 𝑣𝜃 𝜕𝑣𝑧 𝜕𝑣𝑧 𝜕𝑝 1 𝜕 1 𝜕𝜃𝑧 𝜕𝑧𝑧

1 𝜕
𝑑𝑝 = − 𝑟𝑟𝑧 𝑑𝑧 − 𝜌𝑔𝑑𝑧 𝑔𝑧 = −𝑔
𝑟 𝜕𝑟

1𝜕 1 𝜕
න 𝑑𝑝 = − 𝑟𝑟𝑧 න 𝑑𝑧 − 𝜌𝑔 න 𝑑𝑧 𝑝=− 𝑟𝑟𝑧 𝑧 − 𝜌𝑔𝑧 + 𝐶1
𝑟 𝜕𝑟 𝑟 𝜕𝑟

3
 1 𝜕
𝑝=− 𝑟𝑟𝑧 𝑧 − 𝜌𝑔𝑧 + 𝐶1
𝑟 𝜕𝑟

𝑧=0 1𝜕
𝑝0 = − 𝑟𝑟𝑧 0 − 𝜌𝑔0 + 𝐶1
𝑝 = 𝑝0 𝑟 𝜕𝑟

𝐶1 = 𝑝0

𝑧=𝐿 1𝜕
−𝑝𝐿 + 𝑝0 = 𝑟𝑟𝑧 𝐿 + 𝜌𝑔𝐿
𝑝 = 𝑝𝐿 𝑟 𝜕𝑟

1𝜕
𝑝0 − 𝑝𝐿 = 𝑟𝑟𝑧 𝐿 + 𝜌𝑔𝐿
𝑟 𝜕𝑟

4
 El término P se denomina presión modificada y representa el efecto
P = 𝑝 + 𝜌𝑔𝑧 combinado de la presión y la fuerza gravitatoria, siendo z la distancia
medida hacia arriba desde un plano cualquiera, que se toma como
referencia.

P0 = 𝑝0
P0 − PL = 𝑝0 − 𝑝𝐿 − 𝜌𝑔𝐿
P𝐿 = 𝑝𝐿 + 𝜌𝑔𝐿

1𝜕
Reemplazando en −𝑝𝐿 + 𝑝0 = 𝑟𝑟𝑧 𝐿 + 𝜌𝑔𝐿
𝑟 𝜕𝑟

P0 − P𝐿 1 𝜕 P0 − P𝐿
= 𝑟𝑟𝑧 න 𝑟𝑑𝑟 = න 𝑑 𝑟𝑟𝑧
𝐿 𝑟 𝜕𝑟 𝐿

P0 − P𝐿 𝑟 2 𝐶2 𝑟 P0 − P𝐿 𝐶2
+ = 𝑟𝑧 𝑟𝑧 = 𝑟+
2𝐿 𝑟 𝑟 𝑟 2𝐿 𝑟

5
 La cte. de integración C2 no puede ser evaluada en forma inmediata, debido a que no se dispone de la información necesaria
acerca de la densidad de flujo de cantidad de movimiento en ninguna de las dos superficies, en r = KR ni en r = R.
Sin embargo, existirá un máximo de la curva de velocidad en un cierto plano -que hasta ahora no se conoce- siendo este plano
r =  R, para el cual la densidad de flujo de cantidad de movimiento es cero.

r = R P0 − PL 𝐶2
0= 𝑅 +
𝜏𝑟𝑧 = 0 2𝐿 𝑅

P0−PL P0 − PL 𝐶2
𝐶2 = − 2 𝑅 2 Reemplazando en 𝑟𝑧 = 𝑟+ R
2𝐿 2𝐿 𝑟
L
R

P0−PL P0−PL 2𝑅2 P0−PL 2 𝑅 2

P0 − PL 𝑟 𝑅 r
𝜏𝑟𝑧 = 𝑅 − 2 (1)
2𝐿 𝑅 𝑟 6
 𝑑𝑣𝑧
𝜏𝑟𝑧= -𝜇
𝑑𝑟

𝑑𝑣𝑧 P0 − PL 𝑟 𝑅
=− 𝑅 − 2
𝑑𝑟 2µ𝐿 𝑅 𝑟

𝑟 𝑑𝑢 1 1
Integrando: 𝑢= = 𝑑𝑢 = 𝑑𝑟 𝑑𝑟 = 𝑅 𝑑𝑢
𝑅 𝑑𝑟 𝑅 𝑅

P0 − PL 1
න 𝑑𝑣𝑧 = − 𝑅 𝑅 න 𝑢𝑑𝑢 −  𝑅 න 𝑑𝑢
2
2µ𝐿 𝑢

P0 − PL 2 𝑢2
𝑣𝑧 = − 𝑅 − 2 𝑙𝑛𝑢 + 𝐶3 (2)
2µ𝐿 2
7
 P0 − PL 2 1 𝑟 2 𝑟
𝑣𝑧 = − 𝑅 −  𝑙𝑛
2
+ 𝐶3 (2)
2µ𝐿 2 𝑅 𝑅

𝐶𝐶1 𝑟=𝑘𝑅 𝑣𝑧 = 0
Condiciones límites para 𝐶1
𝐶𝐶2 𝑟=𝑅 𝑣𝑧 = 0
Sustituyendo estas dos condiciones límites

P0 − PL 2 2 𝑘 2 𝑅2 𝑘𝑅
Para 𝐶𝐶1: 0=− 𝑅 − 2  𝑙𝑛(
2 ) + 2𝐶3
4µ𝐿 2 𝑅2 𝑅

Para 𝐶𝐶2:

P0 − PL 2 𝑅2 𝑅
0=− 𝑅 − 2  𝑙𝑛
2 + 2𝐶3 𝐶3 = − 1ൗ2
4µ𝐿 𝑅2 𝑅

8
 P0 − PL 2 2
Reemplazando 𝐶1 = − 1Τ2 en la ecuación resultante de aplicar las CC1 0=− 𝑅 𝑘 − 22 𝑙𝑛𝑘 + 2𝐶3
4µ𝐿

1 − 𝑘2 1
0= 𝑘2 − 22 𝑙𝑛𝑘 −1 2 = 𝑙𝑛 = 𝑙𝑛1 − 𝑙𝑛𝑘
1 𝑘
2 𝑙𝑛
𝑘

Sustituyendo en (1) y en (2)

P0 − PL 𝑟 1 − 𝑘2 𝑅 P0 − PL 2 𝑟 2 1 − 𝑘2 𝑟
𝜏𝑟𝑧 = 𝑅 − 𝑣𝑧 = 𝑅 [1 − +( ) 𝑙𝑛( )]
2𝐿 𝑅 1 𝑟 4µ𝐿 𝑅 1 𝑅
2 𝑙𝑛 𝑙𝑛
𝑘 𝑘

Cuando k = 0 las ecuaciones se transforman en las de flujo vertical en un tubo cilíndrico.

9
 R
L

R

z
kR

Observar que en la distribución de 𝑣𝑧 el perfil de velocidad cerca de la superficie del cilindro interior
alcanza más rápidamente el máximo dado que su superficie es menor.
10
La velocidad máxima será :

P0 − PL 2 1 − 𝑘2 1 − 𝑘2
𝑣𝑧 𝑚á𝑥 = 𝑣𝑧 𝑐𝑢𝑎𝑛𝑑𝑜 𝑟=R = 𝑅 1− 1 − 𝑙𝑛
4µ𝐿 1 1
2𝑙𝑛 2𝑙𝑛
𝑘 𝑘

En este caso la fuerza ejercida por el fluido sobre el sólido se obtiene sumando las fuerzas que actúan sobre los cilindros interior y
exterior.
𝐹𝑧 = − 𝑟𝑧ቚ 2𝜋𝑘𝑅𝐿 + 𝑟𝑧ቚ 2𝜋𝑅𝐿
𝑟=𝑘𝑅 𝑟=𝑅

P0 − PL 𝑘𝑅 𝑅 P0 − PL 𝑅 𝑅
𝐹𝑧 = − R [ ( ) −  ( ) ] 2π 𝑘 𝑅 𝐿 +
2 R [ ( ) −  ( ) ]2π 𝑅 𝐿
2
2𝐿 𝑅 𝑘𝑅 2𝐿 𝑅 𝑅

𝐹𝑧 = (P0 − PL)π 𝑅2 [− ( 𝑘 2 − 2 ) + (1 − 2 ) ] 𝐹𝑧 = P0 − PL 𝜋𝑅 2 1 − 𝑘2

11