initial commit / add all books

11 files changed, 295 insertions, 0 deletions

diff --git a/746/CH10/EX10.01/10_01.sce b/746/CH10/EX10.01/10_01.sce
new file mode 100755
index 000000000..5b2e3888e
--- /dev/null
+++ b/746/CH10/EX10.01/10_01.sce
@@ -0,0 +1,14 @@
+//input and power//

+pathname=get_absolute_file_path('10.01.sce')

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

+exec(filename)

+//Impeller exit width b2(in feet):

+b2=Q*12/(2*%pi*R2*Vrb2*7.48*60)

+//Torque of the Shaft, Tshaft(in ft-lbf):

+Tshaft=w*R2^2*p*Q*2*%pi/3600/7.48/144

+//Power, Wm(in hp):

+Wm=w*Tshaft*2*%pi/60/550

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

+printf("\n\nImpeller exit width: %.3f feet\n\n",b2)

+printf("\n\Torque input: %.3f ft-lbf\n\n",Tshaft)

+printf("\n\nPower: %.3f hp\n\n",Wm)

diff --git a/746/CH10/EX10.02/10_02.sce b/746/CH10/EX10.02/10_02.sce
new file mode 100755
index 000000000..214c9e24e
--- /dev/null
+++ b/746/CH10/EX10.02/10_02.sce
@@ -0,0 +1,24 @@
+//volume and power//

+pathname=get_absolute_file_path('10.02.sce')

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

+exec(filename)

+U=0.5*(Dh+Dt)/2*1200*2*%pi/60

+k=tand(alpha1)+cotd(betta1)

+Vn1=U/k

+V1=Vn1/cosd(alpha1)

+Vt1=V1*sind(alpha1)

+Vrb1=Vn1/sind(betta1)

+//Volume flow rate (in m^3/sec):

+Q=%pi/4*Vn1*(Dt^2-Dh^2)

+k=(U-Vn1*cotd(betta2))/Vn1

+alpha2= atand(k)

+V2=Vn1/cosd(alpha2)

+Vt2=V2*sind(alpha2)

+//Rotor Torque (in N-m):

+Tz=p*Q*(Dh+Dt)/4*(Vt2-Vt1)

+//Power required (in W):

+Wm=w*2*%pi/60*Tz

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

+printf("\n\nVolume flow rate: %.3f m^3/sec\n\n",Q)

+printf("\n\nRotor Torque: %.3f N-m\n\n",Tz)

+printf("\n\nPower required: %.3f W\n\n",Wm)

diff --git a/746/CH10/EX10.03/10_03.sce b/746/CH10/EX10.03/10_03.sce
new file mode 100755
index 000000000..bcaf8b526
--- /dev/null
+++ b/746/CH10/EX10.03/10_03.sce
@@ -0,0 +1,34 @@
+//Pump Power//

+pathname=get_absolute_file_path('10.03.sce')

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

+exec(filename)

+[nQ mQ]= size(Q);

+[nps mps]=size(ps);

+[npd mpd]= size(pd);

+[nI mI]= size(I);

+//Correct measured static pressures to he pump centreline p1, p2(in psig):

+ j=1:mps;

+  p1=ps(j)+px*g*zs/144

+   j=1:mpd;

+  p2=pd(j)+px*g*zd/144

+  //The value of Pump head(in feet):

+  j=1:mps;

+  Hp=(p2(j)-p1(j))/(px*g)*144

+  //Values of Hydraulic Power delivered(in hp):

+   j=1:mps;

+   Wh=Q(j).*(p2(j)-p1(j))/7.48/60*144/550

+   //Values of motor power output(in hp):

+  j=1:mI;

+  Pin=Effm*sqrt(3)*PF*E*I(j)/746

+  //Values of Pump Efficiecy:

+   j=1:mI;

+  Effp= Wh(j)./Pin(j)*100

+  //Plotting pump characteristics:

+   plot(Q,Hp,"-o")

+  plot(Q,Pin,"-+")

+  plot(Q,Effp,"-*")

+  xtitle('Pump Characteristics','Volume flow rate(in gpm)',['Pump Efficincy(%)    ','   Pump Head(in feet)   ','    Pump Power input(in hp)   '])

+  legend('Hp','Pin','Effp')

+

+

+

diff --git a/746/CH10/EX10.06/10_06.sce b/746/CH10/EX10.06/10_06.sce
new file mode 100755
index 000000000..e9eea036a
--- /dev/null
+++ b/746/CH10/EX10.06/10_06.sce
@@ -0,0 +1,28 @@
+//Specific and relation//

+pathname=get_absolute_file_path('10.06.sce')

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

+exec(filename)

+//Specific speed in Us customary units:

+Nscu=N*Qus^0.5/Hus^0.75

+//Conversion to SI units:

+w=1170*2*%pi/60;

+Qsi=Qus/7.48/60*0.305^3;

+Hsi=Hus*0.305;

+//Energy per unit mass is:

+h=g*Hsi;

+//Specific speed in SI units:

+Nssi=w*Qsi^0.5/h^0.75

+//Conversion to hertz:

+whz=N/60;

+//Specific speed in European units:

+Nseu=whz*Qsi^0.5/65.5^0.75

+//Relation between specific speeds in Us customary units and European units:

+Conversionfactor1=Nscu/Nseu

+//Relation between specific speeds in Us customary units and SI units:

+Conversionfactor2=Nscu/Nssi

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

+printf("\n\nSpecific speed in US customary units: %.3f \n\n",Nscu)

+printf("\n\nSpecific speed in SI units: %.3f \n\n",Nssi)

+printf("\n\nSpecific speedin European units: %.3f \n\n",Nseu)

+printf("\n\nRelation between specific speeds in Us customary units and European units: %.3f \n\n",Conversionfactor1)

+printf("\n\nRelation between specific speeds in Us customary units and SI units: %.3f \n\n",Conversionfactor2)

diff --git a/746/CH10/EX10.07/10_07.sce b/746/CH10/EX10.07/10_07.sce
new file mode 100755
index 000000000..5202ddee1
--- /dev/null
+++ b/746/CH10/EX10.07/10_07.sce
@@ -0,0 +1,22 @@
+//Comparison of head//

+pathname=get_absolute_file_path('10.07.sce')

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

+exec(filename)

+//Volume flow rate(in gpm) at shut off condition for N2:

+Q2so=N2/N1*Q1so

+//Volume flow(in gpm) rate at best efficiency for N2:

+Q2be=N2/N1*Q1be

+//Relation between pump heads:

+head_relation=(N2/N1)^2

+//Head(in feet) at shut off condition for N2:

+H2so=(N2/N1)^2*H1so

+//Head(in feet) at best efficiency condition for N2:

+H2be=(N2/N1)^2*H1be

+Q1=[Q1so Q1be];

+Q2=[Q2so Q2be];

+H1=[H1so H1be];

+H2=[H2so H2be];

+plot(Q1,H1,"-o")

+plot(Q2,H2,"-*")

+xtitle('Comparison of head for both conditions','Volume Flow Rate','Head') 

+legend('1170','1750')

diff --git a/746/CH10/EX10.08/10_08.sce b/746/CH10/EX10.08/10_08.sce
new file mode 100755
index 000000000..30e6d3aac
--- /dev/null
+++ b/746/CH10/EX10.08/10_08.sce
@@ -0,0 +1,62 @@
+//NPSHA and NPSHR// 

+pathname=get_absolute_file_path('10.08.sce')

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

+exec(filename)

+//Diameter of pipe (in feet):

+Df= Di/12

+//Area of crossection of pipe(in ft^2):

+A=%pi/4*Df^2

+//Velocity of flow(in ft/sec):

+V=Q/7.48/A/60

+//For water at T=80F,viscosity=0.927e-5 ft^2/sec, Reynolds number:

+Re=V*Df/v

+//Friction loss Coefficient for this value of Re:

+f=0.0237;

+//For cast iron, roughness(in feet):

+e=0.00085

+//e/D is:

+e/Df

+//Total head loss(in feet):

+HL=K+f*(SE+OGV)+f*(L/Df)+1

+//The heads are(in feet):

+H1=patm*144/(p*g)

+Vh=V^2/2/g

+//Suction head(in feet):

+Hs=H1+h-HL*Vh

+//NPSHA(in feet):

+NPSHA=Hs+Vh-Hv1

+//For a flow rate of 1000 gpm,NPSHR(in feet) for water at 80 F

+NPSHR=10 

+//PLOTTING NPSHA AND NPSHR VERSUS VOLUME FLOW RATE:

+//For 80 F

+Qp=0:100:1500;

+[nQp mQp]=size(Qp);

+ for j=1:mQp;

+ Vp(j)=Qp(j)/(7.48*A*60);

+ Vhp(j)=(Vp(j))^2/2/g;

+ Hs(j)=H1+h-HL*Vhp(j);

+ end

+

+for j=1:mQp;

+  NPSHAp1(j)=Hs(j)+(Vhp(j))-Hv1;

+end

+

+plot(Qp,NPSHAp1,"-+")

+plot(Qh,NPSHRp,"-o")

+xtitle('Suction head vs Flow rate','Volume flow rate(gpm)','Suction Head(feet)');

+printf("\n\nType (Resume) to continue or (abort) to end\n\n")

+legend('NPSHA','NPSHR')

+pause 

+clf

+

+//For 180 F

+for j=1:mQp;

+  NPSHAp2(j)=Hs(j)+(Vhp(j))-Hv2;

+end

+plot(Qp,NPSHAp2,"-+")

+plot(Qh,NPSHRp,"-o")

+xtitle('Suction head vs Flow rate','Volume flow rate(gpm)','Suction Head(feet)');

+legend('NPSHA','NPSHR')

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

+printf("\n\nNPSHA at Q=1000 gpm of water at 80 F: %.2f ft\n\n",NPSHA)

+printf("\n\nNPSHR at Q=1000 gpm of water at 80 F: %.1f ft\n\n",NPSHR)

diff --git a/746/CH10/EX10.08/10_08.txt b/746/CH10/EX10.08/10_08.txt
new file mode 100755
index 000000000..a971c82e7
--- /dev/null
+++ b/746/CH10/EX10.08/10_08.txt
@@ -0,0 +1,14 @@
+j=0:100:1500;

+NPSHAp2=Hs+(V(j))^2/2/g-Hv2;

+

+

+

+

+

+

+

+

+pause

+clf

+plot(Qp,NPSHAp2)

+xtitle('Volume flow rate(gpm)','Suction Head(feet)','NPSHA when water is at 180 F');
\ No newline at end of file
diff --git a/746/CH10/EX10.11/10_11.sce b/746/CH10/EX10.11/10_11.sce
new file mode 100755
index 000000000..a485336ed
--- /dev/null
+++ b/746/CH10/EX10.11/10_11.sce
@@ -0,0 +1,32 @@
+//Performance curves//

+pathname=get_absolute_file_path('10.11.sce')

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

+exec(filename)

+[nQ mQ]= size(Q1);

+[np mp]= size(p1);

+[nP mP]= size(P1);

+//Volume flow rate for fan 2(in cfm):

+j=1:mQ;

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

+//Pressure values for fan 2(in inches of H2O):

+j=1:mp;

+p2=p1(j)*(d2/d1)*((N2/N1)^2)*((D2/D1)^2)

+//Power values for fan 2(in hp):

+j=1:mP;

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

+plot(Q2,p2)

+xtitle('Performance curves','Volume flow rate(in cfm)','Pressure head(in inches of water)')

+printf("\n\nType (resume) to continue or (abort) to exit\n\n")

+pause 

+clf

+plot(Q2,P2)

+xtitle('Performance curves','Volume flow rate(in cfm)','Power(in hp)')

+printf("\n\nType (resume) to continue or (abort) to exit\n\n")

+pause

+clf

+plot(Q2,Eff)

+xtitle('Performance curves','Volume flow rate(in cfm)','Eficiency(in percentage)')

+//Specific speed of fan(in US customary units) at operating point:

+Nscu= 1150*110000^0.50*0.045^0.75/7.4^0.75

+//Specific speed of fan (in SI units) at operating point:

+Nssi=120*3110^0.5*0.721^0.75/1.84e3^0.75

diff --git a/746/CH10/EX10.12/10_12.sce b/746/CH10/EX10.12/10_12.sce
new file mode 100755
index 000000000..6dfe1e933
--- /dev/null
+++ b/746/CH10/EX10.12/10_12.sce
@@ -0,0 +1,29 @@
+//Power required//

+pathname=get_absolute_file_path('10.12.sce')

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

+exec(filename)

+//From given graph, for maximum delivery condition, Q=48.5gpm.

+//Volume of oil per revolution delivered by the pump(in in^3/rev):

+vc=Qe/N*231

+//Volumetric Effciency of pump at max flow:

+Effv=vc/va

+//Operating point of the pump is found to be at 1500 psig,Q=46.5gpm

+//Power delivered by the fluid(in hp):

+Pf=Qo*po1/7.48/60*144/550

+//Input power(in hp):

+Pi=Pf/Effp

+//The power delivered to the load(in hp):

+Pl=Q*(po1)/7.48/60*144/550

+//Power dissipated by throttling(in hp):

+Pd=Pf-Pl

+//The dissipation with the variable displacement pump(in hp):

+Pvd=Q*(po2-po1)/7.48/60*144/550

+//Power required for te load sensing pump if pump pressure is 100psi above that required by the load(in hp):

+Pls=Q*100/7.48/60*144/550

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

+printf("\n\nVolume of oil per revolution delivered by the pump: %.3f in^3/rev\n\n",vc)

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

+printf("\n\nPower deliverd to the load: %.3f hp\n\n",Pl)

+printf("\n\nPower dissipated by throttling: %.3f hp\n\n",Pd)

+printf("\n\nThe dissipation with the variable displacement pump: %.3f hp\n\n",Pvd)

+printf("\n\nPower required for te load sensing pump if pump pressure is 100psi above that required by the load: %.3f hp\n\n",Pls)

diff --git a/746/CH10/EX10.14/10_14.sce b/746/CH10/EX10.14/10_14.sce
new file mode 100755
index 000000000..efa884132
--- /dev/null
+++ b/746/CH10/EX10.14/10_14.sce
@@ -0,0 +1,16 @@
+//propeller//

+pathname=get_absolute_file_path('10.14.sce')

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

+exec(filename)

+//Propeller Thrust(in MN) :

+Ft=P/V

+//Required power input  to the propeller(in MW):

+Pin=P/Eff

+//Calculating value of D(in m):

+nD=V/J

+D=(Ft*10^6/p/(nD)^2/Cf)^0.5

+//Operating speed (in rpm) is given by:

+n=nD/D*60

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

+printf("\n\nDiameter of the single propeller required to pwer the ship:%.3f m\n\n",D)

+printf("\n\nOperating speed of the propeller: %.3f rpm\n\n",n)

diff --git a/746/CH10/EX10.16/10_16.sce b/746/CH10/EX10.16/10_16.sce
new file mode 100755
index 000000000..69a3995a4
--- /dev/null
+++ b/746/CH10/EX10.16/10_16.sce
@@ -0,0 +1,20 @@
+//Actual//

+pathname=get_absolute_file_path('10.16.sce')

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

+exec(filename)

+//Tip speed ratio of windmill:

+X=N*2*%pi/60*D/2/(V*5/18)

+//Accounting for whirl,max attainable efficiency is:

+Efw=0.53;

+//Kinetic energy flux(in W) is given by:

+KEF=0.5*p*(V*5/18)^3*%pi*(D/2)^2

+//Actual Efficiency:

+Effa=Po/KEF

+//The maximum possible thrust occurs for an interference factor of:

+amax=0.5;

+//Thrust(in W):

+Kx=p*(V*5/18)^2*%pi*(D/2)^2*2*amax*(1-amax)

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

+printf("\n\nTip speed ratio of windmill:%.3f\n\n",X)

+printf("\n\nActual Efficiency: %.3f\n\n",Effa)

+printf("\n\nActual Thrust: %.3f W\n\n",Kx)