1 files changed, 32 insertions, 0 deletions
diff --git a/2063/CH8/EX8.5/8_5.sce b/2063/CH8/EX8.5/8_5.sce
new file mode 100755
index 000000000..5ad0dd1d3
--- /dev/null
+++ b/2063/CH8/EX8.5/8_5.sce
@@ -0,0 +1,32 @@
+clc
+clear
+//Input data
+P1=7;//Inlet pressure of a convergent divergent steam nozzle in bar
+T1=275;//Inlet temperature of the nozzle in degrees celcius
+P2=1;//Discharge pressure of steam in bar
+l=60;//Length of diverging portion of the nozzle in mm
+dt=6;//Diameter of the throat in mm
+f1=10;//Percent of total available enthalpy drop lost in friction in the diverging portion in percentage
+h1=3006.9;//Enthalpy of steam at 7bar pressure and 275 degrees celcius in kJ/kg
+ht=2865.9;//Enthalpy at the throat from Moiller chart in kJ/kg
+he=2616.7;//Enthalpy at the exit from moiller chart in kJ/kg
+vt=0.555;//Specific volume of steam at throat in m^3/kg
+Tt=202.8;//Temperature of steam at throat in degrees celcius from moiller chart
+ve=1.65;//Volume of steam at exit in m^3/kg
+
+//Calculations
+Pt=0.546*P1;//The throat pressure for maximum discharge in bar
+hd=h1-ht;//Enthalpy drop upto throat in kJ/kg
+Vt=44.7*(hd)^(0.5);//Velocity of steam at throat in m/s
+hid=h1-he;//Total isentropic drop from 7 bar,275 degrees celcius to 1 bar in kJ/kg
+hda=(1-(f1/100))*(hid);//Actual heat drop in kJ/kg
+Ve=44.7*(hda)^(0.5);//Velocity at exit in m/s
+At=(3.14/4)*(6/1000)^(2);//Throat area of the nozzle in m^2
+m=(At*Vt)/vt;//Mass flow rate at nozzle throat in kg/s
+Ae=((m*ve)/Ve)*10^4;//Exit area of the nozzle in cm^2
+de=(((Ae*4)/3.14)^(0.5))*10;//Diameter of the nozzle at exit in mm
+alpha=atand((de-dt)/(2*60));//Half of the cone angle of the nozzle in degrees
+alpha1=2*alpha;//Cone angle of the nozzle in degrees
+
+//Output
+printf('(a)Velocity of steam at throat is %3.0f m/s\n (b)Temperature of steam at the throat is %3.1f degrees celcius\n (c)Cone angle of the divergent portion is %3.3f degrees',Vt,Tt,alpha1)