blob: a3d2400afd4a7dfdc7469ba0a193ffb3b9759ba0 (
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
28
29
30
31
32
33
34
35
|
//Chapter-2, Illustration 10, Page 66
//Title: Gas Power Cycles
//=============================================================================
clc
clear
//INPUT DATA
rv=20;//Compression ratio
P1=95;//Pressure at point 1 in kPa
T1=293;//Temperature at point 1 in K
T3=2200;//Temperature at point 3 in K
y=1.4;//Ratio of specific heats
R=287;//Universal gas constant in J/kg-K
Cp=1.005;//Specific heat at constant pressure in kJ/kg-K
//CALCULATIONS
P2=P1*(rv^y);//Pressure at point 2 in kPa
T2=T1*(rv^(y-1));//Temperature at point 2 in K
v2=(R*T2)/(P2*1000);//Specific volume at point 2 in (m^3)/kg
v3=v2*(T3/T2);//Specific volume at point 3 in (m^3)/kg
rc=v3/v2;//Cut-off ratio
nth=(1-(((rc^y)-1)/((rv^(y-1))*y*(rc-1))))*100;//Thermal efficiency
q23=Cp*(T3-T2);//Heat flow between points 2 and 3 in kJ/kg
wnet=(nth*q23)/100;//Net workdone in kJ/kg
MEP=wnet/(v2*(rv-1));//Mean effective pressure in kPa
//OUTPUT
mprintf('Thermal efficiency is %3.1f percent \n Mean effective pressure is %3.2f kPa',nth,MEP)
//==============================END OF PROGRAM=================================
|
