initial commit / add all books

7 files changed, 141 insertions, 0 deletions

diff --git a/812/CH9/EX9.01/9_01.sce b/812/CH9/EX9.01/9_01.sce
new file mode 100755
index 000000000..755aef7b6
--- /dev/null
+++ b/812/CH9/EX9.01/9_01.sce
@@ -0,0 +1,8 @@
+//static pressure//

+pathname=get_absolute_file_path('9.01.sce')

+filename=pathname+filesep()+'9.01-data.sci'

+exec(filename)

+//Change in static pressure between sections 1 and 2:

+C=(((L-2*d1)/(L-2*d2))^4-1)*100;

+printf("\n\nRESULTS\n\n")

+printf("\n\nChange in static pressure between the sections 1 and 2: %.3f percent \n\n",C)

diff --git a/812/CH9/EX9.04/9_04.sce b/812/CH9/EX9.04/9_04.sce
new file mode 100755
index 000000000..db1fd7c5a
--- /dev/null
+++ b/812/CH9/EX9.04/9_04.sce
@@ -0,0 +1,41 @@
+//Displacement thickness and stress//

+pathname=get_absolute_file_path('9.04.sce')

+filename=pathname+filesep()+'9.04-data.sci'

+exec(filename)

+//Reynolds number:

+ReL=U*L/v

+//FOR TURBULENT FLOW

+//Disturbance thickness(in m):

+dL1=0.382/ReL^0.2*L

+//Displacement thickness(in m):

+function y=f(n),y=dL1*(1-n^(1/7))

+endfunction

+dl1=intg(0,1,f)

+//Skin friction coefficient:

+Cf1=0.0594/ReL^0.2

+//Wall shear stress(in N/m^2):

+tw1=Cf1*0.5*d*U^2

+//For LAMINAR FLOW:

+//Disturbance thickness(in m)

+dL2=5/sqrt(ReL)*L

+//Displacement thickness(in m):

+dl2=0.344*dL2

+//Skin friction coefficient:

+Cf2=0.664/sqrt(ReL)

+//Wall shear stress(in N/m^2):

+tw2=Cf2*0.5*d*U^2

+//COMPARISON OF VALUES WITH LAMINAR FLOW

+//Disturbance thickness

+D=dL1/dL2

+//Displacement thickness

+DS=dl1/dl2

+//Wall shear stress

+WSS=tw1/tw2

+printf("\n\nRESULTS\n\n")

+printf("\n\nDisturbace thickness: %.3f m\n\n",dL1)

+printf("\n\nDisplacement thickness: %.3f m\n\n",dl1)

+printf("\n\nWall shear stress: %f N/m^2\n\n",tw1)

+printf("\n\nCOMPARISON WIH LAMINAR FLOW\n\n\n")

+printf("\n\n Disturbance thicknes: %.3f \n\n",D)

+printf("\n\nDisplacement thickness: %.3f\n\n",DS)

+printf("\n\nWall shear stress: %.3f \n\n",WSS)

diff --git a/812/CH9/EX9.05/9_05.sce b/812/CH9/EX9.05/9_05.sce
new file mode 100755
index 000000000..8e146f23f
--- /dev/null
+++ b/812/CH9/EX9.05/9_05.sce
@@ -0,0 +1,19 @@
+//force and power//

+pathname=get_absolute_file_path('9.05.sce')

+filename=pathname+filesep()+'9.05-data.sci'

+exec(filename)

+//Speed in m/s:

+U=s*6076*0.305/3600

+//Reynolds number:

+Re=U*L/v

+//Drag coefficient:

+Cd=0.455/log10(Re)^2.58-1610/Re

+//Area(in m^2):

+A=L*(W+D)

+//Drag force(in N)

+Fd=Cd*A*0.5*d*U^2

+//Power required to overcome skin friction drag(in W):

+P=Fd*U

+printf("\n\nRESULTS\n\n")

+printf("\n\nDrag force: %f N\n\n",Fd)

+printf("\n\nPower required to overcome skin friction drag: %.3f W\n\n",P)

diff --git a/812/CH9/EX9.06/9_06.sce b/812/CH9/EX9.06/9_06.sce
new file mode 100755
index 000000000..3786c2baa
--- /dev/null
+++ b/812/CH9/EX9.06/9_06.sce
@@ -0,0 +1,17 @@
+//Bending moment//

+pathname=get_absolute_file_path('9.06.sce')

+filename=pathname+filesep()+'9.06-data.sci'

+exec(filename)

+//Velocity in m/sec:

+V=s*5/18

+//Reynolds number:

+Re=d*V*D/u

+//Value of Cd is obtained as:

+Cd=0.35;

+//Area(in m^2):

+A=L^2;

+//Moment about the chimney base(in N-m):

+M0=Cd*A*D/4*d*V^2

+printf("\n\nRESULTS\n\n")

+printf("\n\nBending moment at the bottom of the chimney: %.3f N-m\n\n",M0)

+

diff --git a/812/CH9/EX9.07/9_07.sce b/812/CH9/EX9.07/9_07.sce
new file mode 100755
index 000000000..3f1d56c77
--- /dev/null
+++ b/812/CH9/EX9.07/9_07.sce
@@ -0,0 +1,8 @@
+//Time required//

+pathname=get_absolute_file_path('9.07.sce')

+filename=pathname+filesep()+'9.07-data.sci'

+exec(filename)

+//Time required to decelerate to 100 mph(in seconds):

+t=(s1-s2)*2*w/(s1*s2)/Cd/d/A/g*3600/5280

+printf("\n\nRESULTS\n\n")

+printf("\n\nTime required to decelerate to 100 mph: %.3f seconds\n\n",t)

diff --git a/812/CH9/EX9.08/9_08.sce b/812/CH9/EX9.08/9_08.sce
new file mode 100755
index 000000000..54e4c9eca
--- /dev/null
+++ b/812/CH9/EX9.08/9_08.sce
@@ -0,0 +1,28 @@
+//Optimum cruise speed//

+pathname=get_absolute_file_path('9.08.sce')

+filename=pathname+filesep()+'9.08-data.sci'

+exec(filename)

+//Plotting velocity with drag force

+V=175:25:455;

+

+[m n]=size(V);

+for i=1:n

+  CL(i)=2*W/p*(3600/V(i)/5280)^2/A;

+  Cd(i)=Cd0+CL(i)^2/%pi/ar;

+  Fd(i)=Cd(i)/CL(i)*W;

+  FD(i)=Fd(i)/1000;

+end

+plot(V,FD)

+xtitle('Flight speed vs thrust','Flight Speed(in mph)','Drag Force(in 1000lbf)')

+//Optimum cuise speed at speed level is obtained to be 320 mph from the graph.

+Vosl=320;

+//Ratio of speeds at 30000 ft and at sea level is given by:

+r=sqrt(1/0.375);

+//Stall speed at 30000ft is(in mph):

+Vs3=Vssl*r;

+//Optimum Cruise speed at 30000ft(in mph):

+Vo3=Vosl*r;

+printf("\n\nRESULTS\n\n")

+printf("\n\nOptimum cruise speed at sea level: %.3f mph\n\n",Vosl)

+printf("\n\nStall speed at 30000 ft: %.3f mph\n\n",Vs3)

+printf("\n\nOptimum cruise speed at 30000 ft: %.3f\n\n",Vo3)

diff --git a/812/CH9/EX9.09/9_09.sce b/812/CH9/EX9.09/9_09.sce
new file mode 100755
index 000000000..cef9919e3
--- /dev/null
+++ b/812/CH9/EX9.09/9_09.sce
@@ -0,0 +1,20 @@
+//Aerodynamic and Radius//

+pathname=get_absolute_file_path('9.09.sce')

+filename=pathname+filesep()+'9.09-data.sci'

+exec(filename)

+//Reynolds number:

+//Value of wD/2V:

+W=0.5*N*D/1000/V*2*%pi/60

+Red=V*D/v;

+//For this value, CL is obtained as:

+CL=0.3;

+//Aerodynamic lift(in N):

+FL=%pi/8*CL*(D/1000)^2*d*V^2;

+//Radius of curvature of the path in the vertical plane(in m) with topspin:

+Rts=V^2/(g+FL/(m/1000));

+//Radius of curvature without topspin(in m):

+Rwts=V^2/g;

+printf("\n\nRESULTS\n\n")

+printf("\n\nAerodynamic lift acting on the ball:%.3f N\n\n",FL)

+printf("\n\nRadius of curvature of the path when ball has topspin:%.3f m\n\n",Rts)

+printf("\n\nRadius of curvature of the path when ball has topspin: %.3f m\n\n",Rwts)