1 files changed, 34 insertions, 0 deletions
diff --git a/3733/CH4/EX4.16/Ex4_16.sce b/3733/CH4/EX4.16/Ex4_16.sce
new file mode 100644
index 000000000..e04f3f3ad
--- /dev/null
+++ b/3733/CH4/EX4.16/Ex4_16.sce
@@ -0,0 +1,34 @@
+//Example 4_16
+clc;funcprot(0);
+//Given data
+P=5400;// kW
+N=200;// r.p.m
+D=3;// m
+H=240;// m
+n_t=0.82;
+rho=1000;// kg/m^3
+g=9.81;// m/s^2
+
+//Calculation
+//(a)
+Q=(P*1000)/(rho*g*H*n_t);// m^3/sec
+N_u=(N*D/sqrt(H));// Unit speed
+P_u=(P/(D^2*H^(3/2)));// Unit power
+Q_u=(Q/(D^2*(sqrt(H))));// Unit flow
+N_s=(2*%pi*N*sqrt(Q)*sqrt(n_t))/(60*(g*H)^(3/4));// Specific speed
+printf('\n(a)The flow rate,Q=%0.1f m^3/sec \n   The unit speed,N_u=%0.1f \n   The unit power,P_u=%0.3f \n   The unit flow,Q_u=%0.2f \n   The specific speed,N_s=%0.3f',Q,N_u,P_u,Q_u,N_s);
+//(b)
+// When the head is changed to 160 m,the diameter remains same.
+H=160;// m
+N=(N_u*sqrt(H))/D;// rpm
+P_1=(P_u*D^2*H^(3/2));// kW
+Q=(Q_u*D^2*sqrt(H));// m^3/sec
+printf('\n(b)Speed,N=%0.0f r.p.m\n   Power,P=%0.0f kW\n   The flow rate,Q=%0.2f m^3/sec ',N,P_1,Q);
+//(c)
+H=183;// m
+P=2850;// kW
+D_1=sqrt((P/(P_u*H^(3/2))));// m
+N_1=(N_u*sqrt(H))/(D_1);// r.p.m
+Q=Q_u*D_1^2*sqrt(H);// m^3/sec
+printf('\n(c)Diameter,D_1=%0.2f m\n   Speed,N=%0.0f r.p.m\n   The flow rate,Q=%0.2f m^3/sec ',D_1,N_1,Q);
+// The answer vary due to round off error