UNIVERSIDAD NACIONAL DEL CENTRO DEL PERÚ
FACULTAD DE INGENIERÍA QUÍMICA
SEGUNDA EVALUACIÓN PARCIAL ESCRITA DE PROCESOS DE TRANSFERENCIA DE FLUIDOS
APELLIDOS Y NOMBRES: CUESTAS CHAGUA ANA CRISTINA CICLO: VI SECCIÓN: A EAP: IQ
INDICACIONES GENERALES:
- El examen tendrá una duración de horas
- Se podrá emplear, tablas, gráficos, apuntes de clases, separatas, etc.
- El orden, la limpieza, la claridad en las explicaciones, así como los resultados numéricos se tomarán en cuenta en la
calificación.
PROBLEMA 1.- Al sumergir un pequeño tubo en un recipiente que contiene un líquido, se forma
un menisco en la superficie libre debido a la tensión superficial. Los experimentos realizados
señalan que la magnitud de este efecto capilar (∆𝒉) es función del diámetro (D), del peso
específico del líquido (γ) y la tensión superficial (𝝈). Determinar el número de parámetros
repetitivos y hallar los parámetros o números adimensionales correspondientes.
PROBLEMA 2.- Por medio de un sifón ABC se desplaza líquido desde un depósito superior a
otro inferior, calcular el valor máximo de z hasta el punto más alto B, con la condición de que
el líquido que circula permanezca unido.
VOLUMEN DE CONTROL
PROBLEMA 3.- Un venturímetro es situado en una tubería vertical que lleva aceite de dr =
0,90. Los diámetros de entrada y de contracción son de 6 cm y 3 cm respectivamente y la
contracción esta 12 cm por encima de la entrada.
Determinar para una velocidad de 3,0 m/s en la tubería de 6 cm.
a) La diferencia de nivel de mercurio (dr = 13,6) en un manómetro vertical diferencial
conectado entre la entrada y la salida.
b) El caudal que pasa por la tubería
c) El dato de la distancia entre la toma de presión y otra ¿es necesario? y si no es
necesario en qué casos si se requiere.
d) Según el análisis de la parte (c) ¿Qué podemos afirmar, que es independiente de dicha
distancia?
PROBLEMA 4. - Encuéntrese la descarga en el sistema de la figura para las bombas A y B
cada una independientemente, en paralelo o en serie; 1600 m de tubería de hierro forjado
nuevo, D= 300 mm y H = 26.5 m y determínese la potencia necesaria para operar las bombas,
recomiende cuál de las propuestas seria la adecuada para el servicio de bombeo.
BOMBA A BOMBA B
H
Hp (m) Qp ( l/s) e% Hp (m) Qp ( l/s) e %
21,3 0 0 24,4 0 0
18,3 56,6 59 21,3 74 54 P
16,8 72,5 70 18,3 112 70
15,2 85,8 76 15,2 140 80
13,7 97,7 78 12,2 161 73
12,2 108,0 76,3 9,1 174 60
10,7 116,0 72 6,1 177 40
9,1 127,0 65
7,6 130,0 56,5
6,1 134,0 42
Resolución:
El volumen de control es el líquido que se encuentra en el tanque 1, 2 y el líquido que fluye
internamente por la tubería.
H=26.5m
PRIMERA PARTE
Balance de energía total
𝜕𝑄 𝜕𝑊𝑠 𝑃 𝜕 𝜕𝑊𝜇 (1)
− = ∬ (𝑒 + ) ∙ 𝜌(𝑣⃗ ∙ 𝑛⃗⃗)𝑑𝐴 + ∭ 𝑒 ∙ 𝜌 ∙ 𝑑𝑉 +
𝑑𝑡 𝑑𝑡 𝐶.𝑆 𝜌 𝜕𝑡 𝐶.𝑉 𝑑𝑡
De acuerdo a las condiciones mencionadas en el enunciado se tiene:
𝜕𝑄 𝜕𝑊𝑠 𝑃 𝜕 𝜕𝑊𝜇
− = ∬ (𝑒 + ) ∙ 𝜌(𝑣⃗ ∙ 𝑛⃗⃗)𝑑𝐴, 𝑑𝑜𝑛𝑑𝑒, ∭ 𝑒 ∙ 𝜌 ∙ 𝑑𝑉 = 0, =0
𝑑𝑡 𝑑𝑡 𝐶.𝑆 𝜌 𝜕𝑡 𝐶.𝑉 𝑑𝑡
Si se sabe que 𝑒 representa los tipos de energía (cinética, potencial e interna) y que la masa
únicamente entra al volumen de control en la sección (1) y sale en (2), la integral de
superficie se transforma en:
𝑃 𝑣22 𝑃2 𝑣12 𝑃1
∬ (𝑒 + ) ∙ 𝜌(𝑣⃗ ∙ 𝑛⃗⃗)𝑑𝐴 = [ + 𝑔𝑦2 + 𝑢2 + ] (𝜌2 𝑣2 𝐴2 ) − [ + 𝑔𝑦1 + 𝑢1 + ] (𝜌1 𝑣1 𝐴1 )
𝐶.𝑆 𝜌 2 𝜌2 2 𝜌1
Balance de materia, para el equilibrio de masa, se tiene:
𝜌1 𝑣1 𝐴1 = 𝜌2 𝑣2 𝐴2 = 𝑚
Sustituyendo, todos estos valores en la ecuación (1), la expresión que corresponde al
balance de energía, en este problema, será.
𝛿𝑄 𝛿𝑊𝑠 𝑣22 𝑃2 𝑣12 𝑃1
− = [ + 𝑔𝑦2 + 𝑢2 + ] (𝜌2 𝑣2 𝐴2 ) − [ + 𝑔𝑦1 + 𝑢1 + ] (𝜌1 𝑣1 𝐴1 ) (2)
𝑑𝑡 𝑑𝑡 2 𝜌2 2 𝜌1
Si se divide la ecuación (2) entre la rapidez de flujo de masa, se obtiene
𝑞 − 𝑊𝑠 𝑣22 𝑃2 𝑣12 𝑃1
= [ + 𝑔𝑦2 + 𝑢2 + ] − [ + 𝑔𝑦1 + 𝑢1 + ] (3)
𝑚 2 𝜌2 2 𝜌1
𝑞 𝑊𝑠 𝑣22 − 𝑣12 𝑃2 𝑃1 𝑢2 − 𝑢1 𝑞 𝑊𝑠
− = 𝑦2 − 𝑦1 + + − + , 𝑑𝑜𝑛𝑑𝑒, = 0, = 𝐻𝑚
𝑚𝑔 𝑚𝑔 2𝑔 𝜌2 𝑔 𝜌1 𝑔 𝑔 𝑚𝑔 𝑚𝑔
Así, tenemos:
𝑣12 𝑃1 𝑣22 𝑃2 𝐿 𝑣2
𝑦1 + + + (−𝐻𝑚 ) = 𝑦2 + + + ℎ𝐿 , 𝑑𝑜𝑛𝑑𝑒, 𝑢2 − 𝑢1 = ℎ𝐿 𝑔, ℎ𝐿 = 𝑓
2𝑔 𝛾 2𝑔 𝛾 𝐷 2𝑔
Aplicando a las condiciones del problema se tiene:
(o aplicando la ecuación de energía entre los puntos (1) y (2))
(4)
𝑙 𝑣2
0 + 0 + 0 + (−𝐻𝑚 ) = 𝐻 + 0 + 0 + 𝑓
𝐷 2𝑔
𝑙 𝑣2
−𝐻𝑚 = 𝐻 + 𝑓
𝐷 2𝑔
Si se sabe que:
𝜋𝐷 2 4𝑄 8𝑙
𝑄= 𝑣, 𝑣 = , ℎ𝐿 = 𝑓 𝑄2
4 𝜋𝐷 2 𝜋 2 𝐷5 𝑔
Reemplazando en la ecuación (4), se obtiene
8𝑙 8𝑙
−𝐻𝑚 = 𝐻 + 𝑓 𝑄 2 = 𝐻 + 𝑓𝑐1 𝑄2 ; 𝑎𝑠𝑖 𝑠𝑒 𝑡𝑖𝑒𝑛𝑒, 𝑐1 =
𝜋 2 𝐷5 𝑔 𝜋 2 𝐷5 𝑔
−𝐻𝑚 = 𝐻 + 𝑓𝑐1 𝑄2
8𝑙
𝑊𝑝 = −𝐻𝑚 = 𝐻 + 𝑓 𝑄2
𝜋 2 𝐷5 𝑔
Para determinar la curva de operación se tiene de manera general
8𝑙 𝜀
𝑊𝑝 = 𝐻 + 𝑓 𝑄2 , = 𝑐𝑡𝑒
𝜋 2 𝐷5 𝑔 𝐷
SEGUNDA PARTE
DATOS:
Densidad(agua) 998 kg/m3
Viscosidad di. (agua) 0.001 kg/(ms)
Largo(L) 1600 m
Diámetro(D) 0.3 m
Altura(H) 26.5 m
Gravedad(g) 9.81 m/s2
Rugosidad Absoluta 0.000046 m
Hallando la curva de operación en:
BOMBA A y B:
PROGRAMA:
function FLUIDOS
clear all;
clc;
% El presente programa nos permite calcular el wp de una tubería
% con dos bombas en forma individual para determinados caudales
% graficando así una curva de operación y obteniendo un Q0 y H0
% en la gráfica Q vs H y Q vs e del fabricante
disp('==========================================================')
disp(' BOMBA A Y B EN FORMA INDIVIDUAL ')
disp('==========================================================')
% DATOS DEL PROBLEMA
H=26.5; % m
l=1600; % m
D=0.3; % m
vis=1.007*10^-6; % m^2/s
g=9.81; %m/s^2
E=0.0046; % Hierro forjado cm
e=E/(D*100);
fprintf('Para e/D = %2.6f = constante\n',e)
% DATOS DEL FABICANTE
%Bomba A
QA=[0 0.0566 0.0725 0.0858 0.0977 0.108 0.116 0.127 0.130 0.134]; %
m^3/s
MA=[21.3 18.3 16.8 15.2 13.7 12.2 10.7 9.1 7.6 6.1]; % m
EA=[0 59 70 76 78 76.3 72 65 56.5 42];
%Bomba B
QB=[0 0.074 0.112 0.14 0.161 0.174 0.177];
MB=[24.4 21.3 18.3 15.2 12.2 9.1 6.1];
EB=[0 54 70 80 73 60 40];
% PARTE OPERATIVA
% Valores de los caudales para la Bomba A y B
QA1=[0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1];
QB1=[0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1 0.12 0.14
0.15];
k1=length(QA1);
k2=length(QB1);
% Para la Bomba A
disp('==========================================================')
disp(' BOMBA A ')
disp('==========================================================')
disp(' Q(l/s) | Re | f | HL(m) | wp(m)')
disp('==========================================================')
% Determinación de Re, factor fricción, hl y Wp
for i=1:k1
Rea(i)=(QA1(i)*4)/(vis*pi*D);
fa(i)=1/(-2*log10(E/(3.7065*D*100)-
(5.0452/Rea(i))*log10((1/2.8257)*(E/(D*100))^1.1098+(5.8506/Rea(i)^
0.898))))^2; %Ecuación de Chen. FF
HLA(i)=fa(i)*(8*l/(pi^2*D^5*g)*QA1(i)^2);
wpA(i)=H+HLA(i);
fprintf(' %4.0f | %11.4f | %2.4f | %7.4f
|%2.2f\n',QA1(i)*1000,Rea(i),fa(i),HLA(i),wpA(i))
end
% Para la Bomba B
disp('==========================================================')
disp('==========================================================')
disp(' BOMBA B ')
disp('==========================================================')
disp(' Q(l/s) | Re | f | HL(m) | wp(m)')
disp('==========================================================')
% Determinación de Re, factor fricción, hl y Wp
for j=1:k2
Reb(j)=(QB1(j)*4)/(vis*pi*D);
fb(j)=1/(-2*log10(E/(3.7065*D*100)-
(5.0452/Reb(j))*log10((1/2.8257)*(E/(D*100))^1.1098+(5.8506/Reb(j)^
0.898))))^2; %Ecuación de Chen. FF
HLB(j)=fb(j)*(8*l/(pi^2*D^5*g)*QB1(j)^2);
wpB(j)=H+HLB(j);
fprintf(' %4.0f | %11.4f | %2.4f | %7.4f
|%2.2f\n',QB1(j)*1000,Reb(j),fb(j),HLB(j),wpB(j))
end
%GRÁFICAS
% Datos del fabricante (Q - m y Q - e)
subplot(2,2,1),plotyy(QA*1000,MA,QA*1000,EA),grid
title('BOMBA A')
xlabel('Q(L/S)'),ylabel('H(m)')
subplot(2,2,2),plotyy(QB*1000,MB,QB*1000,EB),grid
title('BOMBA B')
xlabel('Q(L/S)'),ylabel('H(m)')
% Datos del fabricante - Curva de Operación
subplot(2,2,3),plot(QA*1000,MA,'-og',QA1*1000,wpA,'-^r'),grid
title('BOMBA A - CURVA DE OPERACION')
legend('Datos del fabricante','Curva de Operación')
xlabel('Q(L/S)'),ylabel('H(m)')
subplot(2,2,4),plot(QB*1000,MB,'-og',QB1*1000,wpB,'-^r'),grid
title('BOMBA B - CURVA DE OPERACION')
legend('Datos del fabricante','Curva de Operación')
xlabel('Q(L/S)'),ylabel('H(m)')
end
Haciendo correr el programa:
==========================================================
BOMBA A Y B EN FORMA INDIVIDUAL
==========================================================
Para e/D = 0.000153 = constante
==========================================================
BOMBA A
==========================================================
Q(l/s) | Re | f | HL(m) | wp(m)
==========================================================
0| 0.0000 | 0.0000 | 0.0000 |26.50
10 | 42146.2941 | 0.0222 | 0.1210 |26.62
20 | 84292.5882 | 0.0194 | 0.4225 |26.92
30 | 126438.8823 | 0.0181 | 0.8855 |27.39
40 | 168585.1764 | 0.0173 | 1.5033 |28.00
50 | 210731.4705 | 0.0167 | 2.2724 |28.77
60 | 252877.7646 | 0.0163 | 3.1905 |29.69
70 | 295024.0587 | 0.0160 | 4.2560 |30.76
80 | 337170.3528 | 0.0157 | 5.4678 |31.97
90 | 379316.6469 | 0.0155 | 6.8252 |33.33
100 | 421462.9410 | 0.0153 | 8.3275 |34.83
==========================================================
BOMBA B
==========================================================
Q(l/s) | Re | f | HL(m) | wp(m)
==========================================================
0| 0.0000 | 0.0000 | 0.0000 |26.50
10 | 42146.2941 | 0.0222 | 0.1210 |26.62
20 | 84292.5882 | 0.0194 | 0.4225 |26.92
30 | 126438.8823 | 0.0181 | 0.8855 |27.39
40 | 168585.1764 | 0.0173 | 1.5033 |28.00
50 | 210731.4705 | 0.0167 | 2.2724 |28.77
60 | 252877.7646 | 0.0163 | 3.1905 |29.69
70 | 295024.0587 | 0.0160 | 4.2560 |30.76
80 | 337170.3528 | 0.0157 | 5.4678 |31.97
90 | 379316.6469 | 0.0155 | 6.8252 |33.33
100 | 421462.9410 | 0.0153 | 8.3275 |34.83
120 | 505755.5292 | 0.0150 | 11.7648 |38.26
140 | 590048.1174 | 0.0148 | 15.7770 |42.28
150 | 632194.4115 | 0.0147 | 17.9981 |44.50
ACOPLANDO LA BOMBA A Y B EN SERIE:
PROGRAMA:
function FLUIDOS_Bomba_serie
clear all;
clc;
% El presente programa nos permite calcular el wp de una tubería
% con dos bombas en serie para determinados caudales
% graficando así una curva de operación y obteniendo un Q0 y H0
% en la gráfica Q vs H y Q vs e del fabricante
disp('==========================================================')
disp(' BOMBA A Y B EN SERIE ')
disp('==========================================================')
% DATOS DEL PROBLEMA
h=26.5; % m
l=1600; % m
D=0.3; % m
vis=1.007*10^-6; % m^2/s
g=9.8; %m/s^2
E=0.0046; % Hierro forjado cm
e=E/(D*100);
% PARTE OPERATIVA
% Valores de los caudales para la Bomba A y B
Q=[0 0.04 0.06 0.08 0.1 0.12];
disp(' Q(l/s)| HA(m) | HB(m) | H(m) ')
disp('=====================================')
% Determinación de la altura a dichos caudales de cada Bomba por la
gráfica
HA=[21.3 19.65 18 15.85 13.35 9.9];
HB=[24.4 23.2 22.15 20.9 19.4 17.5];
H=HA+HB;
n=length(Q);
for i=1:n
fprintf(' %4.0f | %5.2f | %4.2f | %4.2f
\n',Q(i)*1000,HA(i),HB(i),H(i))
end
disp('==========================================================')
disp(' BOMBAS EN SERIE ')
disp('==========================================================')
disp(' Q(l/s) | Re | f | HL(m) | wp(m)')
disp('==========================================================')
% Determinación de Re, factor fricción, hl y Wp
for i=1:n
Re(i)=(Q(i)*4)/(vis*pi*D);
f(i)=1/(-2*log10(E/(3.7065*D*100)-
(5.0452/Re(i))*log10((1/2.8257)*(E/(D*100))^1.1098+(5.8506/Re(i)^0.
898))))^2; %Ecuación de Chen. FF
HL(i)=f(i)*(8*l/(pi^2*D^5*g)*Q(i)^2);
wp(i)=h+HL(i);
fprintf(' %4.0f | %11.4f | %2.4f | %7.4f
|%2.2f\n',Q(i)*1000,Re(i),f(i),HL(i),wp(i))
end
% GRÁFICA
% Datos en serie vs Curva de Operación
plot(Q*1000,H,'-og',Q*1000,wp,'-^r'),grid
title('BOMBAS EN SERIE - CURVA DE OPERACION')
legend('Datos en serie','Curva de Operación')
xlabel('Q(L/S)'),ylabel('H(m)')
end
Haciendo correr el programa:
==========================================================
BOMBA A Y B EN SERIE
==========================================================
Q(l/s)| HA(m) | HB(m) | H(m)
=====================================
0 | 21.30 | 24.40 | 45.70
40 | 19.65 | 23.20 | 42.85
60 | 18.00 | 22.15 | 40.15
80 | 15.85 | 20.90 | 36.75
100 | 13.35 | 19.40 | 32.75
120 | 9.90 | 17.50 | 27.40
==========================================================
BOMBAS EN SERIE
==========================================================
Q(l/s) | Re | f | HL(m) | wp(m)
==========================================================
0| 0.0000 | 0.0000 | 0.0000 |26.50
40 | 168585.1764 | 0.0173 | 1.5049 |28.00
60 | 252877.7646 | 0.0163 | 3.1937 |29.69
80 | 337170.3528 | 0.0157 | 5.4734 |31.97
100 | 421462.9410 | 0.0153 | 8.3360 |34.84
120 | 505755.5292 | 0.0150 | 11.7768 |38.28
Intersección (117,5,27.6)
ACOPLANDO LA BOMBA A Y B EN PARALELO:
PROGRAMA:
function FLUIDOS_Bomba_paralelo
clear all;
clc;
% El presente programa nos permite calcular el wp de una tubería
% con dos bombas en paralwlo para determinados caudales
% graficando así una curva de operación y obteniendo un Q0 y H0
% en la gráfica Q vs H y Q vs e del fabricante
disp('==========================================================')
disp(' BOMBA A Y B EN PARALELO ')
disp('==========================================================')
% DATOS DEL PROBLEMA
h=25.5; % m
l=1500; % m
D=0.25; % m
vis=1.007*10^-6; % m^2/s
g=9.8; %m/s^2
E=0.0046; % Hierro forjado cm
e=E/(D*100);
% PARTE OPERATIVA
% Valores de las alturas para la Bomba A y B
H=[8 10 12 14 16 18 20 21 21.3];
disp(' H(m) | QA | QB | Q(L/s) ')
disp('=======================================')
% Determinación de la altura a dichos caudales de cada Bomba por la
gráfica
QA=[127.5 119.5 109 95 78.5 59.5 34.5 13 0];
QB=[175.5 171 162.5 150 133.5 115 92.5 77.5 74];
Q=(QA+QB)/1000;
n=length(H);
for i=1:n
fprintf(' %4.0f | %6.2f | %6.2f | %4.2f
\n',H(i),QA(i),QB(i),Q(i)*1000)
end
disp('==========================================================')
disp(' BOMBAS EN PARALELO ')
disp('==========================================================')
disp(' Q(l/s) | Re | f | HL(m) | wp(m)')
disp('==========================================================')
% Determinación de Re, factor fricción, hl y Wp
for i=1:n
Re(i)=(Q(i)*4)/(vis*pi*D);
f(i)=1/(-2*log10(E/(3.7065*D*100)-
(5.0452/Re(i))*log10((1/2.8257)*(E/(D*100))^1.1098+(5.8506/Re(i)^0.
898))))^2; %Ecuación de Chen. FF
HL(i)=f(i)*(8*l/(pi^2*D^5*g)*Q(i)^2);
wp(i)=h+HL(i);
fprintf(' %4.0f | %12.4f | %2.4f | %7.4f
|%2.2f\n',Q(i)*1000,Re(i),f(i),HL(i),wp(i))
end
% GRÁFICA
% Datos en serie vs Curva de Operación
plot(Q*1000,H,'-og',Q*1000,wp,'-^r'),grid
title('BOMBAS EN PARALELO - CURVA DE OPERACION')
legend('Datos en paralelo','Curva de Operación')
xlabel('Q(L/S)'),ylabel('H(m)')
end
Haciendo correr el programa:
==========================================================
BOMBA A Y B EN PARALELO
==========================================================
H(m) | QA | QB | Q(L/s)
=======================================
8 | 127.50 | 175.50 | 303.00
10 | 119.50 | 171.00 | 290.50
12 | 109.00 | 162.50 | 271.50
14 | 95.00 | 150.00 | 245.00
16 | 78.50 | 133.50 | 212.00
18 | 59.50 | 115.00 | 174.50
20 | 34.50 | 92.50 | 127.00
21 | 13.00 | 77.50 | 90.50
21 | 0.00 | 74.00 | 74.00
==========================================================
BOMBAS EN PARALELO
==========================================================
Q(l/s) | Re | f | HL(m) | wp(m)
==========================================================
303 | 1532439.2534 | 0.0142 | 165.8254 |191.33
291 | 1469219.8123 | 0.0142 | 152.7195 |178.22
272 | 1373126.2618 | 0.0143 | 133.8261 |159.33
245 | 1239101.0465 | 0.0144 | 109.5452 |135.05
212 | 1072201.7219 | 0.0145 | 82.6817 |108.18
175 | 882543.3984 | 0.0147 | 56.6998 |82.20
127 | 642309.5221 | 0.0150 | 30.7403 |56.24
91 | 457708.7539 | 0.0155 | 16.0859 |41.59
74 | 374259.0916 | 0.0158 | 10.9781 |36.48
De la gráfica no hay ninguna intersección de los datos de las bombas en
paralelo y los datos de la curva de operación.
RESPUESTA:
Lo que se recomienda para la elección de las bombas, es la de mayor de
eficiencia, por lo que se ve en el caso del acoplamiento de las bombas en
serie Q=117.5 L/s y H=27.6 m