blob: 5acd314c1d952dd1d468a57c9c902a03e8962d31 (
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
|
clc
T1=293; //K
T5=1023; //K
T7=T5;
p1=1.5; //bar
p2=6; //bar
n_c=0.82;
n_t=0.82;
e=0.70;
P=350; //kW
cp=1.005; //kJ/kg K
y=1.4;
T3=T1;
px=sqrt(p1*p2);
T2=T1*(px/p1)^((y-1)/y);
T2a=T1+(T2-T1)/n_c;
T4a=T2a;
p5=p2;
T6=T5/(p5/px)^((y-1)/y);
T6a=T5-n_t*(T5-T6);
T8a=T6a;
Ta=T4a+e*(T8a-T4a);
W_net=2*cp*[(T5-T6a)-(T2a-T1)];
Q1=cp*(T5-T4a)+cp*(T7-T6a); //Without regenerator
Q2=cp*(T5-Ta)+cp*(T7-T6a);
disp("n_thermal without regenerator =")
n1=W_net/Q1*100;
disp(n1)
disp("%")
disp("n_thermal woth regenerator =")
n2=W_net/Q2*100;
disp(n2)
disp("%")
disp("(iii) Mass of fluid circulated =")
m=P/W_net;
disp(m)
disp("kg/s")
|
