initial commit / add all books

10 files changed, 166 insertions, 0 deletions

diff --git a/746/CH8/EX8.01/8_01.sce b/746/CH8/EX8.01/8_01.sce
new file mode 100755
index 000000000..cddb2de1f
--- /dev/null
+++ b/746/CH8/EX8.01/8_01.sce
@@ -0,0 +1,15 @@
+//Leakage flow rate//

+pathname=get_absolute_file_path('8.01.sce')

+filename=pathname+filesep()+'8.01-data.sci'

+exec(filename)

+//Leakage flow rate (in mm^3/sec):

+Q=%pi/12*D*a^3*(p1-p2)*10^3/u/L

+//Velocity of flow(in m/sec):

+V=Q/%pi/D/a/1000

+//Specific gravity of SAE 10W oil:

+SG=0.92;

+//Reynolds Number:

+Re=SG*dw*V*a/u/1000

+//As Re<1400, flow is laminar.

+printf("\n\nRESULTS\n\n")

+printf("\n\nLeakage flow rate: %.3f mm^3/sec\n\n",Q)

diff --git a/746/CH8/EX8.02/8_02.sce b/746/CH8/EX8.02/8_02.sce
new file mode 100755
index 000000000..904e37f22
--- /dev/null
+++ b/746/CH8/EX8.02/8_02.sce
@@ -0,0 +1,15 @@
+//Torque and power//

+pathname=get_absolute_file_path('8.02.sce')

+filename=pathname+filesep()+'8.02-data.sci'

+exec(filename)

+//Shear stres (in lbf/ft^2):

+Tyx=u*N*2*%pi/60*D/2/(a/2)

+//Torqe(in inches-lbf):

+T=%pi/2*Tyx*D^2*L/144

+//Power dissipated in the bearing(in hp):

+P=T*N/60*2*%pi/12/550

+//Reynolds number:

+Re=SG*p*N*2*%pi/60*1.5*a/2/u/144

+printf("\n\nRESULTS\n\n")

+printf("\n\nTorque: %.3f inches-lbf\n\n",T)

+printf("\n\nPower dissipated in the bearing: %.3f hp\n\n",P)

diff --git a/746/CH8/EX8.04/8_04.sce b/746/CH8/EX8.04/8_04.sce
new file mode 100755
index 000000000..e65781809
--- /dev/null
+++ b/746/CH8/EX8.04/8_04.sce
@@ -0,0 +1,12 @@
+//Viscosity of fluid//

+pathname=get_absolute_file_path('8.04.sce')

+filename=pathname+filesep()+'8.04-data.sci'

+exec(filename)

+//Viscosity of the liquid(in N-s/m^2):

+u=%pi/128*p*1000*D^4/Q/L/1000

+//Velocity(in m/sec)

+V=Q/(%pi/4*D^2)/1000

+//Reynolds number:

+Re=d*V*D/u/1000

+printf("\n\nRESULTS\n\n")

+printf("\n\nViscosity of fluid %.3f N-s/m^2\n\n",u)

diff --git a/746/CH8/EX8.05/8_05.sce b/746/CH8/EX8.05/8_05.sce
new file mode 100755
index 000000000..7d3d6eac6
--- /dev/null
+++ b/746/CH8/EX8.05/8_05.sce
@@ -0,0 +1,10 @@
+//required//

+pathname=get_absolute_file_path('8.05.sce')

+filename=pathname+filesep()+'8.05-data.sci'

+exec(filename)

+//Reservoir depth required to maintain flow(in m):

+D1=8*Q^2/(%pi)^2/D^4/g*(f*L/D+K+1)

+//Reynolds number:

+Re=4*d*Q/((%pi)*u*D)

+printf("\n\nRESULTS\n\n")

+printf("\n\nReservoir depth required to maintain flow: %.3f m\n\n",D1)

diff --git a/746/CH8/EX8.06/8_06.sce b/746/CH8/EX8.06/8_06.sce
new file mode 100755
index 000000000..9a2b9934a
--- /dev/null
+++ b/746/CH8/EX8.06/8_06.sce
@@ -0,0 +1,13 @@
+//Maximum and power//

+pathname=get_absolute_file_path('8.06.sce')

+filename=pathname+filesep()+'8.06-data.sci'

+exec(filename)

+//Velocity of flow(in ft/sec):

+V=Q/24/3600/(%pi/4*(D/12)^2)*42/7.48

+//Maximum spacing(in ft):

+L=2/f*D/12*(p2-p1)/(SG*d)/V^2*144

+//Power needed at each pump(in hp):

+Win=1/Effp*V*%pi/4*(D/12)^2*(p2-p1)/550*144

+printf("\n\nRESULTS\n\n")

+printf("\n\nMaximum spacing: %.3f feet\n\n",L)

+printf("\n\nPower needed at each pump: %.3f hp\n\n",Win)

diff --git a/746/CH8/EX8.07/8_07.sce b/746/CH8/EX8.07/8_07.sce
new file mode 100755
index 000000000..cc05e2cd7
--- /dev/null
+++ b/746/CH8/EX8.07/8_07.sce
@@ -0,0 +1,10 @@
+//Volume low//

+pathname=get_absolute_file_path('8.07.sce')

+filename=pathname+filesep()+'8.07-data.sci'

+exec(filename)

+//Velocity(in ft/sec):

+V2=sqrt(2*g*l/(f*((L+l)/D*12+8)+1))

+//Volume flow rate(in gpm):

+Q=V2*%pi*(D/12)^2/4*7.48*60

+printf("\n\nRESULTS\n\n")

+printf("\n\nVolume low rate: %.3f\n\n",Q)

diff --git a/746/CH8/EX8.08/8_08.sce b/746/CH8/EX8.08/8_08.sce
new file mode 100755
index 000000000..7d332823e
--- /dev/null
+++ b/746/CH8/EX8.08/8_08.sce
@@ -0,0 +1,25 @@
+//Minimum diameter//

+pathname=get_absolute_file_path('8.08.sce')

+filename=pathname+filesep()+'8.08-data.sci'

+exec(filename)

+//Value of dPmax(in psi):

+dPmax=p1-p2

+//Q in cubic feet/sec:

+Q1=1500/60/7.48;

+//Initially assume diameter to be 4inches:

+D=4;

+//Reynolds number:

+Re=4*Q1/%pi/v/D*12

+//For this value, 

+f=0.012;

+ dP=8*f*L*p*Q1^2/(%pi)^2/D^5*1728;

+while(dP>dPmax)   

+   dP=8*f*L*p*Q1^2/(%pi)^2/D^5*1728;

+   if(dP<dPmax)

+   break

+   else

+   D=D+1;

+ end

+end

+printf("\n\nRESULTS\n\n")

+printf("Minimum diameter that can be used:%.1f inches\n\n",D)

diff --git a/746/CH8/EX8.09/8_09.sce b/746/CH8/EX8.09/8_09.sce
new file mode 100755
index 000000000..f2de249c9
--- /dev/null
+++ b/746/CH8/EX8.09/8_09.sce
@@ -0,0 +1,20 @@
+//Loss Coefficient//

+pathname=get_absolute_file_path('8.09.sce')

+filename=pathname+filesep()+'8.09-data.sci'

+exec(filename)

+//Average velocity (in ft/s):

+V2=4/%pi*Q/D^2*144

+//Reynolds number:

+Re=V2*D/v/12

+//For this value,

+f=0.013;

+//Power law exponent:

+n=-1.7+1.8*log10(Re)

+//Value of V/U:

+v_u=2*n^2/(n+1)/(2*n+1)

+//Value of alpha:

+alpha=(1/v_u)^3*2*n^2/(3+n)/(3+2*n)

+//Loss Coefficient for a square edged entrance:

+K=2*g*h/V2^2-f*L/D*12-alpha;

+printf("\n\nRESULTS\n\n")

+printf("\n\nLoss Coefficient for a square edged entrance: %.3f \n\n",K)

diff --git a/746/CH8/EX8.10/8_10.sce b/746/CH8/EX8.10/8_10.sce
new file mode 100755
index 000000000..29c0d54dd
--- /dev/null
+++ b/746/CH8/EX8.10/8_10.sce
@@ -0,0 +1,20 @@
+//Volume and increase//

+pathname=get_absolute_file_path('8.10.sce')

+filename=pathname+filesep()+'8.10-data.sci'

+exec(filename)

+//Velocity V1(in m/s):

+V1=sqrt(2*g*z0/1.04)

+//Volume flow rate(in m^3/sec):

+Q=V1*%pi*D^2/4

+Kdiff=1-1/A_R^2-Cp

+//For 2nd case:

+//Velocity(in m/s):

+V1=sqrt(2*g*z0/0.59)

+//Volume flow rate(in m^3/s):

+Qd=V1*%pi*D^2/4

+//Increase in discharge after addition of diffuser is:

+dQ=(Qd-Q)/Q*100

+printf("\n\nRESULTS\n\n")

+printf("\n\nVolume flow rate in case1: %.3f m^3/sec\n\n",Q)

+printf("\n\nVolume flow rate in case 2: %.3f m^3/sec\n\n",Qd)

+printf("\n\nIncrease in discharge after addition of diffuser is: %.3f percent\n\n",dQ)

diff --git a/746/CH8/EX8.11/8_11.sce b/746/CH8/EX8.11/8_11.sce
new file mode 100755
index 000000000..1ac0346b6
--- /dev/null
+++ b/746/CH8/EX8.11/8_11.sce
@@ -0,0 +1,26 @@
+//Diameter and head//

+pathname=get_absolute_file_path('8.11.sce')

+filename=pathname+filesep()+'8.11-data.sci'

+exec(filename)

+

+//Value of K*B^2:

+K_B=Q/(%pi/4*D^2)*sqrt(0.5*d1/g/d2/h)

+//Reynods number:

+ReD1=4/%pi*Q/D/v

+//By trial and error method, the value of beta is fixed at:

+betta=0.66;

+//K is then:

+K=K_B/betta^2

+//Diameter of orifice plate(in m):

+Dt=betta*D

+//Value of p3-p2(in N/m^2):

+P1=d1*Q^2/(%pi/4*D^2)^2*(1/0.65/betta^2-1)

+//Value of p1-p2(in N/m^2):

+P2=d2*g*h

+//Head loss between sections 1 and 3(in N-m/kg):

+hLT=(P2-P1)/d1

+//Expressing the permanent pressure as a fractio of the meter differential:

+C=(P2-P1)/P2

+printf("\n\n\nRESULTS\n\n")

+printf("\n\nDiameter of the orifice: %.3f m\n\n",Dt)

+printf("\n\nHead loss between secions 1 and 3: %.3f N-m/kg\n\n",hLT)