initial commit / add all books

1 files changed, 33 insertions, 0 deletions

diff --git a/1955/CH5/EX5.6/example6.sce b/1955/CH5/EX5.6/example6.sce
new file mode 100755
index 000000000..3d24ff7a5
--- /dev/null
+++ b/1955/CH5/EX5.6/example6.sce
@@ -0,0 +1,33 @@
+clc

+clear

+//input data

+R=0.5//Degree of reaction

+P1=14//Initial pressure in bar

+T1=588//Initial temperature in K

+P2=0.14//Final pressure in bar

+Ns=0.75//Stage efficiency 

+RF=1.04//Reheat factor 

+N=20//No. of stages

+W=11770//Total power output in kW

+a1=20//Exit blade angle in degree

+hD=1/12//Ratio of blade height to blade mean diameter 

+

+//calculations

+hs1=3080//Isentropic enthalpy at initial condition from mollier chart in kJ/kg

+hs2=2270//Isentropic enthalpy at final condition from mollier chart in kJ/kg

+dhs=hs1-hs2//Isentropic enthalpy change in kJ/kg

+Nt=Ns*RF//Overall efficiency

+dh=Nt*dhs//Actual enthalpy drop in kJ/kg

+hs=dh/N//Enthalpy drop per stage in kJ/kg

+m=W/dh//Mass flow rate in kg/s

+C11=1.43*1//Velocity of steam at exit from nozzle in m/s in terms of U for 0.5 degree of reaction

+Wm=1*((2*C11*sind(90-a1))-1)//Work done per mass of steam in terms of U^2 in kJ/kg

+U=((hs*10^3)/Wm)^(1/2)//Mean blade velocity in m/s  as work done equals enthalpy drop per stage 

+C1=1.43*U//Velocity of steam at exit from nozzle in m/s 

+Ca=C1*cosd(90-a1)//Flow velocity in m/s

+v=1.618//Specific volume of steam from steam table at 1.05 bar with dry saturated steam in m^3/kg

+D=((m*v)/(hD*3.14*Ca))^(1/2)//Blade mean diameter in m

+N=(U*60)/(3.14*D)//Rotor speed in rpm

+

+//output

+printf('(a)Mass flow rate of steam is %3.2f kg/s\n(b)Mean blade velocity is %3.1f m/s \n(c)Blade mean diameter is %3.3f m \n(d)Rotor speed is %3i rpm',m,U,D,N)