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/CH8/EX8.1/Chapter8_Ex1.sce | 37 +++++++++++++++++++++++++++++++++++++
 1 file changed, 37 insertions(+)
 create mode 100644 1040/CH8/EX8.1/Chapter8_Ex1.sce

(limited to '1040/CH8/EX8.1')

diff --git a/1040/CH8/EX8.1/Chapter8_Ex1.sce b/1040/CH8/EX8.1/Chapter8_Ex1.sce
new file mode 100644
index 000000000..7dd99cd18
--- /dev/null
+++ b/1040/CH8/EX8.1/Chapter8_Ex1.sce
@@ -0,0 +1,37 @@
+//Harriot P., 2003, Chemical Reactor Design (I-Edition), Marcel Dekker, Inc., USA, pp 436.
+//Chapter-8 Ex8.1 Pg No. 323
+//Title:Gas absorption coefficient and fraction of overall resistance
+//============================================================================================================
+clear
+clc
+//INPUT
+rho_oil=0.8;//Density of oil (g/cm3)
+IV_init=130;//Iodine Value initial
+IV_final=80;//Iodine Value final
+P=45;//Pressure of system (psig)
+T_C=204;// Temperature of system (°C)
+t_run=[26 17];//Time required for hydrogenation run 2;
+frac_Ni=[0.005 0.0125]//Fraction of Nickel used for different run
+
+//CALCULATION
+r_ave=((IV_init -IV_final))*(0.039*rho_oil)*(1/60).*(t_run.^(-1));//Relationship between Iodine value and Hydrogen consumption (mol- H2/ L sec)
+H_H2= 4*10^(-3);//Solubility of H2 from Fig8.4 Pg No.322
+P_H2=(P/14.7)+1;//Absolute Pressure in (atm)
+C_H2=P_H2 *H_H2;
+Ci_by_r=C_H2.*(r_ave.^(-1));
+Coeff_R_cat=frac_Ni.^(-1);
+equation=[ones(1,2);Coeff_R_cat]//Simultaneous Equation
+Resistance= Ci_by_r*inv(equation);
+Gas_abs_resistance=(Resistance(1)*100 ).*(Ci_by_r.^(-1));
+Gas_abs_coefficient=(1/Resistance(1));
+
+//OUTPUT
+mprintf('\nThe Gas absorption coefficient is %f sec-1',Gas_abs_coefficient);
+mprintf('\n The Fraction of overall resistance due to gas absorption\n  Run 1 %0.0f%% \n  Run 2 %0.0f%%',Gas_abs_resistance(1),Gas_abs_resistance(2));
+
+//FILE OUTPUT
+fid= mopen('.\Chapter8-Ex1-Output.txt','w');
+mfprintf(fid,'\nThe Gas absorption coefficient is %f sec-1',Gas_abs_coefficient);
+mfprintf(fid,'\n The Fraction of overall resistance due to gas absorption\n  Run 1 %0.0f%% \n  Run 2 %0.0f%%',Gas_abs_resistance(1),Gas_abs_resistance(2));
+mclose(fid);
+//=======================================================END OF PROGRAM=================================================
-- 
cgit