initial commit / add all books

5 files changed, 88 insertions, 0 deletions

diff --git a/914/CH6/EX6.1/ex6_1.sce b/914/CH6/EX6.1/ex6_1.sce
new file mode 100755
index 000000000..11fa71046
--- /dev/null
+++ b/914/CH6/EX6.1/ex6_1.sce
@@ -0,0 +1,14 @@
+clc;

+warning("off");

+printf("\n\n example6.1 - pg200");

+// given

+q=50;  //[gal/min] - volumetric flow rate

+d=2.067/12;  //[ft] - diameter

+A=0.02330;  //[ft^2] - flow area

+p=0.99568*62.43;  //[lb/ft^3] - density of water at 86degF

+mu=0.8007*6.72*10^-4;  //[lb/ft*sec] - viscosity of water at 86degF

+u=q/(60*7.48*A);

+// using the formula Nre=d*u*p/mu;

+Nre=(d*u*p)/mu;

+disp(Nre,"Nre=");

+printf("\n Hence the flow is turbulent.Note also that Nre is dimensionless");

diff --git a/914/CH6/EX6.2/ex6_2.sce b/914/CH6/EX6.2/ex6_2.sce
new file mode 100755
index 000000000..5369e747e
--- /dev/null
+++ b/914/CH6/EX6.2/ex6_2.sce
@@ -0,0 +1,11 @@
+clc;

+warning("off");

+printf("\n\n example6.2 - pg202");

+// given

+p=0.99568*62.43;  //[lb/ft^3] - density of water at 86degF

+mu=0.8007*6.72*10^-4;  //[lb/ft*sec] - viscosity of water at 86degF

+u=4.78;  //[ft/sec] - free stream velocity 

+Nre=5*10^5;  // the lower limit for the transition reynolds number range is substituted

+x=(Nre*mu)/(p*u);

+disp(x,"x");

+printf("\nThus the transition could star at about %fft.The reynolds number at the upper end of the transition range is %e.The value of x at this location is ten times then the value obtained above i.e %fft",x,Nre*10,x*10); 
\ No newline at end of file
diff --git a/914/CH6/EX6.3/ex6_3.sce b/914/CH6/EX6.3/ex6_3.sce
new file mode 100755
index 000000000..71488ba5d
--- /dev/null
+++ b/914/CH6/EX6.3/ex6_3.sce
@@ -0,0 +1,30 @@
+clc;

+warning("off");

+printf("\n\n example6.3 - pg212");

+// given

+t=[0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.10 0.11 0.12];

+Ux=[3.84 3.50 3.80 3.60 4.20 4.00 3.00 3.20 3.40  3.00 3.50 4.30 3.80];

+Uy=[0.43 0.21 0.18 0.30 0.36 0.28 0.35 0.27 0.21 0.22 0.23 0.36 0.35];

+Uz=[0.19 0.16 0.17 0.13 0.09 0.10 0.16 0.15 0.13 0.18 0.17 0.18 0.17];

+// using the formula AREA=(deltat/2)*(U1+U13+2*(U2+U3+U4+U5+U6+U7+U8+U9+U10+U11+U12))

+// for Uxmean

+deltat=0.01;

+T=t(13)-t(1);

+AREA=(deltat/2)*(Ux(1)+Ux(13)+2*(Ux(2)+Ux(3)+Ux(4)+Ux(5)+Ux(6)+Ux(7)+Ux(8)+Ux(9)+Ux(10)+Ux(11)+Ux(12)));

+Uxmean=AREA/T;

+disp(Uxmean,"Uxmean=");

+// for Uymean

+deltat=0.01;

+T=t(13)-t(1);

+AREA=(deltat/2)*(Uy(1)+Uy(13)+2*(Uy(2)+Uy(3)+Uy(4)+Uy(5)+Uy(6)+Uy(7)+Uy(8)+Uy(9)+Uy(10)+Uy(11)+Uy(12)));

+Uymean=AREA/T;

+disp(Uymean,"Uymean=");

+// for Uzmean

+deltat=0.01;

+T=t(13)-t(1);

+AREA=(deltat/2)*(Uz(1)+Uz(13)+2*(Uz(2)+Uz(3)+Uz(4)+Uz(5)+Uz(6)+Uz(7)+Uz(8)+Uz(9)+Uz(10)+Uz(11)+Uz(12)));

+Uzmean=AREA/T;

+disp(Uzmean,"Uzmean=");

+U=(Uxmean^2+Uymean^2+Uzmean^2)^(1/2);

+disp(U,"U=");

+

diff --git a/914/CH6/EX6.5/ex6_5.sce b/914/CH6/EX6.5/ex6_5.sce
new file mode 100755
index 000000000..cec19421e
--- /dev/null
+++ b/914/CH6/EX6.5/ex6_5.sce
@@ -0,0 +1,17 @@
+clc;

+warning('off');

+printf("\n\n example6.5 - pg232");

+// given

+UzmaxbyU=24.83;

+roUbyv=2312;

+Re=100000;

+// using the formula Et/v=95.5*((r/ro)/slope)-1

+// from fig 6.6 at Re=100000

+rbyro=[0 0.040 0.100 0.200 0.300 0.4 0.5 0.6 0.7 0.8 0.9 0.960 1];

+slope=[0 0.105 0.112 0.126 0.144 0.168 0.201 0.252 0.336 0.503 1.007 2.517 94.59];

+for i=2:13

+    Etbyv(i)=95.5*((rbyro(i))/slope(i))-1;

+end

+clf;

+xtitle("eddy viscosity ratio versus dimensionless radius","r/ro","Et/v");

+plot(rbyro,Etbyv);

diff --git a/914/CH6/EX6.9/ex6_9.sce b/914/CH6/EX6.9/ex6_9.sce
new file mode 100755
index 000000000..e4c5b5c05
--- /dev/null
+++ b/914/CH6/EX6.9/ex6_9.sce
@@ -0,0 +1,16 @@
+clc;

+warning("off");

+printf("\n\n example6.9 - pg258");

+// given

+spg=0.84;

+p=0.84*62.4;  //[lbf/ft^3] - density

+dP=80*144;  //[lbf/ft^2] - pressure

+dz=2000;  //[ft] - length of pipe

+gc=32.174;  //[(lbm*ft)/(lbf*sec^2)] - gravitational conversion constant

+dpbydz=-dP/dz;

+do=2.067/12;  //[ft]

+U=2000*(1/24)*(1/3600)*(42)*(1/7.48)*(1/0.02330);

+// using the formula f=((do/2)*(-dp/dz)*gc)/(p*(U)^2)

+f=((do/2)*(-dpbydz)*gc)/(p*(U)^2)

+disp(f,"f=");

+