1 files changed, 33 insertions, 0 deletions
diff --git a/1040/CH4/EX4.1.b/Chapter4_Ex1_b.sce b/1040/CH4/EX4.1.b/Chapter4_Ex1_b.sce
new file mode 100644
index 000000000..d35f481e6
--- /dev/null
+++ b/1040/CH4/EX4.1.b/Chapter4_Ex1_b.sce
@@ -0,0 +1,33 @@
+//Harriot P.,2003,Chemical Reactor Design (I-Edition) Marcel Dekker,Inc., USA,pp 436.
+//Chapter-4 Ex4.1.b Pg No. 135
+//Title: Computation of tortuosity in catalyst pellet
+//===========================================================================================================
+clear
+clc
+//INPUT
+S_g=235;//Total surface per gram (m2/g)
+V_g=0.29;//Pore volume per gram (cm3/g)
+D_AB=0.73;// at 1atm and 298K
+rho_p=1.41;//Density of particle (g/cm3)
+D_He=0.0065;//Effective diffusivity of He (cm2/sec)
+T_ref=293;//Reference temperature (K)
+M_He=4;//Molecular weight of He
+T=298;//Operating temperature
+
+//CALCULATION
+r_bar=2*V_g /(S_g *(10^4));
+D_K=9700*(r_bar)*(T_ref/M_He)^(0.5);//Knudsen flow
+D_AB1=D_AB*(293/298)^(1.7)// at 1.5 atm and 293K
+D_pore=1/((1/D_K)+(1/D_AB1));//pore diffusion
+Epsilon=V_g*rho_p;
+tau=(D_pore*Epsilon)/D_He;//Tortusity
+
+//OUTPUT
+//Console Output
+mprintf('\n The tortusity value = %0.2f',tau);   
+
+//File Output
+fid= mopen('.\Chapter4_Ex1_b_Output.txt','w');
+mfprintf(fid,'\n The tortusity value = %0.2f',tau);
+mclose(fid);
+//========================================================END OF PROGRAM=================================