blob: aa55012870fb769f0d948d77175b5a747678e52b (
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
52
53
54
55
56
|
//example 6.3
clc; funcprot(0);
clf()
//exapple 6.3
// Initialization of Variable
Q=[0 0.01 0.02 0.03 0.04 0.05];//discharge
effi_hyd=[65.4 71 71.9 67.7 57.5 39.2];
effi_over=[0 36.1 56.0 61.0 54.1 37.0];
H_sys=[0 0 0 0 0 0]
d=0.114;//diameter of pipe
d_o=0.096;//diameter of impeller
h=8.75;//elevation
g=9.81;//acc. of gravity
rho=999;//denisity of water
l=60;//length of pipe
theta=0.611;//angle in radians
B=0.0125;//width of blades
pi=3.1412
mu=1.109/1000;//viscosity of water
omega=2*pi*1750/60;
// calculation
for i=1:6
if i==1 then
H_sys(i)=h;
else
H_sys(i)=h+8*Q(i)^2/pi^2/d^4/g*(1+8*l*0.0396/d*(4*rho*Q(i)/pi/d/mu)^-0.25);
end,
end;
H_theor=omega^2*d_o^2/g-omega*Q/2/pi/g/B/tan(theta);
//disp(H_sys"head of system (in m)");
//disp(H_theor);
for i=1:6
H_eff(i)=effi_hyd(i)*H_theor(i)/100;
end
//disp(H_eff);
plot(Q,effi_hyd, 'r--d');
plot(Q,effi_over, 'g');
plot(Q,H_eff,'k');
plot(Q,H_theor);
plot(Q,H_sys ,'c-');
title('system characteritics');
ylabel('Head(m)or Efficiency(%)');
xlabel('volumetric flow rate(m^3/s)');
//calculation of power
//at intersecting point using datatrip b/w H_sys &H_eff
Q=0.0336
effi_over=59.9
H_eff=13.10
P=H_eff*rho*g*Q/effi_over/10;
disp(P ,"Power required to pump fluid at this rate(in KW):")
|
