initial commit / add all books

13 files changed, 269 insertions, 0 deletions

diff --git a/650/CH8/EX8.10/10.sce b/650/CH8/EX8.10/10.sce
new file mode 100755
index 000000000..cd610176e
--- /dev/null
+++ b/650/CH8/EX8.10/10.sce
@@ -0,0 +1,32 @@
+clc

+m=12*10^3/3600; // kg/s

+Rho=815; // kg/m^3

+d=0.05; // m

+e=0.02;

+d1=50; // m

+d2=0.038; // m

+g=9.8; // m

+

+v=4*m/Rho/%pi/d^2;

+

+f1=1/(2*log10(d1/e)+2.28)^2;

+

+L_eq=d1+2*d1*d;

+

+H_50mm=4*f1*L_eq*v^2/(d*2*g);

+

+v=4*m/(Rho*%pi*d2^2);

+

+f2=1/(2*log10(38/e)+2.28)^2;

+

+L_eq=d1+2*d1*d2;

+H_38mm=4*f2*L_eq*v^2/(d2*2*g);

+

+Hr=0.2*v^2/(2*g);

+

+H_L=H_50mm+H_38mm+Hr;

+

+del_p_f=Rho*g*H_L;

+disp("The total pressure drop due to friction through the pipe system is ")

+disp(del_p_f);

+disp("N/m^2")

diff --git a/650/CH8/EX8.11/11.sce b/650/CH8/EX8.11/11.sce
new file mode 100755
index 000000000..8bf60a549
--- /dev/null
+++ b/650/CH8/EX8.11/11.sce
@@ -0,0 +1,14 @@
+clc

+// H_L=1.2*v^2/2/g

+

+// H_L=4*f*L_eq/d*V^2/2/g

+

+// L_eq=60*d

+

+// H_L=240*f*v^2/2/g

+// Combining the two equations for head loss 

+// 1.2*v^2/2/g=240*f*v^2/2/g

+

+f=1.2/240;

+disp("Friction factor =")

+disp(f)
\ No newline at end of file
diff --git a/650/CH8/EX8.12/12.sce b/650/CH8/EX8.12/12.sce
new file mode 100755
index 000000000..30447d572
--- /dev/null
+++ b/650/CH8/EX8.12/12.sce
@@ -0,0 +1,41 @@
+clc

+// dp_AB+dp_BC=dp_AD+dp_DC

+

+// dp_AD=2*f*rho*v^2*L/d

+

+// dp_AD=16600*(3-Q)^2

+// Likewise 

+// dp_AB=16600*Q^2

+// dp_BC=16600*(Q+0.5)^2

+// dp_DC=16600*(2.1-Q)^2

+// By solving above 5 equations, we get

+

+Q=1.175; //litres per second

+

+disp("The rate of flow from B to C =")

+disp(Q+0.5)

+disp("litres per second")

+

+dp_AD=16600*(3-Q)^2;

+dp_AB=16600*Q^2;

+dp_BC=16600*(Q+0.5)^2;

+dp_DC=16600*(2.1-Q)^2;

+

+disp("dp_AD =")

+disp(dp_AD/1000)

+disp("kN/m^2")

+

+disp("dp_AB =")

+disp(dp_AB/1000)

+disp("kN/m^2")

+

+disp("dp_BC =")

+disp(dp_BC/1000)

+disp("kN/m^2")

+

+disp("dp_DC =")

+disp(dp_DC/1000)

+disp("kN/m^2")

+

+

+disp("The lowest pressure drop is in the pipe connecting C and D")
\ No newline at end of file
diff --git a/650/CH8/EX8.13/13.sce b/650/CH8/EX8.13/13.sce
new file mode 100755
index 000000000..86735e71e
--- /dev/null
+++ b/650/CH8/EX8.13/13.sce
@@ -0,0 +1,15 @@
+clc

+H2=0.5; //m

+H1=2; //m

+A=4; //m^2

+f=0.005; 

+L=20; //m

+d=0.025; //m

+g=9.81; // m/s^2

+

+a=%pi*d^2/4;

+

+t=integrate('-A*sqrt((4*f*L/d)+2.5)/a/(sqrt(2*g))*(H)^(-1/2)', 'H', H1, H2);

+disp("Time taken =")

+disp(t)

+disp("s")
\ No newline at end of file
diff --git a/650/CH8/EX8.14/14.sce b/650/CH8/EX8.14/14.sce
new file mode 100755
index 000000000..599beeb26
--- /dev/null
+++ b/650/CH8/EX8.14/14.sce
@@ -0,0 +1,20 @@
+clc

+

+d0=0.15; // m

+d1=0.1; // m

+Q=50/3600; // m^3/s

+f=0.0052;

+Rho=972;

+

+a=%pi/4*((d0)^2-(d1)^2);

+

+P=%pi*((d0)+(d1));

+

+d_eq=4*a/P;

+

+v=Q/a;

+

+del_p_f=2*f*Rho*v^2/d_eq;

+disp("the pressure drop due to friction per metre length of tube is found to be ")

+disp(del_p_f)

+disp("Nm^2/m")
\ No newline at end of file
diff --git a/650/CH8/EX8.15/15.sce b/650/CH8/EX8.15/15.sce
new file mode 100755
index 000000000..a8bc08a96
--- /dev/null
+++ b/650/CH8/EX8.15/15.sce
@@ -0,0 +1,9 @@
+clc

+f=0.005;

+Q=0.07; // m^3/s

+g=9.81; // m/s^2

+

+H_f=integrate('32*f*(Q)^(2)/(%pi)^(2)/g/(0.3-0.0666*L)^(5)', 'L', 0, 3);

+disp("Fractional head loss =")

+disp(H_f)

+disp("m")
\ No newline at end of file
diff --git a/650/CH8/EX8.16/16.sce b/650/CH8/EX8.16/16.sce
new file mode 100755
index 000000000..f5a54c023
--- /dev/null
+++ b/650/CH8/EX8.16/16.sce
@@ -0,0 +1,13 @@
+clc

+g=9.81; // m/s^2

+H=4; // m

+f=0.006;

+L=50; // m

+d=0.1; // m

+

+v1=sqrt(2*g*H/(4*f*L/d + 1.3));

+

+t=integrate('4/(v1^2-v^2)', 'v', 0, 0.99*v1);

+disp("Time taken =")

+disp(t)

+disp("s")
\ No newline at end of file
diff --git a/650/CH8/EX8.4/4.sce b/650/CH8/EX8.4/4.sce
new file mode 100755
index 000000000..d31cc4e39
--- /dev/null
+++ b/650/CH8/EX8.4/4.sce
@@ -0,0 +1,17 @@
+clc

+rho=867; // kg/m^3

+Q=12/3600; // m^3/s

+u=7.5*10^(-4); // Ns/m^2

+L=200; // m

+H=10; // m

+g=9.81; // m/s^2

+

+d=(H*2*g/(4*0.079*(4*rho*Q/%pi/u)^(-1/4)*L*(4*Q/%pi)^2))^(-4/19);

+disp("Internal diameter of the pipeline =")

+disp(d)

+disp("m")

+

+Re=4*rho*Q/%pi/d/u;

+disp("Re =")

+disp(Re)

+disp("The value of Reynolds number lies between 4000 and 10^5, confirming the validity of using the Blasius equation for smooth-walled pipes")
\ No newline at end of file
diff --git a/650/CH8/EX8.5/5.sce b/650/CH8/EX8.5/5.sce
new file mode 100755
index 000000000..cafe99b0a
--- /dev/null
+++ b/650/CH8/EX8.5/5.sce
@@ -0,0 +1,21 @@
+clc

+m=40/60; // kg/s

+rho=873; // kg/m^3

+d=0.025; // m

+u=8.8*10^-4; // Ns/m^2

+dp=55*10^3; //N/m^2

+L=18; // m

+g=9.81; // m/s^2

+

+v2=4*m/rho/%pi/d^2;

+Re=rho*v2*d/u;

+

+// According to this value of Re, Prandtl's equation is satisfied. 

+// 1/sqrt(f)=4*log(Re*sqrt(f))-0.4

+// By trial and error method we get friction factor equal to

+f=0.0055;

+

+H=dp/rho/g + v2^2/2/g + v2^2/2/g*(4*f*L/d+1.5);

+disp("The minimum allowable height =")

+disp(H)

+disp("m")
\ No newline at end of file
diff --git a/650/CH8/EX8.6/6.sce b/650/CH8/EX8.6/6.sce
new file mode 100755
index 000000000..455c96b6b
--- /dev/null
+++ b/650/CH8/EX8.6/6.sce
@@ -0,0 +1,24 @@
+clc;

+Q=15/3600; // m^3/s

+d=0.05; // m

+Rho=780;

+u=1.7*10^(-3); // Ns/m^2

+f=0.0065;

+L=100; // m

+g=9.8; // m^2/s

+

+v=4*Q/%pi/d^2;

+

+del_pf=2*f*Rho*v^2*L/d;

+disp("The pressure drop due to friction is ")

+disp(del_pf);

+disp("kNm-2")

+

+H_f=4*f*L*v^2/(d*2*g);

+H_exit=v^2/2/g;

+H_entrance=v^2/4/g;

+

+H=H_f+H_exit+H_entrance;

+disp("and the difference in levels is")

+disp(H);

+disp("m");
\ No newline at end of file
diff --git a/650/CH8/EX8.7/7.sce b/650/CH8/EX8.7/7.sce
new file mode 100755
index 000000000..0c1fd0c42
--- /dev/null
+++ b/650/CH8/EX8.7/7.sce
@@ -0,0 +1,27 @@
+clc

+f=0.005; 

+L=10; // m

+d=0.025; // m

+g=9.81; // m/s^2

+

+// H_L=4*f*L/d*v^2/2/g+0.5*v^2/2/g

+// H_L=8.5*v^2/2/g

+

+// By Bernoulli equation we get

+// H=2.62+9.5*v2^2/2/g

+

+// Applying the Bernoulli equation between the liquid surface and discharge point 

+// H_L=33.5*v2^2/2/g

+

+// Solving above two we get

+v2=1.9; // m/s

+

+Q=%pi*d^2/4*v2;

+disp("Rate of flow =")

+disp(Q)

+disp("m^3/s")

+

+H=2.62+9.5*v2^2/2/g;

+disp("The minimum allowable height =")

+disp(H)

+disp("m")
\ No newline at end of file
diff --git a/650/CH8/EX8.8/8.sce b/650/CH8/EX8.8/8.sce
new file mode 100755
index 000000000..be40a1922
--- /dev/null
+++ b/650/CH8/EX8.8/8.sce
@@ -0,0 +1,15 @@
+clc;

+d_A=0.025; // m

+v_A=1.21; // m/s

+d_B=0.05; // m

+v_B=1.71; // m/s

+

+Q_A=%pi*d_A^2*v_A/4;

+disp("The rate of flow through parallel pipes A is ")

+disp(Q_A);

+disp("m^3/s")

+

+Q_B=%pi*d_B^2*v_B/4;

+disp("The rate of flow through parallel pipes B is ")

+disp(Q_B);

+disp("m^3/s")
\ No newline at end of file
diff --git a/650/CH8/EX8.9/9.sce b/650/CH8/EX8.9/9.sce
new file mode 100755
index 000000000..330887acd
--- /dev/null
+++ b/650/CH8/EX8.9/9.sce
@@ -0,0 +1,21 @@
+clc

+d2=0.06; // m

+d1=0.12; // m

+k=0.44;

+f=0.05;

+L1=500; // m

+g=9.81; // m/s^2

+

+// v1=d2^2/d1^2*v2

+

+// H_f=4*f*L1/16/d*v2^2/2/g

+// H_c=k*v2^2/2/g

+// H_f=4*f*L2/d*v2^2/2/g

+// H_exit=v2^2/2/g

+

+v2=sqrt(30*2*g/173.4);

+

+Q=%pi*d2^2/4*v2;

+disp("The rate of flow =")

+disp(Q)

+disp("m^3/s")
\ No newline at end of file