initial commit / add all books

6 files changed, 192 insertions, 0 deletions

diff --git a/497/CH11/EX11.1/Chap11_Ex1.sce b/497/CH11/EX11.1/Chap11_Ex1.sce
new file mode 100755
index 000000000..56754f278
--- /dev/null
+++ b/497/CH11/EX11.1/Chap11_Ex1.sce
@@ -0,0 +1,56 @@
+//Kunii D., Levenspiel O., 1991. Fluidization Engineering(II Edition). Butterworth-Heinemann, MA, pp 491

+

+//Chapter-11, Example 1, Page 265

+//Title: Fitting Reported Mass Transfer Data with the Bubbling Bed Model

+//==========================================================================================================

+

+clear

+clc

+

+//INPUT

+db=0.37;//Equilibrium bubble size in cm

+dp=0.028;//Particle size in cm

+rhos=1.06;//Density of solids in g/cc

+ephsilonmf=0.5;//Void fraction at minimum fluidization condition

+phis=0.4;//Sphericity of solids

+gammab=0.005;//Ratio of volume of dispersed solids to that of bubble phase

+rhog=1.18E-3;//Density of air in g/cc

+myu=1.8E-4;//Viscosity of gas in g/cm s

+D=0.065;//Diffusion coefficient of gas in cm^2/s

+Sc=2.35;//Schmidt number

+etad=1;//Adsorption efficiency factor

+y=1;

+umf=1.21;//Velocity at minimum fluidization condition in cm/s

+ut=69;//Terminal velocity in cm/s

+g=980;//Acceleration due to gravity in square cm/s^2

+uo=[10;20;30;40;50];//Superficial gas velocity in cm/s

+

+//CALCULATION

+n=length(uo);

+i=1;

+Rept=(dp*ut*rhog)/myu;

+Shstar=2+(0.6*(Rept^0.5)*(Sc^(1/3)));//Sherwood no. from Eqn.(1)

+Kbc=4.5*(umf/db)+5.85*((D^0.5*g^0.25)/db^(5/4));//Gas interchange coefficient between bubble and cloud from Eqn.(10.27)

+ubr=0.711*(g*db)^0.5;//Rise velocity of the bubble

+while i<=n

+    x(i)=(uo(i)-umf)/(ubr*(1-ephsilonmf));//The term delta/(1-epshilonf) after simplification

+    Shbed(i)=x(i)*[(gammab*Shstar*etad)+((phis*dp^2*y)/(6*D))*Kbc];//Sherwood no. from Eqn.(11)

+    Rep(i)=(dp*uo(i)*rhog)/myu;//Reynolds of the particle

+    i=i+1;

+end

+

+//OUTPUT

+printf('\nThe desired result is the relationship between Shbed and Rep  The points gives a straight line of the form y=mx+c');

+printf('\nRep');

+printf('\t\tShbed');

+i=1;

+while i<=n

+    printf('\n%f',Rep(i));

+    printf('\t%f',Shbed(i));

+    i=i+1;

+end

+plot(Rep,Shbed);

+xlabel("Rep");

+ylabel("Shbed");

+

+//====================================END OF PROGRAM ======================================================
\ No newline at end of file
diff --git a/497/CH11/EX11.1/Chap11_Ex1_R.jpg b/497/CH11/EX11.1/Chap11_Ex1_R.jpg
new file mode 100755
index 000000000..afda114a3
--- /dev/null
+++ b/497/CH11/EX11.1/Chap11_Ex1_R.jpg
diff --git a/497/CH11/EX11.2/Chap11_Ex2.sce b/497/CH11/EX11.2/Chap11_Ex2.sce
new file mode 100755
index 000000000..a942045dc
--- /dev/null
+++ b/497/CH11/EX11.2/Chap11_Ex2.sce
@@ -0,0 +1,40 @@
+//Kunii D., Levenspiel O., 1991. Fluidization Engineering(II Edition). Butterworth-Heinemann, MA, pp 491

+

+//Chapter-11, Example 2, Page 267

+//Title: The Effect of m on Bubble-Emulsion Interchange

+//==========================================================================================================

+

+clear

+clc

+

+//INPUT

+umf=0.12;//Velocity at minimum fluidization condition in cm/s

+uo=40;//Superficial gas velocity in cm/s

+ub=120;//Velocity of the bubble in cm/s

+D=0.7;//Diffusion coefficient of gas in cm^2/s

+abkbe1=1;//Bubble-emuslion interchange coefficient for non absorbing particles(m=0)

+abkbe2=18;//Bubble-emuslion interchange coefficient for highly absorbing particles(m=infinity)

+g=980;//Acceleration due to gravity in square cm/s^2

+

+//CALCULATION

+//For non absorbing particles m=0,etad=0

+Kbc=(ub/uo)*(abkbe1);

+dbguess=2;//Guess value of db

+function[fn]=solver_func(db)//Function defined for solving the system

+    fn=abkbe1-(uo/ub)*(4.5*(umf/db)+5.85*(D^0.5*g^0.25)/(db^(5/4)));//Eqn.(10.27)

+endfunction

+[d]=fsolve(dbguess,solver_func,1E-6);//Using inbuilt function fsolve for solving Eqn.(10.27) for db

+//For highly absorbing particles m=infinity, etad=1

+M=abkbe2-(uo/ub)*Kbc;

+//For intermediate condition

+alpha=100;

+m=10;

+etad=1/(1+(alpha/m));//Fitted adsorption efficiency factor from Eqn.(23)

+abkbe3=M*etad+(uo/ub)*Kbc;

+

+//OUTPUT

+mprintf('\nFor non absorbing particles:\n\tDiameter of bubble=%fcm\n\tBubble-cloud interchange coefficient=%fs^-1',d,Kbc);

+mprintf('\nFor highly absorbing partilces:\n\tM=%f',M);

+mprintf('\nFor intermediate condition:\n\tFitted adsorption efficiency factor:%f\n\tBubble-emuslion interchange coefficient:%fs^-1',etad,abkbe3);

+

+//====================================END OF PROGRAM ======================================================
\ No newline at end of file
diff --git a/497/CH11/EX11.3/Chap11_Ex3.sce b/497/CH11/EX11.3/Chap11_Ex3.sce
new file mode 100755
index 000000000..d972b9fd3
--- /dev/null
+++ b/497/CH11/EX11.3/Chap11_Ex3.sce
@@ -0,0 +1,55 @@
+//Kunii D., Levenspiel O., 1991. Fluidization Engineering(II Edition). Butterworth-Heinemann, MA, pp 491

+

+//Chapter-11, Example 3, Page 273

+//Title: Fitting Reported Heat Transfer Data with the Bubbling Bed Model

+//==========================================================================================================

+

+clear

+clc

+

+//INPUT

+rhos=1.3;//Density of solids in g/cc

+phis=0.806;//Sphericity of solids

+gammab=0.001;//Ratio of volume of dispersed solids to that of bubble phase

+rhog=1.18E-3;//Density of air in g/cc

+Pr=0.69;//Prandtl number

+myu=1.8E-4;//Viscosity of gas in g/cm s

+Cpg=1.00;//Specific heat capacity of gas in J/g K

+ephsilonmf=0.45;//Void fraction at minimum fluidization condition

+kg=2.61E-4;//Thermal concuctivity of gas in W/cm k

+dp=0.036;//Particle size in cm

+umf=6.5;//Velocity at minimum fluidization condition in cm/s

+ut=150;//Terminal velocity in cm/s

+db=0.4;//Equilibrium bubble size in cm

+etah=1;//Efficiency of heat transfer

+uo=[10;20;30;40;50];//Superficial gas velocity in cm/s

+g=980;//Acceleration due to gravity in square cm/s^2

+

+//CALCULATION

+Nustar=2+[((dp*ut*rhog)/myu)^0.5*Pr^(1/3)];//Nusselt no. from Eqn.(25)

+Hbc=4.5*(umf*rhog*Cpg/db)+5.85*((kg*rhog*Cpg)^0.5*g^0.25/db^(5/4));//Total heat interchange across the bubble-cloud boundary from Eqn.(32)

+ubr=0.711*(g*db)^0.5;//Rise velocity of the bubble from Eqn.(6.7)

+n=length(uo);

+i=1;

+while i<=n

+    x(i)=(uo(i)-umf)/(ubr*(1-ephsilonmf));//The term delta/(1-epshilonf) after simplification

+    Nubed(i)=x(i)*[gammab*Nustar*etah+(phis*dp^2/(6*kg))*Hbc];//Nusselt no. from Eqn.(36)

+    Rep(i)=(dp*uo(i)*rhog)/myu;//Reynolds of the particle

+    i=i+1;

+end

+

+//OUTPUT

+printf('\nThe desired result is the relationship between Nubed and Rep which is in the form of a straight line y=mx+c');

+printf('\nRep');

+printf('\t\tNubed');

+i=1;

+while i<=n

+    printf('\n%f',Rep(i));

+    printf('\t%f',Nubed(i));

+    i=i+1;

+end

+plot(Rep,Nubed);

+xlabel("Rep");

+ylabel("Nubed");

+ 

+//====================================END OF PROGRAM ======================================================
\ No newline at end of file
diff --git a/497/CH11/EX11.3/Chap11_Ex3_R.jpg b/497/CH11/EX11.3/Chap11_Ex3_R.jpg
new file mode 100755
index 000000000..719ebfecb
--- /dev/null
+++ b/497/CH11/EX11.3/Chap11_Ex3_R.jpg
diff --git a/497/CH11/EX11.4/Chap11_Ex4.sce b/497/CH11/EX11.4/Chap11_Ex4.sce
new file mode 100755
index 000000000..b3f2ff179
--- /dev/null
+++ b/497/CH11/EX11.4/Chap11_Ex4.sce
@@ -0,0 +1,41 @@
+//Kunii D., Levenspiel O., 1991. Fluidization Engineering(II Edition). Butterworth-Heinemann, MA, pp 491

+

+//Chapter-11, Example 4, Page 274

+//Title: Heating a Particle in a Fluidized Bed

+//==========================================================================================================

+

+clear

+clc

+

+//INPUT

+rhog=1.2;//Density of air in kg/m^3

+myu=1.8E-5;//Viscosity of gas in kg/m s

+kg=2.6E-2;//Thermal concuctivity of gas in W/m k

+dp=1E-4;//Particle size in m

+rhos=8920;//Density of solids in kg/m^3

+Cps=390;//Specific heat capacity of the solid in J/kg K

+ephsilonf=0.5;//Void fraction of the fluidized bed

+umf=0.1;//Velocity at minimum fluidization condition in m/s

+uo=0.1;//Superficial gas velocity in m/s

+pi=3.14

+

+//CALCULATION

+to=0;//Initial temperature of the bed

+T=100;//Temperature of the bed

+t=0.99*T;//Particle temperature i.e. when it approaches 1% of the bed temperature

+mp=(pi/6)*dp^3*rhos;//Mass of the particle

+A=pi*dp^2;//Surface area of the particle

+Rep=(dp*uo*rhog)/myu;//Reynold's no. of the particle

+Nubed=0.0178;//Nusselt no. from Fig.(6)

+hbed1=(Nubed*kg)/dp;//Heat transfer coefficient of the bed

+t1=(mp*Cps/(hbed1*A))*log((T-to)/(T-t));//Time needed for the particle approach 1 percentage of the bed temperature in case(a)

+hbed2=140*hbed1;//Since from Fig.(6) Nup is 140 times Nubed

+t2=(mp*Cps/(hbed2*A))*log((T-to)/(T-t));//Time needed for the particle approach 1 percentage of the bed temperature in case(b)

+

+//OUTPUT

+printf('\nCase(a):Using the whole bed coefficient from Fig.(6)');

+mprintf('\n\tTime needed for the particle approach 1 percentage of the bed temperature is %fs',t1);

+printf('\nCase(b):Uisng the single-particle coefficient of Eqn.(25),also shown in Fig.(6)');

+mprintf('\n\tTime needed for the particle approach 1 percentage of the bed temperature is %fs',t2);

+

+//====================================END OF PROGRAM ======================================================
\ No newline at end of file