1 files changed, 32 insertions, 0 deletions
diff --git a/1040/CH4/EX4.1.a/Chapter4_Ex1_a.sce b/1040/CH4/EX4.1.a/Chapter4_Ex1_a.sce
new file mode 100644
index 000000000..34ff568da
--- /dev/null
+++ b/1040/CH4/EX4.1.a/Chapter4_Ex1_a.sce
@@ -0,0 +1,32 @@
+//Harriot P.,2003,Chemical Reactor Design (I-Edition) Marcel Dekker,Inc., USA,pp 436.
+//Chapter-4 Ex4.1.a Pg No. 135
+//Title:Diffusivity of Chlorine at 573K and 1.5atm
+//===========================================================================================================
+clear
+clc
+//INPUT
+S_g=235;//Total surface per gram (m2/g)
+V_g=0.29E-6;//Pore volume per gram (cm3/g)
+rho_p=1.41;
+D_He=0.0065;//Effective diffusivity of He (cm2/sec)
+M_He=4;//Molecular weight of He
+M_Cl2=70.09;//Molecular weight of Cl2
+T_ref=293;//Reference temperature
+T_degC=300;
+T=T_degC+273;//Reaction temperature(K)
+
+//CALCULATION
+r_bar=2*V_g/S_g;//Mean Pore radius
+D_Cl2=D_He*((M_He/M_Cl2)*(T/T_ref))^(0.5);//Assuming Knudsen flow at 573K
+
+//OUTPUT
+//Console Output
+mprintf('The predicted diffusivity of Chlorine is %0.2E cm2/s ',D_Cl2);
+
+//File Output
+fid= mopen('.\Chapter4_Ex1_a_Output.txt','w');
+mfprintf(fid,'The predicted diffusivity of Chlorine is %0.2E cm2/s ',D_Cl2);
+mclose(fid);
+//============================================END OF PROGRAM=================================================
+
+