1 files changed, 32 insertions, 0 deletions
diff --git a/3685/CH19/EX19.24/Ex19_24.sce b/3685/CH19/EX19.24/Ex19_24.sce
new file mode 100644
index 000000000..e1007ffe5
--- /dev/null
+++ b/3685/CH19/EX19.24/Ex19_24.sce
@@ -0,0 +1,32 @@
+clc
+// Given that
+N0 = 9000 // Rotational speed in RPM
+Q = 6 // Volume flow rate in m^3/s
+p1 = 1 // Initial pressure in bar
+t1 = 25 // Initial temperature in degree centigrade
+p2 = 2.2 // Compressed pressure in bar
+n = 1.33 // Compression index
+Vf = 75 // Velocity of flow in m/s
+beta1 = 30 // Blade angle at inlet in degree
+beta2 = 55 // Blade angle at outlet in degree
+d = 0.75 // Diameter of impeller in m
+cp = 1.005 
+printf("\n Example 19.24\n")
+T1 = t1+273
+T2 = T1*(p2/p1)^((n-1)/n)
+Wc = cp*(T2-T1)
+x = Wc // Where x = Vw2*Vb2
+y = Vf/tand(beta2)// Where y = Vb2-Vw2(Equation 1)
+z = (y^2 +4*x*1000)^(0.5) // Where z = Vw2+Vb2(Equation 2)
+// By solving Equation 1 and Equation 2
+Vb2 = (y+z)/2
+Vw2 = ((z-y)/2)
+N = Vb2*60/(%pi*d)
+Vb1 = Vf/tand(beta1)
+D1 = Vb1*60/(%pi*N)
+b1 = Q/(%pi*D1*Vf)
+Q_ = Q* (1/p2)*(T2/T1)
+b2 = Q_/(%pi*d*Vf)
+printf("\n Speed of impeller = %f RPM,\n Impeller width at inlet = %f cm,\n Impeller width at outlet = %f cm,",N,b1*100,b2*100)
+// The answers given in the book vary because of round off error
+