blob: 54d5d5b06198ee629451a5522c71b42a5d0a4c0b (
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
|
clc
V=0.6; //m^3
p1=12*10^5; //Pa
p2=18*10^5; //Pa
T=298; //K
R0=8.314;
x_O2=0.23;
x_N2=0.77;
n=p1*V/R0/10^3/T;
//Considering 100 kg of air
m_O2=23; //kg
m_N2=77; //kg
M_O2=32;
M_N2=28;
m=100; //kg
R=(m_O2/M_O2 + m_N2/M_N2)*R0/m; //for air
M=R0/R; //for air
m=p1*V/R/T/10^3;
m_O2=x_O2*m;
disp("Mass of O2=")
disp(m_O2)
disp("kg")
m_N2=x_N2*m;
disp("Mass of N2=")
disp(m_N2)
disp("kg")
//After adding CO2 in the vessel
p2=18*10^5; //Pa;
// p_CO2+p_N2+p_O2=18*10^5
// p_N2 + p_O2=12*10^5
p_CO2=6*10^5; //Pa
M_CO2=44;
R_CO2=R0/M_CO2;
m_CO2=p_CO2*V/(R_CO2*10^3*T);
disp("Mass of CO2 = ")
disp(m_CO2)
disp("kg")
|
