From b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b Mon Sep 17 00:00:00 2001 From: priyanka Date: Wed, 24 Jun 2015 15:03:17 +0530 Subject: initial commit / add all books --- 1040/CH8/EX8.1/Chapter8_Ex1_Output.txt | 5 +++++ 1040/CH8/EX8.1/Ex8_1.sce | 37 ++++++++++++++++++++++++++++++++++ 2 files changed, 42 insertions(+) create mode 100644 1040/CH8/EX8.1/Chapter8_Ex1_Output.txt create mode 100644 1040/CH8/EX8.1/Ex8_1.sce (limited to '1040/CH8/EX8.1') diff --git a/1040/CH8/EX8.1/Chapter8_Ex1_Output.txt b/1040/CH8/EX8.1/Chapter8_Ex1_Output.txt new file mode 100644 index 000000000..533ea8272 --- /dev/null +++ b/1040/CH8/EX8.1/Chapter8_Ex1_Output.txt @@ -0,0 +1,5 @@ + +The Gas absorption coefficient is 0.145500 sec-1 + The Fraction of overall resistance due to gas absorption + Run 1 42% + Run 2 65% \ No newline at end of file diff --git a/1040/CH8/EX8.1/Ex8_1.sce b/1040/CH8/EX8.1/Ex8_1.sce new file mode 100644 index 000000000..7dd99cd18 --- /dev/null +++ b/1040/CH8/EX8.1/Ex8_1.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