1 files changed, 41 insertions, 0 deletions
diff --git a/32/CH13/EX13.10/13_10.sce b/32/CH13/EX13.10/13_10.sce
new file mode 100755
index 000000000..5c8770283
--- /dev/null
+++ b/32/CH13/EX13.10/13_10.sce
@@ -0,0 +1,41 @@
+//pathname=get_absolute_file_path('13.10.sce')
+//filename=pathname+filesep()+'13.10-data.sci'
+//exec(filename)
+//Mass flow rate(in kg/s):
+m=5/60
+//Pressure at which steam is discharged(in bar):
+p3=1
+//Initial pressure(in bar):
+p1=10
+//Initial temperature(in K)
+T1=200+273
+//Adiabatic index of compression:
+n=1.3
+//From steam tables:
+h1=2827.9 //kJ/kg
+s1=6.6940 //kJ/kg.K
+v1=0.2060 //m^3/kg
+h2a=2711.23 //kJ/kg
+s2a=6.6749 //kJ/kg.K
+s3=s2a
+h3=2420.08 //kJ/kg
+v3=1.5025 //m^3/kg
+psat=3.44 //bar (at T=138.18 C)
+Tsat=155.12 //C (at p=5.45 bar)
+//Pressure at throat(in bar):
+p2=p1*(2/(n+1))^(n/(n-1))
+//Velocity at exit(in m/s):
+C3=sqrt(2*(h1-h3)*10^3)
+//Exit area(in m^2):
+A3=m*v3/C3
+//Diameter of nozzle at exit(in m):
+d=sqrt(A3*4/%pi)
+//Temperature at throat(in K):
+T2=T1*(p2/p1)^((n-1)/n)
+//Degree of supersaturation:
+d1=p2/psat
+//Amount of undercooling(in C):
+u=Tsat-(T2-273)
+printf("\nRESULT\n")
+printf("\nDegree of supersaturation = %f",d1)
+printf("\nAmount of undercooling = %f C",u)
+\ No newline at end of file