From b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b Mon Sep 17 00:00:00 2001
From: priyanka
Date: Wed, 24 Jun 2015 15:03:17 +0530
Subject: initial commit / add all books

---
 1962/CH9/EX9.1/example9_1.sce   | 27 +++++++++++++++++++++++++++
 1962/CH9/EX9.10/example9_10.sce | 13 +++++++++++++
 1962/CH9/EX9.11/example9_11.sce | 27 +++++++++++++++++++++++++++
 1962/CH9/EX9.12/example9_12.sce | 23 +++++++++++++++++++++++
 1962/CH9/EX9.13/example9_13.sce | 23 +++++++++++++++++++++++
 1962/CH9/EX9.14/example9_14.sce | 17 +++++++++++++++++
 1962/CH9/EX9.2/example9_2.sce   | 16 ++++++++++++++++
 1962/CH9/EX9.3/example9_3.sce   | 17 +++++++++++++++++
 1962/CH9/EX9.4/example9_4.sce   | 10 ++++++++++
 1962/CH9/EX9.5/example9_5.sce   | 16 ++++++++++++++++
 1962/CH9/EX9.6/example9_6.sce   | 16 ++++++++++++++++
 1962/CH9/EX9.7/example9_7.sce   | 17 +++++++++++++++++
 1962/CH9/EX9.8/example9_8.sce   | 24 ++++++++++++++++++++++++
 1962/CH9/EX9.9/example9_9.sce   | 18 ++++++++++++++++++
 14 files changed, 264 insertions(+)
 create mode 100755 1962/CH9/EX9.1/example9_1.sce
 create mode 100755 1962/CH9/EX9.10/example9_10.sce
 create mode 100755 1962/CH9/EX9.11/example9_11.sce
 create mode 100755 1962/CH9/EX9.12/example9_12.sce
 create mode 100755 1962/CH9/EX9.13/example9_13.sce
 create mode 100755 1962/CH9/EX9.14/example9_14.sce
 create mode 100755 1962/CH9/EX9.2/example9_2.sce
 create mode 100755 1962/CH9/EX9.3/example9_3.sce
 create mode 100755 1962/CH9/EX9.4/example9_4.sce
 create mode 100755 1962/CH9/EX9.5/example9_5.sce
 create mode 100755 1962/CH9/EX9.6/example9_6.sce
 create mode 100755 1962/CH9/EX9.7/example9_7.sce
 create mode 100755 1962/CH9/EX9.8/example9_8.sce
 create mode 100755 1962/CH9/EX9.9/example9_9.sce

(limited to '1962/CH9')

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
-- 
cgit