diff options
Diffstat (limited to '914/CH9')
-rwxr-xr-x | 914/CH9/EX9.3/ex9_3.sce | 32 | ||||
-rwxr-xr-x | 914/CH9/EX9.4/ex9_4.sce | 25 | ||||
-rwxr-xr-x | 914/CH9/EX9.5/ex9_5.sce | 17 |
3 files changed, 74 insertions, 0 deletions
diff --git a/914/CH9/EX9.3/ex9_3.sce b/914/CH9/EX9.3/ex9_3.sce new file mode 100755 index 000000000..1f2ffcef9 --- /dev/null +++ b/914/CH9/EX9.3/ex9_3.sce @@ -0,0 +1,32 @@ +clc;
+warning("off");
+printf("\n\n example9.3 - pg389");
+Nblades=4; // no. of blades
+d=9/12; //[ft] - diameter of the impeller
+dt=30/12; //[ft] - diameter of the tank
+Nbaffles=4; // no. of baffles
+h=30; // [inch] - height of unit
+mu=10; //[cP] - viscosity of fluid
+sg=1.1; // specific gravity of fluid
+s=300; //[rpm] - speed of agitator
+CbyT=0.3;
+V=(%pi*dt^3)/4; //volume of tank in ft^3
+V1=V*7.48; //[gal] - volume of tank in gallons
+mu=mu*(6.72*10^-4); //[lb/ft*sec]
+p=sg*62.4; //[lb/ft^3] - density of fluid
+N=s/60; //[rps] - impeller speed in revolutions per second
+Nre=((d^2)*N*p)/mu;
+disp(Nre,"Nre=");
+printf("\n\n Therefore the agitator operates in the turbulent region");
+Npo=1.62;
+gc=32.174;
+P=(Npo*(p*(N^3)*(d^5)))/(gc*550);
+Cf=63025;
+Tq=(P/s)*Cf;
+PbyV=P/V;
+PbyV1=P/V1;
+TqbyV=Tq/V;
+TqbyV1=Tq/V1;
+printf("\n\n The power per unit volume and the torque per unit volume is \n P/V=%f hp/ft^3=%f hp/gal\n Tq/V=%f in*lb/ft^3=%f in*lb/gal",PbyV,PbyV1,TqbyV,TqbyV1);
+
+
diff --git a/914/CH9/EX9.4/ex9_4.sce b/914/CH9/EX9.4/ex9_4.sce new file mode 100755 index 000000000..4648801a5 --- /dev/null +++ b/914/CH9/EX9.4/ex9_4.sce @@ -0,0 +1,25 @@ +clc;
+warning("off");
+printf("\n\n example9.4 - pg391");
+// given
+Tpilot=30;
+Tlab=10;
+N1=690;
+N2=271;
+D2=3;
+D1=1;
+n=(log(N1/N2))/(log(D2/D1));
+V=12000/7.48; //[ft^3]
+T=((4*V)/%pi)^(1/3); //[ft]
+R=12.69/(30/12);
+N3=N2*(1/R)^n; //[rpm] - impeller speed in the reactor
+disp(N3,"impeller speed in rpm=");
+D3=0.75*R; //[ft] - reactor impeller diameter
+disp(D3,"reactor impeller diameter in ft=");
+P=0.1374*((N3/N2)^3)*(R^5);
+disp(P,"power in hp=");
+Cf=63025;
+Tq=(P/N3)*Cf; //[inch*lb]
+disp(Tq,"torque in inch*lb=");
+printf("\n\n At this point, the design is complete. A standard size impeller would be chosen as well as a standard size motor(7.5 hp or 10 hp)");
+
diff --git a/914/CH9/EX9.5/ex9_5.sce b/914/CH9/EX9.5/ex9_5.sce new file mode 100755 index 000000000..1d2e22124 --- /dev/null +++ b/914/CH9/EX9.5/ex9_5.sce @@ -0,0 +1,17 @@ +clc;
+warning("off");
+printf("\n\n example9.5 - pg 393");
+// given
+n=[0.5 0.6 0.7 0.8 0.9 1.0];
+D2=3.806;
+D1=0.25;
+R=D2/D1;
+N1=690;
+N2=N1*((D1/D2)^n);
+P1=9.33*10^-3; //[hp]
+P2=P1*R^(5-3*n);
+printf("\n\n n N,rpm P,hp");
+for i=1:6
+printf("\n %f %f %f",n(i),N2(i),P2(i));
+end
+
|