blob: 63d6ad748443e99fd7ae8bfb597772d363dcf425 (
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
//Initialization of variables
ta=780 //F
tr=80//F
Qa=195 //Btu
pd=14.7 //psia
R=1.986/29
k=1.4
J=778
g=32.174
//calculations
etat=(ta-tr)/(ta+459.7)
W=etat*Qa
vd=R*(tr+460)/pd
va=vd*exp(-(Qa-W)/R/(tr+460))
vb=va*((tr+460)/(ta+460))^(1/(k-1))
vc=vd/va*vb
rv=vd/vb
rv2=vc/vb
rv3=va/vb
imep= W*J/(144*(vd-vb))/5.77
//5.77 is conversion factor
//results
printf("cycle expansion ratio = %.1f ",rv)
printf("\n isothermal expansion ratio = %.1f",rv2)
printf("\n isentropic expansion ratio = %.1f",rv3)
printf("\n imep = %.1f lbf/in^2",imep)
|
