blob: 00dac5414b228933fdb720fd18a3665206b930b5 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
|
clc;
clear;
D=2;//ft
Q=30;//cfs
Dm=3;//in
//Rem=Re; hence (Vm*Dm/kvism)=(V*D/kvis); where kvis is kinematic viscosity
//kvis=kvism; same fluid is used for model and prototype
//(Vm/V)=(D/Dm)
//Q=VA; hence Qm/Q = (Vm*Am)/(V*A)=(Dm/D)
Qm=(Dm/12)*Q/D;//cfs
disp("cfs",Qm,"The required flowrate in the model=")
Drat=0.04:0.01:1;
count=1;
for i=0.04:0.01:1
Vrat(count)=1/i;
count=count+1;
end
plot2d(Drat,Vrat,rect=[0,0,1,25])
xtitle("Vm/V vs Dm/D","Dm/D","Vm/V")
|
