initial commit / add all books

11 files changed, 298 insertions, 0 deletions

diff --git a/1379/CH5/EX5.1.1/example5_1.sce b/1379/CH5/EX5.1.1/example5_1.sce
new file mode 100755
index 000000000..ad93ab2c7
--- /dev/null
+++ b/1379/CH5/EX5.1.1/example5_1.sce
@@ -0,0 +1,25 @@
+

+

+//exapple 5.1 

+clc; funcprot(0);

+// Initialization of Variable

+rho=999.7;

+g=9.81;

+mu=1.308/1000;

+s=1/6950;

+b=0.65;

+h=32.6/100;

+n=0.016;

+//calculation

+//part1

+A=b*h;

+P=b+2*h;

+m=A/P;

+u=s^.5*m^(2/3)/n;

+Q=A*u

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

+C=u/m^0.5/s^0.5;

+disp(C,"chezy coefficient (m^0.5/s):");

+a=-m*rho*g*s/mu;//delu/dely

+disp(a,"velocity gradient in the channel (s^-1):")

+

diff --git a/1379/CH5/EX5.1.10/example5_10.sce b/1379/CH5/EX5.1.10/example5_10.sce
new file mode 100755
index 000000000..020a780f5
--- /dev/null
+++ b/1379/CH5/EX5.1.10/example5_10.sce
@@ -0,0 +1,25 @@
+

+

+//exapple 5.10 

+clc; funcprot(0);

+// Initialization of Variable

+pi=3.14;

+n=0.022;

+B=5.75;

+s=0.15*pi/180;

+Q=16.8;

+function[y]=normal(x)

+    y=Q-B*x/n*(B*x/(B+2*x))^(2/3)*s^0.5;

+endfunction

+x=fsolve(1.33,normal);

+disp(x,"Normal depth in (m):");

+Dc=(Q^2/g/B^2)^(1/3);

+disp(Dc,"Critical depth in (m):");

+delD=.1;

+D=1.55:.1:2.35

+su=0;

+for i=1:9

+    delL=delD/s*(1-(Dc/D(i))^3)/(1-(x/D(i))^3.33);

+    su=su+delL

+end

+disp(su,"distance in (m) from upstream to that place:")

diff --git a/1379/CH5/EX5.1.11/example5_11.sce b/1379/CH5/EX5.1.11/example5_11.sce
new file mode 100755
index 000000000..4547250d2
--- /dev/null
+++ b/1379/CH5/EX5.1.11/example5_11.sce
@@ -0,0 +1,53 @@
+

+

+//exapple 5.11 

+clc; funcprot(0);

+// Initialization of Variable

+g=9.81;

+q=1.49;

+pi=3.14;

+//calculation

+//part1

+Dc=(q^2/g)^.333;

+disp(Dc,"critical depth in (m):");

+//part2

+n=0.021;

+su=1.85*pi/180;//slope upstream

+sd=0.035*pi/180;//slope downstream

+Dnu=(n*q/sqrt(su))^(3/5);

+Dnu=round(Dnu*1000)/1000;

+disp(Dnu,"normal depth upstream in (m):");

+Dnd=(n*q/sqrt(sd))^(3/5);

+disp(Dnd,"normal depth downstream in (m):");

+//part3

+D2u=-0.5*Dnu*(1-sqrt(1+8*q^2/g/Dnu^3));

+D2u=round(D2u*1000)/1000;

+disp(D2u,"conjugate depth for upstream in (m):");

+D1d=-0.5*Dnd*(1-sqrt(1+8*q^2/g/Dnd^3));

+disp(D1d,"conjugate depth for downstream in (m):");

+//part4

+//accurate method

+delD=.022;

+D=0.987:.022:1.141

+dis=0;

+for i=1:8

+    delL=delD/su*(1-(Dc/D(i))^3)/(1-(Dnu/D(i))^3.33);

+    dis=dis+delL

+end

+disp(dis,"distance in (m) of occurence of jump by accurate method:");

+//not so accurate one

+E1=D2u+q^2/2/g/D2u^2;

+E2=Dnd+q^2/2/g/Dnd^2;

+E2=round(E2*1000)/1000;

+E1=round(E1*1000)/1000;

+ahm=(D2u+Dnd)/2;//av. hydraulic mean

+afv=.5*(q/D2u+q/Dnd);//av. fluid velocity

+i=(afv*0.021/ahm^(2/3))^2;

+l=(E2-E1)/(su-i+0.0002);

+disp(l,"distance in (m) of occurence of jump by not so accurate method:")

+//part5

+rho=998;

+Eu=Dnu++q^2/2/g/Dnu^2;

+Eu=round(Eu*1000)/1000;

+P=rho*g*q*(Eu-E1);

+disp(P/1000,"power loss in hydraulic jump per unit width in (kW):")

diff --git a/1379/CH5/EX5.1.2/example5_2.sce b/1379/CH5/EX5.1.2/example5_2.sce
new file mode 100755
index 000000000..107a6c361
--- /dev/null
+++ b/1379/CH5/EX5.1.2/example5_2.sce
@@ -0,0 +1,19 @@
+

+

+//exapple 5.2 

+clc; funcprot(0);

+// Initialization of Variable

+Q=0.885;

+pi=3.1428;

+s=1/960;

+s=round(s*1000000)/1000000;

+b=1.36;

+n=0.014;

+theta=55*pi/180;

+//calculation

+function[y]=flow(x);

+    a=(x*(b+x/tan(theta)))/(b+2*x/sin(theta));

+    y=a^(2/3)*s^(1/2)*(x*(b+x/tan(theta)))/n-Q;

+endfunction

+x=fsolve(0.1,flow);

+disp(x,"depth of water in (m):")

diff --git a/1379/CH5/EX5.1.3/example5_3.sce b/1379/CH5/EX5.1.3/example5_3.sce
new file mode 100755
index 000000000..04660eca2
--- /dev/null
+++ b/1379/CH5/EX5.1.3/example5_3.sce
@@ -0,0 +1,16 @@
+

+

+//exapple 5.3 

+clc; funcprot(0);

+// Initialization of Variable

+n=0.011;

+h=0.12;

+Q=25/10000;

+//calculation

+deff('y=f(x)','y=1/x^2-1');

+x=fsolve(0.1,f);

+theta=2*atan(x);

+A=h*2*h/tan(theta/2)/2;

+P=2*h*sqrt(2);

+s=Q^2*n^2*P^(4/3)/A^(10/3);

+disp(s,"the slope of channel in (radians):")

diff --git a/1379/CH5/EX5.1.4/example5_4.sce b/1379/CH5/EX5.1.4/example5_4.sce
new file mode 100755
index 000000000..f4ff5d23c
--- /dev/null
+++ b/1379/CH5/EX5.1.4/example5_4.sce
@@ -0,0 +1,34 @@
+

+

+//exapple 5.4 

+clc; funcprot(0);

+// Initialization of Variable

+//part1

+//maximizing eqution in theta & get a function

+function[y]=theta(x)

+    y=(x-.5*sin(2*x))/2/x^2-(1-cos(2*x))/2/x;

+endfunction

+x=fsolve(2.2,theta);

+x=round(x*1000)/1000;

+a=(1-cos(x))/2;

+printf("velocity will be maximum when stream depth in times of diameter is %.3f",a);

+//part2

+//maximizing eqution in theta & get a function

+function[y]=theta2(x)

+    y=3*(x-.5*sin(2*x))^2*(1-cos(2*x))/2/x-(x-.5*sin(2*x))^3/2/x^2 ;   

+endfunction

+x1=fsolve(2.2,theta2);

+x1=round(x1*1000)/1000;

+a=(1-cos(x1))/2;

+disp("")

+printf("vlumetric flow will be maximum when stream depth in times of diameter is %.3f",a);

+//part3

+r=1;

+A=1*x-0.5*sin(2*x);

+s=0.35*3.14/180;

+P=2*x*r;

+C=78.6;

+u=C*(A/P)^0.5*s^0.5;

+disp(u,"maximum velocity of obtained fluid (m/s):");

+//part4

+disp(x1,"maximum flow rate obtained at angle in (radians):")

diff --git a/1379/CH5/EX5.1.5/example5_5.sce b/1379/CH5/EX5.1.5/example5_5.sce
new file mode 100755
index 000000000..6621ed73a
--- /dev/null
+++ b/1379/CH5/EX5.1.5/example5_5.sce
@@ -0,0 +1,31 @@
+

+

+//exapple 5.5 

+clc; funcprot(0);

+// Initialization of Variable

+g=9.81;

+h=28/100;

+Cd=0.62;

+B=46/100;

+Q=0.355;

+n=2;//from francis formula

+//calcualtion

+//part1

+u=sqrt(2*g*h);

+disp(u,"velocity of fluid (m/s):");

+//part2a

+H=(3*Q/2/Cd/B/(2*g)^0.5)^(2/3);

+disp(H,"fluid depth over weir in (m):");

+//part2b

+//using francis formula

+function[y]=root(x)

+    y=Q-1.84*(B-0.1*n*x)*x^1.5;

+endfunction

+x=fsolve(0.2,root);

+disp(x,"fluid depth over weir in if SI units uesd in (m):");

+//part3

+H=18.5/100;

+Q=22/1000;

+a=15*Q/8/Cd/(2*g)^0.5/H^2.5;

+theta=2*atan(a);

+disp(theta*180/3.14,"base angle of the notch of weir (degrees)")

diff --git a/1379/CH5/EX5.1.6/example5_6.sce b/1379/CH5/EX5.1.6/example5_6.sce
new file mode 100755
index 000000000..62dbacf5b
--- /dev/null
+++ b/1379/CH5/EX5.1.6/example5_6.sce
@@ -0,0 +1,26 @@
+

+

+//exapple 5.6 

+clc; funcprot(0);

+// Initialization of Variable

+Q=0.675;

+B=1.65;

+D=19.5/100;

+g=9.81;

+//caculation

+u=Q/B/D;

+u=round(u*1000)/1000;

+E=D+u^2/2/g;

+y=poly([8.53/1000 0 -E 1],'x','coeff');

+x=roots(y);

+disp(x(1),"alternative depth in (m)");

+disp("It is shooting flow");

+Dc=2/3*E;

+Qmax=B*(g*Dc^3)^0.5;

+disp(Qmax,"maximum volumetric flow (m^3/s)");

+Fr=u/sqrt(g*D);

+disp(Fr,"Froude no.");

+a=(E-D)/E;

+disp(a*100,"% of kinetic energy in initial system");

+b=(E-x(1))/E;

+disp(b*100,"% of kinetic energy in final system");

diff --git a/1379/CH5/EX5.1.7/example5_7.sce b/1379/CH5/EX5.1.7/example5_7.sce
new file mode 100755
index 000000000..507f6308d
--- /dev/null
+++ b/1379/CH5/EX5.1.7/example5_7.sce
@@ -0,0 +1,19 @@
+

+

+//exapple 5.7 

+clc; funcprot(0);

+// Initialization of Variable

+G=338;//mass flow rate

+rho=998;

+q=G/rho;

+E=0.48;

+n=0.015;

+g=9.81;

+B=0.4;

+y=poly([5.85/1000 0 -E 1],'x','coeff');

+x=roots(y);

+disp(x(1),x(2),"alternate depths (m):");

+s=(G*n/rho/x(2)/(B*x(2)/(B+2*x(2)))^(2/3))^2

+disp(s,"slode when depth is 12.9cm");

+s=(G*n/rho/x(1)/(B*x(1)/(B+2*x(1)))^(2/3))^2

+disp(s,"slode when depth is 45.1cm");

diff --git a/1379/CH5/EX5.1.8/example5_8.sce b/1379/CH5/EX5.1.8/example5_8.sce
new file mode 100755
index 000000000..47945281f
--- /dev/null
+++ b/1379/CH5/EX5.1.8/example5_8.sce
@@ -0,0 +1,29 @@
+

+

+//exapple 5.8 

+clc; funcprot(0);

+// Initialization of Variable

+pi=3.14;

+theta=pi/3;

+h=1/tan(theta);

+B=0.845;

+E=0.375;

+g=9.81;

+//calculation

+//part1

+//deducing a polynomial(quadratic) in Dc 

+a=5*h;

+b=3*B-4*h*E;

+c=-2*E*B;

+y=poly([c b a],'x','coeff');

+x=roots(y);

+disp(x(2),"critical depth in (m):");

+//part2

+Ac=x(2)*(B+x(2)*tan(theta/2));

+Btc=B+x(2)*tan(theta/2)*2;

+Dcbar=Ac/Btc;

+uc=sqrt(g*Dcbar);

+disp(uc,"critical velocity (m/s):");

+//part3

+Qc=Ac*uc;

+disp(Qc,"Critical volumetric flow (m^3/s):");

diff --git a/1379/CH5/EX5.1.9/example5_9.sce b/1379/CH5/EX5.1.9/example5_9.sce
new file mode 100755
index 000000000..119d8f7c6
--- /dev/null
+++ b/1379/CH5/EX5.1.9/example5_9.sce
@@ -0,0 +1,21 @@
+

+

+//exapple 5.9 

+clc; funcprot(0);

+// Initialization of Variable

+B2=1.60;//breadth at 2

+D2=(1-0.047)*1.27;//depth at 2

+g=9.81;

+B1=2.95;//breadth at 1

+D1=1.27;//depth at 1

+Z=0;

+//calculation

+Q=B2*D2*(2*g*(D1-D2-Z)/(1-(B2*D2/B1/D1)^2))^0.5;

+disp(Q,"volumetric flow rate over flat topped weir over rectangular section in non uniform width(m^3/s)");

+//next part

+B2=12.8;

+D1=2.58;

+Z=1.25;

+Q=1.705*B2*(D1-Z)^1.5;

+disp(Q,"volumetric flow rate over flat topped weir over rectangular section in uniform width (m^3/s):")

+