1 files changed, 39 insertions, 0 deletions
diff --git a/1955/CH5/EX5.12/example12.sce b/1955/CH5/EX5.12/example12.sce
new file mode 100755
index 000000000..558f7de57
--- /dev/null
+++ b/1955/CH5/EX5.12/example12.sce
@@ -0,0 +1,39 @@
+clc
+clear
+//input data
+P0=100//Initial pressure of steam in bar
+T0=773//Initial temperature of steam in K
+D=1//Turbine diameter in m
+N=3000//Speed of turbine in rpm
+m=100//Mass flow rate of steam in kg/s
+a1=70//Exit angle of the first stage nozzle in degree
+ns1=0.78//Stage efficiency of first stage 
+ns2=ns1//Stage efficiency of second stage
+
+//calculations
+U=(3.1415*D*N)/60//Mean blade speed in m/s
+C1=(2*U)/sind(a1)//Velocity of steam at exit from nozzle in m/s 
+b11=atand(tand(a1)/2)//Inlet angle of the rotor blade in degree
+b21=b11//Outlet angle of the rotor blade in degree
+b12=b21//Inlet angle of the rotor blade in second stage in degree
+b22=b12//Outlet angle of the rotor blade in second stage in degree
+W=4*m*U^2*10^-6//Total work done in both the stages in MW
+dh02=2*U^2*10^-3//Change in enthalpy in first stage of turbine in kJ/kg
+dh02s=(dh02/ns1)//Change in enthalpy isentropically of turine first stage in kJ/kg
+printf('The values of enthalpy and specific volume are taken from the mollier chart at inlet and exit conditions respectively')
+h0=3370//Enthalpy at beginning of first stage in kJ/kg
+h2=h0-dh02//Enthalpy at the end of first stage in kJ/kg
+h2s=h0-dh02s//Isentropic enthalpy at the end of first stage in kJ/kg
+v2=0.041//Specific volume at the end of first stage in m^3/kg
+dh24=2*U^2*10^-3//Change in enthalpy in second stage of turbine in kJ/kg
+dh24s=dh24/ns2//Change in enthalpy isentropically of turine second stage in kJ/kg
+h4=h2-dh24//Enthalpy at beginning of second stage in kJ/kg
+h4s=h2-dh24s//Isentropic enthalpy at the end of second stage in kJ/kg
+v4=0.05//Specific volume at the end of second stage in m^3/kg
+
+Ca=C1*cosd(a1)//Axial velocity in m/s
+h1r=(m*v2)/(3.1415*D*Ca)//Blade height at first stage rotor exit in m
+h2r=(m*v4)/(3.1415*D*Ca)//Blade height at second stage rotor exit in m
+
+//output
+printf('\n\n(a)rotor blade angles\n    Inlet angle of the rotor blade is %3.2f degree\n    Outlet angle of the rotor blade is %3.2f degree\n    Inlet angle of the rotor blade in second stage is %3.2f degres\n    Outlet angle of the rotor blade in second stage is %3.2f degree\n(b)Total work done or Power developed in both the stages is %3.2f MW\n(c)final state of steam\n    Enthalpy at beginning of first stage is %3i kJ/kg\n    Enthalpy at the end of first stage is %3.2f kJ/kg\n    Isentropic enthalpy at the end of first stage is %3.2f kJ/kg\n    Specific volume at the end of first stage is %3.3f m^3/kg\n    Enthalpy at beginning of second stage is %3.1f kJ/kg\n    Isentropic enthalpy at the end of second stage is %3.2f kJ/kg\n    Specific volume at the end of second stage is %3.2f m^3/kg\n(d)blade height\n    Blade height at first stage rotor exit is %3.4f m\n    Blade height at second stage rotor exit is %3.4f m',b11,b21,b12,b22,W,h0,h2,h2s,v2,h4,h4s,v4,h1r,h2r)