1 files changed, 45 insertions, 0 deletions
diff --git a/1943/CH7/EX7.15/Ex7_15.sce b/1943/CH7/EX7.15/Ex7_15.sce
new file mode 100755
index 000000000..84fdbc19a
--- /dev/null
+++ b/1943/CH7/EX7.15/Ex7_15.sce
@@ -0,0 +1,45 @@
+
+clc
+clear
+//Input data
+p1=35//Pressure in bar
+T1=350+273//Temperature in K
+p2=0.07//Pressure in bar
+x=1/4//Fraction of drop in isentropic enthalpy
+a=20//Nozzle angle in degrees
+nn=88//Nozzle efficiency in percent
+y=0.2//Velocity ratio
+b2=30//Exit blade angle in degrees
+b4=30//Exit blade angle in degrees
+f=0.9//Friction coefficienct
+in=75//Internal efficiency of the turbine in percent
+
+//Calculations
+h1=3106.4//Enthalpy in kJ/kg
+s1=6.6643//Entropy in kJ/kg.K
+x2s=(s1-0.5582)/7.7198//dryness fraction
+h2s=(163.16+x2s*2409.54)//Enthalpy in kJ/kg
+dh=(h1-h2s)//Change in enthalpy in kJ/kg
+h13s=x*dh//Change in enthalpy in kJ/kg
+h13=(nn/100)*h13s//Change in enthalpy in kJ/kg
+V1=(44.72*sqrt(h13))//Velocity in m/s
+Vb=(y*V1)//Velocity in m/s
+b1=atand((V1*sind(a))/((V1*cosd(a))-Vb))//Angle in degrees
+Vr1=((V1*sind(a))/sind(b1))//Velocity in m/s
+Vr2=(f*Vr1)//Velocity in m/s
+dVw1=(Vr1*cosd(b1))+(Vr2*cosd(b2))//Velocity in m/s
+V2=sqrt((Vr2*sind(b2))^2+((Vr2*cosd(b2))-Vb)^2)//Velocity in m/s
+V3=f*V2//Velocity in m/s
+b3=atand((V3*sind(b2))/((V3*cosd(b2))-Vb))//Angle in degrees
+Vr3=((V3*sind(b2))/sind(b3))//Velocity in m/s
+Vr4=f*Vr3//Velocity in m/s
+dVw2=(Vr3*cosd(b3))+(Vr4*cosd(b4))//Velocity in m/s
+dVw=(dVw1+dVw2)//Velocity in m/s
+nb1=((2*dVw*Vb)/V1^2)*100//Efficiency in percent
+ns=(nn*nb1)/100//Efficiency in percent
+ht=(in/100)*dh//Total change in enthalpy in kJ/kg
+hc=(ns/100)*h13s//Total change in enthalpy in kJ/kg
+pp=(hc/ht)*100//Percentage of enthalpy
+
+//Output
+printf('Efficiency of first stage is %3.2f percent \n Percentage of the total power developed by the turbine is %3.2f percent',ns,pp)