Problema 21

Código de programación Matlab

% PROBLEMA 21

clc
clear
format compact

% Agua a 60°F fluye a razón de 1.2 ft^3/s a través de una repentina
% expansión en una tubería horizontal; el diámetro cambia abruptamente de 4
% in a 8 in. Calcular la presión ocasionada por la expansión.

Q=input(‘Valor del flujo volumétrico [ft^3/s] ‘);
%Q=1.2
D1=input(‘Valor del diámetro interno de la tubería 1 [ft] ‘);
%D1=.3355
D2=input(‘Valor del diámetro interno de la tubería 2 [ft] ‘);
%D2=.66513
S1=input(‘Valor del área de la tubería 1 [ft^2] ‘);
%S1=0.0884
S2=input(‘Valor del área de la tubería 2 [ft^2] ‘);
%S2=0.3474
f1=input(‘Valor del factor de fricción 1 ‘);
%f1=.017
f2=input(‘Valor del factor de fricción 2 ‘);
%f2=.014
d=input(‘Valor de la densidad [lb/ft^3] ‘);
%d=62.371
L=input(‘Valor del tramo recto de la tubería [ft] ‘);
%L=1
gc=32.2;

disp(‘Los valores de K son adimensionales’)
kexp=[1-(D1^2/D2^2)]^2
ktr1=f1*(L/D1)
ktr2=f2*(L/D2)
kta=ktr1+kexp
ktb=ktr2+kexp
disp(‘Las unidades de la velocidad son ft/s’)
ua=Q/S1
ub=Q/S2
disp(‘Las unidades de los Hfs son [lbf ft/lb]’)
hfsa=(kta*ua*ua)/(2*gc)
hfsb=(ktb*ub*ub)/(2*gc)
disp(‘La diferencia de presión está dada en [lbf/ft^2]’)
dp=(hfsa+hfsb)*d