Merge pull request #1 from prashantsinalkar/masterHEADmaster

Initial commit

1 files changed, 255 insertions, 0 deletions

diff --git a/Chemical_Reactor_Design_by_P_Harriott/10-Novel_Reactors_.ipynb b/Chemical_Reactor_Design_by_P_Harriott/10-Novel_Reactors_.ipynb
new file mode 100644
index 0000000..d973c25
--- /dev/null
+++ b/Chemical_Reactor_Design_by_P_Harriott/10-Novel_Reactors_.ipynb
@@ -0,0 +1,255 @@
+{
+"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
+}