diff options
Diffstat (limited to '497/CH9')
| -rwxr-xr-x | 497/CH9/EX9.1/Chap9_Ex1.sce | 47 | ||||
| -rwxr-xr-x | 497/CH9/EX9.2/Chap9_Ex2.sce | 68 | ||||
| -rwxr-xr-x | 497/CH9/EX9.3/Chap9_Ex3.sce | 45 |
3 files changed, 160 insertions, 0 deletions
diff --git a/497/CH9/EX9.1/Chap9_Ex1.sce b/497/CH9/EX9.1/Chap9_Ex1.sce new file mode 100755 index 000000000..1f9328035 --- /dev/null +++ b/497/CH9/EX9.1/Chap9_Ex1.sce @@ -0,0 +1,47 @@ +//Kunii D., Levenspiel O., 1991. Fluidization Engineering(II Edition). Butterworth-Heinemann, MA, pp 491
+
+//Chapter-9, Example 1, Page 218
+//Title: Vertical Movement of Solids
+//==========================================================================================================
+
+clear
+clc
+
+//INPUT
+umf=0.015;//Velocity at minimum fluidization condition in m/s
+ephsilonmf=0.5;//Void fraction at minimum fluidization condition
+uo=0.1;//Superficial gas velocity in m/s
+delta=0.2;//Bed fraction in bubbles
+db=0.06;//Equilibrium bubble size in m
+dt=[0.1;0.3;0.6;1.5];//Various vessel sizes in m
+ub=[0.4;0.75;0.85;1.1];//Bubble velocity in m/s
+Dsv=[0.03;0.11;0.14;0.23];//Reported values of vertical dispersion coefficient
+
+//CALCULATION
+n=length(ub);
+i=1;
+fw1=2;//Wake fraction from Hamilton et al.
+fw2=0.32;//Wake fraction from Fig.(5.8)
+fw=(fw1+fw2)*0.5;//Average value of wake fraction
+while i<=n
+ Dsv1(i)=12*((uo*100)^0.5)*((dt(i)*100)^0.9);//Vertical distribution coefficient from Eqn.(3)
+ Dsv2(i)=(fw^2*ephsilonmf*delta*db*ub(i)^2)/(3*umf);//Vertical distribution coefficient from Eqn.(12)
+ i=i+1;
+end
+
+//OUTPUT
+printf('\n\t\tVertical dispersion coefficient(m^2/s)');
+printf('\nVessel Size(m)');
+printf('\tFrom Experiment');
+printf('\tFrom Eqn.(3)');
+printf('\tFrom Eqn.(12)');
+i=1;
+while i<=n
+ mprintf('\n%f',dt(i));
+ mprintf('\t%f',Dsv(i));
+ mprintf('\t%f',Dsv1(i)/10^4);
+ mprintf('\t%f',Dsv2(i));
+ i=i+1;
+end
+
+//====================================END OF PROGRAM ======================================================
\ No newline at end of file diff --git a/497/CH9/EX9.2/Chap9_Ex2.sce b/497/CH9/EX9.2/Chap9_Ex2.sce new file mode 100755 index 000000000..f8eb9da3b --- /dev/null +++ b/497/CH9/EX9.2/Chap9_Ex2.sce @@ -0,0 +1,68 @@ +//Kunii D., Levenspiel O., 1991. Fluidization Engineering(II Edition). Butterworth-Heinemann, MA, pp 491
+
+//Chapter-9, Example 2, Page 222
+//Title: Horizontal Drift Of Solids
+//==========================================================================================================
+
+clear
+clc
+
+//INPUT
+Lmf=0.83;//Length of bed at minimum fluidization condition in m
+dp=450;//Average particle size in micrometer
+ephsilonmf=0.42;//Void fraction at minimum fluidization condition
+umf=0.17;//Velocity at minimum fluidization condition in m/s
+uo=[0.37;0.47;0.57;0.67];//Superficial gas velocity in m/s
+Dsh=[0.0012;0.0018;0.0021;0.0025];//Horizontal Drift Coefficient from Experiment in m^2/s
+db=[0.10;0.14];//Equilibrium bubble size in m
+g=9.81;//Acceleration due to gravity in m/s^2
+
+
+//CALCULATION
+n=length(uo);
+m=length(db);
+j=1;
+i=1;
+k=1;
+alpha=0.77;//Since we are not dealing with Geldart A or AB solids
+uf=umf/ephsilonmf;
+for j = 1:m
+ for i = 1:n
+ ubr(k)=0.711*(db(j)*g)^0.5;//Rise velocity of a single bubble in m/s
+ ub(k)=uo(i)-umf+ubr(k);//Rise velocity of bubbles in a bubbling bed
+ delta(k)=(uo(i)-umf)/(ub(k)+umf);//Bed fraction in bubbles
+ if ubr(i)>uf then Dshc(k)=(3/16)*(delta(k)/(1-delta(k)))*((alpha^2*db(j)*ubr(k)*[(((ubr(k)+2*uf)/(ubr(k)-uf))^(1/3))-1]));//Horizontal Distribution coeff. from Eqn.(14)
+ else Dsh(k)=(3/16)*(delta/(1-delta))*(alpha^2*umf*db/ephsilonmf);//Horizontal Distribution coeff. from Eqn.(15)
+ end
+ Dshc(k)=(3/16)*(delta(k)/(1-delta(k)))*((alpha^2*db(j)*ubr(k)*[(((ubr(k)+2*uf)/(ubr(k)-uf))^(1/3))-1]));//Horizontal Distribution coeff. from Eqn.(14)
+ i=i+1;
+ k=k+1;
+ end
+ i=1;
+ j=j+1;
+end
+
+//OUTPUT
+i=1;
+j=1;
+k=1;
+while k<=m*n
+ mprintf('\nSnce we do not have ub=%fm/s>>uf=%fm/s we use Eqn.(14).',ub(k),uf)
+ printf('\nGas Velocity(m/s)');
+ printf('\tHorizontal Drift Coefficient Calculated(m^2/s)');
+ printf('\tHorizontal Drift Coefficient from Experiment(m^2/s)');
+ while j<=m
+ mprintf('\ndb=%fm',db(j));
+ while i<=n
+ mprintf('\n%f',uo(i));
+ mprintf('\t\t%f',Dshc(k));
+ mprintf('\t\t\t\t\t%f',Dsh(i));
+ i=i+1;
+ k=k+1;
+ end
+ i=1;
+ j=j+1;
+ end
+end
+
+//====================================END OF PROGRAM ======================================================
\ No newline at end of file diff --git a/497/CH9/EX9.3/Chap9_Ex3.sce b/497/CH9/EX9.3/Chap9_Ex3.sce new file mode 100755 index 000000000..0cceb0c04 --- /dev/null +++ b/497/CH9/EX9.3/Chap9_Ex3.sce @@ -0,0 +1,45 @@ +//Kunii D., Levenspiel O., 1991. Fluidization Engineering(II Edition). Butterworth-Heinemann, MA, pp 491
+
+//Chapter-9, Example 3, Page 232
+//Title: Design of Baffle Plates
+//==========================================================================================================
+
+clear
+clc
+
+//INPUT
+Gsup=1.5;//Solid interchange rate in kg/m^2plate s
+dor=19.1;//Orifice diameter in mm
+dp=210;//Particle size in micrometer
+uo=0.4;//Superficial gas velocity in m/s
+fopen=[0.12;0.17;0.26];//Open area fraction
+pi=3.14;
+
+//CALCULATION
+n=length(fopen);
+i=1;
+while i<=n
+ uor(i)=uo/fopen(i);//Gas velocity through the orifice
+ ls1(i)=Gsup/fopen(i);//Flux of solids through the holes
+ i=i+1;
+end
+ls2=[12;20;25];//Flux of solids through holes from Fig.13(c) for different uor values
+fopen1=0.12;//Open area fraction which gives reasonable fit
+lor=sqrt(((pi/4)*dor^2)/fopen1);//Orifice spacing
+
+//OUTPUT
+printf('\nfopen');
+printf('\t\tuor(m/s)');
+printf('\tls from Eqn.(18)');
+printf('\tls from Fig.13(c)');
+i=1;
+while i<=n
+ mprintf('\n%f',fopen(i));
+ mprintf('\t%f',uor(i));
+ mprintf('\t%f',ls1(i));
+ mprintf('\t\t%f',ls2(i));
+ i=i+1;
+end
+mprintf('\n\nFor square pitch, the orifice spacing should be %fmm',lor);
+
+//====================================END OF PROGRAM ======================================================
\ No newline at end of file |