diff options
Diffstat (limited to '746')
161 files changed, 2675 insertions, 0 deletions
diff --git a/746/CH1/EX1.01/1_01.sce b/746/CH1/EX1.01/1_01.sce new file mode 100755 index 000000000..576606142 --- /dev/null +++ b/746/CH1/EX1.01/1_01.sce @@ -0,0 +1,7 @@ +//Heat addition//
+filename=pathname+filesep()+'1.01-data.sci'
+exec(filename)
+//Heat added during the process(in kJ):
+Q12=m*cp*(T2-T1)
+printf("\n\nRESULTS\n\n")
+printf("\n\nHeat added during the process: %f kJ\n\n",Q12/1000)
diff --git a/746/CH1/EX1.02/1_02.sce b/746/CH1/EX1.02/1_02.sce new file mode 100755 index 000000000..5a88f0378 --- /dev/null +++ b/746/CH1/EX1.02/1_02.sce @@ -0,0 +1,13 @@ +//speed and actual speed//
+pathname=get_absolute_file_path('1.02.sce')
+filename=pathname+filesep()+'1.02-data.sci'
+exec(filename)
+//Speed at which the ball hits the ground(in m/sec):
+V=sqrt(m*g/k*(1-%e^(2*k/m*(-y0))))
+//Terminal speed(in m/sec):
+Vt=sqrt(m*g/k)
+//Ratio of actual speed to the terminal speed:
+r=V/Vt;
+printf("\n\nRESULTS\n\n")
+printf("\n\nSpeed at which the ball hits he ground: %f m/sec\n\n",V)
+printf("\n\nRatio of actual speed to the terminal speed: %f\n\n",r)
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)
diff --git a/746/CH11/EX11.01/11_01.sce b/746/CH11/EX11.01/11_01.sce new file mode 100755 index 000000000..073540f56 --- /dev/null +++ b/746/CH11/EX11.01/11_01.sce @@ -0,0 +1,19 @@ +//Change//
+pathname=get_absolute_file_path('11.01.sce')
+filename=pathname+filesep()+'11.01-data.sci'
+exec(filename)
+//Density of air at entry:
+d1=p1*10^3/R/T1
+//Area(in m^2):
+A=m/d1/V1
+//Change in enthalpy of air(in kJ/kg):
+dh=cp*(T2-T1)
+//Change in internal energy of air(in kJ/kg):
+du=cv*(T2-T1)
+//Change in entropy(in kJ/(kg-K)):
+ds=cp*log(T2/T1)-R/1000*log(p2/p1)
+printf("\n\nRESULTS\n\n")
+printf("\n\nDuct Area: %.3f m^2\n\n",A)
+printf("\n\nChange in enthalpy of air: %.3f kJ/kg\n\n",dh)
+printf("\n\nChange in internal energy of air:%.3f kJ/kg\n\n",du)
+printf("\n\nChange in entropy: %.3f kg-K\n\n",ds)
diff --git a/746/CH11/EX11.03/11_03.sce b/746/CH11/EX11.03/11_03.sce new file mode 100755 index 000000000..cc4eee914 --- /dev/null +++ b/746/CH11/EX11.03/11_03.sce @@ -0,0 +1,19 @@ +//Speed of sound// +pathname=get_absolute_file_path('11.03.sce') +filename=pathname+filesep()+'11.03-data.sci' +exec(filename) +//Values of altitude(in m): +Al=0:1000:15000 +[nAl mAl]=size(Al); +//Values of temperature at given altitudes(in K): +T=[288.2 281.7 275.2 268.7 262.2 255.7 249.2 242.7 236.2 229.7 223.3 216.8 216.7 216.7 216.7 216.7]; +[nT mT]=size(T); +//Values of speed of sound at these altitudes(in m/sec): +j=1:mT; +c=sqrt(k*R*T(j)) +//Speed of sound at sea level(in m/sec): +c1=sqrt(k*R*T(1)) +plot(c,Al) +xtitle('Variation of sound speed with altitude','Speed of sound(m/sec)','Altitude(m)') +printf("\n\nRESULTS\n\n") +printf("\n\nSpeed of sound at sea level: %.3f m/sec\n\n",c1) diff --git a/746/CH11/EX11.04/11_04.sce b/746/CH11/EX11.04/11_04.sce new file mode 100755 index 000000000..350f4dc6e --- /dev/null +++ b/746/CH11/EX11.04/11_04.sce @@ -0,0 +1,22 @@ +//pressure and change//
+pathname=get_absolute_file_path('11.04.sce')
+filename=pathname+filesep()+'11.04-data.sci'
+exec(filename)
+//Mach number at entry:
+M1=V1/sqrt(k*R*T1)
+//Stagnation pressure at entry(in kPa):
+p01=p1*(1+(k-1)/2*M1^2)^(k/(k-1))
+//Stagnation temperature at entry(in K):
+T01=T1*(1+(k-1)/2*M1^2)
+//Static pressure at exit(in kPa):
+p2=p02/(1+(k-1)/2*M2^2)^(k/(k-1))
+//Temperature at exit(in K):
+T2=T02/(1+(k-1)/2*M2^2)
+//Change in entropy(in kJ/kg-K):
+ds=cp*log(T2/T1)-R/1000*log(p2/p1)
+printf("\n\nRESULTS\n\n")
+printf("\n\nStagnation pressure at entry: %.3f kPa\n\n",p01)
+printf("\n\nStagnation temperature at entry: %.3f K\n\n",T01)
+printf("\n\nStatic pressure at exit: %.3f kPa\n\n",p2)
+printf("\n\nTemperature at exit: %.3f K\n\n",T2)
+printf("\n\nChange in entropy: %.3f kJ/kg-K\n\n",ds)
diff --git a/746/CH12/EX12.01/12_01.sce b/746/CH12/EX12.01/12_01.sce new file mode 100755 index 000000000..ede9ac151 --- /dev/null +++ b/746/CH12/EX12.01/12_01.sce @@ -0,0 +1,33 @@ +//pressure and area//
+pathname=get_absolute_file_path('12.01.sce')
+filename=pathname+filesep()+'12.01-data.sci'
+exec(filename)
+//Here the stagnation quantities are constant.
+// Stagnation temperature(in K):
+T0=T1*(1+(k-1)/2*M1^2)
+//Stagnation pressure(in kPa):
+p0=p1*((1+(k-1)/2*M1^2)^(k/(k-1)))
+//Finding T2/T1:
+T=t2/t1
+//Temperature at exit(in K):
+T2=T*T1
+//Finding p2/p1:
+P=P2/P1
+//Pressure at exit(in kPa):
+p2=P2*p1
+//Density of air at exit(in kg/m^3):
+d2=p2*10^3/R/T2
+//Velocity of air at exit(in m/sec):
+V2=M2*sqrt(k*R*T2)
+//Finding A2/A1:
+a=a2/a1
+//Area at exit(in m^2):
+A2=a*A1
+printf("\n\nRESULTS\n\n")
+printf("\n\nStagnation temperature: %.3f K\n\n",T0)
+printf("\n\nStagantion pressure: %.3f kPa\n\n",p0)
+printf("\n\nTemperature a exit %.3f K\n\n",T2)
+printf("\n\nPressure at exit: %.3f kPa\n\n",p2)
+printf("\n\nDensity of air at exit: %.3f kg/m^3\n\n",d2)
+printf("\n\nVelocity of air at exit: %.3f m/sec\n\n",V2)
+printf("\n\nArea at exit: %.3f \n\n",A2)
diff --git a/746/CH12/EX12.02/12_02.sce b/746/CH12/EX12.02/12_02.sce new file mode 100755 index 000000000..9410fab06 --- /dev/null +++ b/746/CH12/EX12.02/12_02.sce @@ -0,0 +1,27 @@ +//Mass flow//
+pathname=get_absolute_file_path('12.02.sce')
+filename=pathname+filesep()+'12.02-data.sci'
+exec(filename)
+//Checking for chocking:
+c=pb/p0;
+if(c<=0.528)
+ //choked
+else
+ //Not choked
+ //Therefore pressure at exit = back pressure
+ pe=pb;
+ //Mach number at exit:
+ Me=(((p0/pe)^((k-1)/k)-1)*(2/(k-1)))^0.5
+ //Temperature at exit(in K):
+ Te=T0/(1+(k-1)/2*Me^2)
+ //Velocity at exit(in m/sec):
+ Ve=Me*sqrt(k*R*Te)
+ //Density at exit(in kg/m^3):
+ de=pe*10^3/R/Te
+ //Mass flow rate of air(kg/sec):
+ m=de*Ve*Ae
+end;
+printf("\n\nRESULTS\n\n")
+printf("\n\nMach number at exit: %.3f\n\n",Me)
+printf("\n\nMass flow rate of air: %.3f kg/sec\n\n",m)
+
diff --git a/746/CH12/EX12.03/12_03.sce b/746/CH12/EX12.03/12_03.sce new file mode 100755 index 000000000..23bf28dac --- /dev/null +++ b/746/CH12/EX12.03/12_03.sce @@ -0,0 +1,29 @@ +//mass and area//
+pathname=get_absolute_file_path('12.03.sce')
+filename=pathname+filesep()+'12.03-data.sci'
+exec(filename)
+//Saturation pressure(in psia):
+p0=p1*(1+(k-1)/2*M1^2)^(k/(k-1))
+//Checking for choking:
+x=pb/p0;
+if(x>0.528)
+ //Not choked
+else
+ //choked
+end
+//As there is choking:
+Mt=1;
+//Velocity at entry:
+V1=M1*sqrt(k*R*(T1+460)*32.2)
+//Density at the entry(in lbm/ft^3):
+d1=p1/(R*(T1+460))*144
+//Mass flow rate(in lbm/sec):
+m=d1*V1*A1
+//Finding the valueof A1/A*;
+A=1/M1*((1+(k-1)/2*M1^2)/(1+(k-1)/2))^((k+1)/(2*(k-1)))
+//For choked flow, At=A*
+At=A1/A
+printf("\n\nRESULTS\n\n")
+printf("\n\nMach number at throat: %.3f\n\n",Mt)
+printf("\n\nMass flow rate: %.3f lbm/sec\n\n",m)
+printf("\n\nArea at throat: %.3f ft^2\n\n",At)
diff --git a/746/CH12/EX12.04/12_04.sce b/746/CH12/EX12.04/12_04.sce new file mode 100755 index 000000000..5839d27ec --- /dev/null +++ b/746/CH12/EX12.04/12_04.sce @@ -0,0 +1,32 @@ +//throat//
+pathname=get_absolute_file_path('12.04.sce')
+filename=pathname+filesep()+'12.04-data.sci'
+exec(filename)
+//Temperature at the throat(in K):
+Tt=T0/(1+(k-1)/2*Mt^2)
+//Pressure at throat(in kPa):
+pt=p0*(Tt/T0)^(k/(k-1))
+//Density at throat(in kg/m^3):
+dt=pt*1000/R/Tt
+//Velocity at the throat(in m/s):
+Vt=Mt*sqrt(k*R*Tt)
+//Value of At/A*:
+Ax=1/Mt*((1+(k-1)/2*Mt^2)/(1+(k-1)/2))^((k+1)/(2*(k-1)))
+//Stagnation properties are constant
+//As a result pressure at exit,
+pe=pb;
+//The Mach number at the exit is therefore given by
+Me=sqrt(((p0/pe)^((k-1)/k)-1)*2/(k-1))
+//Calculating the value of Ae/A*:
+Ay=1/Me*((1+(k-1)/2*Me^2)/(1+(k-1)/2))^((k+1)/(2*(k-1)))
+//Value of A*(in m^2):
+A_star=Ae/Ay
+//Area at throat(in m^2):
+At=Ax*A_star
+printf("\n\nRESULTS\n\n")
+printf("\n\nTemperature at the throat: %.3f K\n\n",Tt)
+printf("\n\nPressure at throat: %.3f kPa\n\n",pt)
+printf("\n\nDensity at throat: %.3f kg/m^3\n\n",dt)
+printf("\n\nVelocity at the throat: %.3f m/sec\n\n",Vt)
+printf("\n\nMach number at the exit: %.3f\n\n",Me)
+printf("\n\nArea at throat: %.3f m^2\n\n",At)
diff --git a/746/CH12/EX12.05/12_05.sce b/746/CH12/EX12.05/12_05.sce new file mode 100755 index 000000000..22d125b1b --- /dev/null +++ b/746/CH12/EX12.05/12_05.sce @@ -0,0 +1,13 @@ +//number and flow//
+pathname=get_absolute_file_path('12.05.sce')
+filename=pathname+filesep()+'12.05-data.sci'
+exec(filename)
+//Mach number at the exit:
+Me=sqrt(((p0/pe)^((k-1)/k)-1)*2/(k-1))
+//Temperature at exit(in K):
+Te=T0/(1+(k-1)/2*Me^2)
+//Mass flow rate(in kg/s):
+m=pe*1000*Me*sqrt(k/R/Te)*Ae
+printf("\n\nRESULTS\n\n")
+printf("\n\nMach number at the exit: %.3f \n\n",Me)
+printf("\n\nMass flow rate: %.3f kg/sec\n\n",m)
diff --git a/746/CH12/EX12.06/12_06.sce b/746/CH12/EX12.06/12_06.sce new file mode 100755 index 000000000..f7c00e5cd --- /dev/null +++ b/746/CH12/EX12.06/12_06.sce @@ -0,0 +1,32 @@ +//mass and volume//
+pathname=get_absolute_file_path('12.06.sce')
+filename=pathname+filesep()+'12.06-data.sci'
+exec(filename)
+//Mach umber at section 1:
+M1=sqrt((2/(k-1)*((p0/p1)^((k-1)/k)-1)))
+//Temperature at section 1(in K):
+T1=T0/(1+(k-1)/2*M1^2)
+//Density at section 1(in kg/m^3):
+d1=p1*1000/R/T1
+//Velocity at section1(in m/sec):
+V1=M1*sqrt(k*R*T1)
+//Area at section 1(in m^2):
+A1=%pi/4*D^2
+//Mass flow rate(in kg/sec):
+m=d1*A1*V1
+//Mach number at section 2:
+M2=sqrt((2/(k-1))*((T0/T2)-1))
+//Velocity at section 2(in m/sec):
+V2=M2*sqrt(k*R*T2)
+//Density at section 2(in kg/m^3):
+d2=d1*V1/V2
+//Pressure at section 2(in kPa):
+p2=d2/1000*R*T2
+//Stagnation pressure at section 2(in kPa):
+p02=p2*(1+(k-1)/2*M2^2)^(k/(k-1))
+//Force exerted on control volume by duct wall(in N):
+F=(p2-p1)*1000*A1+m*(V2-V1)
+printf("\n\nRESULTS\n\n")
+printf("\n\nMass flow rate: %.3f kg/sec\n\n",m)
+printf("\n\nLocal isentropic stagnation pressure at section 2:%.3f kPa\n\n",p02)
+printf("\n\nForce exerted on control volume by duct wall:%.3f N\n\n",F)
diff --git a/746/CH12/EX12.07/12_07.sce b/746/CH12/EX12.07/12_07.sce new file mode 100755 index 000000000..44177a095 --- /dev/null +++ b/746/CH12/EX12.07/12_07.sce @@ -0,0 +1,34 @@ +//length//
+pathname=get_absolute_file_path('12.07.sce')
+filename=pathname+filesep()+'12.07-data.sci'
+exec(filename)
+//Mach number at section 1:
+M1= sqrt(2/(k-1)*((p0/(p0+p1))^((k-1)/k)-1))
+//Temperature at section 1(in K):
+T1=T0/(1+(k-1)/2*(M1)^2)
+V1=M1*sqrt(k*R*T1)
+//Pressure at section 1(in kPa):
+p1=g*dHg*(760-18.9)*10^-3
+//Density at section 1(in kg/m^3):
+d1=p1/R/T1
+//At M1=0.190,
+//(p/p*)1:
+P1=5.745
+// (fLmax/Dh)1:
+F1=16.38
+//Value of L13(in m):
+L13=F1*D/f
+//Value of (p/p*)2:
+P2=p2/p1*P1
+//For this value, Value of M2 is obtained as 0.4
+M2=0.4;
+//For M=0.4, fLmX/D=2.309
+F2=2.309
+//Value of L23(in m):
+L23=F2*D/f
+//Length of duct between section 1 and 2(in m):
+L12=L13-L23
+printf("\n\nRESULTS\n\n")
+printf("\n\nLength of duct required for choking from section 1: %3f m\n\n",L13)
+printf("\n\nMach number section 2: %.3f \n\n",M2)
+printf("\n\Length of duct between section 1 and 2: %.3f m\n\n",L12)
diff --git a/746/CH12/EX12.08/12_08.sce b/746/CH12/EX12.08/12_08.sce new file mode 100755 index 000000000..cd86a9749 --- /dev/null +++ b/746/CH12/EX12.08/12_08.sce @@ -0,0 +1,34 @@ +//velocity and entropy//
+pathname=get_absolute_file_path('12.08.sce')
+filename=pathname+filesep()+'12.08-data.sci'
+exec(filename)
+//Density at section 1(in lbm/ft^3):
+d1=p1*144/R/T1
+//Velocity at section 2(in ft/sec):
+V2=(p1-p2)*144/d1/V1*32.2+V1
+//Density at section 2(in lbm/ft3):
+d2=d1*V1/V2
+//Temperature at section 2(in R):
+T2=p2/d2/R*144
+//Mach number at section 2:
+M2=V2/sqrt(k*R*32.16*T2)
+//Stagnation Temperature at section 2(in R):
+T02=T2*(1+(k-1)/2*M2^2)
+//Stagnation pressure at section 2 (in psia):
+p02=p2*(T02/T2)^(k/(k-1))
+//Mach Number at section 1:
+M1=V1/sqrt(k*R*32.16*T1)
+//Stagnation temperature at section 1(in R):
+T01=T1*(1+(k-1)/2*M1^2)
+//Energy added(in Btu/lbm):
+E=Cp*(T02-T01)
+//Change in entropy(in Btu/(lbm-R)):
+dS=Cp*log(T2/T1)-(Cp-Cv)*log(p2/p1)
+printf("\n\nRESULTS\n\n")
+printf("\n\nVelocity at section 2: %.3f ft/sec\n\n",V2)
+printf("\n\nDensity at section 2: %.3f lbm/ft^3\n\n",d2)
+printf("\n\nTemperature at section 2: %.3f R\n\n",T2)
+printf("\n\nStagnation Temperature at section 2: %.3f R\n\n",T02)
+printf("\n\nStagnation pressure at section 2: %.3f psia\n\n",p02)
+printf("\n\nEnergy added: %.3f Btu/lbm\n\n",E)
+printf("\n\nChange in entropy: %.3f Btu/(lbm-R)\n\n",dS)
diff --git a/746/CH12/EX12.09/12_09.sce b/746/CH12/EX12.09/12_09.sce new file mode 100755 index 000000000..d39b89a78 --- /dev/null +++ b/746/CH12/EX12.09/12_09.sce @@ -0,0 +1,58 @@ +//Temperature and entropy//
+pathname=get_absolute_file_path('12.09.sce')
+filename=pathname+filesep()+'12.09-data.sci'
+exec(filename)
+//Mach nuber at section 1:
+M1=V1/sqrt(k*R*T1)
+//For these value of M1 and M2,the following values are obtained:
+//(To/T0*)1:
+t01=0.7934;
+//(T0/T0*)2:
+t02=0.9787;
+//(p0/p0*)1:
+P01=1.503;
+//(p0/p0*)2:
+P02=1.019;
+//(T/T*)1:
+t1=0.5289;
+//(T/T*)2:
+t2=0.9119;
+//(p/p*)1:
+P1=0.3636;
+//(p/p*)2:
+P2=0.7958;
+//(V/V*)1:
+v1=1.455;
+//(V/V*)2:
+v2=1.146;
+//Value of T2/T1:
+t=t2/t1
+//Temperature at section 2(in K):
+T2=t*T1
+//Value of p2/p1:
+p=P2/P1
+//Pressure at section 2(in kPa):
+p2=p*p1
+//Value of V2/V1:
+v=v2/v1
+//Velocity at section 2(in m/sec):
+V2=v*V1
+//Density at section 2(in kg/m^3):
+d2=p2*1000/R/T2
+//At M1, T/T0=0.5556
+T01=T1/0.5556
+//At M2, T/T0=0.7764
+T02=T2/0.7764
+//Heat added(in kJ/kg):
+E=Cp*(T02-T01)
+//Change in entropy(kJ/kg-K):
+dS=Cp*log(T2/T1)-R*log(p2/p1)/1000
+printf("\n\nRESULTS\n\n")
+printf("\n\nTemperature at section 2: %.3f K\n\n",T2)
+printf("\n\nPressure at section 2: %.3f kPa\n\n",p2)
+printf("\n\nVelocity at section 2: %.3f m/sec\n\n",V2)
+printf("\n\nDensity at section 2: %.3f kg/m^3\n\n",d2)
+printf("\n\nStagnation temperature at section 2: %.3f K\n\n",T02)
+printf("\n\nHeat added: %.3f kJ/kg\n\n",E)
+printf("\n\nChange in entropy: %.3f kJ/kg\n\n",dS)
+
diff --git a/746/CH12/EX12.10/12_10.sce b/746/CH12/EX12.10/12_10.sce new file mode 100755 index 000000000..bff927440 --- /dev/null +++ b/746/CH12/EX12.10/12_10.sce @@ -0,0 +1,43 @@ +//Temperature//
+pathname=get_absolute_file_path('12.10.sce')
+filename=pathname+filesep()+'12.10-data.sci'
+exec(filename)
+//Density at section 1(in kg/m^3):
+d1=p1*1000/R/T1
+//Mach number at section 1:
+M1=V1/sqrt(k*R*T1)
+//Stagnation temperature at section 1(in K):
+T01=T1*(1+(k-1)/2*M1^2)
+//Stagnation pressure at section 1(in kPa):
+p01=p1*(1+(k-1)/2*M1^2)^(k/(k-1))
+//The following values are obtained from the appendix:
+//po2/p01:
+p0=0.7209;
+//T2/T1:
+T=1.687;
+//p2/p1:
+p=4.5;
+//V2/V1:
+V=0.3750;
+//Temperature at section 2 (in K):
+T2=T*T1
+//Pressure at section 2(in kPa):
+p2=p*p1
+//Velocity at section 2(in m/sec):
+V2=V*V1
+//Density at section 2 (in kg/m^3):
+d2=p2*1000/R/T2
+//Stagnation pressure at section 2(in kPa):
+p02=p0*p01
+//Stagnation temperature at section 2(in K):
+T02=T01;
+//Change in entropy(in kJ/(kg-K)):
+dS=-R/1000*log(p0)
+printf("\n\nRESULTS\n\n")
+printf("\n\nTemperature at section 2 : %.3f K\n\n",T2)
+printf("\n\nPressure at section 2: %.3f kPa\n\n",p2)
+printf("\n\nVelocity at section 2: %.3f m/sec\n\n",V2)
+printf("\n\nDensity at section 2 : %.3f kg/m^3\n\n",d2)
+printf("\n\nStagnation pressure at section 2: %.3f kPa\n\n",p02)
+printf("\n\nChange in entropy: %.3f kg-K\n\n",dS)
+printf("\n\nStagnation temperature at section 2: %.3f K\n\n",T02)
diff --git a/746/CH2/EX2.02/2_02.sce b/746/CH2/EX2.02/2_02.sce new file mode 100755 index 000000000..b3d57c048 --- /dev/null +++ b/746/CH2/EX2.02/2_02.sce @@ -0,0 +1,17 @@ +pathname=get_absolute_file_path('2.02.sce')
+filename=pathname+filesep()+'2.02-data.sci'
+exec(filename)
+
+//Viscosity in units of lbf-s/ft^2:
+u1=u/100/454/32.2*30.5
+//Kinematic viscosity (in m/sec^2):
+v=u1/SG/d*(0.305)^2
+//Shear stress on the upper plate(lbf/ft^2):
+tu=u1*U/D*1000
+//Shear stress on the lower plate(in Pa)
+tl=tu*4.45/0.305^2
+printf("\n\nRESULTS\n\n")
+printf("\n\nViscosity in units of lbf-s/ft^2: %1.8f $lbf-s/ft^2\n\n",u1)
+printf("\n\nKinematic viscosity: %1.8f m/sec^2\n\n",v)
+printf("\n\nShear stres on the upeer plate: %f lbf/ft^2\n\n",tu)
+printf("\n\nSear stress on the lower plate: %f Pa\n\n",tl)
diff --git a/746/CH3/EX3.01/3_01.jpg b/746/CH3/EX3.01/3_01.jpg Binary files differnew file mode 100755 index 000000000..74d8becc0 --- /dev/null +++ b/746/CH3/EX3.01/3_01.jpg diff --git a/746/CH3/EX3.01/3_01.sce b/746/CH3/EX3.01/3_01.sce new file mode 100755 index 000000000..089d4a614 --- /dev/null +++ b/746/CH3/EX3.01/3_01.sce @@ -0,0 +1,21 @@ +//liquid level//
+pathname=get_absolute_file_path('3.01.sce')
+filename=pathname+filesep()+'3.01-data.sci'
+exec(filename)
+//Tube diameter(in mm):
+D=1:25;
+D1=D/1000
+[m n]=size(D1)
+for i=1:n
+//Change in liquid level for water(in mm):
+dhw(i)=4*STw*cosd(thetaw)/dw/g/D1(i);
+//Change in liquid level for mercury(in mm):
+dhm(i)=4*STm*cosd(thetam)/dm/g/D1(i);
+end;
+
+//Plotting tube daimeter and water level:
+plot(D1*1000,dhw,'-o')
+//Plotting tube daimeter and mercury level:
+plot(D1*1000,dhm,'-*')
+legend(['Water';'Mercury']);
+xtitle('Liquid level vs Tube diameter','Liquid level(in mm)','Tube diameter(in mm)')
diff --git a/746/CH3/EX3.03/3_03.sce b/746/CH3/EX3.03/3_03.sce new file mode 100755 index 000000000..92dc4e6bd --- /dev/null +++ b/746/CH3/EX3.03/3_03.sce @@ -0,0 +1,8 @@ +//pressure difference//
+pathname=get_absolute_file_path('3.03.sce')
+filename=pathname+filesep()+'3.03-data.sci'
+exec(filename)
+//Pressure difference(in lbf/in^2):
+dp=g*d*(-d1+SGm*d2-SGo*d3+SGm*d4+d5)/12/144
+printf("\n\nRESULTS\n\n")
+printf("\n\nPressure difference between A and B: %f lbf/in^2\n\n",dp)
diff --git a/746/CH3/EX3.04/3_04.sce b/746/CH3/EX3.04/3_04.sce new file mode 100755 index 000000000..4030fcd51 --- /dev/null +++ b/746/CH3/EX3.04/3_04.sce @@ -0,0 +1,40 @@ +//temperature and pressure//
+pathname=get_absolute_file_path('3.04.sce')
+filename=pathname+filesep()+'3.04-data.sci'
+exec(filename)
+//Assuming temperature varies linearly with altitude:
+//Temperature gradient(in F/ft):
+m=(T1-T2)/(z2-z1)
+//Value of g/(m*R):
+v=g/m/R/32.2
+//Pressure at Vail Pass(in inches of Hg):
+p12=p1*((T2+460)/(T1+460))^v
+//Percentage change in density:
+pc1=(p12/p1*(T1+460)/(T2+460)-1)*100
+//Assuming density is constant:
+//Pressure at Vail Pass(in inches of Hg):
+p22=p1*(1-(g*(z2-z1)/(R*32.2)/(T1+460)))
+//Percentage change in density:
+pc2=0;
+//Assuming temperature is constant:
+//Pressure at Vail Pass(in inches of Hg):
+p32=p1*%e^(-g*(z2-z1)/(R*32.2)/(T2+460))
+//Percentage change in density:
+pc3=(p32/p1*(T1+460)/(T1+460)-1)*100
+//For an adiabatic atmosphere:
+p42=p1*((62+460)/(80+460))^(k/(k-1))
+//Percentage change in density:
+pc4=(p42/p1*(T1+460)/(T2+460)-1)*100
+printf("\n\nRESULTS\n\n")
+printf("\n\n1) If temperature varies linearly with altitude\n\n")
+printf("\n\n\tAtmospheric pressure at Vail Pass: %f inches of Hg\n\n",p12)
+printf("\n\n\tPercentage change in density wrt Denver: %f percent\n\n",pc1)
+printf("\n\n2) If density is constant\n\n")
+printf("\n\n\tAtmospheric pressure at Vail Pass: %f inches of Hg\n\n",p22)
+printf("\n\n\tPercentage change in density wrt Denver: %f percent\n\n",pc2)
+printf("\n\n3) If temperature is constant\n\n")
+printf("\n\n\tAtmospheric pressure at Vail Pass: %f inches of Hg\n\n",p32)
+printf("\n\n\tPercentage change in density wrt Denver: %f percent\n\n",pc3)
+printf("\n\n4) For an adiabatic atmosphere\n\n")
+printf("\n\n\tAtmospheric pressure at Vail Pass: %f inches of Hg\n\n",p42)
+printf("\n\n\tPercentage change in density wrt Denver: %f percent\n\n",pc4)
diff --git a/746/CH3/EX3.05/3_05.sce b/746/CH3/EX3.05/3_05.sce new file mode 100755 index 000000000..4cea618b9 --- /dev/null +++ b/746/CH3/EX3.05/3_05.sce @@ -0,0 +1,25 @@ +//force and pressure//
+pathname=get_absolute_file_path('3.05.sce')
+filename=pathname+filesep()+'3.05-data.sci'
+exec(filename)
+//Net force on the gate(in kN):
+Fr=d*g*w*(D*L+L^2/2*sind(theta))
+//Centre of pressure:
+//Calculation for y coordinate:
+ yc=D/sind(theta)+L/2
+ //Area(in m^2):
+ A=L*w
+ //Moment of inertia of rectangular gate(in m^4):
+ Ixx=w*L^3/12
+ //y coordinate(in m):
+ y=yc+Ixx/A/yc
+//Calculation for x coordinate:
+ Ixy=0
+ xc=w/2
+ //x coordinate(in m):
+ x=xc+Ixy/A/xc
+printf("\n\nRESULTS\n\n")
+printf("\n\nNet force on the gate: %f kN\n\n",Fr/1000)
+printf("\n\nCoordinate of centre of pressure:(%0.1f,%0.1f)\n\n",x,y)
+
+
diff --git a/746/CH3/EX3.06/3_06.sce b/746/CH3/EX3.06/3_06.sce new file mode 100755 index 000000000..89f2daf73 --- /dev/null +++ b/746/CH3/EX3.06/3_06.sce @@ -0,0 +1,9 @@ +//force//
+pathname=get_absolute_file_path('3.06.sce')
+filename=pathname+filesep()+'3.06-data.sci'
+exec(filename)
+//Force required to keep the door shut(in lbf):
+function y=f(z),y=b/L*p0*z+d*b/L*(L*z-z^2),endfunction
+Ft=intg(0,L,f)
+printf("\n\nRESULTS\n\n")
+printf("\n\nForce requiredto kep the door shut: %.1f lbf\n\n",Ft)
diff --git a/746/CH3/EX3.07/3_07.sce b/746/CH3/EX3.07/3_07.sce new file mode 100755 index 000000000..8624fca4a --- /dev/null +++ b/746/CH3/EX3.07/3_07.sce @@ -0,0 +1,18 @@ +//force at equilibrium//
+pathname=get_absolute_file_path('3.07.sce')
+filename=pathname+filesep()+'3.07-data.sci'
+exec(filename)
+//Horizontal component of resultant force(in kN):
+Frh=0.5*d*g*w*D^2
+//Line of action of Frh(in m):
+y1=0.5*D+w*D^3/12/(0.5*D)/(w*D)
+//Vertical component of resultant force(in kN):
+function y=q(x), y=d*g*w*(D-sqrt(a*x)),endfunction
+Frv=intg(0,D^2/a,q)
+//Line of acion of Frv(in m):
+function k=f(x), k=d*g*w/Frv*x*(D-sqrt(a*x)),endfunction
+xa=intg(0,D^2/a,f)
+//Force required to keep the gate in equilibrium(in kN):
+Fa=1/l*(xa*Frv+(D-y1)*Frh)
+printf("\n\nRESULTS\n\n")
+printf("\n\nForce required to keep the gate at equilibrium: %f kN\n\n",Fa/1000)
diff --git a/746/CH4/EX4.01/4_01.sce b/746/CH4/EX4.01/4_01.sce new file mode 100755 index 000000000..e50c6df62 --- /dev/null +++ b/746/CH4/EX4.01/4_01.sce @@ -0,0 +1,19 @@ +//Velocity//
+pathname=get_absolute_file_path('4.01.sce')
+filename=pathname+filesep()+'4.01-data.sci'
+exec(filename)
+//If I=integral of(pV.dA):
+//For system: Ics=IA1+IA2+IA3+IA4.
+//For area 1
+IA1=-d*V1*A1
+//For area 3: IA2=d*V3*A3=m3
+IA3=m3
+//For area 4: IA4=-d*V4*A4=-d*Q4
+IA4=-d*Q4
+//For area 2:
+IA2=-IA1-IA3-IA4
+//Velocity at section 2(in ft/sec):
+V2=IA2/d/A2
+//V2 is in the negative y direction
+printf("\n\nRESULTS\n\n")
+printf("\n\nVelocity at section 2: -%.0fj ft/sec\n\n",V2)
diff --git a/746/CH4/EX4.02/4_02.sce b/746/CH4/EX4.02/4_02.sce new file mode 100755 index 000000000..a209b74e0 --- /dev/null +++ b/746/CH4/EX4.02/4_02.sce @@ -0,0 +1,20 @@ +//Mass flow//
+pathname=get_absolute_file_path('4.02.sce')
+filename=pathname+filesep()+'4.02-data.sci'
+exec(filename)
+//If I=integral of(pV.dA):
+//For system: ICS=Iab+Ibc+Icd+Ida
+//But ICS=0
+
+//For Aab:
+function p=f(y),p=-d*U*w*y^0,endfunction
+IAab=intg(0,t,f)
+
+//For Acd:
+function q=g(y),q=d*U*w*(2*y/t-(y/t)^2),endfunction
+IAcd=intg(0,t,g)
+
+//Mass flow rate across surface bc(in kg/sec):
+mbc=(-IAab-IAcd)/1000
+printf("\n\nRESULTS\n\n")
+printf("\n\nMass flow rate across surface bc: %.4f kg/sec\n\n",mbc)
diff --git a/746/CH4/EX4.03/4_03.sce b/746/CH4/EX4.03/4_03.sce new file mode 100755 index 000000000..d9ad2ace3 --- /dev/null +++ b/746/CH4/EX4.03/4_03.sce @@ -0,0 +1,9 @@ +//density//
+pathname=get_absolute_file_path('4.03.sce')
+filename=pathname+filesep()+'4.03-data.sci'
+exec(filename)
+//Rate of change of air density in tank(in (kg/m^3)/s):
+r=-d*v*A/V/10^6
+printf("\n\nRESULTS\n\n")
+printf("\n\nRate of change of air density in tank: %.3f kg/m^3\n\n",r)
+printf("\n\nThe density decreases as is indicated by the negative sign\n\n")
diff --git a/746/CH4/EX4.04/4_04.sce b/746/CH4/EX4.04/4_04.sce new file mode 100755 index 000000000..40e6d7721 --- /dev/null +++ b/746/CH4/EX4.04/4_04.sce @@ -0,0 +1,24 @@ +//Horizontal force//
+pathname=get_absolute_file_path('4.04.sce')
+filename=pathname+filesep()+'4.04-data.sci'
+exec(filename)
+//1) Control Volume selected so that area of left surface is equal to the area of the right surface
+u1=15;
+//Force of support on control volum(in kN):
+function y=f(A),y=-u1*d*V,endfunction
+Rx1=intg(0,0.01,f)
+//Horizontal force on support(in kN):
+Kx=-Rx1
+//2) Control volumes are selected do that the area of the left and right surfaces are equial to the area of the plate
+
+function z=g(A),z=-u1*d*V,endfunction
+Fsx=intg(0,0.01,g)
+//Net force on plate:Fx=0=-Bx-pa*Ap+Rx
+// Rx=pa*Ap+Bx
+//From the above, it is obtained that:
+Rx2=-2.25
+//Horizontal force on support(in kN):
+Kx2=-Rx2
+printf("\n\nRESULTS\n\n")
+printf("\n\nHorizontal force on support: %.3f kN\n\n",Kx/1000)
+
diff --git a/746/CH4/EX4.05/4_05.sce b/746/CH4/EX4.05/4_05.sce new file mode 100755 index 000000000..27c799ef0 --- /dev/null +++ b/746/CH4/EX4.05/4_05.sce @@ -0,0 +1,16 @@ +//Scale//
+pathname=get_absolute_file_path('4.05.sce')
+filename=pathname+filesep()+'4.05-data.sci'
+exec(filename)
+//Weight of water in the tank(in lbf):
+d1=62.4;
+WH2O=d1*A*h
+v=-5;
+//Total body force in negative y direction(lbf):
+function y=f(A),y=-v*d2*V1,endfunction
+F=intg(0,A1,f)
+//Force of scale on control volume(in kN):
+Ry=W+WH2O-F
+printf("\n\nRESULTS\n\n")
+printf("\n\nScale Reading: %.3f lbf\n\n",Ry)
+
diff --git a/746/CH4/EX4.06/4_06.sce b/746/CH4/EX4.06/4_06.sce new file mode 100755 index 000000000..8220dc73e --- /dev/null +++ b/746/CH4/EX4.06/4_06.sce @@ -0,0 +1,10 @@ +//force exerted per unt//
+pathname=get_absolute_file_path('4.06.sce')
+filename=pathname+filesep()+'4.06-data.sci'
+exec(filename)
+//X-component of reaction force per unit width of the gate(in N/m):
+Rxw=(d*(V2^2*D2-V1^2*D1))-(d*g/2*(D1^2-D2^2))
+//Horizontal force exerted per unt width on the gate(in N/m):
+Kxw=-Rxw
+printf("\n\nRESULTS\n\n")
+printf("\n\nHorizontal force exerted per unt width on the gate: %.3f kN/m\n\n",Kxw/1000)
diff --git a/746/CH4/EX4.07/4_07.sce b/746/CH4/EX4.07/4_07.sce new file mode 100755 index 000000000..6f0000427 --- /dev/null +++ b/746/CH4/EX4.07/4_07.sce @@ -0,0 +1,16 @@ +//Force to hold//
+pathname=get_absolute_file_path('4.07.sce')
+filename=pathname+filesep()+'4.07-data.sci'
+exec(filename)
+//Velocity at section 1(in m/sec):
+V1=V2*A2/A1
+//Gauge pressure(in kPa):
+p1g=p1-patm
+u1=V1;u2=-V2;
+//Reaction force component in the x direction(in N):
+Rx=-p1g*A1-u1*d*V1*A1
+//Reaction force component in the y direction(in N):
+Ry=u2*d*V2*A2
+printf("\n\nRESULTS\n\n")
+printf("\n\nForce to hold elbow acting to the left: %.3f kN\n\n",Rx/1000)
+printf("\n\nForce to hold elbow acting downwards: %.3f N\n\n",Ry)
diff --git a/746/CH4/EX4.08/4_08.sce b/746/CH4/EX4.08/4_08.sce new file mode 100755 index 000000000..fd8e5eddf --- /dev/null +++ b/746/CH4/EX4.08/4_08.sce @@ -0,0 +1,7 @@ +//Tension//
+filename=pathname+filesep()+'4.08-data.sci'
+exec(filename)
+//Tension required to pull the belt(in lbf):
+T=Vbelt*m/32.2
+printf("\n\nRESULTS\n\n")
+printf("\n\nTension required to pull the belt: %.3f lbf\n\n",T)
diff --git a/746/CH4/EX4.09/4_09.sce b/746/CH4/EX4.09/4_09.sce new file mode 100755 index 000000000..3415f9e67 --- /dev/null +++ b/746/CH4/EX4.09/4_09.sce @@ -0,0 +1,8 @@ +//pressure required//
+pathname=get_absolute_file_path('4.09.sce')
+filename=pathname+filesep()+'4.09-data.sci'
+exec(filename)
+//Minimum gauge pressure required(in lbf/in^2):
+pg=8/%pi^2*d/D1^4*Q^2*((D1/D2)^4-1)*144
+printf("\n\nRESULTS\n\n")
+printf("Minimum gauge pressure required: %.3f lbf/in^2",pg)
diff --git a/746/CH4/EX4.10/4_10.sce b/746/CH4/EX4.10/4_10.sce new file mode 100755 index 000000000..39073d8e6 --- /dev/null +++ b/746/CH4/EX4.10/4_10.sce @@ -0,0 +1,19 @@ +//Net force//
+pathname=get_absolute_file_path('4.10.sce')
+filename=pathname+filesep()+'4.10-data.sci'
+exec(filename)
+u1=V-U
+u2=(V-U)*cosd(theta)
+v2=(V-U)*sind(theta)
+V1=V-U
+V2=V1
+//X component of moment equation(in N):
+function y=f(A),y=u1*-(d*V1),endfunction
+function z=g(A),z=u2*d*V2,endfunction
+Rx=intg(0,A,f)+intg(0,A,g)
+
+//Y component of moment equation(in N):
+function a=h(A),a=v2*d*V1,endfunction
+Ry=intg(0,A,h) //This is after neglecting weight of vane and the water.
+printf("\n\nRESULTS\n\n")
+printf("\n\nNet force on the vane: %.3f i+%.2f j kN\n\n",Rx/1000,Ry/1000)
diff --git a/746/CH4/EX4.11/4_11.sce b/746/CH4/EX4.11/4_11.sce new file mode 100755 index 000000000..718c78372 --- /dev/null +++ b/746/CH4/EX4.11/4_11.sce @@ -0,0 +1,16 @@ +//PLOTTING//
+pathname=get_absolute_file_path('4.11.sce')
+filename=pathname+filesep()+'4.11-data.sci'
+exec(filename)
+//Evaluating the value of Vb:
+Vb=V*(1-cosd(theta))*d*A/M
+//Value of U/V for various values of t
+t=0:20;
+[m n]=size(t)
+for i=1:n
+ U_V(i)=Vb*t(i)/(1+Vb*t(i));
+end
+
+//Plotting U/V vs t:
+plot(t,U_V)
+xtitle('U/V vs t','t (in sec)','U/V')
diff --git a/746/CH4/EX4.12/4_12.sce b/746/CH4/EX4.12/4_12.sce new file mode 100755 index 000000000..86464c549 --- /dev/null +++ b/746/CH4/EX4.12/4_12.sce @@ -0,0 +1,11 @@ +//Velocity of rocket//
+pathname=get_absolute_file_path('4.12.sce')
+filename=pathname+filesep()+'4.12-data.sci'
+exec(filename)
+//Acceleration of rocket at t=0(in m/sec^2):
+Ve*me/M0-g
+//Velocity of rocket at t=10 (in m/sec):
+function y=f(t),y=Ve*me/(M0-me*t)-g,endfunction
+Vcv=intg(0,t,f)
+printf("\n\nRESULTS\n\n")
+printf("\n\nVelocity of rocket at t=10: %.1f m/sec\n\n",Vcv)
diff --git a/746/CH4/EX4.14/4_14.sce b/746/CH4/EX4.14/4_14.sce new file mode 100755 index 000000000..755aa6047 --- /dev/null +++ b/746/CH4/EX4.14/4_14.sce @@ -0,0 +1,15 @@ +//Relative speed and friction//
+pathname=get_absolute_file_path('4.14.sce')
+filename=pathname+filesep()+'4.14-data.sci'
+exec(filename)
+//Area of jet(in mm^2):
+Ajet=%pi/4*D^2
+//Jet speed relative to the nozzle(in m/sec):
+Vrel=Q/2/Ajet*10^6/60/1000
+//Value of w*R in m/sec:
+wR=w*R*2*%pi/60/1000
+//Friction torque at pivot(in N-m):
+Tf=R*(Vrel*cosd(alpha)-wR)*d*Q/1000/60/1000
+printf("\n\nRESULTS\n\n")
+printf("\n\nJet speed relative to each nozzle: %.2f m/sec\n\n",Vrel)
+printf("\n\nFriction torque at pivot: %.5f N-m\n\n",Tf)
diff --git a/746/CH4/EX4.16/4_16.sce b/746/CH4/EX4.16/4_16.sce new file mode 100755 index 000000000..4feb7092b --- /dev/null +++ b/746/CH4/EX4.16/4_16.sce @@ -0,0 +1,11 @@ +//Rate of heat//
+pathname=get_absolute_file_path('4.16.sce')
+filename=pathname+filesep()+'4.16-data.sci'
+exec(filename)
+//Velocity at exit(in ft/sec):
+V2=m*R*(T2+460)/A2/p2/144
+//As power input is to CV, Ws=-600
+//Rate of heat transfer(in Btu/sec):
+Q=Ws*550/778+m*cp*(T2-T1)+m*V2^2/2/32.2/778
+printf("\n\nRESULTS\n\n")
+printf("\n\nRate of heat transfer: %.3f Btu/sec\n\n",Q)
diff --git a/746/CH4/EX4.17/4_17.sce b/746/CH4/EX4.17/4_17.sce new file mode 100755 index 000000000..38b8f55c8 --- /dev/null +++ b/746/CH4/EX4.17/4_17.sce @@ -0,0 +1,10 @@ +//Mass flow rate//
+pathname=get_absolute_file_path('4.17.sce')
+filename=pathname+filesep()+'4.17-data.sci'
+exec(filename)
+//Density of tank(in kg/m^3):
+d=(p1+patm)/R/T
+//Mass flow rate of air into the tank(in kg/sec):
+m=d*V*cv*r/R/T*1000
+printf("\n\nRESULTS\n\n")
+printf("\n\nMass flow rate of air into the tank: %.3f g/sec\n\n",m)
diff --git a/746/CH5/EX5.02/5_02.sce b/746/CH5/EX5.02/5_02.sce new file mode 100755 index 000000000..bbf900b84 --- /dev/null +++ b/746/CH5/EX5.02/5_02.sce @@ -0,0 +1,8 @@ +//Rate of change//
+pathname=get_absolute_file_path('5.02.sce')
+filename=pathname+filesep()+'5.02-data.sci'
+exec(filename)
+//Rate of change of density with time(in kg/m^3-s):
+r=-d*V/L
+printf("\n\nRESULTS\n\n")
+printf("\n\nRate of change of density with time: %.1f kg/m^3-s\n\n",r)
diff --git a/746/CH5/EX5.07/5_07.sce b/746/CH5/EX5.07/5_07.sce new file mode 100755 index 000000000..186042f6f --- /dev/null +++ b/746/CH5/EX5.07/5_07.sce @@ -0,0 +1,15 @@ +//angular and rotation//
+pathname=get_absolute_file_path('5.07.sce')
+filename=pathname+filesep()+'5.07-data.sci'
+exec(filename)
+//At point b, u=3 mm/sec
+u=3;
+//Displacemet of b(in mm):
+xb=u*t
+//Rate of angular deformation(in s^-1):
+def=U/h
+//Rate of rotation(in s^-1):
+rot=-0.5*U/h
+printf("\n\nRSULTS\n\n")
+printf("\n\nRate of angular deformation: %.1f /sec\n\n",def)
+printf("\n\nRate of rotation: %.1f /sec\n\n",rot)
diff --git a/746/CH5/EX5.08/5_08.sce b/746/CH5/EX5.08/5_08.sce new file mode 100755 index 000000000..f01cb0ced --- /dev/null +++ b/746/CH5/EX5.08/5_08.sce @@ -0,0 +1,31 @@ +//Rates and area//
+pathname=get_absolute_file_path('5.08.sce')
+filename=pathname+filesep()+'5.08-data.sci'
+exec(filename)
+//Value of T:
+T=log(3/2)/A
+x0=1:2;
+y0=1:2;
+for i=1:2
+ for j=1:2
+ //For X coordinate:
+ X(i)(j)=x0(i)*%e^(A*T)
+ //For Y coordinate:
+ Y(i)(j)=y0(j)*%e^(-A*T)
+ end
+end
+plot(X,Y)
+//Rates of linear deformation in X direction:
+Ax=0.3;
+//Rate of linear deformation in the y direction:
+Ay=-0.3;
+//Rate of volume dilation(s^-1):
+v=A-A
+//Area of abcd:
+A1=1;
+//Area of a'b'c'd':
+A2=(3-3/2)*(4/3-2/3)
+printf("\n\nRESULTS\n\n")
+printf("\n\nRates of linear deformation in X and Y direction: %.1f /s , %.1f /s\n\n",Ax,Ay)
+printf("\n\nRate of volume dilation: %.0f /sec\n\n",v)
+printf("\n\nArea of abcd and a,b,c,d:%.1f m^2, %.1f m^\n\n",A1,A2)
diff --git a/746/CH5/EX5.09/5_09.sce b/746/CH5/EX5.09/5_09.sce new file mode 100755 index 000000000..f13da2f7a --- /dev/null +++ b/746/CH5/EX5.09/5_09.sce @@ -0,0 +1,8 @@ +//Volume flow rate//
+pathname=get_absolute_file_path('5.09.sce')
+filename=pathname+filesep()+'5.09-data.sci'
+exec(filename)
+//Volume flow rate(in m^3/sec):
+Q=d*g*sind(theta)*b*(h/1000)^3*1000/u/3
+printf("RESULTS")
+printf("\n\nVolume flow rate: %.4f m^3/sec\n\n",Q)
diff --git a/746/CH6/EX6.01/6_01.sce b/746/CH6/EX6.01/6_01.sce new file mode 100755 index 000000000..d0d412fd5 --- /dev/null +++ b/746/CH6/EX6.01/6_01.sce @@ -0,0 +1,10 @@ +//Volume flow rate//
+pathname=get_absolute_file_path('06.01.sce')
+filename=pathname+filesep()+'06.01-data.sci'
+exec(filename)
+//Velocity of flow(in m/sec):
+V=sqrt(dw/log((r+w)/r)*g/da*p/1000)
+//Volume flow rate(in m^3/sec):
+Q=V*(d*w)
+printf("\n\nRESULTS\n\n")
+printf("\n\nVolume flow rate: %.3f m^3/sec\n\n",Q)
diff --git a/746/CH6/EX6.02/6_02.sce b/746/CH6/EX6.02/6_02.sce new file mode 100755 index 000000000..bf235e779 --- /dev/null +++ b/746/CH6/EX6.02/6_02.sce @@ -0,0 +1,8 @@ +//Velocity of flow//
+pathname=get_absolute_file_path('06.02.sce')
+filename=pathname+filesep()+'06.02-data.sci'
+exec(filename)
+//Velocity of flow(in m/sec):
+V=sqrt(2*dw*g*p/1000*SG/da)
+printf("\n\nRESULTS\n\n")
+printf("\n\nVelocity of flow: %.3f m/sec\n\n",V)
diff --git a/746/CH6/EX6.03/6_03.sce b/746/CH6/EX6.03/6_03.sce new file mode 100755 index 000000000..5f4c26300 --- /dev/null +++ b/746/CH6/EX6.03/6_03.sce @@ -0,0 +1,9 @@ +//prssure required//
+filename=pathname+filesep()+'06.03-data.sci'
+exec(filename)
+//Velocity of flwat the inlet(in m/sec):
+V1=Ae/Ai*V2
+//Gauge pressure required at the inlet(in kPa):
+p=0.5*da*(V2^2-V1^2)
+printf("\n\nRESULTS\n\n")
+printf("\n\nGauge prssure required at the inlet: %.3f kPa\n\n",p/1000)
diff --git a/746/CH6/EX6.04/6_04.sce b/746/CH6/EX6.04/6_04.sce new file mode 100755 index 000000000..81e13bb2d --- /dev/null +++ b/746/CH6/EX6.04/6_04.sce @@ -0,0 +1,11 @@ +//Speed and pressure//
+pathname=get_absolute_file_path('06.04.sce')
+filename=pathname+filesep()+'06.04-data.sci'
+exec(filename)
+//Speed of water at exit(in m/sec):
+V2=sqrt(2*g*z)
+//Pressure at point A in the flow(kPa):
+pA=p1+d*g*(0-l)-0.5*d*V2^2
+printf("\n\nRESULTS\n\n")
+printf("\n\nSpeed of water at exit: %.3f m/sec\n\n",V2)
+printf("\n\nPressure at point A in the flow: %3f kPa\n\n",pA/1000)
diff --git a/746/CH6/EX6.05/6_05.sce b/746/CH6/EX6.05/6_05.sce new file mode 100755 index 000000000..c72c372cf --- /dev/null +++ b/746/CH6/EX6.05/6_05.sce @@ -0,0 +1,11 @@ +//flow//
+pathname=get_absolute_file_path('06.05.sce')
+filename=pathname+filesep()+'06.05-data.sci'
+exec(filename)
+//Velocity of flow at the exit(in ft/sec):
+V2=sqrt(2*g*(Du-Dd/12))
+//Volume flow rate/width(ft^2/sec):
+Q=V2*Dd/12
+printf("\n\nRESULTS\n\n")
+printf("\n\nVelocity of flow at the exit: %.3f ft/sec\n\n",V2)
+printf("\n\nVolume flow rate/width: %.3f ft^2/sec\n\n",Q)
diff --git a/746/CH6/EX6.06/6_06.sce b/746/CH6/EX6.06/6_06.sce new file mode 100755 index 000000000..3b954de53 --- /dev/null +++ b/746/CH6/EX6.06/6_06.sce @@ -0,0 +1,16 @@ +//pressure//
+pathname=get_absolute_file_path('06.06.sce')
+filename=pathname+filesep()+'06.06-data.sci'
+exec(filename)
+//Pressure of air at 1000 m(in N/m^2):
+p=P1*pa
+//Density of air at 1000m(in kg/m^3):
+d=D1*da
+//Stagnation pressure at A(in kPa):
+p0A=p+0.5*d*(V*1000/3600)^2
+//Static pressure at B(in kPa):
+pB=p+d/2*((V*1000/3600)^2-Vb^2)
+printf("\n\nRESULTS\n\n")
+printf("\n\nStagnation pressure at A: %.3f kPa\n\n",p0A/1000)
+printf("\n\nStatic pressure at B: %.3f kPa\n\n",pB/1000)
+
diff --git a/746/CH6/EX6.08/6_08.sce b/746/CH6/EX6.08/6_08.sce new file mode 100755 index 000000000..1b322d0bd --- /dev/null +++ b/746/CH6/EX6.08/6_08.sce @@ -0,0 +1,12 @@ +//temperature//
+pathname=get_absolute_file_path('06.08.sce')
+filename=pathname+filesep()+'06.08-data.sci'
+exec(filename)
+//Velocity of flow at exit(in ft/sec):
+V4=sqrt(2*g*(z3-0))
+//Mass flow rate of water(in slug/sec):
+m=d*V4*A4/144
+//Rise in temperature between points 1 and 2(in R):
+T=Q*3413/3600/m/32.2
+printf("\n\nRESULTS\n\n")
+printf("\n\nRise in temperature between points 1 and 2: %.3f R\n\n",T)
diff --git a/746/CH6/EX6.09/6_09.sce b/746/CH6/EX6.09/6_09.sce new file mode 100755 index 000000000..ba6a27459 --- /dev/null +++ b/746/CH6/EX6.09/6_09.sce @@ -0,0 +1,15 @@ +//Streamline flow//
+pathname=get_absolute_file_path('06.09.sce')
+filename=pathname+filesep()+'06.09-data.sci'
+exec(filename)
+t=0:5
+//Value of sqrt(2gh):
+x=sqrt(2*g*h)
+//Value of 1/2L*sqrt(2gh):
+y=1/2/L*x
+[m n]=size(t)
+i=1:n;
+//Velocity(in m/sec):
+V2=x*tanh(y*t(i))
+plot(t,V2);
+xtitle('Streamline flow from 1 to 2','Time(in s)','V2(in m/sec)')
diff --git a/746/CH7/EX7.04/7_04.sce b/746/CH7/EX7.04/7_04.sce new file mode 100755 index 000000000..b6e3ab182 --- /dev/null +++ b/746/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/746/CH7/EX7.05/7_05.sce b/746/CH7/EX7.05/7_05.sce new file mode 100755 index 000000000..e0e2e8ef6 --- /dev/null +++ b/746/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/746/CH7/EX7.06/7_06.sce b/746/CH7/EX7.06/7_06.sce new file mode 100755 index 000000000..ea7c51970 --- /dev/null +++ b/746/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)
diff --git a/746/CH8/EX8.01/8_01.sce b/746/CH8/EX8.01/8_01.sce new file mode 100755 index 000000000..cddb2de1f --- /dev/null +++ b/746/CH8/EX8.01/8_01.sce @@ -0,0 +1,15 @@ +//Leakage flow rate//
+pathname=get_absolute_file_path('8.01.sce')
+filename=pathname+filesep()+'8.01-data.sci'
+exec(filename)
+//Leakage flow rate (in mm^3/sec):
+Q=%pi/12*D*a^3*(p1-p2)*10^3/u/L
+//Velocity of flow(in m/sec):
+V=Q/%pi/D/a/1000
+//Specific gravity of SAE 10W oil:
+SG=0.92;
+//Reynolds Number:
+Re=SG*dw*V*a/u/1000
+//As Re<1400, flow is laminar.
+printf("\n\nRESULTS\n\n")
+printf("\n\nLeakage flow rate: %.3f mm^3/sec\n\n",Q)
diff --git a/746/CH8/EX8.02/8_02.sce b/746/CH8/EX8.02/8_02.sce new file mode 100755 index 000000000..904e37f22 --- /dev/null +++ b/746/CH8/EX8.02/8_02.sce @@ -0,0 +1,15 @@ +//Torque and power//
+pathname=get_absolute_file_path('8.02.sce')
+filename=pathname+filesep()+'8.02-data.sci'
+exec(filename)
+//Shear stres (in lbf/ft^2):
+Tyx=u*N*2*%pi/60*D/2/(a/2)
+//Torqe(in inches-lbf):
+T=%pi/2*Tyx*D^2*L/144
+//Power dissipated in the bearing(in hp):
+P=T*N/60*2*%pi/12/550
+//Reynolds number:
+Re=SG*p*N*2*%pi/60*1.5*a/2/u/144
+printf("\n\nRESULTS\n\n")
+printf("\n\nTorque: %.3f inches-lbf\n\n",T)
+printf("\n\nPower dissipated in the bearing: %.3f hp\n\n",P)
diff --git a/746/CH8/EX8.04/8_04.sce b/746/CH8/EX8.04/8_04.sce new file mode 100755 index 000000000..e65781809 --- /dev/null +++ b/746/CH8/EX8.04/8_04.sce @@ -0,0 +1,12 @@ +//Viscosity of fluid//
+pathname=get_absolute_file_path('8.04.sce')
+filename=pathname+filesep()+'8.04-data.sci'
+exec(filename)
+//Viscosity of the liquid(in N-s/m^2):
+u=%pi/128*p*1000*D^4/Q/L/1000
+//Velocity(in m/sec)
+V=Q/(%pi/4*D^2)/1000
+//Reynolds number:
+Re=d*V*D/u/1000
+printf("\n\nRESULTS\n\n")
+printf("\n\nViscosity of fluid %.3f N-s/m^2\n\n",u)
diff --git a/746/CH8/EX8.05/8_05.sce b/746/CH8/EX8.05/8_05.sce new file mode 100755 index 000000000..7d3d6eac6 --- /dev/null +++ b/746/CH8/EX8.05/8_05.sce @@ -0,0 +1,10 @@ +//required//
+pathname=get_absolute_file_path('8.05.sce')
+filename=pathname+filesep()+'8.05-data.sci'
+exec(filename)
+//Reservoir depth required to maintain flow(in m):
+D1=8*Q^2/(%pi)^2/D^4/g*(f*L/D+K+1)
+//Reynolds number:
+Re=4*d*Q/((%pi)*u*D)
+printf("\n\nRESULTS\n\n")
+printf("\n\nReservoir depth required to maintain flow: %.3f m\n\n",D1)
diff --git a/746/CH8/EX8.06/8_06.sce b/746/CH8/EX8.06/8_06.sce new file mode 100755 index 000000000..9a2b9934a --- /dev/null +++ b/746/CH8/EX8.06/8_06.sce @@ -0,0 +1,13 @@ +//Maximum and power//
+pathname=get_absolute_file_path('8.06.sce')
+filename=pathname+filesep()+'8.06-data.sci'
+exec(filename)
+//Velocity of flow(in ft/sec):
+V=Q/24/3600/(%pi/4*(D/12)^2)*42/7.48
+//Maximum spacing(in ft):
+L=2/f*D/12*(p2-p1)/(SG*d)/V^2*144
+//Power needed at each pump(in hp):
+Win=1/Effp*V*%pi/4*(D/12)^2*(p2-p1)/550*144
+printf("\n\nRESULTS\n\n")
+printf("\n\nMaximum spacing: %.3f feet\n\n",L)
+printf("\n\nPower needed at each pump: %.3f hp\n\n",Win)
diff --git a/746/CH8/EX8.07/8_07.sce b/746/CH8/EX8.07/8_07.sce new file mode 100755 index 000000000..cc05e2cd7 --- /dev/null +++ b/746/CH8/EX8.07/8_07.sce @@ -0,0 +1,10 @@ +//Volume low//
+pathname=get_absolute_file_path('8.07.sce')
+filename=pathname+filesep()+'8.07-data.sci'
+exec(filename)
+//Velocity(in ft/sec):
+V2=sqrt(2*g*l/(f*((L+l)/D*12+8)+1))
+//Volume flow rate(in gpm):
+Q=V2*%pi*(D/12)^2/4*7.48*60
+printf("\n\nRESULTS\n\n")
+printf("\n\nVolume low rate: %.3f\n\n",Q)
diff --git a/746/CH8/EX8.08/8_08.sce b/746/CH8/EX8.08/8_08.sce new file mode 100755 index 000000000..7d332823e --- /dev/null +++ b/746/CH8/EX8.08/8_08.sce @@ -0,0 +1,25 @@ +//Minimum diameter//
+pathname=get_absolute_file_path('8.08.sce')
+filename=pathname+filesep()+'8.08-data.sci'
+exec(filename)
+//Value of dPmax(in psi):
+dPmax=p1-p2
+//Q in cubic feet/sec:
+Q1=1500/60/7.48;
+//Initially assume diameter to be 4inches:
+D=4;
+//Reynolds number:
+Re=4*Q1/%pi/v/D*12
+//For this value,
+f=0.012;
+ dP=8*f*L*p*Q1^2/(%pi)^2/D^5*1728;
+while(dP>dPmax)
+ dP=8*f*L*p*Q1^2/(%pi)^2/D^5*1728;
+ if(dP<dPmax)
+ break
+ else
+ D=D+1;
+ end
+end
+printf("\n\nRESULTS\n\n")
+printf("Minimum diameter that can be used:%.1f inches\n\n",D)
diff --git a/746/CH8/EX8.09/8_09.sce b/746/CH8/EX8.09/8_09.sce new file mode 100755 index 000000000..f2de249c9 --- /dev/null +++ b/746/CH8/EX8.09/8_09.sce @@ -0,0 +1,20 @@ +//Loss Coefficient//
+pathname=get_absolute_file_path('8.09.sce')
+filename=pathname+filesep()+'8.09-data.sci'
+exec(filename)
+//Average velocity (in ft/s):
+V2=4/%pi*Q/D^2*144
+//Reynolds number:
+Re=V2*D/v/12
+//For this value,
+f=0.013;
+//Power law exponent:
+n=-1.7+1.8*log10(Re)
+//Value of V/U:
+v_u=2*n^2/(n+1)/(2*n+1)
+//Value of alpha:
+alpha=(1/v_u)^3*2*n^2/(3+n)/(3+2*n)
+//Loss Coefficient for a square edged entrance:
+K=2*g*h/V2^2-f*L/D*12-alpha;
+printf("\n\nRESULTS\n\n")
+printf("\n\nLoss Coefficient for a square edged entrance: %.3f \n\n",K)
diff --git a/746/CH8/EX8.10/8_10.sce b/746/CH8/EX8.10/8_10.sce new file mode 100755 index 000000000..29c0d54dd --- /dev/null +++ b/746/CH8/EX8.10/8_10.sce @@ -0,0 +1,20 @@ +//Volume and increase//
+pathname=get_absolute_file_path('8.10.sce')
+filename=pathname+filesep()+'8.10-data.sci'
+exec(filename)
+//Velocity V1(in m/s):
+V1=sqrt(2*g*z0/1.04)
+//Volume flow rate(in m^3/sec):
+Q=V1*%pi*D^2/4
+Kdiff=1-1/A_R^2-Cp
+//For 2nd case:
+//Velocity(in m/s):
+V1=sqrt(2*g*z0/0.59)
+//Volume flow rate(in m^3/s):
+Qd=V1*%pi*D^2/4
+//Increase in discharge after addition of diffuser is:
+dQ=(Qd-Q)/Q*100
+printf("\n\nRESULTS\n\n")
+printf("\n\nVolume flow rate in case1: %.3f m^3/sec\n\n",Q)
+printf("\n\nVolume flow rate in case 2: %.3f m^3/sec\n\n",Qd)
+printf("\n\nIncrease in discharge after addition of diffuser is: %.3f percent\n\n",dQ)
diff --git a/746/CH8/EX8.11/8_11.sce b/746/CH8/EX8.11/8_11.sce new file mode 100755 index 000000000..1ac0346b6 --- /dev/null +++ b/746/CH8/EX8.11/8_11.sce @@ -0,0 +1,26 @@ +//Diameter and head//
+pathname=get_absolute_file_path('8.11.sce')
+filename=pathname+filesep()+'8.11-data.sci'
+exec(filename)
+
+//Value of K*B^2:
+K_B=Q/(%pi/4*D^2)*sqrt(0.5*d1/g/d2/h)
+//Reynods number:
+ReD1=4/%pi*Q/D/v
+//By trial and error method, the value of beta is fixed at:
+betta=0.66;
+//K is then:
+K=K_B/betta^2
+//Diameter of orifice plate(in m):
+Dt=betta*D
+//Value of p3-p2(in N/m^2):
+P1=d1*Q^2/(%pi/4*D^2)^2*(1/0.65/betta^2-1)
+//Value of p1-p2(in N/m^2):
+P2=d2*g*h
+//Head loss between sections 1 and 3(in N-m/kg):
+hLT=(P2-P1)/d1
+//Expressing the permanent pressure as a fractio of the meter differential:
+C=(P2-P1)/P2
+printf("\n\n\nRESULTS\n\n")
+printf("\n\nDiameter of the orifice: %.3f m\n\n",Dt)
+printf("\n\nHead loss between secions 1 and 3: %.3f N-m/kg\n\n",hLT)
diff --git a/746/CH9/EX9.01/9_01.sce b/746/CH9/EX9.01/9_01.sce new file mode 100755 index 000000000..755aef7b6 --- /dev/null +++ b/746/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/746/CH9/EX9.04/9_04.sce b/746/CH9/EX9.04/9_04.sce new file mode 100755 index 000000000..db1fd7c5a --- /dev/null +++ b/746/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/746/CH9/EX9.05/9_05.sce b/746/CH9/EX9.05/9_05.sce new file mode 100755 index 000000000..8e146f23f --- /dev/null +++ b/746/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/746/CH9/EX9.06/9_06.sce b/746/CH9/EX9.06/9_06.sce new file mode 100755 index 000000000..3786c2baa --- /dev/null +++ b/746/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/746/CH9/EX9.07/9_07.sce b/746/CH9/EX9.07/9_07.sce new file mode 100755 index 000000000..3f1d56c77 --- /dev/null +++ b/746/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/746/CH9/EX9.08/9_08.sce b/746/CH9/EX9.08/9_08.sce new file mode 100755 index 000000000..54e4c9eca --- /dev/null +++ b/746/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/746/CH9/EX9.09/9_09.sce b/746/CH9/EX9.09/9_09.sce new file mode 100755 index 000000000..cef9919e3 --- /dev/null +++ b/746/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)
diff --git a/746/DEPENDENCIES/10_01.sci b/746/DEPENDENCIES/10_01.sci new file mode 100755 index 000000000..fd0ee70a5 --- /dev/null +++ b/746/DEPENDENCIES/10_01.sci @@ -0,0 +1,10 @@ +//Volume flow rate in gpm:
+Q= 150;
+//Value of Vrb2 in ft/sec:
+Vrb2=10;
+//Radius of outter impeller in inches:
+R2=2;
+//Impeller Speed in rpm:
+w=3450;
+//Density of air in slug/ft^3
+p=1.94;
diff --git a/746/DEPENDENCIES/10_02.sci b/746/DEPENDENCIES/10_02.sci new file mode 100755 index 000000000..e612973fa --- /dev/null +++ b/746/DEPENDENCIES/10_02.sci @@ -0,0 +1,14 @@ +//Tip Diameter in metres: +Dt=1.1; +//Hub Diameter in metres: +Dh=0.8; +//Operating Speed in rpm: +w=1200; +//Absolute inlet angle in degrees: +alpha1=30; +//Blade inlet angle in degrees: +betta1=30; +//Blade outlet angle in degrees: +betta2=60; +//Density of air in kg/m^3 +p=1.23;
\ No newline at end of file diff --git a/746/DEPENDENCIES/10_03.sci b/746/DEPENDENCIES/10_03.sci new file mode 100755 index 000000000..4418e53e3 --- /dev/null +++ b/746/DEPENDENCIES/10_03.sci @@ -0,0 +1,24 @@ +//Rate of flow in gm:
+Q=[0 500 800 1000 1100 1200 1400 1500];
+//Suction pressure in psig:
+ps=[ 0.65 0.25 -0.35 -0.92 -1.24 -1.62 -2.42 -2.89];
+//Discharge pressure in psig:
+pd=[53.3 48.3 42.3 36.9 33 27.8 15.3 7.3];
+//Motor Current in amps:
+I=[18 26.2 31 33.9 35.2 36.3 38 39];
+//Acceleration due to gravity in ft/s^2:
+g=32.2;
+//Value of Zs in feet
+zs=1;
+//Density of air in slug/ft^3:
+px=1.94;
+//Value of ZD in feet:
+zd=3;
+//Density of fluid in slug/ft^3:
+py=1000;
+//Motor Efficiency:
+Effm=0.9;
+//Motor Supply in volts:
+E=460;
+//Power Factor:
+PF=0.875;
diff --git a/746/DEPENDENCIES/10_06.sci b/746/DEPENDENCIES/10_06.sci new file mode 100755 index 000000000..9d36e828b --- /dev/null +++ b/746/DEPENDENCIES/10_06.sci @@ -0,0 +1,8 @@ +//Head in Us customary units:
+Hus=21.9;
+//Volume flow rate in US customary units:
+Qus=300;
+//Working seed in rpm:
+N=1170;
+//Aceleration due to graviy in m/s^2
+g=9.81;
\ No newline at end of file diff --git a/746/DEPENDENCIES/10_07.sci b/746/DEPENDENCIES/10_07.sci new file mode 100755 index 000000000..d1e66c4a9 --- /dev/null +++ b/746/DEPENDENCIES/10_07.sci @@ -0,0 +1,12 @@ +//Volume flow rate(in gpm) at shut off condition for N1:
+Q1so=0;
+//Volume flow(in gpm) rate at best efficiency for N1:
+Q1be=300;
+//Head(in feet) at shut off condition for N1:
+H1so=25;
+//Head(in feet) at best efficiency condition for N1:
+H1be=21.9
+//Operation Speed 1:
+N1=1170;
+//Operation speed 2:
+N2=1750;
\ No newline at end of file diff --git a/746/DEPENDENCIES/10_08.sci b/746/DEPENDENCIES/10_08.sci new file mode 100755 index 000000000..6b0113905 --- /dev/null +++ b/746/DEPENDENCIES/10_08.sci @@ -0,0 +1,32 @@ +//For 5 inch nominal pipe line, diameter D:
+Di=5.047;
+//Length of pipeline(in feet):
+L=6;
+//Operatng spped (in rpm):
+N=1750;
+//Water level abovepump centreline(in feet):
+h=3.5;
+//Temperature 1 of water(in Farenheit):
+T1=80;
+//Temperature 2 of water (in Farenheit):
+T2=180;
+//Volume flow rate of water(in gpm):
+Q=1000;
+//Minor loss Coefficients:
+K=0.5; SE=30; OGV=8;
+//Atmospheric pressure(in lbf/in^2):
+patm=14.7;
+//Density of air(slug/ft^3):
+p=1.93;
+//Acceleration due to gravity(in ft/sec^2):
+g=32.2;
+//Head(in feet) due to vapor pressure of water for T =80F:
+Hv1=1.17;
+//Head(in feet) due to vapor pressure of water for T =180F:
+Hv2=17.3;
+//Kinematic viscosity of water at 80F:
+v=0.927e-5;
+//Value of discharges for plotting NPSHR(in gpm):
+Qh=[500 700 900 1100 1300]
+//Values of NPSHR obtained from Fig. D3 of appendix D:
+NPSHRp=[7 8 9.5 12 16]
\ No newline at end of file diff --git a/746/DEPENDENCIES/10_11.sci b/746/DEPENDENCIES/10_11.sci new file mode 100755 index 000000000..562d62a80 --- /dev/null +++ b/746/DEPENDENCIES/10_11.sci @@ -0,0 +1,21 @@ +//Diameter of fan 1 (in inches):
+D1=36;
+//Operating speed of fan 1(in rpm):
+N1=600
+//Density of air used in fan 1(in lbm/ft^3):
+d1=0.075;
+//Diameter of fan 2(in inches):
+D2=42;
+//Operating speed of fan 2(in rpm):
+N2=1150;
+//Density of aif usd in fan 2(in lbm/ft^3):
+d2= 0.045;
+//The following values are obtained from the given graph
+//Values of volume flow rate(in cfm) through fan 1:
+Q1= [0 10000 20000 30000 40000 50000 60000];
+//Values of pressure( in inches of H2O):
+p1=[ 3.68 3.75 3.50 2.96 2.12 1.02 0];
+//Values of power(in hp):
+P1=[ 11.1 15.1 18.6 21.4 23.1 23.1 21.0];
+//Efficiency (in %):
+Eff=[0 37 59 65 57 34 0];
diff --git a/746/DEPENDENCIES/10_12.sce b/746/DEPENDENCIES/10_12.sce new file mode 100755 index 000000000..6dfe1e933 --- /dev/null +++ b/746/DEPENDENCIES/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/DEPENDENCIES/10_12.sci b/746/DEPENDENCIES/10_12.sci new file mode 100755 index 000000000..5310deba1 --- /dev/null +++ b/746/DEPENDENCIES/10_12.sci @@ -0,0 +1,18 @@ +//Operation speed(in rpm):
+N=2000;
+//Volume flow rate(in gpm):
+Q=20;
+//Pressure(in psig):
+p=1500;
+//Actual Pump Displacement(in^3/rev):
+va=5.9;
+//Volume flow rate at operating condition(in gpm):
+Qo=46.5;
+//Volume flow rate at maximum delivery(in gpm):
+Qe=48.5;
+//Pressure at operation condition(in psi):
+po1=1500;
+//Efficiency of pump at operating condition:
+Effp=0.84;
+//Pressure at operating condition case 2(in psig):
+po2=3000;
\ No newline at end of file diff --git a/746/DEPENDENCIES/10_14.sci b/746/DEPENDENCIES/10_14.sci new file mode 100755 index 000000000..e6a9085d1 --- /dev/null +++ b/746/DEPENDENCIES/10_14.sci @@ -0,0 +1,15 @@ +//Total propulsion power requirement(in MW):
+P=11.4;
+//From the given curves,Value of coefficients atoptimum efficiency are as follows:
+//Speed of advance coefficient:
+J=0.85;
+//Thrust Coefficient:
+Cf=0.1;
+//Torque Coefficint:
+Ct=0.02;
+//Efficiency:
+Eff=0.66;
+//Velocity of ship(in m/sec):
+V=6.69;
+//Density of water(in kg/m^3):
+p=1025;
\ No newline at end of file diff --git a/746/DEPENDENCIES/10_16.sci b/746/DEPENDENCIES/10_16.sci new file mode 100755 index 000000000..11f28ee33 --- /dev/null +++ b/746/DEPENDENCIES/10_16.sci @@ -0,0 +1,15 @@ +//Diameter of windmill(in m):
+D=26;
+//Operating speed(in rpm):
+N=20;
+//Wind speed(in km/hr):
+V=36;
+//Power Output(in W):
+Po=41000;
+//Maximum efficiency occurs in following conditions:
+//Efficiency:
+Eff=0.593;
+//Inteference Factor:
+a=1/3;
+//Density of air(in kg/m^3):
+p=1.23;
\ No newline at end of file diff --git a/746/DEPENDENCIES/11_01.sci b/746/DEPENDENCIES/11_01.sci new file mode 100755 index 000000000..8917647ce --- /dev/null +++ b/746/DEPENDENCIES/11_01.sci @@ -0,0 +1,19 @@ +//Temperature of air entering the cold section(in K):
+T1=440;
+//Absolute pressure of air entering the cold section(in kPa):
+p1=188;
+//Velocity of air entering the cold section(in m/sec):
+V1=210;
+//Temperature of air at outlet:(in K)
+T2=351;
+//Absolute pressure of air at outlet(in kPa):
+p2=213;
+//Rate of heat loss in the section(in kJ/sec):
+//Gas Constant(in N-m):
+R= 287;
+//Mass flow rate of air(in kg/sec):
+m=0.15;
+//Specific heat at constant pressue(in kJ/(kg-K)):
+cp=1;
+//Specific energy at constant volume(in kJ/(kg-K)):
+cv=0.717;
\ No newline at end of file diff --git a/746/DEPENDENCIES/11_03.sci b/746/DEPENDENCIES/11_03.sci new file mode 100755 index 000000000..89ce0534b --- /dev/null +++ b/746/DEPENDENCIES/11_03.sci @@ -0,0 +1,6 @@ +//Value of k:
+k=1.4;
+//Gas Constant(in Kj/(kg-K)):
+R=287;
+
+
diff --git a/746/DEPENDENCIES/11_04.sci b/746/DEPENDENCIES/11_04.sci new file mode 100755 index 000000000..deccf6ea9 --- /dev/null +++ b/746/DEPENDENCIES/11_04.sci @@ -0,0 +1,18 @@ +//Pressure at entry(in kPa):
+p1=350;
+//Temperature at entry(in K)
+T1=333;
+//Velocity at entry(in m/s):
+V1=183;
+//Mach no. at exit:
+M2=1.3;
+//Stagnation pressure at exit(in kPa):
+p02=385;
+//Stagnation temperature at exit(in K):
+T02=350;
+//Value of k:
+k=1.4;
+//Gas constant(in N-m/kg-K)
+R=287;
+//Specific heat at constant pressure(kJ/(kg-K):
+cp=1;
\ No newline at end of file diff --git a/746/DEPENDENCIES/12_01.sci b/746/DEPENDENCIES/12_01.sci new file mode 100755 index 000000000..aee1d5c02 --- /dev/null +++ b/746/DEPENDENCIES/12_01.sci @@ -0,0 +1,32 @@ +//Mach number at entry:
+M1=0.3;
+//Temperature at entry(in K):
+T1=335;
+//Pressure at entry(in kPa):
+p1=650;
+//Area at entry(in m^2):
+A1=0.001;
+//Mach number at exit:
+M2=0.8;
+///Value of k:
+k=1.4;
+//For the Mach no:0.3:
+//T/T0:
+t1=0.9823,
+//p/p0:
+P1=0.9395;
+//d/d0:
+den1=0.9564;
+//A/A*:
+a1=2.035;
+//For the Mach no:0.8:
+// T/T0:
+t2=0.8865;
+//p/p0:
+P2=0.6560;
+//d/d0:
+den2=0.7400;
+//A/A*:
+a2=1.038;
+//Gas Constant(in N-m/kg-K):
+R=287;
\ No newline at end of file diff --git a/746/DEPENDENCIES/12_02.sci b/746/DEPENDENCIES/12_02.sci new file mode 100755 index 000000000..20510fa79 --- /dev/null +++ b/746/DEPENDENCIES/12_02.sci @@ -0,0 +1,12 @@ +//Throat area of nozzle(in m^2):
+Ae=0.001;
+//Back pressure of air(in kPa):
+pb=591;
+//Stagnation pressure(in kPa):
+p0=1000;
+//Stagnation temperature(in K):
+T0=333;
+//Gas Constant(in N-m/kg-K):
+R=287;
+//Value of k:
+k=1.4;
diff --git a/746/DEPENDENCIES/12_03.sci b/746/DEPENDENCIES/12_03.sci new file mode 100755 index 000000000..faaa10025 --- /dev/null +++ b/746/DEPENDENCIES/12_03.sci @@ -0,0 +1,14 @@ +//Mach number at entry:
+M1=0.52;
+//Temperature at entry(in F):
+T1=40;
+//Pressure at entry(in psia):
+p1=60;
+//Area at entry(in ft^2):
+A1=0.013;
+//Back pressure(in psia):
+pb=30;
+//Gas Consant(in ft-lbf/lbm-R)
+R=53.3;
+//Value of k:
+k=1.4;
diff --git a/746/DEPENDENCIES/12_04.sci b/746/DEPENDENCIES/12_04.sci new file mode 100755 index 000000000..9fa20131a --- /dev/null +++ b/746/DEPENDENCIES/12_04.sci @@ -0,0 +1,14 @@ +//Stagnation temperature(in K):
+T0=350;
+//Stagnation pressure(in kPa):
+p0=1000;
+//Back Pressure(in kPa):
+pb=954;
+//Mach number at throat:
+Mt=0.68;
+//Area at exit(in m^2):
+Ae=0.001;
+//Value of k:
+k=1.4;
+//Gas Constant(in N-m/kg-K):
+R=287;
\ No newline at end of file diff --git a/746/DEPENDENCIES/12_05.sci b/746/DEPENDENCIES/12_05.sci new file mode 100755 index 000000000..9c90eefd1 --- /dev/null +++ b/746/DEPENDENCIES/12_05.sci @@ -0,0 +1,15 @@ +//Stagnation temperature(in K):
+T0=350;
+//Stagnation pressure(in kPa)\
+p0=1000;
+//Pressure at exit(in kPa)
+pe=87.5;
+//Back Pressure(in kPa):
+pb=50;
+//Area at exit(in m^2):
+Ae=0.001;
+//Gas Constant(in N-m/kg-K)
+R=287;
+//Value of k:
+k=1.4;
+
diff --git a/746/DEPENDENCIES/12_06.sci b/746/DEPENDENCIES/12_06.sci new file mode 100755 index 000000000..bde4e1e59 --- /dev/null +++ b/746/DEPENDENCIES/12_06.sci @@ -0,0 +1,14 @@ +//Diameter of pipe(in m):
+D=7.16*10^-3;
+//Stagnation pressure(in kPa):
+p0=101;
+//Stagnation temperature(in K):
+T0=296;
+//Pressure at section 1(in kPa):
+p1=98.5;
+//Temperature at section 2(in K):
+T2=287;
+//Gas constant(in N-m/kg-K):
+R=287;
+//Value of k:
+k=1.4;
diff --git a/746/DEPENDENCIES/12_07.sci b/746/DEPENDENCIES/12_07.sci new file mode 100755 index 000000000..030f1977a --- /dev/null +++ b/746/DEPENDENCIES/12_07.sci @@ -0,0 +1,22 @@ +//Stagnation temperature(in K):
+T0=296;
+//Stagnation pressure(in mm of Hg):
+p0=760;
+//Gauge pressure at section 1(in mm of Hg):
+p1=-18.9;
+//Gauge pressure at section2(in mm of Hg):
+p2=-412;
+//Mach number at 3:
+M3=1;
+//Gas constant:
+R=287;
+//Density of mercury(kg/m^3):
+dHg=13500;
+//Acceleration due to gravity(in m/sec^2):
+g=9.8;
+//Friction factor:
+f=0.0235;
+//Diameter of tube(in m):
+D=7.16*10^-3;
+//Value of k:
+k=1.4;
diff --git a/746/DEPENDENCIES/12_08.sci b/746/DEPENDENCIES/12_08.sci new file mode 100755 index 000000000..2e23c70f6 --- /dev/null +++ b/746/DEPENDENCIES/12_08.sci @@ -0,0 +1,19 @@ +//Temperature at section 1(in R):
+T1=600;
+//Pressure at section 1(in psia):
+p1=20;
+//Pressure at section 2(in psia):
+p2=10;
+//Velocity at section 1(in ft/s):
+V1=360;
+//Cross-sectional area of the duct(in ft^2):
+A=0.25;
+//Gas Constant(in ft-lbf/lbm-R):
+R=53.3;
+//Value of k:
+k=1.4;
+//Specific heat at constant presure(in Btu/lbm-R):
+Cp=0.24;
+//Specific heat at constant volume(in Btu/lbm-R):
+Cv=0.171;
+
diff --git a/746/DEPENDENCIES/12_09.sci b/746/DEPENDENCIES/12_09.sci new file mode 100755 index 000000000..f2d16264f --- /dev/null +++ b/746/DEPENDENCIES/12_09.sci @@ -0,0 +1,14 @@ +//Temperature at section 1(in K):
+T1=333;
+//Pressure at section 1(in kPa):
+p1=135;
+//Velocity at section 1(in m/sec):
+V1=732;
+//Mach number at section 2:
+M2=1.2;
+//Gas constant(in N-m/kg-K):
+R=287;
+//Value of k:
+k=1.4;
+//Specific heat at constant presure(in N-m/kg-K):
+Cp=1;
\ No newline at end of file diff --git a/746/DEPENDENCIES/12_10.sci b/746/DEPENDENCIES/12_10.sci new file mode 100755 index 000000000..aab81766f --- /dev/null +++ b/746/DEPENDENCIES/12_10.sci @@ -0,0 +1,10 @@ +//Temperature at section 1(in K):
+T1=278;
+//Pressure at section 1(in kPa):
+p1=65;
+//Velocity at section 1(in m/sec):
+V1=668;
+//Value of gas constant(in N-m/kg-K):
+R=287;
+//Value of k:
+k=1.4;
diff --git a/746/DEPENDENCIES/1_01.sci b/746/DEPENDENCIES/1_01.sci new file mode 100755 index 000000000..20b290d7e --- /dev/null +++ b/746/DEPENDENCIES/1_01.sci @@ -0,0 +1,10 @@ +//Mass of oxygen present(in kg):
+m=0.95;
+//Initial temperatur(in K):
+T1=300;
+//Final temperature of oxygen(in K):
+T2=900;
+//Pressure of oxygen(in kPa):
+p=150;
+//Specific heat at constant pressure(in J/kg-K):
+cp=909.4;
diff --git a/746/DEPENDENCIES/1_02.sci b/746/DEPENDENCIES/1_02.sci new file mode 100755 index 000000000..013a4b783 --- /dev/null +++ b/746/DEPENDENCIES/1_02.sci @@ -0,0 +1,8 @@ +//Mass of ball(in kg):
+m=0.2;
+//Height fom which ball is dropped(in m):
+y0=500;
+//Value of k:
+k=2*10^-4;
+//Accleration due to gravity(in m/sec^2):
+g=9.81;
\ No newline at end of file diff --git a/746/DEPENDENCIES/2_02.sci b/746/DEPENDENCIES/2_02.sci new file mode 100755 index 000000000..bdb2db7dd --- /dev/null +++ b/746/DEPENDENCIES/2_02.sci @@ -0,0 +1,11 @@ +//Liquid Viscosity(in cp):
+u=0.65;
+//Specific gravity:
+SG=0.88;
+//Density of water(in slug/ft^3):
+d=1.94;
+//Velocity with which plate is moved(in m/sec):
+U=0.3;
+//Distance between the plates(in mm):
+D=0.3;
+
diff --git a/746/DEPENDENCIES/3_01.sci b/746/DEPENDENCIES/3_01.sci new file mode 100755 index 000000000..903102918 --- /dev/null +++ b/746/DEPENDENCIES/3_01.sci @@ -0,0 +1,14 @@ +//Surface tension of water(in mN/m):
+STw=72.8*10^-3;
+//Surface Tension of mercury(in mN/m):
+STm=375*10^-3;
+//Contact angle for water:
+thetaw=0;
+//COntact angle for mercury:
+thetam=140;
+//Density of water(in kg/m^3):
+dw=1;
+//Density of mercury(in kg/m^3):
+dm=13.6;
+//Acceleration de to gravity(in m/sec):
+g=9.81;
diff --git a/746/DEPENDENCIES/3_03.sci b/746/DEPENDENCIES/3_03.sci new file mode 100755 index 000000000..c2fd37685 --- /dev/null +++ b/746/DEPENDENCIES/3_03.sci @@ -0,0 +1,16 @@ +//Acceleration due to gravity(in ft/sec^2):
+g=32.2;
+//Specific gravity of mercury:
+SGm=13.6;
+//Specific gravity of oil:
+SGo=0.88;
+//Specific gravity of water:
+SGw=1;
+//Density of water(in slug/ft^3):
+d=1.94;
+//Heights of liquid in various tubes(in inches):
+d1=10;
+d2=3;
+d3=4;
+d4=5;
+d5=8;
diff --git a/746/DEPENDENCIES/3_04.sci b/746/DEPENDENCIES/3_04.sci new file mode 100755 index 000000000..61f2a35d0 --- /dev/null +++ b/746/DEPENDENCIES/3_04.sci @@ -0,0 +1,17 @@ +//Elevation of Denver(in ft):
+z1=5280;
+//Pressure at Denver(in mm of Hg):
+p1=24.8;
+//Temperature at Denver(in F):
+T1=80;
+//Elevation at Vail Pass(in ft):
+z2=10600;
+//Temperature at Vsil Pass(in F):
+T2=62;
+//Value of R in ft-lbf/lbm-R):
+R=53.3;
+//Acceleration due togravity(in ft/sec^2):
+g=32.2;
+//Value of adiabatic constant:
+k=1.4;
+
diff --git a/746/DEPENDENCIES/3_05.sci b/746/DEPENDENCIES/3_05.sci new file mode 100755 index 000000000..4388b9227 --- /dev/null +++ b/746/DEPENDENCIES/3_05.sci @@ -0,0 +1,12 @@ +//Length of gate(in m):
+L=4;
+//Width of gate(in m):
+w=5;
+//Depth of gate under water(in m):
+D=2;
+//Density of water(in kg/m^3:
+d=999;
+//Acceleration due to gravity(in m/sec^2):
+g=9.81;
+//Angle of gate with horizontal:
+theta=30;
diff --git a/746/DEPENDENCIES/3_06.sci b/746/DEPENDENCIES/3_06.sci new file mode 100755 index 000000000..2bfa3f6c0 --- /dev/null +++ b/746/DEPENDENCIES/3_06.sci @@ -0,0 +1,9 @@ +//Pressure apllied on the door(in psfg):
+p0=100;
+//Length of door(in feet):
+L=3;
+//Breadth of the door(in feet):
+b=2;
+//Density of liqiuid(in lbf/ft^3):
+d=100;
+
diff --git a/746/DEPENDENCIES/3_07.sci b/746/DEPENDENCIES/3_07.sci new file mode 100755 index 000000000..0a44b164a --- /dev/null +++ b/746/DEPENDENCIES/3_07.sci @@ -0,0 +1,12 @@ +//Width of gate(in m):
+w=5;
+//Depth of water(in m):
+D=4;
+//Density of water(in kg/m^3);
+d=999;
+//Accelration deto gravity(in m/sec^2):
+g=9.81;
+//Value of a (in m):
+a=4;
+//Point where force acts(in m):
+l=5;
diff --git a/746/DEPENDENCIES/4_01.sci b/746/DEPENDENCIES/4_01.sci new file mode 100755 index 000000000..fee4e0aa4 --- /dev/null +++ b/746/DEPENDENCIES/4_01.sci @@ -0,0 +1,16 @@ +//Area of 1 (in ft^2):
+A1=0.2;
+//Area of 2 (in ft^2):
+A2=0.5;
+//Area of 3 (in ft^2):
+A3=0.4;
+//Area of 4 (in ft^2):
+A4=0.4;
+//Density of water (in slug/ft^3):
+d=1.94;
+//Mass flow rate out of section 3(in slug/sec):
+m3=3.88;
+//Volme flow rate in section 4 (in ft^3/sec):
+Q4=1;
+//Velocity at 1(in ft/sec):
+V1=10;
diff --git a/746/DEPENDENCIES/4_02.sci b/746/DEPENDENCIES/4_02.sci new file mode 100755 index 000000000..d5e051ffe --- /dev/null +++ b/746/DEPENDENCIES/4_02.sci @@ -0,0 +1,8 @@ +//Flow velocity ahead of the plate(in m/sec):
+U=30;
+//Boundary layer tckness at location d(in mm):
+t=5;
+//Density of fluid air(in k/m^3):
+d=1.24;
+//Plate wdth perpendicular to the plate(in m):
+w=0.6;
diff --git a/746/DEPENDENCIES/4_03.sci b/746/DEPENDENCIES/4_03.sci new file mode 100755 index 000000000..fe4d65c83 --- /dev/null +++ b/746/DEPENDENCIES/4_03.sci @@ -0,0 +1,12 @@ +//Volume of tank(in m^3):
+V=0.05;
+//Pressure of air(In kPa):
+p=800;
+//Temperature of tank(in C):
+T=15;
+//Velocity of leavig air(in m/sec):
+v=311;
+//Density of air(in kg/m^3):
+d=6.13;
+//Area ofvalve exit(in mm^2):
+A=65;
diff --git a/746/DEPENDENCIES/4_04.sci b/746/DEPENDENCIES/4_04.sci new file mode 100755 index 000000000..e0fad90ee --- /dev/null +++ b/746/DEPENDENCIES/4_04.sci @@ -0,0 +1,6 @@ +//Velocity of water leaving the nozle(in m/sec):
+V=15;
+//Area of nozzle(in m^2):
+A=0.01;
+//Density of water(in kg/m^3):
+d=999;
diff --git a/746/DEPENDENCIES/4_05.sci b/746/DEPENDENCIES/4_05.sci new file mode 100755 index 000000000..d5a532a75 --- /dev/null +++ b/746/DEPENDENCIES/4_05.sci @@ -0,0 +1,18 @@ +//Height of the container(in ft):
+l=2;
+//Area of cross section(in ft^2):
+A=1;
+//Weight of container(in lbf):
+W=5;
+//Water depth (in ft):
+h=1.9;
+//Area of opening 1(in ft^2):
+A1=0.1;
+//Velocity at opening 1(in ft/sec):
+V1=-5;
+//Area of opening 2(in ft^2):
+A2=0.1;
+//Area of opening 1(in ft^2):
+A3=0.1;
+//Density of water(in slug/f^3):
+d2=1.94;
diff --git a/746/DEPENDENCIES/4_06.sci b/746/DEPENDENCIES/4_06.sci new file mode 100755 index 000000000..1546d639d --- /dev/null +++ b/746/DEPENDENCIES/4_06.sci @@ -0,0 +1,13 @@ +//Diameter of channel(in m):
+D1=1.5;
+//Velcity of flow in channel(in m/sec):
+V1=0.2;
+//Diameter at section 2(in m):
+D2=0.0563;
+//Velocity a section 2(in m/sec):
+V2=5.33;
+//Density of water(in kg/m^3):
+d=999;
+//Acceleration due to gravity(in m/sec2):
+g=9.81;
+
diff --git a/746/DEPENDENCIES/4_07.sci b/746/DEPENDENCIES/4_07.sci new file mode 100755 index 000000000..33b028eac --- /dev/null +++ b/746/DEPENDENCIES/4_07.sci @@ -0,0 +1,10 @@ +//Pressure at inlet tothe elbow(in N/m^2):
+p1=2.21*10^5;
+//Area of crosssection(in m^2):
+A1=0.01;
+//Velocity at secton 2(in m/sec):
+V2=16;
+//Area of cross section of section 2(in m^2):
+A2=0.0025;
+//Atmospheric pressure(in kPa):
+patm=1.012*10^5;
diff --git a/746/DEPENDENCIES/4_08.sci b/746/DEPENDENCIES/4_08.sci new file mode 100755 index 000000000..bcbb1f39b --- /dev/null +++ b/746/DEPENDENCIES/4_08.sci @@ -0,0 +1,7 @@ +//Velocity of conveyor belt(in ft/sec):
+Vbelt=3;
+//Velocity of sand alling onto belt(in ft/sec):
+Vsand=5;
+//Flow rate(in lbm/sec):
+m=500;
+
\ No newline at end of file diff --git a/746/DEPENDENCIES/4_09.sci b/746/DEPENDENCIES/4_09.sci new file mode 100755 index 000000000..9cc7ed26f --- /dev/null +++ b/746/DEPENDENCIES/4_09.sci @@ -0,0 +1,8 @@ +//Nozzle inlet diameter(in inchess):
+D1=3;
+//Nozzle exit diameter(in inches):
+D2=1;
+//Desired volume flow rate(in ft^3/sec):
+Q=0.7;
+//Density of water(in slug/ft^3):
+d=1.94;
diff --git a/746/DEPENDENCIES/4_10.sci b/746/DEPENDENCIES/4_10.sci new file mode 100755 index 000000000..822dc60b4 --- /dev/null +++ b/746/DEPENDENCIES/4_10.sci @@ -0,0 +1,10 @@ +//Vane turning angle:
+theta=60;
+//Speed of vane(in m/sec):
+U=10;
+//Area of nozzle(in m2):
+A=0.003;
+//Flow velocity of water(in m/sec):
+V=30;
+//Density of water(in kg/m^3):
+d=999;
diff --git a/746/DEPENDENCIES/4_11.sci b/746/DEPENDENCIES/4_11.sci new file mode 100755 index 000000000..4fe857acb --- /dev/null +++ b/746/DEPENDENCIES/4_11.sci @@ -0,0 +1,10 @@ +//Mass of vane and cart(in kg):
+M=75;
+//Turning angle of vane:
+theta=60;
+//Speed of water leaving nozzle horizontally(in m/sec):
+V=35;
+//Exit area of nozzle(in m^):
+A=0.003;
+//Density of water(in kg/m^3):
+d=999;
diff --git a/746/DEPENDENCIES/4_12.sci b/746/DEPENDENCIES/4_12.sci new file mode 100755 index 000000000..1137b7719 --- /dev/null +++ b/746/DEPENDENCIES/4_12.sci @@ -0,0 +1,10 @@ +//Initial mass of th rocket(in kg):
+M0=400;
+//Rate of fuel consumption(in kg/sec):
+me=5;
+//Exhaust velocity(in m/sec):
+Ve=3500;
+//Acceleration due to gravity(in m/sec^2):
+g=9.81;
+//Time after which velocity is to be calculated(in sec):
+t=10;
diff --git a/746/DEPENDENCIES/4_14.sci b/746/DEPENDENCIES/4_14.sci new file mode 100755 index 000000000..e06524131 --- /dev/null +++ b/746/DEPENDENCIES/4_14.sci @@ -0,0 +1,17 @@ +//Inlet gauge pressure(in kPa):
+p=20;
+//Volume flow rate of water through the sprinkler(in l/min):
+Q=7.5;
+//Speed of rotstion of sprinkler(in rpm):
+w=30;
+//Diameter of jet f sprinkle(in mm):
+D=4;
+//Radius of sprinkler(in mm):
+R=150;
+//Supply pressure to sprinkler(in kPa):
+p=20;
+//Angle at which jet is sprayed wrt horizontal:
+alpha=30;
+//Density of water(in kg/m^):
+d=999;
+
diff --git a/746/DEPENDENCIES/4_16.sci b/746/DEPENDENCIES/4_16.sci new file mode 100755 index 000000000..3bc2cc9a1 --- /dev/null +++ b/746/DEPENDENCIES/4_16.sci @@ -0,0 +1,19 @@ +//Pressure at entry(in psia):
+p1=14.7;
+//Temperature at entry(in F):
+T1=70;
+//Pressure at exit(in psia):
+p2=50;
+//Temprature a exit(in F):
+T2=100;
+//Cross sectional area of the pipe at exit(in ft^2):
+A2=1;
+//Mass flow rate(in lbf/sec):
+m=20;
+//Power input to the compressor(in hp):
+Ws=-600;
+//Value of cp(in Btu/lbm-R):
+cp=0.24;
+//Value of gas constant(in ft-lbf/(lbm-R))
+R=53.3;
+
diff --git a/746/DEPENDENCIES/4_17.sci b/746/DEPENDENCIES/4_17.sci new file mode 100755 index 000000000..0a7be40c9 --- /dev/null +++ b/746/DEPENDENCIES/4_17.sci @@ -0,0 +1,17 @@ +//Volume of tak(in m^3):
+V=0.1;
+//Temperature of line and tank(in K):
+T=293;
+//Initial tank gauge pressure(in N/m^2):
+p1=1*10^5;
+//Absolute line pressure(in N/m^2):
+p=2*10^6;
+//Rate of rise of temperature after opening of the valve(in C/sec):
+r=0.05;
+//Atmospheric pressure(in N/m^2):
+patm=1.01*10^5;
+//Gas Constant(in N-m/(kg-K)):
+R=287;
+//Value of cv(in N-m/kg-K):
+cv=717;
+
diff --git a/746/DEPENDENCIES/5_02.sci b/746/DEPENDENCIES/5_02.sci new file mode 100755 index 000000000..f0eef8b04 --- /dev/null +++ b/746/DEPENDENCIES/5_02.sci @@ -0,0 +1,6 @@ +//Distance f piston from closed end of the cylinder at the give instant(in m):
+L=0.15;
+//Density of gas(in kg/m^3):
+d=18;
+//Velocity of piston(in m/sec):
+V=12;
diff --git a/746/DEPENDENCIES/5_07.sci b/746/DEPENDENCIES/5_07.sci new file mode 100755 index 000000000..c59a89585 --- /dev/null +++ b/746/DEPENDENCIES/5_07.sci @@ -0,0 +1,7 @@ +//Value of (in mm/sec):
+U=4;
+//Value of h(in mm):
+h=4;
+//Tme at which to find position(in sec):
+t=1.5;
+
diff --git a/746/DEPENDENCIES/5_08.sci b/746/DEPENDENCIES/5_08.sci new file mode 100755 index 000000000..c44082b87 --- /dev/null +++ b/746/DEPENDENCIES/5_08.sci @@ -0,0 +1,2 @@ +//Value of A(in sec^-1):
+A=0.3;
diff --git a/746/DEPENDENCIES/5_09.sci b/746/DEPENDENCIES/5_09.sci new file mode 100755 index 000000000..38a19d43b --- /dev/null +++ b/746/DEPENDENCIES/5_09.sci @@ -0,0 +1,12 @@ +//Thickness of water film(in mm):
+h=1;
+//Width of surface(in m):
+b=1;
+//Angle of inclination of surface:
+theta=15;
+//Density of water(in kg/m^3):
+d=999;
+//Acceleration du to gravity(in m/sec^2):
+g=9.81;
+//Viscosity(kg/m-s):
+u=10^-3;
\ No newline at end of file diff --git a/746/DEPENDENCIES/6_01.sci b/746/DEPENDENCIES/6_01.sci new file mode 100755 index 000000000..c2acec12b --- /dev/null +++ b/746/DEPENDENCIES/6_01.sci @@ -0,0 +1,15 @@ +//Depth of the duct(in m):
+d=0.3;
+//Width of the duct(in m):
+w=0.1;
+//Inner radius of the bend(in m):
+r=0.25;
+//Pressure difference between the taps(in mm of Hg):
+p=40;
+//Density of water(in kg/m^3):
+dw=999;
+//Acceleration due to gravity(in m/sec^2):
+g=9.8;
+//Density of air(in kg/m^3):
+da=1.23;
+
diff --git a/746/DEPENDENCIES/6_02.sci b/746/DEPENDENCIES/6_02.sci new file mode 100755 index 000000000..ccbd18b0f --- /dev/null +++ b/746/DEPENDENCIES/6_02.sci @@ -0,0 +1,10 @@ +//Pressure diference(in mm of mecury):
+p=30;
+//Density of water(in kg/m^3):
+dw=1000;
+//Aceleration due to gravity(in m/sec^2):
+g=9.81;
+//Density of air(in kg/m^3):
+da=1.23;
+//Specific gravity of mercury:
+SG=13.6;
diff --git a/746/DEPENDENCIES/6_03.sci b/746/DEPENDENCIES/6_03.sci new file mode 100755 index 000000000..2384ded21 --- /dev/null +++ b/746/DEPENDENCIES/6_03.sci @@ -0,0 +1,8 @@ +//Area of nozzle at input(in m^2):
+Ai=0.1;
+//Area of nozzle at exit(in m^2):
+Ae=0.02;
+//Outlet velocity of flow(in m/sec):
+V2=50;
+//Density of air(in kg/m^3):
+da=1.23;
diff --git a/746/DEPENDENCIES/6_04.sci b/746/DEPENDENCIES/6_04.sci new file mode 100755 index 000000000..8e11ea28e --- /dev/null +++ b/746/DEPENDENCIES/6_04.sci @@ -0,0 +1,12 @@ +//Length of tube above surface(in m):
+l=1;
+//Depth of exit below water surface(in m):
+z=7;
+//Acceleration due to gravity(in m/sec^2):
+g=9.81;
+//Density of water(in kg/m^3):
+d=999;
+//Atmospheric pressure(in N/m^2):
+p1=1.01*10^5;
+
+
diff --git a/746/DEPENDENCIES/6_05.sci b/746/DEPENDENCIES/6_05.sci new file mode 100755 index 000000000..9def2ba8e --- /dev/null +++ b/746/DEPENDENCIES/6_05.sci @@ -0,0 +1,6 @@ +//Depth of water at the upstream(on feet):
+Du=1.5;
+//Depth of water at the vena contracta downstream from the gate(in inches):
+Dd=2;
+//Acceleration due to gravity(in ft/sec^2):
+g=32.2;
diff --git a/746/DEPENDENCIES/6_06.sci b/746/DEPENDENCIES/6_06.sci new file mode 100755 index 000000000..6a0015e66 --- /dev/null +++ b/746/DEPENDENCIES/6_06.sci @@ -0,0 +1,14 @@ +//Speed of plane(in km/hr):
+V=150;
+//Speed at point B relative to the wing(in m/sec):
+Vb=60;
+//Density of air(in kg/m^3):
+da=1.23;
+//Atmospheris pressure(in N/m^2):
+pa=1.01*10^5;
+//At 1000m,
+//p/pSL:
+P1=0.8870;
+//d/dSL:
+D1=0.9075;
+
diff --git a/746/DEPENDENCIES/6_08.sci b/746/DEPENDENCIES/6_08.sci new file mode 100755 index 000000000..1f2dfbe41 --- /dev/null +++ b/746/DEPENDENCIES/6_08.sci @@ -0,0 +1,10 @@ +//Area of cross section of the nozzle(in in^2):
+A4=0.864;
+//Capacity of heater(in kW):
+Q=10
+//Acceleration due to gravity(in ft/sec^2):
+g=32.2;
+//Water level in reservoir above datum line(in ft):
+z3=10;
+//Density of water(In slug/ft^3):
+d=1.94;
diff --git a/746/DEPENDENCIES/6_09.sci b/746/DEPENDENCIES/6_09.sci new file mode 100755 index 000000000..c852adc49 --- /dev/null +++ b/746/DEPENDENCIES/6_09.sci @@ -0,0 +1,8 @@ +//Depth to which water is filled(in m):
+h=3;
+//Length of pipe(in m):
+L=6;
+//Diameter of pipe (in mm):
+D=150;
+//Acceleration due to gravity(in m/sec^2):
+g=9.81;
diff --git a/746/DEPENDENCIES/7_04.sci b/746/DEPENDENCIES/7_04.sci new file mode 100755 index 000000000..c7fe1e512 --- /dev/null +++ b/746/DEPENDENCIES/7_04.sci @@ -0,0 +1,16 @@ +//Diameter of the prototype(in ft):
+Dp=1;
+//Speed of towing of prototype(in knots):
+Vp=5;
+//Diameter of model(in inches):
+Dm=6;
+//Drag for model at test condition(in lbf):
+Fm=5.58;
+//Density of seawater at 5 C for prototype(in slug/ft^3):
+dp=1.99;
+//Kinematic viscosity at 5 C for prototype(in ft^2/sec):
+vp=1.69*10^-5;
+//Density of air at STP for model(in slug/ft^3):
+dm=0.00238;
+//Kinematic viscosity of air at STP for model(in ft^2/sec):
+vm=1.57*10^-4;
diff --git a/746/DEPENDENCIES/7_05.sci b/746/DEPENDENCIES/7_05.sci new file mode 100755 index 000000000..a3990f7ad --- /dev/null +++ b/746/DEPENDENCIES/7_05.sci @@ -0,0 +1,18 @@ +//Width of the prototype(in ft):
+wp=8;
+//Frontal area of the prototype(in ft^2):
+Ap=84;
+//Model Scale:
+S=1/16;
+//Density of air(in kg/m^3):
+d=1.23;
+//Air speed in wind tunnel(in m/sec):
+V=[18 21.8 26 30.1 35 38.5 40.9 44.1 46.7];
+//Drag force(in N):
+Fd=[3.1 4.41 6.09 7.97 10.7 12.9 14.7 16.9 18.9];
+//Kinematic viscosity(in m^2/sec):
+v=1.46*10^-5;
+//Density of air(in kg/m^3):
+d=1.23;
+//Speed of prototype(in km/hr):\
+Vp=100;
diff --git a/746/DEPENDENCIES/7_06.sci b/746/DEPENDENCIES/7_06.sci new file mode 100755 index 000000000..34626f0e4 --- /dev/null +++ b/746/DEPENDENCIES/7_06.sci @@ -0,0 +1,16 @@ +//Efficinc of pump:
+Effp=0.8;
+//Design specific speed(in rpm):
+Nscu1=2000;
+//Impeller diameter(in inches):
+D1=8;
+//Opertion sped at esign point flow condition(in rpm):
+N1=1170;
+//Flow rate at design point flow condition(in gpm):
+Q1=300;
+//Density of water (in slug/ft^3):
+d1=1.94;
+//Acceleration due to gravity(in ft^2/sec):
+g=32.2;
+//Working speed 2(in rpm):
+N2=1750;
diff --git a/746/DEPENDENCIES/8_01.sci b/746/DEPENDENCIES/8_01.sci new file mode 100755 index 000000000..2d4620b38 --- /dev/null +++ b/746/DEPENDENCIES/8_01.sci @@ -0,0 +1,16 @@ +//Operation pressure of hydraulic system(in kPa):
+p1=20000;
+//Operation temperature of hydraulic system(in C):
+T=55;
+//Piston diameter(in mm):
+D=25;
+//Viscosity of SAE 10W at 55C(in kg/(m-s):
+u=0.018;
+//Mean radial clearance of a cylinder(in mm):
+a=0.005;
+//Gauge pressure on lower pressure side of piston(in kPa):
+p2=1000;
+//Lenth of piston(in mm):
+L=15;
+//Denity of water(in kg/m^3):
+dw=1000;
\ No newline at end of file diff --git a/746/DEPENDENCIES/8_02.sci b/746/DEPENDENCIES/8_02.sci new file mode 100755 index 000000000..fa2011727 --- /dev/null +++ b/746/DEPENDENCIES/8_02.sci @@ -0,0 +1,16 @@ +//temperature fo operation(in F):
+T=210;
+//Diameter of te bearing(in inches):
+D=3;
+//Diametral clearance(in inches):
+a=0.0025;
+//Length of shaft(in inhes):
+L=1.25;
+//Speed of rotation of the shaft(in rpm):
+N=3600;
+//Viscosity of the oil(in lbf-s/ft^2):
+u=2.01*10^-4;
+//Specific gravity of SAE 10W:
+SG=0.92;
+//Density of water (in slug/ft^3)
+p=1.94;
\ No newline at end of file diff --git a/746/DEPENDENCIES/8_04.sci b/746/DEPENDENCIES/8_04.sci new file mode 100755 index 000000000..230ba6ea5 --- /dev/null +++ b/746/DEPENDENCIES/8_04.sci @@ -0,0 +1,10 @@ +//Flow rate through capilarry viscometer(in mm^3/sec):
+Q=880;
+//Tube length(in m):
+L=1;
+//Tube diameter(in mm):
+D=0.5;
+//Pressure drop(in kPa):
+p=1000;
+//Density of oil(in kg/m^3):
+d=999;
diff --git a/746/DEPENDENCIES/8_05.sci b/746/DEPENDENCIES/8_05.sci new file mode 100755 index 000000000..50ad20760 --- /dev/null +++ b/746/DEPENDENCIES/8_05.sci @@ -0,0 +1,16 @@ +//Volme flow rate of water(in m^3/sec):
+Q=0.0084;
+//Length of horizontal pipe(in m):
+L=100;
+//Diameter of pipe(in m):
+D=0.075;
+//Density of water(in kg/m^3):
+d=999;
+//Friction factor:
+f=0.017;
+//Minor lossses coefficient:
+K=0.5;
+//Viscosity(in kg/m-s):
+u=10^-3;
+//Acceleration due to gravity(in /sec^2):
+g=9.8;
\ No newline at end of file diff --git a/746/DEPENDENCIES/8_06.sci b/746/DEPENDENCIES/8_06.sci new file mode 100755 index 000000000..87e0e9dda --- /dev/null +++ b/746/DEPENDENCIES/8_06.sci @@ -0,0 +1,20 @@ +//Flow rate of crude oil(in bbl):
+Q=1.6*10^6;
+//Inside diamete of pipe(i inches):
+D=48;
+//Maximum allowable pressure(in psi):
+p2=1200;
+//Minimum pressure required to keep gases dissolves(in psi):
+p1=50;
+//Specific gravity of crde oil:
+SG=0.93;
+//Viscosity at 140 F(in lbf-s/ft^2):
+u=3.5*10^-4;
+//Efficincy of pump:
+Effp=0.85;
+//Density(in slug/ft^3):
+d=1.94;
+//Viscosity (in lbf-sec):
+u=3.5*10^-4;
+//Friction factor:
+f=0.017;
diff --git a/746/DEPENDENCIES/8_07.sci b/746/DEPENDENCIES/8_07.sci new file mode 100755 index 000000000..915a806e7 --- /dev/null +++ b/746/DEPENDENCIES/8_07.sci @@ -0,0 +1,10 @@ +//Height of standpipe(in ft):
+l=80;
+//Length of longest pipe(in ft):
+L=600;
+//Diameter of pipe(in inches):
+D=4;
+//Friction factor:
+f=0.031;
+//Acceleration due to gravityin ft/sec^2):
+g=32.2;
diff --git a/746/DEPENDENCIES/8_08.sci b/746/DEPENDENCIES/8_08.sci new file mode 100755 index 000000000..7682c363f --- /dev/null +++ b/746/DEPENDENCIES/8_08.sci @@ -0,0 +1,12 @@ +//Length of Al tubing (in ft):
+L=500;
+//Volume flow rate of pump output(in gpm):
+Q=1500;
+//Discharge pressure(in psig):
+p1=65;
+//Sprinkler pressure(in psig):
+p2=30;
+//Kinematic viscosity(in ft^2/sec):
+v=1.21*10^-5;
+//Density(in slug/ft^3):
+p=1.94;
diff --git a/746/DEPENDENCIES/8_09.sci b/746/DEPENDENCIES/8_09.sci new file mode 100755 index 000000000..afd919d32 --- /dev/null +++ b/746/DEPENDENCIES/8_09.sci @@ -0,0 +1,12 @@ +//Length of copper wire(in ft):
+L=10;
+//Inner diameer of pipe(in inches):
+D=1.5;
+//Dischare(in ft^3/sec):
+Q=0.566;
+//Level of reservoir above pipe centreline(inn feet):
+h=85.1;
+//Kinematic viscosity at 70 F(in ft^2/s):
+v=1.05*10^-5;
+//Acceleration due to gravity(in ft/sec^2):
+g=32.2;
\ No newline at end of file diff --git a/746/DEPENDENCIES/8_10.sci b/746/DEPENDENCIES/8_10.sci new file mode 100755 index 000000000..4d8f2fbb0 --- /dev/null +++ b/746/DEPENDENCIES/8_10.sci @@ -0,0 +1,12 @@ +//Nozzle exit diameter(in mm):
+D=25;
+//N/R1 value value:
+N_R=3;
+//AR value:
+A_R=2;
+//Static head available from the main(in m):
+z0=1.5;
+//Acceleration due to gravity(in m/sec^2):
+g=9.8;
+//Value of Cp:
+Cp=0.45;
diff --git a/746/DEPENDENCIES/8_11.sci b/746/DEPENDENCIES/8_11.sci new file mode 100755 index 000000000..33248b034 --- /dev/null +++ b/746/DEPENDENCIES/8_11.sci @@ -0,0 +1,15 @@ +//Volume flw rate of ai(in m^3/sec):
+Q=1;
+//Diameter of pipe(in m):
+D=0.25;
+//Density of air (in kg/m^3):
+d1=1.23;
+//Acceleration due to gravity(in m/s^2):
+g=9.8;
+//Density of water(in kg/m^3):
+d2=999;
+//Maxmum range of manometer(in m):
+h=0.3;
+//Kinematic viscosity(in m^2/s):
+v=1.46*10^-5;
+
diff --git a/746/DEPENDENCIES/9_01.sci b/746/DEPENDENCIES/9_01.sci new file mode 100755 index 000000000..ac62714ef --- /dev/null +++ b/746/DEPENDENCIES/9_01.sci @@ -0,0 +1,8 @@ +//Lengh of side of the test section(in mm):
+L=305;
+//Freesteam speed at section 1(in m/sec):
+U1=26;
+//Displacement thickness at section 1(in mm):
+d1=1.5;
+//Displacment thickness at section 2(in mm):
+d2=2.1;
diff --git a/746/DEPENDENCIES/9_04.sci b/746/DEPENDENCIES/9_04.sci new file mode 100755 index 000000000..05ba26f5a --- /dev/null +++ b/746/DEPENDENCIES/9_04.sci @@ -0,0 +1,8 @@ +//Veocity of flow(in m/sec):
+U=1;
+//Length of flat plate(in m):
+L=1;
+//Density of water(in kg/m^3):
+d=999;
+//Kinematic viscosity of water(in m^2/sec):
+v=10^-6;
diff --git a/746/DEPENDENCIES/9_05.sci b/746/DEPENDENCIES/9_05.sci new file mode 100755 index 000000000..80abff69f --- /dev/null +++ b/746/DEPENDENCIES/9_05.sci @@ -0,0 +1,12 @@ +//Length of the supertanker(in m):
+L=360;
+//Width of supertanker(in m):
+W=70;
+//Draft of the supertanker(in m):
+D=50;
+//Cruising speed in water(in knots):
+s=13;
+//Kinematic viscosity at 10 C
+v=1.37*10^-6;
+//Density of sea water(in kg/m^3):
+d=1020;
diff --git a/746/DEPENDENCIES/9_06.sci b/746/DEPENDENCIES/9_06.sci new file mode 100755 index 000000000..40b30d6a4 --- /dev/null +++ b/746/DEPENDENCIES/9_06.sci @@ -0,0 +1,13 @@ +//Diameter of chimney(in m):
+D=1;
+//Height of chimney(in m):
+L=25;
+//Speed of wind(in kmph):
+s=50;
+//Density of air(in kg/m^3):
+d=1.23;
+//Viscosity of air(in kg/(m-s)):
+u=1.79*10^-5;
+//Pressure(in kPa):
+p=101;
+
diff --git a/746/DEPENDENCIES/9_07.sci b/746/DEPENDENCIES/9_07.sci new file mode 100755 index 000000000..a6add1b08 --- /dev/null +++ b/746/DEPENDENCIES/9_07.sci @@ -0,0 +1,14 @@ +//Weight of the dragster(in lbf):
+w=1600;
+//Speed of dragster(in mph):
+s1=270;
+//Area of drag chute(in ft^2):
+A=25;
+//Speed of dragster after deceleraton(in mph):
+s2=100;
+//Acceleration due to gravity(in ft/sec^2):
+g=32.2;
+//Density of air(in slug/ft^3):
+d=0.00238;
+//Value of coefficient of drag:
+Cd=1.42;
diff --git a/746/DEPENDENCIES/9_08.sci b/746/DEPENDENCIES/9_08.sci new file mode 100755 index 000000000..a6004efa4 --- /dev/null +++ b/746/DEPENDENCIES/9_08.sci @@ -0,0 +1,15 @@ +//Wing area(in ft^2):
+A=1600;
+//Aspect ratio:
+ar=6.5;
+//Groos weight of aircraft(in lbf):
+W=150000;
+//Coefficient of drag at zero lift :
+Cd0=0.0182;
+//Sonic speed at sea level(in mph):
+c=759;
+//Density of air(in slug/ft^3):
+p=0.00238;
+//Srall speed at sea level(in mph):
+Vssl=175;
+
diff --git a/746/DEPENDENCIES/9_09.sci b/746/DEPENDENCIES/9_09.sci new file mode 100755 index 000000000..3c491ff54 --- /dev/null +++ b/746/DEPENDENCIES/9_09.sci @@ -0,0 +1,14 @@ +//Mass of the tennis ball(in grams):
+m=57;
+//Diameter ofthe ball (in mm):
+D=64;
+//Velocity with which te ball is hit(in m/s):\
+V=25;
+//Topspin given on the ball(in rpm):
+N=7500;
+//Acceleration due to gravity(in m/s^2):
+g=9.81;
+//Kinematic viscosity(in m^2/s):
+v=1.46*10^-5
+//Desity of air(in kg/m^3):
+d=1.23;
|
