1 files changed, 51 insertions, 0 deletions
diff --git a/1040/CH7/EX7.1.b/Chapter7_Ex1_b.sce b/1040/CH7/EX7.1.b/Chapter7_Ex1_b.sce
new file mode 100644
index 000000000..882f08df8
--- /dev/null
+++ b/1040/CH7/EX7.1.b/Chapter7_Ex1_b.sce
@@ -0,0 +1,51 @@
+//Harriot P.,2003,Chemical Reactor Design (I-Edition) Marcel Dekker,Inc.,USA,pp 436. 
+//Chapter-7 Ex7.1.b Pg No.260
+//Title:Reaction Volume and Reactor Size 
+//===========================================================================================================
+clear
+clc
+//INPUT
+k2=8.5;//Reaction rate constant (L/mol-sec)
+T=50;//Reaction condition temperature(°C)
+P=2;//Reaction Pressure (atm)
+H_O2=8*10^4;// Solubility (atm/mol fraction)
+F=17000//Feed rate (L/hr)
+C_B_feed=1.6;//Feed concentration(M)
+C_B_product=0.8;//Product concentration(M)
+k_L_a=900;//Liquid film mass transfer coefficient(hr-1)
+k_g_a=80;//Gas film mass transfer coefficient(mol/hr L atm)
+Epsilon=0.1;//Porosity
+Kg_a=0.596;//Refer the overall reaction rate calculated in Ex7.1.a
+percent_inc=0.2;//Percentage excess required for reactor volume
+
+//CALCULATION
+delta_C_B=C_B_feed-C_B_product;
+mol_O2_needed=F*delta_C_B/4;
+N_air=100;//Assuming 100 mole of feed air
+f_O2=0.209;//Fraction of O2
+f_N2=1-f_O2;//Fraction of N2
+N_O2_in=N_air*f_O2;
+N_N2_in=N_air*f_N2;
+N_O2_out=N_O2_in/2;//Half of O2 fed
+N_N2_out=N_N2_in;
+N_air_out=N_N2_out+N_O2_out;
+P_O2_out=P*(N_O2_out/N_air_out);
+P_O2_in=P*(N_O2_in/N_air);
+P_O2_bar=(P_O2_in-P_O2_out)/(log(P_O2_in/P_O2_out));//Log mean Pressure
+volume=mol_O2_needed/(Kg_a*P_O2_bar);
+reactor_vol=volume+volume*percent_inc;
+volume_gal=volume*0.264;
+reactor_vol_gal=reactor_vol*0.264;
+
+//OUTPUT
+//Console Output
+mprintf('\n Reaction volume  calculated : %0.0f L ',volume );
+mprintf('\n Reactor size to be chosen : %0.0f L',reactor_vol);
+//File Output
+fid= mopen('.\Chapter7_Ex1_b_Output.txt','w');
+mfprintf(fid,'\n Reaction volume  calculated : %0.0f L ',volume );
+mfprintf(fid,'\n Reactor size to be chosen : %0.0f L',reactor_vol);
+mclose(fid);
+//=============================================END OF PROGRAM============================================================
+// Disclaimer : The numerically calculated value of reaction volume is 18008 L not 18000 L as mentioned in the textbook
+