From b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b Mon Sep 17 00:00:00 2001
From: priyanka
Date: Wed, 24 Jun 2015 15:03:17 +0530
Subject: initial commit / add all books

---
 1457/CH15/EX15.1/15_1.sce    | 22 ++++++++++++++++++++++
 1457/CH15/EX15.2/15_2.sce    | 35 +++++++++++++++++++++++++++++++++++
 1457/CH15/EX15.3.a/15_3.sce  | 18 ++++++++++++++++++
 1457/CH15/EX15.3.b/15_3b.sce | 22 ++++++++++++++++++++++
 4 files changed, 97 insertions(+)
 create mode 100755 1457/CH15/EX15.1/15_1.sce
 create mode 100755 1457/CH15/EX15.2/15_2.sce
 create mode 100755 1457/CH15/EX15.3.a/15_3.sce
 create mode 100755 1457/CH15/EX15.3.b/15_3b.sce

(limited to '1457/CH15')

diff --git a/1457/CH15/EX15.1/15_1.sce b/1457/CH15/EX15.1/15_1.sce
new file mode 100755
index 000000000..8617221d9
--- /dev/null
+++ b/1457/CH15/EX15.1/15_1.sce
@@ -0,0 +1,22 @@
+clc
+//Initialization of variables
+z2=500//ft
+z1=300//ft
+D=[1 1.5 2 2.5 3 4 6]
+g=32.2
+gam=62.4
+//calculations
+Dj=D/12
+Vj=sqrt((z2-z1)*2*g./(1.04 + 640.*Dj.^4))
+Aj=%pi/4 *Dj.^2
+Q=Aj.*Vj
+Pjet=gam*Q.*Vj.^2 /(2*g) /550
+Pj=max(Pjet)
+for i=1:length(Pjet)
+    if(Pjet(i) ==Pj)
+        break
+    end
+end
+diameter=D(i)
+//results
+printf("Thus a pipe of %d in will be the optimum",diameter)
diff --git a/1457/CH15/EX15.2/15_2.sce b/1457/CH15/EX15.2/15_2.sce
new file mode 100755
index 000000000..4cb75d10a
--- /dev/null
+++ b/1457/CH15/EX15.2/15_2.sce
@@ -0,0 +1,35 @@
+clc
+//Initialization of variables
+phi=0.46
+g=32.2
+k=0.44
+cv=0.98
+d=10 //in
+A=0.545 //ft^2
+beta=160 //degrees
+//calculations
+u=phi*sqrt(2*g)
+V1=cv*sqrt(2*g)
+gQ=62.4*A*V1
+T=d/2 *gQ/g *(1 -cosd(beta) /sqrt(1+k) )*sqrt(2*g)*(cv-phi)
+Power=T*2*u/d
+//results
+printf("Torque required = %d ft lb",T)
+printf("\n Power transferred = %d ft lb/s",Power)
+Pi=gQ
+He=Power/Pi
+printf("\n Hydraulic efficiency = %.2f",He)
+v1=V1-u
+v2=v1/(sqrt(1+k))
+hl=k*v2^2 /(2*g)
+printf("Head loss in bucket friction = %.4f",hl*100)
+Hn=(1/cv^2 -1)*V1^2 /(2*g)
+printf("\n Head loss in  nozzle = %.4f",Hn*100)
+V2cos=u+v2*cosd(beta)
+V2sin=v2*sind(beta)
+alpha=acot(V2cos/V2sin)
+V2=V2sin/sin(alpha)
+Hd=V2^2/(2*g)
+printf("\n Head loss at discharge = %.2f",Hd*100)
+Htotal=Hd+Hn+hl
+printf("\n Total head loss = %.2f",Htotal*100)
diff --git a/1457/CH15/EX15.3.a/15_3.sce b/1457/CH15/EX15.3.a/15_3.sce
new file mode 100755
index 000000000..275e081f5
--- /dev/null
+++ b/1457/CH15/EX15.3.a/15_3.sce
@@ -0,0 +1,18 @@
+clc
+//Initialization of variables
+cv=0.98
+g=32.2
+h=1320 //ft
+A=0.196 //ft^2
+eta=0.85
+ne=400
+phi=0.45
+//calculations
+V=cv*sqrt(2*g*h)
+Q=A*V
+bhp=eta*62.4*Q*h/550
+ns=ne*sqrt(bhp) /h^(5/4)
+u=phi*sqrt(2*g*h)
+D=u*60/%pi/ne
+//results
+printf("Pitch diameter = %.2f ft",D)
diff --git a/1457/CH15/EX15.3.b/15_3b.sce b/1457/CH15/EX15.3.b/15_3b.sce
new file mode 100755
index 000000000..4f0e1d41f
--- /dev/null
+++ b/1457/CH15/EX15.3.b/15_3b.sce
@@ -0,0 +1,22 @@
+clc
+//Initialization of variables
+cv=0.98
+g=32.2
+h=1320 //ft
+A=0.196 //ft^2
+eta=0.85
+ne=400
+phi=0.45
+//calculations
+V=cv*sqrt(2*g*h)
+Q=A*V/3
+bhp=eta*62.4*Q*h/550
+ne2=600
+ns1=ne2*sqrt(bhp) /h^(5/4)
+D=2500/ne2
+Dj=sqrt(Q*4/V/%pi)
+//results
+printf("Jet diameter = %.3f ft",Dj)
+printf("\n Specific speed = %.2f ",ns1)
+printf("\n Pitch Diameter = %.2f ft",D)
+printf("\n Operating speed = %d rpm",ne2)
-- 
cgit