1 files changed, 22 insertions, 0 deletions
diff --git a/51/CH12/EX12.6/12_6.sce b/51/CH12/EX12.6/12_6.sce
new file mode 100755
index 000000000..2d949a627
--- /dev/null
+++ b/51/CH12/EX12.6/12_6.sce
@@ -0,0 +1,22 @@
+clc;
+clear;
+exec("C:\Program Files\scilab-5.3.0\bin\TCP\12_6data.sci");
+//energy equation between a point on surface of lake and the nozzle outlet
+//z0=(V1^2)/(2*32.2) + hL
+//from continuity equation, V=(A1*V1/A)=((D1/D)^2)*V1
+//neglecting minor losses, 
+//z0=(1+(f*l/D)*((D1/D)^4))*(V1^2)/(2*32.2)
+//Wshaft=d*Q*u*(U-V1)*(1-cos(B))
+//The maximum power occurs at U=V1/2 and dWshaft/dD1=0
+a=(2*32.2*z0)^0.5;//113.5
+b=f*(l/(D/12))*(1/(D/12))^4;//152
+c=a*%pi*1.94*(1-cos(B*%pi/180))/4;//323
+d=(c*a*a/4);//1.04*(10^6)
+//by the above conditions, and applying Q=(%pi*(D1^2)*V1/4)
+D1=(1/(2*b))^(1/4);//ft
+disp("ft",D1,"The nozzle diameter for maximum power output=")
+Wshaft=-((d*D1^2)/(1+(b*D1^4))^(3/2))/550;//hp
+disp("hp",Wshaft,"The maximum power output=")
+V1=a/(1+(b*(D1^4)))^0.5;//ft/sec
+omega=(V1/(2*R))*60/(2*%pi);//rpm
+disp("rpm",omega,"The angular velocity of the rotor=")