diff options
| author | prashantsinalkar | 2017-10-10 12:27:19 +0530 |
|---|---|---|
| committer | prashantsinalkar | 2017-10-10 12:27:19 +0530 |
| commit | 7f60ea012dd2524dae921a2a35adbf7ef21f2bb6 (patch) | |
| tree | dbb9e3ddb5fc829e7c5c7e6be99b2c4ba356132c /1040/CH7/EX7.1/Chapter7_Ex1.sce | |
| parent | b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (diff) | |
| download | Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.tar.gz Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.tar.bz2 Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.zip | |
initial commit / add all books
Diffstat (limited to '1040/CH7/EX7.1/Chapter7_Ex1.sce')
| -rw-r--r-- | 1040/CH7/EX7.1/Chapter7_Ex1.sce | 78 |
1 files changed, 78 insertions, 0 deletions
diff --git a/1040/CH7/EX7.1/Chapter7_Ex1.sce b/1040/CH7/EX7.1/Chapter7_Ex1.sce new file mode 100644 index 000000000..d3f4bfd49 --- /dev/null +++ b/1040/CH7/EX7.1/Chapter7_Ex1.sce @@ -0,0 +1,78 @@ +//Harriot P.,2003,Chemical Reactor Design (I-Edition) Marcel Dekker,Inc.,USA,pp 436.
+//Chapter-7 Ex7.1 Pg No.260
+//Title:Overall Reaction Rate Coefficient, Percent Resistance, Reaction Volume and Reactor Size
+//===========================================================================================================
+clear
+clc
+// COMMON 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
+percent_inc=0.2;//Percentage excess required for reactor volume
+
+
+//CALCULATION (Ex7.1.a)
+H_O2_conv=H_O2*18/1000;// Convert (atm L/mole O2)
+k_L_a_by_H=k_L_a/H_O2_conv;
+reaction_resistance=H_O2_conv/(k2*C_B_product*(1-Epsilon)*3600);
+Kg_a=1/((1/k_g_a)+(1/k_L_a_by_H)+(reaction_resistance));//Refer equation7.10
+gasfilm_resistance_per=((1/k_g_a)/(1/Kg_a))*100;
+liq_film_resistance_per=((1/k_L_a_by_H)/(1/Kg_a))*100;
+reaction_resistance_per=((reaction_resistance)/(1/Kg_a))*100;
+
+//CALCULATION (Ex7.1.b)
+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 (Ex7.1.a)
+mprintf('\n OUTPUT Ex7.1.a');
+mprintf('\n==========================================================');
+mprintf('\nThe percentage gas-film resistance : %0.1f%%',gasfilm_resistance_per);
+mprintf('\nThe percentage liquid-film resistance: %0.1f%%',liq_film_resistance_per);
+mprintf('\nThe percentage chemical reaction resistance: %0.1f%%',reaction_resistance_per);
+
+//OUTPUT (Ex7.1.b)
+mprintf('\n\n\n OUTPUT Ex7.1.b');
+mprintf('\n==========================================================');
+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-Output.txt','w');
+mfprintf(fid,'\n OUTPUT Ex7.1.a');
+mfprintf(fid,'\n==========================================================');
+mfprintf(fid,'\nThe percentage gas-film resistance : %0.1f%%',gasfilm_resistance_per);
+mfprintf(fid,'\nThe percentage liquid-film resistance: %0.1f%%',liq_film_resistance_per);
+mfprintf(fid,'\nThe percentage chemical reaction resistance: %0.1f%%',reaction_resistance_per);
+mfprintf(fid,'\n\n\n OUTPUT Ex7.1.b');
+mfprintf(fid,'\n==========================================================');
+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======================================================
+
+
|