pro14

Problema 14

Código de programación Matlab

% PROBLEMA 14

clc
clear
format compact

% Una tubería vieja de 2m de diámetro tiene una rugosidad de 30mm. Si se le
% aplica un revestimiento de 12mm de espesor, la rugosidad se reducirá a 1 mm.
% ¿Qué tanto se ahorra en el costo de bombeo por cada año y por cada km de
% tubería para un gasto de 6 m^3/s de agua a 20°C?. La bomba y los motores
% tiene eficiencia del 80% y el costo de la energía es de 1 centavo por cada
% kilowatt-hora.

dz=input(‘Valor de la diferencia de alturas= ‘);
%dz=0
u1=input(‘Valor de la velocidad del fluido= ‘);
%u1=0
Q=input(‘Valor del gasto [m^3/s]’);
%Q=6
d=input(‘Valor de la densidad del fluido []= ‘);
%d=998.2
v=input(‘Valor de la viscosidad= ‘);
%v=9.8e-4

disp(‘Procedimiento para calcular el costo antes del revestimiento’)

D=input(‘Valor del diámetro de la tuberías en [m]= ‘);
%D=2
e=input(‘Valor de E= ‘);
%e=.03
s=pi*(D/2)^2
U2=Q/s
Nre=(D*U2*d)/v
ed=e/D
f=input(‘Valor del factor de fricción= ‘);
%f=.044
gc=9.81;
L=input(‘Valor del tramo recto= ‘);
%L=1000
disp(‘El valor de Wf esta dado en [kgfm/kg]’)
Wf=((U2*U2)/(2*gc))+(f*U2*U2*L)/(2*gc*D)
disp(‘El valor de wf esta dado en [lbfm/lb]’)
wf=Wf/(.3048)
w=Q*d
w1=w/(.4536)
Ph=wf*w1
ph1=(Ph*745.7)/550
E=input(‘Valor de la eficiencia global = ‘);
%E=.8
Pr=ph1/E
disp(‘kw es el valor del costo antes del revestimiento ‘)
kw1=(Pr*3600)/3.6e6
disp(‘Procedimiento para calcular el costo después del revestimiento’)
r=input(‘Valor del revestimiento [m]= ‘);
%r=.012
D2=D-r
e2=input(‘Valor de E= ‘);
%e2=1e-3
s2=pi*(D/2)^2
U3=Q/s2
Nre=(D2*U3*d)/v
ed=e2/D2
f2=input(‘Valor del factor de fricción= ‘);
%f2=.0168
gc=9.81;
disp(‘El valor de Wf esta dado en [kgfm/kg]’)
Wf2=U3/(2*gc)+(f2*U3*U3*L)/(2*gc*D2)
disp(‘El valor de wf esta dado en [lbfm/lb]’)
wf2=Wf2/(.3048)
w2=Q*d;
w3=w2/(.4536)
Ph2=wf2*w3;
ph3=(Ph2*745.7)/550
pr2=ph3/E
kw2=pr2*3600
disp(‘kw3 es el valor del costo después del revestimiento ‘)
kw3=kw2/3.6e6
disp(‘Por lo tanto el ahorro en el costo del bombeo por kilómetro es’)
ahorro=kw1 – kw3
disp(‘Por lo tanto el ahorro en el costo del bombeo por kilómetro por año=’)
a=ahorro*24*365