2 files changed, 34 insertions, 0 deletions
diff --git a/3814/CH9/EX9.1/EX9_1.sce b/3814/CH9/EX9.1/EX9_1.sce
new file mode 100644
index 000000000..7c508cb33
--- /dev/null
+++ b/3814/CH9/EX9.1/EX9_1.sce
@@ -0,0 +1,16 @@
+//  to find hydralic power 
+// ex 9.1 pgno.215
+clc;
+p= 200e3 // pressure of water
+g=9800 
+zs=0.1
+pd=600e3
+sp=0.85 // efficiency of the pump
+h=((p/g)+zs)+((pd/g)+zs)
+printf(" %f  m",h)
+q=0.2 // increase the pressure
+h1=81.6
+hp=g*q*h1 // efficiency 
+printf("\n %e W",hp)
+n=hp/sp // electrical power i e shaft power
+printf("\n electrical power = %5e  W",n)
diff --git a/3814/CH9/EX9.2/EX9_2.sce b/3814/CH9/EX9.2/EX9_2.sce
new file mode 100644
index 000000000..19e232698
--- /dev/null
+++ b/3814/CH9/EX9.2/EX9_2.sce
@@ -0,0 +1,18 @@
+//determine the total head of the pump
+// ex 9.2  page no 215
+clc;
+q=37.5e-3*4 // water flow rate
+A=3.14*0.15^2// area  of suction in 15cm in meter
+V=q/A
+printf("velocity at the suction Q/As=  %2.2f m/s",V)
+Ad=3.14*0.125^2 // area of suction in 12.5cm
+Vd=q/Ad
+printf("\n  velocity at the discharge side Q/Ad=  %2.2f m/s",Vd)
+ps=54e3
+gamma1=9800// constant gamma
+g=9.8
+vs=2 // velocity of suction
+pd=160e3 // power density
+vd=3 // velocity of discharge side
+H=((ps/gamma1)+(vs^2/(2*g)))+((pd/gamma1)+(vd^2/(2*g)))
+disp(H)