1 files changed, 24 insertions, 0 deletions
diff --git a/51/CH5/EX5.23/5_23.sce b/51/CH5/EX5.23/5_23.sce
new file mode 100755
index 000000000..1d3fcb915
--- /dev/null
+++ b/51/CH5/EX5.23/5_23.sce
@@ -0,0 +1,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")