ACTIVIDAD ENTREGABLE TUBERIAS
Comportamiento de los fluidos al interior de las tuberías
R. Aguas, D. Hernández, L. Iguaran, A. Villacís
[email protected] ,
[email protected],
Iviguaranuninorte.edu.co ,
[email protected] EJERCICIO 1
En la figura mostrada a continuación se puede observar el comportamiento del flujo
en las válvulas de salida (agujas) antes de impactar la turbina en la central hidroeléctrica
Hidro Sogamoso.
Si el perfil de velocidad tiene forma parabólica. Si la velocidad máxima, V=25m/s,
y el perfil está descrito por:
𝑢 = 𝑘𝑏 (𝑏𝑜 − 𝑏) 𝑘: constante
Si el diámetro D es de 15mm, calcular la velocidad máxima, sabiendo que V es la velocidad
uniforme antes de sentir el efecto de la aguja en la distribución de velocidad.
Desarrollo del Problema
Dada la información proporcionada en el problema, donde se establece que la
velocidad uniforme antes de que la aguja tenga efecto es V = 25 m/s, se puede inferir que
esta velocidad también representa la velocidad promedio del sistema, según lo ilustrado en
la Figura 1.
Figura.1. Desarrollo de la capa límite de velocidad en una tubería
Fuente: Mecánica de fluidos Cengel 4ta Edición
En particular, tal y como se ilustra en la Figura 1, se observa que la velocidad
promedio permanece constante a lo largo de toda la sección de la tubería. Sin embargo, se
producen variaciones en las velocidades dependiendo de la ubicación del fluido con
respecto al tubo. En este contexto, las paredes del tubo muestran velocidades mínimas
(aproximadamente cero), las cuales aumentan gradualmente hacia el centro del tubo, donde
se espera que se alcance la velocidad máxima. Este fenómeno se describe mediante el perfil
de velocidades del flujo. Para este caso, dicho perfil se modela con la ecuación (1), que
describe una parábola. A continuación, se detalla el cálculo asociado con la determinación
de la ubicación donde la velocidad alcanza su máximo, empleando un proceso de
maximización que implica la derivada de la ecuación que modela el perfil.
𝑢 = 𝑘𝑏 (𝑏𝑜 − 𝑏) 𝑘: constante (1)
Para hallar la ubicación del punto de mayor velocidad se tiene:
du
=0
db
d
db
( kb ( b0−b ) )
d
k
db
( b b 0−b2 )
k b 0−2 bk=0
k b 0=2 bk
b0
b= (2)
2
Así, se determina el punto donde se albergan las velocidades máximas en todo el
sistema, lo cual se destaca por coincidir particularmente con el eje centroidal del tubo.
Seguidamente, una vez identificada esta ubicación máxima, es posible sustituir su valor
encontrado en la ecuación (2) para obtener el valor asociado a la velocidad máxima del
fluido. Reemplazando (2) en (1) se tiene:
( ) ( )
2
b0 b0 b0 b0 b0
U max =k b 0− =k =k (3)
2 2 2 2 4
2
b0
Con esto se conoce que la velocidad máxima equivale a k veces y sabiendo que
4
K es una constante es posible utilizar una ecuación para calcular la velocidad promedio y
obtener el valor de K en función de b 0.
b0
1 25 m
V promedio = ∫
b0 −0 0
kb ( b0−b ) db=
s
(4 )
Por efectos prácticos siendo bo = a
a
1
a0
∫ kb ( a−b ) db=25 m/s
k
¿
a
a a
k
V = a ∫ bdb−∫ b db
2
a 0 0
[ ] [ ]
a a
b2 b3
k −
a 0 3 0
2 3
a k a
k[ ]− [ ]
2 a 3
2 2
k a 2k a
3 − =25 m/s
2∗3 3∗2
2 2 2
3 k a −2 k a k a
=
6 6
2
ka
=25 m/s
6
(25 m/s )(6)
k=
a2
Entonces restaurando la variable se tiene que:
150 m/ s
k= 2
(5)
bo
Reemplazando (5) en (3):
2
b0
U max =k
4
( )
m
150
s 2
2
∗ b0
bo
U max =
4
m
U max =37 , 5
s
Es crucial destacar que el resultado de velocidad obtenido induce a pesar que el
flujo es turbulento, ya que, si el flujo estuviera en un régimen completamente laminar y
asumiendo que está completamente desarrollado, es decir, ha superado la región de entrada,
la velocidad máxima del perfil correspondería a dos veces la velocidad promedio, en este
caso 50 m/s. Para obtener más información al respecto (Véase capítulo sobre flujo en
tuberías del libro de Cengel 4ta edición)
Por lo tanto, asumiendo que el flujo está completamente desarrollado y es
turbulento, la velocidad obtenida es menor que la velocidad máxima que podría obtenerse
en un flujo laminar y en particular de mayor magnitud al promedio, como se evidencia en
los cálculos anteriores.
Figura.2. Región de entrada
Fuente: Mecánica de fluidos Cengel 4ta Edición
Es fundamental mencionar que en los cálculos realizados se asume que la longitud
del tubo es lo suficientemente larga para superar la región de entrada. Por lo tanto, la
velocidad que se está calculando corresponde a la velocidad máxima que posee el fluido
después de haber superado esta región inicial y justo antes de ser perturbada por la aguja.
Velocidad máxima
menor a 2 veces
velocidad
promedio.
Figura.3. Perfil desarrollado en flujos laminares y turbulentos
Fuente: Mecánica de fluidos Cengel 4ta Edición
EJERCICIO 2
Si existe un potencial de velocidades dado por:
𝑢 = 𝑎(𝑥 2 − 𝑦 2 )
𝑣 = −2𝑎𝑥𝑦
𝑤=0
Verificar si el flujo es rotacional y dibujarlo
V = (u, v, w) donde: u, v y w son (x, y, z)
i j k
∂ ∂ ∂
∇∗V =
∂x ∂y ∂x
u v w
u=a ( x 2− y 2 )
v=−2 axy
w=0
i j k
∂ ∂ ∂
∇∗V =
∂x ∂y ∂x
a ( x − y ) −2 axy
2 2
0
Los componentes del rotor se simplifican a:
∇∗V = ( ∂∂ wy − ∂∂ vz , ∂u∂z − ∂∂wx , ∂∂ vx − ∂∂ uy )
(
∇∗V = 0−0 , 0−0 ,
∂
∂x
∂
(−2 axy ) − ( a ( x 2− y 2 ) )
∂y )
El rotacional del campo de velocidades es (0,0,0) lo que indica que el flujo es irrotacional.
Figura.4. Grafica vectorial flujo de rotación
Fuente: Elaboración propia
EJERCICIO 3
¿Cuál de las siguientes expresiones satisface el principio de conservación de la masa para
un fluido incompresible?
a. 𝑢 = −𝑥 ; 𝑣 = 𝑦
b. 𝑢 = 3𝑦 ; 𝑣 = 3𝑥
c. 𝑢 = 4𝑥 ; 𝑣 = −4𝑦
d. 𝑢 = 3𝑥𝑡 ; 𝑣 = 3𝑦𝑡
e. 𝑢 = 𝑥𝑦 + 𝑦 2 𝑡 ; 𝑣 = 𝑥𝑦 + 𝑥 4 𝑡
f. 𝑢 = 4𝑥 2𝑦 3 ; 𝑣 = −2𝑥𝑦 4
Para determinar cuál de las expresiones satisface el principio de conservación de la
masa para un fluido incompresible, se hace uso de la ecuación de continuidad en
coordenadas cartesianas:
∂ ρ ∂( ρu) ∂ ( ρ u) ∂(ρ u)
+ + + =0
∂t ∂x ∂y ∂z
Como es un flujo incompresible, la densidad (ρ) es constante por lo que sale de la
derivada parcial temporal. Entonces, quedaría:
ρ ( +
∂x ∂ y )
∂ ( u) ∂( u )
=0
Se busca entonces satisfacer la ecuación haciendo que la operación de los 2 términos
dentro del paréntesis sea 0:
1) u=−x , v = y
∂ (u) ∂(v)
=−1 =1
∂x ∂y
−1+1=0 | Satisface el principio de conservación de la masa.
2) u=3 y , v=3 x
∂ (u) ∂( v )
=0 =0
∂x ∂y
0+ 0=0 | Satisface el principio de conservación de la masa.
3) u=4 x , v=−4 y
∂ (u) ∂ (v )
=4 =−4
∂x ∂y
4−4=0 | Satisface el principio de conservación de la masa.
4) u=3 xt , v=3 yt
∂ (u) ∂( v )
=3 t =3 t
∂x ∂y
3 t+3 t=6 t | No satisface la ecuación de continuidad a menos que t sea 0, lo cual no es una
condición realista para un flujo de fluido.
5) u=xy+ y 2 t , v=xy + x 4 t
∂ (u) ∂ (v )
=y =x
∂x ∂y
y + x | No satisface el principio de conservación de la masa.
6) u=4 x 2 y 3 , v =−2 x y 4
∂ (u) ∂( v )
=8 x y 3 =−8 x y 3
∂x ∂y
3 3
8 x y −8 x y =0| Satisface el principio de conservación de la masa
EJERCICIO 4
Si el movimiento de un fluido está descrito por un campo de velocidades
𝑉⃗ = 𝑢𝑖̂+ 𝑣𝑗̂+ 𝑤𝑘̂ 𝑢, 𝑣,𝑗: componentes de velocidad en las direcciones 𝑖,̂𝑗̂𝑦 𝑘̂.
Obtener la ecuación que describe el campo de aceleraciones para una partícula de fluido.
Para obtener la ecuación que describe el campo de aceleraciones para una partícula
de fluido, primero debemos considerar la derivada total de la velocidad con respecto al
tiempo.
∂𝑉⃗/∂t = ∂𝑢/∂t 𝑖̂ + ∂𝑣/∂t 𝑗̂ + ∂𝑤/∂t 𝑘̂
Donde:
∂𝑢/∂t: Tasa de cambio de la componente 𝑢 de la velocidad en el tiempo.
∂𝑣/∂t: Tasa de cambio de la componente 𝑣 de la velocidad en el tiempo.
∂𝑤/∂t: Tasa de cambio de la componente 𝑤 de la velocidad en el tiempo.
La derivada convectiva de la velocidad describe cómo cambia la velocidad de una
partícula a medida que se mueve a través de un campo de velocidades espacialmente
variable. Se calcula como:
(𝑉⃗ ⋅ ∇)𝑉⃗ = (𝑢𝑖̂ + 𝑣𝑗̂ + 𝑤𝑘̂) ⋅ (∂𝑖̂/∂𝑥 + ∂𝑗̂/∂𝑦 + ∂𝑘̂/∂𝑧)
Desarrollando el producto escalar y simplificando, obtenemos:
(𝑉⃗ ⋅ ∇)𝑉⃗ = 𝑢(∂𝑢/∂𝑥) + 𝑣(∂𝑢/∂𝑦) + 𝑤(∂𝑢/∂𝑧) + 𝑢(∂𝑣/∂𝑥) + 𝑣(∂𝑣/∂𝑦) + 𝑤(∂𝑣/∂𝑧) +
𝑢(∂𝑤/∂𝑥) + 𝑣(∂𝑤/∂𝑦) + 𝑤(∂𝑤/∂𝑧)
Sustituyendo las derivadas parciales temporales y convectivas en la ecuación de la
derivada total, obtenemos la ecuación que describe el campo de aceleraciones para una
partícula de fluido:
𝐷𝑉⃗/Dt = ∂𝑢/∂t 𝑖̂ + ∂𝑣/∂t 𝑗̂ + ∂𝑤/∂t 𝑘̂
+ 𝑢(∂𝑢/∂𝑥) + 𝑣(∂𝑢/∂𝑦) + 𝑤(∂𝑢/∂𝑧)
+ 𝑢(∂𝑣/∂𝑥) + 𝑣(∂𝑣/∂𝑦) + 𝑤(∂𝑣/∂𝑧)
+ 𝑢(∂𝑤/∂𝑥) + 𝑣(∂𝑤/∂𝑦) + 𝑤(∂𝑤/∂𝑧)
Esta ecuación se conoce como la ecuación de aceleración de Euler y es
fundamental para la descripción del movimiento de los fluidos.
