blob: 449342b19582d10e142ed521075db255556079cf (
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
|
clc
clear
//Initial conditions at 7 bar pressure
P1=7; //in bars
Vg1=0.273; //in m^3/kg
V1=Vg1; //in m^3/kg
Hg1=2763.5; //in kJ/kg
H1=Hg1;
Tsat=165+273; //in K
Sf=1.992; //in kJ/kg K
Sfg=4.716; //in kJ/kg K
Sg=6.708; //in kJ/kg K
n=1.1;
//Final conditions at 0.5 bar
P2=0.5; //in bars
V2=((P1*(V1^1.1))/P2)^(1/1.1); //using P(V)^1.1=Constant
W=((P1*100*V1)-(P2*100*V2))/(n-1);
printf('Work Done: %3.2f kJ',W);
printf('\n');
Hf2=340.6; //in kJ/kg
Hfg2=2305.4; //in kJ/kg
Vg2=3.24; //in m^3/kg
x2=V2/Vg2; //Dryness Fraction
H2=Hf2+(x2*Hfg2);
U1=H1-(P1*100*V1);
U2=H2-(P2*100*V2);
U=U2-U1;
printf('Change in Internal Energy: %3.2f kJ/kg',U);
printf('\n');
Q=U+W; //From First law of Thermodynamics
printf('Heat Transferred: %3.2f kJ/kg',Q);
printf('\n');
S1=Sg;
//At 0.5 bar
Sf2=1.091; //in kJ/kg K
Sfg2=6.503; //in kJ/kg K
Sg2=7.594; //in kJ/kg K
S2=Sf2+(x2*Sfg2);
S=S2-S1;
printf('Change in Entropy: %3.2f kJ/kg K',S);
printf('\n');
|
