diff options
Diffstat (limited to 'Chemical_Reactor_Design_by_P_Harriott/10-Novel_Reactors_.ipynb')
| -rw-r--r-- | Chemical_Reactor_Design_by_P_Harriott/10-Novel_Reactors_.ipynb | 255 |
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 +} |