4 files changed, 152 insertions, 0 deletions
diff --git a/497/CH10/EX10.1/Chap10_Ex1.sce b/497/CH10/EX10.1/Chap10_Ex1.sce
new file mode 100755
index 000000000..95be23a0e
--- /dev/null
+++ b/497/CH10/EX10.1/Chap10_Ex1.sce
@@ -0,0 +1,60 @@
+//Kunii D., Levenspiel O., 1991. Fluidization Engineering(II Edition). Butterworth-Heinemann, MA, pp 491
+
+//Chapter-10, Example 1, Page 253
+//Title: Estimate Interchange Coefficients in Bubbling Beds
+//==========================================================================================================
+
+clear
+clc
+
+//INPUT
+umf=[0.01;0.045];//Velocity at minimum fluidization condition in m/s
+ephsilonmf=[0.5;0.5];//Void fraction at minimum fluidization condition
+D=[2E-5;7E-5];//Diffusion coefficient of gas in m^2/s
+g=9.81;//Acceleration due to gravity in m/s^2
+
+//CALCULATION
+db=[5;10;15;20];
+n=length(umf);
+m=length(db)'
+for i = 1:n
+    for j = 1:m
+            Kbc(i,j)=4.5*(umf(i)/db(j))+5.85*((D(i)^0.5*g^0.25)/db(j)^(5/4));//Gas interchange coefficient between bubble and cloud from Eqn.(27)
+            Kce(i,j)=6.77*((D(i)*ephsilonmf(i)*0.711*(g*db(j))^0.5)/db(j)^3)^0.5;//Gas interchange coefficient between emulsion and cloud from Eqn.(34)
+            Kbe(i,j)=(Kbc(i,j)*Kce(i,j))/(Kbc(i,j)+Kce(i,j));//Gas interchange coefficient between bubble and emulsion from Eqn.(14)
+    end;
+end
+
+//OUTPUT
+i=1;
+j=1;
+k=1;
+while k<=m*n
+    printf('\n\t\tKbc for fine particles and He');
+    printf('\tKbc for coarse particles and ozone');
+    printf('\tKbe for fine particles and He');
+    printf('\tKbe for coarse particles and ozone');
+    while j<=m
+        mprintf('\ndb=%fm',db(j)*10^-2);
+        while i<=n
+            mprintf('\t%f',Kbc(k));
+            mprintf('\t\t\t%f',Kbe(k));
+            i=i+1;   
+            k=k+1;
+            printf('\t\t\t');
+        end
+    i=1;
+    j=j+1;
+    end
+end
+Kbe=Kbe';
+Kbc=Kbc';
+plot2d("ll",db,[Kbc Kbe]);
+xtitle('Plot of Kbc,Kbe vs db','db',['Kbc','Kbe']);
+printf('\nComparing the points with the plot of Kbc,Kbe vs db in Fig.(12), we can conlcude the following:');
+printf('\nKbc for fine particles and helium: line 2 in Fig.(12)');
+printf('\nKbc for coarser particles and ozone: line 3 in Fig.(12)');
+printf('\nKbe for fine particles and helium: line 4 in Fig.(12)');
+printf('\nKbe for coarser particles and ozone: line 5 in Fig.(12)');
+
+//====================================END OF PROGRAM ======================================================
+\ No newline at end of file
diff --git a/497/CH10/EX10.1/Chap10_Ex1_R.jpg b/497/CH10/EX10.1/Chap10_Ex1_R.jpg
new file mode 100755
index 000000000..e71881544
--- /dev/null
+++ b/497/CH10/EX10.1/Chap10_Ex1_R.jpg
diff --git a/497/CH10/EX10.2/Chap10_Ex2.sce b/497/CH10/EX10.2/Chap10_Ex2.sce
new file mode 100755
index 000000000..bb6414d36
--- /dev/null
+++ b/497/CH10/EX10.2/Chap10_Ex2.sce
@@ -0,0 +1,44 @@
+//Kunii D., Levenspiel O., 1991. Fluidization Engineering(II Edition). Butterworth-Heinemann, MA, pp 491
+
+//Chapter-10, Example 2, Page 254
+//Title: Compare the Relative Importance of Kbc and Kce
+//==========================================================================================================
+
+clear
+clc
+
+//INPUT
+D=0.69;//Diffusion coefficient of gas in cm^2/s
+umf=1.0;//Velocity at minimum fluidization condition in cm/s
+ephsilonmf=0.5;//Void fraction at minimum fluidization condition
+db=[5;15];//Equilibrium bubble size in cm
+g=980;//Acceleration due to gravity in cm/s^2
+
+//CALCULATION
+n=length(db);
+i=1;
+while i<=n
+    Kbc(i)=4.5*(umf/db(i))+5.85*((D^0.5*g^0.25)/db(i)^(5/4));//Gas interchange coefficient between bubble and cloud from Eqn.(27)
+    Kce(i)=6.77*((D*ephsilonmf*0.711*(g*db(i))^0.5)/db(i)^3)^0.5;//Gas interchange coefficient between emulsion and cloud from Eqn.(34)
+    Kbe(i)=(Kbc(i)*Kce(i))/(Kbc(i)+Kce(i));//Gas interchange coefficient between bubble and emulsion from Eqn.(14)
+    e(i)=(Kce(i)-Kbe(i))/Kbe(i);//Error when minor resistance is ignored
+    i=i+1;
+end
+
+//OUTPUT
+printf('\ndb(cm)');
+printf('\t\tCalculated Kbc');
+printf('\tCalculated Kce');
+printf('\t\tKbe from Eqn.(14)');
+printf('\tErron when minor resistance is ignored (in percentage)');
+i=1;
+while i<=n
+    mprintf('\n%f',db(i));
+    mprintf('\t%f',Kbc(i));
+    mprintf('\t%f',Kce(i));
+    mprintf('\t\t%f',Kbe(i));
+    mprintf('\t\t%f',e(i)*100);
+    i=i+1;   
+end
+
+//====================================END OF PROGRAM ======================================================
+\ No newline at end of file
diff --git a/497/CH10/EX10.3/Chap10_Ex3.sce b/497/CH10/EX10.3/Chap10_Ex3.sce
new file mode 100755
index 000000000..3b5572357
--- /dev/null
+++ b/497/CH10/EX10.3/Chap10_Ex3.sce
@@ -0,0 +1,48 @@
+//Kunii D., Levenspiel O., 1991. Fluidization Engineering(II Edition). Butterworth-Heinemann, MA, pp 491
+
+//Chapter-10, Example 3, Page 255
+//Title: Compare Interchange Rates for Adsorbed and Nonadsorbed Gases
+//==========================================================================================================
+
+clear
+clc
+
+//INPUT
+Kbe=[0.028;0.05];//Reported range for gas interchange coefficient between bubble and emulsion
+uo=0.30;//Superficial gas velocity in m/s
+db=0.13;//Equilibrium bubble size in m
+m=7;
+ephsilonmf=0.5;//Void fraction at minimum fluidization condition
+umf=0.0018;//Velocity at minimum fluidization condition in m/s
+D=[9E-6;22E-6];//Diffusion coefficient of gas in m^2/s
+g=9.81;//Acceleration due to gravity in m/s^2
+
+//CALCULATION
+n=length(Kbe);
+i=1;
+while i<=n
+    Kbem(i)=(6/db)*Kbe(i);//Gas interchange coefficient between bubble and emulsion from Eqn.(19)
+    Kbc(i)=4.5*(umf/db)+5.85*((D(i)^0.5*g^0.25)/db^(5/4));//Gas interchange coefficient between bubble and cloud from Eqn.(27)
+    Kce(i)=6.77*((D(i)*ephsilonmf*0.711*(g*db)^0.5)/db^3)^0.5;//Gas interchange coefficient between emulsion and cloud from Eqn.(34)
+    Kbe(i)=(Kbc(i)*Kce(i))/(Kbc(i)+Kce(i));//Gas interchange coefficient between bubble and emulsion from Eqn.(14)
+    c(i)=(Kbem(i)/Kbe(i));
+    i=i+1;
+end
+
+//OUTPUT
+printf('\nKbe from Eqn.(19)');
+printf('\tKbc from Eqn.(27)');
+printf('\tKce from Eqn.(34)');
+printf('\tKbe from Eqn.(14)');
+printf('\tComparison of Kbe from Eqn.(19) and that from Eqn.(14)');
+i=1;
+while i<=n
+    mprintf('\n%f',Kbem(i));
+    mprintf('\t\t%f',Kbc(i));
+    mprintf('\t\t%f',Kce(i));
+    mprintf('\t\t%f',Kbe(i));
+    mprintf('\t\t%f',c(i));
+    i=i+1;   
+end
+
+//====================================END OF PROGRAM ======================================================
+\ No newline at end of file