1 files changed, 40 insertions, 0 deletions
diff --git a/1775/CH4/EX4.7/Chapter4_Example7.sce b/1775/CH4/EX4.7/Chapter4_Example7.sce
new file mode 100755
index 000000000..22c03f53c
--- /dev/null
+++ b/1775/CH4/EX4.7/Chapter4_Example7.sce
@@ -0,0 +1,40 @@
+//Chapter-4, Illustration 7, Page 170
+//Title: Steam Nozzles and Steam Turbines
+//=============================================================================
+clc
+clear
+
+//INPUT DATA
+P1=10;//Pressure at point 1 in bar
+P2=0.5;//Pressure at point 2 in bar
+h1=3050;//Enthalpy at point 1 in kJ/kg
+h2s=2480;//Enthalpy at point 2s in kJ/kg
+ht=2910;//Enthalpy at throat in kJ/kg
+n=1.3;//Adiabatic gas constant
+r=0.1;//Total available heat drop
+v1=0.258;//Specific volume at point 1 in (m^3)/kg
+h2f=340.6;//Enthalpy for exit pressure from steam tables in kJ/kg
+hfg=2305.4;//Enthalpy for exit pressure from steam tables in kJ/kg
+m=0.5;//Mass flow rate in kg/s
+
+//CALCULATIONS
+x=n/(n-1);//Ratio
+Pt=((2/(n+1))^x)*P1;//Temperature at throat in bar
+h2=h2s+(r*(h1-h2s));//Enthalpy at point 2 in kJ/kg
+vt=((P1/Pt)^(1/n))*v1;//Specific volume at throat in (m^3)/kg
+v2=((P1/P2)^(1/n))*v1;//Specific volume at point 2 in (m^3)/kg
+Ct=sqrt(2000*(h1-ht));//Throat velocity in m/s
+At=((m*vt)/Ct)*(10^6);//Throat area in (mm^2)
+C2=sqrt(2000*(h1-h2));//Exit velocity in m/s
+A2=((m*v2)/C2)*(10^6);//Exit area in (mm^2)
+x2=((h2-h2f)/hfg)*100;//Steam quality at exit
+
+//OUTPUT
+mprintf('Throat area is %3.0f (mm^2) \n Exit area is %3.0f (mm^2) \n Steam quality at exit is %3.0f percent',At,A2,x2)
+
+
+
+
+
+
+//==============================END OF PROGRAM=================================