diff options
Diffstat (limited to '3685/CH13/EX13.4/Ex13_4.sce')
| -rw-r--r-- | 3685/CH13/EX13.4/Ex13_4.sce | 30 |
1 files changed, 30 insertions, 0 deletions
diff --git a/3685/CH13/EX13.4/Ex13_4.sce b/3685/CH13/EX13.4/Ex13_4.sce new file mode 100644 index 000000000..964a442b1 --- /dev/null +++ b/3685/CH13/EX13.4/Ex13_4.sce @@ -0,0 +1,30 @@ +clc
+T1 = 50 // Temperature before compression stroke in degree Celsius
+rk = 16 // Compression ratio
+g = 1.4 // Heat capacity ratio
+P3 = 70 // Maximum cycle pressure in bar
+cv = 0.718 // Constant volume heat addition capacity
+cp = 1.005 // Constant pressure heat addition capacity
+R = 0.287 // Gas constant
+T2 = (T1+273)*((rk^(g-1))) //Temperature after compression stroke
+P1 = 1 // Pressure before compression in bar
+P2 = P1*(rk)^g // Pressure after compression
+T3 = T2*(P3/P2) // Temperature after constant volume heat addition
+Q23 = cv*(T3-T2) // Constant volume heat added
+T4 = (Q23/cp)+T3 // Temperature after constant pressure heat addition
+v43 = T4/T3 // cut off ratio
+v54 = rk/v43 // Expansion ratio
+T5 = T4*(1/v54)^(g-1) // Temperature after expansion
+P5 = P1*(T5/(T1+273)) // Pressure after expansion
+Q1 = cv*(T3-T2)+cp*(T4-T3) // Total heat added
+Q2 = cv*(T5-T1-273) // Heat rejected
+n_cycle = 1-(Q2/Q1) // Cycle efficiency
+v1 = (R*(T1+273))/(P1*1e2) // Volume before compression
+v2 = (1/16)*v1 // Swept volume
+Wnet = Q1*n_cycle // Net work done
+Pm = Wnet/(v1-v2) // Mean pressure
+printf("\n Example 13.4\n")
+printf("\n Efficiency of the cycle is %f percent",n_cycle*100)
+printf("\n Mean effective pressure is %f bar",Pm/100)
+//The answers vary due to round off error
+
|