initial commit / add all books

2 files changed, 134 insertions, 0 deletions

diff --git a/1040/CH7/EX7.5/Chapter7_Ex5_Output.txt b/1040/CH7/EX7.5/Chapter7_Ex5_Output.txt
new file mode 100644
index 000000000..b342fabe5
--- /dev/null
+++ b/1040/CH7/EX7.5/Chapter7_Ex5_Output.txt
@@ -0,0 +1,20 @@
+

+ OUTPUT Ex7.5.a

+==========================================================

+	The maximum rate of CO absorption at 15 atm : 0.030172 (mol/L s)

+	The kinetic rate of CO absorption at 180(°C) : 0.003003 (mol/L s)

+	The predicted value of k_L_a : 0.29 (s-1)

+

+

+ OUTPUT Ex7.5.b

+==========================================================

+	The Dimensions of the reactor are 

+	Diameter:2 m

+	Height:14.34 m

+

+

+ OUTPUT Ex7.5.c

+==========================================================

+	The new dimensions of the reactor

+	Diameter:2.8 m

+	Height:6.3 m
\ No newline at end of file
diff --git a/1040/CH7/EX7.5/Ex7_5.sce b/1040/CH7/EX7.5/Ex7_5.sce
new file mode 100644
index 000000000..b987b4f34
--- /dev/null
+++ b/1040/CH7/EX7.5/Ex7_5.sce
@@ -0,0 +1,114 @@
+//Harriot P.,2003,Chemical Reactor Design (I-Edition) Marcel Dekker,Inc.,USA,pp 436

+//Chapter-7 Ex7.5 Pg No.293

+//Title:Maximum rate of CO absorption  and Dimensions of Bubble Column Reactor

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

+clear

+clc

+// COMMON INPUT

+P_dash=5;//Partial pressure of acetic acid (atm)

+P_total=20;//Total Pressure (atm)

+myu=0.19;// Viscosity of acetic acid

+T_C=180;//Temperature in (°C)

+T_K=T_C+273;//Temperature in (K)

+sigma_20=28;//Surface Tension(Dynes/cm) at 20 (°C)

+sigma_180=20;//Surface Tension (Dynes/cm)at 180 (°C)

+M_CO=28;//Molecular weight of CO

+M_B=60.05;//Molecular weight acetic acid

+V_A= 30.7;//Molar volume

+S_CO=7*10^(-3);//Solubility of CO (mol/L atm)

+f_CO=0.75;//Fraction of CO in feed

+f_acetic_acid=1-f_CO;//Fraction of Acetic acid

+R=82.056*(10^-3);//(cm3 atm/ K  mol)

+rho_air=1.21;//(kg/m3)density of air at 20 (°C)

+sigma_H2O=72;//Surface tension (Dynes/cm)

+myu_H2O=1;//Viscosity of water

+k_L_a_air_water=0.051;//(sec-1)

+D_O2_water=2.4*(10^-5);//(cm2/sec)diffusivity for oxygen in waterat 20(°C) 

+Conc_Rh=4*10^(-3);//Concentration of Rohdium(M)

+Conc_CH3I=1;//Concentration of Methyl Iodide(M)

+F_product_acetic_acid=0.1;// Rate of acetic acid produced (kmol/sec)

+f_CO_reacted=0.8;//80% of CO reacted

+u_g=0.1;//(m/sec)

+Epsilon_air_water_new=0.07;//At velocity  3(cm/sec)

+Epsilon_air_water_old= 0.12;//At velocity 6(cm/sec)

+u_g_c=5*(10^(-2));//Gas Velocity  Ex7.5.c(m/sec)

+

+

+

+//CALCUATION (Ex7.5.a)

+D_CO=(7.4*10^(-8)*M_B^(1/2)*T_K)/(myu*V_A^(0.6));//Diffusivity of CO (Wilke–Chang equation Eq4.17)

+M_ave=f_CO*M_CO+M_B*f_acetic_acid;//Average Molecular weight

+rho_g=M_ave*P_total/(R*T_K);//From ideal gas law

+epsilon_air_water= 0.12;//At velocity 6(cm/sec)

+epsilon=epsilon_air_water*(sigma_H2O/sigma_180)^(0.4)*(myu/myu_H2O)^(0.2)*(rho_g/rho_air)^(0.2);//From equation 7.64

+u_G=6;//From figure 7.12(cm/sec)

+k_L_a=k_L_a_air_water*(D_CO/D_O2_water)^(0.5)*(epsilon/epsilon_air_water);//From equation 7.69

+P_CO=P_total-P_dash;

+C_CO_Star=S_CO*P_CO;

+r_max=C_CO_Star*k_L_a;//Rate of CO absorption at 15 atm

+r_test=158.8*(10^(6))*exp(-8684/T_K)*(Conc_Rh)*(Conc_CH3I);//Kinetic rate at 180 (°C) 

+

+//CALCULATION(Ex7.5.b)

+F_feed_CO=F_product_acetic_acid/f_CO_reacted;//Rate of flow of CO (kmol/sec)

+F_total=F_feed_CO/f_CO;

+Q=F_total*R*T_K/(P_total);

+S=Q/u_g;

+D_t=sqrt(4*S/%pi);

+r_test_b=(158.8*(10^(6))*exp(-8684/T_K)*(Conc_Rh)*(Conc_CH3I))*(10^(-3));//Kinetic rate at 180 (°C) 

+liquid_vol= (F_product_acetic_acid/r_test_b)*(10^(-3));//liquid volume (m3)

+h0=liquid_vol/S;//clear liquid

+h=h0/(1-epsilon);//aerated liquid

+

+//CALCULATION(Ex7.5.c)

+Q=F_total*R*T_K/(P_total);

+S=Q/u_g_c;

+D_t_c=sqrt(4*S/%pi);

+Epsilon_new=(Epsilon_air_water_new/Epsilon_air_water_old)*epsilon;

+liquid_vol= (F_product_acetic_acid/r_test_b)*(10^(-3));//liquid volume (m3)

+h0=liquid_vol/S;//clear liquid

+h_new=h0/(1-Epsilon_new);//aerated liquid

+

+//OUTPUT (Ex7.5.a)

+mprintf('\n OUTPUT Ex7.5.a');

+mprintf('\n==========================================================');

+mprintf('\n\tThe maximum rate of CO absorption at 15 atm : %f (mol/L s)',r_max);

+mprintf('\n\tThe kinetic rate of CO absorption at 180(°C) : %f (mol/L s)',r_test);

+mprintf('\n\tThe predicted value of k_L_a : %0.2f (s-1)',k_L_a);

+

+//OUTPUT (Ex7.5.b)

+mprintf('\n\n\n OUTPUT Ex7.5.b');

+mprintf('\n==========================================================');

+mprintf('\n\tThe Dimensions of the reactor are ');

+mprintf('\n\tDiameter:%0.0f m',D_t); 

+mprintf('\n\tHeight:%0.2f m',h); 

+

+//OUTPUT (Ex7.5.c)

+mprintf('\n\n\n OUTPUT Ex7.5.c');

+mprintf('\n==========================================================');

+mprintf('\n\tThe new dimensions of the reactor');

+mprintf('\n\tDiameter:%0.1f m',D_t_c);

+mprintf('\n\tHeight:%0.1f m',h_new);

+

+//FILE OUTPUT

+fid= mopen('.\Chapter7-Ex5-Output.txt','w');

+mfprintf(fid,'\n OUTPUT Ex7.5.a');

+mfprintf(fid,'\n==========================================================');

+mfprintf(fid,'\n\tThe maximum rate of CO absorption at 15 atm : %f (mol/L s)',r_max);

+mfprintf(fid,'\n\tThe kinetic rate of CO absorption at 180(°C) : %f (mol/L s)',r_test);

+mfprintf(fid,'\n\tThe predicted value of k_L_a : %0.2f (s-1)',k_L_a);

+mfprintf(fid,'\n\n\n OUTPUT Ex7.5.b');

+mfprintf(fid,'\n==========================================================');

+mfprintf(fid,'\n\tThe Dimensions of the reactor are ');

+mfprintf(fid,'\n\tDiameter:%0.0f m',D_t); 

+mfprintf(fid,'\n\tHeight:%0.2f m',h); 

+mfprintf(fid,'\n\n\n OUTPUT Ex7.5.c');

+mfprintf(fid,'\n==========================================================');

+mfprintf(fid,'\n\tThe new dimensions of the reactor');

+mfprintf(fid,'\n\tDiameter:%0.1f m',D_t_c);

+mfprintf(fid,'\n\tHeight:%0.1f m',h_new);

+mclose(fid);

+

+//=================================================END OF PROGRAM===========================================================

+

+

+