summaryrefslogtreecommitdiff
path: root/1040/CH7/EX7.4/Chapter7_Ex4.sce
diff options
context:
space:
mode:
Diffstat (limited to '1040/CH7/EX7.4/Chapter7_Ex4.sce')
-rw-r--r--1040/CH7/EX7.4/Chapter7_Ex4.sce52
1 files changed, 52 insertions, 0 deletions
diff --git a/1040/CH7/EX7.4/Chapter7_Ex4.sce b/1040/CH7/EX7.4/Chapter7_Ex4.sce
new file mode 100644
index 000000000..7b8932178
--- /dev/null
+++ b/1040/CH7/EX7.4/Chapter7_Ex4.sce
@@ -0,0 +1,52 @@
+//Harriot P.,2003,Chemical Reactor Design (I-Edition) Marcel Dekker,Inc.,USA,pp 436
+//Chapter-7 Ex7.4 Pg No.279
+//Title:Local selectivity due to mass transfer limitations
+//===========================================================================================================
+clear
+clc
+//INPUT
+C_Ai=0.02;//(M)
+C_B0=3;//(M)
+D_A=10^(-5);//(cm2/sec)
+D_B=D_A;//(cm2/sec)
+D_C=D_B;//(cm2/sec)
+k_1=10^(4);//(L/mol sec)
+k_star_l=0.015;//(cm/sec)
+n=1;
+C_c0=[0 1.4];
+X=[0 0.5]// Conversion
+Phi=[33 23];//From figure 7.7
+
+
+//CALCULATION
+k_2=0.09*k_1;
+for i=1:2
+ C_B(i)=(1-X(i))*C_B0;
+sqrt_M(i)=sqrt(C_B(i)*k_1*D_A)/k_star_l;
+Phi_a(i)=(1+(C_B(i)/(n*C_Ai))*(D_B/D_A));//Refer equation 7.51
+C_Bbar_by_C_B(i)=(Phi(i)/sqrt_M(i))^2;//Refer equation 7.59
+delta_C_B(i)=(1-C_Bbar_by_C_B(i))*C_B(i);//Refer equation 7.60
+delta_C_c(i)=delta_C_B(i);
+C_cbar(i)=delta_C_c(i)+C_c0(i);
+C_Bbar(i)=C_Bbar_by_C_B(i)*(C_B(i));
+S(i)=(1-(k_2*C_cbar(i)/(C_Bbar(i)*k_1)))*100;//Refer equation 7.56
+end
+
+//OUTPUT
+mprintf('\n\tLocal selectivity due to mass transfer limitations ');
+mprintf('\n\tThe local selectivity for Zero Conversion : %0.0f%%',S(1));
+mprintf('\n\tThe local selectivity for 50%% Conversion : %0.0f%%',S(2));
+
+//FILE OUTPUT
+fid= mopen('.\Chapter7-Ex4-Output.txt','w');
+mfprintf(fid,'\n\tLocal selectivity due to mass transfer limitations ');
+mfprintf(fid,'\n\tThe local selectivity for Zero Conversion is %0.0f%%',S(1));
+mfprintf(fid,'\n\tThe local selectivity for 50%% Conversion is %0.0f%%',S(2));
+mclose(fid);
+//======================================================END OF PROGRAM===================================================
+
+
+
+
+
+