From b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b Mon Sep 17 00:00:00 2001
From: priyanka
Date: Wed, 24 Jun 2015 15:03:17 +0530
Subject: initial commit / add all books

---
 914/CH6/EX6.1/ex6_1.sce | 14 ++++++++++++++
 914/CH6/EX6.2/ex6_2.sce | 11 +++++++++++
 914/CH6/EX6.3/ex6_3.sce | 30 ++++++++++++++++++++++++++++++
 914/CH6/EX6.5/ex6_5.sce | 17 +++++++++++++++++
 914/CH6/EX6.9/ex6_9.sce | 16 ++++++++++++++++
 5 files changed, 88 insertions(+)
 create mode 100755 914/CH6/EX6.1/ex6_1.sce
 create mode 100755 914/CH6/EX6.2/ex6_2.sce
 create mode 100755 914/CH6/EX6.3/ex6_3.sce
 create mode 100755 914/CH6/EX6.5/ex6_5.sce
 create mode 100755 914/CH6/EX6.9/ex6_9.sce

(limited to '914/CH6')

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=");
+
-- 
cgit