blob: 7b390666b9dd63424720dc4392b1ce9ef9f0726e (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
|
clc;
warning("off");
printf("\n\n example14.1 - pg726");
// given
T=40+273.15; //[K] - temperature
P=1; //[atm] - pressure
sigma=3.711*10^-10; //[m]
etadivkb=78.6; //[K]
A=1.16145;
B=0.14874;
C=0.52487;
D=0.77320;
E=2.16178;
F=2.43787;
Tstar=T/(etadivkb);
// using the formula si=(A/(Tstar^B))+(C/exp(D*Tstar))+(E/exp(F*Tstar)
si=(A/(Tstar^B))+(C/exp(D*Tstar))+(E/exp(F*Tstar));
M=28.966; //[kg/mole] - molecular weight
// using the formula mu=(2.6693*(10^-26))*(((M*T)^(1/2))/((sigma^2)*si))
mu=(2.6693*(10^-26))*(((M*T)^(1/2))/((sigma^2)*si));
printf("\n\n The viscosity of air is \n mu=%eNs/m^2=%fcp",mu,mu*10^3);
|
