initial commit / add all books

4 files changed, 73 insertions, 0 deletions

diff --git a/3785/CH7/EX7.1/Ex7_1.sce b/3785/CH7/EX7.1/Ex7_1.sce
new file mode 100644
index 000000000..d65b2e7a8
--- /dev/null
+++ b/3785/CH7/EX7.1/Ex7_1.sce
@@ -0,0 +1,27 @@
+// Example 7_1

+clc;funcprot(0);

+// Given data

+D=6;// The diameter of a steel pipe in inch

+Q=2000;// Volume flow rate in gpm

+L=1.0;// Length in km

+nu=1.0*10^-6;// Kinematic viscosity in m^2/s

+rho=1*10^3;// The density of water in kg/m^3

+

+// Calculation

+// (a)

+D=D*2.54*10^-2;// m

+Q=(Q*3.782*10^-3)/60;// m^3/s

+Vbar=(4*Q)/(%pi*D^2);// m/s

+Re_D=(Vbar*D)/nu;// Reynolds number

+// (b)

+epsilon=5*10^-5;// physical height in m

+function[X]=frictionfactor(y)

+    X(1)=-(2.0*log10(((epsilon/D)/3.7)+(2.51/(Re_D*sqrt(y(1))))))-(1/sqrt(y(1)));

+endfunction

+// Guessing a value of f=1*10^-2;

+y=[1*10^-2];

+f=fsolve(y,frictionfactor);

+dp=f*((1/2)*rho*Vbar^2)*((L*10^3)/D);// The pressure drop in  Pa

+P=dp*Q;// The power required to maintain the flow in W

+printf("\n(a)Re_D=%1.3e.The How is turbulent since the Reynolds number exceeds the transition value of 2300. \n(b)The pressure drop,deltap=%1.3e Pa \n(c)The power required to maintain the flow,P=%1.3e W",Re_D,dp,P);

+// The answer is varied due to round off error

diff --git a/3785/CH7/EX7.3/Ex7_3.sce b/3785/CH7/EX7.3/Ex7_3.sce
new file mode 100644
index 000000000..c59def7cf
--- /dev/null
+++ b/3785/CH7/EX7.3/Ex7_3.sce
@@ -0,0 +1,19 @@
+// Example 7_3

+clc;funcprot(0);

+// Given data

+L=100;// The length of the ship in m

+A=3*10^3;// Surface area in m^2

+rho=1.03*10^3;// The density of sea water in kg/m^3

+V=8;// Speed in m/s

+epsilon=1*10^-4;// The surface roughness in m

+nu=1*10^-6;// The kinematic viscosity in m^2/s

+

+// Calculation

+Re_L=(V*L)/nu;// The length Reynolds number Re_L

+// If the ship surface were smooth,

+C_D_fp=0.455/(log10(Re_L))^2.58;// The drag coefficient

+// For a rough surface,

+C_D_fp=0.30/(log10(14.7*(L/epsilon))^2.5);// The drag coefficient for a rough surface

+D=((1/2)*rho*V^2)*A*C_D_fp;// The ship's frictional drag force in N

+P=D*V;// The power in MW

+printf("\nThe ships frictional drag force,D=%1.4e N \nThe power required to overcome drag force,DV=%1.3f MW",D,P/10^6);

diff --git a/3785/CH7/EX7.4/Ex7_4.sce b/3785/CH7/EX7.4/Ex7_4.sce
new file mode 100644
index 000000000..7702c9253
--- /dev/null
+++ b/3785/CH7/EX7.4/Ex7_4.sce
@@ -0,0 +1,14 @@
+// Example 7_4

+clc;funcprot(0);

+// Given data

+z_1=1;// m

+z_2=10;// m

+k=0.4;// The von Karman constant 

+ubar_1=6;// m/s

+ubar_2=9;// m/s

+

+// Calculation

+ustar=(ubar_2-ubar_1)/(2.5*log(10));// m/s

+y_0=10/exp(ubar_2/(2.5*ustar));// m

+C_f=(2*ustar^2)/ubar_2^2;// The friction coefficient

+printf("\nu_*=%0.3f m/s \ny_0=%1.2e m \nThe friction coefficient,C_f=%1.2e",ustar,y_0,C_f);

diff --git a/3785/CH7/EX7.6/Ex7_6.sce b/3785/CH7/EX7.6/Ex7_6.sce
new file mode 100644
index 000000000..04609301b
--- /dev/null
+++ b/3785/CH7/EX7.6/Ex7_6.sce
@@ -0,0 +1,13 @@
+// Example 7_6

+clc;funcprot(0);

+// Given data

+x=40;// Fixed distance in m

+V_v=100;// Vehicular speed in m/s

+C_D=1.0;// The truck drag coefficient

+A=9;// The trucks frontal area in m^2

+alpha_w=0.05;

+

+// Calculation

+V_w=V_v*((C_D*A)/(%pi*(2*alpha_w*x)^2));// km/h

+dV=V_v-V_w;// The relative air speed in km/h

+printf("\nThe relative air speed,V_v-V_w=%2.1f km/h",dV)