1 files changed, 40 insertions, 0 deletions
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)