(************** Content-type: application/mathematica ************** CreatedBy='Mathematica 5.0' Mathematica-Compatible Notebook This notebook can be used with any Mathematica-compatible application, such as Mathematica, MathReader or Publicon. The data for the notebook starts with the line containing stars above. To get the notebook into a Mathematica-compatible application, do one of the following: * Save the data starting with the line of stars above into a file with a name ending in .nb, then open the file inside the application; * Copy the data starting with the line of stars above to the clipboard, then use the Paste menu command inside the application. Data for notebooks contains only printable 7-bit ASCII and can be sent directly in email or through ftp in text mode. Newlines can be CR, LF or CRLF (Unix, Macintosh or MS-DOS style). NOTE: If you modify the data for this notebook not in a Mathematica- compatible application, you must delete the line below containing the word CacheID, otherwise Mathematica-compatible applications may try to use invalid cache data. For more information on notebooks and Mathematica-compatible applications, contact Wolfram Research: web: http://www.wolfram.com email: info@wolfram.com phone: +1-217-398-0700 (U.S.) Notebook reader applications are available free of charge from Wolfram Research. *******************************************************************) (*CacheID: 232*) (*NotebookFileLineBreakTest NotebookFileLineBreakTest*) (*NotebookOptionsPosition[ 34078, 994]*) (*NotebookOutlinePosition[ 34744, 1017]*) (* CellTagsIndexPosition[ 34700, 1013]*) (*WindowFrame->Normal*) Notebook[{ Cell[CellGroupData[{ Cell[TextData[{ StyleBox["Computing of NonIsothermal Effectiveness Factor ", "Subsubtitle", FontSize->18, FontSlant->"Plain", FontVariations->{"CompatibilityType"->0}], StyleBox["\n", FontWeight->"Plain", FontVariations->{"CompatibilityType"->0}], StyleBox["Author's Data", FontSize->14, FontWeight->"Bold"], StyleBox[": ", FontSize->14], StyleBox["Housam BINOUS\nDepartment of Chemical Engineering\nNational \ Institute of Applied Sciences and Technology\nTunis, TUNISIA\nEmail: \ binoushousam@yahoo.com ", FontSize->14, FontWeight->"Plain"] }], "Title", TextAlignment->Center, Background->RGBColor[0.605478, 0.996109, 0.605478]], Cell[TextData[{ StyleBox["References:", FontWeight->"Bold"], "\nM. B. Cutlip and M. Shacham, Problem Solving in Chemical Engineering \ with Numerical Methods, Prentice Hall, 1999." }], "Subsubtitle", Background->RGBColor[0.996109, 0.753918, 0.996109]] }, Open ]], Cell[CellGroupData[{ Cell[TextData[StyleBox["Simultaneous diffusion and nonitothermal first order \ reaction in spherical catalyst pellets.", "Subsubtitle", FontVariations->{"CompatibilityType"->0}]], "Title", Background->RGBColor[0.765637, 0.996109, 0.996109]], Cell[CellGroupData[{ Cell["Integration of Heat Transfer Equation", "Subsubtitle", Background->RGBColor[0.996109, 0.996109, 0.68751]], Cell[CellGroupData[{ Cell[BoxData[ \(T[r] = Ts + \[CapitalDelta]H\_R\ De/ ke\ \((CA[r] - CAs)\) /. \((\[CapitalDelta]H\_R\ De/ ke \[Rule] \(-\[Beta]\)\ Ts/CAs)\)\)], "Input"], Cell[BoxData[ \(Ts - \(Ts\ \[Beta]\ \((\(-CAs\) + CA[r])\)\)\/CAs\)], "Output"] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell["Rate constant", "Subsubtitle", Background->RGBColor[0.996109, 0.996109, 0.68751]], Cell[CellGroupData[{ Cell[BoxData[ \(ka = ks\ a\ Exp[\(-\[Epsilon]\)\ \((Ts/T[r] - 1)\)] /. a\ ks \[Rule] \ \[Phi]^2\ De/Rs^2\)], "Input"], Cell[BoxData[ \(\(De\ \[ExponentialE]\^\(\(-\[Epsilon]\)\ \((\(-1\) + Ts\/\(Ts - \(Ts\ \ \[Beta]\ \((\(-CAs\) + CA[r])\)\)\/CAs\))\)\)\ \[Phi]\^2\)\/Rs\^2\)], "Output"] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell["Numerical values of the parameters of the problem", "Subsubtitle", Background->RGBColor[0.996109, 0.996109, 0.68751]], Cell[BoxData[ \(Rs = 0.01; \[Phi] = 1; \[Epsilon] = 30; CAs = 0.01; Ts = 400; \[CapitalDelta]H\_R = \(-8\)\ 10^7; \[Beta] = 0.2; De = 10^\(-6\);\)], "Input"] }, Closed]], Cell[CellGroupData[{ Cell["\<\ Split Boundary Problem solved using the Shooting Method. Must solve \ simultaneously for nonisothermal effectiveness factor.\ \>", "Subsubtitle", Background->RGBColor[0.996109, 0.996109, 0.68751]], Cell[BoxData[ \(myODEsoln[\[CapitalOmega]_] := NDSolve[{\(-De\)\ D[ r^2\ D[CA[r], {r, 1}], {r, 1}] \[Equal] \ \(-ka\)\ CA[ r]\ r^2, D[\[Eta][r], {r, 1}] \[Equal] \ 3\ Exp[\(-\[Epsilon]\)\ \((Ts/T[r] - 1)\)]\ CA[ r]\ r^2/\((CAs\ Rs^3)\), \[Eta][10^\(-5\)] \[Equal] \ 0, \(CA\ '\)[10^\(-5\)] \[Equal] \ 0, CA[10^\(-5\)] \[Equal] \ \[CapitalOmega]\ \ }, {CA[r], \[Eta][ r]}, {r, 10^\(-5\), Rs}, MaxSteps \[Rule] 1000]\)], "Input"], Cell[BoxData[{ \(yend[\[CapitalOmega]_?NumericQ] := Flatten[\((CA[r] /. myODEsoln[\[CapitalOmega]])\) /. r \[Rule] Rs]\), "\n", \(\(bc = \(FindRoot[ yend[\[CapitalOmega]] \[Equal] CAs, {\[CapitalOmega], 0, 0.000001}, MaxIterations \[Rule] \ 1000]\)[\([1, 2]\)];\)\)}], "Input"], Cell[CellGroupData[{ Cell[BoxData[ \(bc\)], "Input"], Cell[BoxData[ \(0.000023367246172997336`\)], "Output"] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell["Plot of concentration profile", "Subsubtitle", Background->RGBColor[0.996109, 0.996109, 0.68751]], Cell[CellGroupData[{ Cell[BoxData[ \(Plot[Evaluate[CA[r] /. myODEsoln[bc]], {r, 0.00001, Rs}]\)], "Input"], Cell[GraphicsData["PostScript", "\<\ %! %%Creator: Mathematica %%AspectRatio: .61803 MathPictureStart /Mabs { Mgmatrix idtransform Mtmatrix dtransform } bind def /Mabsadd { Mabs 3 -1 roll add 3 1 roll add exch } bind def %% Graphics %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10 scalefont setfont % Scaling calculations 0.0238095 95.2381 0.0147151 58.8604 [ [.21429 .00222 -15 -9 ] [.21429 .00222 15 0 ] [.40476 .00222 -15 -9 ] [.40476 .00222 15 0 ] [.59524 .00222 -15 -9 ] [.59524 .00222 15 0 ] [.78571 .00222 -15 -9 ] [.78571 .00222 15 0 ] [.97619 .00222 -12 -9 ] [.97619 .00222 12 0 ] [.01131 .13244 -30 -4.5 ] [.01131 .13244 0 4.5 ] [.01131 .25016 -30 -4.5 ] [.01131 .25016 0 4.5 ] [.01131 .36788 -30 -4.5 ] [.01131 .36788 0 4.5 ] [.01131 .4856 -30 -4.5 ] [.01131 .4856 0 4.5 ] [.01131 .60332 -24 -4.5 ] [.01131 .60332 0 4.5 ] [ 0 0 0 0 ] [ 1 .61803 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath 0 g .25 Mabswid [ ] 0 setdash .21429 .01472 m .21429 .02097 L s [(0.002)] .21429 .00222 0 1 Mshowa .40476 .01472 m .40476 .02097 L s [(0.004)] .40476 .00222 0 1 Mshowa .59524 .01472 m .59524 .02097 L s [(0.006)] .59524 .00222 0 1 Mshowa .78571 .01472 m .78571 .02097 L s [(0.008)] .78571 .00222 0 1 Mshowa .97619 .01472 m .97619 .02097 L s [(0.01)] .97619 .00222 0 1 Mshowa .125 Mabswid .07143 .01472 m .07143 .01847 L s .11905 .01472 m .11905 .01847 L s .16667 .01472 m .16667 .01847 L s .2619 .01472 m .2619 .01847 L s .30952 .01472 m .30952 .01847 L s .35714 .01472 m .35714 .01847 L s .45238 .01472 m .45238 .01847 L s .5 .01472 m .5 .01847 L s .54762 .01472 m .54762 .01847 L s .64286 .01472 m .64286 .01847 L s .69048 .01472 m .69048 .01847 L s .7381 .01472 m .7381 .01847 L s .83333 .01472 m .83333 .01847 L s .88095 .01472 m .88095 .01847 L s .92857 .01472 m .92857 .01847 L s .25 Mabswid 0 .01472 m 1 .01472 L s .02381 .13244 m .03006 .13244 L s [(0.002)] .01131 .13244 1 0 Mshowa .02381 .25016 m .03006 .25016 L s [(0.004)] .01131 .25016 1 0 Mshowa .02381 .36788 m .03006 .36788 L s [(0.006)] .01131 .36788 1 0 Mshowa .02381 .4856 m .03006 .4856 L s [(0.008)] .01131 .4856 1 0 Mshowa .02381 .60332 m .03006 .60332 L s [(0.01)] .01131 .60332 1 0 Mshowa .125 Mabswid .02381 .04415 m .02756 .04415 L s .02381 .07358 m .02756 .07358 L s .02381 .10301 m .02756 .10301 L s .02381 .16187 m .02756 .16187 L s .02381 .1913 m .02756 .1913 L s .02381 .22073 m .02756 .22073 L s .02381 .27959 m .02756 .27959 L s .02381 .30902 m .02756 .30902 L s .02381 .33845 m .02756 .33845 L s .02381 .39731 m .02756 .39731 L s .02381 .42674 m .02756 .42674 L s .02381 .45617 m .02756 .45617 L s .02381 .51503 m .02756 .51503 L s .02381 .54446 m .02756 .54446 L s .02381 .57389 m .02756 .57389 L s .25 Mabswid .02381 0 m .02381 .61803 L s 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath .5 Mabswid .02476 .01609 m .02594 .01609 L .027 .01609 L .02824 .01609 L .0294 .01609 L .03148 .01609 L .03373 .01609 L .03884 .0161 L .04355 .0161 L .05317 .01612 L .06336 .01615 L .0737 .01618 L .08339 .01623 L .10545 .01635 L .12568 .01651 L .14418 .01669 L .16476 .01695 L .18385 .01726 L .22201 .01808 L .24251 .01868 L .2611 .01934 L .28195 .02026 L .30113 .02129 L .33964 .02407 L .35844 .02587 L .37909 .02829 L .40021 .03133 L .41947 .03472 L .45834 .04378 L .4774 .04959 L .49814 .05714 L .53889 .07647 L .57811 .1018 L .61583 .13321 L .65599 .17439 L .69464 .22076 L .73574 .2755 L .77532 .33135 L .81339 .38601 L .85391 .44355 L .89292 .4972 L .93437 .55158 L .97432 .60107 L .97619 .60332 L s % End of Graphics MathPictureEnd \ \>"], "Graphics", ImageSize->{288, 177.938}, ImageMargins->{{43, 0}, {0, 0}}, ImageRegion->{{0, 1}, {0, 1}}, ImageCache->GraphicsData["Bitmap", "\<\ CF5dJ6E]HGAYHf4PAg9QL6QYHgOol0025oo`03001oogoo06io o`@008Yoo`008Goo00<007ooOol0LWoo0P00R7oo000QOol00`00Oomoo`1dOol20026Ool0025oo`03 001oogoo07Ioo`8008Aoo`008Goo00<007ooOol0N7oo0P00PWoo000QOol00`00Oomoo`1jOol20020 Ool0025oo`8007eoo`8007ioo`008Goo00<007ooOol0OWoo0P00O7oo000QOol00`00Oomoo`20Ool2 001jOol0025oo`03001oogoo089oo`8007Qoo`008Goo00<007ooOol0Q7oo00<007ooOol0MGoo000Q Ool00`00Oomoo`25Ool2001eOol0025oo`03001oogoo08Moo`8007=oo`008Goo00<007ooOol0RGoo 00<007ooOol0L7oo000QOol2002;Ool00`00Oomoo`1_Ool0025oo`03001oogoo08]oo`8006moo`00 8Goo00<007ooOol0SGoo00<007ooOol0K7oo000QOol00`00Oomoo`2>Ool00`00Oomoo`1[Ool0025o o`03001oogoo08moo`03001oogoo06Yoo`000Woo0P0017oo0P0017oo0P0017oo0P000goo100017oo 00<007ooOol0T7oo0P00JWoo00001Goo001oogoo000027oo00@007ooOol0009oo`04001oogoo0002 Ool00`00Oomoo`05Ool00`00Oomoo`2BOol00`00Oomoo`1WOol00005Ool007ooOol00008Ool01000 Oomoo`000Woo00@007ooOol000=oo`03001oogoo00Aoo`8009Aoo`03001oogoo06Ioo`0000Eoo`00 Oomoo`0000Qoo`04001oogoo0002Ool01000Oomoo`0017oo00<007ooOol00goo00<007ooOol0U7oo 00<007ooOol0IGoo00001Goo001oogoo000027oo00@007ooOol0009oo`04001oogoo0002Ool01000 Oomoo`0017oo00<007ooOol0UGoo00<007ooOol0I7oo0002Ool2000:Ool20004Ool20004Ool20005 Ool00`00Oomoo`2FOol00`00Oomoo`1SOol0025oo`03001oogoo09Moo`03001oogoo069oo`008Goo 00<007ooOol0V7oo00<007ooOol0HGoo000QOol00`00Oomoo`2HOol00`00Oomoo`1QOol0025oo`03 001oogoo09Uoo`03001oogoo061oo`008Goo0P00Vgoo00<007ooOol0Ggoo000QOol00`00Oomoo`2K Ool00`00Oomoo`1NOol0025oo`03001oogoo09aoo`03001oogoo05eoo`008Goo00<007ooOol0WGoo 00<007ooOol0G7oo000QOol00`00Oomoo`2NOol00`00Oomoo`1KOol0025oo`03001oogoo09moo`03 001oogoo05Yoo`008Goo00<007ooOol0X7oo00<007ooOol0FGoo000QOol2002ROol00`00Oomoo`1H Ool0025oo`03001oogoo0:9oo`03001oogoo05Moo`008Goo00<007ooOol0XWoo00<007ooOol0Egoo 000QOol00`00Oomoo`2SOol00`00Oomoo`1FOol0025oo`03001oogoo0:Aoo`03001oogoo05Eoo`00 8Goo00<007ooOol0YGoo00<007ooOol0E7oo000QOol00`00Oomoo`2VOol00`00Oomoo`1COol0025o o`800:Qoo`03001oogoo059oo`008Goo00<007ooOol0Z7oo00<007ooOol0DGoo000QOol00`00Oomo o`2XOol00`00Oomoo`1AOol0025oo`03001oogoo0:Uoo`03001oogoo051oo`008Goo00<007ooOol0 ZWoo00<007ooOol0Cgoo000QOol00`00Oomoo`2[Ool00`00Oomoo`1>Ool0009oo`8000Aoo`8000Ao o`8000Aoo`8000Eoo`<000=oo`03001oogoo0:]oo`03001oogoo04ioo`0000Eoo`00Oomoo`0000Qo o`04001oogoo0002Ool01000Oomoo`001Goo00<007ooOol00Woo00<007ooOol0[7oo00<007ooOol0 CGoo00001Goo001oogoo000027oo00@007ooOol0009oo`04001oogoo0002Ool50003Ool2002^Ool0 0`00Oomoo`1Goo0002Ool20004Ool20004Ool20004Ool20003Ool30005Ool00`00Oomoo`30Ool00`00Oomoo`0i Ool00005Ool007ooOol00008Ool01000Oomoo`000Woo00@007ooOol0009oo`04001oogoo0004Ool0 0`00Oomoo`31Ool00`00Oomoo`0hOol00005Ool007ooOol00008Ool01000Oomoo`000Woo00@007oo Ool0009oo`04001oogoo0004Ool20033Ool00`00Oomoo`0gOol00005Ool007ooOol00008Ool01000 Oomoo`000Woo00@007ooOol0009oo`<000Eoo`03001oogoo0<=oo`03001oogoo03Ioo`0000Eoo`00 Oomoo`0000Qoo`04001oogoo0002Ool01000Oomoo`000goo00<007ooOol017oo00<007ooOol0`goo 00<007ooOol0=Woo0002Ool2000:Ool20004Ool20004Ool30004Ool00`00Oomoo`34Ool00`00Oomo o`0eOol0025oo`03001oogoo0Ool00`00 Oomoo`0[Ool0025oo`03001oogoo05oo`03 001oogoo01Qoo`008Goo00<007ooOol0hWoo00<007ooOol05goo000QOol2003TOol00`00Oomoo`0F Ool0025oo`03001oogoo0>Aoo`03001oogoo01Eoo`008Goo00<007ooOol0iGoo00<007ooOol057oo 000QOol00`00Oomoo`3VOol00`00Oomoo`0COol0025oo`03001oogoo0>Ioo`03001oogoo01=oo`00 8Goo00<007ooOol0igoo00<007ooOol04Woo000QOol00`00Oomoo`3XOol00`00Oomoo`0AOol0025o o`800>Yoo`03001oogoo011oo`008Goo00<007ooOol0jGoo00<007ooOol047oo000QOol00`00Oomo o`3ZOol00`00Oomoo`0?Ool0025oo`03001oogoo0>]oo`03001oogoo00ioo`008Goo00<007ooOol0 k7oo00<007ooOol03Goo000QOol00`00Oomoo`3/Ool00`00Oomoo`0=Ool000Qoo`8000Aoo`8000Ao o`8000=oo`D000=oo`03001oogoo0>eoo`03001oogoo00aoo`001goo00@007ooOol000Qoo`04001o ogoo0004Ool00`00Oomoo`03Ool00`00Oomoo`3^Ool00`00Oomoo`0;Ool000Moo`04001oogoo0008 Ool01000Oomoo`0017oo00<007ooOol00goo0P00l7oo00<007ooOol02Woo0007Ool01000Oomoo`00 27oo00@007ooOol000Aoo`03001oogoo00=oo`03001oogoo0?aoo`001goo00@007ooOol000Qoo`04 001oogoo0004Ool00`00Oomoo`03Ool00`00Oomoo`3lOol000Qoo`8000Yoo`8000Aoo`8000Eoo`03 001oogoo0?aoo`00ogoo8Goo003oOolQOol00?moob5oo`00ogoo8Goo003oOolQOol00?moob5oo`00 \ \>"], ImageRangeCache->{{{0, 287}, {176.938, 0}} -> {-0.00137187, -0.00117609, \ 4.13524*^-05, 6.69096*^-05}}], Cell[BoxData[ TagBox[\(\[SkeletonIndicator] Graphics \[SkeletonIndicator]\), False, Editable->False]], "Output"] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell["Plot of temperature profile", "Subsubtitle", Background->RGBColor[0.996109, 0.996109, 0.68751]], Cell[CellGroupData[{ Cell[BoxData[ \(Plot[Evaluate[T[r] /. myODEsoln[bc]], {r, 0.00001, Rs}]\)], "Input"], Cell[GraphicsData["PostScript", "\<\ %! %%Creator: Mathematica %%AspectRatio: .61803 MathPictureStart /Mabs { Mgmatrix idtransform Mtmatrix dtransform } bind def /Mabsadd { Mabs 3 -1 roll add 3 1 roll add exch } bind def %% Graphics %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10 scalefont setfont % Scaling calculations 0.0238095 95.2381 -2.9352 0.00737478 [ [.21429 .00222 -15 -9 ] [.21429 .00222 15 0 ] [.40476 .00222 -15 -9 ] [.40476 .00222 15 0 ] [.59524 .00222 -15 -9 ] [.59524 .00222 15 0 ] [.78571 .00222 -15 -9 ] [.78571 .00222 15 0 ] [.97619 .00222 -12 -9 ] [.97619 .00222 12 0 ] [.01131 .16221 -18 -4.5 ] [.01131 .16221 0 4.5 ] [.01131 .30971 -18 -4.5 ] [.01131 .30971 0 4.5 ] [.01131 .4572 -18 -4.5 ] [.01131 .4572 0 4.5 ] [.01131 .6047 -18 -4.5 ] [.01131 .6047 0 4.5 ] [ 0 0 0 0 ] [ 1 .61803 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath 0 g .25 Mabswid [ ] 0 setdash .21429 .01472 m .21429 .02097 L s [(0.002)] .21429 .00222 0 1 Mshowa .40476 .01472 m .40476 .02097 L s [(0.004)] .40476 .00222 0 1 Mshowa .59524 .01472 m .59524 .02097 L s [(0.006)] .59524 .00222 0 1 Mshowa .78571 .01472 m .78571 .02097 L s [(0.008)] .78571 .00222 0 1 Mshowa .97619 .01472 m .97619 .02097 L s [(0.01)] .97619 .00222 0 1 Mshowa .125 Mabswid .07143 .01472 m .07143 .01847 L s .11905 .01472 m .11905 .01847 L s .16667 .01472 m .16667 .01847 L s .2619 .01472 m .2619 .01847 L s .30952 .01472 m .30952 .01847 L s .35714 .01472 m .35714 .01847 L s .45238 .01472 m .45238 .01847 L s .5 .01472 m .5 .01847 L s .54762 .01472 m .54762 .01847 L s .64286 .01472 m .64286 .01847 L s .69048 .01472 m .69048 .01847 L s .7381 .01472 m .7381 .01847 L s .83333 .01472 m .83333 .01847 L s .88095 .01472 m .88095 .01847 L s .92857 .01472 m .92857 .01847 L s .25 Mabswid 0 .01472 m 1 .01472 L s .02381 .16221 m .03006 .16221 L s [(420)] .01131 .16221 1 0 Mshowa .02381 .30971 m .03006 .30971 L s [(440)] .01131 .30971 1 0 Mshowa .02381 .4572 m .03006 .4572 L s [(460)] .01131 .4572 1 0 Mshowa .02381 .6047 m .03006 .6047 L s [(480)] .01131 .6047 1 0 Mshowa .125 Mabswid .02381 .05159 m .02756 .05159 L s .02381 .08846 m .02756 .08846 L s .02381 .12534 m .02756 .12534 L s .02381 .19908 m .02756 .19908 L s .02381 .23596 m .02756 .23596 L s .02381 .27283 m .02756 .27283 L s .02381 .34658 m .02756 .34658 L s .02381 .38345 m .02756 .38345 L s .02381 .42033 m .02756 .42033 L s .02381 .49408 m .02756 .49408 L s .02381 .53095 m .02756 .53095 L s .02381 .56782 m .02756 .56782 L s .25 Mabswid .02381 0 m .02381 .61803 L s 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath .5 Mabswid .02476 .60332 m .02594 .60332 L .027 .60332 L .02824 .60332 L .0294 .60332 L .03148 .60332 L .03373 .60332 L .03884 .60331 L .04355 .6033 L .05317 .60329 L .06336 .60326 L .0737 .60322 L .08339 .60318 L .10545 .60306 L .12568 .6029 L .14418 .60271 L .16476 .60245 L .18385 .60215 L .22201 .60132 L .24251 .60072 L .2611 .60006 L .28195 .59914 L .30113 .59811 L .33964 .59532 L .35844 .59351 L .37909 .59109 L .40021 .58804 L .41947 .58465 L .45834 .57556 L .4774 .56974 L .49814 .56217 L .53889 .5428 L .57811 .51741 L .61583 .48593 L .65599 .44464 L .69464 .39817 L .73574 .3433 L .77532 .28733 L .81339 .23253 L .85391 .17485 L .89292 .12108 L .93437 .06657 L .97432 .01697 L .97619 .01472 L s % End of Graphics MathPictureEnd \ \>"], "Graphics", ImageSize->{288, 177.938}, ImageMargins->{{43, 0}, {0, 0}}, ImageRegion->{{0, 1}, {0, 1}}, ImageCache->GraphicsData["Bitmap", "\<\ CF5dJ6E]HGAYHf4PAg9QL6QYHgWoo00@007ooOol000Qoo`04001oogoo0002 Ool01000Oomoo`000Woo00<007ooOol05goo00@007ooOol000Qoo`04001oogoo0002Ool01000Oomo o`001Goo00<007ooOol05Goo00@007ooOol000Qoo`04001oogoo0002Ool01000Oomoo`000Woo00@0 07ooOol001Moo`04001oogoo0008Ool01000Oomoo`000Woo00@007ooOol0009oo`04001oogoo000J Ool01000Oomoo`0027oo00@007ooOol000Aoo`03001oogoo005oo`00>Woo00@007ooOol000Qoo`04 001oogoo0002Ool01000Oomoo`000goo00<007ooOol05Woo00@007ooOol000Qoo`04001oogoo0002 Ool01000Oomoo`000Woo1@005Woo00@007ooOol000Qoo`04001oogoo0002Ool01000Oomoo`000Woo 00@007ooOol001Moo`04001oogoo0008Ool01000Oomoo`000Woo00@007ooOol0009oo`04001oogoo 000JOol01000Oomoo`0027oo00@007ooOol000Aoo`03001oogoo005oo`00>Woo00@007ooOol000Qo o`04001oogoo0002Ool01000Oomoo`0017oo00<007ooOol05Goo00@007ooOol000Qoo`04001oogoo 0002Ool01000Oomoo`000Woo00@007ooOol001Moo`04001oogoo0008Ool01000Oomoo`000Woo00@0 07ooOol0009oo`<001Qoo`04001oogoo0008Ool01000Oomoo`000Woo00@007ooOol000=oo`8001]o o`04001oogoo0008Ool01000Oomoo`0017oo00<007ooOol00Goo000jOol01000Oomoo`0027oo00@0 07ooOol0009oo`04001oogoo0002Ool01000Oomoo`005Woo00@007ooOol000Qoo`04001oogoo0002 Ool01000Oomoo`000goo00<007oo00005goo00@007ooOol000Qoo`04001oogoo0002Ool01000Oomo o`000goo00<007ooOol05goo00@007ooOol000Qoo`04001oogoo0002Ool01000Oomoo`000Woo00@0 07ooOol001Yoo`04001oogoo0008Ool01000Oomoo`0017oo00<007ooOol00Goo000kOol2000:Ool2 0004Ool20004Ool2000HOol2000:Ool20004Ool20005Ool2000HOol2000:Ool20004Ool20004Ool3 000HOol2000:Ool20004Ool20004Ool2000LOol2000:Ool20004Ool20003Ool00?moob5oo`005Goo 00<007ooOol0ogoo2Goo000EOol00`00Oomoo`3oOol9Ool001Eoo`03001oogoo0?moo`Uoo`005Goo 00<007ooOol0ogoo2Goo000>Oooo000=0006Ool001Eoo`03001oogoo00Yoo`03001oogoo00Uoo`03 001oogoo00Yoo`03001oogoo00Yoo`03001oogoo00Uoo`03001oogoo00Yoo`03001oogoo00Yoo`03 001oogoo00Uoo`03001oogoo00Yoo`03001oogoo00Yoo`03001oogoo00Uoo`03001oogoo00Yoo`03 001oogoo00Yoo`03001oogoo00Uoo`03001oogoo00Yoo`03001oogoo00Yoo`03001oogoo00Uoo`03 001oogoo00Yoo`03001oogoo00Yoo`03001oogoo00Uoo`8000aoo`005Goo00<007ooOol0<7oo00<0 07ooOol0;goo00<007ooOol0<7oo00<007ooOol0<7oo00<007ooOol0;Woo00<007oo000037oo000E Ool00`00Oomoo`3hOol00`00Oomoo`0=Ool001Eoo`03001oogoo0?Moo`03001oogoo00ioo`005Goo 00<007ooOol0mgoo00<007ooOol03Woo000EOol00`00Oomoo`3fOol00`00Oomoo`0?Ool001Eoo`03 001oogoo0?Eoo`03001oogoo011oo`005Goo00<007ooOol0m7oo00<007ooOol04Goo000EOol00`00 Oomoo`3cOol00`00Oomoo`0BOol001Eoo`800?=oo`03001oogoo01=oo`005Goo00<007ooOol0lWoo 00<007ooOol04goo000EOol00`00Oomoo`3aOol00`00Oomoo`0DOol001Eoo`03001oogoo0?1oo`03 001oogoo01Eoo`005Goo00<007ooOol0kgoo00<007ooOol05Woo000EOol00`00Oomoo`3^Ool00`00 Oomoo`0GOol001Eoo`03001oogoo0>ioo`03001oogoo01Moo`005Goo00<007ooOol0kGoo00<007oo Ool067oo000EOol00`00Oomoo`3/Ool00`00Oomoo`0IOol001Eoo`03001oogoo0>]oo`03001oogoo 01Yoo`005Goo0P00k7oo00<007ooOol06Woo000EOol00`00Oomoo`3ZOol00`00Oomoo`0KOol001Eo o`03001oogoo0>Uoo`03001oogoo01aoo`005Goo00<007ooOol0j7oo00<007ooOol07Goo000EOol0 0`00Oomoo`3XOol00`00Oomoo`0MOol001Eoo`03001oogoo0>Moo`03001oogoo01ioo`005Goo00<0 07ooOol0iWoo00<007ooOol07goo000EOol00`00Oomoo`3UOol00`00Oomoo`0POol001Eoo`03001o ogoo0>Eoo`03001oogoo021oo`005Goo00<007ooOol0i7oo00<007ooOol08Goo000EOol2003TOol0 0`00Oomoo`0ROol001Eoo`03001oogoo0>=oo`03001oogoo029oo`005Goo00<007ooOol0hWoo00<0 07ooOol08goo000EOol00`00Oomoo`3QOol00`00Oomoo`0TOol001Eoo`03001oogoo0>1oo`03001o ogoo02Eoo`005Goo00<007ooOol0h7oo00<007ooOol09Goo000EOol00`00Oomoo`3OOol00`00Oomo o`0VOol001Eoo`03001oogoo0=ioo`03001oogoo02Moo`000goo0`0000=oo`0000000P000goo0P00 1Goo00<007ooOol0gWoo00<007ooOol09goo0004Ool01000Oomoo`001Goo00@007ooOol000Aoo`03 001oogoo0=eoo`03001oogoo02Qoo`0000=oo`0000000`000Woo00<007ooOol00Woo00@007ooOol0 00Aoo`800=eoo`03001oogoo02Uoo`0000Eoo`00Oomoo`0000Aoo`05001oogooOol00002Ool00`00 Oomoo`02Ool00`00Oomoo`3KOol00`00Oomoo`0ZOol0009oo`03001oo`00009oo`04001oogoo0002 Ool01000Oomoo`0017oo00<007ooOol0fgoo00<007ooOol0:Woo0003Ool20003Ool20004Ool20005 Ool00`00Oomoo`3JOol00`00Oomoo`0[Ool001Eoo`03001oogoo0=Uoo`03001oogoo02aoo`005Goo 00<007ooOol0fGoo00<007ooOol0;7oo000EOol00`00Oomoo`3HOol00`00Oomoo`0]Ool001Eoo`03 001oogoo0=Moo`03001oogoo02ioo`005Goo00<007ooOol0eWoo00<007ooOol0;goo000EOol2003G Ool00`00Oomoo`0_Ool001Eoo`03001oogoo0=Eoo`03001oogoo031oo`005Goo00<007ooOol0e7oo 00<007ooOol0Goo000EOol0 0`00Oomoo`3;Ool00`00Oomoo`0jOol001Eoo`03001oogoo0goo000EOol00`00Oomoo`39Ool00`00Oomoo`0lOol001Eoo`03001o ogoo0Ool2001hOol001Eoo`03001oogoo 08eoo`03001oogoo07Qoo`005Goo00<007ooOol0S7oo00<007ooOol0NGoo000EOol2002;Ool2001l Ool001Eoo`03001oogoo08Uoo`03001oogoo07aoo`005Goo00<007ooOol0Qgoo0P00Ogoo000EOol0 0`00Oomoo`25Ool20021Ool001Eoo`03001oogoo08=oo`8008=oo`005Goo00<007ooOol0PGoo0P00 QGoo000EOol00`00Oomoo`1oOol20027Ool001Eoo`03001oogoo07eoo`8008Uoo`005Goo00<007oo Ool0NWoo0`00Rgoo000EOol00`00Oomoo`1hOol2002>Ool001Eoo`8007Ioo`<0091oo`005Goo00<0 07ooOol0Lgoo0P00Tgoo000EOol00`00Oomoo`1_Ool4002EOol001Eoo`03001oogoo06Uoo`H009Uo o`005Goo00<007ooOol0I7oo1@00Wgoo000EOol00`00Oomoo`1OOol5002TOol001Eoo`03001oogoo 05Uoo`H00:Uoo`005Goo00<007ooOol0CGoo3000[goo0003Ool30002Ool20004Ool20005Ool00`00 Oomoo`0eOolH002kOol000Aoo`04001oogoo0002Ool01000Oomoo`000Woo00<007ooOol00Woo>000 dgoo00000goo0000000300001Goo001oogoo00000Woo00@007ooOol000Aoo`800?moo`Yoo`0000Eo o`00Oomoo`0000=oo`8000=oo`04001oogoo0004Ool00`00Oomoo`3oOol9Ool0009oo`03001oo`00 009oo`04001oogoo0002Ool01000Oomoo`0017oo00<007ooOol0ogoo2Goo0003Ool20003Ool20004 Ool20005Ool00`00Oomoo`3oOol9Ool00?moob5oo`00ogoo8Goo0000\ \>"], ImageRangeCache->{{{0, 287}, {176.938, 0}} -> {-0.000840883, 392.969, \ 3.94216*^-05, 0.509091}}], Cell[BoxData[ TagBox[\(\[SkeletonIndicator] Graphics \[SkeletonIndicator]\), False, Editable->False]], "Output"] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell["Checking Boundary conditions", "Subsubtitle", Background->RGBColor[0.996109, 0.996109, 0.68751]], Cell[CellGroupData[{ Cell[BoxData[ \(\(T[r] /. myODEsoln[bc]\) /. r \[Rule] Rs\)], "Input"], Cell[BoxData[ \({400.0000012317936`}\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(\(CA[r] /. myODEsoln[bc]\) /. r \[Rule] Rs\)], "Input"], Cell[BoxData[ \({0.009999999846025804`}\)], "Output"] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell["\<\ Temperature inside catalyst pellet is 479.813 K. This temperature is higher \ than Ts=400 K (temperature at surface of the catalyst pellet). This is due to \ the fact that reaction is exothermic (\[Beta] =0.2>0). Thus, we expect \ \[Eta], the nonisothermal effectiveness factor, to be higher than 1. For \ \[Phi] = 1 and \[Epsilon] = 30, we find \[Eta] = 5.81 which agrees with \ figure 7-10 of book by M.B.Cutlip and M.Shacham,Problem Solving in Chemical \ Engineering with Numerical Methods,Prentice Hall,1999.\ \>", "Subsubtitle", Background->RGBColor[0.996109, 0.996109, 0.68751]], Cell[CellGroupData[{ Cell[BoxData[ \(\(T[r] /. myODEsoln[bc]\) /. r \[Rule] 10^\(-5\)\)], "Input"], Cell[BoxData[ \({479.81306203061604`}\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(\(\[Eta][r] /. myODEsoln[bc]\) /. r \[Rule] Rs\)], "Input"], Cell[BoxData[ \({5.81672303848871`}\)], "Output"] }, Open ]] }, Closed]] }, Open ]] }, FrontEndVersion->"5.0 for Microsoft Windows", ScreenRectangle->{{0, 1024}, {0, 695}}, CellGrouping->Manual, WindowSize->{622, 671}, WindowMargins->{{16, Automatic}, {Automatic, 1}} ] (******************************************************************* Cached data follows. If you edit this Notebook file directly, not using Mathematica, you must remove the line containing CacheID at the top of the file. The cache data will then be recreated when you save this file from within Mathematica. *******************************************************************) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[CellGroupData[{ Cell[1776, 53, 684, 20, 229, "Title"], Cell[2463, 75, 261, 6, 87, "Subsubtitle"] }, Open ]], Cell[CellGroupData[{ Cell[2761, 86, 244, 3, 118, "Title"], Cell[CellGroupData[{ Cell[3030, 93, 113, 1, 46, "Subsubtitle"], Cell[CellGroupData[{ Cell[3168, 98, 194, 4, 30, "Input"], Cell[3365, 104, 83, 1, 42, "Output"] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell[3497, 111, 89, 1, 48, "Subsubtitle"], Cell[CellGroupData[{ Cell[3611, 116, 137, 3, 30, "Input"], Cell[3751, 121, 172, 2, 71, "Output"] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell[3972, 129, 125, 1, 48, "Subsubtitle"], Cell[4100, 132, 171, 3, 50, "Input"] }, Closed]], Cell[CellGroupData[{ Cell[4308, 140, 208, 4, 68, "Subsubtitle"], Cell[4519, 146, 546, 10, 110, "Input"], Cell[5068, 158, 344, 7, 70, "Input"], Cell[CellGroupData[{ Cell[5437, 169, 35, 1, 30, "Input"], Cell[5475, 172, 58, 1, 29, "Output"] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell[5582, 179, 105, 1, 48, "Subsubtitle"], Cell[CellGroupData[{ Cell[5712, 184, 89, 1, 30, "Input"], Cell[5804, 187, 13411, 371, 186, 3599, 245, "GraphicsData", "PostScript", \ "Graphics"], Cell[19218, 560, 130, 3, 29, "Output"] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell[19397, 569, 103, 1, 48, "Subsubtitle"], Cell[CellGroupData[{ Cell[19525, 574, 88, 1, 30, "Input"], Cell[19616, 577, 12827, 351, 186, 3371, 230, "GraphicsData", "PostScript", \ "Graphics"], Cell[32446, 930, 130, 3, 29, "Output"] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell[32625, 939, 104, 1, 48, "Subsubtitle"], Cell[CellGroupData[{ Cell[32754, 944, 74, 1, 30, "Input"], Cell[32831, 947, 54, 1, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[32922, 953, 75, 1, 30, "Input"], Cell[33000, 956, 57, 1, 29, "Output"] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell[33106, 963, 596, 9, 168, "Subsubtitle"], Cell[CellGroupData[{ Cell[33727, 976, 81, 1, 30, "Input"], Cell[33811, 979, 55, 1, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[33903, 985, 79, 1, 30, "Input"], Cell[33985, 988, 53, 1, 29, "Output"] }, Open ]] }, Closed]] }, Open ]] } ] *) (******************************************************************* End of Mathematica Notebook file. *******************************************************************)