7 files changed, 124 insertions, 0 deletions
diff --git a/3814/CH8/EX8.1/Ex8_1.sce b/3814/CH8/EX8.1/Ex8_1.sce
new file mode 100644
index 000000000..5adbd33ce
--- /dev/null
+++ b/3814/CH8/EX8.1/Ex8_1.sce
@@ -0,0 +1,15 @@
+// select pump deliver 1890 l/min 
+clc
+Q=(1890e-3/60)
+disp(Q)
+p=448e3
+N=3600 //rev/min
+w=(2*N*%pi)/60
+g=9.8
+gammma=9800
+mprintf('\n speed in rad/s = w= %f rad/s',w)
+H=p/gammma
+mprintf('\n head in meters H = %f m',H)
+wp=((w*(sqrt(Q)))/((g*H)^(3/4)))
+mprintf('\n specific speed of the pump giveb by Wp=%f',wp)
+mprintf('Wp<1 therefore radial pump selected')
diff --git a/3814/CH8/EX8.2/Ex8_2.sce b/3814/CH8/EX8.2/Ex8_2.sce
new file mode 100644
index 000000000..1a91bc644
--- /dev/null
+++ b/3814/CH8/EX8.2/Ex8_2.sce
@@ -0,0 +1,9 @@
+// elevation that the diameter pump can be situated above the water surface of suction
+clc
+patm=101e3
+pv=1666
+g=9800
+npsh=7.4
+z1=((patm-pv)/g)-npsh
+mprintf('\n Z1= %f m',z1)
+mprintf('\n the pump must be place at approximately %f m above the suction reservoir of water surface',z1)
diff --git a/3814/CH8/EX8.3/Ex8_3.sce b/3814/CH8/EX8.3/Ex8_3.sce
new file mode 100644
index 000000000..a51ce2e1e
--- /dev/null
+++ b/3814/CH8/EX8.3/Ex8_3.sce
@@ -0,0 +1,10 @@
+// radial flow pump characteristic is given by
+clc
+a=196
+b=-10.7
+c=7.9
+a1=(1/(2*a))
+Q=a1*(-b+sqrt((b^2)+(4*a*c)))
+mprintf('\n Operating point at  Q = %f m3/s',Q)
+H=15+(85*Q^2)
+mprintf('\n H = %f m',H)
diff --git a/3814/CH8/EX8.4/Ex8_4.sce b/3814/CH8/EX8.4/Ex8_4.sce
new file mode 100644
index 000000000..3290b62fe
--- /dev/null
+++ b/3814/CH8/EX8.4/Ex8_4.sce
@@ -0,0 +1,40 @@
+// determine flow rate,theoretical head required power,pressure across the impeller
+clc
+b1=44
+r1=21e-3 // mm
+B=11e-3 // mm
+r2=66e-3 //mm
+b2=5e-3//mm
+N=2500  //rpm
+h1=0
+g=9.8
+alpha=90 // degree
+D1=2
+D2=2
+u1=(2*%pi*N*r1)/60
+gamm1=9800
+p1=1000
+mprintf('\n peripherial velocity at inlet u1=wR1 =%f m/s',u1)
+u2=(2*%pi*N*r2)/60
+mprintf('\n peripherial velocity at exit  u2= wR2=%f m/s',u2)
+V1=tand(b1)*u1
+mprintf('\n V1f = %f m/s',V1)
+Q=%pi*2*r1*B*V1
+mprintf('\n Q = %f m3/s',Q)
+V2f=Q/(2*%pi*r2*b2)
+mprintf('\n V2f =%f m/s',V2f)
+V2w=V2f/(tand(30))
+mprintf('\n u2-V2w = %f ',V2w)
+v2w=u2-V2w
+mprintf('\n V2w = %f m/s',v2w)
+alpha2=atand(V2f/v2w)
+mprintf('\n  alpha2 = %f degree',alpha2)
+v2=v2w/cosd(18.9)
+mprintf('\n  V2= %f m/s',v2)
+H1=(u2*v2w)/g
+mprintf('\n H1 = %f m',H1) 
+p=gamm1*Q*H1
+mprintf('\n H1 = %f watt',p)
+P2=(p1*g*H1)-((p1/2)*(v2^2-V1^2))
+mprintf('\n p2-p1 = %e Pa',P2)
+mprintf('\n p2-p1 = %f bar',P2/10^5)
diff --git a/3814/CH8/EX8.5/Ex8_5.sce b/3814/CH8/EX8.5/Ex8_5.sce
new file mode 100644
index 000000000..20eaab5ab
--- /dev/null
+++ b/3814/CH8/EX8.5/Ex8_5.sce
@@ -0,0 +1,17 @@
+// determine  discharge and head with two identail
+clc
+f=0.025
+l=70
+D=0.3
+k=2.5
+g=9.8
+m=((f*l/D)+k)/(2*g*(((%pi*D*D)/4)^2))
+disp(m)
+mprintf('\n H1 =15 +%d Q^2',m)
+b=5.35
+a=112.8
+c=7.9
+Q=(1/(2*a))*(b+sqrt((b^2)+(4*a*c)))
+mprintf('\n Q= %f m3/s',Q)
+H1=15+85*Q^2
+mprintf('\n H1 = %f m',H1)
diff --git a/3814/CH8/EX8.6/Ex8_6.sce b/3814/CH8/EX8.6/Ex8_6.sce
new file mode 100644
index 000000000..3e265430b
--- /dev/null
+++ b/3814/CH8/EX8.6/Ex8_6.sce
@@ -0,0 +1,16 @@
+// Determine the velocity of flow theoretical head and power required to drive the pump
+clc
+Q=150e-3
+d1=300e-3
+d2=150e-3
+n=500
+g=9.8
+Vf=(Q)/((%pi/4)*(d1^2-d2^2))
+mprintf('\n Velocity of flow Vf=Q/A %f m/s',Vf)
+D=(d1+d2)/2
+mprintf('\n Peripherial velocity is calculated on the mean diameter D =%f m',D)
+u=((2*%pi*n)/60)*(D/2)
+H1=(u^2/g)-((u*Vf)/g)*(cotd(75)+cotd(70))
+mprintf('\n Theoretical Head H = %f m',H1)
+P=g*Q*H1
+mprintf('\n Required power P = %f kw',P)
diff --git a/3814/CH8/EX8.7/Ex8_7.sce b/3814/CH8/EX8.7/Ex8_7.sce
new file mode 100644
index 000000000..34947701c
--- /dev/null
+++ b/3814/CH8/EX8.7/Ex8_7.sce
@@ -0,0 +1,17 @@
+// maximum height that the pump
+clc
+Q=3.5e-3
+nphs=4.5
+t=15
+d=0.1
+patm=101e3
+g=9.8
+k=20
+V=(4*Q)/(%pi*d^2)
+mprintf('\n Velocity in the suction pipe V =Q/A  %f m/s',V)
+h1=(k*V^2)/(2*g)
+mprintf('\n h1= %f m',h1)
+pv=1666
+z1=(patm/9800)-(h1)-(pv/9800)-nphs
+mprintf('\n z1 = %f m',z1)
+mprintf('\n the pump should be located higher than %f m above the water surface ',z1)