initial commit / add all books

18 files changed, 357 insertions, 0 deletions

diff --git a/2066/CH7/EX7.1.b/7_1b.sce b/2066/CH7/EX7.1.b/7_1b.sce
new file mode 100755
index 000000000..23e75ddf6
--- /dev/null
+++ b/2066/CH7/EX7.1.b/7_1b.sce
@@ -0,0 +1,16 @@
+clc

+clear

+//Initialization of variables

+z1=2 //ft

+Q=0.1 //gal/min

+alpha=2

+g=32.2 //ft/s^2

+L=4 //ft

+D=1/96 //ft

+//calculations

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

+hl=z1-alpha*v2^2 /(2*g)

+Nr=64/hl *L/D *v2^2 /(2*g)

+mu=v2*D/Nr

+//results

+printf("Kinematic viscosity = %.2e ft^2/s",mu)

diff --git a/2066/CH7/EX7.1.c/7_1c.sce b/2066/CH7/EX7.1.c/7_1c.sce
new file mode 100755
index 000000000..fad25124e
--- /dev/null
+++ b/2066/CH7/EX7.1.c/7_1c.sce
@@ -0,0 +1,16 @@
+clc

+clear

+//Initialization of variables

+z1=2 //ft

+Q=0.1 //gal/min

+alpha=2

+g=32.2 //ft/s^2

+L=4 //ft

+D=1/96 //ft

+//calculations

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

+hl=z1-alpha*v2^2 /(2*g)

+Nr=64/hl *L/D *v2^2 /(2*g)

+Ld=0.058*Nr*D

+//results

+printf("Theoretical entrance transistion length = %.3f ft",Ld)

diff --git a/2066/CH7/EX7.1/7_1.sce b/2066/CH7/EX7.1/7_1.sce
new file mode 100755
index 000000000..bc45119cb
--- /dev/null
+++ b/2066/CH7/EX7.1/7_1.sce
@@ -0,0 +1,16 @@
+clc

+clear

+//Initialization of variables

+z1=2 //ft

+Q=0.1 //gal/min

+alpha=2

+g=32.2 //ft/s^2

+L=4 //ft

+D=1/96 //ft

+//calculations

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

+hl=z1-alpha*v2^2 /(2*g)

+Nr=64/hl *L/D *v2^2 /(2*g)

+//results

+printf("Reynolds number is %d.  Hence the flow is laminar",Nr)

+//The answer is a bit different due to rounding off error in textbook

diff --git a/2066/CH7/EX7.10/7_10.sce b/2066/CH7/EX7.10/7_10.sce
new file mode 100755
index 000000000..0e097bcf9
--- /dev/null
+++ b/2066/CH7/EX7.10/7_10.sce
@@ -0,0 +1,24 @@
+clc

+clear

+//Initialization of variables

+e0=0.00085 //ft

+alpha=0.25 ///year

+t=15 //years

+r0=3 //in

+Q=500 //gal/min

+d=6 //in

+mu=2.04e-5 //lb-sec/ft^2

+rho=1.94 //slugs/ft^3

+g=32.2 //ft/s^2

+L=1 //ft

+gam=62.4

+//calculations

+e15=e0*(1+ alpha*t)

+ratio=r0/(12*e15)

+V=Q/(7.48*60*%pi/4 *(d/12)^2)

+Nr=V*d*rho/(mu*12)

+f=0.036

+hl=f*L/(d/12) *V^2 /(2*g)

+dp=gam*hl

+//results

+printf("Pressure difference = %.2f lb/ft^2 per foot of horizontal pipe",dp)

diff --git a/2066/CH7/EX7.11/7_11.sce b/2066/CH7/EX7.11/7_11.sce
new file mode 100755
index 000000000..7328292b6
--- /dev/null
+++ b/2066/CH7/EX7.11/7_11.sce
@@ -0,0 +1,24 @@
+clc

+clear

+//Initialization of variables

+d2=4 //in

+d1=3 //in

+e=0.0005 //ft

+mu=3.75e-5 //lb-sec/ft^2

+rho=1.94 //slugs/ft^3

+Q=100 //gal/min

+L=100 //ft

+g=32.2 //ft/s^2

+gam=62.4

+//calculations

+A=%pi/4 *((d2/12)^2 -(d1/12)^2)

+WP=%pi*(d1+d2)/12

+R=A/WP

+RR= 2*R/e

+V= Q/(7.48*60*A)

+Nr=V*4*R*rho/mu

+f=0.035

+hl=f*L/(4*R) *V^2 /(2*g)

+hp=hl*Q/(7.48*60) *gam/550

+//results

+printf("horsepower required = %.2f hp/100 ft",hp)

diff --git a/2066/CH7/EX7.12/7_12.sce b/2066/CH7/EX7.12/7_12.sce
new file mode 100755
index 000000000..7e7fb336a
--- /dev/null
+++ b/2066/CH7/EX7.12/7_12.sce
@@ -0,0 +1,18 @@
+clc

+clear

+//Initialization of variables

+p1=25 //psig

+p2=20 //psig

+d1=18 //in

+d2=12 //in

+Cl=0.25

+gam=62.4

+g=32.2 //ft/s^2

+//calculations

+Vr=(d2/d1)^2

+xv=(p2-p1)*144/gam

+V22=xv/(-1-Cl+Vr^2) *2*g

+V2=sqrt(V22)

+Q=V2*%pi/4 *(d2/12)^2

+//results

+printf("Discharge = %.1f ft^3/s",Q)

diff --git a/2066/CH7/EX7.13/7_13.sce b/2066/CH7/EX7.13/7_13.sce
new file mode 100755
index 000000000..cecc01551
--- /dev/null
+++ b/2066/CH7/EX7.13/7_13.sce
@@ -0,0 +1,19 @@
+clc

+clear

+//Initialization of variables

+V61=10.8 //ft/s

+V81=6.05 //ft/s

+r0=3 //in

+e=0.00015

+d1=6 //in

+rho=1.94 //slugs/ft^3

+mu=2.34e-5 //ft-lb/s^2

+//calculations

+roe=r0/(12*e)

+Nr1=V61*(d1/12)*rho/mu

+f6=0.0165

+V6=11.6 //ft/s

+V8=6.52 //ft/s

+Q=V6*%pi/4 *(d1/12)^2

+//results

+printf("Discharge = %.2f ft^3/s",Q)

diff --git a/2066/CH7/EX7.14/7_14.sce b/2066/CH7/EX7.14/7_14.sce
new file mode 100755
index 000000000..dd1894f7a
--- /dev/null
+++ b/2066/CH7/EX7.14/7_14.sce
@@ -0,0 +1,24 @@
+clc

+clear

+//Initialization of variables

+L=1000 //ft

+Q=2000/(7.48*60) //ft63/s

+g=32.2 //ft/s^2

+p=5 //psi/1000 ft

+gam=62.4

+sp=0.7

+f=0.02

+r0=0.904/2

+e=0.00015

+mu=7e-6 //lb-ft/s^2

+L=1000 //ft

+//calculations

+hl=p*144/(sp*gam)

+D5=f*8*L*Q^2 /(%pi^2 *g*hl)

+D=D5^(1/5)

+Nr=4*Q*sp*gam/(g*(%pi*D*mu))

+f2=0.0145

+D5=f2*8*L*Q^2 /(%pi^2 *g*hl)

+D1=D5^(1/5)

+//results

+printf("Diameter of steel pipe = %.3f ft",D1)

diff --git a/2066/CH7/EX7.2/7_2.sce b/2066/CH7/EX7.2/7_2.sce
new file mode 100755
index 000000000..7fac1ee5e
--- /dev/null
+++ b/2066/CH7/EX7.2/7_2.sce
@@ -0,0 +1,18 @@
+clc

+clear

+//Initialization of variables

+Q=350 //gal/min

+D=6 //in

+rho=0.84

+gam=62.4

+g=32.2 //ft/s^2

+mu=9.2e-5 //lb-sec/ft^2

+L=5280 //ft

+//calculations

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

+Nr=V*D/12 *rho*gam/g /mu

+f=0.3164/(Nr)^0.25

+hl=f*L*12/D *V^2 /(2*g)

+hp=hl*gam*Q*rho/(550*7.48*60)

+//resu;ts

+printf("Horsepower required = %.2f hp/mile",hp)

diff --git a/2066/CH7/EX7.3/7_3.sce b/2066/CH7/EX7.3/7_3.sce
new file mode 100755
index 000000000..f474a0ca5
--- /dev/null
+++ b/2066/CH7/EX7.3/7_3.sce
@@ -0,0 +1,10 @@
+clc

+clear

+//Initialization of variables

+n=7

+//calculations

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

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

+//results

+printf("alpha = %.2f",alpha)

+printf("\n beta = %.2f",bet)

diff --git a/2066/CH7/EX7.5/7_5.sce b/2066/CH7/EX7.5/7_5.sce
new file mode 100755
index 000000000..d5e3f0520
--- /dev/null
+++ b/2066/CH7/EX7.5/7_5.sce
@@ -0,0 +1,25 @@
+clc

+clear

+//Initialization of variables

+spg=0.84

+z=1 //in

+gam=62.4

+patm=14.7 //psia

+T=459.6+85 //R

+R=53.3

+g=32.2 //ft/s^2

+D=3 //ft

+mu=3.88e-7 //lb-sec/ft^2

+//calculations

+dp=spg*z/12 *gam

+rho=patm*144/(R*T*g)

+umax=sqrt(2*dp/rho)

+V=0.8*umax

+Nr=V*D*rho/mu

+V2=0.875*umax

+mass=rho*%pi/4 *D^2 *V2

+emf=V2^2 /(2*g)

+hp=emf*mass*g/550

+//results

+printf("Mass flow rate = %.2f slug/sec",mass)

+printf("\n Horsepower input of the fan = %.2f hp",hp)

diff --git a/2066/CH7/EX7.7.a/7_7a.sce b/2066/CH7/EX7.7.a/7_7a.sce
new file mode 100755
index 000000000..752cf1561
--- /dev/null
+++ b/2066/CH7/EX7.7.a/7_7a.sce
@@ -0,0 +1,26 @@
+clc

+clear

+//Initialization of variables

+D=36 //in

+rho=0.00226  //slug/ft^3

+mu=3.88e-7 //lb-sec/ft^2

+umax=62.2 //ft/s

+V=54.5 //ft/s

+Nr=9.5e5

+r0=18 //in

+r=12 //in

+n=8.8

+k=0.4

+//calculations

+f=0.0032 + 0.221/(Nr^0.237)

+Vs=sqrt(f/8) *V

+y=r0-r

+u1=umax*(y/r0)^(1/n)

+u2=umax+ 2.5*Vs*log(y/r0)

+u3=umax+ Vs/k *(sqrt(1-y/r0) + log(1-sqrt(1-y/r0)))

+u4=Vs*(5.5+ 5.75*log10(Vs*y/12 *rho/mu))

+//results

+printf("Using equation 7-13, velocity = %.1f ft/s",u1)

+printf("\n Using equation 7-18, velocity = %.1f ft/s",u2)

+printf("\n Using equation 7-25, velocity = %.1f ft/s",u3)

+printf("\n Using equation 7-34a, velocity = %.1f ft/s",u4)

diff --git a/2066/CH7/EX7.7.b/7_7b.sce b/2066/CH7/EX7.7.b/7_7b.sce
new file mode 100755
index 000000000..dbf46a956
--- /dev/null
+++ b/2066/CH7/EX7.7.b/7_7b.sce
@@ -0,0 +1,23 @@
+clc

+clear

+//Initialization of variables

+D=36 //in

+rho=0.00226  //slug/ft^3

+mu=3.88e-7 //lb-sec/ft^2

+umax=62.2 //ft/s

+V=54.5 //ft/s

+Nr=9.5e5

+r0=18 //in

+r=12 //in

+n=8.8

+k=0.4

+//calculations

+f=0.0032 + 0.221/(Nr^0.237)

+Vs=sqrt(f/8) *V

+y=r0-r

+delta1=D*5*sqrt(8) /(Nr*sqrt(f))

+vss=70

+thick=13*delta1

+//results

+printf("Outer edge of buffer zone is at %d",vss)

+printf("\n Thickness of buffer zone = %.4f in",thick)

diff --git a/2066/CH7/EX7.7.c/7_7c.sce b/2066/CH7/EX7.7.c/7_7c.sce
new file mode 100755
index 000000000..6b894f00a
--- /dev/null
+++ b/2066/CH7/EX7.7.c/7_7c.sce
@@ -0,0 +1,23 @@
+clc

+clear

+//Initialization of variables

+D=36 //in

+rho=0.00226  //slug/ft^3

+mu=3.88e-7 //lb-sec/ft^2

+umax=62.2 //ft/s

+V=54.5 //ft/s

+Nr=9.5e5

+r0=18 //in

+r=12 //in

+n=8.8

+k=0.4

+//calculations

+f=0.0032 + 0.221/(Nr^0.237)

+Vs=sqrt(f/8) *V

+delta1=D*5*sqrt(8) /(Nr*sqrt(f))

+y=delta1

+u2=Vs^2 *delta1/12 *rho/mu

+u1=62.2 *(delta1/18)^(1/n)

+//results

+printf("using equation 7-13, velocity = %.1f ft/s",u1)

+printf("\n using equation 7-30, velocity = %.1f ft/s",u2)

diff --git a/2066/CH7/EX7.7.d/7_7d.sce b/2066/CH7/EX7.7.d/7_7d.sce
new file mode 100755
index 000000000..629c53265
--- /dev/null
+++ b/2066/CH7/EX7.7.d/7_7d.sce
@@ -0,0 +1,23 @@
+clc

+clear

+//Initialization of variables

+D=36 //in

+rho=0.00226  //slug/ft^3

+mu=3.88e-7 //lb-sec/ft^2

+umax=62.2 //ft/s

+V=54.5 //ft/s

+Nr=9.5e5

+r0=18 //in

+r=12 //in

+n=8.8

+k=0.4

+//calculations

+f=0.0032 + 0.221/(Nr^0.237)

+Vs=sqrt(f/8) *V

+delta1=D*5*sqrt(8) /(Nr*sqrt(f))

+y=14*delta1

+u2=62.2*(y/18)^(1/n)

+u3=Vs*(5.50 + 5.75*log10(Vs*y/12 *rho/mu))

+//results

+printf("Using equation 7-13, velocity = %.1f ft/s",u2)

+printf("\n using equation 7-34a, velocity = %.1f ft/s",u3)

diff --git a/2066/CH7/EX7.7.e/7_7e.sce b/2066/CH7/EX7.7.e/7_7e.sce
new file mode 100755
index 000000000..cc48028d0
--- /dev/null
+++ b/2066/CH7/EX7.7.e/7_7e.sce
@@ -0,0 +1,24 @@
+clc

+clear

+//Initialization of variables

+D=36 //in

+rho=0.00226  //slug/ft^3

+mu=3.88e-7 //lb-sec/ft^2

+umax=62.2 //ft/s

+V=54.5 //ft/s

+Nr=9.5e5

+r0=18 //in

+r=12 //in

+n=8.8

+k=0.4

+//calculations

+f=0.0032 + 0.221/(Nr^0.237)

+Vs=sqrt(f/8) *V

+delta1=D*5*sqrt(8) /(Nr*sqrt(f))

+u2=Vs^2 *delta1/12 *rho/mu

+T0=rho*Vs^2

+T02=mu*u2/delta1 *12

+//results

+printf("Using equation 7-9a, shearing stress = %.5f lb/ft^2",T0)

+printf("\n Using equation 7-28, shearing stress = %.5f lb/ft^2",T02)

+disp("The answers are a bit different due to rounding off error in textbook")

diff --git a/2066/CH7/EX7.8/7_8.sce b/2066/CH7/EX7.8/7_8.sce
new file mode 100755
index 000000000..7bbe6e55e
--- /dev/null
+++ b/2066/CH7/EX7.8/7_8.sce
@@ -0,0 +1,12 @@
+clc

+clear

+//Initialization of variables

+umax=62.2 //ft/s

+r0=18 //in

+e=0.0696 //in

+r=6 //in

+//calculations

+Vs=umax/(8.5 + 5.75*log10(r0/e))

+u=Vs*(8.5 + 5.75*log10(r/e))

+//results

+printf("Velocity = %.1f ft/s",u)

diff --git a/2066/CH7/EX7.9/7_9.sce b/2066/CH7/EX7.9/7_9.sce
new file mode 100755
index 000000000..f9988e2d4
--- /dev/null
+++ b/2066/CH7/EX7.9/7_9.sce
@@ -0,0 +1,16 @@
+clc

+clear

+//Initialization of variables

+d=8 //in

+V=3.65 //ft/s

+u1=4.75 //ft/s

+r0=4 //in

+//calculations

+f=0.0449

+Q=V*%pi/4 *(d/12)^2

+Vs=(u1-V)/3.75

+r0e=10^((u1/Vs - 8.5)/5.75)

+e=r0/r0e

+//results

+printf("Flow rate = %.2f ft^3/s",Q)

+printf("\n roughness factor = %.3f in",e)