initial commit / add all books

14 files changed, 264 insertions, 0 deletions

diff --git a/1962/CH9/EX9.1/example9_1.sce b/1962/CH9/EX9.1/example9_1.sce
new file mode 100755
index 000000000..e6824b550
--- /dev/null
+++ b/1962/CH9/EX9.1/example9_1.sce
@@ -0,0 +1,27 @@
+

+//example 9.1

+//page 308

+clc; funcprot(0);

+//initialisation of variable

+S=1.26;

+mu=0.862//N.s/m^2

+rhow=998;//density

+rhog=1.26*998;//density of glycerine

+Q=0.1;

+g=9.81;

+L=100;

+D=0.2;

+pi=3.14;

+f=0.0688;//coeff of friction

+A=pi*0.2^2/4;

+V=Q/A;

+R=V*D*rhog/mu;

+hf=f*L/D*V^2/2/g;

+disp(hf,"head loss if glucerine is flowing(m of glycerine)");

+//part2

+mu=1.005*10^-3;

+R=V*D*rhow/mu;

+f=0.021;//coeff. of friction

+hf=f*L/D*V^2/2/g;

+disp(hf,"head loss if water is flowing(m of water)");

+clear

diff --git a/1962/CH9/EX9.10/example9_10.sce b/1962/CH9/EX9.10/example9_10.sce
new file mode 100755
index 000000000..0b9bd44d8
--- /dev/null
+++ b/1962/CH9/EX9.10/example9_10.sce
@@ -0,0 +1,13 @@
+

+//example 9.10

+//page 333

+clc; funcprot(0);

+//initialisation of variable

+f1=0.021//friction in pipe 1

+d1=0.2;

+k=f1*(10+50+30+20+100)/d1^5+2*0.9/d1^4+10/d1^4+2*1.2/d1^4;//k=(HL)1/(16Q^2/2pi^2g)

+f2=0.023//friction pipe 2

+d2=0.15;

+L2=k*d2^5/f2;

+disp(L2,"length of pipe2 (m)=")

+clear

diff --git a/1962/CH9/EX9.11/example9_11.sce b/1962/CH9/EX9.11/example9_11.sce
new file mode 100755
index 000000000..245f74487
--- /dev/null
+++ b/1962/CH9/EX9.11/example9_11.sce
@@ -0,0 +1,27 @@
+

+//example 9.11

+//page 3337

+clc; funcprot(0);

+//initialisation of variable

+Q=2;

+nu=1.007*10^-56;

+d2=0.3;

+d1=0.5;

+d3=0.15;

+L1=500;

+L2=200;

+L3=300;

+pi=3.14;

+g=9.81;

+f=0.018//(solved using iteration from moody's chart)

+A1=pi*d1^2/4;

+A2=pi*d2^2/4;

+A3=pi*d3^2/4;

+//assumption

+q1=0.5;

+q2=sqroot(L1*q1^2/d1^5*d2^5/L2);

+q3=sqroot(L1*q1^2/d1^5*d3^5/L3);

+Q1=2*q1/(q1+q2+q3);

+hf=f*L1*Q1^2/A1^2/2/g/d1;

+disp(hf,"head loss (m)=");

+clear

diff --git a/1962/CH9/EX9.12/example9_12.sce b/1962/CH9/EX9.12/example9_12.sce
new file mode 100755
index 000000000..449987fbb
--- /dev/null
+++ b/1962/CH9/EX9.12/example9_12.sce
@@ -0,0 +1,23 @@
+

+//example 9.12

+//page 340

+clc; funcprot(0);

+//initialisation of variable

+L1=5000;

+D1=1.5;

+EL1=135;//energy level 1

+L2=800;

+D2=0.5;

+EL2=120;

+EL3=112;

+L3=1500;

+D3=0.75;

+//using iteration junction height=131.6;

+

+Q1=sqroot((EL1-131.6)/0.7594);

+disp(Q1,"flow in pipe 1(m^3/s)");

+Q2=sqroot((131.6-EL2)/36.33);

+disp(Q2,"flow in pipe 2(m^3/s)");

+Q3=sqroot((131.6-EL3)/8.168);

+disp(Q3,"flow in pipe 3(m^3/s)");

+clear

diff --git a/1962/CH9/EX9.13/example9_13.sce b/1962/CH9/EX9.13/example9_13.sce
new file mode 100755
index 000000000..7738d4d44
--- /dev/null
+++ b/1962/CH9/EX9.13/example9_13.sce
@@ -0,0 +1,23 @@
+

+//example 9.13

+//page 342

+clc; funcprot(0);

+//initialisation of variable

+L1=5000;

+D1=1.5;

+EL1=136;//energy level 1

+L2=800;

+D2=0.5;

+EL2=122;

+EL3=110;

+L3=1500;

+D3=0.75;

+//using iteration junction height=131;

+j=131;

+Q1=sqroot((EL1-j)/0.7594);

+disp(Q1,"flow in pipe 1(m^3/s)");

+Q2=sqroot((j-EL2)/36.33);

+disp(Q2,"flow in pipe 2(m^3/s)");

+Q3=sqroot((j-EL3)/8.168);

+disp(Q3,"flow in pipe 3(m^3/s)");

+clear

diff --git a/1962/CH9/EX9.14/example9_14.sce b/1962/CH9/EX9.14/example9_14.sce
new file mode 100755
index 000000000..cdf3575f9
--- /dev/null
+++ b/1962/CH9/EX9.14/example9_14.sce
@@ -0,0 +1,17 @@
+

+//example 9.14

+//page 349

+clc; funcprot(0);

+//initialisation of variable

+z=100;//elevation

+f=0.0125;//using iteration

+L=10000;

+L2=1000;

+D=2;

+g=32.2;

+pi=3.14;

+V=sqrt(2*g*z/f/L*D);

+Q=pi*D^2/4*V;

+x=34-10-f*L2*V^2/2/g/D-V^2/2/g;

+disp(30-x,"minimum depth of ridge(ft)=");

+clear

diff --git a/1962/CH9/EX9.2/example9_2.sce b/1962/CH9/EX9.2/example9_2.sce
new file mode 100755
index 000000000..f328c9f12
--- /dev/null
+++ b/1962/CH9/EX9.2/example9_2.sce
@@ -0,0 +1,16 @@
+

+//example 9.2

+//page 311

+clc; funcprot(0);

+//initialisation of variable

+L=100;//length

+D=0.2;//diameter

+hf=5.43;//head loss

+g=9.81

+pi=3.14;

+f=0.021//friction coeff.

+A=sqroot(2*g*D*hf/L);//area

+V=A/sqroot(f);

+Q=V*pi*D^2/4;

+disp(Q,"flow rate (m^3/s)=");

+clear

diff --git a/1962/CH9/EX9.3/example9_3.sce b/1962/CH9/EX9.3/example9_3.sce
new file mode 100755
index 000000000..426e1928f
--- /dev/null
+++ b/1962/CH9/EX9.3/example9_3.sce
@@ -0,0 +1,17 @@
+

+//example 9.3

+//page 311

+clc; funcprot(0);

+//initialisation of variable

+Q=0.1;

+hf=5.43;//head loss

+L=100;

+nu=1.007*10^-6;//kinematic viscosity

+pi=3.14;

+g=9.81;

+A=8*L*Q^2/hf/g/pi^2;

+//using moody's chart

+f=0.021;

+D=(A*f)^0.2;

+disp(D,"diameter of the pipe(m)=");

+clear

diff --git a/1962/CH9/EX9.4/example9_4.sce b/1962/CH9/EX9.4/example9_4.sce
new file mode 100755
index 000000000..3e4a46483
--- /dev/null
+++ b/1962/CH9/EX9.4/example9_4.sce
@@ -0,0 +1,10 @@
+

+//example 9.4

+//page 318

+clc; funcprot(0);

+//initialisation of variable

+L=500;

+S=0.004;//slope of slope line

+hf=S*L;

+disp(hf,"head loss (ft)=");

+clear

diff --git a/1962/CH9/EX9.5/example9_5.sce b/1962/CH9/EX9.5/example9_5.sce
new file mode 100755
index 000000000..ce9ff27b3
--- /dev/null
+++ b/1962/CH9/EX9.5/example9_5.sce
@@ -0,0 +1,16 @@
+

+//example 9.5

+//page 319

+clc; funcprot(0);

+//initialisation of variable

+epsilon=0.025;//roughness

+L=500;

+pi=3.14;

+g=9.81;

+D=30/100;

+Q=0.1;

+K=pi/4*sqroot(2*g)*(2*log10(D*100/epsilon)+1.14)*D^2.5;

+S=(Q/K)^2;

+hf=S*L;

+disp(hf,"head loss(m)=");

+clear

diff --git a/1962/CH9/EX9.6/example9_6.sce b/1962/CH9/EX9.6/example9_6.sce
new file mode 100755
index 000000000..d15b62920
--- /dev/null
+++ b/1962/CH9/EX9.6/example9_6.sce
@@ -0,0 +1,16 @@
+

+//example 9.6

+//page 319

+clc; funcprot(0);

+//initialisation of variable

+epsilon=0.025;//roughness

+L=500;

+pi=3.14;

+g=9.81;

+D=20/100;

+Q=0.1;

+S=5.43/100;

+K=pi/4*sqroot(2*g)*(2*log10(D*100/epsilon)+1.14)*D^2.5;

+Q=K*sqroot(S);

+disp(Q,"flow rate(m^3/s)=");

+clear

diff --git a/1962/CH9/EX9.7/example9_7.sce b/1962/CH9/EX9.7/example9_7.sce
new file mode 100755
index 000000000..7072f5d91
--- /dev/null
+++ b/1962/CH9/EX9.7/example9_7.sce
@@ -0,0 +1,17 @@
+

+//example 9.7

+//page 320

+clc; funcprot(0);

+//initialisation of variable

+epsilon=0.025;//roughness

+L=500;

+pi=3.14;

+g=9.81;

+Q=0.1;

+S=5.43/100;

+K=Q/sqroot(S);

+//solving for D

+deff('y=f(D)','y=3.14/4*sqroot(2*9.81)*(2*log10(D/0.025)+1.14)*D^2.5-.3');

+[x]=fsolve(0.1,f);

+disp(x,"diameter(m):");

+clear

diff --git a/1962/CH9/EX9.8/example9_8.sce b/1962/CH9/EX9.8/example9_8.sce
new file mode 100755
index 000000000..765d9ec61
--- /dev/null
+++ b/1962/CH9/EX9.8/example9_8.sce
@@ -0,0 +1,24 @@
+

+//example 9.8

+//page 326

+clc; funcprot(0);

+//initialisation of variable

+Q=0.075;

+L=30;

+D=0.1;

+pi=3.14;

+k=0.5;

+K=10;

+g=9.81;

+nu=1.007*10^-6//kinematic viscosity

+A=pi*D^2/4;

+V=Q/A;

+R=V*D/nu;

+//using moody's chart

+f=0.025;

+hf1=f*L*V^2/2/g/D;//head loss by friction

+hf2=k*V^2/2/g;//head loss due to contraction

+hf3=K*V^2/2/g;//head loss due to expansion

+hf=hf1+hf2+hf3;

+disp(hf,"total head loss (m of H2O)=");

+clear

diff --git a/1962/CH9/EX9.9/example9_9.sce b/1962/CH9/EX9.9/example9_9.sce
new file mode 100755
index 000000000..17efb5d14
--- /dev/null
+++ b/1962/CH9/EX9.9/example9_9.sce
@@ -0,0 +1,18 @@
+

+//example 9.9

+//page 329

+clc; funcprot(0);

+//initialisation of variable

+nu=1.007*10^-6;

+g=9.81;

+hf=100;

+pi=3.14;

+//on solving V1 and V2 in terms of V3 and

+// iterate for f1 and f2

+//we get 

+f1=0.019;

+f2=0.022;

+V3=sqroot(2*g*hf/(8.4*f1+268.85*f2+4.85));

+Q3=V3*pi*0.08^2/4;

+disp(Q3,"flow in pipe 3(m^3/s)=");

+clear