diff options
author | prashantsinalkar | 2017-10-10 12:27:19 +0530 |
---|---|---|
committer | prashantsinalkar | 2017-10-10 12:27:19 +0530 |
commit | 7f60ea012dd2524dae921a2a35adbf7ef21f2bb6 (patch) | |
tree | dbb9e3ddb5fc829e7c5c7e6be99b2c4ba356132c /3685/CH19/EX19.24/Ex19_24.sce | |
parent | b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (diff) | |
download | Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.tar.gz Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.tar.bz2 Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.zip |
initial commit / add all books
Diffstat (limited to '3685/CH19/EX19.24/Ex19_24.sce')
-rw-r--r-- | 3685/CH19/EX19.24/Ex19_24.sce | 32 |
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
+
|