initial commit / add all books

3 files changed, 75 insertions, 0 deletions

diff --git a/812/CH7/EX7.04/7_04.sce b/812/CH7/EX7.04/7_04.sce
new file mode 100755
index 000000000..b6e3ab182
--- /dev/null
+++ b/812/CH7/EX7.04/7_04.sce
@@ -0,0 +1,18 @@
+//speed and force//

+pathname=get_absolute_file_path('7.04.sce')

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

+exec(filename)

+//Velocity of prototype in ft/sec

+Vp1=Vp*6080/3600

+//Reynolds number of prototype:

+Rep=Vp1*Dp/vp

+//Rep=Rem

+//Therefore:

+Rem=Rep;

+//Velocity of air for wind tunnel(in ft/sec):

+Vm=Rem*vm/(Dm/12)

+//Drag force on prototype(in lbf):

+Fp=Fm*(dp/dm)*(Vp1/Vm)^2*(Dp/(Dm/12))^2

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

+printf("\n\nTest speed in air: %.3f ft/sec\n\n",Vm)

+printf("\n\nDrag force on prototype: %.3f lbf\n\n",Fp)

diff --git a/812/CH7/EX7.05/7_05.sce b/812/CH7/EX7.05/7_05.sce
new file mode 100755
index 000000000..e0e2e8ef6
--- /dev/null
+++ b/812/CH7/EX7.05/7_05.sce
@@ -0,0 +1,29 @@
+//speed force and power//

+pathname=get_absolute_file_path('7.05.sce')

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

+exec(filename)

+//Width of the model(in m):

+wm=S*wp*0.3048

+//Area of model (in m^2):

+Am=S^2*Ap*0.305^2

+[m n]=size(V)

+i=1:n

+//Aerodynamic drag coefficient():

+Cd=2.*Fd(i)/d./(V(i))^2/0.0305

+//Reynolds number:

+Re=V(i)*wm/v

+plot(Re,Cd);

+a=gca()

+a.data_bounds=[100000,0.4;500000,0.6]

+xtitle('Aerodynamic drag coefficient vs drag force','Reynolds number','Model Drag Coeff.')

+//It is seen that drag coefficient becomes constant at CD=0.46above Re=4*10^5 at which speed of air is 40m/s

+CDc=0.46;

+Va=40;

+//Drag force (in N):

+FDp=CDc/2*d*(Vp*5/18)^2*Ap*0.305^2

+//Power required to pull prototype at 100 kmph(in W)

+Pp=FDp*Vp*5/18

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

+printf("\n\nSpeed above which Cd is constant: %.3f m/sec\n\n",Va)

+printf("\n\nDrag Force: %.3f kN\n\n",FDp/1000)

+printf("\n\nPower required to pull prototype at 100 kmph: %.3f kW\n\n",Pp/1000)

diff --git a/812/CH7/EX7.06/7_06.sce b/812/CH7/EX7.06/7_06.sce
new file mode 100755
index 000000000..ea7c51970
--- /dev/null
+++ b/812/CH7/EX7.06/7_06.sce
@@ -0,0 +1,28 @@
+//power and speed//

+pathname=get_absolute_file_path('7.06.sce')

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

+exec(filename)

+//The same pump is used for both the conditions.Hence:

+D2=D1;

+//The same water is used for both the conditions. Hence:

+d2=d1;

+//Flow rate at condition 2(in gpm):

+Q2=Q1*N2/N1*(D2/D1)^3

+//Head at condition 1(in ft):

+H1=(N1*sqrt(Q1)/Nscu1)^(4/3)

+//Head at condition 1(in ft):

+H2=H1*(N2/N1)^2*(D2/D1)^2

+//Pump output power at condition 1(in hp):

+P1=d1*g*Q1*H1/7.48/60/550

+//Pump output power at condition 2(in hp):

+P2=P1*(d2/d1)*(N2/N1)^3*(D2/D1)^5

+//Required input power(in hp):

+Pin=P2/Effp

+//Specific speed at condition 2:

+Nscu2=N2*sqrt(Q2)/H2^(3/4)

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

+printf("\n\nVolume flow rate at condition 2: %.3f gpm\n\n\n",Q2)

+printf("\n\nHead at condition: %.3f ft\n\n\n",H2)

+printf("\n\nPump output power at condition: %.3f hp\n\n\n",P2)

+printf("\n\nRequired input power: %.3f hp\n\n\n",Pin)

+printf("\n\nSpecific speed at condition 2: %.3f\n\n\n",Nscu2)