From 7f60ea012dd2524dae921a2a35adbf7ef21f2bb6 Mon Sep 17 00:00:00 2001 From: prashantsinalkar Date: Tue, 10 Oct 2017 12:27:19 +0530 Subject: initial commit / add all books --- 1040/CH4/EX4.2.b/Chapter4_Ex2_b.sce | 39 ++++++++++++++++++++++++++++++ 1040/CH4/EX4.2.b/Chapter4_Ex2_b_Output.txt | 4 +++ 2 files changed, 43 insertions(+) create mode 100644 1040/CH4/EX4.2.b/Chapter4_Ex2_b.sce create mode 100644 1040/CH4/EX4.2.b/Chapter4_Ex2_b_Output.txt (limited to '1040/CH4/EX4.2.b') diff --git a/1040/CH4/EX4.2.b/Chapter4_Ex2_b.sce b/1040/CH4/EX4.2.b/Chapter4_Ex2_b.sce new file mode 100644 index 000000000..3af409740 --- /dev/null +++ b/1040/CH4/EX4.2.b/Chapter4_Ex2_b.sce @@ -0,0 +1,39 @@ +//Harriot P.,2003,Chemical Reactor Design (I-Edition) Marcel Dekker,Inc., USA,pp 436 +//Chapter-4 Ex4.2.b Pg No.140 +//Title:Computation of surface mean pore radius, Diffusivity and tortusity +//=========================================================================================================== +clear +clc +//INPUT +S_g=150;//Total surface per gram (m2/g) +V_g=0.45;//Pore volume per gram (cm3/g) +V_i=0.30;//Micropore volume per gram (cm3/g) +V_a=0.15;// Macropore volume per gram (cm3/g) +rho_P=1.2;//Density of particle (g/cm3) +tau=2.5; +r_bar_i=40*(10^(-8));//Micropore radius +r_bar_a=2000*(10^(-8));//Macropore radius +D_AB=0.49;//For N2–O2 at 1 atm +M_O2=32;//Molecular weight of O2 +T=493;//Opereating Temperature +D_e=0.0235;//Refer Ex4.2a (cm2/s) Pg. No. 141 + +//CALCULATION +Epsilon=V_g*rho_P; +r_bar=2*V_g/(S_g*10^4); +D_K=9700*(r_bar)*sqrt(T/M_O2);//Knudsen Flow +D_Pore=1/((1/D_K)+(1/D_AB)); +tau=D_Pore*Epsilon/D_e; + +//OUTPUT +//Console Output +mprintf('\n The calculated surface mean pore radius = %.0e cm',r_bar); +mprintf('\n The predicted pore diffusivity = %0.2e cm2/sec',D_Pore); +mprintf('\n The corresponding tortusity = %0.2f',tau); +//File Output +fid= mopen('.\Chapter4_Ex2_b_Output.txt','w'); +mfprintf(fid,'\n The calculated surface mean pore radius = %.0e cm',r_bar); +mfprintf(fid,'\n The predicted pore diffusivity = %0.2e cm2/sec',D_Pore); +mfprintf(fid,'\n The corresponding tortusity = %0.2f',tau); +mclose(fid); +//===========================================END OF PROGRAM================================================== diff --git a/1040/CH4/EX4.2.b/Chapter4_Ex2_b_Output.txt b/1040/CH4/EX4.2.b/Chapter4_Ex2_b_Output.txt new file mode 100644 index 000000000..7dccc0924 --- /dev/null +++ b/1040/CH4/EX4.2.b/Chapter4_Ex2_b_Output.txt @@ -0,0 +1,4 @@ + + The calculated surface mean pore radius = 6e-07 cm + The predicted pore diffusivity = 2.18e-02 cm2/sec + The corresponding tortusity = 0.50 \ No newline at end of file -- cgit