summaryrefslogtreecommitdiff
path: root/650/CH8
diff options
context:
space:
mode:
Diffstat (limited to '650/CH8')
-rwxr-xr-x650/CH8/EX8.10/10.sce32
-rwxr-xr-x650/CH8/EX8.11/11.sce14
-rwxr-xr-x650/CH8/EX8.12/12.sce41
-rwxr-xr-x650/CH8/EX8.13/13.sce15
-rwxr-xr-x650/CH8/EX8.14/14.sce20
-rwxr-xr-x650/CH8/EX8.15/15.sce9
-rwxr-xr-x650/CH8/EX8.16/16.sce13
-rwxr-xr-x650/CH8/EX8.4/4.sce17
-rwxr-xr-x650/CH8/EX8.5/5.sce21
-rwxr-xr-x650/CH8/EX8.6/6.sce24
-rwxr-xr-x650/CH8/EX8.7/7.sce27
-rwxr-xr-x650/CH8/EX8.8/8.sce15
-rwxr-xr-x650/CH8/EX8.9/9.sce21
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
generated by cgit (git 2.34.1) at 2024-12-22 21:10:37 +0000