summaryrefslogtreecommitdiff
path: root/551/CH4/EX4.53/53.sce
blob: ab3be111f80be8c2d0675bbf3fd4b51ee510c6b8 (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
m=1; //kg
p=8; //bar
s1=6.55; //kJ/kg.K
T=200; //0C
s_f1=2.0457; //kJ/kg.K
s_fg1=4.6139; //kJ/kg.K
h_f1=720.9; //kJ/kg
h_fg1=2046.5; //kJ/kg
h2=2839.3; //kJ/kg

x1=(s1-s_f1)/s_fg1;
h1=h_f1+x1*h_fg1;
Q=h2-h1;
disp("Heat supplied=")
disp(Q)
disp("kJ/kg")

// For T-s diagram

s=0:0.01:10;
T=(-(s-5)^2+298);
plot(s,T)

T=[295.44 295.44];
s=[6.6 3.45];
plot(s,T,'g')

s=[6.6 7];
T=[295.44 300];
plot(s,T,'g')

s=[6.55 6.55];
T=[270 295.44];
plot(s,T,'r')

s=[6.6 6.6];
T=[270 295.44];
plot(s,T,'--r')

s=[6.66 6.66];
T=[270 295.44];
plot(s,T,'r')

xtitle("T-s diagram", "s(kJ/kg K)", "T(K)")
//The area in red represents the heat flow and it goes upto x-axis