1 files changed, 57 insertions, 0 deletions
diff --git a/1040/CH4/EX4.4/Chapter4_Ex4.sce b/1040/CH4/EX4.4/Chapter4_Ex4.sce
new file mode 100644
index 000000000..942eae0cb
--- /dev/null
+++ b/1040/CH4/EX4.4/Chapter4_Ex4.sce
@@ -0,0 +1,57 @@
+//Harriot P.,2003,Chemical Reactor Design (I-Edition) Marcel Dekker,Inc. USA,pp 436.
+//Chapter-4 Ex4.4 Pg No.157
+//Title: Effectiveness factor for solid catalyzed reaction
+//======================================================================================================================
+clear
+clc
+//INPUT
+D_e_A=0.02;//(cm2/s)
+D_e_B=0.03;//(cm2/s)
+D_e_C=0.015;//(cm2/s)
+X_f_A=0.3;
+X_f_B=(1-X_f_A);
+eta_assumed=0.68;//Effectiveness factor from Fig.4.8 for first order reaction
+T=150;//(deg C)
+T_K=T+273;//(K)
+r=0.3;//(cm)Radius of catalyst sphere
+P_opt=4;//(atm)Operating Pressure 
+R=82.056;//(cm3 atm/K mol)Gas constant 
+
+
+//CALCULATION
+//Kinetic equation r= (2.5*10^-5*P_A*P_B)/(1+0.1*P_A+2*P_C)^2
+P_A=X_f_A*P_opt;
+P_B=X_f_B*P_opt;
+r_star=(2.5*10^-5*P_A*P_B)/(1+0.1*P_A)^2;
+C_A=P_A/(R*T_K);
+k=r_star/C_A;
+Phi= r*(k/D_e_A)^(0.5);
+P_A_bar=eta_assumed*P_A;
+delta_P_A=P_A*(1-eta_assumed);
+delta_P_B=delta_P_A*(D_e_A/D_e_B);
+P_B_bar=P_B-delta_P_B;
+delta_P_C=delta_P_A*(D_e_A/D_e_C);
+P_C_bar=delta_P_C;
+r_calc=(2.5*10^-5*P_A_bar*P_B_bar)/(1+0.1*P_A_bar+2*P_C_bar)^2
+eta_calc=r_calc/r_star;
+eta_approx=(eta_calc+eta_assumed)/2;
+
+//OUTPUT
+//Console Output
+mprintf('\tBased on average pressures calculated Rate and Effectiveness factor');
+mprintf('\n\t r : %0.2E (mol/s cm3)',r_calc);
+mprintf('\n\t eta_calc : %0.3f ',eta_calc);
+mprintf('\n The actual value of Effectiveness factor eta_actual :%0.1f',eta_approx);
+
+//File Output
+fid= mopen('.\Chapter4-Ex4-Output.txt','w');
+mfprintf(fid,'\tBased on average pressures calculated Rate and Effectiveness factor');
+mfprintf(fid,'\n\t r : %0.2E (mol/s cm3)',r_calc);
+mfprintf(fid,'\n\t eta_calc : %0.3f ',eta_calc);
+mfprintf(fid,'\n The actual value of Effectiveness factor eta_actual :%0.1f',eta_approx);
+mclose(fid);
+//================================================END OF PROGRAM==================================================================================
+
+
+
+