initial commit / add all books

6 files changed, 174 insertions, 0 deletions

diff --git a/3814/CH5/EX5.1/EX5_1.sce b/3814/CH5/EX5.1/EX5_1.sce
new file mode 100644
index 000000000..e8a2fcfd5
--- /dev/null
+++ b/3814/CH5/EX5.1/EX5_1.sce
@@ -0,0 +1,25 @@
+

+// to find efficiency

+// ex 5.1 pgno.115

+clc

+p=67.5*1000 // 67.5 kw to develop wheel

+no=0.83 // efficiency of installation 

+gamma1=9800 // constant gamma

+g=9.8 //gravitational acceleration 

+N=400 // rotates

+H=60 // head of water 60 m

+Q=p/(no*gamma1*H)// volume flow rate

+printf(" Q= %.3f m3/s",Q)

+v1=sqrt(2*g*H) // velocity of the jet

+printf("\n  V1 = %f m/s",v1)

+d=sqrt((0.138*4)/(3.14*v1))

+printf("\n %e m",d)

+r=0.46 //  ratio of bucket speed in rev/min

+u=v1*r //velocity

+printf("\n %f m/s",u)

+D=(H*u)/(%pi*N)

+printf("\n %f m",D)

+w=(2*N*%pi)/(H) //specific speed of trubine

+mprintf('\n specific speed of turbine %f rad/s',w)

+wt=(w*((p/1000)^(0.5)))/((g*H)^((5)/(4)))

+mprintf('\n  wt  = %f',wt)

diff --git a/3814/CH5/EX5.2/Ex5_2.sce b/3814/CH5/EX5.2/Ex5_2.sce
new file mode 100644
index 000000000..4b7bd4cea
--- /dev/null
+++ b/3814/CH5/EX5.2/Ex5_2.sce
@@ -0,0 +1,21 @@
+// to find flow rate and shaft power develpoed by the turbine

+// ex 5.2 pgno 116

+clc

+g=9.8 // gravitional acceleration

+H=400 // head

+hf=23.6 // penstock and nozzle

+d=80e-3 // diameter of the jet

+u=40 // bucket speed

+k=.85 // ratio of heat

+deg=165 // degree

+n1=0.9 // rotational speed

+V1=sqrt(2*g*(H-hf)) // velocity of jet

+mprintf('\n  velocity of jet v1= %f m/s',V1)

+E=u/g*((V1-u)*(1-(k*cosd(deg)))) // eulers head

+mprintf('\n eulers head E = %f m',E)

+Q=(%pi/4)*d^2*V1 // flow rate

+mprintf('\n Flow Rate Q = %f m3/s',Q)

+Pe=g*Q*E // power developed by the runner

+mprintf(' \n power developed by the runner =Pe= %f kw',Pe)

+P=Pe*n1 // nint

+mprintf('\n nint = %f kw',P)

diff --git a/3814/CH5/EX5.3/Ex5_3.sce b/3814/CH5/EX5.3/Ex5_3.sce
new file mode 100644
index 000000000..c105f1805
--- /dev/null
+++ b/3814/CH5/EX5.3/Ex5_3.sce
@@ -0,0 +1,35 @@
+

+// a pelton wheel hydraulic efficiency and over all efficiency

+// ex 5.3 pgno. 117

+clc

+D=1.45 // diameter of the wheel

+N=375 // shaft running

+u=(%pi*D*N)/60 // peripheral velocity

+k=0.9 // coefficient of the bucket

+p=3750 //peripherial velocity

+hf=200*0.1 // head availabe

+mprintf('\n  peripherial velocity u =%f m/s',u)

+mprintf('\n Total Head = %d m',hf)

+

+h1=200  // total head

+l=20 // losses

+H=h1-l  //effective head

+g=9.8 // gravity

+mprintf('\n effective head H = %d m',H)

+V1=sqrt(2*g*H) // velocity of the jet

+mprintf('\n velocity of the jet V1= %f m/s',V1)

+S=u/V1 // speed ratio

+mprintf('\n Speed Ratio =u/V1= %f',S)

+nh=2*((S)*(1-S)*(1-k*cosd(165))) // hydraulic efficiency 

+mprintf('\n  Hydraulic efficiency nh= %f percentage',(nh*100))

+E=(u/g)*(V1-u)*(1-(k*cosd(165))) // euler's head

+mprintf('\n E =%f m',E)

+no=k*nh // realation between

+mprintf('\n Relation between n0= %f',no)

+hp=p/no // hydraulic power

+mprintf('\n hydraulic power = %d kw',hp)

+gamma1=9800 // constant gamma

+Q=(1000*hp)/(2*gamma1*H) // flow rate

+mprintf('\n Flow rate Q = %f m3/s',Q)

+d=sqrt((4*Q)/(%pi*V1)) // diameter 

+mprintf('\n d = %f m',d)

diff --git a/3814/CH5/EX5.4/Ex5_5.sce b/3814/CH5/EX5.4/Ex5_5.sce
new file mode 100644
index 000000000..8e4cd5fd7
--- /dev/null
+++ b/3814/CH5/EX5.4/Ex5_5.sce
@@ -0,0 +1,26 @@
+// to find coefficient of velocity speed ratio,jet diameter

+//ex 5.5 pgno119

+clc

+cv=0.98//velocity of volume

+g=9.8//gravity

+h=130//head loss

+V1=cv*(sqrt(2*g*h))//velocity of jet 

+mprintf('\n velocity of the jet = %f m/s',V1)

+s=0.46//specific gravity

+u=s*V1//velocity 

+mprintf('\n u = %f m/s',u)

+N=200//Rotational speed

+D=(60*u)/(%pi*N)//Diameter

+mprintf('\n  peripherial velocity u =%f m',D)

+d=D/9

+mprintf('\n d =%f m',d)

+p=8e6//petlon turbine 

+no=0.87//eficiency

+gamma1=9800//constant gamma

+Q=(p/(no*gamma1*h))//volume flow rate

+mprintf('\n Q = %f m3/s',Q)

+n=(Q*4)/(%pi*d*d*V1)//

+mprintf('\n number of jets n =%d',n)

+Z=(D/(2*d))+15

+mprintf('\n number of buckets %f ',Z)

+

diff --git a/3814/CH5/EX5.6/Ex5_6.sce b/3814/CH5/EX5.6/Ex5_6.sce
new file mode 100644
index 000000000..c3d291b2f
--- /dev/null
+++ b/3814/CH5/EX5.6/Ex5_6.sce
@@ -0,0 +1,39 @@
+// to find pelton turbine completely 

+// ex 5.6 pgno.120

+clc

+P=100/4 //power each unit

+mprintf('\n power output of each unit P = %d MW',P)

+gammma=9800 //constant gammma

+Q=6.85 //flow rate

+H=580 //head

+g1=9.8

+N=428 // speed

+t1=60 // temperature

+n=2 // types of turbine

+k=0.95 //ratio of head

+Hp=(gammma*Q*H)/(1000*1000) // hydraulic efficiency

+mprintf('\n hydraclic power = %f MW',Hp)

+on=P/Hp // overall efficiency

+mprintf('\n Overall efficiency = %f',on)

+sp=0.46 // assuming speed ratio

+V1=sqrt(2*g1*H) // velocity of jet

+mprintf('\n velocity of the jet V1 =%f m/s',V1)

+u=V1*sp // peripherial velocity

+mprintf('\n u =%f m/s',u)

+D=(t1*u)/(%pi*N)// peripherial velocity

+mprintf('\n  peripherial velocity %f m',D)

+d=sqrt(((Q)/((%pi/4)*V1*n)))

+mprintf('\n %f m',d)

+Z=((D)/(2*d))+15 //number of buckets

+mprintf('\n number of bukets Z =%f m',Z)

+m=D/d // jet ratio

+mprintf('\n jet ratio = m= %f',m)

+L=2.5*d // length

+mprintf('\n Radial length of bucket L = %f m',L)

+B=4*d // width

+mprintf('\n width of bucket B =%f m',B)

+mprintf('\n Depth of bucket hyrauclic efficiency %f m',d)

+nb=2*(u/V1)*(1-(u/V1))*(1-cosd(160))

+mprintf('\n %f',nb)

+nm=on/nb

+mprintf('\n nm =%f',nm)

diff --git a/3814/CH5/EX5.7/Ex5_7.sce b/3814/CH5/EX5.7/Ex5_7.sce
new file mode 100644
index 000000000..57da720ce
--- /dev/null
+++ b/3814/CH5/EX5.7/Ex5_7.sce
@@ -0,0 +1,28 @@
+

+// to find flow rate wheel diameter of each jet

+// ex 5.7 pgno.121

+clc

+p=4.5e6 //pelton wheel develop 

+no=0.8 // wheel diameter

+g=9800 //gravitional acceleration

+h=120 //head loss

+g1=9.8 

+p1=1000 // over all efficiency

+N=200 //  rotational speed

+Q=p/(no*g*h) // flow rate 

+mprintf('\n overall efficiency no= %f m3/s',Q)

+v1=sqrt(2*9.8*h) // velocity of the jet

+mprintf('\n velocity of jet =%f m/s',v1)

+u=v1*0.42 // peripherial velocity

+mprintf('\n u =%f m/s',u)

+N=200 // speed

+D=(60*u)/(%pi*N) // diameter of the jet

+mprintf('\n D =%f m',D)

+d=D/8 // to find diameter

+mprintf('\n diameter of jet = %f meter',d) // to display diameter

+n=(Q*4)/(%pi*d*d*v1) // to calculate jets

+mprintf('\n number of jets n = %d',n) // to display number of jets

+w=(2*N*%pi)/60 // to calculate speed

+wt=w*(((p/p1)^0.5)/((g1*h)^(5/4))) // specific speed

+mprintf('\n wt =%f RPM',wt)

+