blob: d44290385b4087388c5d6838e21ab267faebfcc7 (
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
49
50
51
|
clear;
clc;
printf("\n Example 3.4");
v_particle = 6; //velocity of the particle in m/sec
v_water = 1.2; //veloicity of the water in m/sec
v_rel = v_particle - v_water;//relative velocity of particles relative to the fluid in m/sec
Re1 = 6*10^(-3)*v_rel*1000/(1*10^(-3));//Re1 is the reynold's no.
printf("\nReynold no. is %d",Re1);
//When the particle has been retarded to a velocity such that Re=500
ydot = (v_rel*500)/Re1;
printf("\nParticle velocity is %.3f m/sec\n",ydot);
c = 0.33/(6*10^(-3))*(1000/2500);
f = sqrt((3*6*10^(-3)*(2500-1000)*9.81)/1000);
function[y]=Fa(t)
y = (-1/22)*(log(cos(0.517*22*t) + 4.8/0.517*sin(0.517*22*t)));
funcprot(0);
endfunction
function[yd]=deriv(t)
yd= -0.083+(0.517*(9.28*cos(11.37*t) - sin(11.37*t))/(cos(11.37*t) + 9.28*sin(11.37*t)));
funcprot(0);
endfunction
function[ydd]=double_deriv(t)
ydd = -0.517*(11.37)^2*(9.28*cos(11.37*t) - sin(11.37*t))/(cos(11.37*t) + 9.28*sin(11.37*t));
funcprot(0);
endfunction
told = 0;
while 1
tnew = told - deriv(told)/double_deriv(told);
if (tnew == told) then
y = Fa(told);
d = y;
printf("\nThe distance moved with speed less than 0.083m/sec is %.3fm",d);
t=told;
printf("\n The time taken by particle to move this distance is %.3fsec",t);
break;
end
told = tnew;
end
printf("\nThe distance moved by the particle relative to the walls of the plant %.3fm",1.2*t - d);
|
