summaryrefslogtreecommitdiff
path: root/551/CH13/EX13.26/26.sce
blob: bcf833d65fe7aa9842ef1b1b40841f664c49f599 (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
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
clc
p1=1; //bar
T1=363; //K
r=9; 
p3=68; //bar
p4=68; //bar
Q=1750; //kJ/kg
y=1.4;
cv=0.71;
cp=1.0;

disp("(i) Pressures and temperatures at salient points")
p2=p1*(r)^y;
disp("p2=")
disp(p2)
disp("bar")

T2=T1*r^(y-1);
disp("T2=")
disp(T2)
disp("K")

disp("p3=")
disp(p3)
disp("bar")

disp("p4=")
disp(p4)
disp("bar")

T3=T2*(p3/p2);
disp("T3=")
disp(T3)
disp("K")

Q1=cv*(T3-T2); //heat added at constant volume
Q2=Q-Q1; //heat added at constant pressure

T4=Q2/cp+T3;
disp("T4=")
disp(T4)
disp("K")

rho=T4/T3; //V4/V3=T4/T3

p5=p4*(rho/r)^y;
disp("p5=")
disp(p5)
disp("bar")

T5=T4*(rho/r)^(y-1);
disp("T5=")
disp(T5)
disp("K")


disp("(ii) Air standard efficiency =")
Qr=cv*(T5-T1);
n_airstandard=(Q-Qr)/Q;
disp(n_airstandard)


disp("(iii) Mean effective pressure =")
pm=1/(r-1)*(p3*(rho-1) + (p4*rho-p5*r)/(y-1) - (p2-p1*r)/(y-1));
disp(pm)
disp("bar")