8 files changed, 227 insertions, 0 deletions
diff --git a/3814/CH7/EX7.1/Ex7_1.sce b/3814/CH7/EX7.1/Ex7_1.sce
new file mode 100644
index 000000000..082157785
--- /dev/null
+++ b/3814/CH7/EX7.1/Ex7_1.sce
@@ -0,0 +1,27 @@
+// determine inlet and exit angles of blades mean diameter
+// ex 7.1 pgno. 169
+clc
+H=21.8 // head of turbine
+P=21 // MW power
+N=140 // number runs of rpm
+D1=4.5 // diameter in m
+Dh=2.0 // in meter
+nh=0.94 // efficiency
+nn= 0.88// efficiency in exit angles
+g=9.8
+M=(D1+Dh)/2 //mean diameter in m
+mprintf('Mean Diameter  = %f m',M)
+w=(2*N*%pi)/60
+mprintf('\n w= %f rad/s',w)
+u=(w*M)/2
+mprintf('\n u=wr = %f m/s',u)
+Vlw=(nh*g*H)/u
+mprintf('\n Vlw = %f m/s',Vlw)
+Q=(P*1000*1000)/(nn*H*9800) // flow rate
+mprintf('\n Q=%f m3/s',Q)
+vf=(4*Q)/(%pi*((D1^2)-(Dh^2)))
+mprintf('\n Vf = %f m/s',vf)// velocity in m/s
+b1=vf/(u-Vlw) // tan b1
+mprintf(' \n  tanb1 =%f degree ',atand(b1)) // inlet angles
+b2=vf/u // tan b2
+mprintf('\n b2= %f degree',atand(b2))// exit angles
diff --git a/3814/CH7/EX7.2/Ex7_2.sce b/3814/CH7/EX7.2/Ex7_2.sce
new file mode 100644
index 000000000..f0ccdf985
--- /dev/null
+++ b/3814/CH7/EX7.2/Ex7_2.sce
@@ -0,0 +1,33 @@
+// calculate axial velocity ,flow rate ,exit blade angle,eulers power
+// ex 7.2 pgno. 170
+clc
+D1=2 // m
+Dh=0.8 // m
+N=250 //rpm
+alpha1 =42 // degree
+beta1=148// degree
+D=(D1+Dh)/2 // diameter
+g=9.8 
+mprintf('\n D= %f m',D)
+u=(%pi*D*N)/60 // peripherical velocity of blade
+mprintf('\n u =%f m/s',u)
+a=(180-148) //area
+disp(a)
+d=a*18.3 // diameter
+disp(tan(d))
+
+vlw=(tand(42)+tand(32))  
+disp(vlw)
+Vlw=tand(d)/vlw
+disp(Vlw)
+vlw=7.5 
+vf=vlw*tand(alpha1) // inlet trangle of velocities
+mprintf('\n Vf = %f m/s',vf)
+Q=(%pi/4)*((D1^2-Dh^2)*vf) // flow rate
+mprintf('\n Q = %f m3/s',Q)
+E=(u*vlw)/g // euler's head
+mprintf('\n Euler s Head E =%f m',E)
+Ep=(9800*Q*E)/1000 // eulers power
+mprintf('\n Eulers power = %e W',Ep)
+b2=vf/u 
+mprintf('\n b2 = %f Degree',atand(b2))
diff --git a/3814/CH7/EX7.3/Ex7_3.sce b/3814/CH7/EX7.3/Ex7_3.sce
new file mode 100644
index 000000000..8e643066e
--- /dev/null
+++ b/3814/CH7/EX7.3/Ex7_3.sce
@@ -0,0 +1,68 @@
+// determine blade angles at hub,mean and tip diameters
+clc
+Dt=4.5 // meter
+Dh=2 // meter
+p=20e6 //watts
+N=150 // rpm
+nh=0.94 // hydraulic efficiency 
+n0=0.88 // overall efficiency
+h=21 // head
+g=9.8
+Q=(p)/(n0*h*9800) // 
+mprintf('\n Q =%f m3/s',Q)
+vf=(4*Q)/(%pi*(Dt^2-Dh^2)) // velocity of flow
+mprintf('\n Vf = %f m/s',vf)
+Vw=(60*g*h*nh)/(2*N*%pi) // velocity of whirl
+mprintf('\n rVw = %f',Vw)
+D=(Dt+Dh)/2 // diameters
+mprintf('\n D =%f',D)
+rm=D/2 
+mprintf('\n rm =%f',rm)
+Vm1=Vw/rm
+mprintf('\n Vm1= %f m/s',Vm1)
+rt=rm+0.625
+Vm2=Vw/rt
+mprintf('\n Vm2 =%f m/s',Vm2)
+rt=rm-0.625
+Vm3=Vw/rt
+mprintf('\n Vm2 =%f m/s',Vm3)
+uh=(%pi*Dh*N)/60
+mprintf('\n uh=%f m/s',uh)
+um=(%pi*D*N)/60
+mprintf('\n um=%f m/s',um)
+ut=(%pi*Dt*N)/60
+mprintf('\n ut=%f m/s',ut)
+b1h=vf/(uh-Vm3)
+mprintf('     hub       : ')
+mprintf('\n tan(pi-beta1h)=  %f',b1h)
+be=atand(b1h)
+B1h=180-be
+mprintf('\n B1h = %f degree',B1h)
+B2=vf/uh
+mprintf('\n B2h =%f degree ',atand(B2))
+mprintf('     mean       : ')
+b2h=vf/(um-Vm1)
+mprintf('\n tan(pi-beta1h)=  %f',b2h)
+be1=atand(b2h)
+B2h=180-be1
+mprintf('\n B1m = %f degree',B2h)
+B2m=vf/um
+mprintf('\n B1m =%f degree ',atand(B2m))
+mprintf('\    ')
+mprintf('\    ')
+mprintf('\    ')
+mprintf('         Tip:        ')
+b3h=vf/(ut-Vm2)
+mprintf('\n                      tan(pi-beta1h)=  %f',b3h)
+be2=atand(b3h)
+B3h=180-be2
+mprintf('\n B1m = %f degree',B3h)
+B3m=vf/ut
+mprintf('\n B1m =%f degree ',atand(B3m))
+
+
+
+
+
+
+
diff --git a/3814/CH7/EX7.4/Ex7_4.sce b/3814/CH7/EX7.4/Ex7_4.sce
new file mode 100644
index 000000000..90bd4878e
--- /dev/null
+++ b/3814/CH7/EX7.4/Ex7_4.sce
@@ -0,0 +1,23 @@
+//determine the mean speed of turbine
+// ex 7.4 pgno.174 
+clc
+p=60e6 // power
+h=40// meter
+no=0.85 // overall efficiency
+x=0.5 // flow ratio 
+g=9.8 
+Ku=1.6 // speed ratio
+Q=p/(9800*no*h) // flow in kaplan turbine
+mprintf('\n Q = %d m3/s',Q)
+Vf=x*sqrt(2*g*h) // velocity of flow
+mprintf('\n Vf = %d m/s',Vf)
+d1=(%pi*(1-(0.35^2)))/4
+d=d1*Vf // diameter
+df=sqrt(180/d)
+mprintf('\n D1 = %f m',df)
+dh=0.35*df 
+mprintf('\n Dh = %f m',dh)
+D=(dh+df)/2 // mean diameter 
+mprintf('\n mean diameter D= %f m',D)
+N=(Ku*sqrt(2*g*h)*60)/(%pi*D) //rotational speed 
+mprintf('\n N = %f rev/min',N)
diff --git a/3814/CH7/EX7.5/Ex7_5.sce b/3814/CH7/EX7.5/Ex7_5.sce
new file mode 100644
index 000000000..fc11e5109
--- /dev/null
+++ b/3814/CH7/EX7.5/Ex7_5.sce
@@ -0,0 +1,26 @@
+// determine blade inlet angle and eulers head
+// ex 7.5 pgno.175
+clc
+h=35 //meter
+D=2 //m
+N=145 //rev/min
+alpa1=30// degree
+g=9.8 //
+b1=28// degree
+H=32.6// 0.93*35 head
+V1=sqrt(2*g*H) // inlet velocity
+mprintf('\n inlet velocity V1 = %f m/s',V1)
+u=(%pi*D*N)/60 // 
+mprintf('\n u = %f m/s',u)
+Vr1=sqrt(V1^2+u^2-(2*u*V1*cosd(alpa1))) // inlet triangle of velocity
+mprintf('\n Vr1 = %f m/s',Vr1)
+v=(-u^2+Vr1^2+V1^2)
+v1=2*Vr1*V1
+V=v/v1
+s=acosd(V)
+B1=(180-s) // Beta
+mprintf('\n Beta B1 = %f degree',B1)
+Vlw=V1*cosd(alpa1)
+mprintf('\n Vlw =% f m/s',Vlw)
+E=(u*Vlw)/g // eulers head
+mprintf('\n E =%f m',E)
diff --git a/3814/CH7/EX7.6/Ex7_6.sce b/3814/CH7/EX7.6/Ex7_6.sce
new file mode 100644
index 000000000..af44e2738
--- /dev/null
+++ b/3814/CH7/EX7.6/Ex7_6.sce
@@ -0,0 +1,17 @@
+// design of bersia power station
+// ex 7.6 pgno.178 
+clc
+p=24.7e6// power watt
+h=26.5 //m
+N=187.5 // rev/min
+Q=104 //m3/s
+w=(2*N*%pi)/60
+g=9.8
+mprintf('\n  w= %f rad/s',w)
+wt=(w*sqrt(p/10^3))/(g*h)^(5/4)
+mprintf('\n wt =%f',wt)
+Ns=(N*sqrt(p/10^3))/(h^(5/4)) // speed
+mprintf('\n Ns =%f',Ns)
+n0=p/(9800*Q*h) // overall efficiency
+mprintf('\n Overall efficiency n0= %f percentage',n0*100)
+mprintf('\n Based on specific speed values obtained kaplan turbine is selected with an overall efficiency of %f percentage',n0*100)
diff --git a/3814/CH7/EX7.7/EX7_7.sce b/3814/CH7/EX7.7/EX7_7.sce
new file mode 100644
index 000000000..52f5252d2
--- /dev/null
+++ b/3814/CH7/EX7.7/EX7_7.sce
@@ -0,0 +1,16 @@
+// design of termengor power stations
+// ex 7.7 pgno.178
+ clc
+ Q=125.4 // flow rate at m3/s
+ H=101 // m
+ N=214.3 // speed in rev/min
+ p=90e6 // power to turbine in wattsr
+ g=9.8
+ w=(2*N*%pi)/60
+ mprintf('\n Wt =%f rad/s',w)
+ wt=(w*(sqrt(p/10^3)))/((g*H)^(5/4))
+ mprintf('\n wt= %f',wt)
+ Ns=(N*sqrt(90*1000))/((H)^(5/4))// specific speed
+ mprintf('\n Ns = %f ',Ns)
+ n0=p/(9800*Q*H) // overall efficiency
+ mprintf('\n Over all efficiency n0= %f percentage',n0*100)
diff --git a/3814/CH7/EX7.8/Ex7_8.sce b/3814/CH7/EX7.8/Ex7_8.sce
new file mode 100644
index 000000000..0f2b32b45
--- /dev/null
+++ b/3814/CH7/EX7.8/Ex7_8.sce
@@ -0,0 +1,17 @@
+// to find JOR power station
+// ex 7.8 pgno.179 
+clc
+P=26.1e6 // power in mega watts
+H=587.3 // m
+N=428 // revloution /minutes
+Q=6.85 // m3/s
+w=(2*N*%pi)/60
+g=9.8
+mprintf('\n W= %f rad/s',w)
+wt=(w*(sqrt(P/10^3)))/((g*H)^(5/4))
+mprintf('\n wt =%f',wt)
+Ns=(N*(sqrt(26.1e6)))/(H^(5/4)) // speed
+mprintf('\n Ns = %f',Ns)
+mprintf('\n error in text book they have taken p=26.1e3 instead of 26.1 e6')
+n0=P/(9800*Q*H) // overall efficiency
+mprintf(' \nOverall efficiency n0= %f percentage',n0*100)