Fluid HW involving MATLAB cource

example_program.pdf

clc;clear V1 = input('input guess for velocity in 100mm piping '); Error = 1; while Error >0.001 rr1 = 0.15/100; rr2 = 0.15/50; V2 = -4*V1+15.279;

%calculate reynolds number for each pipeline

Re1 = V1*0.1/(1e-06); Re2 = V2*0.05/(1e-06);

%calculate friction factor in each pipeline

f1 = @(ff)1/ff^0.5+2*log10(rr1/3.7+2.51/((Re1)*ff^0.5)); Friction1 = fzero(f1,[0.008 0.1]);

f2 = @(ff)1/ff^0.5+2*log10(rr2/3.7+2.51/((Re2)*ff^0.5)); Friction2 = fzero(f2,[0.008 0.1]);

%calculate head loss in each pipeline

hl1 = Friction1*(3/0.1)*V1^2/2 + 2*Friction1*(20)*V1^2/2+ 1.6*V1^2/2 +

Friction1*340*V1^2/2;

hl2 = Friction2*(7/0.05)*V2^2/2 + 2*Friction2*(60)*V2^2/2 +

Friction2*(8)*V2^2/2 + 2*Friction2*(30)*V2^2/2;

Error = abs(hl1 - hl2); V1 = V1 + 0.00001; end