{
"cells": [
 {
		   "cell_type": "markdown",
	   "metadata": {},
	   "source": [
       "# Chapter 10: Novel Reactors "
	   ]
	},
{
		   "cell_type": "markdown",
		   "metadata": {},
		   "source": [
			"## Example 10.1: Fraction_unconverted_naphthalene_based_on_model_II.sce"
		   ]
		  },
  {
"cell_type": "code",
	   "execution_count": null,
	   "metadata": {
	    "collapsed": true
	   },
	   "outputs": [],
"source": [
"//Harriot P., 2003, Chemical Reactor Design (I-Edition), Marcel Dekker, Inc., USA, pp 436.\n",
"//Chapter-10 Ex10.1 Pg No. 408\n",
"//Title:Fraction unconverted naphthalene based on model II\n",
"//===========================================================================================================\n",
"clear\n",
"clc\n",
"//INPUT\n",
"T_ref=273;//Reference Temperature\n",
"T_feed=300+T_ref;//Temperature in (K)\n",
"SV_STP=[60000 120000];//Space velocity(Hr-1)\n",
"t_cell=0.04;//Thickness(cm)\n",
"cell_unit_area=100/(2.54^2);//No of cells per unit area(cells/cm2)\n",
"L_inch=6;// Length of monolithic converter (Inches)\n",
"Epsilon=0.68;//Porosity\n",
"myu=0.0284*(10^-2);//Viscosity of air(Poise)\n",
"rho=6.17*10^(-4);//Density of air (g/cm3)\n",
"\n",
"\n",
"//CALCULATION\n",
"d=sqrt(1/cell_unit_area)- t_cell;\n",
"Epsilon=(d^2/(d+t_cell)^2);\n",
"\n",
"//Assume the wash coating lowers d to 0.21 cm and Epsilon to 0.68:\n",
"d_new=0.21;\n",
"Epsilon_new =0.68\n",
"a=4*Epsilon_new/d_new;\n",
"SV=SV_STP.*(T_feed/(T_ref*3600));//Refer equation 10.13\n",
"L_cm=L_inch*2.54;\n",
"u0=SV.*(L_cm);\n",
"u=u0.*(1/Epsilon);\n",
"Nu=myu/rho;//Kinematic viscosity\n",
"D_CO_N2_1=0.192;//Diffusion coefficients for binary gas mixtures(cm2/sec) at 288K\n",
"D_CO_N2_2=D_CO_N2_1*(T_feed/288)^(1.7);////Diffusion coefficients for binary gas mixtures(cm2/sec) at 573K\n",
"Sc=Nu/D_CO_N2_2;\n",
"for i=1:2\n",
"Re(i)=d_new*u(i)/Nu;\n",
"Re_Sc_d_by_L(i)=Re(i)*Sc*(d_new/L_cm);\n",
"Sh(i) = 3.66 *(1+0.095*Re_Sc_d_by_L(i))^(0.45);//Refer equation 10.7\n",
"k_c(i)=Sh(i)*D_CO_N2_2/d_new;\n",
"X(i)=1-exp((-k_c(i)*a*L_cm*u0(i)^(-1)));//Refer equation10.12\n",
"Percent_X(i)=X(i)*100;\n",
"end\n",
"\n",
"//OUTPUT\n",
"mprintf('\n The Conversion expected for the given space velocities ');\n",
"mprintf(' \n Space Velocity (hr-1)\t \t Conversion (%%)');\n",
"mprintf('\n ======================================================');\n",
"for i=1:2\n",
"    mprintf('\n %.0f \t \t \t \t %.1f',SV_STP(i),Percent_X(i));\n",
"end\n",
"\n",
"//FILE OUTPUT\n",
"fid= mopen('.\Chapter10-Ex1-Output.txt','w');\n",
"mfprintf(fid,'\n The Conversion expected for the given space velocities ');\n",
"mfprintf(fid,' \n Space Velocity (hr-1)\t \t Conversion (%%)');\n",
"mfprintf(fid,'\n ======================================================');\n",
"for i=1:2\n",
"    mfprintf(fid,'\n %.0f \t \t \t \t %.1f',SV_STP(i),Percent_X(i));\n",
"end\n",
"mclose(fid);\n",
"\n",
"\n",
"//================================================END OF PROGRAM=========================================================\n",
"\n",
"\n",
"\n",
"\n",
"\n",
"\n",
"\n",
"\n",
"\n",
"\n",
"\n",
"\n",
"\n",
"\n",
"\n",
"\n",
""
   ]
   }
,
{
		   "cell_type": "markdown",
		   "metadata": {},
		   "source": [
			"## Example 10.2: Conversion_as_a_function_of_No_of_Gauzes.sce"
		   ]
		  },
  {
"cell_type": "code",
	   "execution_count": null,
	   "metadata": {
	    "collapsed": true
	   },
	   "outputs": [],
"source": [
"//Harriot P., 2003, Chemical Reactor Design (I-Edition), Marcel Dekker, Inc., USA, pp 436.\n",
"//Chapter-10 Ex10.2 Pg No. 414\n",
"//Title:Conversion as a function of No. of Gauzes\n",
"//===========================================================================================================\n",
"clear\n",
"clc\n",
"// COMMON INPUT\n",
"M_NH3=17;//Molecular weight NH3\n",
"M_air=29;//Molecular weight air\n",
"f_air=0.9;//Fraction of air in feed\n",
"f_NH3=(1-f_air);//Fraction of NH3 in feed\n",
"myu_air=0.0435*(10^-2);//Viscosity of air (Poise)\n",
"P_atm=(100+14.7)/14.7;//Pressure of the system\n",
"P_ref=1;//Reference Pressure\n",
"T_ref=273;//Reference temperature\n",
"T_inlet=300+T_ref;//Inlet Temperature\n",
"V_ref=22400;\n",
"T_surf=700+T_ref;//Surface Temperature\n",
"u0=1.8;//Velocity at 300 °C (m/sec)\n",
"d=0.076*(10^-1);//Size of wire (cm)\n",
"D_NH3_N2=0.23;//Diffusivity at 298 K 1 atm(cm2/s)\n",
"N=32;//Gauzes (wires/cm)\n",
"frac_N2 = 0.25*(10^(-2));//Fraction of NH3 fed into N2 (Byproduct reaction)\n",
"n =[1 2 5 10 15 20];//No. of Gauzes\n",
"\n",
"\n",
"//CALCULATION (Ex 10.2.a)\n",
"M_ave =f_air*M_air+f_NH3*M_NH3;\n",
"rho =(M_ave*T_ref*P_atm)/(V_ref*T_surf*P_ref);\n",
"u0_surf = u0*(T_surf/T_inlet);\n",
"Re = rho*u0_surf*100*d/myu_air;\n",
"Gamma = [1-32*(d)]^2;//From equation 10.5\n",
"Re_Gamma = Re/Gamma;\n",
"D_NH3 = 0.23*(T_surf/298)^(1.7)*(1/7.8);// at 7.8 atm 700 °C\n",
"Sc =(myu_air*P_ref)/(rho*D_NH3);\n",
"j_D = 0.644*(Re_Gamma)^(-0.57);//Refer equation 10.14\n",
"k_c = j_D*(u0_surf*100/Gamma)*(1/(Sc)^(2/3));\n",
"a_dash = 2*(%pi)*(d)*N\n",
"k_c_a_dash_u0 =(k_c*a_dash)/(u0_surf*100);\n",
"m = length(n)\n",
"for i = 1:m\n",
"    X(i) = (1-exp(-k_c_a_dash_u0*n(i)));\n",
"end\n",
"//CALCULATION (Ex 10.2.b)\n",
"for i = 1:m\n",
"    X(i) = (1-exp(-k_c_a_dash_u0*n(i)));\n",
"    Yield(i) = X(i)-frac_N2*n(i);\n",
"end\n",
"\n",
"\n",
"//OUTPUT(Ex 10.2.a)\n",
"mprintf('\n OUTPUT Ex10.2.a');\n",
"mprintf('\n=====================================');\n",
"mprintf('\n \tThe Ammonia Conversion');\n",
"mprintf('\n=====================================');\n",
"mprintf('\n\t Gauzes        Conversion');\n",
"mprintf('\n\t  (n)             (X)');\n",
"mprintf('\n=====================================');\n",
"for i=1:m\n",
"    mprintf('\n\t  %.0f  \t \t %.3f',n(i),X(i));\n",
"end\n",
"\n",
"//OUTPUT(Ex 10.2.b)\n",
"mprintf('\n\n\n OUTPUT Ex10.2.b');\n",
"mprintf('\n==========================================');\n",
"mprintf('\n \tThe Ammonia Yield');\n",
"mprintf('\n==========================================');\n",
"mprintf('\n\t Gauzes        Yield');\n",
"mprintf('\n\t  (n)             (X-%fn)',frac_N2);\n",
"mprintf('\n==========================================');\n",
"for i=1:m\n",
"    mprintf('\n\t  %.0f \t \t %.3f',n(i),Yield(i));\n",
"end\n",
"//FILE OUTPUT\n",
"fid= mopen('.\Chapter10-Ex2-Output.txt','w');\n",
"mfprintf(fid,'\n OUTPUT Ex10.2.a');\n",
"mfprintf(fid,'\n=====================================');\n",
"mfprintf(fid,'\n \tThe Ammonia Conversion');\n",
"mfprintf(fid,'\n=====================================');\n",
"mfprintf(fid,'\n\t Gauzes        Conversion');\n",
"mfprintf(fid,'\n\t  (n)             (X)');\n",
"mfprintf(fid,'\n=====================================');\n",
"for i=1:m\n",
"    mfprintf(fid,'\n\t  %.0f  \t \t %.3f',n(i),X(i));\n",
"end\n",
"mfprintf(fid,'\n\n\n OUTPUT Ex10.2.b');\n",
"mfprintf(fid,'\n==========================================');\n",
"mfprintf(fid,'\n \tThe Ammonia Yield');\n",
"mfprintf(fid,'\n==========================================');\n",
"mfprintf(fid,'\n\t Gauzes        Yield');\n",
"mfprintf(fid,'\n\t  (n)             (X-%fn)',frac_N2);\n",
"mfprintf(fid,'\n==========================================');\n",
"for i=1:m\n",
"    mfprintf(fid,'\n\t  %.0f \t \t %.3f',n(i),Yield(i));\n",
"end\n",
"mclose(fid);\n",
"\n",
"//====================================================END OF PROGRAM====================================================\n",
"\n",
"\n",
"\n",
"\n",
"\n",
"\n",
"\n",
"\n",
"\n",
""
   ]
   }
],
"metadata": {
		  "kernelspec": {
		   "display_name": "Scilab",
		   "language": "scilab",
		   "name": "scilab"
		  },
		  "language_info": {
		   "file_extension": ".sce",
		   "help_links": [
			{
			 "text": "MetaKernel Magics",
			 "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md"
			}
		   ],
		   "mimetype": "text/x-octave",
		   "name": "scilab",
		   "version": "0.7.1"
		  }
		 },
		 "nbformat": 4,
		 "nbformat_minor": 0
}