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
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
|
printf("\t example 13.1 \n");
// at atmospheric pressure,Pt=760 mm Hg
printf("\t approximate values are mentioned in the book \n");
x(1)=0.077; // mole fraction of C4
x(2)=0.613; // mole fraction of C5
x(3)=0.310; // mole fraction of C6
printf("\t for T 100 F \n");
Pp(1)=3170; // vapour pressure of C4, from fig 13.3
Pp(2)=790; // vapour pressure of C5,from fig 13.3
Pp(3)=250; // vapour pressure of C6,from fig 13.3
i=1;
while(i<4)
p(i)=(Pp(i)*x(i));
printf(" \n x(i) Pp(i) p(i) \n "+string(x(i))+" "+string(Pp(i))+" "+string(p(i))+" \n");
i=i+1;
end
pt=p(1)+p(2)+p(3);
printf("\t total pressure is : %.1f mm Hg \n",pt);
printf("\t pressure is too high \n");
printf("\t for T 96 F \n");
Pp(1)=2990; // vapour pressure of C4, from fig 13.3
Pp(2)=725; // vapour pressure of C5,from fig 13.3
Pp(3)=229; // vapour pressure of C6,from fig 13.3
i=1;
while(i<4)
p(i)=(Pp(i)*x(i));
printf(" \n x(i) Pp(i) p(i) \n "+string(x(i))+" "+string(Pp(i))+" "+string(p(i))+" \n");
i=i+1;
end
pt=p(1)+p(2)+p(3);
printf("\t total pressure is : %.1f mm Hg \n",pt);
printf("\t pressure is too low \n");
printf("\t for T 97 F \n");
Pp(1)=3040; // vapour pressure of C4, from fig 13.3
Pp(2)=740; // vapour pressure of C5,from fig 13.3
Pp(3)=234; // vapour pressure of C6,from fig 13.3
i=1;
while(i<4)
p(i)=(Pp(i)*x(i));
printf(" \n x(i) Pp(i) p(i) \n "+string(x(i))+" "+string(Pp(i))+" "+string(p(i))+" \n");
i=i+1;
end
pt=p(1)+p(2)+p(3);
printf("\t total pressure is : %.1f mm Hg \n",pt);
i=1;
while(i<4)
y(i)=(Pp(i)*x(i)/pt);
printf("\n x(i) y(i) \n "+string(x(i))+" "+string(y(i))+" \n");
i=i+1;
end
printf("\t solution for b \n");
// Similarly at what temperature will the mixture start to boil if the system is under a pressure of 35 psia
printf("\t for T 150 F \n");
Pp(1)=6100; // vapour pressure of C4, from fig 13.3
Pp(2)=1880; // vapour pressure of C5,from fig 13.3
Pp(3)=680; // vapour pressure of C6,from fig 13.3
i=1;
while(i<4)
p(i)=(Pp(i)*x(i));
printf(" \n x(i) Pp(i) p(i) \n "+string(x(i))+" "+string(Pp(i))+" "+string(p(i))+" \n");
i=i+1;
end
pt=p(1)+p(2)+p(3);
printf("\t total pressure is : %.0f mm Hg \n",pt);
printf("\t pressure is too high \n");
printf("\t for T 149F \n");
Pp(1)=6050; // vapour pressure of C4, from fig 13.3
Pp(2)=1850; // vapour pressure of C5,from fig 13.3
Pp(3)=670; // vapour pressure of C6,from fig 13.3
i=1;
while(i<4)
p(i)=(Pp(i)*x(i));
printf(" \n x(i) Pp(i) p(i) \n "+string(x(i))+" "+string(Pp(i))+" "+string(p(i))+" \n");
i=i+1;
end
pt=p(1)+p(2)+p(3);
printf("\t total pressure is : %.0f mm Hg \n",pt);
i=1;
while(i<4)
y(i)=(Pp(i)*x(i)/pt);
printf("\n x(i) y(i) \n "+string(x(i))+" "+string(y(i))+" \n");
i=i+1;
end
printf("\t solution for c \n");
printf("\t for T 95F \n");
K(1)=3.13; // fig 7
K(2)=0.92; // fig 7
K(3)=0.30; // fig 7
i=1;
while(i<4)
y(i)=(K(i)*x(i));
printf("\n x(i) K(i) y(i) \n "+string(x(i))+" "+string(K(i))+" "+string(y(i))+" \n");
i=i+1;
end
yt=y(1)+y(2)+y(3);
printf("\t yt is : %.3f \n",yt);
printf("\t yt is too low \n");
printf("\t for T 100F \n");
K(1)=3.35; // fig 7
K(2)=1; // fig 7
K(3)=0.335; // fig 7
i=1;
while(i<4)
y(i)=(K(i)*x(i));
printf("\n x(i) K(i) y(i) \n "+string(x(i))+" "+string(K(i))+" "+string(y(i))+" \n");
i=i+1;
end
yt=y(1)+y(2)+y(3);
printf("\t yt is : %.3f \n",yt);
printf("\t yt is too low \n");
printf("\t for T 102F \n");
K(1)=3.45; // fig 7
K(2)=1.02; // fig 7
K(3)=0.35; // fig 7
i=1;
while(i<4)
y(i)=(K(i)*x(i));
printf("\n x(i) K(i) y(i) \n "+string(x(i))+" "+string(K(i))+" "+string(y(i))+" \n");
i=i+1;
end
yt=y(1)+y(2)+y(3);
printf("\t yt is : %.3f \n",yt);
printf("\t solution for d \n");
// The use of K values gives y, directly and permits use of the total mol fraction of yt = 1.00 as the criterion for equilibrium
printf("\t for T 150F \n");
K(1)=2.8; // fig 7
K(2)=1.01; // fig 7
K(3)=0.4; // fig 7
i=1;
while(i<4)
y(i)=(K(i)*x(i));
printf("\n x(i) K(i) y(i) \n "+string(x(i))+" "+string(K(i))+" "+string(y(i))+" \n");
i=i+1;
end
yt=y(1)+y(2)+y(3);
printf("\t yt is : %.3f \n",yt);
printf("\t yt is too low \n");
printf("\t for T 153F \n");
K(1)=2.90; // fig 7
K(2)=1.06; // fig 7
K(3)=0.415; // fig 7
i=1;
while(i<4)
y(i)=(K(i)*x(i));
printf("\n x(i) K(i) y(i) \n "+string(x(i))+" "+string(K(i))+" "+string(y(i))+" \n");
i=i+1;
end
yt=y(1)+y(2)+y(3);
printf("\t yt is : %.3f \n",yt);
printf("\t solution for e at pt=760mm Hg \n");
y(1)=0.077; // mole fraction of C4
y(2)=0.613; // mole fraction of C5
y(3)=0.310; // mole fraction of C6
printf("\t for T 130F \n");
K(1)=5; // fig 7
K(2)=1.65; // fig 7
K(3)=0.62; // fig 7
i=1;
while(i<4)
x(i)=(y(i)/K(i));
printf("\n y(i) K(i) x(i) \n "+string(y(i))+" "+string(K(i))+" "+string(x(i))+" \n");
i=i+1;
end
xt=x(1)+x(2)+x(3);
printf("\t xt is : %.3f \n",xt);
printf("\t xt is too low \n");
printf("\t for T 120F \n");
K(1)=4.4; // fig 7
K(2)=1.4; // fig 7
K(3)=0.51; // fig 7
i=1;
while(i<4)
x(i)=(y(i)/K(i));
printf("\n y(i) K(i) x(i) \n "+string(y(i))+" "+string(K(i))+" "+string(x(i))+" \n");
i=i+1;
end
xt=x(1)+x(2)+x(3);
printf("\t xt is : %.3f \n",xt);
printf("\t xt is high \n");
printf("\t for T 123F \n");
K(1)=4.6; // fig 7
K(2)=1.49; // fig 7
K(3)=0.545; // fig 7
i=1;
while(i<4)
x(i)=(y(i)/K(i));
printf("\n y(i) K(i) x(i) \n "+string(y(i))+" "+string(K(i))+" "+string(x(i))+" \n");
i=i+1;
end
xt=x(1)+x(2)+x(3);
printf("\t xt is : %.3f \n",xt);
printf("\t dew point at 760mm is 123F \n");
printf("\t dew point at 35psia \n");
printf("\t for T 174F \n");
K(1)=3.7; // fig 7
K(2)=1.38; // fig 7
K(3)=0.58; // fig 7
i=1;
while(i<4)
x(i)=(y(i)/K(i));
printf("\n y(i) K(i) x(i) \n "+string(y(i))+" "+string(K(i))+" "+string(x(i))+" \n");
i=i+1;
end
xt=x(1)+x(2)+x(3);
printf("\t xt is : %.3f \n",xt);
printf("\t dew point is 174F \n");
// end
|
