summaryrefslogtreecommitdiff
path: root/1040/CH7/EX7.5.a/Chapter7_Ex5_a.sce
diff options
context:
space:
mode:
Diffstat (limited to '1040/CH7/EX7.5.a/Chapter7_Ex5_a.sce')
-rw-r--r--1040/CH7/EX7.5.a/Chapter7_Ex5_a.sce57
1 files changed, 57 insertions, 0 deletions
diff --git a/1040/CH7/EX7.5.a/Chapter7_Ex5_a.sce b/1040/CH7/EX7.5.a/Chapter7_Ex5_a.sce
new file mode 100644
index 000000000..1ad78a946
--- /dev/null
+++ b/1040/CH7/EX7.5.a/Chapter7_Ex5_a.sce
@@ -0,0 +1,57 @@
+//Harriot P.,2003,Chemical Reactor Design (I-Edition) Marcel Dekker,Inc.,USA,pp 436
+//Chapter-7 Ex7.5.a Pg No.293
+//Title:Maximum rate of CO absorption
+//===========================================================================================================
+clear
+clc
+//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
+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 water at 20(°C)
+Conc_Rh=4*10^(-3);//Concentration of Rohdium(M)
+Conc_CH3I=1;//Concentration of Methyl Iodide(M)
+
+//CALCUATION
+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)
+
+//OUTPUT
+//Console Output
+mprintf('\n\tThe maximum rate of CO absorption at 15 atm : %0.3f (mol/L s)',r_max);
+mprintf('\n\tThe kinetic rate of CO absorption at 180(°C) : %0.3f (mol/L s)',r_test);
+mprintf('\n\tThe predicted value of k_L_a : %0.2f (s-1)',k_L_a);
+//File Output
+fid= mopen('.\Chapter7_Ex5_a_Output.txt','w');
+mfprintf(fid,'\n\tThe maximum rate of CO absorption at 20 atm : %0.3f (mol/L s)',r_max);
+mfprintf(fid,'\n\tThe kinetic rate of CO absorption at 180(°C) : %0.3f (mol/L s)',r_test);
+mfprintf(fid,'\n\tThe predicted value of k_L_a : %0.2f (s-1)',k_L_a);
+mclose(fid);
+//=================================================END OF PROGRAM===========================================================
+
+
+