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/CH9/EX9.1/Chapter9_Ex1.sce | 48 +++++++++++++++++++++++++++++++++++++++++
 1 file changed, 48 insertions(+)
 create mode 100644 1040/CH9/EX9.1/Chapter9_Ex1.sce

(limited to '1040/CH9/EX9.1')

diff --git a/1040/CH9/EX9.1/Chapter9_Ex1.sce b/1040/CH9/EX9.1/Chapter9_Ex1.sce
new file mode 100644
index 000000000..a43fd6b50
--- /dev/null
+++ b/1040/CH9/EX9.1/Chapter9_Ex1.sce
@@ -0,0 +1,48 @@
+//Harriot P.,2003,Chemical Reactor Design (I-Edition) Marcel Dekker,Inc.,USA,pp 436.
+//Chapter-9 Ex9.1 Pg No.376
+//Title:Model II- Volumetric Mass Transfer Coefficient (K)
+//============================================================================================================
+
+clear
+clc
+//INPUT
+u0=[ 0.1 0.3 0.5 0.75 0.95 1.15];//Fluid Velocities (m/sec)
+X=[0.923 0.872 0.846 0.775 0.728 0.664];//Conversion
+h_by_h0=[1.26 1.44 1.66 2.0 2.3 2.7];//Height of bed under fluidized condition by height of packed bed
+Epsilon_m=0.456;//Fraction of voids in packed bed
+h0=0.75;//Height of packed bed (m)
+k_r=4.45 ;//Reaction rate constant(sec-1)
+W=5;//Weight of the bed (kg)
+
+
+//CALCULATION
+n=length(X);
+for i=1:n
+    K0_L_by_u0(i)=log(1/(1-X(i)));//Refer equation 9.21 Pg No.371
+    L(i)=h_by_h0(i)*h0;
+    one_minus_epsilon(i)=(1-Epsilon_m)/h_by_h0(i); 
+    k_rhob(i)=k_r*one_minus_epsilon(i); 
+    K0(i)=K0_L_by_u0(i)*u0(i)/L(i);
+    K(i)=1/((K0(i).^(-1))-(1/k_rhob(i))); //Refer equation 9.19 Pg No.371
+end
+
+
+//OUTPUT
+mprintf('\nThe values of K for given velocities')
+mprintf('\n u (m/sec) \t K (sec-1) ');
+mprintf('\n===================================================================');
+for i=1:n
+    mprintf('\n %.3g \t \t %0.3f',u0(i),K(i));            
+end
+
+//FILE OUTPUT
+fid= mopen('.\Chapter9-Ex1-Output.txt','w');
+mfprintf(fid,'\nThe values of K for given velocities')
+mfprintf(fid,'\n u (m/sec) \t K (sec-1) ');
+mfprintf(fid,'\n===================================================================');
+for i=1:n
+    mfprintf(fid,'\n %.3g \t \t %0.3f',u0(i),K(i));            
+end    
+
+//==============================================END OF PROGRAM=============================================    
+    
-- 
cgit