blob: 40d145d54bdd99f7e3bf6c9f4ad25cfed73af67e (
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
38
39
40
41
42
43
44
45
46
|
clc
r1=15.3; //V1/V2
r2=7.5; //V4/V3
p1=1; //bar
T1=300; //K
n_mech=0.8;
C=42000; //kJ/kg
y=1.4;
R=287;
cp=1.005;
cv=0.718;
V2=1; ////Assuming V2=1 m^3
T2=T1*r1^(y-1);
p2=p1*r1^y;
T3=r1/r2*T2;
m=p2*10^5*V2/R/T2;
T4=T3/r2^(y-1);
Q_added=m*cp*(T3-T2);
Q_rejected=m*cv*(T4-T1);
W=Q_added-Q_rejected;
pm=W/(r1-1)/V2/100;
disp("Mean effective pressure =")
disp(pm)
disp("bar")
ratio=p2/pm;
disp("Ratio of maximum pressure to mean effective pressure =")
disp(ratio)
n_cycle=W/Q_added;
disp("Cycle efficiency =")
disp(n_cycle)
n_thI=0.5;
n_cycle1=n_thI*n_cycle;
n_thB=n_mech*n_cycle1;
BP=1;
mf=BP/C/n_thB*3600;
disp("Fuel consumption per kWh =")
disp(mf)
disp("kg/kWh")
|
