blob: 12e7d01bde9ae41de023579f0af2e3dfee9bc511 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
|
clc
clear
//Input data
p1=1.05//Inlet pressure in kg/cm^2
T1=15+273//Inlet temperature in K
p2=33.4//Pressure at the end of adiabatic compression in kg/cm^2
r=5//The ratio of expansion
Cp=0.238//Specific heat at constant pressure in kJ/kg.K
Cv=0.17//Specific heat at constant volume in kJ/kg.K
g=1.4//Ratio of specific heats
//Calculations
r1=(p2/p1)^(1/g)//Compression ratio
k=r1/r//Cutoff ratio
T2=(p2/p1)^((g-1)/g)*T1//Temperature at the end of adiabatic compression in K
T3=T2*k//Temperature at the end of constant pressure heat addition in K
T4=T3*(1/r)^(g-1)//Temperature at the end of adiabatic expansion in K
qa=(Cp*(T3-T2))//Heat added in kcal/kg of air
qre=(Cv*(T4-T1))//Heat rejected in kcal/kg of air
nt=((qa-qre)/qa)*100//Ideal thermal efficiency in percent
//Output
printf('The ideal thermal efficiency is %3.1f percent',nt)
|
