1 files changed, 72 insertions, 0 deletions
diff --git a/1040/CH6/EX6.1/Ex6_1.sce b/1040/CH6/EX6.1/Ex6_1.sce
new file mode 100644
index 000000000..9af0bfd8e
--- /dev/null
+++ b/1040/CH6/EX6.1/Ex6_1.sce
@@ -0,0 +1,72 @@
+//Harriot P., 2003, Chemical Reactor Design (I-Edition), Marcel Dekker, Inc., USA, pp 436.
+//Chapter-6 Ex6.1 Pg No.236
+//Title:Power Consumption at 300 rpm,speed of stirrer and blending time
+//====================================================================================================================
+clear
+clc
+// COMMON INPUT
+D_a=0.1;
+D_t=0.3;
+H=0.3;
+N_P=5.5;
+rho=1000;
+n=5;
+S_f=6;//Scale up factor in diameter
+P_by_V_limit=10;//Pressure per unit volume (HP/1000gal)
+n1=5;
+Da_by_Dt1=D_a/D_t;
+Da_by_Dt2=0.5;
+
+//CALCULATION (Ex6.1.a)
+P_unit_vol=(N_P*n^3*D_a^5)/(%pi*(1/4)*D_t^2*H);
+P_thousand_gal=P_unit_vol*5.067;
+t=(4/n)*(D_t/D_a)^2*(H/D_t);
+P_unit_vol_new=S_f^2*P_thousand_gal;
+
+//CALCULATION (Ex6.1.b)
+n_limit=(P_by_V_limit/P_unit_vol_new)^(1/3) *n1;//Pressure per unit vol propotional to n3
+t_inc_factor=n1/n_limit;//t inversely propotional to n
+rotational_speed=n_limit*60;//Speed in rpm
+
+//CALCULATION (Ex6.1.c)
+n2=(Da_by_Dt1/Da_by_Dt2)^(5/3)*n_limit;
+rotaional_speed=n2*60;
+t1=4*(1/Da_by_Dt1)^2*(H/D_t)*(1/n_limit);
+t2=4*(1/Da_by_Dt2)^2*(H/D_t)*(1/n2);
+
+//OUTPUT (Ex6.1.a)
+mprintf('\n OUTPUT Ex6.1.a');
+mprintf('\n==========================================================');
+mprintf('\n The Power consumption per unit volume at 300rpm = %.2f HP/1000 gal',P_thousand_gal);
+mprintf('\n\ The Power consumption scaling up sixfold in diameter = %.0f HP/1000 gal',P_unit_vol_new); 
+
+
+//OUTPUT (Ex6.1.b)
+mprintf('\n\n\n OUTPUT Ex6.1.b');
+mprintf('\n==========================================================');
+mprintf('\n The speed of the stirrer  = %.2f sec-1 or %.0f rpm',n_limit,rotational_speed);
+mprintf('\n Blending time increases by factor of %.2f ',t_inc_factor); 
+
+//OUTPUT(Ex6.1.c)
+mprintf('\n\n\n OUTPUT Ex6.1.c');
+mprintf('\n==========================================================');
+mprintf('\n The new stirrer speed = %.2f sec-1 or %.0f rpm',n2,rotaional_speed); 
+mprintf('\n The new blending time  for Da/Dt ratio of 0.5 = %.1f sec',t2); 
+
+//FILE OUTPUT
+fid= mopen('.\Chapter6-Ex1-Output.txt','w');
+mfprintf(fid,'\n OUTPUT Ex6.1.a');
+mfprintf(fid,'\n==========================================================');
+mfprintf(fid,'\n The Power consumption per unit volume at 300rpm = %.2f HP/1000 gal',P_thousand_gal);
+mfprintf(fid,'\n\ The Power consumption scaling up sixfold in diameter = %.0f HP/1000 gal',P_unit_vol_new);
+mfprintf(fid,'\n\n\n OUTPUT Ex6.1.b');
+mfprintf(fid,'\n==========================================================');
+mfprintf(fid,'\n The speed of the stirrer  = %.2f sec-1 or %.0f rpm',n_limit,rotational_speed);
+mfprintf(fid,'\n Blending time increases by factor of %.2f ',t_inc_factor); 
+mfprintf(fid,'\n\n\n OUTPUT Ex6.1.c');
+mfprintf(fid,'\n==========================================================');
+mfprintf(fid,'\n The new stirrer speed = %.2f sec-1 or %.0f rpm',n2,rotaional_speed); 
+mfprintf(fid,'\n The new blending time  for Da/Dt ratio of 0.5 = %.1f sec',t2);
+mclose(fid);
+//======================================================END OF PROGRAM=================================================
+//Disclaimer: In Ex6.1.c there is an arithematic error in the value of D_a/D_t. The value of  D_a/D_t should be 11.4 instead of the value reported in the textbook for D_a/D_t=11.1.