Problemas de Mecánica de Fluidos
Problemas de Mecánica de Fluidos
Problemas de Mecánica de Fluidos
1. (Puntaje 5) Un sistema de bombeo está diseñado para transferir 225 gal/min de agua a 50
°F desde un depósito con inferior a un tanque elevado como se muestra en la figura. La
tuberia es de acero comercial de D = 4 in y tiene una longitud total de tubería L=368 ft.
Seleccionar del diagrama de curvas de operacion de la bomba adjunto: el diametro del
impulsor y determinar el punto de operación para entregar por lo menos 225 gal/min. Del
mismo diagrama especificar la eficiencia y la potencia de la bomba.
𝑔𝑎𝑙 𝑓𝑡 3
𝑄 = 225 = .5013
𝑚𝑖𝑛 𝑠
4
𝐷 = 4 𝑖𝑛 = 𝑓𝑡
12
𝐿 = 368 𝑓𝑡
De la ecuación de la energía
∝1 𝑣12 𝑃1 ∝2 𝑣22 2
+ + 𝑧1 + ℎ𝑏𝑜𝑚𝑏𝑎 = + + 𝑧2 + ℎ𝑡𝑢𝑟𝑏𝑖𝑛𝑎 + ℎ𝐿𝑡𝑜𝑡𝑎𝑙
2𝑔 𝜌𝑔 2𝑔 𝜌𝑔
Tenemos que
𝑣1 = 𝑣2 = 𝑃1 = 𝑃2 = ℎ𝑡𝑢𝑟𝑏𝑖𝑛𝑎 = 0
Y tenemos que
𝑣2 𝐿
ℎ𝐿𝑡𝑜𝑡𝑎𝑙 =
2𝑔
(𝑓 𝐷 + ∑ 𝑘) ec.2
𝑓𝑡 3
4 (0.5013
𝑄 4𝑄 𝑠 ) 𝑓𝑡
𝑣=𝜋 = = 2 = 5.74
2 𝜋𝐷2 4 𝑠
4𝐷 𝜋 ( 𝑓𝑡)
12
𝑠𝑙𝑢𝑔 𝑓𝑡 4
𝜌𝑣𝐷 1.940 𝑓𝑡 3 (5.74 𝑠 ) (12 𝑓𝑡)
𝑁𝑅𝑒 = = = 135.965 × 103
𝜇 𝑙𝑏 ∗ 𝑠
2.73 × 10−5
𝑓𝑡 2
Tenemos que ∑ 𝑘 = 𝑘𝑣á𝑙𝑣𝑢𝑙𝑎 𝑎𝑏𝑖𝑒𝑟𝑡𝑎 + 𝑘𝑣á𝑙𝑣𝑢𝑙𝑎 𝑐ℎ𝑒𝑐𝑘 + 𝑘𝑐𝑜𝑑𝑜 𝑟𝑒𝑔𝑢𝑙𝑎𝑟 90° + 𝑘𝑐𝑜𝑑𝑜 𝑟𝑒𝑔𝑢𝑙𝑎𝑟 90° + 𝑘𝑠𝑎𝑙𝑖𝑑𝑎 +
𝑘𝑒𝑛𝑡𝑟𝑎𝑑𝑎 = .11 + 2 + 2 ∗ .23 + 1 + .5 = 4.07
Sustituyendo en la ecuación 2
𝑓𝑡 2
(5.74 ) 368𝑓𝑡
𝑠2
ℎ𝐿𝑇𝑜𝑡𝑎𝑙 = ((0.0175) ( ) + 4.07) = 13.03 𝑓𝑡
𝑓𝑡 4 𝑖𝑛
2 (32.2 2 )
𝑠 12 𝑖𝑛 𝑓𝑡
Con esto y con ayuda del gráfico, observamos que con 225 gal/min y una carga de 183.03 ft,
la potencia sería de 20 HP y con una eficiencia del 54%, además de que también se observa
que el diámetro del impulsor es de 8 in
Punto de operación
𝑣2 368 𝑓𝑡
ℎ𝑏𝑜𝑚𝑏𝑎 = 170 𝑓𝑡 + 2𝑔 (𝑓 4 + 4.07)
𝑓𝑡
12
𝑄 4𝑄
Donde 𝑣 = 𝜋 = 𝜋𝐷2
𝐷2
4
Con ayuda de Matlab, tomamos diferentes valores de Q desde 0 hasta 250 ft, de modo que
obtenemos la siguiente curva
Q=(0:5:250);
%Conversion gal/min a ft3/s
Qf = Q*0.00222801;
%Cálculo de la velocidad:
V = Qf / ((pi/4)*(16/144));
%Cálculo de número de Reynolds
NRE = (1.94*V*(4/12))/(2.73e-5);
%Cálculo de factor de fricción.
f=0.25./(log10((0.00045/3.7)+(5.74./(NRE.^0.9)))).^2;
%Ecuación característica de la bomba
hbomba = 170+(2.039)*Qf.^2.*(f*1104+4.07);
hold on;
plot(Q,hbomba,'linewidth',2);
grid minor;
xlim([0 400]);
ylim([50 475]);
La superponemos con la gráfica de la bomba y colocamos un punto de operación a los 225
gal/min
Diámetro impulsor = 8 in
Potencia= 20 HP
Carga= 183.03 ft
Datos conocidos:
𝐿 =7𝑚
𝐾𝑓𝑖𝑙𝑡𝑟𝑜 = 8
𝜖 = 0.002 𝑚𝑚
𝐿 𝑚3
𝑄 = 150 = 2.5 × 10−3
𝑚𝑖𝑛 𝑠
De la ecuación de energía
∝1 𝑣12 𝑃1 ∝2 𝑣22 2
+ + 𝑧1 + ℎ𝑏𝑜𝑚𝑏𝑎 = + + 𝑧2 + ℎ𝑡𝑢𝑟𝑏𝑖𝑛𝑎 + ℎ𝑡𝑜𝑡𝑎𝑙
2𝑔 𝜌𝑔 2𝑔 𝜌𝑔
𝑣2 𝐿
ℎ𝑡𝑜𝑡𝑎𝑙 = (𝑓 + ∑ 𝑘)
2𝑔 𝐷
donde ∑ 𝑘 = 𝑘𝑠𝑎𝑙𝑖𝑑𝑎 + 𝑘𝑓𝑖𝑙𝑡𝑟𝑜 + 𝑘𝑒𝑛𝑡𝑟𝑎𝑑𝑎 = 0.68 + 8 + 1 = 9.68
𝑣2 7𝑚
5𝑚= 𝑚 (𝑓 𝐷 + 9.68)
2(9.81 2 )
𝑠
𝑄 4𝑄
la velocidad la obtenemos de la ecuación de continuidad 𝑣 = 𝜋 = 𝜋𝐷2
𝐷2
4
4𝑄 2
( 2) 7𝑚
5 𝑚 = 𝜋𝐷 𝑚 (𝑓 + 9.68)
2(9.81 2 ) 𝐷
𝑠
obtenemos
10𝑔𝐷 5 𝜋2 9.68
𝑓= − 𝐷 …ec1
112𝑄 2 7
Con ayuda de la rugosidad podemos obtener la siguiente ecuación, dado que la rugosidad
para el acero inoxidable es 0.002 × 10−3
𝜀 0.002×10−3
= …ec2
𝐷 𝐷
Otra de las ecuaciones de las que podemos auxiliarnos es la del NRe con 𝜌 =
786 (𝑝𝑎𝑟𝑎 𝑎𝑙𝑐𝑜ℎ𝑜𝑙 𝑖𝑠𝑜𝑝𝑟𝑜𝑝í𝑙𝑖𝑐𝑜 𝑎 25 𝑔𝑟𝑎𝑑𝑜𝑠)𝑦 𝜇 = 1.945 × 10−3
786𝑘𝑔4𝑄𝐷
𝜌𝑣𝐷 𝑚3 𝜋𝐷2 786(4𝑄)
𝑁𝑅𝑒 = = = 𝜋(1.945×10−3 )𝐷 …ec3
𝜇 𝜇
Con estas ecuaciones podemos suponer un valor de 𝐷 con el que obtendremos un valor de 𝑓
que usaremos en las ecuaciones 2 y 3 para de esta forma, obtener el coeficiente de rugosidad
relativa y el NRe, para así usar el diagrama de Moody para encontrar un 𝑓 que coincida con el
calculado en la ecuación 1, o para ser más exactos, usamos la ecuación de Swameey- Jain
0.25
𝑓= 2
𝜖
𝐷 + 5.74 )]
[log ( 3.7
𝑁𝑅𝑒 0.9
Iniciamos con valores inferiores a .05
format long
syms D
%ecuación para obtener factor de fricción
e=((10*9.81*(pi^2)*(D^5))/(112*(0.0025^2)))-(9.68*D/7);
a= subs (e,0.0348)
f2= 0.000002/D;
b= subs (f2,0.0348)
f3= (786*4*.0025)/(pi*0.00194*D);
c= subs(f3,0.0348)
%Ecuación de Swamee-Jain
x=0.25/(log10((b/3.7)+(5.74/(c^0.9))))^2
𝝐
D propuesto Valor 𝒇 Valor 𝑫 Valor Nre Valor
obtenido con la obtenido en obtenido en ec. 𝒇 obtenido con
ecuación 1 ecuación 2 3 ecuación de
Swameey- Jain
0.04 m 0.086 5 × 10−5 3.25 × 104 0.023
0.045 m 0.19 4.44 × 10−5 2.89 × 104 0.023
0.048 m 0.28 4.16 × 10−5 2.71 × 104 0.024
0.035 m 0.0242 5.71 × 10−5 3.68 × 104 0.0224
0.0349 m 0.0233 5.73 × 10−5 3.69 × 104 0.0224
0.03495 m 0.023 5.72 × 10−5 3.68 × 104 0.0224
0.0348 m 0.02247 5.74 × 10−5 3.705 × 104 0.02243
Como podemos ver en la tabla, con D=0.0348 m obtenemos un valor 𝑓 que satisface las
ecuaciones necesarias para que ese valor sea considerado como el correcto.