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
|
clc;
clear all;
disp("composition determination")
T=273+15;// K
G=8314;// kJ/(kgmole*K)
p=1.1*10^(5);// N/m^2
pO2=0.21*p;
pN2=0.79*p;
CO2=pO2/(G*T);// kg mole/m^3
disp("kg mole/m^3",CO2,"molar concentration of O2, CO2 =")
CN2=pN2/(G*T);// kg mole/m^3
disp("kg mole/m^3",CN2,"molar concentration of N2, CN2 =")
MO2=32;
MN2=28;
rhoO2=MO2*CO2;//kg/m^3
rhoN2=MN2*CN2;//kg/m^3
disp("kg/m^3",rhoO2,"molar density of O2 rhoO2 =")
disp("kg/m^3",rhoN2,"molar density of N2, rhoN2 =")
rho=rhoO2+rhoN2;
xO2=rhoO2/rho;
disp(xO2,"mass fraction of O2 xO2 =")
xN2=rhoN2/rho;
disp(xN2,"mass fraction of N2, xN2 =")
C=CO2+CN2
mO2=CO2/C;
disp(xO2,"mole fraction of O2, mO2 =")
mN2=CN2/C;
disp(xN2,"mole fraction of N2, mN2 =")
|
