summaryrefslogtreecommitdiff
path: root/3751/CH11/EX11.13/Ex11_13.sce
diff options
context:
space:
mode:
Diffstat (limited to '3751/CH11/EX11.13/Ex11_13.sce')
-rw-r--r--3751/CH11/EX11.13/Ex11_13.sce56
1 files changed, 56 insertions, 0 deletions
diff --git a/3751/CH11/EX11.13/Ex11_13.sce b/3751/CH11/EX11.13/Ex11_13.sce
new file mode 100644
index 000000000..0dde95b26
--- /dev/null
+++ b/3751/CH11/EX11.13/Ex11_13.sce
@@ -0,0 +1,56 @@
+//Fluid Systems - By - Shiv Kumar
+//Chapter 11- Centrifugal Pumps
+//Example 11.13
+//To Find (i)The Discharge of the Pump (ii)The Pressure at Suction and Delivery side of the Pump.
+
+ clc
+ clear
+
+//Given Data:-
+ h_st=35; //Static Head, m
+ h_s=4; //Suction Head, m
+ D=150; //Diameter of Pipes, mm
+ Ds=D; //Diameter of Suction Pipe, mm
+ Dd=D; //Diameter of Delivery Pipe, mm
+ h_fs=1.6; //Head loss in Suction pipe, m
+ h_fd=6.5; //Head loss in Delivery Pipe, m
+ Do=380; //Diameter of Impeller at Outlet, mm
+ bo=25; //Width of Impeller at Outlet, mm
+ N=1200; //Speed, rpm
+ beta_o=35; //Ezxit Blade Angle, degrees
+ eta_man=80/100; //Manometric Efficiency
+
+
+//Data Used:-
+ g=9.81; //Acceleration due to gravity, m/s^2
+
+//Computations:-
+ Do=Do/1000; //m
+ D=D/1000; //m
+ Ds=Ds/1000; //m
+ Dd=Dd/1000; //m
+ bo=bo/1000; //m
+
+ Hm=h_st+h_fs+h_fd; // Manometric Head, m
+ uo=%pi*Do*N/60; // Tangential velocity of Impeller at Outlet, m/s
+ Vwo=g*Hm/(uo*eta_man); //m/s
+ Vfo=(uo-Vwo)*tand(beta_o); //m/s
+
+ //(i)The Discharge of the Pump, Q
+ Q=%pi*Do*bo*Vfo*1000; //litres/s
+
+ // (ii)The Pressure at Suction and Delivery side of the Pump
+
+ A=(%pi/4)*D^2; //m^2
+ Vd=Q*10^-3/A; //m/s
+ Vs=Vd; //m/s
+ Hs=h_s+h_fs+Vs^2/(2*g); //Pressure on Suction Side, m of water
+ h_d=h_st-h_s; //m
+ Hd=h_d+h_fd+Vd^2/(2*g); //Pressure on Delivery Side, m of water
+
+
+//Result:-
+ printf(" (i)The Discharge of the Pump, Q =%.2f litres/s\n",Q) //The answer vary due to round off error
+ printf(" (ii) Pressure on Suction Side, Hs =-%.3f m of water \n",Hs) //The answer vary due to round off error
+ printf(" Pressure on Delivery Side, Hd =%.2f m of water \n",Hd) //The answer vary due to round off error
+