blob: 32ec2ccaee1e163d1daa8fbb0865c85f0dd4a837 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
|
clc;
warning("off");
printf("\n\n example12.12 - pg 594");
// given
pp=1.13*10^4; //[kg/m^3] - density of lead particle
p=1.22; //[kg/m^3] - density of air
g=9.80; //[m/sec^2] - acceleration due to gravity
d=2*10^-3; //[m] - diameter of particle
mu=1.81*10^-5; //[kg/m*sec] - viscosity of air
// let us assume
Cd=0.44;
Ut=((4*d*g*(pp-p))/(3*p*Cd))^(1/2);
disp(Ut)
Nre=(Ut*d*p)/(mu);
// from fig 12,16 value of Cd is
Cd=0.4;
Ut=((4*d*g*(pp-p))/(3*p*Cd))^(1/2);
Nre=(Ut*d*p)/(mu);
// Within the readibility of the chart Cd is unchanged and therefore the above obtained Cd is the final answer
printf("\n\n The terminal velocity is \n Ut = %f m/sec",Ut);
|
