From b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b Mon Sep 17 00:00:00 2001 From: priyanka Date: Wed, 24 Jun 2015 15:03:17 +0530 Subject: initial commit / add all books --- 497/CH14/EX14.1/Chap14_Ex1.sce | 36 ++++++++++++++++++++ 497/CH14/EX14.2/Chap14_Ex2.sce | 72 +++++++++++++++++++++++++++++++++++++++ 497/CH14/EX14.2/Chap14_Ex2_R.jpg | Bin 0 -> 23087 bytes 497/CH14/EX14.3/Chap14_Ex3.sce | 22 ++++++++++++ 497/CH14/EX14.4/Chap14_Ex4.sce | 30 ++++++++++++++++ 497/CH14/EX14.5/Chap14_Ex5.sce | 52 ++++++++++++++++++++++++++++ 6 files changed, 212 insertions(+) create mode 100755 497/CH14/EX14.1/Chap14_Ex1.sce create mode 100755 497/CH14/EX14.2/Chap14_Ex2.sce create mode 100755 497/CH14/EX14.2/Chap14_Ex2_R.jpg create mode 100755 497/CH14/EX14.3/Chap14_Ex3.sce create mode 100755 497/CH14/EX14.4/Chap14_Ex4.sce create mode 100755 497/CH14/EX14.5/Chap14_Ex5.sce (limited to '497/CH14') diff --git a/497/CH14/EX14.1/Chap14_Ex1.sce b/497/CH14/EX14.1/Chap14_Ex1.sce new file mode 100755 index 000000000..8417234c0 --- /dev/null +++ b/497/CH14/EX14.1/Chap14_Ex1.sce @@ -0,0 +1,36 @@ +//Kunii D., Levenspiel O., 1991. Fluidization Engineering(II Edition). Butterworth-Heinemann, MA, pp 491 + +//Chapter-14, Example 1, Page 343 +//Title: Flow with Elutriation +//========================================================================================================== + +clear +clc + +//INPUT +Fo=2.7;//Feed rate in kg/min +Fof=0.9;//Feed rate of fines in feed in kg/min +Foc=1.8;//Feed rate of coarse in feed in kg/min +W=17;//Bed weight in kg +kf=0.8;//Elutriation of fines in min^-1 +kc=0.0125;//Elutriation of coarse in min^-1 + +//CALCULATION +F1guess=1;//Guess value of F1 +function[fn]=solver_func(F1)//Function defined for solving the system + fn=F1-(Fof/(1+(W/F1)*kf))-(Foc/(1+(W/F1)*kc));//Eqn.(17) +endfunction +[F1]=fsolve(F1guess,solver_func,1E-6);//Inbuilt function fsolve to solve for F1 +F1f=Fof/(1+(W/F1)*kf);//Flow rate of fines in entrained streams from Eqn.(16) +F1c=Foc/(1+(W/F1)*kc);//Flow rate of coarse in entrained streams from Eqn.(16) +F2f=Fof-F1f;//Flow rate of fines in overflow streams from Eqn.(9) +F2c=Foc-F1c;//Flow rate of coarse in overflow streams from Eqn.(9) +tbarf=1/((F1/W)+kf);//Mean residence time of fines from Eqn.(12) +tbarc=1/((F1/W)+kc);//Mean residence time of coarse from Eqn.(12) + +//OUTPUT +mprintf('\nFlow rate in entrained stream:\n\tFines:%fkg/min\n\tCoarse:%fkg/min',F1f,F1c); +mprintf('\nFlow rate in overflow stream:\n\tFines:%fkg/min\n\tCoarse:%fkg/min',F2f,F2c); +mprintf('\nMean residence time:\n\tFines:%fmins\n\tCoarse:%fmins',tbarf,tbarc); + +//====================================END OF PROGRAM ====================================================== \ No newline at end of file diff --git a/497/CH14/EX14.2/Chap14_Ex2.sce b/497/CH14/EX14.2/Chap14_Ex2.sce new file mode 100755 index 000000000..a1cf0488c --- /dev/null +++ b/497/CH14/EX14.2/Chap14_Ex2.sce @@ -0,0 +1,72 @@ +//Kunii D., Levenspiel O., 1991. Fluidization Engineering(II Edition). Butterworth-Heinemann, MA, pp 491 + +//Chapter-14, Example 2, Page 344 +//Title: Flow with Elutriation and Change in Density of Solids +//========================================================================================================== + +clear +clc + +//INPUT +dt=4;//Diameter of reactor in m +ephsilonm=0.4;//Void fraction of static bed +rhos=2500;//Density of solid in the bed in kg/m^3 +Lm=1.2;//Height of static bed in m +Fo=3000;//Feed rate in kg/hr +beta1=1.2;//Increase in density of solids +dp=[3;4;5;6;7;8;9;10;11;12;3;14;16;18;20;22;24;26;28;30]*10^-2;//Size of particles in mm +po=[0;0.3;0.8;1.3;1.9;2.6;3.5;4.4;5.7;6.7;7.5;7.8;7.5;6.3;5.0;3.6;2.4;1.3;0.5;0];//Size distribution of solids in mm^-1 +k=[0;10;9.75;9.5;8.75;7.5;6.0;4.38;2.62;1.20;0.325;0;0;0;0;0;0;0;0;0]*10^-4;//Elutriation constant in s^-1 +pi=3.14; + +//CALCULATION +W=(pi/4*dt^2)*Lm*(1-ephsilonm)*rhos;//Weight of solids in bed +n=length(dp); +i=1; +F1guess=1000;//Guess value for F1 +F1c=2510:10:2700; +while i<=n + function[fn]=solver_func(F1)//Function defined for solving the system + if k(i)==0 then x(i)=0; break + else x(i)=(po(i)/(W*k(i)/F1))*log(1+(W*k(i)/F1)); + end + fn=F1/(Lm*Fo)-x(i); + endfunction + [F1(i)]=fsolve(F1guess,solver_func,1E-6);//Using inbuilt function fsolve for solving Eqn.(20) for F1 + c(i)=F1c(i)/(Lm*Fo); + if F1(i)==0 then a(i)=0; + else a(i)=(po(i)/(W*k(i)/F1(i)))*log(1+(W*k(i)/F1(i))); + end + i=i+1; +end +plot(F1,a,F1,c); +xtitle('F1 vs a,c','F1','a,c'); +F1n=2500;//The point were both the curves meet +F2=beta1*Fo-F1n;//Flow rate of the second leaving stream +j=1; +m=length(dp); +while j<=m + p1(j)=(1/F1n)*((Fo*po(j))/(1+(W/F1n)*k(j)));//Size distribution of stream 1 in mm^-1 from Eqn.(16) + p2(j)=k(j)*W*p1(j)/F2;//Size distribution of stream 2 in mm^-1 from Eqn.(7) + if p1(j)==0 & p2(j)==0 then tbar(j)=0; + else if p1(j)==0 then tbar(j)=(W*p1(j))/(F2*p2(j)); + else if p2(j)==0 then tbar(j)=(W*p1(j))/(F1n*p1(j)); + else tbar(j)=(W*p1(j))/(F1n*p1(j)+F2*p2(j));//Average time in hr from Eqn.(11) + end + end + end + j=j+1; +end + +//OUTPUT +printf('\nFlow rate of stream 1:%fkg/hr',F1n); +printf('\nFlow rate of stream 2:%fkg/hr',F2); +j=1; +mprintf('\ntbar(hr)'); +while j<=m + mprintf('\n%f',tbar(j)); + j=j+1; +end + +//====================================END OF PROGRAM ====================================================== +//DISCLAIMER: The value obtained for tbar is deviating highly form the one given in textbook. However, the value obtained by manual calculation is close to the ones obtained from the program. \ No newline at end of file diff --git a/497/CH14/EX14.2/Chap14_Ex2_R.jpg b/497/CH14/EX14.2/Chap14_Ex2_R.jpg new file mode 100755 index 000000000..e38a4223f Binary files /dev/null and b/497/CH14/EX14.2/Chap14_Ex2_R.jpg differ diff --git a/497/CH14/EX14.3/Chap14_Ex3.sce b/497/CH14/EX14.3/Chap14_Ex3.sce new file mode 100755 index 000000000..4e1a107a1 --- /dev/null +++ b/497/CH14/EX14.3/Chap14_Ex3.sce @@ -0,0 +1,22 @@ +//Kunii D., Levenspiel O., 1991. Fluidization Engineering(II Edition). Butterworth-Heinemann, MA, pp 491 + +//Chapter-14, Example 3, Page 351 +//Title: Single-Size Feed of Shrinking Particles +//========================================================================================================== + +clear +clc + +//INPUT +dp=1;//Particle size in mm +Fo=10;//Feed rate in kg/min +k=0.1;//Particle shrinkage rate in mm/min + +//CALCULATION +R=k/2;//Particle shrinkage rate in terms of radius +W=(Fo*dp/2)/(4*R);//Bed weight from Eqn.(42) + +//OUTPUT +printf('\nWeight of bed:%fkg',W); + +//====================================END OF PROGRAM ====================================================== \ No newline at end of file diff --git a/497/CH14/EX14.4/Chap14_Ex4.sce b/497/CH14/EX14.4/Chap14_Ex4.sce new file mode 100755 index 000000000..ccc99bb2f --- /dev/null +++ b/497/CH14/EX14.4/Chap14_Ex4.sce @@ -0,0 +1,30 @@ +//Kunii D., Levenspiel O., 1991. Fluidization Engineering(II Edition). Butterworth-Heinemann, MA, pp 491 + +//Chapter-14, Example 4, Page 352 +//Title: Wide Size Distribution of Shrinking Particle +//========================================================================================================== + +clear +clc + +//INPUT +dpi=[1.05;0.95;0.85;0.75;0.65;0.55;0.45;0.35;0.25;0.15;0.05];//Mean size in mm +Fo=[0;0.5;3.5;8.8;13.5;17.0;18.2;17.0;13.5;7.3;0]*10^-2;//Feed rate in kg/s +k=[0;0;0;0;0;0;0;0;2.0;12.5;62.5]*10^-5;//Elutriation constant in s^-1 +R=-1.58*10^-5;//Rate of particle shrinkage in mm/s +deldpi=0.1;//Size intervals in mm + +//CALCULATION +n=length(dpi); +m=2;//Starting with the largest value size interval that contains solids +W(m-1)=0; +while m<=n + W(m)=(Fo(m)-R*W(m-1)/deldpi)/(k(m)-R/deldpi-3*R/dpi(m));//From Eqn.(33) + m=m+1; +end +Wt=sum(W);//Total sum + +//OUTPUT +printf('\nTotal mass in the bed:%fkg',Wt); + +//====================================END OF PROGRAM ====================================================== \ No newline at end of file diff --git a/497/CH14/EX14.5/Chap14_Ex5.sce b/497/CH14/EX14.5/Chap14_Ex5.sce new file mode 100755 index 000000000..74a98ae99 --- /dev/null +++ b/497/CH14/EX14.5/Chap14_Ex5.sce @@ -0,0 +1,52 @@ +//Kunii D., Levenspiel O., 1991. Fluidization Engineering(II Edition). Butterworth-Heinemann, MA, pp 491 + +//Chapter-14, Example 5, Page 353 +//Title: Elutriation and Attrition of Catalyst +//========================================================================================================== + +clear +clc + +//INPUT +dpi=[0.17;0.15;0.13;0.11;0.09;0.07;0.05;0.03;0.01];//Mean size of particles in mm +a=[0;0.95;2.45;5.2;10.1;23.2;35.65;20.0;2.45]*10^-2;//Feed composition Fo(dpi)/Fo +y=[0;0;0;0;0;0;0.625;10.225;159.25]*10^-6;//Elutriation and cyclone efficiency k(dpi)(1-eta(dpi)) +F=0.01;//Rate at which solids are withdrawn in kg/s +W=40000;//Weight of bed in kg +dp1=0.11//Initial size in mm +dp2=0.085;//Size after shrinking in mm +dpmin=0.01;//Minimum size in mm +deldpi=2*10^-2;//Size inerval in mm +t=20.8;//Time in days +si=1; + +//CALCULATION +kdash=log((dp1-dpmin)/(dp2-dpmin))/(t*24*3600);//Rate of particle shrinkage from Eqn.(24) +n=length(dpi); +m=2; +Fo=0.05;//Initial value of Fo +F1(m-1)=0; +s=0; +c=0; +t=1E-6; +while m<=n + R(m)=-kdash*(dpi(m)-dpmin);//Rate of size change + x(m)=(a(m)*Fo-W*R(m-1)*F1(m-1)/deldpi)/(F+(W*y(m))-(W*R(m)/deldpi)-3*W*R(m)/dpi(m));//Eqn.(34) + F1(m)=x(m)*F; + c=c+x(m); + m=m+1; + if abs(c-1)