diff options
Diffstat (limited to '3685/CH13/EX13.3/Ex13_3.sce')
| -rw-r--r-- | 3685/CH13/EX13.3/Ex13_3.sce | 27 |
1 files changed, 27 insertions, 0 deletions
diff --git a/3685/CH13/EX13.3/Ex13_3.sce b/3685/CH13/EX13.3/Ex13_3.sce new file mode 100644 index 000000000..5109df3e6 --- /dev/null +++ b/3685/CH13/EX13.3/Ex13_3.sce @@ -0,0 +1,27 @@ +clc
+rk = 16 // Compression ratio
+T1 = 15 // Air inlet temperature in degree Celsius
+P1 = 0.1 // Air inlet pressure in MPa
+T3 = 1480 // Highest temperature in cycle in degree Celsius
+g = 1.4 // Heat capacity ratio
+R = 0.287 // Gas constant
+T2 = (T1+273)*(rk^(g-1)) // Temperature after compression
+rc = (T3+273)/T2 // cut off ratio
+cp = 1.005 // Constant pressure heat constant
+cv = 0.718 // Constant volume heat constant
+Q1 = cp*(T3+273-T2) // Heat addition
+T4 = (T3+273)*((rc/rk)^(g-1)) // Temperature after heat addition
+Q2 = cv*(T4-T1-273) // Heat rejection
+n = 1-(Q2/Q1) // cycle efficiency
+n_ = 1-((1/g))*(1/rk^(g-1))*((rc^(g-1))/(rc-1)) // cycle efficiency from another formula
+Wnet = Q1*n // Net work
+v1 = (R*(T1+273))/(P1*1e3) // Volume before compression
+v2 = v1/rk // Volume after compression
+Pm = Wnet/(v1-v2) // Mean pressure
+printf("\n Example 13.3\n")
+printf("\n Cut-off ratio is %f ",rc)
+printf("\n Heat supplied per kg of air is %f kJ/kg",Q1)
+printf("\n Cycle efficiency is %f percent",n*100)
+printf("\n Mean effective pressure is %f kPa",Pm)
+//The answers vary due to round off error
+
|