blob: a16575690759c78ff789d8c16dbc0331d2ac6380 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
|
clear ;
clc;
printf("\t Example 11.3\n");
T1=276; //temp.of air,K
aa=3.711; //lennard jones constant or collision diameter,A
ab=2.827; //lennard jones constant or collision diameter,A
b1=78.6; //lennard jones constant,K
b2=59.7; //lennard jones constant,K
a=(aa+ab)/2; //effective molecular diameter for collisions of hydrogen and air,m
b=(b1*b2)^0.5; //effective potential well depth,K
c=T1/b;
d=0.8822; //potential well function
Dab=1.8583*10^-7*T1^1.5/(a^2*d)*(1/2.016+1/28.97)^0.5; //diffusion coefficient of hydrogen in air,m^2/s
printf("\t diffusion coefficient of hydrogen in air is :% -5e m^2/s an experimental value from table is 6.34*10^-5 m^2/s,so the prediction is high by 5 percent.\n",Dab);
//end
|
