1 files changed, 33 insertions, 0 deletions
diff --git a/629/CH14/EX14.10/example14_10.sce b/629/CH14/EX14.10/example14_10.sce
new file mode 100644
index 000000000..448ee4aee
--- /dev/null
+++ b/629/CH14/EX14.10/example14_10.sce
@@ -0,0 +1,33 @@
+clear
+clc
+//Example 14.10 IMPULSE TURBINE
+z1=1670; //[m]
+z2=1000; //[m]
+g=9.81; //[m/s^2]
+L=6000; //length[m]
+f=0.015;
+Dj=0.18; //diameter [m]
+Aj=%pi*Dj^2/4 //area[m^2]
+Dp=1; //diameter of penstock[m]
+Ap=%pi*Dp^2/4 //area[m^2
+//Energy equation, (p1/gamma)+(V1^2/2g)+z1=(p2/gamma)+(Vj^2/2g)+z2+hL
+//p1=0,p2=0 ,Vp=Vj*(Aj/Ap), hL=f*L*Vp^2/(D*2*g)
+Vj=sqrt(2*g*(z1-z2)/(1+f*L*(Aj/Ap)^2/Dp)) //jet velocity[m/s]
+Gamma=9810; 
+Q=Aj*Vj //[m^3/s])
+P=Q*Gamma*Vj^2/(2*g*10^3) //gross power[kW]
+eta=0.85; //efficiency
+Pd=P*eta //power developed [kW]
+printf("\nThe power developed by the turbine = %.f kW.\n",Pd)
+
+//Angular speed of wheel
+r=1.5; //radius[m]
+Vb=Vj/2 //[m/s]
+w=Vb/r //[rad/s]
+//Wheel speed
+N=w*60/(2*%pi) //in rpm
+printf("\nAngular speed of wheel for maximum efficiency = %.f rpm.\n",N)
+
+//Torque
+T=Pd/w //[kN.m]
+printf("\nThe torque on the turbine shaft = %.f kN.m\n",T)
+\ No newline at end of file