Asignación Fenómenos de Transporte

Ejemplo 2.4>

-Presunciones

Estado estacionario

Flujo laminar

Fluido Newtoniano e incompresible

Coordenadas cilíndricas (r, 𝜃, 𝑧)

-Velocidades

𝑣𝑟 = 0
𝑣𝜃 = 0
𝑣𝑧 = 𝑣𝑧 (𝑟, 𝜃, 𝑧, 𝑡)
-Volumen de Control

Vc=2𝜋𝑟∆𝑟𝐿

-Balance de momentum

2𝜋𝑟𝐿(𝜏𝑟𝑧 |𝑟) − 2𝜋𝑟𝐿(𝜏𝑟𝑧 |𝑟 + ∆𝑟) + 2𝜋𝑟∆𝑟((𝜌𝑣𝑧 𝑣𝑧 |𝑧 = 0) − (𝜌𝑣𝑧 𝑣𝑧 |𝑧 = 𝐿)) + 2𝜋𝑟∆𝑟(𝑝0 − 𝑝𝐿 )

− 2𝜋𝑟∆𝑟𝐿 𝜌𝑔 = 0
Se cancela el termino de convección

2𝜋𝑟𝐿(𝜏𝑟𝑧 |𝑟) − 2𝜋𝑟𝐿(𝜏𝑟𝑧 |𝑟 + ∆𝑟) + 2𝜋𝑟∆𝑟((𝜌𝑣𝑧 𝑣𝑧 |𝑧 = 0) − (𝜌𝑣𝑧 𝑣𝑧 |𝑧 = 𝐿)) + 2𝜋𝑟∆𝑟(𝑝0 − 𝑝𝐿 )

− 2𝜋𝑟∆𝑟𝐿 𝜌𝑔 = 0
Se divide la ecuación entre las constantes y el diferencial del Volumen de control
2𝜋𝑟𝐿 2𝜋𝑟𝐿 2𝜋𝑟∆𝑟 2𝜋𝑟∆𝑟𝐿
(𝜏𝑟𝑧 |𝑟) − (𝜏𝑟𝑧 |𝑟 + ∆𝑟) + (𝑝0 − 𝑝𝐿 ) − 𝜌𝑔 = 0
2𝜋∆𝑟𝐿 2𝜋∆𝑟𝐿 2𝜋∆𝑟𝐿 2𝜋∆𝑟𝐿
(𝑟𝜏𝑟𝑧 |𝑟) − (𝑟𝜏𝑟𝑧 |𝑟 + ∆𝑟) 𝑟(𝑝0 − 𝑝𝐿 )
+ − 𝑟𝜌𝑔 = 0
∆𝑟 𝐿
(𝑟𝜏𝑟𝑧 |𝑟) − (𝑟𝜏𝑟𝑧 |𝑟 + ∆𝑟) (𝑝0 − 𝑝𝐿 )
+𝑟[ − 𝜌𝑔] = 0
∆𝑟 𝐿

(𝑟𝜏𝑟𝑧 |𝑟) − (𝑟𝜏𝑟𝑧 |𝑟 + ∆𝑟) (𝑝0 − 𝑝𝐿 ) − 𝜌𝑔𝐿

+𝑟[ ]=0
∆𝑟 𝐿

Se inserta el termino de presión modificada

(𝑟𝜏𝑟𝑧 |𝑟) − (𝜏𝑟𝑟𝑧 |𝑟 + ∆𝑟) 𝑃0 − 𝑃𝐿
+𝑟[ ]=0
∆𝑟 𝐿
Se aplica el lim para obtener definición de derivada
(𝑟𝜏𝑟𝑧 |𝑟) − (𝑟𝜏𝑟𝑧 |𝑟 + ∆𝑟) 𝑃0 − 𝑃𝐿
lim = −𝑟 [ ]
∆𝑟→0 ∆𝑟 𝐿
𝑑(𝑟𝜏𝑟𝑧 ) 𝑃0 − 𝑃𝐿
− = −𝑟 [ ]
𝑑𝑟 𝐿
𝑑(𝑟𝜏𝑟𝑧 ) 𝑃0 − 𝑃𝐿
= 𝑟[ ]
𝑑𝑟 𝐿
Se aplica el método de variables separables para solucionar la EDO
𝑃0 − 𝑃𝐿
∫ 𝑑(𝑟𝜏𝑟𝑧 ) = [ ] ∫ 𝑟𝑑𝑟
𝐿
𝑃0 − 𝑃𝐿 𝑟 2
𝑟𝜏𝑟𝑧 = + 𝐶1
𝐿 2
𝑃0 − 𝑃𝐿 𝑟 2 𝐶1
𝜏𝑟𝑧 = +
2𝐿 𝑟 𝑟
𝑃0 − 𝑃𝐿 𝐶1
𝜏𝑟𝑧 = 𝑟+ (1)
2𝐿 𝑟
Se reemplaza la primera condición limite

𝑟 = 𝜆𝑅 , 𝜏𝑟𝑧 = 0
𝑃0 − 𝑃𝐿 𝐶1
0= 𝜆𝑅 +
2𝐿 𝜆𝑅
 𝑃0 − 𝑃𝐿 𝐶1
− 𝜆𝑅 =
2𝐿 𝜆𝑅
𝑃0 − 𝑃𝐿
− (𝜆𝑅)2 = 𝐶1
2𝐿
Se reemplaza en (1)

𝑃0 − 𝑃𝐿 𝑃0 − 𝑃𝐿 (𝜆𝑅)2
𝜏𝑟𝑧 = 𝑟−
2𝐿 2𝐿 𝑟
𝑃0 − 𝑃𝐿 𝑃0 − 𝑃𝐿 (𝜆𝑅)2
𝜏𝑟𝑧 = 𝑟−
2𝐿 2𝐿 𝑟
𝑃0 − 𝑃𝐿 (𝜆𝑅)2
𝜏𝑟𝑧 = (𝑟 − )
2𝐿 𝑟

𝑃0 − 𝑃𝐿 𝜆2 𝑅2
𝜏𝑟𝑧 = (𝑟 − )
2𝐿 𝑟
(𝑃0 − 𝑃𝐿 )𝑅 𝑟 𝑅
𝜏𝑟𝑧 = ( − 𝜆2 )
2𝐿 𝑅 𝑟
Se reemplaza la ley de viscosidad de Newton
𝑑𝑣𝑧 (𝑃0 − 𝑃𝐿 )𝑅 𝑟 𝑅
−𝜇 = ( − 𝜆2 )
𝑑𝑟 2𝐿 𝑅 𝑟
Se aplica variables separables nuevamente
(𝑃0 − 𝑃𝐿 )𝑅 𝑟 𝑅
𝑑𝑣𝑧 = − ( − 𝜆2 ) 𝑑𝑟
2𝜇𝐿 𝑅 𝑟
(𝑃0 − 𝑃𝐿 )𝑅 𝑟 𝑅
∫ 𝑑𝑣𝑧 = − ∫ ( − 𝜆2 ) 𝑑𝑟
2𝜇𝐿 𝑅 𝑟
(𝑃0 − 𝑃𝐿 )𝑅 𝑟 𝑅
𝑣𝑧 = − [∫ 𝑑𝑟 − ∫ 𝜆2 𝑑𝑟]
2𝜇𝐿 𝑅 𝑟
(𝑃0 − 𝑃𝐿 )𝑅 𝑟 𝑅
𝑣𝑧 = − [∫ 𝑑𝑟 − ∫ 𝜆2 𝑑𝑟]
2𝜇𝐿 𝑅 𝑟
(𝑃0 − 𝑃𝐿 )𝑅 1 1
𝑣𝑧 = − [ ∫ 𝑟𝑑𝑟 − 𝜆2 𝑅 ∫ 𝑑𝑟]
2𝜇𝐿 𝑅 𝑟
(𝑃0 − 𝑃𝐿 )𝑅 𝑟 2
𝑣𝑧 = − [ − 𝜆2 𝑅 𝑙𝑛𝑟 + 𝐶2 ] (2)
2𝜇𝐿 2𝑅

Se plantea las condiciones limite para plantear el cálculo de 𝜆

𝑟 = 𝑘𝑅 , 𝑣𝑧 = 0
𝑟=𝑅 , 𝑣𝑧 = 0
Se reemplazan las condiciones limite en (2)

(𝑃0 − 𝑃𝐿 )𝑅 (𝑘𝑅 )2
0=− [ − 𝜆2 𝑅 ln(𝑘𝑅) + 𝐶2 ] (3)
2𝜇𝐿 2𝑅

(𝑃0 − 𝑃𝐿 )𝑅 𝑅 2
0=− [ − 𝜆2 𝑅 ln(𝑅) + 𝐶2 ] (4)
2𝜇𝐿 2𝑅
(𝑃0 −𝑃𝐿 )𝑅
Se pasa a dividir el termino − 2𝜇𝐿
en las 2 ecuaciones

(𝑘𝑅 )2
0= − 𝜆2 𝑅 ln(𝑘𝑅) + 𝐶2 (5)
2𝑅
𝑅2
0= − 𝜆2 𝑅 ln(𝑅) + 𝐶2 (6)
2𝑅
Se resta (5) – (6)

(𝑘𝑅 )2 𝑅 2
0= − − 𝜆2 𝑅 ln(𝑘𝑅) + 𝜆2 𝑅 ln(𝑅)
2𝑅 2𝑅
A partir de aquí se desarrolla la ecuación para conocer la función de 𝜆

𝑘 2𝑅 𝑅
0= − − 𝜆2 𝑅 ln(𝑘𝑅) + 𝜆2 𝑅 ln(𝑅)
2 2
𝑘2 1
0 = 𝑅( − − 𝜆2 ln(𝑘𝑅) + 𝜆2 ln(𝑅))
2 2

𝑘2−1
0= − 𝜆2 ln(𝑘𝑅) + 𝜆2 ln(𝑅)
2
𝑘2−1
0= + 𝜆2 (ln(𝑅) − ln(𝑘𝑅))
2
𝑘2−1 𝑅
0= + 𝜆2 (ln ( ))
2 𝑘𝑅
𝑘2−1 1
0= + 𝜆2 ln ( )
2 𝑘
1
0 = 𝑘 2 − 1 + 2𝜆2 ln ( )
𝑘
1
1 − 𝑘 2 = +2𝜆2 ln ( )
𝑘
1−𝑘2
= 2𝜆2
1
ln ( )
𝑘
1−𝑘2
2𝜆2 =
1
ln ( )
𝑘
 1−𝑘2
𝜆2 =
1
2 ln ( )
𝑘

1−𝑘2
𝜆=√
1
2 ln ( )
𝑘

Este 𝜆 es una constante la cual tiene un valor entre 0 y 1, siendo esta mayor que k. Dicha
constante al ser multiplicada por el Radio del tubo exterior representa el punto donde la
densidad de flujo de momentum toma un valor de 0.

Ahora, si se retira el cilindro interior k=0, entonces se plantea

1−𝑘2
𝜆=√
1
2 ln ( )
𝑘

1−𝑘2
lim 𝜆 = √
𝑘→0 1
2 ln ( )
𝑘

1 − 02
lim 𝜆 = √
𝑘→0 1
2 ln ( )
0

1
lim 𝜆 = √
𝑘→0 2∞

1
lim 𝜆 = √
𝑘→0 ∞

lim 𝜆 = √0
𝑘→0

lim 𝜆 = 0
𝑘→0

Por tanto 𝜆 = 0 cuando el cilindro interior no existe

𝑘+1
Ahora se plantea si 𝜆 esta mas cerca del cilindro interior o el exterior, donde es la relación
2
para la distancia media entre los cilindros interior y exterior

𝑘+1 1−𝑘2
>√
2 1
2 ln ( )
𝑘
Como 0<k<1 se le dan valores a k para ver si se cumple esta inecuación, sino se cumple
𝜆 estará mas cerca del cilindro exterior

K=0,5

0,5 + 1 1 − 0,52
>√
2 1
2 ln ( )
0,5

0,75 > 0,7355

Esta inecuación se cumple, por tanto 𝜆 esta mas cerca del cilindro interior.

Ejemplo 2.5 modificado

Al plantear el ejercicio se ve que la densidad no afecta las distribuciones de velocidad ni de

momentum, por tanto, todo el ejercicio será igual su planteamiento solo que al cambiar las
viscosidades el perfil de esfuerzos cambiara y el de velocidades también

-Presunciones

Estado estacionario

Flujo laminar

Fluido Newtoniano e incompresible

Coordenadas rectangulares ( 𝑥, 𝑦, 𝑧)

-Velocidades

𝑣𝑥 = 0
𝑣𝑦 = 0
𝑣𝑧 = 𝑣𝑧 ( 𝑥, 𝑦, 𝑧, 𝑡)
-Volumen de Control

Vc=𝐿𝑊∆𝑥

-Balance de momentum

𝐿𝑊(𝜏𝑥𝑧 |𝑟) − 𝐿𝑊(𝜏𝑥𝑧 |𝑟 + ∆𝑟) + 𝑊∆𝑥((𝜌𝑣𝑧 𝑣𝑧 |𝑧 = 0) − (𝜌𝑣𝑧 𝑣𝑧 |𝑧 = 𝐿)) + 𝑊∆𝑥(𝑝0 − 𝑝𝐿 ) = 0

Se cancela el termino de convección

𝐿𝑊(𝜏𝑥𝑧 |𝑥) − 𝐿𝑊(𝜏𝑥𝑧 |𝑥 + ∆𝑥) + 𝑊∆𝑥((𝜌𝑣𝑧 𝑣𝑧 |𝑧 = 0) − (𝜌𝑣𝑧 𝑣𝑧 |𝑧 = 𝐿)) + 𝑊∆𝑥(𝑝0 − 𝑝𝐿 ) = 0

Se divide la ecuación entre las constantes y el diferencial del Volumen de control

𝐿𝑊 𝐿𝑊 𝑊∆𝑥 𝐿𝑊∆𝑥
(𝜏𝑥𝑧 |𝑥) − (𝜏𝑥𝑧 |𝑥 + ∆𝑥) + (𝑝0 − 𝑝𝐿 ) − 𝜌𝑔 = 0
𝐿𝑊∆𝑥 𝐿𝑊∆𝑥 𝐿𝑊∆𝑥 𝐿𝑊∆𝑥
(𝜏𝑥𝑧 |𝑥) − (𝜏𝑥𝑧 |𝑥 + ∆𝑥) (𝑝0 − 𝑝𝐿 )
+ =0
∆𝑥 𝐿
Se aplica el lim para obtener definición de derivada
(𝜏𝑥𝑧 |𝑥) − (𝜏𝑥𝑧 |𝑥 + ∆𝑥) (𝑝0 − 𝑝𝐿 )
lim =−
∆𝑥→0 ∆𝑥 𝐿
𝑑(𝜏𝑥𝑧 ) (𝑝0 − 𝑝𝐿 )
− =−
𝑑𝑥 𝐿
𝑑(𝜏𝑥𝑧 ) (𝑝0 − 𝑝𝐿 )
=
𝑑𝑥 𝐿
Se aplica el método de variables separables para solucionar la EDO
(𝑝0 − 𝑝𝐿 )
∫ 𝑑(𝜏𝑥𝑧 ) = ∫ 𝑑𝑥
𝐿
(𝑝0 − 𝑝𝐿 )
𝜏𝑥𝑧 = 𝑥 + 𝐶1 (1)
𝐿
Como se tienen dos fases se obtienen dos expresiones para 𝜏𝑥𝑧
(𝑝0 − 𝑝𝐿 )
𝜏𝑥𝑧 𝐼 = 𝑥 + 𝐶1 𝐼 (2)
𝐿
(𝑝0 − 𝑝𝐿 )
𝜏𝑥𝑧 𝐼𝐼 = 𝑥 + 𝐶1 𝐼𝐼 (3)
𝐿
Se plantea la condición limite entre la interfase liq-liq

𝑥=0 , 𝜏𝑥𝑧 𝐼𝐼 = 𝜏𝑥𝑧 𝐼

Por tanto
 (𝑝0 − 𝑝𝐿 ) (𝑝0 − 𝑝𝐿 )
0 + 𝐶1 𝐼𝐼 = 0 + 𝐶1 𝐼
𝐿 𝐿
𝐶1 𝐼𝐼 = 𝐶1 𝐼
Se puede plantear

𝐶1 = 𝐶1 𝐼𝐼 = 𝐶1 𝐼
Se continua con el planteamiento, se reemplaza la ley de viscosidad de Newton
𝑑𝑣𝑧 (𝑝0 − 𝑝𝐿 )
−𝜇 = 𝑥 + 𝐶1
𝑑𝑥 𝐿
𝑑𝑣𝑧 (𝑝0 − 𝑝𝐿 ) 𝐶1
=− 𝑥−
𝑑𝑥 𝜇𝐿 𝜇
Se aplica variables separables nuevamente
(𝑝0 − 𝑝𝐿 ) 𝐶1
∫ 𝑑𝑣𝑧 = − ∫ 𝑥 𝑑𝑥 − ∫ 𝑑𝑥
𝜇𝐿 𝜇
(𝑝0 − 𝑝𝐿 ) 2 𝐶1
𝑣𝑧 = − 𝑥 − 𝑥 + 𝐶2 (4)
2𝜇𝐿 𝜇
Para cada fase se obtiene
(𝑝0 − 𝑝𝐿 ) 2 𝐶1
𝑣𝑧 𝐼 = − 𝐼
𝑥 − 𝐼 𝑥 + 𝐶2 𝐼 (5)
2𝜇 𝐿 𝜇
(𝑝0 − 𝑝𝐿 ) 2 𝐶1
𝑣𝑧 𝐼𝐼 = − 𝑥 − 𝐼𝐼 𝑥 + 𝐶2 𝐼𝐼 (6)
2𝜇𝐼𝐼 𝐿 𝜇

Se plantea las condiciones límites para hallar c1, 𝐶2 𝐼𝐼 , 𝐶2 𝐼

𝑥=0 , 𝑣𝑧 𝐼 = 𝑣𝑧 𝐼𝐼

𝑥 = −𝑏 , 𝑣𝑧 𝐼 = 0

𝑥=𝑏 , 𝑣𝑧 𝐼𝐼 = 0
Se reemplazan las condiciones límite:
(𝑝0 − 𝑝𝐿 ) 2 𝐶1 (𝑝0 − 𝑝𝐿 ) 2 𝐶1
− 𝐼𝐼
0 − 𝐼𝐼 0 + 𝐶2 𝐼𝐼 = − 0 − 𝐼 0 + 𝐶2 𝐼
2𝜇 𝐿 𝜇 2𝜇𝐼 𝐿 𝜇

𝐶2 𝐼𝐼 = 𝐶2 𝐼
Se plantea igual que C1

𝐶2 = 𝐶2 𝐼𝐼 = 𝐶2 𝐼
 (𝑝0 − 𝑝𝐿 ) 2 𝐶1
0=− 𝑏 + 𝐼 𝑏 + 𝐶2 (7)
2𝜇𝐼 𝐿 𝜇
(𝑝0 − 𝑝𝐿 ) 2 𝐶1
0=− 𝑏 − 𝐼𝐼 𝑏 + 𝐶2 (8)
2𝜇𝐼𝐼 𝐿 𝜇
Se restan (7) – (8) y se realizan despejes para calcular C1
(𝑝0 − 𝑝𝐿 ) 2 (𝑝0 − 𝑝𝐿 ) 2 𝐶1 𝐶1
0=− 𝐼
𝑏 + 𝐼𝐼
𝑏 + 𝐼 𝑏 + 𝐼𝐼 𝑏 + 𝐶2 − 𝐶2
2𝜇 𝐿 2𝜇 𝐿 𝜇 𝜇
(𝑝0 − 𝑝𝐿 ) 2 1 1 1 1
0= 𝑏 ( 𝐼𝐼 − 𝐼 ) + 𝐶1 𝑏 ( 𝐼 + 𝐼𝐼 )
2𝐿 𝜇 𝜇 𝜇 𝜇
(𝑝0 − 𝑝𝐿 ) 2 𝜇𝐼 − 𝜇𝐼𝐼 𝜇𝐼𝐼 + 𝜇𝐼
0= 𝑏 ( 𝐼𝐼 𝐼 ) + 𝐶1 𝑏 ( 𝐼 𝐼𝐼 )
2𝐿 𝜇 𝜇 𝜇𝜇

(𝑝0 − 𝑝𝐿 ) 2 𝜇𝐼 − 𝜇𝐼𝐼 𝜇𝐼𝐼 + 𝜇𝐼

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

Se reemplaza C1 en (7)

(𝑝0 − 𝑝𝐿 ) 2 (𝑝0 − 𝑝𝐿 ) 2 𝜇𝐼 − 𝜇𝐼𝐼

0=− 𝑏 − 𝑏 ( 𝐼𝐼 ) + 𝐶2
2𝜇𝐼 𝐿 2𝜇𝐼 𝐿 𝜇 + 𝜇𝐼

Se calcula C2

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

(𝑝0 − 𝑝𝐿 ) 2 𝜇𝐼𝐼 + 𝜇𝐼 𝜇𝐼 − 𝜇𝐼𝐼

0=− 𝑏 ( 𝐼𝐼 + ) + 𝐶2
2𝜇𝐼 𝐿 𝜇 + 𝜇𝐼 𝜇𝐼𝐼 + 𝜇𝐼

(𝑝0 − 𝑝𝐿 ) 2 𝜇𝐼𝐼 + 𝜇𝐼 + 𝜇𝐼 − 𝜇𝐼𝐼

0=− 𝑏 ( ) + 𝐶2
2𝜇𝐼 𝐿 𝜇𝐼𝐼 + 𝜇𝐼

(𝑝0 − 𝑝𝐿 ) 2 2𝜇𝐼
0=− 𝑏 ( ) + 𝐶2
2𝜇𝐼 𝐿 𝜇𝐼𝐼 + 𝜇𝐼

(𝑝0 − 𝑝𝐿 ) 2 2𝜇𝐼
𝑏 ( ) = 𝐶2
2𝜇𝐼 𝐿 𝜇𝐼𝐼 + 𝜇𝐼

Se reemplazan en las ecuaciones (5) y (6)

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

(𝑝0 − 𝑝𝐿 ) 2 (𝑝0 − 𝑝𝐿 ) 𝜇𝐼 − 𝜇𝐼𝐼 (𝑝0 − 𝑝𝐿 ) 2 𝜇𝐼𝐼

𝑣𝑧 𝐼𝐼 = − 𝑥 + 𝑏 ( ) 𝑥 + 𝑏 ( )
2𝜇𝐼𝐼 𝐿 2𝜇𝐼𝐼 𝐿 𝜇𝐼𝐼 + 𝜇𝐼 2𝜇𝐼𝐼 𝐿 𝜇𝐼𝐼 + 𝜇𝐼

Se reordenan las ecuaciones

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

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

Como C1 es igual para las dos fases el perfil de esfuerzos es el siguiente

(𝑝0 − 𝑝𝐿 ) (𝑝0 − 𝑝𝐿 ) 𝜇𝐼 − 𝜇𝐼𝐼

𝜏𝑥𝑧 = 𝑥− 𝑏 ( 𝐼𝐼 )
𝐿 2𝐿 𝜇 + 𝜇𝐼

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

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

Pregunta para discusión 11

No es posible, debido a que si se disponen los fluidos como se ve en la figura ambos tienen picos de
velocidad esto provocara la existencia de un punto intermedio en el cual se reduce la velocidad y este
planteamiento impediría la relación lineal demostrada del esfuerzo cortante en el anterior punto. Por
tanto, en fluidos inmiscibles solo existe un pico en el perfil de velocidad para ambas capas y no como
se presenta aquí.

Problema 2B-4
 V0

-Presunciones

Estado estacionario

Flujo laminar

Fluido Newtoniano e incompresible

Coordenadas rectangulares ( 𝑥, 𝑦, 𝑧)

-Velocidades

𝑣𝑥 = 0
𝑣𝑦 = 0

𝑣𝑧 = 𝑣𝑧 ( 𝑥, 𝑦, 𝑧, 𝑡)
-Volumen de Control

Vc=𝐿𝑊∆𝑥

-Balance de momentum

𝐿𝑊(𝜏𝑥𝑧 |𝑟) − 𝐿𝑊(𝜏𝑥𝑧 |𝑟 + ∆𝑟) + 𝑊∆𝑥((𝜌𝑣𝑧 𝑣𝑧 |𝑧 = 0) − (𝜌𝑣𝑧 𝑣𝑧 |𝑧 = 𝐿)) + 𝑊∆𝑥(𝑝0 − 𝑝𝐿 )

− 𝐿𝑊∆𝑥 𝜌𝑔 = 0
Se cancela el termino de convección

𝐿𝑊(𝜏𝑥𝑧 |𝑥) − 𝐿𝑊(𝜏𝑥𝑧 |𝑥 + ∆𝑥) + 𝑊∆𝑥((𝜌𝑣𝑧 𝑣𝑧 |𝑧 = 0) − (𝜌𝑣𝑧 𝑣𝑧 |𝑧 = 𝐿)) + 𝑊∆𝑥(𝑝0 − 𝑝𝐿 )

− 𝐿𝑊∆𝑥 𝜌𝑔 = 0
Se divide la ecuación entre las constantes y el diferencial del Volumen de control
𝐿𝑊 𝐿𝑊 𝑊∆𝑥 𝐿𝑊∆𝑥
(𝜏𝑥𝑧 |𝑥) − (𝜏𝑥𝑧 |𝑥 + ∆𝑥) + (𝑝0 − 𝑝𝐿 ) − 𝜌𝑔 = 0
𝐿𝑊∆𝑥 𝐿𝑊∆𝑥 𝐿𝑊∆𝑥 𝐿𝑊∆𝑥
(𝜏𝑥𝑧 |𝑥) − (𝜏𝑥𝑧 |𝑥 + ∆𝑥) (𝑝0 − 𝑝𝐿 )
+ − 𝜌𝑔 = 0
∆𝑥 𝐿
(𝜏𝑥𝑧 |𝑥) − (𝜏𝑥𝑧 |𝑥 + ∆𝑥) (𝑝0 − 𝑝𝐿 )
+[ − 𝜌𝑔] = 0
∆𝑥 𝐿

(𝜏𝑥𝑧 |𝑥) − (𝜏𝑥𝑧 |𝑥 + ∆𝑥) (𝑝0 − 𝑝𝐿 ) − 𝜌𝑔𝐿

+[ ]=0
∆𝑥 𝐿

Se inserta el termino de presión modificada

(𝜏𝑥𝑧 |𝑥) − (𝜏𝑥𝑧 |𝑥 + ∆𝑥) 𝑃0 − 𝑃𝐿
+[ ]=0
∆𝑥 𝐿
Se aplica el lim para obtener definición de derivada
(𝜏𝑥𝑧 |𝑥) − (𝜏𝑥𝑧 |𝑥 + ∆𝑥) 𝑃0 − 𝑃𝐿
lim = −[ ]
∆𝑥→0 ∆𝑥 𝐿
𝑑(𝜏𝑥𝑧 ) 𝑃0 − 𝑃𝐿
− = −[ ]
𝑑𝑥 𝐿
𝑑(𝜏𝑥𝑧 ) 𝑃0 − 𝑃𝐿
=[ ]
𝑑𝑥 𝐿
Se aplica el método de variables separables para solucionar la EDO
𝑃0 − 𝑃𝐿
∫ 𝑑(𝜏𝑥𝑧 ) = [ ] ∫ 𝑑𝑥
𝐿
𝑃0 − 𝑃𝐿
𝜏𝑥𝑧 = 𝑥 + 𝐶1 (1)
𝐿
Se reemplaza la ley de viscosidad de Newton
𝑑𝑣𝑧 (𝑃0 − 𝑃𝐿 )
−𝜇 = 𝑥 + 𝐶1
𝑑𝑥 𝐿
𝑑𝑣𝑧 (𝑃0 − 𝑃𝐿 ) 𝐶1
=− 𝑥−
𝑑𝑥 𝜇𝐿 𝜇
Se plantea que
𝐶1
− = 𝑘1
𝜇
Por tanto
𝑑𝑣𝑧 (𝑃0 − 𝑃𝐿 )
=− 𝑥 + 𝑘1
𝑑𝑥 𝜇𝐿
Se aplica variables separables nuevamente
(𝑃0 − 𝑃𝐿 )
∫ 𝑑𝑣𝑧 = − ∫ 𝑥 𝑑𝑥 + 𝑘1 ∫ 𝑑𝑥
𝜇𝐿
(𝑃0 − 𝑃𝐿 ) 2
𝑣𝑧 = − 𝑥 + 𝑘1 𝑥 + 𝐶2 (2)
2𝜇𝐿
Se plantea las condiciones límites para hallar k1 y c2

𝑥=𝐵 , 𝑣𝑧 = 𝑣0
𝑥 = −𝐵 , 𝑣𝑧 = 0
Se reemplazan en (2)
(𝑃0 − 𝑃𝐿 ) 2
𝑣0 = − 𝐵 + 𝑘1 𝐵 + 𝐶2 (3)
2𝜇𝐿
(𝑃0 − 𝑃𝐿 ) 2
0=− 𝐵 − 𝑘1 𝐵 + 𝐶2 (4)
2𝜇𝐿
Restamos (3)-(4)
(𝑃0 − 𝑃𝐿 ) 2 (𝑃0 − 𝑃𝐿 ) 2
𝑣0 = − 𝐵 + 𝐵 + 𝑘1 𝐵 + 𝑘1 𝐵 + 𝐶2 − 𝐶2
2𝜇𝐿 2𝜇𝐿
𝑣0 = 2𝑘1 𝐵
𝑣0
= 𝑘1
2𝐵
Se reemplaza k1 en (4)
(𝑃0 − 𝑃𝐿 ) 2 𝑣0
0=− 𝐵 − 𝐵 + 𝐶2 (4)
2𝜇𝐿 2𝐵
(𝑃0 − 𝑃𝐿 ) 2 𝑣0
0=− 𝐵 − + 𝐶2
2𝜇𝐿 2
(𝑃0 − 𝑃𝐿 ) 2 𝑣0
𝐵 + = 𝐶2
2𝜇𝐿 2
Se reemplaza k1 y C2 en (2)
(𝑃0 − 𝑃𝐿 ) 2 𝑣0 (𝑃0 − 𝑃𝐿 ) 2 𝑣0
𝑣𝑧 = − 𝑥 + 𝑥+ 𝐵 +
2𝜇𝐿 2𝐵 2𝜇𝐿 2
Agrupamos términos
(𝑃0 − 𝑃𝐿 ) 2 (𝑃0 − 𝑃𝐿 ) 2 𝑣0 𝑥 𝑣0
𝑣𝑧 = + 𝐵 − 𝑥 + +
2𝜇𝐿 2𝜇𝐿 2𝐵 2
(𝑃0 − 𝑃𝐿 ) 2 𝑣0 𝑥
𝑣𝑧 = [𝐵 − 𝑥 2 ] + ( + 1)
2𝜇𝐿 2 𝐵
 (𝑃0 − 𝑃𝐿 )𝐵2 𝑥2 𝑣0 𝑥
𝑣𝑧 = [1 − 2 ] + (1 + )
2𝜇𝐿 𝐵 2 𝐵

(𝑃0 − 𝑃𝐿 )𝐵2 𝑥 2 𝑣0 𝑥
𝑣𝑧 = [1 − ( ) ] + (1 + )
2𝜇𝐿 𝐵 2 𝐵
Ahora k1 se reemplaza en (1)
𝑃0 − 𝑃𝐿
𝜏𝑥𝑧 = 𝑥 + 𝐶1
𝐿
Teniendo en cuenta que
𝐶1
− = 𝑘1
𝜇
𝐶1 = −𝑘1 𝜇
𝑣0
𝐶1 = − 𝜇
2𝐵

Entonces
𝑃0 − 𝑃𝐿 𝜇 𝑣0
𝜏𝑥𝑧 = ( ) 𝑥−
𝐿 2𝐵
-Presunciones

Estado estacionario

Flujo laminar

Fluido Newtoniano e incompresible

Coordenadas cilíndricas (r, 𝜃, 𝑧)

-Velocidades

𝑣𝑟 = 0
𝑣𝜃 = 0
𝑣𝑧 = 𝑣𝑧 (𝑟, 𝜃, 𝑧, 𝑡)
-Volumen de Control

Vc=2𝜋𝑟∆𝑟𝐿

-Balance de momentum

2𝜋𝑟𝐿(𝜏𝑟𝑧 |𝑟) − 2𝜋𝑟𝐿(𝜏𝑟𝑧 |𝑟 + ∆𝑟) + 2𝜋𝑟∆𝑟((𝜌𝑣𝑧 𝑣𝑧 |𝑧 = 0) − (𝜌𝑣𝑧 𝑣𝑧 |𝑧 = 𝐿)) + 2𝜋𝑟∆𝑟𝐿 𝜌𝑔 = 0

Se cancela el termino de convección

2𝜋𝑟𝐿(𝜏𝑟𝑧 |𝑟) − 2𝜋𝑟𝐿(𝜏𝑟𝑧 |𝑟 + ∆𝑟) + 2𝜋𝑟∆𝑟((𝜌𝑣𝑧 𝑣𝑧 |𝑧 = 0) − (𝜌𝑣𝑧 𝑣𝑧 |𝑧 = 𝐿)) + 2𝜋𝑟∆𝑟𝐿 𝜌𝑔 = 0

Se divide la ecuación entre las constantes y el diferencial del Volumen de control

2𝜋𝑟𝐿 2𝜋𝑟𝐿 2𝜋𝑟∆𝑟𝐿
(𝜏𝑟𝑧 |𝑟) − (𝜏𝑟𝑧 |𝑟 + ∆𝑟) + 𝜌𝑔 = 0
2𝜋∆𝑟𝐿 2𝜋∆𝑟𝐿 2𝜋∆𝑟𝐿
(𝑟𝜏𝑟𝑧 |𝑟) − (𝑟𝜏𝑟𝑧 |𝑟 + ∆𝑟)
+ 𝑟𝜌𝑔 = 0
∆𝑟
Se aplica el lim para obtener definición de derivada
(𝑟𝜏𝑟𝑧 |𝑟) − (𝑟𝜏𝑟𝑧 |𝑟 + ∆𝑟)
lim = −𝑟𝜌𝑔
∆𝑟→0 ∆𝑟
𝑑(𝑟𝜏𝑟𝑧 )
− = −𝑟𝜌𝑔
𝑑𝑟
𝑑(𝑟𝜏𝑟𝑧 )
= 𝑟𝜌𝑔
𝑑𝑟
Se aplica el método de variables separables para solucionar la EDO

∫ 𝑑(𝑟𝜏𝑟𝑧 ) = 𝜌𝑔 ∫ 𝑟𝑑𝑟

𝑟2
𝑟𝜏𝑟𝑧 = 𝜌𝑔 + 𝐶1
2
 𝜌𝑔 𝑟 2 𝐶1
𝜏𝑟𝑧 = +
2 𝑟 𝑟
𝜌𝑔 𝐶1
𝜏𝑟𝑧 = 𝑟+ (1)
2 𝑟
Se plantea la primera condición limite

𝑟 = 𝑎𝑅 , 𝜏𝑟𝑧 = 0
Se reemplaza en (1)
𝜌𝑔 𝐶1
0= 𝑎𝑅 +
2 𝑎𝑅
𝜌𝑔 𝐶1
− =
2 𝑎𝑅
𝜌𝑔
− (𝑎𝑅)2 = 𝐶1
2
Se reemplaza en (1)

𝜌𝑔 𝜌𝑔 (𝑎𝑅)2
𝜏𝑟𝑧 = 𝑟−
2 2 𝑟
𝜌𝑔 (𝑎𝑅)2
𝜏𝑟𝑧 = (𝑟 − )
2 𝑟

𝜌𝑔 𝑎2 𝑅 2
𝜏𝑟𝑧 = (𝑟 − )
2 𝑟
𝜌𝑔𝑅 𝑟 𝑅
𝜏𝑟𝑧 = ( − 𝑎2 )
2 𝑅 𝑟
Se reemplaza la ley de viscosidad de Newton
𝑑𝑣𝑧 𝜌𝑔𝑅 𝑟 𝑅
−𝜇 = ( − 𝑎2 )
𝑑𝑟 2 𝑅 𝑟
Se aplica variables separables nuevamente
𝜌𝑔𝑅 𝑟 𝑅
𝑑𝑣𝑧 = − ( − 𝑎2 ) 𝑑𝑟
2𝜇 𝑅 𝑟
𝜌𝑔𝑅 𝑟 𝑅
∫ 𝑑𝑣𝑧 = − ∫ ( − 𝑎2 ) 𝑑𝑟
2𝜇 𝑅 𝑟
𝜌𝑔𝑅 𝑟 𝑅
𝑣𝑧 = − [∫ 𝑑𝑟 − ∫ 𝑎2 𝑑𝑟]
2𝜇 𝑅 𝑟
𝜌𝑔𝑅 1 1
𝑣𝑧 = − [ ∫ 𝑟𝑑𝑟 − 𝑎2 𝑅 ∫ 𝑑𝑟]
2𝜇 𝑅 𝑟
 𝜌𝑔𝑅 𝑟 2
𝑣𝑧 = − [ − 𝑎2 𝑅 ln(𝑟) + 𝐶2 ] (2)
2𝜇 2𝑅

Se plantea la segunda condición limite

𝑟=𝑅 , 𝑣𝑧 = 0
Se reemplaza en (2)

𝜌𝑔𝑅 𝑅 2
0=− [ − 𝑎2 𝑅 ln(𝑅) + 𝐶2 ]
2𝜇 2𝑅

Se despeja C2

𝑅2
0= − 𝑎2 𝑅 ln(𝑅) + 𝐶2
2𝑅
𝑅
0= − 𝑎2 𝑅 ln(𝑅) + 𝐶2
2
𝑅
𝑎2 𝑅 ln(𝑅) − = 𝐶2
2
Se reemplaza c2 en (2)

𝜌𝑔𝑅 𝑟 2 𝑅
𝑣𝑧 = − [ − 𝑎2 𝑅 ln(𝑟) + 𝑎2 𝑅 ln(𝑅) − ]
2𝜇 2𝑅 2

Se saca factor común R/2 nuevamente:

𝜌𝑔𝑅 2 𝑟 2
𝑣𝑧 = − [ − 2𝑎2 ln(𝑟) + 2𝑎2 ln(𝑅) − 1]
4𝜇 𝑅 2

El menos distribuye en el corchete:

𝜌𝑔𝑅 2 𝑟 2
𝑣𝑧 = − [ − 2𝑎2 ln(𝑟) + 2𝑎2 ln(𝑅) − 1]
4𝜇 2𝑅 2

𝜌𝑔𝑅 2 𝑟2
𝑣𝑧 = [1 − 2 + 2𝑎2 ln(𝑟) − 2𝑎2 ln(𝑅)]
4𝜇 𝑅

Se saca en el corchete un factor común para 2𝑎2

𝜌𝑔𝑅 2 𝑟2
𝑣𝑧 = [1 − 2 + 2𝑎2 (ln(𝑟) − ln(𝑅))] (3)
4𝜇 𝑅

Se aplica propiedad de resta de logaritmos

𝜌𝑔𝑅 2 𝑟2 𝑟
𝑣𝑧 = [1 − 2 + 2𝑎2 ln ( )]
4𝜇 𝑅 𝑅

𝜌𝑔𝑅 2 𝑟 2 𝑟
𝑣𝑧 = [1 − ( ) + 2𝑎2 ln ( )] 𝑎)
4𝜇 𝑅 𝑅
Para calcular el flujo masico se plantea
2𝜋 𝑎𝑅
𝜔 =𝜌∗∫ ∫ 𝑣𝑧 𝑟𝑑𝑟𝑑𝜃
0 𝑅

Se reemplaza vz pero la ecuación (3) para desarrollar las dobles integrales:

2𝜋 𝑎𝑅
𝜌𝑔𝑅 2 𝑟2
𝜔 =𝜌∗∫ ∫ [ [1 − 2 + 2𝑎2 (ln(𝑟) − ln(𝑅))]] 𝑟𝑑𝑟𝑑𝜃
0 𝑅 4𝜇 𝑅

𝜌𝑔𝑅2
Se saca 4𝜇
de las integrales ya que es un factor constante

𝜌𝑔𝑅 2 2𝜋 𝑎𝑅 𝑟2
𝜔=𝜌∗ ∫ ∫ [1 − 2 + 2𝑎2 (ln(𝑟) − ln(𝑅))] 𝑟𝑑𝑟𝑑𝜃
4𝜇 0 𝑅 𝑅

Se distribuye el rdr que cada uno de los factores del corchete

𝜌𝑔𝑅 2 2𝜋 𝑎𝑅 𝑟2
𝜔=𝜌∗ ∫ ∫ [𝑟𝑑𝑟 − 2 𝑟𝑑𝑟 + 2𝑎2 (ln(𝑟)𝑟𝑑𝑟 − ln(𝑅)𝑟𝑑𝑟)] 𝑑𝜃
4𝜇 0 𝑅 𝑅

Se distribuye la integral para cada factor

𝜌𝑔𝑅 2 2𝜋 𝑎𝑅 𝑎𝑅 2
𝑟 𝑎𝑅
𝜔=𝜌∗ ∫ [∫ 𝑟𝑑𝑟 ∫ 2
𝑟𝑑𝑟 + ∫ 2𝑎2 (ln(𝑟)𝑟𝑑𝑟 − ln(𝑅)𝑟𝑑𝑟)] 𝑑𝜃
4𝜇 0 𝑅 𝑅 𝑅 𝑅

𝜌𝑔𝑅 2 2𝜋 𝑎𝑅 1 𝑎𝑅 𝑎𝑅
𝜔=𝜌∗ ∫ [∫ 𝑟𝑑𝑟 − 2 ∫ 𝑟 3 𝑑𝑟 + 2𝑎2 ∫ (ln(𝑟)𝑟𝑑𝑟 − ln(𝑅)𝑟𝑑𝑟)] 𝑑𝜃
4𝜇 0 𝑅 𝑅 𝑅 𝑅

𝜌𝑔𝑅 2 2𝜋 𝑎𝑅 1 𝑎𝑅 𝑎𝑅 𝑎𝑅
𝜔=𝜌∗ ∫ [∫ 𝑟𝑑𝑟 − 2 ∫ 𝑟 3 𝑑𝑟 + 2𝑎2 [∫ ln(𝑟)𝑟𝑑𝑟 − ln(𝑅) ∫ 𝑟𝑑𝑟]] 𝑑𝜃
4𝜇 0 𝑅 𝑅 𝑅 𝑅 𝑅

Se realiza cada integral del corchete:

𝑎𝑅 𝑎𝑅 𝑎𝑅
𝜌𝑔𝑅 2 2𝜋 𝑟 2 1 𝑟4 2
𝑎𝑅
𝑟2
𝜔=𝜌∗ ∫ [[ ] − 2 [ ] + 2𝑎 [∫ ln(𝑟)𝑟𝑑𝑟 − ln(𝑅) [ ] ]] 𝑑𝜃
4𝜇 0 2 𝑅 𝑅 4 𝑅 𝑅 2 𝑅

Esta integral que queda es por partes por tanto hay que plantearla

∫ 𝑟 ln(𝑟)𝑑𝑟 = 𝑢𝑣 − ∫ 𝑣 𝑑𝑢

Donde
𝑑𝑟
𝑢 = ln(𝑟) 𝑑𝑢 =
𝑟
𝑟2
𝑑𝑣 = 𝑟 𝑑𝑟 𝑣=
2
Por tanto:
 𝑟2 𝑟 2 𝑑𝑟
∫ 𝑟 ln(𝑟)𝑑𝑟 = ln(𝑟) − ∫
2 2 𝑟
𝑟2 1
∫ 𝑟 ln(𝑟)𝑑𝑟 = ln(𝑟) − ∫ 𝑟𝑑𝑟
2 2
𝑟2 1 𝑟2
∫ 𝑟 ln(𝑟)𝑑𝑟 = ln(𝑟) −
2 22
𝑟2 𝑟2
∫ 𝑟 ln(𝑟)𝑑𝑟 = ln(𝑟) −
2 4
Se reemplaza y se evalúa para cumplir la integral definida:

𝑎𝑅 𝑎𝑅 𝑎𝑅 𝑎𝑅
𝜌𝑔𝑅 2 2𝜋 𝑟 2 1 𝑟4 𝑟2 𝑟2 𝑟2
𝜔=𝜌∗ ∫ [[ ] − 2 [ ] + 2𝑎2 [[ ln(𝑟) − ] − ln(𝑅) [ ] ]] 𝑑𝜃
4𝜇 0 2 𝑅 𝑅 4 𝑅 2 4 𝑅 2 𝑅

Se evalúan todos factores del corchete

𝜌𝑔𝑅 2 2𝜋 (𝑎𝑅)2 𝑅 2 1 (𝑎𝑅)4 𝑅 4 (𝑎𝑅)2 (𝑎𝑅)2 𝑅 2 𝑅2 (𝑎𝑅)2 𝑅 2
𝜔=𝜌∗ ∫ [( − )− 2( − ) + 2𝑎2 [( ln(𝑎𝑅) − − ln(𝑅) + ) − ln(𝑅) ( − )]] 𝑑𝜃
4𝜇 0 2 2 𝑅 4 4 2 4 2 4 2 2

Se realiza la segunda integral

𝜌𝑔𝑅 2 (𝑎𝑅)2 𝑅 2 1 (𝑎𝑅)4 𝑅 4 (𝑎𝑅)2 (𝑎𝑅)2 𝑅 2 𝑅2 (𝑎𝑅)2 𝑅 2 2𝜋
𝜔 =𝜌∗ [( − )− 2( − ) + 2𝑎2 [( ln(𝑎𝑅) − − ln(𝑅) + ) − ln(𝑅) ( − )]] ∫ 𝑑𝜃
4𝜇 2 2 𝑅 4 4 2 4 2 4 2 2 0

𝜌𝑔𝑅 2 (𝑎𝑅)2 𝑅 2 1 (𝑎𝑅)4 𝑅 4 (𝑎𝑅)2 (𝑎𝑅)2 𝑅 2 𝑅2 (𝑎𝑅)2 𝑅 2

𝜔=𝜌∗ [( − )− 2( − ) + 2𝑎2 [( ln(𝑎𝑅) − − ln(𝑅) + ) − ln(𝑅) ( − )]] [𝜃]2𝜋
0
4𝜇 2 2 𝑅 4 4 2 4 2 4 2 2

𝜌𝑔𝑅 2 (𝑎𝑅)2 𝑅 2 1 (𝑎𝑅)4 𝑅 4 (𝑎𝑅)2 (𝑎𝑅)2 𝑅 2 𝑅2 (𝑎𝑅)2 𝑅 2

𝜔 =𝜌∗ [( − )− 2( − ) + 2𝑎2 [( ln(𝑎𝑅) − − ln(𝑅) + ) − ln(𝑅) ( − )]] ∗ 2𝜋
4𝜇 2 2 𝑅 4 4 2 4 2 4 2 2

Hasta aquí ya se puede dejar, pero la idea es simplificar la ecuación

𝜌 2 𝑔𝑅 2 𝑎2 𝑅2 𝑅 2 1 𝑎4 𝑅 4 𝑅 4 𝑎2 𝑅 2 𝑎2 𝑅2 𝑅 2 𝑅2 𝑎2 𝑅2 𝑅 2
𝜔 = 2𝜋 [( − )− 2( − ) + 2𝑎2 [( ln(𝑎𝑅) − − ln(𝑅) + ) − ln(𝑅) ( − )]]
4𝜇 2 2 𝑅 4 4 2 4 2 4 2 2

𝜌 2 𝑔𝑅 2 𝑎2 𝑅 2 𝑅 2 𝑎4 𝑅4 𝑅 4 𝑎2 𝑅2 𝑎2 𝑅 2 𝑅 2 𝑅2 𝑎2 𝑅 2 𝑅2
𝜔 = 2𝜋 [ − − 2
+ 2 + 2𝑎2 [ ln(𝑎𝑅) − − ln(𝑅) + − ln(𝑅) + ln(𝑅) ]]
4𝜇 2 2 4𝑅 4𝑅 2 4 2 4 2 2

𝜌 2 𝑔𝑅 2 𝑎2 𝑅 2 𝑅 2 𝑎4 𝑅2 𝑅 2 𝑎2 𝑅 2 𝑎2 𝑅2 𝑅2 𝑅2 𝑎2 𝑅2 𝑅2
𝜔 = 2𝜋 [ − − + + 2𝑎2 ln(𝑎𝑅) − 2𝑎2 − 2𝑎2 ln(𝑅) + 2𝑎2 − ln(𝑅)2𝑎2 + ln(𝑅)2𝑎2 ]
4𝜇 2 2 4 4 2 4 2 4 2 2

𝜌 2 𝑔𝑅 2 𝑎2 𝑅2 𝑅 2 𝑎4 𝑅2 𝑅 2 𝑎4 𝑅2 𝑎2 𝑅 2
𝜔 = 2𝜋 [ − − + + 𝑎4 𝑅 2 ln(𝑎𝑅) − − 𝑎2 𝑅 2 ln(𝑅) + − 𝑎4 𝑅2 ln(𝑅) + 𝑎2 𝑅2 ln(𝑅)]
4𝜇 2 2 4 4 2 2

Se suman y se restan los valores semejantes

𝜌 2 𝑔𝑅 2 𝑎2 𝑅2 𝑅 2 𝑎4 𝑅2 𝑅 2 𝑎4 𝑅2 𝑎2 𝑅 2
𝜔 = 2𝜋 [ − − + + 𝑎4 𝑅 2 ln(𝑎𝑅) − − 𝑎2 𝑅 2 ln(𝑅) + − 𝑎4 𝑅2 ln(𝑅) + 𝑎2 𝑅2 ln(𝑅)]
4𝜇 2 2 4 4 2 2
 𝜌2 𝑔𝑅 2 𝑎2 𝑅 2 𝑎2 𝑅 2 𝑅 2 𝑎4 𝑅 2 𝑅 2 𝑎4 𝑅 2
𝜔 = 2𝜋 [ + − − + + 𝑎4 𝑅2 ln(𝑎𝑅) − − 𝑎4 𝑅 2 ln(𝑅)]
4𝜇 2 2 2 4 4 2

𝜌2 𝑔𝑅 2 2 2 𝑅 2 𝑎4 𝑅 2 𝑅 2 𝑎4 𝑅 2
𝜔 = 2𝜋 [𝑎 𝑅 − − + + 𝑎4 𝑅 2 ln(𝑎𝑅) − − 𝑎4 𝑅2 ln(𝑅)]
4𝜇 2 4 4 2

Se saca factor común 𝑅 2

𝜌2 𝑔𝑅 4 2 1 𝑎4 1 𝑎4
𝜔 = 2𝜋 [𝑎 − − + + 𝑎4 ln(𝑎𝑅) − − 𝑎4 ln(𝑅)]
4𝜇 2 4 4 2

Se agrupan los términos

𝜌2 𝑔𝑅 4 1 1 𝑎4 𝑎4
𝜔 = 2𝜋 [ − + 𝑎2 − − + 𝑎4 ln(𝑎𝑅) − 𝑎4 ln(𝑅)]
4𝜇 4 2 4 2

𝜌2 𝑔𝑅 4 2 1 3𝑎4
𝜔=𝜋 [𝑎 − − + 𝑎4 ln(𝑎𝑅) − 𝑎4 ln(𝑅)]
2𝜇 4 4

Se saca factor común ¼

𝜌2 𝑔𝑅 4
𝜔=𝜋 [4𝑎2 − 1 − 3𝑎4 + 4𝑎4 (ln(𝑎𝑅) − ln(𝑅))]
8𝜇
Se aplica la propiedad de resta de logaritmos

𝜌2 𝑔𝑅 4 𝑎𝑅
𝜔=𝜋 [4𝑎2 − 1 − 3𝑎4 + 4𝑎4 ln ( )]
8𝜇 𝑅
𝜌2 𝑔𝑅 4
𝜔=𝜋 [4𝑎2 − 1 − 3𝑎4 + 4𝑎4 ln(𝑎)]
8𝜇
𝜌2 𝑔𝑅 4
𝜔=𝜋 [((2𝑎)2 − 1) + 𝑎4 (4 ln(𝑎) − 3)] (𝑏)
8𝜇
Hasta aquí se terminaría de simplificar

c) Ahora se plantea que 𝑎 debe ser casi 1 ya que en este punto tendría el valor del radio del
tubo R, como debe ser un valor cercano a 1 se plantea que 𝑎 = 1 + 𝜎 , donde 𝜎 es un valor
muy pequeño ya que 𝑎 debe ser casi 1. Por tanto, se reemplaza:

𝜌2 𝑔𝑅 4
𝜔=𝜋 [4(1 + 𝜎)2 − 1 − 3(1 + 𝜎)4 + 4(1 + 𝜎)4 ln(1 + 𝜎)]
8𝜇
Se expanden todos los binomios y el Ln por series de Taylor

𝜌2 𝑔𝑅 4
𝜔=𝜋 [4(1 + 2𝜎 + 𝜎 2 ) − 1 − 3(𝜎 4 + 4𝜎 3 + 6𝜎 2 + 4𝜎 + 1) + 4(𝜎 4 + 4𝜎 3 + 6𝜎 2 + 4𝜎
8𝜇
𝜎2 𝜎3
+ 1)(𝜎 − + + ⋯ . . )]
2 3
 𝜌2 𝑔𝑅 4
𝜔=𝜋 [4 + 8𝜎 + 4𝜎 2 − 1 − 3𝜎 4 − 12𝜎 3 − 18𝜎 2 − 12𝜎 − 3 + (4𝜎 4 + 16𝜎 3 + 24𝜎 2
8𝜇
𝜎2 𝜎3
+ 16𝜎 + 4)(𝜎 − + + ⋯ . . )]
2 3

De la multiplicación del ultimo termino solo se mantienen los que tengan potencia menor a
4, ya que el orden de magnitud es muy pequeño aun, de igual forma los ya presentes en el
paréntesis se despreciaran por lo dicho con anterioridad:

𝜌2 𝑔𝑅 4 16
𝜔=𝜋 [4 + 8𝜎 + 4𝜎 2 − 1 − 12𝜎 3 − 18𝜎 2 − 12𝜎 − 3 + 24𝜎 3 + 16𝜎 2 + 4𝜎 − 𝜎 3
8𝜇 2
4 2 4 3
− 𝜎 + 𝜎 ]
2 3
Se suman o restan los semejantes:

𝜌2 𝑔𝑅 4 4
𝜔=𝜋 [−12𝜎 3 + 24𝜎 3 − 8𝜎 3 + 𝜎 3 ]
8𝜇 3
𝜌2 𝑔𝑅 4 16 3
𝜔=𝜋 [ 𝜎 ]
8𝜇 3
𝜌2 𝑔𝑅 4 16 3
𝜔=𝜋 [ 𝜎 ]
8𝜇 3
2𝜋𝜌2 𝑔𝑅 4 3
𝜔= 𝜎
3𝜇
Para que esta ecuación se semejante a:

Se debe plantear que 𝑊 = 2𝜋𝑅, 𝛿 3 = 𝑅 3 𝜎 3 , cos 𝛽 = 1. Planteando todo lo ultimo las

ecuaciones serán semejantes.

Ejercicio 2B-7
-Presunciones

Estado estacionario

Flujo laminar

Fluido Newtoniano e incompresible

Coordenadas cilíndricas (r, 𝜃, 𝑧)

-Velocidades

𝑣𝑟 = 0
𝑣𝜃 = 0
𝑣𝑧 = 𝑣𝑧 (𝑟, 𝜃, 𝑧, 𝑡)
-Volumen de Control

Vc=2𝜋𝑟∆𝑟𝐿

-Balance de momentum

2𝜋𝑟𝐿(𝜏𝑟𝑧 |𝑟) − 2𝜋𝑟𝐿(𝜏𝑟𝑧 |𝑟 + ∆𝑟) + 2𝜋𝑟∆𝑟((𝜌𝑣𝑧 𝑣𝑧 |𝑧 = 0) − (𝜌𝑣𝑧 𝑣𝑧 |𝑧 = 𝐿)) + 2𝜋𝑟∆𝑟(𝑝0 − 𝑝𝐿 )

− 2𝜋𝑟∆𝑟𝐿 𝜌𝑔 = 0
Se cancela el termino de convección

2𝜋𝑟𝐿(𝜏𝑟𝑧 |𝑟) − 2𝜋𝑟𝐿(𝜏𝑟𝑧 |𝑟 + ∆𝑟) + 2𝜋𝑟∆𝑟((𝜌𝑣𝑧 𝑣𝑧 |𝑧 = 0) − (𝜌𝑣𝑧 𝑣𝑧 |𝑧 = 𝐿)) + 2𝜋𝑟∆𝑟(𝑝0 − 𝑝0 )

=0
Se divide la ecuación entre las constantes y el diferencial del Volumen de control
2𝜋𝑟𝐿 2𝜋𝑟𝐿 2𝜋𝑟∆𝑟
(𝜏𝑟𝑧 |𝑟) − (𝜏𝑟𝑧 |𝑟 + ∆𝑟) + (𝑝 − 𝑝0 ) = 0
2𝜋∆𝑟𝐿 2𝜋∆𝑟𝐿 2𝜋∆𝑟𝐿 0
 (𝑟𝜏𝑟𝑧 |𝑟) − (𝑟𝜏𝑟𝑧 |𝑟 + ∆𝑟)
=0
∆𝑟
Se aplica el lim para obtener definición de derivada
(𝑟𝜏𝑟𝑧 |𝑟) − (𝑟𝜏𝑟𝑧 |𝑟 + ∆𝑟)
lim =0
∆𝑟→0 ∆𝑟
𝑑(𝑟𝜏𝑟𝑧 )
− =0
𝑑𝑟
𝑑(𝑟𝜏𝑟𝑧 )
=0
𝑑𝑟
Se aplica el método de variables separables para solucionar la EDO

∫ 𝑑(𝑟𝜏𝑟𝑧 ) = ∫ 0 𝑑𝑟

𝑟𝜏𝑟𝑧 = 𝐶1
𝐶1
𝜏𝑟𝑧 = (1)
𝑟
Se reemplaza la ley de viscosidad de Newton
𝑑𝑣𝑧 𝐶1
−𝜇 =
𝑑𝑟 𝑟
Se aplica variables separables nuevamente
𝐶1
𝑑𝑣𝑧 = − 𝑑𝑟
𝜇𝑟
𝐶1
Se plantea que 𝑘1 = −
𝜇

𝑘1
𝑑𝑣𝑧 = 𝑑𝑟
𝑟
𝑑𝑟
∫ 𝑑𝑣𝑧 = 𝑘1 ∫
𝑟
𝑣𝑧 = 𝑘1 ln 𝑟 + 𝐶2
𝑣𝑧 = 𝑘1 ln 𝑟 + 𝐶2 (2)
Se plantea las condiciones límite para calcular las constante k1 y C2

𝑟 = 𝑘𝑅 , 𝑣𝑧 = 𝑣0
𝑟=𝑅 , 𝑣𝑧 = 0
Se reemplazan ambas condiciones limites en (2)

𝑣0 = 𝑘1 ln 𝑘𝑅 + 𝐶2 (3)
 0 = 𝑘1 ln 𝑅 + 𝐶2 (4)
Se restan (3)-(4)

𝑣0 = 𝑘1 ln 𝑘𝑅 − 𝑘1 ln 𝑅 + 𝐶2 − 𝐶2
𝑣0 = 𝑘1 ln 𝑘𝑅 − 𝑘1 ln 𝑅
𝑣0 = 𝑘1 (ln 𝑘𝑅 − ln 𝑅)
Se aplica propiedad de resta de logaritmos
𝑘𝑅
𝑣0 = 𝑘1 ln ( )
𝑅
𝑣0 = 𝑘1 ln 𝑘
𝑣0
= 𝑘1
ln 𝑘
Se reemplaza k1 en (4)
𝑣0
0= ln 𝑅 + 𝐶2
ln 𝑘
𝑣0
− ln 𝑅 = 𝐶2
ln 𝑘
Se reemplazan C2 y k1 en (2)
𝑣0 𝑣0
𝑣𝑧 = ln 𝑟 − ln 𝑅
ln 𝑘 ln 𝑘
Se ajusta la ecuación
𝑣0
𝑣𝑧 = (ln 𝑟 − ln 𝑅)
ln 𝑘
𝑣0 𝑟
𝑣𝑧 = ln
ln 𝑘 𝑅
𝑟
𝑣𝑧 ln 𝑅
= 𝑎)
𝑣0 ln 𝑘
Se calcula el flujo masico planteando
2𝜋 𝑅
𝜔 =𝜌∗∫ ∫ 𝑣𝑧 𝑟𝑑𝑟𝑑𝜃
0 𝑘𝑅

Se reemplaza vz
2𝜋 𝑅
𝑣0
𝜔 =𝜌∗∫ ∫ [ (ln 𝑟 − ln 𝑅)] 𝑟𝑑𝑟𝑑𝜃
0 𝑘𝑅 ln 𝑘

𝑣0 2𝜋 𝑅
𝜔=𝜌∗ ∫ ∫ [(ln 𝑟 − ln 𝑅)] 𝑟𝑑𝑟𝑑𝜃
ln 𝑘 0 𝑘𝑅
Se distribuye rdr

𝑣0 2𝜋 𝑅
𝜔=𝜌∗ ∫ ∫ [(𝑟 ln 𝑟 𝑑𝑟 − 𝑟 ln 𝑅 𝑑𝑟)] 𝑑𝜃
ln 𝑘 0 𝑘𝑅

Se distribuye la primera integral

𝑣0 2𝜋 𝑅 𝑅
𝜔=𝜌∗ ∫ [∫ 𝑟 ln 𝑟 𝑑𝑟 − ∫ 𝑟 ln 𝑅 𝑑𝑟] 𝑑𝜃
ln 𝑘 0 𝑘𝑅 𝑘𝑅

𝑣0 2𝜋 𝑅 𝑅
𝜔=𝜌∗ ∫ [∫ 𝑟 ln 𝑟 𝑑𝑟 − ln 𝑅 ∫ 𝑟 𝑑𝑟] 𝑑𝜃
ln 𝑘 0 𝑘𝑅 𝑘𝑅

La primera integral ya se desarrollo en el ejercicio anterior por tanto la reemplazamos aquí

también:
𝑅 𝑅
𝑣0 2𝜋 𝑟 2 𝑟2 𝑟2
𝜔 =𝜌∗ ∫ [[ ln(𝑟) − ] − ln 𝑅 [ ] ] 𝑑𝜃
ln 𝑘 0 2 4 𝑘𝑅 2 𝑘𝑅

𝑣0 2𝜋 𝑅 2 𝑅 2 (𝑘𝑅)2 (𝑘𝑅)2 𝑅2 (𝑘𝑅)2

𝜔=𝜌∗ ∫ [( ln(𝑅) − − ln(𝑘𝑅) + ) − ln 𝑅 ( − )] 𝑑𝜃
ln 𝑘 0 2 4 2 4 2 2

𝜌𝑣0 2𝜋 𝑅 2 𝑅 2 (𝑘𝑅)2 (𝑘𝑅)2 𝑅2 (𝑘𝑅)2

𝜔= ∫ [( ln(𝑅) − − ln(𝑘𝑅) + ) − ln 𝑅 ( − )] 𝑑𝜃
ln 𝑘 0 2 4 2 4 2 2

Se realiza la segunda integral

𝜌𝑣0 𝑅 2 𝑅 2 (𝑘𝑅)2 (𝑘𝑅)2 𝑅2 (𝑘𝑅)2

𝜔= [( ln(𝑅) − − ln(𝑘𝑅) + ) − ln 𝑅 ( − )] [𝜃]2𝜋
0
ln 𝑘 2 4 2 4 2 2

𝜌𝑣0 𝑅2 𝑅 2 (𝑘𝑅)2 (𝑘𝑅)2 𝑅 2 (𝑘𝑅)2

𝜔= [( ln(𝑅) − − ln(𝑘𝑅) + ) − ln 𝑅 ( − )] ∗ 2𝜋
ln 𝑘 2 4 2 4 2 2

𝜌𝑣0 𝑅 2 𝑅 2 (𝑘𝑅)2 (𝑘𝑅)2 𝑅2 (𝑘𝑅)2

𝜔 = 2𝜋 [( ln(𝑅) − − ln(𝑘𝑅) + ) − ln 𝑅 ( − )]
ln 𝑘 2 4 2 4 2 2

Se distribuye el ln en el corchete:

𝜌𝑣0 𝑅 2 𝑅 2 (𝑘𝑅)2 (𝑘𝑅)2 𝑅 2 (𝑘𝑅)2

𝜔 = 2𝜋 [ ln(𝑅) − − ln(𝑘𝑅) + − ln 𝑅 + ln 𝑅]
ln 𝑘 2 4 2 4 2 2

Se restan o suman los semejantes

𝜌𝑣0 𝑅 2 𝑅 2 (𝑘𝑅)2 (𝑘𝑅)2 𝑅 2 (𝑘𝑅)2

𝜔 = 2𝜋 [ ln(𝑅) − − ln(𝑘𝑅) + − ln 𝑅 + ln 𝑅]
ln 𝑘 2 4 2 4 2 2

𝜌𝑣0 𝑅 2 (𝑘𝑅)2 (𝑘𝑅)2 (𝑘𝑅)2

𝜔 = 2𝜋 [− − ln(𝑘𝑅) + + ln 𝑅]
ln 𝑘 4 2 4 2
 𝜌𝑣0 𝑅 2 𝑘 2 𝑅2 𝑘 2 𝑅2 𝑘 2 𝑅2
𝜔 = 2𝜋 [− − ln(𝑘𝑅) + + ln 𝑅]
ln 𝑘 4 2 4 2
𝑅2
Se saca factor común 4

𝜌𝑣0 𝑅 2
𝜔 = 2𝜋 [−1 − 2𝑘 2 ln(𝑘𝑅) + 𝑘 2 + 2𝑘 2 ln 𝑅]
ln 𝑘 4
𝜋𝜌𝑣0 𝑅 2 2
𝜔= [𝑘 − 1 − 2𝑘 2 ln(𝑘𝑅) + 2𝑘 2 ln 𝑅]
2 ln 𝑘
Se ajusta la ecuación

𝜋𝜌𝑣0 𝑅 2 2
𝜔= [𝑘 − 1 + 2𝑘 2 (ln 𝑅 − ln(𝑘𝑅))]
2 ln 𝑘
𝜋𝜌𝑣0 𝑅 2 2 𝑅
𝜔= [𝑘 − 1 + 2𝑘 2 ln ]
2 ln 𝑘 𝑘𝑅
𝜋𝜌𝑣0 𝑅2 2 1
𝜔= [𝑘 − 1 + 2𝑘 2 ln ]
2 ln 𝑘 𝑘
1
Se plantea que ln 𝑘 = ln 1 − ln 𝑘 = − ln 𝑘

𝜋𝜌𝑣0 𝑅 2 2
𝜔= [𝑘 − 1 − 2𝑘 2 ln 𝑘]
2 ln 𝑘
Se distribuye el ln k que se encuentra fuera del corchete

𝜋𝜌𝑣0 𝑅 2 𝑘 2 − 1 2𝑘 2 ln 𝑘
𝜔= [ − ]
2 ln 𝑘 ln 𝑘

𝜋𝜌𝑣0 𝑅 2 𝑘 2 − 1
𝜔= [ − 2𝑘 2 ]
2 ln 𝑘

Se ajusta la ecuación para que quede como se nos pide

𝜋𝜌𝑣0 𝑅 2 1 − 𝑘 2
𝜔= [− − 2𝑘 2 ]
2 ln 𝑘

𝜋𝜌𝑣0 𝑅 2 1 − 𝑘 2
𝜔= [ − 2𝑘 2 ]
2 − ln 𝑘

Por la relación del logaritmo que se planteó arriba se obtiene:

𝜋𝜌𝑣0 𝑅2 1 − 𝑘 2
𝜔= [ − 2𝑘 2 ]
2 1
ln
𝑘

𝜋𝑅 2 𝑣0 𝜌 1 − 𝑘 2
𝜔= [ − 2𝑘 2 ] 𝑏)
2 1
ln
𝑘
Ahora se hallará la función para el esfuerzo
𝐶1
𝜏𝑟𝑧 =
𝑟
Donde
𝐶1
− = 𝑘1
𝜇
𝐶1 = −𝜇𝑘1
𝜇𝑣0
𝐶1 = −
ln 𝑘
Por tanto
𝜇𝑣0 1
𝜏𝑟𝑧 = −
ln 𝑘 𝑟
Para calcular la fuerza que ejerce el fluido sobre la varilla se plantea:

𝐹𝑧 = 𝐴(−𝜏𝑟𝑧 |𝑟 = 𝑘𝑅)
El esfuerzo es negativo debido a que al moverse la varilla el momentum de esta supera al del
fluido, por tanto, el fluido ejercera un esfuerzo negativo sobre esta:
𝜇𝑣0 1
𝐹𝑧 = −2𝜋𝑘𝑅𝐿 (− )
ln 𝑘 𝑘𝑅
2𝜋𝐿𝜇𝑣0
𝐹𝑧 =
ln 𝑘
1
Se plantea que ln = − ln 𝑘
𝑘

2𝜋𝐿𝜇𝑣0
𝐹𝑧 =
1
− ln
𝑘
2𝜋𝐿𝜇𝑣0
𝐹𝑧 = − 𝑐)
1
ln
𝑘
Se plantea que 𝑘 = 1 − 𝜖, y reemplazándolo en la fórmula de Fuerza
2𝜋𝐿𝜇𝑣0
𝐹𝑧 =
ln 𝑘
2𝜋𝐿𝜇𝑣0
𝐹𝑧 =
ln 1 − 𝜖
Se extiende el denominado por series de Taylor:
 2𝜋𝐿𝜇𝑣0
𝐹𝑧 =
𝜖2 𝜖3 𝜖4
(−𝜖 − 2 − 3 − 4 − ⋯ . . )

En el ejercicio 2B-5 se recomienda solo usar los primeros 4 términos de la serie

2𝜋𝐿𝜇𝑣0
𝐹𝑧 =
𝜖2 𝜖3 𝜖4
(−𝜖 − 2 − 3 − 4 )

2𝜋𝐿𝜇𝑣0
𝐹𝑧 =
𝜖 𝜖2 𝜖3
−𝜖(1 + 2 + 3 + 4 )

2𝜋𝐿𝜇𝑣0
𝐹𝑧 = −
𝜖 𝜖2 𝜖3
𝜖(1 + 2 + 3 + 4 )

2𝜋𝐿𝜇𝑣0 1
𝐹𝑧 = −
𝜖 𝜖 𝜖2 𝜖3
1+2+ 3 + 4

Ahora si se expande este polinomio de la derecha por series de Taylor a 𝜖 = 0, se obtiene

2𝜋𝐿𝜇𝑣0 𝜖 𝜖2 𝜖3
𝐹𝑧 = − (1 − − − ) 𝑑)
𝜖 2 12 24

Ejercicio 2B-11

-Presunciones

Estado estacionario

Flujo laminar

Fluido Newtoniano e incompresible

Coordenadas rectangulares ( 𝑥, 𝑦, 𝑧)

-Velocidades
𝑣𝑥 = 𝑣𝑧 ( 𝑥, 𝑦, 𝑧, 𝑡)
𝑣𝑦 = 0

𝑣𝑧 = 0
-Volumen de Control

Vc=𝐴∆𝑦

-Balance de momentum

𝐴(𝜏𝑦𝑥 |𝑦) − 𝐴(𝜏𝑦𝑥 |𝑦 + ∆𝑦) + 𝑑𝑧∆𝑦((𝜌𝑣𝑥 𝑣𝑥 |𝑥 = 0) − (𝜌𝑣𝑥 𝑣𝑥 |𝑥 = 𝐿)) = 0

Se cancela el termino de convección

𝐴(𝜏𝑦𝑥 |𝑦) − 𝐴(𝜏𝑦𝑥 |𝑦 + ∆𝑦) + 𝑑𝑧∆𝑦((𝜌𝑣𝑥 𝑣𝑥 |𝑥 = 0) − (𝜌𝑣𝑥 𝑣𝑥 |𝑥 = 𝐿)) = 0

Se divide la ecuación entre las constantes y el diferencial del Volumen de control

𝐴 𝐴
(𝜏𝑦𝑥 |𝑦) − (𝜏 |𝑦 + ∆𝑦) = 0
𝐴∆𝑦 𝐴∆𝑦 𝑦𝑥

(𝜏𝑦𝑥 |𝑦) − (𝜏𝑦𝑥 |𝑦 + ∆𝑦)

=0
∆𝑦
Se aplica el lim para obtener definición de derivada

(𝜏𝑦𝑥 |𝑥) − (𝜏𝑦𝑥 |𝑥 + ∆𝑥)

lim =0
∆𝑦→0 ∆𝑦

𝑑(𝜏𝑦𝑥 )
− =𝑜
𝑑𝑦

𝑑(𝜏𝑦𝑥 )
=0
𝑑𝑦
Se aplica el método de variables separables para solucionar la EDO

∫ 𝑑(𝜏𝑦𝑥 ) = ∫ 0 𝑑𝑦

𝜏𝑦𝑥 = 𝐶1 (1)

Se reemplaza la ley de viscosidad de Newton

𝑑𝑣𝑥
−𝜇 = 𝐶1
𝑑𝑦
𝑑𝑣𝑥 𝐶1
=−
𝑑𝑦 𝜇
Se plantea que
 𝐶1
− = 𝑘1
𝜇
Por tanto
𝑑𝑣𝑥
= −𝑘1
𝑑𝑦
Se aplica variables separables nuevamente

∫ 𝑑𝑣𝑥 = 𝑘1 ∫ 𝑑𝑦

𝑣𝑥 = 𝑘1 𝑦 + 𝐶2 (2)
Se plantea las condiciones límites para hallar k1 y c2

𝑦=0 , 𝑣𝑥 = 0
𝑥=𝑏 , 𝑣𝑥 = 𝑉
Se reemplazan las condiciones limite en (2)

0 = 𝑘1 0 + 𝐶2
De este primer reemplazo se obtiene que C2 =0

𝑉 = 𝑘1 𝑏
Por tanto
𝑉
= 𝑘1
𝑏
Se reemplazan k1 y c2 en (2)
𝑉
𝑣𝑥 = 𝑦 (3)
𝑏
Se reemplazan los valores dichos por el problema para obtener la expresión de la velocidad
en coordenadas esféricas:

𝑣𝑥 = 𝑣𝜙

𝑉 = Ω𝑟
𝑏 = 𝑟 sin 𝜓0 = 𝑟𝜓0
𝜋
𝑦 = 𝑟 sin 𝜓 = 𝑟𝜓 = 𝑟 ( − 𝜃)
2
Reemplazando todos los valores en(3) obtenemos
Ω𝑟 𝜋
𝑣𝜙 = 𝑟 ( − 𝜃)
𝑟𝜓0 2
 Ω 𝜋
𝑣𝜙 = 𝑟 ( − 𝜃)
𝜓0 2
𝑣𝜙 Ω 𝜋
= ( − 𝜃) 𝑎)
𝑟 𝜓0 2
En el apéndice B1 se plantea que
sin 𝜃 𝑑 𝑣𝜙 1 𝑑
𝜏𝜃𝜙 = −𝜇 [ ( ) ] ≈ −𝜇 [ (𝑣 ) ]
𝑟 𝑑𝜃 sin 𝜃 𝑟 𝑑𝜃 𝜙
Se calcula la derivada de este término:
𝑑𝑣𝜙 𝑑 Ω 𝜋
= ( 𝑟 ( − 𝜃))
𝑑𝜃 𝑑𝜃 𝜓0 2

𝑑𝑣𝜙 Ω 𝑑 1 𝜋
= 𝑟 ( ( − 𝜃))
𝑑𝜃 𝜓0 𝑑𝜃 2
𝑑𝑣𝜙 Ω 𝑑 𝜋 𝑑𝜃
= 𝑟 [ ( )− ]
𝑑𝜃 𝜓0 𝑑𝜃 2 𝑑𝜃
𝑑𝑣𝜙 Ω
= 𝑟 [0 − 1]
𝑑𝜃 𝜓0
𝑑𝑣𝜙 Ω
=− 𝑟
𝑑𝜃 𝜓0
Se reemplaza en el esfuerzo
1 𝑑
𝜏𝜃𝜙 = −𝜇 [ (𝑣 ) ]
𝑟 𝑑𝜃 𝜙
1 Ω
𝜏𝜃𝜙 = −𝜇 [ (− 𝑟 ) ]
𝑟 𝜓0
Ω
𝜏𝜃𝜙 = 𝜇 𝑏)
𝜓0
Se calcula el torque planteando
2𝜋 𝑅
𝜋
𝑇𝑧 = ∫ ∫ (𝜏𝜃𝜙 ∙ 𝑟|𝜃 = ) 𝑟𝑑𝑟𝑑𝜙
0 0
2
2𝜋 𝑅
Ω 2
𝑇𝑧 = ∫ ∫ 𝜇 𝑟 𝑑𝑟𝑑𝜙
0 0 𝜓0

Ω 2𝜋 𝑅 2
𝑇𝑧 = 𝜇 ∫ ∫ 𝑟 𝑑𝑟𝑑𝜙
𝜓0 0 0
𝑅
Ω 2𝜋 𝑟 3
𝑇𝑧 = 𝜇 ∫ [ ] 𝑑𝜙
𝜓0 0 3 0
 Ω 2𝜋 𝑅 3
𝑇𝑧 = 𝜇 ∫ 𝑑𝜙
𝜓0 0 3

Ω 𝑅3
𝑇𝑧 = 𝜇 [𝜙]2𝜋
0
𝜓0 3
Ω 𝑅3
𝑇𝑧 = 𝜇 2𝜋
𝜓0 3
2 Ω𝑅 3
𝑇𝑧 = 𝜋𝜇 𝑐)
3 𝜓0
d) Se tienen como datos:
𝑟𝑎𝑑 9,5493°
Ω = 10 =
𝑚𝑖𝑛 𝑠
𝑘𝑔
𝜇 = 100 𝑐𝑝 = 0,1
𝑚𝑠

𝑅 = 10 𝑐𝑚 = 0,1 𝑚
𝜓0 = 0,5 °
Se reemplazan en la formula del torque y se obtiene
9,5493°
2 𝑘𝑔 (0,1 𝑚)3
𝑇𝑧 = 𝜋 ∙ 0,1 𝑠 = 0,004 𝑁𝑚 = 40000 𝑑𝑖𝑛𝑎 − 𝑐𝑚 𝑑)
3 𝑚𝑠 0,5 °