blob: 5a3dbd287bb7140eae2cecf154feb70dd5ca50bd (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
|
//clc()
P = 101.3;//kPa
T = [283 293 303 313];
for i=1:4
Ps(i) = exp(13.8858 - 2788.51/(T(i)-52.36));
end
for j =1:4
Ys(j) = Ps(j) * 78.048 / ((P - Ps(j))*29);
end
disp(Ps)
disp(Ys)
plot(T,Ys,rect = [270,0,323,0.9]);
//Tas = Tg - L *(Y1as - Y1) / Cs
//Cs = Cb + Y1 *Ca = 1.005 + Y1 * 1.2,
L = 435.4;//kJ/kgK
//for different value of Tg and Y1 tried, we have the following set of values
Tg = [283 290.4 300 310.1 320.8];
Y1 = [0.1701150 0.15 0.125 0.1 0.075];
plot(Tg,Y1);
xtitle('Saturation curve and adiabatic cooling line','Temperature, K','Y, kg benzene/kg dry air');
|
