blob: 5109df3e6ae19a8be5d3e94574c952843e6ae14a (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
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
|
