initial commit / add all books

15 files changed, 259 insertions, 0 deletions

diff --git a/2789/CH9/EX9.1/Ex9_1.sce b/2789/CH9/EX9.1/Ex9_1.sce
new file mode 100755
index 000000000..e9f83bd71
--- /dev/null
+++ b/2789/CH9/EX9.1/Ex9_1.sce
@@ -0,0 +1,13 @@
+clear;
+clc;
+
+//page no.281
+d = 6;//inches
+v = 15;//fps
+l = 100;//ft
+h_L = 17.5;//ft
+f = h_L*(d/(12*l))*(2*32.2/v^2);
+V_f = v*sqrt(f/8);
+printf('The friction velocity = %.2f fps',V_f);
+
+//there is an error in the answer given in textbook
diff --git a/2789/CH9/EX9.10/Ex9_10.sce b/2789/CH9/EX9.10/Ex9_10.sce
new file mode 100755
index 000000000..b443c691c
--- /dev/null
+++ b/2789/CH9/EX9.10/Ex9_10.sce
@@ -0,0 +1,15 @@
+clear;

+clc;

+

+//page no. 305

+Q = 90;//gpm

+d = 3;//in

+l = 3000;//ft

+V = Q/(60*7.48*0.25*%pi*(d/12)^2);

+R_h = (d/12)/4;

+C_hw = 140;

+S = (V/(1.318*140*R_h^0.63))^(1/0.54);

+h_L = S*l;

+printf('The loss of head = %.1f ft of water',h_L);

+

+//there is a minute error in the answer given in textbook

diff --git a/2789/CH9/EX9.11/Ex9_11.sce b/2789/CH9/EX9.11/Ex9_11.sce
new file mode 100755
index 000000000..9971f758a
--- /dev/null
+++ b/2789/CH9/EX9.11/Ex9_11.sce
@@ -0,0 +1,20 @@
+clear;

+clc;

+

+//page no. 307

+

+G = 40;// lb/min

+d = 3;// in

+T = 100;// degreeF

+p = 50;// psia

+l = 2000;//ft

+Re = ((G/60)*(d/12))/(0.0491*32.2*4*10^-7);

+f = 0.015;

+gam1 = p*(144/(53.3*(T+460)));

+V1 = (G/60)/(gam1*0.0491);

+a = sqrt(1.4*32.2*53.3*(T+460));

+M1 = V1/a;

+M2_limit = sqrt(1/1.4);

+l = (((1-(M1/M2_limit)^2)/(1.4*M1^2)) - 2*log(M2_limit/M1))*(0.25/0.015);

+p2 = 38.9;//psia, from trial and error method 

+printf('p2 = %.1f psia',p2);

diff --git a/2789/CH9/EX9.12/Ex9_12.sce b/2789/CH9/EX9.12/Ex9_12.sce
new file mode 100755
index 000000000..189384681
--- /dev/null
+++ b/2789/CH9/EX9.12/Ex9_12.sce
@@ -0,0 +1,15 @@
+clear;

+clc;

+

+//page no. 312

+

+d = 12;// in

+D = 24;//in

+theta = 20;//degrees

+G = 10;//cfs

+p = 20;//psi

+V12 = G/(0.25*%pi);

+V24 = V12/4;

+K_L = 0.43;

+p24 = ((p*144/62.4) + (V12^2 /(2*32.2)) - ((V24^2)/(2*32.2)) - K_L*(V12-V24)^2 /(2*32.2))/2.314;

+printf('Pressure in the larger pipe = %.1f psi',p24);

diff --git a/2789/CH9/EX9.13/Ex9_13.sce b/2789/CH9/EX9.13/Ex9_13.sce
new file mode 100755
index 000000000..f1a8303c9
--- /dev/null
+++ b/2789/CH9/EX9.13/Ex9_13.sce
@@ -0,0 +1,18 @@
+clear;

+clc;

+

+//page no. 322

+

+d = 12;// in

+l = 1000;//ft

+h1 = 200;//elevation

+h2 = 250;//elevation

+T = 50;//degreeF

+f1 = 0.030;

+V1 = sqrt((h2-h1)*2*32.2/(0.5+f1*l +1));

+R1 = V1/0.00000141;

+f2 = 0.019;

+V2 = sqrt((h2-h1)*2*32.2/(0.5+f2*l +1));

+R2 = V1/0.00000141;

+Q = 0.25*%pi*(d/12)^2 *V2; 

+printf('Velocity = %.1f fps\n flow rate = %.1f cfs',V2,Q);

diff --git a/2789/CH9/EX9.14/Ex9_14.sce b/2789/CH9/EX9.14/Ex9_14.sce
new file mode 100755
index 000000000..05952413c
--- /dev/null
+++ b/2789/CH9/EX9.14/Ex9_14.sce
@@ -0,0 +1,14 @@
+clear;

+clc;

+

+//page no. 322

+

+l = 200;//ft

+Q = 0.1;//cfs

+del_h = 5;//ft

+T = 50;//degreeF

+d = 0.187;//ft

+V = Q/(0.25*%pi*d^2);

+R = V*d/0.0000141;

+f = (del_h*2*32.2/V^2 -(1+0.5))*(d/l);

+printf('Required diameter of the pipe = %.2f in.',d*12);

diff --git a/2789/CH9/EX9.15/Ex9_15.sce b/2789/CH9/EX9.15/Ex9_15.sce
new file mode 100755
index 000000000..23a5cc14a
--- /dev/null
+++ b/2789/CH9/EX9.15/Ex9_15.sce
@@ -0,0 +1,24 @@
+clear;

+clc;

+

+//page no. 324

+

+Q = 2.5;//cfs

+T = 50;//degreeF

+d1 = 8;//in

+d2 = 6;//in

+l1 = 1000;//ft

+l2 = 2000;//ft

+V8 = Q/(0.25*%pi*(d1/12)^2);

+V6 = Q/(0.25*%pi*(d2/12)^2);

+R8 = V8*0.667/0.0000141;

+f8 = 0.020;

+R6 = V6*0.5/0.0000141;

+f6 = 0.019;

+h_L8 = f8*(l1/0.667)*(V8^2 /(2*32.2));

+h_L6 = f6*(l2/0.5)*(V6^2 /(2*32.2));

+Ep = 100+h_L8+h_L6;

+n = Q*62.4*(Ep)/550;

+V8 = sqrt((30/f8)*2*32.2/(l1/0.667));

+Q_max = V8*0.25*%pi*(d1/12)^2;

+printf('Maximum reliable flow that can be pumped = %.1f cfs',Q_max);

diff --git a/2789/CH9/EX9.16/Ex9_16.sce b/2789/CH9/EX9.16/Ex9_16.sce
new file mode 100755
index 000000000..18e9d2f76
--- /dev/null
+++ b/2789/CH9/EX9.16/Ex9_16.sce
@@ -0,0 +1,19 @@
+clear;

+clc;

+

+//page no. 327

+

+Q = 5;//cfs

+d = 12;//in

+l = 5000;//ft

+h = 70;//ft

+L = 2000;//ft

+K = (h/Q^1.85);

+a = (L/l)*K;

+b = ((l-L)/l)*K;

+Q_ = (h/((b+a*(0.5^(1.85)))))^(1/1.85);

+Q_A = Q_/2;

+Q_B = Q_/2;

+del = Q_-Q;//gain capcaity

+percent = (del/Q)*100;//gain percentage

+printf('The gain of capacity by looping the pipe is %.1f cfs or %d percentage',del,percent);

diff --git a/2789/CH9/EX9.2/Ex9_2.sce b/2789/CH9/EX9.2/Ex9_2.sce
new file mode 100755
index 000000000..7d8d9947e
--- /dev/null
+++ b/2789/CH9/EX9.2/Ex9_2.sce
@@ -0,0 +1,17 @@
+clear;

+clc;

+

+//page no. 285

+

+T = 100;// degreeF

+d = 3;// inches

+Re = 80000;// Reynolds number

+e = 0.006;// inches

+l = 1000;//feet

+f1 = 0.021;//friction factor

+nu = 0.729*10^-5;// sqft/sec

+V = Re*nu/0.25;

+h_L1 = f1*(l/0.25)*(V^2 /(2*32.2));

+f = 0.316/Re^0.25;

+h_L = (f/f1)*h_L1;

+printf('Head loss expected = %.1f ft\n and head loss expected if the pipe were smooth = %.2f ft',h_L1,h_L);

diff --git a/2789/CH9/EX9.3/Ex9_3.sce b/2789/CH9/EX9.3/Ex9_3.sce
new file mode 100755
index 000000000..f6dfe2936
--- /dev/null
+++ b/2789/CH9/EX9.3/Ex9_3.sce
@@ -0,0 +1,14 @@
+clear;

+clc;

+

+//page no. 288

+

+T = 100;//degreeF

+d = 3;// inches

+Re = 80000;// Reynolds number

+e = 0.006;//inches

+l = 1000;//ft

+f = 0.0255;//friction factor

+V = 2.33;//fps

+h_L = f*(l/0.25)*(V^2 /(2*32.2));

+printf('Head loss expected = %.1f ft',h_L);

diff --git a/2789/CH9/EX9.4/Ex9_4.sce b/2789/CH9/EX9.4/Ex9_4.sce
new file mode 100755
index 000000000..3d22da8cd
--- /dev/null
+++ b/2789/CH9/EX9.4/Ex9_4.sce
@@ -0,0 +1,17 @@
+clear;

+clc;

+

+//page no. 290

+

+Q = 100;//gallons per minute

+sg = 0.90;

+nu = 0.0012;// lb-sec/sqft

+d = 3;// in

+l = 1000;//ft

+r = 1;//in

+V = 4.53;//fps

+Re = V*(d/12)*sg*1.935/nu;

+h_L = (64/Re)*(12*l/d)*(V^2 /(2*32.2));

+v = 2*V*(1 - (2/d)^2);

+tau = 62.4*sg*h_L/(2*l*12);

+printf('v = %.2f fps\n h_L = %.1f ft of oil\n tau = %.3f psf',v,h_L,tau);

diff --git a/2789/CH9/EX9.6/Ex9_6.sce b/2789/CH9/EX9.6/Ex9_6.sce
new file mode 100755
index 000000000..c3ed3b8ed
--- /dev/null
+++ b/2789/CH9/EX9.6/Ex9_6.sce
@@ -0,0 +1,22 @@
+clear;

+clc;

+

+//page no. 295

+

+Q = 90;// gallons per minute

+T = 68;//degreeF

+d = 3;// in

+l = 3000;//ft

+r = 1;// in

+f = 0.018;

+V = Q/(60*7.48*0.25*%pi*(d/12)^2);

+Re = V*(d/12)*1.935/(0.000021);

+h_L = f*(l/0.25)*(V^2 /(2*32.2));

+tau_0 = f*1.935*V^2 /8;

+tau1 = 2*tau_0/d;

+v_c = V*(1+4.07*sqrt(f/8));

+v_ = sqrt(tau_0/1.935);

+v1 = v_*(5.50+5.75*log10(v_*(r/(2*12))/0.00001085));

+v1_ = v_c-v_*5.75*log10(0.5*d/(r/2));

+delta = d*32.8/(Re*sqrt(f));

+printf('Head lost = %.1f ft of water\n Wall shear stress = %.3f psf\n the center velocity = %.2f fps\n shearing stress = %.3f psf\n v1 = %.2f fps\n delta = %.4f in.',h_L,tau_0,v_c,tau1,v1_,delta);

diff --git a/2789/CH9/EX9.7/Ex9_7.sce b/2789/CH9/EX9.7/Ex9_7.sce
new file mode 100755
index 000000000..f3ef610d1
--- /dev/null
+++ b/2789/CH9/EX9.7/Ex9_7.sce
@@ -0,0 +1,19 @@
+clear;

+clc;

+

+//page no. 298

+

+d = 12;// in

+v = 10;//fps

+e = 2;//in

+k = 0.002;//relative roughness

+l = 1000;//ft

+f = (1/(1.14+2*log10(1/k)))^2;

+v_c = v*(1+4.07*sqrt(f/8));

+tau_0 = f*1.935*v^2 /8;

+v2 = v_c - tau_0*5.75*log10(0.5*d/e);

+v2_ = 8.48*tau_0 + tau_0*5.75*log10(e/(12*k));

+h_L = f*(l)*v^2 /(2*32.2); 

+printf('f = %.4f\n v_c = %.2f fps\n v2 = %.1f fps\n h_L = %.1f ft of water',f,v_c,v2_,h_L);

+

+//there are small errors in the answer given in textbook

diff --git a/2789/CH9/EX9.8/Ex9_8.sce b/2789/CH9/EX9.8/Ex9_8.sce
new file mode 100755
index 000000000..89b1bdae1
--- /dev/null
+++ b/2789/CH9/EX9.8/Ex9_8.sce
@@ -0,0 +1,15 @@
+clear;

+clc;

+

+//page no. 300

+

+V = 4.08;// fps

+Re = 93800;//Reynolds number

+r = 1;//in

+m = 1/7;

+R = 3;//in

+f = 0.316/(Re^0.25);

+v_c = V/(2/((m+1)*(m+2)));

+v1 = v_c*(r/R)^(1/7);

+tau_0 = f*1.935*V^2 /8;

+printf('f = %.3f\n v_c = %.2f fps\n v1 = %.2f fps\n wall shear = %.3f ps',f,v_c,v1,tau_0);

diff --git a/2789/CH9/EX9.9/Ex9_9.sce b/2789/CH9/EX9.9/Ex9_9.sce
new file mode 100755
index 000000000..cff2ffda0
--- /dev/null
+++ b/2789/CH9/EX9.9/Ex9_9.sce
@@ -0,0 +1,17 @@
+clear;

+clc;

+

+//page no. 302

+

+p = 14.7;//psia

+T = 60;// degreeF

+l = 2000;//ft

+b = 18;//in

+h = 12;// in

+v = 10;// fps

+R_h = (b*h)/(2*12*(b+h));

+Re = v*4*R_h*0.0763/(32.2*0.000000375);

+f = 0.019;

+h_L = f*(l/(4*R_h))*v^2 /(2*32.2);

+del_p = 0.0763*h_L;

+printf('loss of head = %.1f ft of air\n and the pressure drop = %.2f psf = %.3f psi',h_L,del_p,del_p*0.0069);