blob: 3002c430c4451f6d54bb2da11c93aa3dab43d84c (
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
36
37
|
//Chapter-2, Illustration 4, Page 58
//Title: Gas Power Cycles
//=============================================================================
clc
clear
//INPUT DATA
rv=9.5;//Compression ratio
P1=100;//Air pressure in kPa
T1=290;//Air temperature in K
V1=600*(10^-6);//Volume of air in m^3
T4=800;//Final temperature in K
R=287;//Universal gas constan in J/kg.K
Cv=0.718;//Specific heat at constant volume in kJ/kg.K
y=1.4;//Ratio of specific heats
//CALCULATIONS
T3=T4*(rv^(y-1));//Temperature at the end of constant volume heat addition in K
P2=P1*(rv^y);//Pressure at point 2 in kPa
T2=T1*(rv^(y-1));//Temperature at point 2 in K
P3=P2*(T3/T2);//Pressure at point 3 in kPa
m=(P1*1000*V1)/(R*T1);//Specific mass in kg/s
Q=m*Cv*(T3-T2);//Heat transferred in kJ
n=(1-(1/rv^(y-1)))*100;//Thermal efficiency
Wnet=(n*Q)/100;//Net workdone in kJ
MEP=Wnet/(V1*(1-(1/rv)));//Mean effective pressure in kPa
//OUTPUT
mprintf('Maximum pressure of the cycle is %3.2f kPa \n Maximum temperature of the cycle is %3.1f K \n Amount of heat transferred is %3.2f kJ \n Thermal efficiency is %3.1f percent \n Mean effective pressure is %3.1f kPa',P3,T3,Q,n,MEP)
//==============================END OF PROGRAM=================================
|
