blob: c64e5a791edff042e782d84d6e2ba16c75ba99be (
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 5, Page 60
//Title: Gas Power Cycles
//=============================================================================
clc
clear
//INPUT DATA
rv=8;//Compression ratio
P1=95;//Pressure at point 1 in kPa
T1=300;//Temperature at point 1 in K
q23=750;//Heat transferred during constant volume heat addition process in kJ/kg
y=1.4;//Ratio of specific heats
Cv=0.718;//Specific heat at constant volume in kJ/kg-K
R=287;//Universal gas constant in J/kg-K
//CALCULATIONS
T2=T1*(rv^(y-1));//Temperature at point 2 in K
P2=P1*(rv^y);//Pressure at point 2 in kPa
T3=(q23/Cv)+T2;//Temperature at point 3 in K
P3=P2*(T3/T2);//Pressure at point 3 in kPa
nth=(1-(1/(rv^(y-1))))*100;//Thermal efficiency
Wnet=(nth*q23)/100;//Net work output in kJ/kg
v1=(R*T1)/(P1*1000);//Speific volume at point 1 in (m^3)/kg
MEP=Wnet/(v1*(1-(1/rv)));//Mean effective pressure in kPa
//OUTPUT
mprintf('Pressure at the end of heat addition process is %3.1f kPa \n Temperature at the end of heat addition process is %3.1f K \n Net work output is %3.2f kJ/kg \n Thermal efficiency is %3.2f percent \n Mean effective pressure is %3.0f kPa',P3,T3,Wnet,nth,MEP)
//==============================END OF PROGRAM=================================
|
