VEXPO
VEXPO
Ciclo: 5
Alumnos:
𝑘𝑅
𝜏𝑟𝑧 |𝑟=𝑘𝑅 = 0
𝑅
- Coordenadas:
𝐶𝑖𝑙í𝑛𝑑𝑟𝑖𝑐𝑎𝑠 (𝑧, 𝑟, 𝜃)
- Velocidad:
𝑣⃗(𝑣𝑧 , 𝑣𝑟 , 𝑣𝜃 )
- Variación de la velocidad:
𝑣𝑧 (𝑧, 𝑟, 𝜃)
a. Determine el perfil de velocidades 𝑣𝑧 .
- Información:
𝑣𝑟 = 𝑣𝜃 = 0
𝜕𝑣𝑧 𝜕𝑣𝑧
= =0
𝜕𝜃 𝜕𝑧
𝜕 2 𝑣𝑧 𝜕 2 𝑣𝑧
= =0
𝜕𝜃 2 𝜕𝑧 2
𝜕𝑣𝑧 𝜕 2 𝑣𝑧
≠0→ ≠0
𝜕𝑟 𝜕𝑟 2
- Ecuación de movimiento en términos de 𝜏:
1𝑑
𝜌𝑔𝑧 = − (𝑟𝜏𝑟𝑧 )
𝑟 𝑑𝑟
(−𝑟𝜌𝑔𝑧 )𝑑𝑟 = 𝑑(𝑟𝜏𝑟𝑧 )
𝑟 2 𝜌𝑔𝑧
𝑟𝜏𝑟𝑧 = + 𝑐1
2
𝑟𝜌𝑔𝑧 𝑐1
𝜏𝑟𝑧 = +
2 𝑟
- Condiciones de frontera 𝜏𝑟𝑧 |𝑟=𝑘𝑅 = 0:
𝑘𝑅𝜌𝑔𝑧 𝑐1
𝜏𝑟𝑧 (𝑘𝑅) = + =0
2 𝑘𝑅
(𝑘𝑅)2 𝜌𝑔𝑧
𝑐1 = −
2
𝑟𝜌𝑔𝑧 (𝑘𝑅)2 𝜌𝑔𝑧 1
𝜏𝑟𝑧 = −
2 2 𝑟
𝜌𝑔𝑧 (𝑘𝑅)2
𝜏𝑟𝑧 = [𝑟 − ]
2 𝑟
- Por el 𝜏𝑟𝑧 para 𝑣𝑧 :
𝜌𝑔𝑧 (𝑘𝑅)2
− [𝑟 − ] 𝑑𝑟 = 𝑑𝑣𝑧
2𝜇 𝑟
𝜌𝑔𝑧 𝜌𝑔𝑧 2
𝑣𝑧 = (𝑘𝑅)2 ln 𝑟 − 𝑟 + 𝑐2
2𝜇 4𝜇
- Condición de frontera 𝑣𝑧 |𝑟=𝑅 = 0:
𝜌𝑔𝑧 𝜌𝑔𝑧 2
𝑣𝑧 (𝑅) = (𝑘𝑅)2 ln 𝑅 − 𝑅 + 𝑐2 = 0
2𝜇 4𝜇
𝜌𝑔𝑧 𝑅 2 1
𝑐2 = [ − 𝑘 2 ln 𝑅]
2𝜇 2
𝜌𝑔𝑧 𝜌𝑔𝑧 2 𝜌𝑔𝑧 𝑅 2 1
𝑣𝑧 = (𝑘𝑅)2 ln 𝑟 − 𝑟 + [ − 𝑘 2 ln 𝑅]
2𝜇 4𝜇 2𝜇 2
𝜌𝑔𝑧 𝑟2 1
𝑣𝑧 = [(𝑘𝑅)2 ln 𝑟 − + 𝑅 2 ( − 𝑘 2 ln 𝑅)]
2𝜇 2 2
𝜌𝑔𝑧 2 2 𝑟 𝑟 2 + 𝑅2
𝑣𝑧 = [𝑅 𝑘 (ln ) − ]
2𝜇 𝑅 2
b. Determine cuál es la fuerza tangencial que hace el líquido contra la pared.
𝜌𝑔𝑧 (𝑘𝑅)2
𝐹𝑧 = [ [𝑟 − ]] ∗ 2𝜋𝑅𝐿
2 𝑟
𝜌𝑔𝑧 (𝑘𝑅)2
𝐹𝑧 = [𝑅 − ] ∗ 2𝜋𝑅𝐿
2 𝑅
𝜌𝑔𝑧
𝐹𝑧 = [𝑅 − 𝑘 2 𝑅] ∗ 2𝜋𝑅𝐿
2
𝐹𝑧 = 𝜋𝐿𝜌𝑔𝑧 𝑅 2 [1 − 𝑘 2 ]
2. Un líquido fluye entre dos tubos coaxiales de radio 𝑘𝑅 (el interior) y 𝑅 (el exterior).
El cilindro interior asciende de tal forma que ejerce una fuerza 𝐹𝑧 igual a – 𝜙 (estamos
considerando el sentido del eje z hacia abajo). No hay diferencia de presiones
apreciable que genere o consuma cantidad de movimiento.
𝜙
𝐹𝑧 = −𝜙 𝜏𝑟𝑧 (𝑘𝑅) = −
2𝜋𝑘𝑅𝐿
𝐹𝑧
𝜏𝑟𝑧 =
2𝜋𝑟𝐿
𝑣𝑧 (𝑅) = 0
- Coordenadas:
𝐶𝑖𝑙í𝑛𝑑𝑟𝑖𝑐𝑎𝑠 (𝑧, 𝜃, 𝑟)
- Velocidad:
𝑣⃗(𝑣𝑟 , 𝑣𝜃 , 𝑣𝑧 )
- Variación de la velocidad:
𝑣𝑧 (𝑧, 𝜃, 𝑟)
- Derivadas:
𝑣𝑟 = 𝑣𝜃 = 0
𝜕𝑣𝑧 𝜕𝑣𝑧
= =0
𝜕𝜃 𝜕𝑧
𝜕 2 𝑣𝑧 𝜕 2 𝑣𝑧
= =0
𝜕𝜃 2 𝜕𝑧 2
𝜕𝑣𝑧 𝜕 2 𝑣𝑧
≠0→ ≠0
𝜕𝑟 𝜕𝑟 2
a. Determine el perfil del esfuerzo cortante 𝜏𝑟𝑧
1𝑑
(𝑟𝜏𝑟𝑧 ) = 𝜌𝑔𝑧
𝑟 𝑑𝑟
𝑑(𝑟𝜏𝑟𝑧 ) = 𝑟𝜌𝑔𝑧 𝑑𝑟
𝑟 2 𝜌𝑔𝑧
𝑟𝜏𝑟𝑧 = + 𝑐1
2
𝑟𝜌𝑔𝑧 𝑐1
𝜏𝑟𝑧 = +
2 𝑟
- Por la condición de frontera:
𝑘𝑅𝜌𝑔𝑧 𝑐1 𝜙
𝜏𝑟𝑧 (𝑘𝑅) = + =−
2 𝑘𝑅 2𝜋𝑘𝑅𝐿
𝜙 (𝑘𝑅)2 𝜌𝑔𝑧
𝑐1 = − −
𝜋𝐿 2
𝑟𝜌𝑔𝑧 𝜙 (𝑘𝑅)2 𝜌𝑔𝑧 1
𝜏𝑟𝑧 = + [− − ]
2 𝜋𝐿 2 𝑟
𝜌𝑔𝑧 (𝑘𝑅)2 𝜙
𝜏𝑟𝑧 = [𝑟 − ]−
2 𝑟 𝜋𝐿𝑟
𝜙 𝜌𝑔𝑧 𝑟2 1
𝑣𝑧 = [ln 𝑟/𝑅] + [(𝑘𝑅) ln 𝑟 − + 𝑅 2 ( − 𝑘 2 ln 𝑅)]
2
𝜋𝐿𝜇 2𝜇 2 2
3. Un líquido fluye entre dos tubos coaxiales de radio kR (el interior) y R (el exterior).
a. Determine el perfil del esfuerzo cortante 𝜏𝑟𝑧
b. Determine el perfil de velocidades 𝑣𝑧
𝑧=0 𝐿 𝑧=𝐿
𝑣𝑧 (𝑅) = 0
𝑃0 𝑃1 𝑅
𝑘𝑅 𝑣𝑧 (𝑘𝑅) = 0
- Coordenadas:
𝐶𝑖𝑙í𝑛𝑑𝑟𝑖𝑐𝑎𝑠 (𝑧, 𝜃, 𝑟)
- Velocidad:
𝑣⃗(𝑣𝑟 , 𝑣𝜃 , 𝑣𝑧 )
- Variación de la velocidad:
𝑣𝑧 (𝑧, 𝜃, 𝑟)
- Derivadas:
𝑣𝑟 = 𝑣𝜃 = 0
𝜕𝑣𝑧 𝜕𝑣𝑧
= =0
𝜕𝜃 𝜕𝑧
𝜕 2 𝑣𝑧 𝜕 2 𝑣𝑧
= =0
𝜕𝜃 2 𝜕𝑧 2
𝜕𝑣𝑧 𝜕 2 𝑣𝑧
≠0→ ≠0
𝜕𝑟 𝜕𝑟 2
𝑝0 − 𝑝1 𝑅2 1 − 𝑘 2
𝜏𝑟𝑧 =( ) [𝑟 + ( )]
2𝐿 2𝑟 ln 𝑘
𝑝0 − 𝑝1 𝑐1 𝑑𝑣𝑧
( ) 𝑟 + = −𝜇 [ ]
2𝐿 𝑟 𝑑𝑟
𝑝0 − 𝑝1 𝑐1
𝑑𝑣𝑧 = [− ( ) 𝑟 − ] 𝑑𝑟
2𝜇𝐿 𝑟𝜇
𝑝0 − 𝑝1 2 𝑐1
𝑣𝑧 = − ( ) 𝑟 − ln 𝑟 + 𝑐2
4𝜇𝐿 𝜇
- Condiciones de frontera:
𝑝0 − 𝑝1 2 𝑐1
𝑣𝑧 (𝑅) = − ( ) 𝑅 − ln 𝑅 + 𝑐2 = 0
4𝜇𝐿 𝜇
𝑝0 − 𝑝1 𝑐1
𝑣𝑧 (𝑘𝑅) = − ( ) (𝑘𝑅)2 − ln 𝑘𝑅 + 𝑐2 = 0
4𝜇𝐿 𝜇
𝑐1 𝑝0 − 𝑝1
[ln 𝑘𝑅 − ln 𝑅] + ( ) [(𝑘𝑅)2 − 𝑅 2 ] = 0
𝜇 4𝜇𝐿
𝑐1 𝑝0 − 𝑝1 [𝑅 2 − (𝑘𝑅)2 ]
=( )
𝜇 4𝜇𝐿 ln 𝑘
𝑝0 − 𝑝1 𝑝0 − 𝑝1 [𝑅 2 − (𝑘𝑅)2 ]
𝑐2 = ( ) (𝑘𝑅)2 + ( ) ln 𝑘𝑅
4𝜇𝐿 4𝜇𝐿 ln 𝑘
𝑝0 − 𝑝1 (1 − 𝑘 2 ) ln 𝑘𝑅
𝑐2 = 𝑅 2 ( ) [𝑘 2 + ]
4𝜇𝐿 ln 𝑘
𝑝0 − 𝑝1 2 𝑝0 − 𝑝1 [𝑅 2 − (𝑘𝑅)2 ] 2
𝑝0 − 𝑝1 2
(1 − 𝑘 2 ) ln 𝑘𝑅
𝑣𝑍 = − ( )𝑟 − ( ) ln 𝑟 + 𝑅 ( ) [𝑘 + ]
4𝜇𝐿 4𝜇𝐿 ln 𝑘 4𝜇𝐿 ln 𝑘
𝑝0 − 𝑝1 (1 − 𝑘 2 ) ln 𝑘𝑅 [1 − 𝑘 2 ]
𝑣𝑧 = ( ) [𝑅 2 (𝑘 2 + − ) − 𝑟 2]
4𝜇𝐿 ln 𝑘 ln 𝑘
𝒑𝟎 − 𝒑𝟏 𝟐 𝟐
(𝟏 − 𝒌𝟐 )
𝒗𝒛 = ( ) [𝑹 (𝒌 + [𝐥𝐧 𝒌𝑹 − 𝟏]) − 𝒓𝟐 ]
𝟒𝝁𝑳 𝐥𝐧 𝒌
4. Se tiene un líquido newtoniano entre dos cilindros coaxiales de radio kR (el interior) y KR
(el exterior). Éstos giran a velocidades angulares opuestas e iguales a Ω el cilindro interior y
-2 Ω el cilindro exterior.
a. ¿En qué punto la velocidad V del fluido es igual a cero?
b. Determine el valor del flujo volumétrico (caudal Q).
Resolución.
-Analizamos
𝑑 1𝑑
𝜇 ( (𝑟𝑣𝜃 )) = 0
𝑑𝑟 𝑟 𝑑𝑟
Integramos
1𝑑
𝑑( (𝑟𝑣𝜃 )) = 0
𝑟 𝑑𝑟
1𝑑
(𝑟𝑣𝜃 ) = 𝑐1
𝑟 𝑑𝑟
𝑑(𝑟𝑣𝜃 ) = 𝑐1 𝑟𝑑𝑟
Volvemos integrar
𝑐1 2
𝑟𝑣𝜃 = 𝑟 + 𝑐2
2
𝐶1 𝐶2
𝑉𝜃 = 𝑟+
2 𝑟
Condiciones de frontera
𝑟 = 𝐾𝑅 → 𝑣𝜃 = 𝛺𝐾𝑅
𝑟 = 𝑅 → 𝑣𝜃 = −2𝛺𝐾𝑅
𝐶1 𝐶2
𝛺𝑘𝑅 = 𝑘𝑅 + . . . . (1)
2 𝑘𝑅
𝐶1 𝐶2 −1
(−2𝛺𝑅 = 𝑅+ )∗ . . . (2)
2 𝑅 𝐾
2𝛺𝑅 𝐶1 𝐶1
𝛺𝑘𝑅 + = 𝑘𝑅 − 𝑅
𝑘 2 2
𝛺𝐾 2 𝑅 + 2𝛺𝑅 𝐶1 𝑅
= (𝑘𝑅 − )
𝐾 2 𝑘
𝛺𝐾 2 𝑅 + 2𝛺𝑅 𝐶1 𝑘 2 𝑅 − 𝑅
= ( )
𝐾 2 𝐾
2[𝛺𝑅(𝑘 2 + 2)] = 𝐶1 (𝑘 2 𝑅 − 𝑅)
2𝛺𝑅𝑘 2 + 4𝛺𝑅
= 𝐶1
(𝐾 2 𝑅 − 𝑅)
2𝛺𝑅(𝐾 2 + 2)
= 𝐶1
𝑅(𝑘 2 − 1)
𝟐𝜴(𝒌𝟐 + 𝟐)
= 𝑪𝟏
(𝒌𝟐 − 𝟏)
𝐶1 𝐶2
𝛺𝑘𝑅 = 𝑘𝑅 + . . . . (1)
2 𝑘𝑅
𝟐𝜴(𝒌𝟐 + 𝟐)
= 𝑪𝟏
(𝒌𝟐 − 𝟏)
2
2𝛺(𝑘 + 2)
2
(𝑘 − 1) 𝐶
𝛺𝑘𝑅 = 𝑘𝑅 + 2
2 𝑘𝑅
2
𝛺(𝑘 + 2)𝑘𝑅 𝐶2
𝛺𝑘𝑅 − =
2
(𝑘 − 1) 𝑘𝑅
2 2 𝐶2
𝛺𝑘𝑅(𝑘 − 1) − (𝛺𝑘 + 2𝛺)𝑘𝑅 =
𝑘𝑅
3 3 𝐶2
𝛺𝑘 𝑅 − 𝛺𝑘𝑅 − 𝛺𝑅𝑘 − 2𝛺𝑘𝑅 =
𝑘𝑅
2
−3𝛺(𝑘𝑅)
[ 2
] = 𝐶2
(𝑘 − 1)
𝟑𝜴(𝒌𝑹)𝟐
= 𝑪𝟐
(𝟏 − 𝒌𝟐 )
Por lo tanto, el perfil de velocidad queda de la siguiente manera , reemplazando los valores de las
constantes 1 y 2 :
2𝛺(𝑘 2 + 2) 3𝛺(𝑘𝑅)2
(𝑘 2 − 1) (1 − 𝑘 2 )
𝑉𝜃 = 𝑟+
2 𝑟
𝜴(𝒌𝟐 + 𝟐) 𝟑𝜴𝒌𝟐 𝑹𝟐
𝑽𝜽 = ∗𝒓+
(𝒌𝟐 − 𝟏) (𝟏 − 𝒌𝟐 ) ∗ 𝒓
Igualamos 𝑉𝜃 = 0
𝛺(𝑘 2 + 2) 3𝛺𝑘 2 𝑅2
𝑉𝜃 = ∗𝑟+ =0
(𝑘 2 − 1) (1 − 𝑘 2 ) ∗ 𝑟
𝑟𝛺(𝑘 2 + 2) 3𝛺𝑘 2 𝑅2
= 2
(𝑘 2 − 1) (𝑘 − 1)𝑟
2
3 𝑘 2 𝑅2
𝑟 = 2
(𝑘 + 2)
𝟑
𝒓 = 𝒌𝑹 √
(𝒌𝟐 + 𝟐)
b. Determine el valor del flujo volumétrico (caudal Q).
5. Se tiene un líquido newtoniano entre dos cilindros coaxiales de longitud L y de radios kR
(el interior) y R (el exterior). El interior gira y ejerce sobre el fluido un esfuerzo cortante
𝜏𝑟θ constante e igual a 𝜙
a. Determine el perfil del esfuerzo cortante 𝜏𝑟θ
b. Determine el perfil de velocidades 𝑉θ
Resolución
Velocidad:
𝑣 = (𝑣𝑟 , 𝑣𝜃 , 𝑣𝑧 )
Variación de la Velocidad
𝑣𝜃 = 𝑣𝜃 (𝑟, 𝜃, 𝑧)
Derivadas
𝟏 𝐝 𝟐
𝟎=[ (𝐫 𝛕𝐫𝛉 )]
𝐫 𝟐 𝐝𝐫
1 d 2
[ (r τrθ )] = 0
r 2 dr
𝑟 2 𝜏𝑟𝜃 = 𝑐1
𝑐1
𝜏𝑟𝜃 =
𝑟2
𝑐1
𝜏𝑟𝜃 (𝑘𝑅) = =𝜙
(𝑘𝑅)2
𝑐1 = 𝜙(𝑘𝑅)2
𝜙(𝑘𝑅)2
𝜏𝑟𝜃 =
𝑟2
L 𝐹𝑧 = 𝜙
W
𝑦 B
𝑟
=𝐴
𝑧
𝑉⃗
- Coordenadas:
𝑅𝑒𝑐𝑡𝑎𝑛𝑔𝑢𝑙𝑎𝑟𝑒𝑠 (𝑥, 𝑦, 𝑧)
- Velocidad:
𝑣⃗(𝑣𝑥 , 𝑣𝑦 , 𝑣𝑧 )
- Variación de la velocidad:
𝑣𝑧 (𝑥, 𝑦, 𝑧)
- Derivadas:
𝑣𝑥 = 𝑣𝑦 = 0
𝜕𝑣𝑧 𝜕𝑣𝑧
= =0
𝜕𝑦 𝜕𝑧
𝜕 2 𝑣𝑧 𝜕 2 𝑣𝑧
= =0
𝜕𝑦 2 𝜕𝑧 2
𝜕𝑣𝑧 𝜕 2 𝑣𝑧
≠0→ ≠0
𝜕𝑥 𝜕𝑥 2
𝑑 2 𝑣𝑍
𝜇[ 2 ] = 0
𝑑𝑥
𝑑𝑣𝑧
= 𝑐1
𝑑𝑥
𝑣𝑧 = 𝑐1 𝑥 + 𝑐2
- Por la condición de frontera:
𝑣𝑧 (0) = 𝑐1 (0) + 𝑐2 = 𝑉
𝑐2 = 𝑉
𝑣𝑧 = 𝑐1 𝑥 + 𝑉
- En la ecuación de Newton:
𝑑𝑣𝑧
𝜏𝑥𝑧 = −𝜇 [ ]
𝑑𝑥
𝜏𝑥𝑧 = −𝜇𝑐1
- Condición de frontera
𝜙
𝜏𝑥𝑧 (𝐵) = = −𝜇𝑐1
𝑊𝐿
𝜙
𝑐1 = −
𝜇𝑊𝐿
𝜙
𝑣𝑧 = − 𝑥+𝑉
𝜇𝑊𝐿
- Para el caudal:
𝐵 𝑊
𝑄 = ∫ ∫ 𝑣𝑧 ∙ 𝑑𝑦 𝑑𝑥
0 0
𝐵 𝑊
𝜙
𝑄 = ∫ ∫ (− 𝑥 + 𝑉) ∙ 𝑑𝑦 𝑑𝑥
0 0 𝜇𝑊𝐿
𝐵
𝜙
𝑄 = ∫ (− 𝑥𝑊 + 𝑉𝑊) 𝑑𝑥
0 𝜇𝑊𝐿
𝐵2𝜙
𝑄=− + 𝑉𝐵𝑊
𝜇2𝐿
b. Determine el valor de la velocidad promedio.
𝑄
< 𝑣𝑧 > = 𝐵 𝑊
∫0 ∫0 𝑑𝑥𝑑𝑦
𝐵2𝜙
− 𝜇2𝐿 + 𝑉𝐵𝑊
< 𝑣𝑧 > =
𝑊𝐵
𝐵𝜙
< 𝑣𝑧 > = 𝑉 −
𝜇2𝑊𝐿