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 --- 746/CH10/EX10.01/10_01.sce | 14 +++++++++++ 746/CH10/EX10.02/10_02.sce | 24 ++++++++++++++++++ 746/CH10/EX10.03/10_03.sce | 34 +++++++++++++++++++++++++ 746/CH10/EX10.06/10_06.sce | 28 +++++++++++++++++++++ 746/CH10/EX10.07/10_07.sce | 22 ++++++++++++++++ 746/CH10/EX10.08/10_08.sce | 62 ++++++++++++++++++++++++++++++++++++++++++++++ 746/CH10/EX10.08/10_08.txt | 14 +++++++++++ 746/CH10/EX10.11/10_11.sce | 32 ++++++++++++++++++++++++ 746/CH10/EX10.12/10_12.sce | 29 ++++++++++++++++++++++ 746/CH10/EX10.14/10_14.sce | 16 ++++++++++++ 746/CH10/EX10.16/10_16.sce | 20 +++++++++++++++ 11 files changed, 295 insertions(+) create mode 100755 746/CH10/EX10.01/10_01.sce create mode 100755 746/CH10/EX10.02/10_02.sce create mode 100755 746/CH10/EX10.03/10_03.sce create mode 100755 746/CH10/EX10.06/10_06.sce create mode 100755 746/CH10/EX10.07/10_07.sce create mode 100755 746/CH10/EX10.08/10_08.sce create mode 100755 746/CH10/EX10.08/10_08.txt create mode 100755 746/CH10/EX10.11/10_11.sce create mode 100755 746/CH10/EX10.12/10_12.sce create mode 100755 746/CH10/EX10.14/10_14.sce create mode 100755 746/CH10/EX10.16/10_16.sce (limited to '746/CH10') 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) -- cgit