summaryrefslogtreecommitdiff
path: root/1379/CH9/EX9.1.7/example9_7.sce
blob: 50f11f30f9fe50c1eb2c2a6a0fa6b9f68e5033d6 (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


//exapple 9.7 
clc; funcprot(0);
// Initialization of Variable
rhog=1200;//density of glycerol
mu=1.45;//viscosity of glycerol
pi=3.1428;
g=9.81;
rhos=2280;//density of sphere
d=8/1000;
s=0;
uf=0.8*0.026;
//calculation
function[a]=intre()
    u=linspace(0,uf,1000);
    for i=1:1000
        y=((pi/6*d^3*rhos*g-pi*d^3/6*rhog*g-0.5*pi*d^2/4*24*mu/d/rhog*rhog*u(i))/pi*6/d^3/rhos)^(-1)*uf/1000;
        s=s+y;
    end
    a=s;
endfunction
[t]=intre();
disp(t,"Time taken by particle to reach 80% of its velocity in (s):");