blob: 8df170c9952d6d5259a93ce6d9492f5e10ac501e (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
|
clc;
disp("Example 4.11")
d=0.1 // diameter in m
l=25 // length in m
density=1000 // density in kg/m^3
delP= 14700 // in N/m^2
mew= 0.001 // in kg/ms
ka= d*density*((delP*d)^0.5)/(((2*density*l)^0.5)*mew)
Re= (ka/0.281)^(8/7)
if(Re<50000)
v1=Re*mew/(d*density)
disp(v1)
end
if(Re>50000)
Re=(ka/0.2145)^(10/9)
v1=Re*mew/(d*density)
disp(v1)
end
|
