1 files changed, 25 insertions, 0 deletions
diff --git a/3685/CH13/EX13.1/Ex13_1.sce b/3685/CH13/EX13.1/Ex13_1.sce
new file mode 100644
index 000000000..5a50b95bf
--- /dev/null
+++ b/3685/CH13/EX13.1/Ex13_1.sce
@@ -0,0 +1,25 @@
+clc
+T1 = 35 // Air inlet temperature in degree Celsius
+P1 = 0.1 // Air inlet pressure in MPa
+Q1 = 2100 // Heat supply in kJ/kg
+R = 0.287 // gas constant
+rk = 8 // Compression ratio
+g = 1.4 // Heat capacity ratio
+n_cycle = 1-(1/rk^(g-1)) // cycle efficiency 
+v1 = (R*(T1+273))/(P1*1e3) // Initial volume
+v2 = v1/8 // Volume after compression
+T2 = (T1+273)*(v1/v2)^(g-1) // Temperature after compression
+cv = 0.718 // Constant volume heat capacity in kJ/kg
+T3 = Q1/cv + T2 // Temperature at after heat addition
+P21 = (v1/v2)^g // Pressure ratio
+P2 = P21*P1 // Pressure after compression
+P3 = P2*(T3/T2) // Pressure after heat addition
+Wnet = Q1*n_cycle // Net work output
+Pm = Wnet/(v1-v2) // Mean pressure
+printf("\n Example 13.1\n")
+printf("\n Cycle efficiency is %f percent",n_cycle*100)
+printf("\n Maximum temperature in the cycle is %d K",T3)
+printf("\n Maximum pressure in the cycle is %f MPa",P3)
+printf("\n Mean effective pressure is %f MPa",Pm/1e3)
+//The answers vary due to round off error
+