blob: 1d3fcb915efaf9a407e933b4ffab7c2ce70612bc (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
|
clc;
clear;
dia=120;//mm
p=1.0;//kPa
//using energy equation
//Q=A2*v2=A2*((p1-p2)/(d*(1+Kl)/2)); d =density, Kl= loss coefficient
Kl1=0.05;;
Kl2=0.5;
//for rounded entrance cyliindrical vent
Q1=(%pi*((dia/1000)^2)/4)*(p*1000*2/(1.23*(1+Kl1)))^0.5;
//for cylindrical vent
Q2=(%pi*((dia/1000)^2)/4)*(p*1000*2/(1.23*(1+Kl2)))^0.5;
disp("m^3/sec",Q1,"The volume fowrate associated with the rounded entrance cylindrical vent configuration =")
disp("m^3/sec",Q2,"The volume fowrate associated with the cylindrical vent configuration =")
KLoss=0:0.01:0.5;
count=1;
for i=0:0.01:0.5
flow(count)=(%pi*((dia/1000)^2)/4)*(p*1000*2/(1.23*(1+i)))^0.5;
count=count+1;
end
plot2d(KLoss,flow,rect=[0,0,0.5,0.5])
xtitle("Q vs KL","KL","Q, (m^3)/sec")
|
