diff options
Diffstat (limited to 'Mass_Transfer_Operations_by_R_E_Treybal/13-Leaching.ipynb')
| -rw-r--r-- | Mass_Transfer_Operations_by_R_E_Treybal/13-Leaching.ipynb | 356 |
1 files changed, 356 insertions, 0 deletions
diff --git a/Mass_Transfer_Operations_by_R_E_Treybal/13-Leaching.ipynb b/Mass_Transfer_Operations_by_R_E_Treybal/13-Leaching.ipynb new file mode 100644 index 0000000..1f5fc0e --- /dev/null +++ b/Mass_Transfer_Operations_by_R_E_Treybal/13-Leaching.ipynb @@ -0,0 +1,356 @@ +{ +"cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# Chapter 13: Leaching" + ] + }, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 13.1: Unsteady_State_Operation.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clear;\n", +"clc;\n", +"\n", +"// Illustration 13.1\n", +"// Page: 722\n", +"\n", +"printf('Illustration 13.1 - Page: 722\n\n');\n", +"\n", +"// Solution \n", +"\n", +"//***Data***//\n", +"Density_L = 1137;// [kg/cubic m]\n", +"Density_S = 960;// [kg/cubic m]\n", +"Density_p = 1762;// [kg/cubic m]\n", +"A_prime = 16.4;// [square m/kg]\n", +"g = 9.81;// [square m/s]\n", +"sigma = 0.066;// [N/m]\n", +"Z = 3;// [m]\n", +"dia = 1;// [m]\n", +"//**********//\n", +"\n", +"e = 1-(Density_S/Density_p);// [fraction void]\n", +"ap = A_prime*Density_S;// [square m/cubic m]\n", +"// By Eqn. 6.67:\n", +"dp = 6*(1-e)/ap;// [m]\n", +"// By Eqn. 13.6:\n", +"K = dp^2*e^3*g/(150*(1-e)^2);// [cubic m/s]\n", +"check = K*Density_L*g/(g*sigma);\n", +"if check<0.02\n", +" // By Eqn. 13.3: \n", +" So = 0.075;\n", +"else\n", +" // By Eqn. 13.4:\n", +" So = 0.0018/(check)\n", +"end\n", +"// By Eqn. 13.2:\n", +"ZD = (0.275/g)/((K/g)^0.5*(Density_L/sigma));// [m]\n", +"// By Eqn. 13.1:\n", +"Sav = ((Z-ZD)*So/Z)+(ZD/Z);\n", +"// VolRatio=Vol liquid retained/Vol bed.\n", +"VolRatio = Sav*e;\n", +"printf('Vol liquid retained/Vol bed : %f cubic m/cubic m\n',VolRatio);\n", +"Mass = VolRatio*%pi*dia^2*Z*Density_L/4;// [kg]\n", +"// Mass ratio=Mass Liquid/Mass dry solid\n", +"MassRatio = VolRatio*Density_L/(Density_S);\n", +"printf('Mass liquid/Mass dry solid: %f kg/kg\n',MassRatio);" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 13.2: Multistage_Crosscurrent_Leaching.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clear;\n", +"clc;\n", +"\n", +"// Illustration 13.2\n", +"// Page: 749\n", +"\n", +"printf('Illustration 13.2 - Page: 749\n\n');\n", +"\n", +"// Solution \n", +"\n", +"//***Data***//\n", +"// Eqb=[x(Wt fraction NaOH in clear solution) N(kg CaCO3/kg soln in settled sludge) y*(wt fraction NaOH in soln of settled sludge)]\n", +"// a=H2O b=CaCO3 c=NaOH\n", +"Eqb = [0.090 0.495 0.0917;0.0700 0.525 0.0762;0.0473 0.568 0.0608;0.0330 0.600 0.0452;0.0208 0.620 0.0295;0.01187 0.650 0.0204;0.00710 0.659 0.01435;0.00450 0.666 0.01015];\n", +"deff('[y]=f80(x)','y=0');\n", +"x = 0:0.01:0.12;\n", +"Mass_c = 0.1;// [kg]\n", +"Mass_b = 0.125;// [kg]\n", +"Mass_a = 0.9;// [kg]\n", +"//**************//\n", +"\n", +"scf(42);\n", +"plot(x,f80,Eqb(:,3),Eqb(:,2));\n", +"xgrid();\n", +"xlabel('x,y Wt. fraction of NaOH in loquid');\n", +"ylabel('N kg CaCO3 / kg solution');\n", +"legend('N Vs x','N Vs Y');\n", +"title('Equilibrium Plot')\n", +"// Basis: 1 kg soln in original mixture.\n", +"// As in Fig. 13.27 (Pg 750)\n", +"// The original mixture corresponds to M1:\n", +"NM1 = 0.125;// [kg CaCO3/kg soln]\n", +"yM1 = 0.1;// [kg NaOH/kg solution]\n", +"// The tie line through M1 is drawn. At point E1 representing the settled sludge:\n", +"N1 = 0.47;// [kg CaCO3/kg soln]\n", +"y1 = 0.100;// [kg NaOH/kg solution]\n", +"E1 = Mass_b/N1;// [kg soln. in sludge]\n", +"Ro = 1-E1;// [kg clear soln drawn]\n", +"\n", +"// Stage 2:\n", +"xo = 0;// [kg NaOH/kg soln]\n", +"// By Eqn. 13.11:\n", +"M2 = E1+Ro;// [kg liquid]\n", +"// By Eqn. 13.12:\n", +"NM2 = Mass_b/(E1+Ro);// [kg CaCO3/kg soln]\n", +"// M2 is located on line RoE1. At this value of N, and the tie line through M2 is drawn. At E2:\n", +"N2 = 0.62;// [kg CaCO3/kg soln]\n", +"y2 = 0.035;// [kg NaOH/kg solution]\n", +"E2 = Mass_b/N2;// [kg soln. in sludge]\n", +"Ro = 1-E2;// [kg clear soln drawn]\n", +"\n", +"// Stage 3:\n", +"xo = 0;// [kg NaOH/kg soln]\n", +"// By Eqn. 13.11:\n", +"M3 = E2+Ro;// [kg liquid]\n", +"// By Eqn. 13.12:\n", +"NM3 = Mass_b/M3;// [kg CaCO3/kg soln]\n", +"// Tie line E3R3 is located through M3.At E3:\n", +"N3 = 0.662;// [kg CaCO3/kg soln]\n", +"y3 = 0.012;// [kg NaOH/kg solution]\n", +"// By Eqn. 13.8:\n", +"E3 = Mass_b/N3;// [kg soln. in sludge]\n", +"printf('The fraction of original NaOH in the slurry: %f \n',E3*y3/Mass_c);" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 13.3: Multistage_Countercurrent_Leaching.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clear;\n", +"clc;\n", +"\n", +"// Illustration 13.3\n", +"// Page: 754\n", +"\n", +"printf('Illustration 13.3 - Page: 754\n\n');\n", +"\n", +"// Solution (a)\n", +"\n", +"//***Data***//\n", +"// a=H2O b=CaCO3 c=NaOH \n", +"mass_c = 400;// [kg/h]\n", +"x1 = 0.1;// [wt fraction NaOH in overflow]\n", +"//**********//\n", +"\n", +"Mb = 100;// [kg/kmol]\n", +"Mc = 40;// [kg/kmol]\n", +"rate_c = mass_c/Mc;// [kmol/h]\n", +"rate_b = rate_c/2;// [kmol/h]\n", +"mass_b = rate_b*Mb;// [kg/h]\n", +"// After trial calculations:\n", +"y3 = 0.01;// [kg NaOH/kg solution]\n", +"N3 = 0.666;// [kg CaCO3/kg solution]\n", +"E3 = mass_b/N3;// [kg/h]\n", +"lost_c = E3*y3;// [kg/h]\n", +"sludge_a = E3-lost_c;// [kg/h]\n", +"overflow_c = mass_c-lost_c;// [kg NaOH/kg solution]\n", +"R1 = overflow_c/x1;// [kg overflow/h]\n", +"R1_a = R1-overflow_c;// [kg/h]\n", +"RNpPlus1 = R1_a+sludge_a;// [kg/h]\n", +"// For purpose of calculation, it may be imagined that agitators are not present in the flowsheet and the first thickner is fed with the dry mixture of the reaction products, CaCO3 and NaOH, together with overflow from the second thickner.\n", +"F = 400;// [kg NaOH/h]\n", +"NF = mass_b/F;// [kg CaCO3/kg NaOH]\n", +"yF = 1;// [wt fraction NaOH in dry solid, CaCO3 free basis]\n", +"// Points R1, E3, RNpPlus1 and F are plotted as in Fig 13.30 (Pg 755) and locate the point deltaR at the intersection of lines FR1 and E3RNpPlus1 extended. The coordinates of point deltaR are NdeltaR=-0.1419, ydeltaR=-0.00213. Further computation must be done on enlarged section of the equilibrium diagram (Fig 13.31 (Pg 755)). Point deltaR is plotted and the stages stepped off in a usual manner. The construction are projected on the xy diagram. Three stages produce a value: y3=0.001\n", +"printf('The NaOH lost in sludge: %f%%\n',(lost_c/mass_c)*100);\n", +"printf('\n');\n", +"\n", +"// Solution (b)\n", +"//*** Data***//\n", +"lost_c = 0.001*mass_c;// [kg/h]\n", +"//***********//\n", +"\n", +"NNp_by_yNp = mass_b/lost_c;// [kg CaCO3/kg NaOH in final sludge]\n", +"// In order to determine the liquid content of the final sludge:\n", +"// Eqb=[N y_star]\n", +"Eqb = [0.659 0.01435;0.666 0.01015;0.677 0.002;0.679 0.001;0.680 0.0005];\n", +"N_by_ystar = zeros(5);\n", +"for i = 1:5\n", +" N_by_ystar(i) = Eqb(i,1)/(Eqb(i,2));\n", +"end\n", +"scf(43);\n", +"plot(Eqb(:,1),Eqb(:,2));\n", +"xgrid();\n", +"xlabel('x Wt fraction of NaOH');\n", +"ylabel('N kg CaCO3 / kg solution');\n", +"title('Equilibrium plot')\n", +"// By Interpolation, for N_by_ystar=NNp_by_yNp:\n", +"NNp = interpln([(N_by_ystar)';Eqb(:,1)'],NNp_by_yNp);// [kg CaCO3/kg soln]\n", +"yNp = NNp/NNp_by_yNp;// [wt fraction NaOH in the liquid of the final sludge]\n", +"ENp = mass_b/NNp;// [kg/h]\n", +"ENp_a = ENp-lost_c;// [kg/h]\n", +"overflow_c = mass_c-lost_c;// [kg/h]\n", +"R1 = overflow_c/0.1;// [kg/h]\n", +"R1_a = R1-overflow_c;// [kg/h]\n", +"RNpPlus1 = R1_a+sludge_a;// [kg/h]\n", +"// On the operating diagram (Fig 13.32 (Pg 757)) point deltaR is located and stages were constructed. \n", +"// Beyond the fourth stage, the ratio of the overflow to the liquid in the sludge become substantially constant.\n", +"R_by_E = RNpPlus1/ENp;\n", +"// This is the initial slope of the operating line on the lower part of the figure.\n", +"// From Illustration 13.2:\n", +"m = 0.01015/0.00450;\n", +"Value1 = R_by_E/m;\n", +"xNpPlus1 = 0;// [kg NaOH/kg solution]\n", +"y4 = 0.007;// [wt fraction NaOH in the liquid]\n", +"Value2 = (yNp-(m*xNpPlus1))/(y4-(m*xNpPlus1));\n", +"// From Fig 5.16: (Pg 129):\n", +"// An Additional 2.3 stages beyond 4 are computed graphically are required.\n", +"// An additional two stage will make yNp/y4=0.099:\n", +"yNp = 0.099*y4;// [wt fraction NaOH in the liquid]\n", +"printf('%f kg NaOH was lost if 6 thickners were used\n',yNp*ENp);\n", +"// An additional three stage will make yNp/y4=0.0365:\n", +"yNp = 0.0365*y4;// [wt fraction NaOH in the liquid]\n", +"printf('%f kg NaOH was lost if 7 thickners were used\n',yNp*ENp);" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 13.4: Multistage_Countercurrent_Leaching.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clear;\n", +"clc;\n", +"\n", +"// Illustration 13.4\n", +"// Page: 758\n", +"\n", +"printf('Illustration 13.4 - Page: 758\n\n');\n", +"\n", +"// Solution\n", +"\n", +"//***Data***//\n", +"// a:oil b:soyabean c:hexane\n", +"// Data=[100y*(Wt % oil in soln) 1/N(kg soln retained/kg insoluble solid)]\n", +"Data = [0 0.58;20 0.66;30 0.70];\n", +"// Soyabean feed:\n", +"percent_b = 20;// [soluble]\n", +"yF = 1;// [mass fraction oil,solid free basis]\n", +"// Solvent:\n", +"RNpPlus1 = 1;// [hexane,kg]\n", +"xNpPlus1 = 0;// [mass fraction oil]\n", +"// Leached Solids:\n", +"leached = 0.005;// [fraction of oil to be leached]\n", +"// Miscella:\n", +"percent_miscella = 10;// [percent of insoluble solid]\n", +"//**********//\n", +"\n", +"N = zeros(3);\n", +"ystar_By_N = zeros(3);\n", +"for i = 1:3\n", +" N(i) = 1/Data(i,2);// [kg insoluble solid/kg soln retained]\n", +" ystar_By_N(i) = Data(i,1)/(100*N(i));// [kg oil/kg insoluble solid]\n", +"end\n", +"// Basis: 1 kg flakes introduced\n", +"// Soyabean feed:\n", +"mass_b = 1-(percent_b/100);// [insoluble,kg]\n", +"F = 1-mass_b;// [kg]\n", +"NF = mass_b/F;// [kg insoluble solid/kg oil]\n", +"\n", +"// Leached Solids:\n", +"Ratio = leached/(1-leached);// [kg oil/kg insoluble solid]\n", +"// By interpolation:\n", +"Np = interpln([ystar_By_N';N'],Ratio);\n", +"miscella_b = (percent_miscella/100)*mass_b;// [Insoluble solid lost to miscella,kg]\n", +"leached_b = (1-(percent_miscella/100))*mass_b;// [Insoluble solid in miscella,kg]\n", +"ENp = leached_b/Np;// [kg soln retained]\n", +"retained_a = Ratio*leached_b;// [oil retained,kg]\n", +"retained_c = ENp-retained_a;// [Hexane retained,kg]\n", +"yNp = retained_a/ENp;// [mass fraction of oil in retained liquid]\n", +"\n", +"// Miscella:\n", +"mass_c = 1-retained_c;// [kg]\n", +"mass_a = F-retained_a;// [kg]\n", +"R1 = mass_c+mass_a;// [clear miscella,kg]\n", +"x1 = mass_a/R1;// [mass fraction of oil in the liquid]\n", +"NR1 = miscella_b/R1;// [kg insoluble solid/kg soln]\n", +"\n", +"// The operating diagram is shown in Fig 13.33 (Pg 759).\n", +"// Point R1 represents the cloudy miscella and is therefore is displaced from the axis of he graph at NR1. Point deltaR is located as usual and the stages determined with the N=0 axis for all the stages but the first.\n", +"printf('Between 4 and 5 stages are required\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 +} |