(*********************************************************************** Mathematica-Compatible Notebook This notebook can be used on any computer system with Mathematica 4.0, MathReader 4.0, or any compatible application. 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[ 178672, 5231]*) (*NotebookOutlinePosition[ 195529, 5833]*) (* CellTagsIndexPosition[ 195097, 5816]*) (*WindowFrame->Normal*) Notebook[{ Cell["\<\ Solution of the natural convection boundary-layer flow near a \ heated flat plate\ \>", "Title"], Cell[TextData[{ "This notebook has been written in ", StyleBox["Mathematica ", FontSlant->"Italic"], "by \n\n", StyleBox["Mark J. McCready\nProfessor and Chair of Chemical Engineering\n\ University of Notre Dame\nNotre Dame IN 46556\nUSA", FontSize->14], "\n\nmjm@nd.edu\nhttp://www.nd.edu/~mjm/\n\n\nIt is copyrighted to the \ extent allowed by whatever laws pertain to the World Wide Web and the \ Internet.\n\nI would hope that as a professional courtesy, that if you use \ it, that this notice remain visible to other users. \nThere is no charge for \ copying and dissemination \n\nVersion: 5/25/00\nMore recent versions of \ this notebook may be available at the web site:\n", ButtonBox["http://www.nd.edu/~mjm/thermal_boundarylayer.nb", ButtonData:>{ URL[ "http://www.nd.edu/~mjm/thermal_boundarylayer.nb"], None}, ButtonStyle->"Hyperlink"] }], "Text"], Cell["\<\ Reference: F. P. Incropera and D. P. DeWitt, Fundamentals of Heat \ and Mass Transfer, 4th ed. Wiley, 1996. \ \>", "Subsubsection"], Cell[CellGroupData[{ Cell["Summary", "Subtitle"], Cell[TextData[{ "This notebook is intended to give a first introduction to natural \ convection heat transfer from a solid body. The application of interest is \ heat transfer from a heated sphere (a light bulb) that is enclosed in a clear \ plastic box. The physical experiment is done at Notre Dame in the \ undergraduate chemical engineering lab (", ButtonBox["http://www.nd.edu/~chegdept/classes/cheg358.html", ButtonData:>{ URL[ "http://www.nd.edu/~chegdept/classes/cheg358.html"], None}, ButtonStyle->"Hyperlink"], "). The surface temperature of the light bulb is measured as a function of \ position around the light bulb. These data show the development of a thermal \ boundary layer.\n\nThis notebook starts with the boundary layer equations for \ a heated, vertical flat plate under conditions of high Grashof number. Using \ a similarity variable that involves the local Grashof number, the PDE's are \ reduced to ODE's. The ODE's are solved for different Prandtl numbers to show \ the velocity and temperature profiles. " }], "Text", Evaluatable->False, AspectRatioFixed->True], Cell[CellGroupData[{ Cell[TextData[{ StyleBox["Mathematica", FontSlant->"Italic"], " aside" }], "Subsection"], Cell[TextData[{ "In ", StyleBox["Mathematica", FontSlant->"Italic"], ", it is convenient to give all expressions a \"name\". I try to pick ones \ that are consistent with what is being done (but sometimes \"temp\" is used). \ This assignment is done with an \"=\" sign. To make an equation, a \"==\" \ is used. This distinction is very useful in computer algebra and is employed \ in all of the packages with which I am familiar ." }], "Text"], Cell[CellGroupData[{ Cell["Input notation", "Subsubsection"], Cell["\<\ I would enter the dynamic boundary condition from a key pad \ as\ \>", "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(dc = \(-\ \[Gamma]\)\ D[\[Eta][x, t], {x, 2}]\ - \n\t\t\(\[Rho]\_1\) \((\(-D[\[Phi]\_1[x, y, t], t]\) - U1\ D[\[Phi]\_1[x, y, t], x] - g\ \[Eta][x, t])\) + \n\t\t\t\(\[Rho]\_2\) \((\(-D[\[Phi]\_2[x, y, t], t]\) - \n\t\t\t\t\tU2\ D[\[Phi]\_2[x, y, t], x] - g\ \[Eta][x, t])\)\)], "Input"], Cell[BoxData[ FormBox[ RowBox[{ RowBox[{\(-\[Gamma]\), " ", RowBox[{ SuperscriptBox["\[Eta]", TagBox[\((2, 0)\), Derivative], MultilineFunction->None], "(", \(x, t\), ")"}]}], "-", RowBox[{\(\[Rho]\_1\), " ", RowBox[{"(", RowBox[{\(\(-g\)\ \(\[Eta](x, t)\)\), "-", RowBox[{ SubsuperscriptBox["\[Phi]", "1", TagBox[\((0, 0, 1)\), Derivative], MultilineFunction->None], "(", \(x, y, t\), ")"}], "-", RowBox[{"U1", " ", RowBox[{ SubsuperscriptBox["\[Phi]", "1", TagBox[\((1, 0, 0)\), Derivative], MultilineFunction->None], "(", \(x, y, t\), ")"}]}]}], ")"}]}], "+", RowBox[{\(\[Rho]\_2\), " ", RowBox[{"(", RowBox[{\(\(-g\)\ \(\[Eta](x, t)\)\), "-", RowBox[{ SubsuperscriptBox["\[Phi]", "2", TagBox[\((0, 0, 1)\), Derivative], MultilineFunction->None], "(", \(x, y, t\), ")"}], "-", RowBox[{"U2", " ", RowBox[{ SubsuperscriptBox["\[Phi]", "2", TagBox[\((1, 0, 0)\), Derivative], MultilineFunction->None], "(", \(x, y, t\), ")"}]}]}], ")"}]}]}], TraditionalForm]], "Output"] }, Open ]], Cell[TextData[{ "In ", StyleBox["Standard form", FontSlant->"Italic"], ", which is a little shorter, you would need the typeset window to make \ this practical. " }], "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(dc = \(-\((\[Gamma]\ \[PartialD]\_{x, 2}\[Eta][x, t])\)\) - \[Rho]\_1\ \((\(-\[PartialD]\_t \[Phi]\_1[x, y, t]\) - U1\ \[PartialD]\_x \[Phi]\_1[x, y, t] - g\ \[Eta][x, t])\) + \[Rho]\_2\ \((\(-\[PartialD]\_t \[Phi]\_2[x, y, t]\) - U2\ \[PartialD]\_x \[Phi]\_2[x, y, t] - g\ \[Eta][x, t])\)\)], "Input", CellTags->"standard"], Cell[BoxData[ FormBox[ RowBox[{ RowBox[{\(-\[Gamma]\), " ", RowBox[{ SuperscriptBox["\[Eta]", TagBox[\((2, 0)\), Derivative], MultilineFunction->None], "(", \(x, t\), ")"}]}], "-", RowBox[{\(\[Rho]\_1\), " ", RowBox[{"(", RowBox[{\(\(-g\)\ \(\[Eta](x, t)\)\), "-", RowBox[{ SubsuperscriptBox["\[Phi]", "1", TagBox[\((0, 0, 1)\), Derivative], MultilineFunction->None], "(", \(x, y, t\), ")"}], "-", RowBox[{"U1", " ", RowBox[{ SubsuperscriptBox["\[Phi]", "1", TagBox[\((1, 0, 0)\), Derivative], MultilineFunction->None], "(", \(x, y, t\), ")"}]}]}], ")"}]}], "+", RowBox[{\(\[Rho]\_2\), " ", RowBox[{"(", RowBox[{\(\(-g\)\ \(\[Eta](x, t)\)\), "-", RowBox[{ SubsuperscriptBox["\[Phi]", "2", TagBox[\((0, 0, 1)\), Derivative], MultilineFunction->None], "(", \(x, y, t\), ")"}], "-", RowBox[{"U2", " ", RowBox[{ SubsuperscriptBox["\[Phi]", "2", TagBox[\((1, 0, 0)\), Derivative], MultilineFunction->None], "(", \(x, y, t\), ")"}]}]}], ")"}]}]}], TraditionalForm]], "Output"] }, Open ]], Cell[TextData[ButtonBox["Back to start of calculations.", ButtonData:>"start", ButtonStyle->"Hyperlink"]], "Text"], Cell[TextData[{ "Finally it may be more easily ", StyleBox["read", FontVariations->{"Underline"->True}], " in ", StyleBox["Traditional form", FontSlant->"Italic"], " as shown here. Note that this particular cell is inactive and will not \ be evaluated. You can change this, if you really want to, by selecting the \ cell, going into the \"preferences\" under the Edit menu and under Cell \ Options, Evaluation Options, set Evaluatable to True. " }], "Text", CellTags->"traditional"], Cell[CellGroupData[{ Cell[BoxData[ FormBox[ RowBox[{"dc", "=", RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"\[Gamma]", " ", FractionBox[\(\[PartialD]\^2\( \[Eta](x, t)\)\), \(\[PartialD]x\^2\), MultilineFunction->None]}], ")"}]}], "-", RowBox[{\(\[Rho]\_1\), " ", RowBox[{"(", RowBox[{ RowBox[{"-", FractionBox[\(\[PartialD]\[Phi]\_1[x, y, t]\), \(\[PartialD]t\), MultilineFunction->None]}], "-", RowBox[{"U1", " ", FractionBox[\(\[PartialD]\[Phi]\_1[x, y, t]\), \(\[PartialD]x\), MultilineFunction->None]}], "-", \(g\ \(\[Eta](x, t)\)\)}], ")"}]}], "+", RowBox[{\(\[Rho]\_2\), " ", RowBox[{"(", RowBox[{ RowBox[{"-", FractionBox[\(\[PartialD]\[Phi]\_2[x, y, t]\), \(\[PartialD]t\), MultilineFunction->None]}], "-", RowBox[{"U2", " ", FractionBox[\(\[PartialD]\[Phi]\_2[x, y, t]\), \(\[PartialD]x\), MultilineFunction->None]}], "-", \(g\ \(\[Eta](x, t)\)\)}], ")"}]}]}]}], TraditionalForm]], "Input", Evaluatable->False, CellLabelAutoDelete->True], Cell[BoxData[ FormBox[ RowBox[{ RowBox[{\(-\[Gamma]\), " ", RowBox[{ SuperscriptBox["\[Eta]", TagBox[\((2, 0)\), Derivative], MultilineFunction->None], "(", \(x, t\), ")"}]}], "-", RowBox[{\(\[Rho]\_1\), " ", RowBox[{"(", RowBox[{\(\(-g\)\ \(\[Eta](x, t)\)\), "-", RowBox[{ SubsuperscriptBox["\[Phi]", "1", TagBox[\((0, 0, 1)\), Derivative], MultilineFunction->None], "(", \(x, y, t\), ")"}], "-", RowBox[{"U1", " ", RowBox[{ SubsuperscriptBox["\[Phi]", "1", TagBox[\((1, 0, 0)\), Derivative], MultilineFunction->None], "(", \(x, y, t\), ")"}]}]}], ")"}]}], "+", RowBox[{\(\[Rho]\_2\), " ", RowBox[{"(", RowBox[{\(\(-g\)\ \(\[Eta](x, t)\)\), "-", RowBox[{ SubsuperscriptBox["\[Phi]", "2", TagBox[\((0, 0, 1)\), Derivative], MultilineFunction->None], "(", \(x, y, t\), ")"}], "-", RowBox[{"U2", " ", RowBox[{ SubsuperscriptBox["\[Phi]", "2", TagBox[\((1, 0, 0)\), Derivative], MultilineFunction->None], "(", \(x, y, t\), ")"}]}]}], ")"}]}]}], TraditionalForm]], "Output", Evaluatable->False, CellLabelAutoDelete->True] }, Open ]], Cell[TextData[ButtonBox["Back to start of calculations.", ButtonData:>"start", ButtonStyle->"Hyperlink"]], "Text"], Cell["\<\ Note that you can convert any \"cryptic\" input expression by \ selecting the cell, going up to the cell menu and selecting \"convert to \ traditional form\". This is \"shift curly t\" on the keypad. \ \>", "Text"] }, Open ]] }, Closed]] }, Closed]], Cell[CellGroupData[{ Cell["Keywords", "Subtitle"], Cell["\<\ Natural Convection Heat Transfer, Thermal Boundary Layer, \ Boundary-Layer Flow, Similarity Variable, High Grashof Number Heat \ Transfer\ \>", "Text"] }, Closed]], Cell[CellGroupData[{ Cell["Problem overview", "Subtitle"], Cell[CellGroupData[{ Cell["Physical situation of interest", "Subsection"], Cell[CellGroupData[{ Cell[TextData[{ "The problem of interest is a heated sphere in a fluid of infinite extent \ that is otherwise not flowing. We expect that if the sphere is heated it \ will be at a higher temperature than the surrounding fluid, energy will be \ transferred to the fluid and the resulting temperature difference within the \ fluid will lead to a density gradient. If the sphere is in a uniform gravity \ field, you did not want pay to put your sphere on ", StyleBox["Mir", FontSlant->"Italic"], ", or you could not convince NASA to give you money to put it on the space \ shuttle, then the density gradient will lead to a body force acting on the \ fluid around the sphere. The region of low density fluid will be forced \ upward and replaced by higher density fluid. The question is, will the fluid \ flow?\n\nDimensional analysis for this problem would identify a ", StyleBox["Grashof", FontSlant->"Italic"], " number that represents the ratio of buoyancy force that would be causing \ flow over the viscous force that is resisting the flow. For slightly-heated \ road tar, we expect a low Grashof number and thus no flow. For a lighted \ light bulb in air, the Grashof number can be calculated as ", Cell[BoxData[ \(TraditionalForm\`\(\(10\^7\ or\ more\)\(,\)\)\)]], " this will certainly cause a strong flow.\n\nFor the case of low ", StyleBox["Gr", FontSlant->"Italic"], ", heat removal will be by conduction and radiation. If you define a \ Nusselt number in the usual way, hD/k, the value will be two. This is the \ pure condition limit. Since the Nusselt number gives the ratio of total heat \ transfer to pure condition (we are not considering radiation), we expect that \ as the Grashof number gets large, there will be a significant fluid flow and \ thus significant convective enhancement. " }], "Text", CellTags->"Grashof"], Cell[BoxData[ ButtonBox[\(return\ to\ conclusions\), ButtonData:>"conclusions", ButtonStyle->"Hyperlink"]], "Input"] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell["High Grashof number heat transfer", "Subsection"], Cell["\<\ For a Grashof number that is large compared to unity, we expect \ that the resistance to heat transfer, alternatively the region where most of \ the temperature change occurs, will be confined to a thin region close to the \ surface. If this region is thin compared to the macroscopic dimensions, we \ expect that boundary-layer approximations will give valid simplifications to \ the governing equations. This will be the case. Thus the boundary-layer \ equations will serve as the basic theory for the convection/conduction heat \ transfer from the heated sphere. \ \>", "Text", CellTags->"boundary_layer"], Cell[BoxData[ ButtonBox[\(return\ to\ conclusions\), ButtonData:>"conclusions", ButtonStyle->"Hyperlink"]], "Input"] }, Closed]], Cell[CellGroupData[{ Cell["Boundary-layer physics", "Subsection"], Cell[CellGroupData[{ Cell[TextData[{ "We recall from fluid dynamics that we expect boundary-layers to occur \ within the fluid close to solid surfaces at high Reynolds numbers. The \ current interest where the Grashof number is large is analogous. In a \ momentum boundary-layer the inertia and viscous forces are about the same \ order of magnitude. The governing equation thus contains the dominant \ inertia and viscous terms. As the Reynolds number is increased, the boundary \ layer gets thinner as \[Delta]/L ~ 1/Re^1/2. This scaling for this can be \ determined from the form of the equations with out solving any differential \ equations. The drag on the plate can be determined to order one accuracy \ simply by knowing the boundary layer thickness, \[Delta], and calculating \ \[Tau] = \[Mu] du/dy =~ \[Mu] U", Cell[BoxData[ \(TraditionalForm\`\_\[Infinity]\)]], "/\[Delta]. \n\nYou may not be familiar with heat transfer boundary layer \ behavior, but a lot of similarities occur. First now the conduction and \ convection terms are the same order of magnitude. Thus ", StyleBox["convective", FontSlant->"Italic"], " heat transfer involves both conduction and convection. You will notice \ below that the convection-conduction boundary-layer equation for energy, \ written in terms of T, looks quite similar to the flow direction \ boundary-layer equation for momentum. The primary difference is that the \ energy equation does not have a term like dp/dx, unless there is a heat \ source within the boundary layer. This is a source term and in the case of \ fluid flow is the reason the fluid is flowing (a momentum source). \n\n\ Another similarity between the present problem and fluid flow boundary layers \ is that we can get the basic scaling of the solution from the form of the \ equations. Thus once we have the equations, can figure out how to \ nondimensionalize them and finally turn the PDE's into ODE's, we will know \ how the Nusselt number scales with Grashof number without solving the \ equations! \n\nStrictly speaking boundary-layer theory is the first order \ result of some clever \"perturbation\" expansions of the differential \ equations. Leal shows this in his book, ( L. G. Leal (1992) ", StyleBox["Laminar Flow and Convective Transport Processes", FontSlant->"Italic"], ", Butterworth.) He also shows that for a flow around a sphere, or heat \ transfer from a sphere, the governing equations at the lowest order (which is \ what we are solving) are the same for a sphere or a flat plate. The main \ reason for this is because the layer is very thin compared to the radius of \ the sphere. Thus \"the earth is flat\"! -- or at least curvature does not \ matter. \n\nIf we derive the continuity, momentum and energy equations for \ flow past a heated flat plate we get " }], "Text"], Cell[BoxData[ ButtonBox[\(return\ to\ conclusions\), ButtonData:>"conclusions", ButtonStyle->"Hyperlink"]], "Input"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["\<\ The basic mass, momentum and energy equations for a free convection \ boundary - layer\ \>", "Subsection", CellTags->"govern_eqs"], Cell["Continuity", "Text"], Cell[BoxData[ FormBox[ RowBox[{ RowBox[{ FractionBox[\(\[PartialD]\(u(x, y)\)\), \(\[PartialD]x\), MultilineFunction->None], "+", FractionBox[\(\[PartialD]\(v(x, y)\)\), \(\[PartialD]y\), MultilineFunction->None]}], " ", "=", " ", "0"}], TraditionalForm]], "Input", Evaluatable->False], Cell["Momentum", "Text"], Cell[BoxData[ FormBox[ RowBox[{ RowBox[{ RowBox[{\(u(x, y)\), FractionBox[\(\[PartialD]\(u(x, y)\)\), \(\[PartialD]x\), MultilineFunction->None]}], " ", "+", " ", RowBox[{\(v(x, y)\), " ", FractionBox[\(\[PartialD]\(u(x, y)\)\), \(\[PartialD]y\), MultilineFunction->None]}]}], "=", " ", RowBox[{\(g\ \((T(x, y) - T0)\)\ \[Beta]\), " ", "+", RowBox[{"\[Nu]", " ", FractionBox[\(\[PartialD]\^2\( u(x, y)\)\), \(\[PartialD]y\^2\), MultilineFunction->None]}]}]}], TraditionalForm]], "Input", Evaluatable->False], Cell["Energy", "Text"], Cell[BoxData[ FormBox[ RowBox[{ RowBox[{ RowBox[{\(u(x, y)\), FractionBox[\(\[PartialD]\(T(x, y)\)\), \(\[PartialD]x\), MultilineFunction->None]}], " ", "+", " ", RowBox[{\(v(x, y)\), FractionBox[\(\[PartialD]\(T(x, y)\)\), \(\[PartialD]y\), MultilineFunction->None]}]}], " ", "=", " ", RowBox[{"\[Alpha]", " ", FractionBox[\(\[PartialD]\^2\( T(x, y)\)\), \(\[PartialD]y\^2\), MultilineFunction->None], " "}]}], TraditionalForm]], "Input", Evaluatable->False], Cell[BoxData[ ButtonBox[\(return\ to\ conclusions\), ButtonData:>"conclusions", ButtonStyle->"Hyperlink"]], "Input"], Cell["\<\ These equations are coupled, meaning you cannot solve any of them \ without simultaneously solving the other two. Further, the momentum equation \ is nonlinear. \ \>", "Text"] }, Closed]] }, Closed]], Cell[CellGroupData[{ Cell["Derivation of the ODE's from the original PDE's.", "Subtitle", Evaluatable->False, AspectRatioFixed->True], Cell[CellGroupData[{ Cell["Do background work to set up the problem", "Subsection"], Cell[CellGroupData[{ Cell["Existence of a similarity variable. ", "Subsubsection"], Cell[TextData[{ "A good reference for this section is L. G. Leal (1992) ", StyleBox["Laminar Flow and Convective Transport Processes", FontSlant->"Italic"], ", Butterworth.) \n\nIf you encounter a fluid flow problem in a \ rectangular channel, you expect the characteristic length to be the channel \ height. For a circular pipe, it will be the diameter or radius. For flow \ past a sphere, the radius is a reasonable length. Thus we always expect to \ find a characteristic length, that is contained in the physics of the flow. \ However, for a boundary-layer flow, where the boundary-layer is much smaller \ than the scale of the container, or a boundary-layer on the outside of an \ object where again the layer is much (say ~1/100) smaller than the device, \ what now is the relevant length scale? Would nature leave us with no \ characteristic length? \n\nThe answer is that there is a ", StyleBox["natural", FontSlant->"Italic"], " length scale. It is just not a ", StyleBox["geometric", FontSlant->"Italic"], " length scale. In the absence of a reasonable geometric length scale we \ look to define a \"similarity\" variable. The implication is that the \ velocity or temperature profiles are \"self-similar\", and with this proper \ scaling can be collapsed onto the same curve. The length scale for the \ normal direction within a boundary - layer is constructed from the physical \ properties of the flow and any other length or time variables that come to \ play in the problem.\n\nFor example, for the \"suddenly started\" plate \ problem, where a long and wide flat plate is started in motion, the \ boundary-layer thickness scales as ", Cell[BoxData[ \(TraditionalForm\`y\)]], "/", Cell[BoxData[ \(TraditionalForm\`\@\(\[Nu]\ t\)\)]], ", where ", StyleBox["t", FontSlant->"Italic"], " is time and \[Nu] is the kinematic viscosity. (The boundary layer gets \ thicker in time as ", Cell[BoxData[ \(TraditionalForm\`\@\(\[Nu]\ t\)\)]], ".) \n\nTo wrap this up, you always want to investigate the scaling of your \ problems. If there is no geometric length scale, (i.e., problem has a \ infinite spatial domain), you expect the existence of a similarity variable \ which is a combination of the other variables and physical properties. \n\n\ Use of a similarity variable can lead to transformations that turn the PDE's \ into ODE's. This can work for both linear and nonlinear equations. " }], "Text"] }, Open ]], Cell[CellGroupData[{ Cell["\<\ The basic mass, momentum and energy equations for a free convection \ boundary - layer\ \>", "Subsubsection"], Cell[CellGroupData[{ Cell[BoxData[ \(conteq\ = \ \(\(D[u[x, y], x]\)\(\ \)\(+\)\(D[v[x, y], y]\)\(\ \)\)\)], "Input"], Cell[BoxData[ FormBox[ RowBox[{ RowBox[{ SuperscriptBox["v", TagBox[\((0, 1)\), Derivative], MultilineFunction->None], "(", \(x, y\), ")"}], "+", RowBox[{ SuperscriptBox["u", TagBox[\((1, 0)\), Derivative], MultilineFunction->None], "(", \(x, y\), ")"}]}], TraditionalForm]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(umomeq\ = \ u[x, y]\ D[u[x, y], x]\ + \ v[x, y]\ D[u[x, y], y]\ - \ g\ \[Beta]\ \((T[x, y] - \ T0)\)\ - \ \[Nu]\ \ \ D[ u[x, y], {y, 2}]\)], "Input"], Cell[BoxData[ FormBox[ RowBox[{\(\(-g\)\ \[Beta]\ \((T(x, y) - T0)\)\), "+", RowBox[{\(v(x, y)\), " ", RowBox[{ SuperscriptBox["u", TagBox[\((0, 1)\), Derivative], MultilineFunction->None], "(", \(x, y\), ")"}]}], "-", RowBox[{"\[Nu]", " ", RowBox[{ SuperscriptBox["u", TagBox[\((0, 2)\), Derivative], MultilineFunction->None], "(", \(x, y\), ")"}]}], "+", RowBox[{\(u(x, y)\), " ", RowBox[{ SuperscriptBox["u", TagBox[\((1, 0)\), Derivative], MultilineFunction->None], "(", \(x, y\), ")"}]}]}], TraditionalForm]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(Teq\ = \ u[x, y]\ D[T[x, y], x]\ + \ v[x, y]\ D[T[x, y], y]\ - \ \[Alpha]\ D[T[x, y], {y, 2}]\)], "Input"], Cell[BoxData[ FormBox[ RowBox[{ RowBox[{\(v(x, y)\), " ", RowBox[{ SuperscriptBox["T", TagBox[\((0, 1)\), Derivative], MultilineFunction->None], "(", \(x, y\), ")"}]}], "-", RowBox[{"\[Alpha]", " ", RowBox[{ SuperscriptBox["T", TagBox[\((0, 2)\), Derivative], MultilineFunction->None], "(", \(x, y\), ")"}]}], "+", RowBox[{\(u(x, y)\), " ", RowBox[{ SuperscriptBox["T", TagBox[\((1, 0)\), Derivative], MultilineFunction->None], "(", \(x, y\), ")"}]}]}], TraditionalForm]], "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["Stream function formulation to satisfy continuity", "Subsubsection"], Cell["\<\ We want to turn these three coupled PDE's into 2 coupled ODE's that \ are nondimensionalized. We will define a stream function from the continuity \ equation, define a similarity variable involving x and y. Then transform the \ equations. \ \>", "Text"], Cell[TextData[{ "Note that continuity equation is satisfied automatically by a stream \ function formulation,\n u = ", Cell[BoxData[ \(TraditionalForm\`\[PartialD]\[Psi][x, y]\/\[PartialD]y\)]], " , v =- ", Cell[BoxData[ \(TraditionalForm\`\[PartialD]\[Psi][x, y]\/\[PartialD]x\)]], " " }], "Text"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"conteq", "/.", RowBox[{"{", RowBox[{\(u[x, y] \[Rule] D[\[Psi][x, y], y]\), ",", \(v[x, y] \[Rule] \(-D[\[Psi][x, y], x]\)\), ",", RowBox[{ RowBox[{ SuperscriptBox["u", TagBox[\((a1_, a2_)\), Derivative], MultilineFunction->None], "[", \(x, y\), "]"}], "\[RuleDelayed]", " ", \(\[PartialD]\_\({x, a1}, {y, a2}\)D[\[Psi][x, y], y]\)}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ SuperscriptBox["v", TagBox[\((a1_, a2_)\), Derivative], MultilineFunction->None], "[", \(x, y\), "]"}], "\[RuleDelayed]", " ", \(\[PartialD]\_\({x, a1}, {y, a2}\)D[\(-\[Psi][x, y]\), x]\)}]}], "}"}]}]], "Input"], Cell[BoxData[ \(TraditionalForm\`0\)], "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["Similarity variable formulation", "Subsubsection"], Cell["\<\ Now we need to do the variable transformations to make these last \ two equations dimensionless. For the simplest boundary - layer problem, \ start up of a Couette flow, this is relatively easy. For this problem it is \ a lot more complicated. \ \>", "Text"], Cell["\<\ Start by defining a similarity variable. We will just choose the \ correct one, there are formal procedures for finding these. (See Leal's book \ for how to do this.)\ \>", "Text"], Cell[BoxData[ \(eta = y/x \((Gr[x]\/4)\)\^\(1/4\)\)], "Input"], Cell[BoxData[ \(TraditionalForm\`\(y\ \@\(Gr(x)\)\%4\)\/\(\@2\ x\)\)], "Output", CellTags->"similarity_variable"], Cell[TextData[ButtonBox["Return to conclusions", ButtonData:>"conclusions", ButtonStyle->"Hyperlink"]], "Text"], Cell["Likewise for the stream function,", "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(psi = f[\[Eta][x, y]]\ 4\ \[Nu]\ \((Gr[x]\/4)\)\^\(1/4\)\)], "Input"], Cell[BoxData[ \(TraditionalForm\`2\ \@2\ \[Nu]\ \(f(\[Eta](x, y))\)\ \@\(Gr(x)\)\%4\)], "Output"] }, Open ]], Cell["\<\ We use the Greek, in the equations, but for substitution purposes, \ it is convenient to have a different name, eta and psi.\ \>", "Text"] }, Open ]], Cell[CellGroupData[{ Cell["\<\ Define the velocities in terms of the stream function and \ similarity variable\ \>", "Subsubsection"], Cell[TextData[{ "We need to use the chain rule to find all of the derivative substitutions \ and to find u and v. \n\nFirst find u[x,y], u = ", Cell[BoxData[ \(TraditionalForm\`\[PartialD]\[Psi][x, y]\/\[PartialD]y\)]], " = ", Cell[BoxData[ \(TraditionalForm\`\[PartialD]\[Psi]\/\[PartialD]\[Eta]\)]], Cell[BoxData[ \(TraditionalForm\`\[PartialD]\[Eta]\/\[PartialD]y\)]] }], "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(\(\(eu1\)\(=\)\(D[psi, y]\)\(\ \)\)\)], "Input"], Cell[BoxData[ FormBox[ RowBox[{"2", " ", \(\@2\), " ", "\[Nu]", " ", \(\@\(Gr(x)\)\%4\), " ", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "(", \(\[Eta](x, y)\), ")"}], " ", RowBox[{ SuperscriptBox["\[Eta]", TagBox[\((0, 1)\), Derivative], MultilineFunction->None], "(", \(x, y\), ")"}]}], TraditionalForm]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(eu1 = D[psi, \[Eta][x, y]] D[\[Eta][x, y], y]\)], "Input"], Cell[BoxData[ FormBox[ RowBox[{"2", " ", \(\@2\), " ", "\[Nu]", " ", \(\@\(Gr(x)\)\%4\), " ", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "(", \(\[Eta](x, y)\), ")"}], " ", RowBox[{ SuperscriptBox["\[Eta]", TagBox[\((0, 1)\), Derivative], MultilineFunction->None], "(", \(x, y\), ")"}]}], TraditionalForm]], "Output"] }, Open ]], Cell["\<\ Now get the y derivative of \[Eta] and substitute. \ \>", "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(eu2\ = \ eu1 /. D[\[Eta][x, y], y] \[Rule] D[eta, y]\)], "Input"], Cell[BoxData[ FormBox[ FractionBox[ RowBox[{"2", " ", "\[Nu]", " ", \(\@\(Gr(x)\)\), " ", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "(", \(\[Eta](x, y)\), ")"}]}], "x"], TraditionalForm]], "Output"] }, Open ]], Cell[TextData[{ "Now we need v.\n v =- ", Cell[BoxData[ \(TraditionalForm\`\[PartialD]\[Psi][x, y]\/\[PartialD]x\)]], " = - ", Cell[BoxData[ \(TraditionalForm\`\[PartialD]\[Psi]\/\[PartialD]\[Eta]\)]], Cell[BoxData[ \(TraditionalForm\`\[PartialD]\[Eta]\/\[PartialD]x\)]] }], "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(vee1 = \(-D[psi, x]\)\)], "Input"], Cell[BoxData[ FormBox[ RowBox[{ RowBox[{"-", FractionBox[ RowBox[{"\[Nu]", " ", \(f(\[Eta](x, y))\), " ", RowBox[{ SuperscriptBox["Gr", "\[Prime]", MultilineFunction->None], "(", "x", ")"}]}], \(\@2\ \(Gr(x)\)\^\(3/4\)\)]}], "-", RowBox[{ "2", " ", \(\@2\), " ", "\[Nu]", " ", \(\@\(Gr(x)\)\%4\), " ", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "(", \(\[Eta](x, y)\), ")"}], " ", RowBox[{ SuperscriptBox["\[Eta]", TagBox[\((1, 0)\), Derivative], MultilineFunction->None], "(", \(x, y\), ")"}]}]}], TraditionalForm]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(vee2 = vee1 /. D[\[Eta][x, y], x] \[Rule] D[eta, x]\)], "Input"], Cell[BoxData[ FormBox[ RowBox[{ RowBox[{"-", FractionBox[ RowBox[{"\[Nu]", " ", \(f(\[Eta](x, y))\), " ", RowBox[{ SuperscriptBox["Gr", "\[Prime]", MultilineFunction->None], "(", "x", ")"}]}], \(\@2\ \(Gr(x)\)\^\(3/4\)\)]}], "-", RowBox[{ "2", " ", \(\@2\), " ", "\[Nu]", " ", \(\@\(Gr(x)\)\%4\), " ", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "(", \(\[Eta](x, y)\), ")"}], " ", RowBox[{"(", RowBox[{ FractionBox[ RowBox[{"y", " ", RowBox[{ SuperscriptBox["Gr", "\[Prime]", MultilineFunction->None], "(", "x", ")"}]}], \(4\ \@2\ x\ \(Gr(x)\)\^\(3/4\)\)], "-", \(\(y\ \@\(Gr(x)\)\%4\)\/\(\@2\ x\^2\)\)}], ")"}]}]}], TraditionalForm]], "Output"] }, Open ]] }, Closed]] }, Closed]], Cell[CellGroupData[{ Cell["\<\ Now derive the 2 nonlinear ODE's that describe momentum and energy \ transfer\ \>", "Subsection", CellTags->"pde_to_ode"], Cell["\<\ It is easiest to do the substitutions starting with the original \ equations. We will see if it works. \ \>", "Text"], Cell[TextData[ButtonBox["Return to conclusions", ButtonData:>"conclusions", ButtonStyle->"Hyperlink"]], "Text"], Cell[CellGroupData[{ Cell["Here is the momentum equation", "Subsubsection"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"ueq1", "=", RowBox[{"umomeq", "/.", RowBox[{"{", RowBox[{\(u[x, y] \[Rule] eu2\), ",", \(v[x, y] \[Rule] vee2\), ",", RowBox[{ RowBox[{ SuperscriptBox["u", TagBox[\((a1_, a2_)\), Derivative], MultilineFunction->None], "[", \(x, y\), "]"}], "\[RuleDelayed]", " ", \(\[PartialD]\_\({x, a1}, {y, a2}\)eu2\)}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ SuperscriptBox["v", TagBox[\((a1_, a2_)\), Derivative], MultilineFunction->None], "[", \(x, y\), "]"}], "\[RuleDelayed]", " ", \(\[PartialD]\_\({x, a1}, {y, a2}\)vee2\)}]}], "}"}]}]}]], "Input"], Cell[BoxData[ FormBox[ RowBox[{ RowBox[{"-", FractionBox[ RowBox[{"2", " ", \(\@\(Gr(x)\)\), " ", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{ SuperscriptBox["f", TagBox[\((3)\), Derivative], MultilineFunction->None], "(", \(\[Eta](x, y)\), ")"}], " ", SuperscriptBox[ RowBox[{ SuperscriptBox["\[Eta]", TagBox[\((0, 1)\), Derivative], MultilineFunction->None], "(", \(x, y\), ")"}], "2"]}], "+", RowBox[{ RowBox[{ SuperscriptBox["f", "\[Prime]\[Prime]", MultilineFunction->None], "(", \(\[Eta](x, y)\), ")"}], " ", RowBox[{ SuperscriptBox["\[Eta]", TagBox[\((0, 2)\), Derivative], MultilineFunction->None], "(", \(x, y\), ")"}]}]}], ")"}], " ", \(\[Nu]\^2\)}], "x"]}], "+", RowBox[{\(1\/x\), RowBox[{"(", RowBox[{"2", " ", \(\@\(Gr(x)\)\), " ", RowBox[{"(", RowBox[{ RowBox[{"-", FractionBox[ RowBox[{"\[Nu]", " ", \(f(\[Eta](x, y))\), " ", RowBox[{ SuperscriptBox["Gr", "\[Prime]", MultilineFunction->None], "(", "x", ")"}]}], \(\@2\ \(Gr(x)\)\^\(3/4\)\)]}], "-", RowBox[{ "2", " ", \(\@2\), " ", "\[Nu]", " ", \(\@\(Gr(x)\)\%4\), " ", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "(", \(\[Eta](x, y)\), ")"}], " ", RowBox[{"(", RowBox[{ FractionBox[ RowBox[{"y", " ", RowBox[{ SuperscriptBox["Gr", "\[Prime]", MultilineFunction->None], "(", "x", ")"}]}], \(4\ \@2\ x\ \(Gr(x)\)\^\(3/4\)\)], "-", \(\(y\ \@\(Gr(x)\)\%4\)\/\(\@2\ x\^2\)\)}], ")"}]}]}], ")"}], " ", RowBox[{ SuperscriptBox["f", "\[Prime]\[Prime]", MultilineFunction->None], "(", \(\[Eta](x, y)\), ")"}], " ", RowBox[{ SuperscriptBox["\[Eta]", TagBox[\((0, 1)\), Derivative], MultilineFunction->None], "(", \(x, y\), ")"}], " ", "\[Nu]"}], ")"}]}], "+", RowBox[{\(1\/x\), RowBox[{"(", RowBox[{"2", " ", \(\@\(Gr(x)\)\), " ", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "(", \(\[Eta](x, y)\), ")"}], " ", RowBox[{"(", RowBox[{ FractionBox[ RowBox[{"\[Nu]", " ", RowBox[{ SuperscriptBox["Gr", "\[Prime]", MultilineFunction->None], "(", "x", ")"}], " ", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "(", \(\[Eta](x, y)\), ")"}]}], \(x\ \@\(Gr(x)\)\)], "-", FractionBox[ RowBox[{"2", " ", "\[Nu]", " ", \(\@\(Gr(x)\)\), " ", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "(", \(\[Eta](x, y)\), ")"}]}], \(x\^2\)], "+", FractionBox[ RowBox[{"2", " ", "\[Nu]", " ", \(\@\(Gr(x)\)\), " ", RowBox[{ SuperscriptBox["f", "\[Prime]\[Prime]", MultilineFunction->None], "(", \(\[Eta](x, y)\), ")"}], " ", RowBox[{ SuperscriptBox["\[Eta]", TagBox[\((1, 0)\), Derivative], MultilineFunction->None], "(", \(x, y\), ")"}]}], "x"]}], ")"}], " ", "\[Nu]"}], ")"}]}], "-", \(g\ \[Beta]\ \((T(x, y) - T0)\)\)}], TraditionalForm]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(ueq2 = \ ueq1 /. {D[\[Eta][x, y], x] \[Rule] D[eta, x], D[\[Eta][x, y], y] \[Rule] D[eta, y], \[IndentingNewLine]D[\[Eta][x, y], {x, 2}] \[Rule] D[eta, {x, 2}], D[\[Eta][x, y], {y, 2}] \[Rule] D[eta, {y, 2}]}\)], "Input"], Cell[BoxData[ FormBox[ RowBox[{ RowBox[{"-", FractionBox[ RowBox[{\(Gr(x)\), " ", RowBox[{ SuperscriptBox["f", TagBox[\((3)\), Derivative], MultilineFunction->None], "(", \(\[Eta](x, y)\), ")"}], " ", \(\[Nu]\^2\)}], \(x\^3\)]}], "+", RowBox[{\(1\/x\^2\), RowBox[{"(", RowBox[{\(\@2\), " ", \(\(Gr(x)\)\^\(3/4\)\), " ", RowBox[{"(", RowBox[{ RowBox[{"-", FractionBox[ RowBox[{"\[Nu]", " ", \(f(\[Eta](x, y))\), " ", RowBox[{ SuperscriptBox["Gr", "\[Prime]", MultilineFunction->None], "(", "x", ")"}]}], \(\@2\ \(Gr(x)\)\^\(3/4\)\)]}], "-", RowBox[{ "2", " ", \(\@2\), " ", "\[Nu]", " ", \(\@\(Gr(x)\)\%4\), " ", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "(", \(\[Eta](x, y)\), ")"}], " ", RowBox[{"(", RowBox[{ FractionBox[ RowBox[{"y", " ", RowBox[{ SuperscriptBox["Gr", "\[Prime]", MultilineFunction->None], "(", "x", ")"}]}], \(4\ \@2\ x\ \(Gr(x)\)\^\(3/4\)\)], "-", \(\(y\ \@\(Gr(x)\)\%4\)\/\(\@2\ x\^2\)\)}], ")"}]}]}], ")"}], " ", RowBox[{ SuperscriptBox["f", "\[Prime]\[Prime]", MultilineFunction->None], "(", \(\[Eta](x, y)\), ")"}], " ", "\[Nu]"}], ")"}]}], "+", RowBox[{\(1\/x\), RowBox[{"(", RowBox[{"2", " ", \(\@\(Gr(x)\)\), " ", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "(", \(\[Eta](x, y)\), ")"}], " ", RowBox[{"(", RowBox[{ FractionBox[ RowBox[{"\[Nu]", " ", RowBox[{ SuperscriptBox["Gr", "\[Prime]", MultilineFunction->None], "(", "x", ")"}], " ", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "(", \(\[Eta](x, y)\), ")"}]}], \(x\ \@\(Gr(x)\)\)], "-", FractionBox[ RowBox[{"2", " ", "\[Nu]", " ", \(\@\(Gr(x)\)\), " ", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "(", \(\[Eta](x, y)\), ")"}]}], \(x\^2\)], "+", FractionBox[ RowBox[{"2", " ", "\[Nu]", " ", \(\@\(Gr(x)\)\), " ", RowBox[{"(", RowBox[{ FractionBox[ RowBox[{"y", " ", RowBox[{ SuperscriptBox["Gr", "\[Prime]", MultilineFunction->None], "(", "x", ")"}]}], \(4\ \@2\ x\ \(Gr(x)\)\^\(3/4\)\)], "-", \(\(y\ \@\(Gr(x)\)\%4\)\/\(\@2\ x\^2\)\)}], ")"}], " ", RowBox[{ SuperscriptBox["f", "\[Prime]\[Prime]", MultilineFunction->None], "(", \(\[Eta](x, y)\), ")"}]}], "x"]}], ")"}], " ", "\[Nu]"}], ")"}]}], "-", \(g\ \[Beta]\ \((T(x, y) - T0)\)\)}], TraditionalForm]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(ueq3 = ueq2 /. \[Eta][x, y] \[Rule] \[Eta]\)], "Input"], Cell[BoxData[ FormBox[ RowBox[{ RowBox[{"-", FractionBox[ RowBox[{\(Gr(x)\), " ", RowBox[{ SuperscriptBox["f", TagBox[\((3)\), Derivative], MultilineFunction->None], "(", "\[Eta]", ")"}], " ", \(\[Nu]\^2\)}], \(x\^3\)]}], "+", FractionBox[ RowBox[{\(\@2\), " ", \(\(Gr(x)\)\^\(3/4\)\), " ", RowBox[{"(", RowBox[{ RowBox[{"-", FractionBox[ RowBox[{"\[Nu]", " ", \(f(\[Eta])\), " ", RowBox[{ SuperscriptBox["Gr", "\[Prime]", MultilineFunction->None], "(", "x", ")"}]}], \(\@2\ \(Gr(x)\)\^\(3/4\)\)]}], "-", RowBox[{ "2", " ", \(\@2\), " ", "\[Nu]", " ", \(\@\(Gr(x)\)\%4\), " ", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "(", "\[Eta]", ")"}], " ", RowBox[{"(", RowBox[{ FractionBox[ RowBox[{"y", " ", RowBox[{ SuperscriptBox["Gr", "\[Prime]", MultilineFunction->None], "(", "x", ")"}]}], \(4\ \@2\ x\ \(Gr(x)\)\^\(3/4\)\)], "-", \(\(y\ \@\(Gr(x)\)\%4\)\/\(\@2\ x\^2\)\)}], ")"}]}]}], ")"}], " ", RowBox[{ SuperscriptBox["f", "\[Prime]\[Prime]", MultilineFunction->None], "(", "\[Eta]", ")"}], " ", "\[Nu]"}], \(x\^2\)], "+", RowBox[{\(1\/x\), RowBox[{"(", RowBox[{"2", " ", \(\@\(Gr(x)\)\), " ", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "(", "\[Eta]", ")"}], " ", RowBox[{"(", RowBox[{ FractionBox[ RowBox[{"\[Nu]", " ", RowBox[{ SuperscriptBox["Gr", "\[Prime]", MultilineFunction->None], "(", "x", ")"}], " ", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "(", "\[Eta]", ")"}]}], \(x\ \@\(Gr(x)\)\)], "-", FractionBox[ RowBox[{"2", " ", "\[Nu]", " ", \(\@\(Gr(x)\)\), " ", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "(", "\[Eta]", ")"}]}], \(x\^2\)], "+", FractionBox[ RowBox[{"2", " ", "\[Nu]", " ", \(\@\(Gr(x)\)\), " ", RowBox[{"(", RowBox[{ FractionBox[ RowBox[{"y", " ", RowBox[{ SuperscriptBox["Gr", "\[Prime]", MultilineFunction->None], "(", "x", ")"}]}], \(4\ \@2\ x\ \(Gr(x)\)\^\(3/4\)\)], "-", \(\(y\ \@\(Gr(x)\)\%4\)\/\(\@2\ x\^2\)\)}], ")"}], " ", RowBox[{ SuperscriptBox["f", "\[Prime]\[Prime]", MultilineFunction->None], "(", "\[Eta]", ")"}]}], "x"]}], ")"}], " ", "\[Nu]"}], ")"}]}], "-", \(g\ \[Beta]\ \((T(x, y) - T0)\)\)}], TraditionalForm]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(ueq4\ = \ Expand[\ ueq3\ x^3/\[Nu]^2]\)], "Input"], Cell[BoxData[ FormBox[ RowBox[{\(-\(\(g\ \[Beta]\ \(T(x, y)\)\ x\^3\)\/\[Nu]\^2\)\), "+", \(\(g\ T0\ \[Beta]\ x\^3\)\/\[Nu]\^2\), "+", RowBox[{"2", " ", SuperscriptBox[ RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "(", "\[Eta]", ")"}], "2"], " ", RowBox[{ SuperscriptBox["Gr", "\[Prime]", MultilineFunction->None], "(", "x", ")"}], " ", "x"}], "-", RowBox[{\(f(\[Eta])\), " ", RowBox[{ SuperscriptBox["Gr", "\[Prime]", MultilineFunction->None], "(", "x", ")"}], " ", RowBox[{ SuperscriptBox["f", "\[Prime]\[Prime]", MultilineFunction->None], "(", "\[Eta]", ")"}], " ", "x"}], "-", RowBox[{"4", " ", \(Gr(x)\), " ", SuperscriptBox[ RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "(", "\[Eta]", ")"}], "2"]}], "-", RowBox[{\(Gr(x)\), " ", RowBox[{ SuperscriptBox["f", TagBox[\((3)\), Derivative], MultilineFunction->None], "(", "\[Eta]", ")"}]}]}], TraditionalForm]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(ueq5\ = \ \ ueq4 /. \(Gr'\)[x]\ \[Rule] 3\ \ \(\(Gr[x]/x\)\(\ \ \ \)\)\)], "Input"], Cell[BoxData[ FormBox[ RowBox[{\(-\(\(g\ \[Beta]\ \(T(x, y)\)\ x\^3\)\/\[Nu]\^2\)\), "+", \(\(g\ T0\ \[Beta]\ x\^3\)\/\[Nu]\^2\), "+", RowBox[{"2", " ", \(Gr(x)\), " ", SuperscriptBox[ RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "(", "\[Eta]", ")"}], "2"]}], "-", RowBox[{"3", " ", \(f(\[Eta])\), " ", \(Gr(x)\), " ", RowBox[{ SuperscriptBox["f", "\[Prime]\[Prime]", MultilineFunction->None], "(", "\[Eta]", ")"}]}], "-", RowBox[{\(Gr(x)\), " ", RowBox[{ SuperscriptBox["f", TagBox[\((3)\), Derivative], MultilineFunction->None], "(", "\[Eta]", ")"}]}]}], TraditionalForm]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(ueq6\ = \ ueq5 /. T[x, y] \[Rule] Tstar[\[Eta]]\ *\((\ Ts - T0)\)\ + \ T0\)], "Input"], Cell[BoxData[ FormBox[ RowBox[{\(-\(\(g\ \[Beta]\ \((T0 + \((Ts - T0)\)\ \(Tstar(\[Eta])\))\)\ x\^3\)\/\[Nu]\^2\)\), "+", \(\(g\ T0\ \[Beta]\ x\^3\)\/\[Nu]\^2\), "+", RowBox[{"2", " ", \(Gr(x)\), " ", SuperscriptBox[ RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "(", "\[Eta]", ")"}], "2"]}], "-", RowBox[{"3", " ", \(f(\[Eta])\), " ", \(Gr(x)\), " ", RowBox[{ SuperscriptBox["f", "\[Prime]\[Prime]", MultilineFunction->None], "(", "\[Eta]", ")"}]}], "-", RowBox[{\(Gr(x)\), " ", RowBox[{ SuperscriptBox["f", TagBox[\((3)\), Derivative], MultilineFunction->None], "(", "\[Eta]", ")"}]}]}], TraditionalForm]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(ufinal\ = \ Simplify[ueq6 /. {\[Beta] -> \ Gr[x]\ \(\(\[Nu]^2/g\)/\((Ts - T0)\)\)/x^3}\ ]/ Gr[x]\)], "Input"], Cell[BoxData[ FormBox[ RowBox[{ RowBox[{"2", " ", SuperscriptBox[ RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "(", "\[Eta]", ")"}], "2"]}], "-", \(Tstar(\[Eta])\), "-", RowBox[{"3", " ", \(f(\[Eta])\), " ", RowBox[{ SuperscriptBox["f", "\[Prime]\[Prime]", MultilineFunction->None], "(", "\[Eta]", ")"}]}], "-", RowBox[{ SuperscriptBox["f", TagBox[\((3)\), Derivative], MultilineFunction->None], "(", "\[Eta]", ")"}]}], TraditionalForm]], "Output"] }, Open ]] }, Closed]], Cell[TextData[ButtonBox["Return to conclusions", ButtonData:>"conclusions", ButtonStyle->"Hyperlink"]], "Text"], Cell[CellGroupData[{ Cell["Here is the energy equation", "Subsubsection"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Teq1", "=", RowBox[{"Teq", "/.", RowBox[{"{", RowBox[{\(u[x, y] \[Rule] eu2\), ",", \(v[x, y] \[Rule] vee2\), ",", RowBox[{ RowBox[{ SuperscriptBox["u", TagBox[\((a1_, a2_)\), Derivative], MultilineFunction->None], "[", \(x, y\), "]"}], "\[RuleDelayed]", " ", \(\[PartialD]\_\({x, a1}, {y, a2}\)eu2\)}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ SuperscriptBox["v", TagBox[\((a1_, a2_)\), Derivative], MultilineFunction->None], "[", \(x, y\), "]"}], "\[RuleDelayed]", " ", \(\[PartialD]\_\({x, a1}, {y, a2}\)vee2\)}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ SuperscriptBox["T", TagBox[\((a1_, a2_)\), Derivative], MultilineFunction->None], "[", \(x, y\), "]"}], "\[RuleDelayed]", " ", \(\[PartialD]\_\({x, a1}, {y, a2}\)Tstar[\[Eta][x, y]]\ \((Ts - T0)\)\)}]}], "}"}]}]}]], "Input"], Cell[BoxData[ FormBox[ RowBox[{ RowBox[{\((Ts - T0)\), " ", RowBox[{"(", RowBox[{ RowBox[{"-", FractionBox[ RowBox[{"\[Nu]", " ", \(f(\[Eta](x, y))\), " ", RowBox[{ SuperscriptBox["Gr", "\[Prime]", MultilineFunction->None], "(", "x", ")"}]}], \(\@2\ \(Gr(x)\)\^\(3/4\)\)]}], "-", RowBox[{ "2", " ", \(\@2\), " ", "\[Nu]", " ", \(\@\(Gr(x)\)\%4\), " ", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "(", \(\[Eta](x, y)\), ")"}], " ", RowBox[{"(", RowBox[{ FractionBox[ RowBox[{"y", " ", RowBox[{ SuperscriptBox["Gr", "\[Prime]", MultilineFunction->None], "(", "x", ")"}]}], \(4\ \@2\ x\ \(Gr(x)\)\^\(3/4\)\)], "-", \(\(y\ \@\(Gr(x)\)\%4\)\/\(\@2\ x\^2\)\)}], ")"}]}]}], ")"}], " ", RowBox[{ SuperscriptBox["Tstar", "\[Prime]", MultilineFunction->None], "(", \(\[Eta](x, y)\), ")"}], " ", RowBox[{ SuperscriptBox["\[Eta]", TagBox[\((0, 1)\), Derivative], MultilineFunction->None], "(", \(x, y\), ")"}]}], "-", RowBox[{\((Ts - T0)\), " ", "\[Alpha]", " ", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{ SuperscriptBox["Tstar", "\[Prime]\[Prime]", MultilineFunction->None], "(", \(\[Eta](x, y)\), ")"}], " ", SuperscriptBox[ RowBox[{ SuperscriptBox["\[Eta]", TagBox[\((0, 1)\), Derivative], MultilineFunction->None], "(", \(x, y\), ")"}], "2"]}], "+", RowBox[{ RowBox[{ SuperscriptBox["Tstar", "\[Prime]", MultilineFunction->None], "(", \(\[Eta](x, y)\), ")"}], " ", RowBox[{ SuperscriptBox["\[Eta]", TagBox[\((0, 2)\), Derivative], MultilineFunction->None], "(", \(x, y\), ")"}]}]}], ")"}]}], "+", FractionBox[ RowBox[{ "2", " ", \((Ts - T0)\), " ", "\[Nu]", " ", \(\@\(Gr(x)\)\), " ", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "(", \(\[Eta](x, y)\), ")"}], " ", RowBox[{ SuperscriptBox["Tstar", "\[Prime]", MultilineFunction->None], "(", \(\[Eta](x, y)\), ")"}], " ", RowBox[{ SuperscriptBox["\[Eta]", TagBox[\((1, 0)\), Derivative], MultilineFunction->None], "(", \(x, y\), ")"}]}], "x"]}], TraditionalForm]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(Teq2 = \ Teq1 /. {D[\[Eta][x, y], x] \[Rule] D[eta, x], D[\[Eta][x, y], y] \[Rule] D[eta, y], \[IndentingNewLine]D[\[Eta][x, y], {x, 2}] \[Rule] D[eta, {x, 2}], D[\[Eta][x, y], {y, 2}] \[Rule] D[eta, {y, 2}]}\)], "Input"], Cell[BoxData[ FormBox[ RowBox[{ FractionBox[ RowBox[{ "2", " ", \((Ts - T0)\), " ", "\[Nu]", " ", \(\@\(Gr(x)\)\), " ", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "(", \(\[Eta](x, y)\), ")"}], " ", RowBox[{"(", RowBox[{ FractionBox[ RowBox[{"y", " ", RowBox[{ SuperscriptBox["Gr", "\[Prime]", MultilineFunction->None], "(", "x", ")"}]}], \(4\ \@2\ x\ \(Gr(x)\)\^\(3/4\)\)], "-", \(\(y\ \@\(Gr(x)\)\%4\)\/\(\@2\ x\^2\)\)}], ")"}], " ", RowBox[{ SuperscriptBox["Tstar", "\[Prime]", MultilineFunction->None], "(", \(\[Eta](x, y)\), ")"}]}], "x"], "+", RowBox[{\(1\/\(\@2\ x\)\), RowBox[{"(", RowBox[{\((Ts - T0)\), " ", \(\@\(Gr(x)\)\%4\), " ", RowBox[{"(", RowBox[{ RowBox[{"-", FractionBox[ RowBox[{"\[Nu]", " ", \(f(\[Eta](x, y))\), " ", RowBox[{ SuperscriptBox["Gr", "\[Prime]", MultilineFunction->None], "(", "x", ")"}]}], \(\@2\ \(Gr(x)\)\^\(3/4\)\)]}], "-", RowBox[{ "2", " ", \(\@2\), " ", "\[Nu]", " ", \(\@\(Gr(x)\)\%4\), " ", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "(", \(\[Eta](x, y)\), ")"}], " ", RowBox[{"(", RowBox[{ FractionBox[ RowBox[{"y", " ", RowBox[{ SuperscriptBox["Gr", "\[Prime]", MultilineFunction->None], "(", "x", ")"}]}], \(4\ \@2\ x\ \(Gr(x)\)\^\(3/4\)\)], "-", \(\(y\ \@\(Gr(x)\)\%4\)\/\(\@2\ x\^2\)\)}], ")"}]}]}], ")"}], " ", RowBox[{ SuperscriptBox["Tstar", "\[Prime]", MultilineFunction->None], "(", \(\[Eta](x, y)\), ")"}]}], ")"}]}], "-", FractionBox[ RowBox[{\((Ts - T0)\), " ", "\[Alpha]", " ", \(\@\(Gr(x)\)\), " ", RowBox[{ SuperscriptBox["Tstar", "\[Prime]\[Prime]", MultilineFunction->None], "(", \(\[Eta](x, y)\), ")"}]}], \(2\ x\^2\)]}], TraditionalForm]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(Teq3 = Teq2 /. \[Eta][x, y] \[Rule] \[Eta]\)], "Input"], Cell[BoxData[ FormBox[ RowBox[{ FractionBox[ RowBox[{ "2", " ", \((Ts - T0)\), " ", "\[Nu]", " ", \(\@\(Gr(x)\)\), " ", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "(", "\[Eta]", ")"}], " ", RowBox[{"(", RowBox[{ FractionBox[ RowBox[{"y", " ", RowBox[{ SuperscriptBox["Gr", "\[Prime]", MultilineFunction->None], "(", "x", ")"}]}], \(4\ \@2\ x\ \(Gr(x)\)\^\(3/4\)\)], "-", \(\(y\ \@\(Gr(x)\)\%4\)\/\(\@2\ x\^2\)\)}], ")"}], " ", RowBox[{ SuperscriptBox["Tstar", "\[Prime]", MultilineFunction->None], "(", "\[Eta]", ")"}]}], "x"], "+", RowBox[{\(1\/\(\@2\ x\)\), RowBox[{"(", RowBox[{\((Ts - T0)\), " ", \(\@\(Gr(x)\)\%4\), " ", RowBox[{"(", RowBox[{ RowBox[{"-", FractionBox[ RowBox[{"\[Nu]", " ", \(f(\[Eta])\), " ", RowBox[{ SuperscriptBox["Gr", "\[Prime]", MultilineFunction->None], "(", "x", ")"}]}], \(\@2\ \(Gr(x)\)\^\(3/4\)\)]}], "-", RowBox[{ "2", " ", \(\@2\), " ", "\[Nu]", " ", \(\@\(Gr(x)\)\%4\), " ", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "(", "\[Eta]", ")"}], " ", RowBox[{"(", RowBox[{ FractionBox[ RowBox[{"y", " ", RowBox[{ SuperscriptBox["Gr", "\[Prime]", MultilineFunction->None], "(", "x", ")"}]}], \(4\ \@2\ x\ \(Gr(x)\)\^\(3/4\)\)], "-", \(\(y\ \@\(Gr(x)\)\%4\)\/\(\@2\ x\^2\)\)}], ")"}]}]}], ")"}], " ", RowBox[{ SuperscriptBox["Tstar", "\[Prime]", MultilineFunction->None], "(", "\[Eta]", ")"}]}], ")"}]}], "-", FractionBox[ RowBox[{\((Ts - T0)\), " ", "\[Alpha]", " ", \(\@\(Gr(x)\)\), " ", RowBox[{ SuperscriptBox["Tstar", "\[Prime]\[Prime]", MultilineFunction->None], "(", "\[Eta]", ")"}]}], \(2\ x\^2\)]}], TraditionalForm]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(Teq4\ = \ ExpandAll[\ Teq3\ \ ]\)], "Input"], Cell[BoxData[ FormBox[ RowBox[{ FractionBox[ RowBox[{"T0", " ", "\[Nu]", " ", \(f(\[Eta])\), " ", RowBox[{ SuperscriptBox["Gr", "\[Prime]", MultilineFunction->None], "(", "x", ")"}], " ", RowBox[{ SuperscriptBox["Tstar", "\[Prime]", MultilineFunction->None], "(", "\[Eta]", ")"}]}], \(2\ x\ \@\(Gr(x)\)\)], "-", FractionBox[ RowBox[{"Ts", " ", "\[Nu]", " ", \(f(\[Eta])\), " ", RowBox[{ SuperscriptBox["Gr", "\[Prime]", MultilineFunction->None], "(", "x", ")"}], " ", RowBox[{ SuperscriptBox["Tstar", "\[Prime]", MultilineFunction->None], "(", "\[Eta]", ")"}]}], \(2\ x\ \@\(Gr(x)\)\)], "+", FractionBox[ RowBox[{"T0", " ", "\[Alpha]", " ", \(\@\(Gr(x)\)\), " ", RowBox[{ SuperscriptBox["Tstar", "\[Prime]\[Prime]", MultilineFunction->None], "(", "\[Eta]", ")"}]}], \(2\ x\^2\)], "-", FractionBox[ RowBox[{"Ts", " ", "\[Alpha]", " ", \(\@\(Gr(x)\)\), " ", RowBox[{ SuperscriptBox["Tstar", "\[Prime]\[Prime]", MultilineFunction->None], "(", "\[Eta]", ")"}]}], \(2\ x\^2\)]}], TraditionalForm]], "Output"] }, Open ]], Cell[TextData[{ "Up to this point, we have not specified Gr[x]. To make x vanish from the \ equation, we need the following form for Gr[x]. Gr[x] ~", Cell[BoxData[ \(TraditionalForm\`\(\(\ \)\(x\^3\)\)\)]], ". " }], "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(Teq5\ = \ \ Teq4 /. \(Gr'\)[x]\ \[Rule] 3\ \ \(\(Gr[x]/x\)\(\ \ \ \)\)\)], "Input"], Cell[BoxData[ FormBox[ RowBox[{ FractionBox[ RowBox[{ "3", " ", "T0", " ", "\[Nu]", " ", \(f(\[Eta])\), " ", \(\@\(Gr(x)\)\), " ", RowBox[{ SuperscriptBox["Tstar", "\[Prime]", MultilineFunction->None], "(", "\[Eta]", ")"}]}], \(2\ x\^2\)], "-", FractionBox[ RowBox[{ "3", " ", "Ts", " ", "\[Nu]", " ", \(f(\[Eta])\), " ", \(\@\(Gr(x)\)\), " ", RowBox[{ SuperscriptBox["Tstar", "\[Prime]", MultilineFunction->None], "(", "\[Eta]", ")"}]}], \(2\ x\^2\)], "+", FractionBox[ RowBox[{"T0", " ", "\[Alpha]", " ", \(\@\(Gr(x)\)\), " ", RowBox[{ SuperscriptBox["Tstar", "\[Prime]\[Prime]", MultilineFunction->None], "(", "\[Eta]", ")"}]}], \(2\ x\^2\)], "-", FractionBox[ RowBox[{"Ts", " ", "\[Alpha]", " ", \(\@\(Gr(x)\)\), " ", RowBox[{ SuperscriptBox["Tstar", "\[Prime]\[Prime]", MultilineFunction->None], "(", "\[Eta]", ")"}]}], \(2\ x\^2\)]}], TraditionalForm]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(Teq6 = FullSimplify[ ExpandAll[2\ Teq5\ \ \(x^2/Sqrt[Gr[x]]\)\ /\[Nu]]]\)], "Input"], Cell[BoxData[ FormBox[ FractionBox[ RowBox[{\((T0 - Ts)\), " ", RowBox[{"(", RowBox[{ RowBox[{"3", " ", "\[Nu]", " ", \(f(\[Eta])\), " ", RowBox[{ SuperscriptBox["Tstar", "\[Prime]", MultilineFunction->None], "(", "\[Eta]", ")"}]}], "+", RowBox[{"\[Alpha]", " ", RowBox[{ SuperscriptBox["Tstar", "\[Prime]\[Prime]", MultilineFunction->None], "(", "\[Eta]", ")"}]}]}], ")"}]}], "\[Nu]"], TraditionalForm]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(Teq7 = Teq6/\((T0 - Ts)\)\ /. \[Alpha] \[Rule] \[Nu]/\((\ Pr)\)\)], "Input"], Cell[BoxData[ FormBox[ FractionBox[ RowBox[{ RowBox[{"3", " ", "\[Nu]", " ", \(f(\[Eta])\), " ", RowBox[{ SuperscriptBox["Tstar", "\[Prime]", MultilineFunction->None], "(", "\[Eta]", ")"}]}], "+", FractionBox[ RowBox[{"\[Nu]", " ", RowBox[{ SuperscriptBox["Tstar", "\[Prime]\[Prime]", MultilineFunction->None], "(", "\[Eta]", ")"}]}], "Pr"]}], "\[Nu]"], TraditionalForm]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(Tfinal = ExpandAll[Teq7\ \ \ \ Pr]\)], "Input"], Cell[BoxData[ FormBox[ RowBox[{ RowBox[{"3", " ", "Pr", " ", \(f(\[Eta])\), " ", RowBox[{ SuperscriptBox["Tstar", "\[Prime]", MultilineFunction->None], "(", "\[Eta]", ")"}]}], "+", RowBox[{ SuperscriptBox["Tstar", "\[Prime]\[Prime]", MultilineFunction->None], "(", "\[Eta]", ")"}]}], TraditionalForm]], "Output"] }, Open ]] }, Open ]], Cell[TextData[ButtonBox["Return to conclusions", ButtonData:>"conclusions", ButtonStyle->"Hyperlink"]], "Text"] }, Open ]], Cell[CellGroupData[{ Cell["The final boundary-Layer ODE's are", "Subsection"], Cell[CellGroupData[{ Cell[BoxData[ \(ufinal\)], "Input"], Cell[BoxData[ FormBox[ RowBox[{ RowBox[{"2", " ", SuperscriptBox[ RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "(", "\[Eta]", ")"}], "2"]}], "-", \(Tstar(\[Eta])\), "-", RowBox[{"3", " ", \(f(\[Eta])\), " ", RowBox[{ SuperscriptBox["f", "\[Prime]\[Prime]", MultilineFunction->None], "(", "\[Eta]", ")"}]}], "-", RowBox[{ SuperscriptBox["f", TagBox[\((3)\), Derivative], MultilineFunction->None], "(", "\[Eta]", ")"}]}], TraditionalForm]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(Tfinal\)], "Input"], Cell[BoxData[ FormBox[ RowBox[{ RowBox[{"3", " ", "Pr", " ", \(f(\[Eta])\), " ", RowBox[{ SuperscriptBox["Tstar", "\[Prime]", MultilineFunction->None], "(", "\[Eta]", ")"}]}], "+", RowBox[{ SuperscriptBox["Tstar", "\[Prime]\[Prime]", MultilineFunction->None], "(", "\[Eta]", ")"}]}], TraditionalForm]], "Output"] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell["Boundary conditions", "Subsection"], Cell[TextData[{ "The boundary conditions are that there is no slip on the surface, and that \ the temp and velocity relax to 0 far away and that the temperature on the \ surface is fixed\n\nThe are in terms of the newest notation:\n\n\[Eta] = 0, f \ = f' = 0, T* = 1,\n\[Eta]->\[Infinity], f' = 0, T* = 0. \n\nNote that f is \ the normal velocity, f' is the tangential velocity and T* = ", Cell[BoxData[ \(TraditionalForm\`\(\(\ \)\((T - T0)\)\)\/\((Ts - T0)\)\)]], ", where T0 is the far away temperature. " }], "Text"] }, Closed]], Cell[CellGroupData[{ Cell["Can we conclude anything at this point?", "Subsection"], Cell[TextData[ButtonBox["Return to conclusions", ButtonData:>"conclusions", ButtonStyle->"Hyperlink"]], "Text"], Cell[TextData[{ "We have nondimensional equations for the boundary-layer flow and \ temperature profile that depend only on the fluid property, the Prandtl \ number. Thus we expect to be able to calculate the dimensionless heat flux \ which is only a function of Pr. \n\nThe dimensional heat flux, must then be \ given by \"unscaling\" the equations. Note that doing this will involve ", Cell[BoxData[ \(TraditionalForm\`\(\(Gr\^\(1/4\)\)\(.\)\)\)]], " We thus conclude that this is the key scaling for convective heat \ transfer and this is how the data should scale (in the simplest case.) " }], "Text", CellTags->"gr_1/4"] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell["Correlation of the experimental data", "Subtitle"], Cell[TextData[ButtonBox["Return to conclusions", ButtonData:>"conclusions", ButtonStyle->"Hyperlink"]], "Text"], Cell[TextData[{ "There is a recommended correlation (from Churchill) for free convection \ heat transfer from a sphere:\n\n", Cell[BoxData[ \(TraditionalForm\`Nu\_D\)]], " = 2 + ", Cell[BoxData[ \(TraditionalForm\`\(0.589\ Ra\_D\^\(1/4\)\)\/\([1 + \ \((0.469/Pr)\)\^\ \(9/16\)]\)\^\(4/9\)\)]], ".\n\n(See page 504 of I &D.) \n\nThe Nusselt and Rayleigh numbers are \ based on the sphere diameter. Thus this correlation does have ", Cell[BoxData[ \(TraditionalForm\`Gr\^\(1/4\)\)]], " in agreement with the scaling of the boundary-layer equations. Note that \ Ra\[Congruent] Gr Pr. " }], "Text", CellTags->"correlation"] }, Closed]], Cell[CellGroupData[{ Cell["\<\ Numerical solution of the coupled ODE's that result after a \ similarity variable is introduced to the natural convection boundary - layer \ equations for flow near a flat plate.\ \>", "Subtitle"], Cell[CellGroupData[{ Cell["Set up the problem as a system of first order ODEs.", "Subsubsection"], Cell[BoxData[{ RowBox[{ StyleBox[\(\(From\ page\ 488\ in\ I\ &\) D, \ the\ \ boundary\ layer\ equations\ \ \(are : \[IndentingNewLine]\ \[IndentingNewLine]f''' + \ 3\ f\ f'' - \ 2\ f'\^2 + \ \ \(T\^*\)\), \[IndentingNewLine]\(T\^*\)''\ + \ 3\ Pr\ \ f\ \(T\^*\)'\ = \ 0\), FontFamily->"Times"], "\[IndentingNewLine]", StyleBox[" ", FontFamily->"Times"]}], "\[IndentingNewLine]", RowBox[{ RowBox[{ StyleBox["These", FontFamily->"Times"], StyleBox[" ", FontFamily->"Times"], StyleBox["nonlinear", FontFamily->"Times"], StyleBox[" ", FontFamily->"Times"], StyleBox["odes", FontFamily->"Times"], StyleBox[" ", FontFamily->"Times"], StyleBox["do", FontFamily->"Times"], StyleBox[" ", FontFamily->"Times"], StyleBox["not", FontFamily->"Times"], StyleBox[" ", FontFamily->"Times"], StyleBox["have", FontFamily->"Times"], StyleBox[" ", FontFamily->"Times"], StyleBox["a", FontFamily->"Times"], StyleBox[" ", FontFamily->"Times"], StyleBox["known", FontFamily->"Times"], StyleBox[" ", FontFamily->"Times"], StyleBox["analytical", FontFamily->"Times"], StyleBox[" ", FontFamily->"Times"], StyleBox["solution", FontFamily->"Times"], StyleBox[" ", FontFamily->"Times"], RowBox[{ StyleBox["(", FontFamily->"Times"], StyleBox[" ", FontFamily->"Times"], RowBox[{ StyleBox[\(for\ 2\ particular\ values\ of\ Pr\), FontFamily->"Times"], StyleBox[",", FontFamily->"Times"], StyleBox[" ", FontFamily->"Times"], StyleBox[\(a\ similar\ set\ of\ ODE' s\ does\ have\ an\ analytic\ solution\), FontFamily->"Times"], StyleBox[",", FontFamily->"Times"], StyleBox[" ", FontFamily->"Times"], RowBox[{ StyleBox["See", FontFamily->"Times"], StyleBox[":", FontFamily->"Times"], StyleBox[" ", FontFamily->"Times"], RowBox[{ StyleBox[\(L . \ G . \ Leal\), "SmallText", FontFamily->"Times"], StyleBox[" ", "SmallText", FontFamily->"Times"], StyleBox[\((1992)\), "SmallText", FontFamily->"Times"], StyleBox[" ", "SmallText", FontFamily->"Times"], StyleBox["Laminar", "SmallText", FontFamily->"Times", FontWeight->"Plain", FontSlant->"Italic"], StyleBox[" ", "SmallText", FontFamily->"Times", FontWeight->"Plain", FontSlant->"Italic"], StyleBox["flow", "SmallText", FontFamily->"Times", FontWeight->"Plain", FontSlant->"Italic"], StyleBox[" ", "SmallText", FontFamily->"Times", FontWeight->"Plain", FontSlant->"Italic"], StyleBox["and", "SmallText", FontFamily->"Times", FontWeight->"Plain", FontSlant->"Italic"], StyleBox[" ", "SmallText", FontFamily->"Times", FontWeight->"Plain", FontSlant->"Italic"], StyleBox["Convective", "SmallText", FontFamily->"Times", FontWeight->"Plain", FontSlant->"Italic"], StyleBox[" ", "SmallText", FontFamily->"Times", FontWeight->"Plain", FontSlant->"Italic"], StyleBox["Transport", "SmallText", FontFamily->"Times", FontWeight->"Plain", FontSlant->"Italic"], StyleBox[" ", "SmallText", FontFamily->"Times", FontWeight->"Plain", FontSlant->"Italic"], StyleBox["Processes", "SmallText", FontFamily->"Times", FontWeight->"Plain", FontSlant->"Italic"]}]}], StyleBox[",", "SmallText", FontFamily->"Times", FontWeight->"Plain"], StyleBox[" ", "SmallText", FontFamily->"Times", FontWeight->"Plain"], StyleBox[\(\(Butterworth\)\(.\)\), "SmallText", FontFamily->"Times", FontWeight->"Plain"]}], StyleBox[")", FontFamily->"Times"]}]}], StyleBox[" ", FontFamily->"Times"]}], "\[IndentingNewLine]", StyleBox[\(However, they\ can\ be\ easily\ solved\ numerically . \ \ \ \ Numerical\ \ solutions\ for\ ode' s\ are\ often\ done\ by\ creating\ \ a\ system\ of\ first\ order\ \ odes . This\ is\ done\ by\ defining\), FontFamily->"Times"], "\[IndentingNewLine]", StyleBox[\(y1 = f, \[IndentingNewLine]y2 = f'\), FontFamily->"Times"], "\[IndentingNewLine]", StyleBox[\(y3 = f''\), FontFamily->"Times"], "\[IndentingNewLine]", StyleBox[\(y4 = T\), FontFamily->"Times"], "\[IndentingNewLine]", StyleBox[\(y5 = T', \[IndentingNewLine]\[IndentingNewLine]These\ have\ to\ be\ \ related . \ \ Thus\ we\ have\), FontFamily->"Times"], "\[IndentingNewLine]", StyleBox[\(y1' = y2\ \((by\ definition)\)\), FontFamily->"Times"], "\[IndentingNewLine]", StyleBox[\(y2' = y3\ \((by\ definition)\)\), FontFamily->"Times"], "\[IndentingNewLine]", StyleBox[\(y3' = \(-3\)\ y1*y3 + 2\ y2\^2 - y4\ \((which\ is\ from\ the\ original\ \(\(ODE\)\(.\)\))\)\), FontFamily->"Times"], "\[IndentingNewLine]", StyleBox[\(y4' = y5\), FontFamily->"Times"], "\[IndentingNewLine]", StyleBox[\(y5' = \(-\((y1\ y5)\)\) 3\ Pr \((\ which\ is\ the\ other\ \(\(ODE\)\(.\)\))\)\), FontFamily->"Times"]}], "Text", CellMargins->{{18, Inherited}, {Inherited, Inherited}}, TextAlignment->Left] }, Closed]], Cell[CellGroupData[{ Cell[TextData[{ "Numerical solution uses a shooting method and a Runge - Kutta Integration. \ \nReference: Numerical Recipes in Fortran, See also the notes from ChEg \ 258,(D. T. Leighton). ", ButtonBox["http://www.nd.edu/~dtl/cheg258/cheg258-1999/notes/l37/overheads.\ html", ButtonData:>{ URL[ "http://www.nd.edu/~dtl/cheg258/cheg258-1999/notes/l37/overheads.html"], None}, ButtonStyle->"Hyperlink"] }], "Subsubsection"], Cell[TextData[ButtonBox["Return to conclusions", ButtonData:>"conclusions", ButtonStyle->"Hyperlink"]], "Text"], Cell[TextData[{ "The key equations for the shooting method are:\n\n", StyleBox["\[Alpha]", FontWeight->"Bold"], " \[CenterDot] \[Delta] ", StyleBox["V", FontWeight->"Bold"], " = -", StyleBox["F", FontWeight->"Bold"], "\n", StyleBox["V", FontWeight->"Bold"], "new = ", StyleBox["V", FontWeight->"Bold"], "old + \[Delta]", StyleBox["V\n\n", FontWeight->"Bold"], "where ", StyleBox["\[Alpha]", FontWeight->"Bold"], " is the Jacobian Matrix obtained by varying the initial guesses for the \ unknown initial conditions, ", StyleBox["F", FontWeight->"Bold"], " is the error in the solution produced from the current initial guesses, \ ", StyleBox["V", FontWeight->"Bold"], "old, and \[Delta]", StyleBox["V", FontWeight->"Bold"], " is the correction to the initial guesses for the next step. \n\nThis \ turns out to be a touchy calculation. Here are a set of initial conditions \ that give a solution with some trouble. You will find that higher ", StyleBox["Pr", FontSlant->"Italic"], " is harder and that if youter is too large, the calculation does not work.\ \n\nfinit = {{.8},{-.25}}; \npr=10; (*Prandtl number*)\nyouter=5.5; (* outer \ value of \[Eta] *)\ndeltv = .0006;(* increment on the initial guesses used \ to generate the Jacobian*)" }], "Text", CellTags->"shooting"] }, Closed]], Cell[CellGroupData[{ Cell["Here is the numerical solution", "Subsection"], Cell[CellGroupData[{ Cell["Initialize some variables.", "Subsubsection"], Cell[BoxData[{ RowBox[{\(finit\ = \ {{ .8}, {\(- .25\)}};\), " ", StyleBox[\( (*initial\ condition\ vector\ for\ the\ 2\ unknown\ \ derivatives, \ \(f''\)[0]\ and\ \(T'\)[0]*) \), FontFamily->"Comic Sans MS", FontWeight->"Plain"], StyleBox[" ", FontFamily->"Comic Sans MS", FontWeight->"Plain"]}], "\[IndentingNewLine]", RowBox[{\(pr = 1;\), " ", StyleBox[\( (*Prandtl\ number*) \), FontFamily->"Comic Sans MS", FontWeight->"Plain"]}], "\[IndentingNewLine]", RowBox[{\(youter = 5;\), " ", StyleBox[\( (*\ outer\ value\ of\ \[Eta]\ *) \), FontFamily->"Comic Sans MS", FontWeight->"Plain"]}], "\[IndentingNewLine]", RowBox[{\(deltv\ = \ .0001;\), StyleBox[\( (*\ increment\ on\ the\ initial\ guesses\ used\ to\ generate\ the\ \ Jacobian*) \), FontFamily->"Comic Sans MS", FontWeight->"Plain"]}]}], "Input", AspectRatioFixed->True], Cell[CellGroupData[{ Cell["\<\ If you like to know how to pick out a number from a vector, look \ here. \ \>", "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(finit\)], "Input"], Cell[BoxData[ FormBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"0.8`"}, {\(-0.25`\)} }, ColumnAlignments->{Decimal}], "\[NoBreak]", ")"}], TraditionalForm]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(finit[\([1, 1]\)]\)], "Input"], Cell[BoxData[ \(TraditionalForm\`0.8`\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(finit[\([2, 1]\)]\)], "Input"], Cell[BoxData[ \(TraditionalForm\`\(-0.25`\)\)], "Output"] }, Open ]] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell["\<\ We need to run through the integration once to get a first values \ for the error. \ \>", "Subsubsection", Evaluatable->False, AspectRatioFixed->True], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"zz", "=", RowBox[{"NDSolve", "[", "\[IndentingNewLine]", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ SuperscriptBox["y1", "\[Prime]", MultilineFunction->None], "[", "x", "]"}], "==", \(y2[x]\)}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ SuperscriptBox["y2", "\[Prime]", MultilineFunction->None], "[", "x", "]"}], "==", \(y3[x]\)}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ SuperscriptBox["y3", "\[Prime]", MultilineFunction->None], "[", "x", "]"}], "==", \(\(-3\)\ \ y1[x]\ y3[x] + 2\ \ y2[x]^2\ - y4[x]\)}], ",", "\[IndentingNewLine]", \(\(y4'\)[x]\ \[Equal] \ y5[x]\), ",", "\[IndentingNewLine]", \(\(y5'\)[x]\ \[Equal] \ 3\ pr\ \((\ \(-\ y1[x]\)\ y5[x])\)\), ",", \(y1[0] == 0\), ",", \(y2[0] \[Equal] 0\), ",", \(y3[0] \[Equal] finit[\([1, 1]\)]\), ",", "\[IndentingNewLine]", \(y4[0] \[Equal] 1\), ",", \(y5[0]\ \[Equal] finit[\([2, 1]\)]\)}], "}"}], ",", "\[IndentingNewLine]", \({y1[x], y2[x], y3[x], y4[x], y5[x]}\), ",", \({x, 0, youter}\)}], "]"}]}], ";"}], "\n", \(eps1 = \(-\((zz /. x \[Rule] youter)\)\[LeftDoubleBracket]1, 2, 2\[RightDoubleBracket]\);\), "\[IndentingNewLine]", \(eps2 = \ \(-\((zz /. x \[Rule] youter)\)\[LeftDoubleBracket]1, 4, 2\[RightDoubleBracket]\);\), "\n", \(Print[\ "\", eps1, "\< f''[0]= \>", finit[\([1, 1]\)], "\< Terror= \>", \ eps2, "\< Tprime[0] = \>", \ finit[\([2, 1]\)]];\)}], "Input", AspectRatioFixed->True], Cell[BoxData[ FormBox[ InterpretationBox[\("Uerror= "\[InvisibleSpace]0.6768205491772687`\ \[InvisibleSpace]" f''[0]= "\[InvisibleSpace]0.8`\[InvisibleSpace]" \ Terror= "\[InvisibleSpace]\(-0.6210401410465517`\)\[InvisibleSpace]" \ Tprime[0] = "\[InvisibleSpace]\(-0.25`\)\), SequenceForm[ "Uerror= ", 0.67682054917726875, " f''[0]= ", 0.80000000000000004, " Terror= ", -0.6210401410465517, " Tprime[0] = ", -0.25], Editable->False], TraditionalForm]], "Print"] }, Open ]], Cell["\<\ The streamwise velocity and dimensionless temperature should be 0 \ as \[Eta]--> \[Infinity]. Of course we cannot not calculate out to very high \ values of \[Eta]. We need just high enough \[Eta]. \ \>", "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(eps1\)], "Input"], Cell[BoxData[ \(TraditionalForm\`0.6768205491772687`\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(eps2\)], "Input"], Cell[BoxData[ \(TraditionalForm\`\(-0.6210401410465517`\)\)], "Output"] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell["Here is the main iteration loop.", "Subsubsection"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"While", "[", RowBox[{\(Abs[eps1] + Abs[eps2] > .0001\), ",", "\[IndentingNewLine]", StyleBox[\( (*calculate\ the\ change\ due\ to\ changing\ the\ \ velocity\ gradient\ guess*) \), FontFamily->"Comic Sans MS", FontWeight->"Plain"], "\[IndentingNewLine]", RowBox[{ RowBox[{"zz", "=", RowBox[{"NDSolve", "[", "\[IndentingNewLine]", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ SuperscriptBox["y1", "\[Prime]", MultilineFunction->None], "[", "x", "]"}], "==", \(y2[x]\)}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ SuperscriptBox["y2", "\[Prime]", MultilineFunction->None], "[", "x", "]"}], "==", \(y3[x]\)}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ SuperscriptBox["y3", "\[Prime]", MultilineFunction->None], "[", "x", "]"}], "==", \(\(-3\)\ \ y1[x]\ y3[x] + 2\ \ y2[x]^2\ - y4[x]\)}], ",", "\[IndentingNewLine]", \(\(y4'\)[x]\ \[Equal] \ y5[x]\), ",", "\[IndentingNewLine]", \(\(y5'\)[x]\ \[Equal] \ 3\ pr\ \((\ \(-\ y1[x]\)\ y5[x])\)\), ",", \(y1[0] == 0\), ",", \(y2[0] \[Equal] 0\), ",", \(y3[0] \[Equal] finit[\([1, 1]\)] + deltv\), ",", "\[IndentingNewLine]", \(y4[0] \[Equal] 1\), ",", \(y5[0]\ \[Equal] finit[\([2, 1]\)]\)}], "}"}], ",", "\[IndentingNewLine]", \({y1[x], y2[x], y3[x], y4[x], y5[x]}\), ",", \({x, 0, youter}\)}], "]"}]}], ";", "\[IndentingNewLine]", \(eps12 = \(-\((zz /. x \[Rule] youter)\)\[LeftDoubleBracket]1, 2, 2\[RightDoubleBracket]\)\), ";", "\[IndentingNewLine]", \(eps22 = \(-\((zz /. x \[Rule] youter)\)\[LeftDoubleBracket]1, 4, 2\[RightDoubleBracket]\)\), ";", "\[IndentingNewLine]", StyleBox[\( (*calculate\ the\ change\ due\ to\ changing\ the\ \ temperature\ gradient\ guess*) \), FontFamily->"Comic Sans MS", FontWeight->"Plain"], "\[IndentingNewLine]", RowBox[{"zz", "=", RowBox[{"NDSolve", "[", "\[IndentingNewLine]", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ SuperscriptBox["y1", "\[Prime]", MultilineFunction->None], "[", "x", "]"}], "==", \(y2[x]\)}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ SuperscriptBox["y2", "\[Prime]", MultilineFunction->None], "[", "x", "]"}], "==", \(y3[x]\)}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ SuperscriptBox["y3", "\[Prime]", MultilineFunction->None], "[", "x", "]"}], "==", \(\(-3\)\ \ y1[x]\ y3[x] + 2\ \ y2[x]^2\ - y4[x]\)}], ",", "\[IndentingNewLine]", \(\(y4'\)[x]\ \[Equal] \ y5[x]\), ",", "\[IndentingNewLine]", \(\(y5'\)[x]\ \[Equal] \ 3\ pr\ \((\ \(-\ y1[x]\)\ y5[x])\)\), ",", \(y1[0] == 0\), ",", \(y2[0] \[Equal] 0\), ",", \(y3[0] \[Equal] finit[\([1, 1]\)]\), ",", "\[IndentingNewLine]", \(y4[0] \[Equal] 1\), ",", \(y5[0]\ \[Equal] finit[\([2, 1]\)] + deltv\)}], "}"}], ",", "\[IndentingNewLine]", \({y1[x], y2[x], y3[x], y4[x], y5[x]}\), ",", \({x, 0, youter}\)}], "]"}]}], ";", "\[IndentingNewLine]", \(eps13 = \(-\((zz /. x \[Rule] youter)\)\[LeftDoubleBracket]1, 2, 2\[RightDoubleBracket]\)\), ";", "\[IndentingNewLine]", \(eps23 = \(-\((zz /. x \[Rule] youter)\)\[LeftDoubleBracket]1, 4, 2\[RightDoubleBracket]\)\), ";", "\[IndentingNewLine]", StyleBox[\( (*calculate\ the\ Jacobian*) \), FontFamily->"Comic Sans MS", FontWeight->"Plain"], "\[IndentingNewLine]", \(alf = {{\((eps12 - eps1)\)/ deltv, \((eps13 - eps1)\)/ deltv}, \ \ \ \ \ \ \[IndentingNewLine]\ \ \ \ \ \ \ \ \ \ \ {\((eps22 - eps2)\)/deltv, \((eps23 - eps2)\)/deltv}}\), ";", "\[IndentingNewLine]", StyleBox[\( (*calculate\ the\ new\ guesses*) \), FontFamily->"Comic Sans MS", FontWeight->"Plain"], "\[IndentingNewLine]", \(deltvnew\ = \ \(-Inverse[ alf]\ . \ {{eps1}, {eps2}}\)\), ";", "\[IndentingNewLine]", \(finitnew\ = \ finit\ + \ deltvnew\), ";", "\[IndentingNewLine]", StyleBox[\( (*\(u pdate\)\ everything\ and\ save\ the\ old\ numbers*) \), FontFamily->"Comic Sans MS", FontWeight->"Plain"], "\[IndentingNewLine]", \(finitold = \ finit\), ";", "\[IndentingNewLine]", \(eps1old\ = \ eps1\), ";", "\[IndentingNewLine]", \(eps2old\ = \ eps2\), ";", "\[IndentingNewLine]", \(finit = finitnew\), ";", "\[IndentingNewLine]", StyleBox[\( (*run\ the\ new\ numbers\ to\ see\ if\ the\ solution\ \ has\ converged*) \), FontFamily->"Comic Sans MS", FontWeight->"Plain"], "\[IndentingNewLine]", RowBox[{"zz", "=", RowBox[{"NDSolve", "[", "\[IndentingNewLine]", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ SuperscriptBox["y1", "\[Prime]", MultilineFunction->None], "[", "x", "]"}], "==", \(y2[x]\)}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ SuperscriptBox["y2", "\[Prime]", MultilineFunction->None], "[", "x", "]"}], "==", \(y3[x]\)}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ SuperscriptBox["y3", "\[Prime]", MultilineFunction->None], "[", "x", "]"}], "==", \(\(-3\)\ \ y1[x]\ y3[x] + 2\ \ y2[x]^2\ - y4[x]\)}], ",", "\[IndentingNewLine]", \(\(y4'\)[x]\ \[Equal] \ y5[x]\), ",", "\[IndentingNewLine]", \(\(y5'\)[x]\ \[Equal] \ 3\ pr\ \((\ \(-\ y1[x]\)\ y5[x])\)\), ",", \(y1[0] == 0\), ",", \(y2[0] \[Equal] 0\), ",", \(y3[0] \[Equal] finit[\([1, 1]\)]\), ",", "\[IndentingNewLine]", \(y4[0] \[Equal] 1\), ",", \(y5[0]\ \[Equal] finit[\([2, 1]\)]\)}], "}"}], ",", "\[IndentingNewLine]", \({y1[x], y2[x], y3[x], y4[x], y5[x]}\), ",", \({x, 0, youter}\)}], "]"}]}], ";", "\n", \(eps1 = \(-\((zz /. x \[Rule] youter)\)\[LeftDoubleBracket]1, 2, 2\[RightDoubleBracket]\)\), ";", "\[IndentingNewLine]", \(eps2 = \(-\((zz /. x \[Rule] youter)\)\[LeftDoubleBracket]1, 4, 2\[RightDoubleBracket]\)\), ";", "\n", \(Print[\ "\", eps1, "\< f''[0]= \>", finit[\([1, 1]\)], "\< Terror= \>", \ eps2, "\< Tprime[0] = \>", \ finit[\([2, 1]\)]]\), ";"}]}], "]"}]], "Input"], Cell[BoxData[ FormBox[ InterpretationBox[\("Uerror= \ "\[InvisibleSpace]\(-0.5060173373935412`\)\[InvisibleSpace]" f''[0]= "\ \[InvisibleSpace]0.657017452666272`\[InvisibleSpace]" Terror= "\ \[InvisibleSpace]0.058349432930894035`\[InvisibleSpace]" Tprime[0] = "\ \[InvisibleSpace]\(-0.638140869477595`\)\), SequenceForm[ "Uerror= ", -0.50601733739354116, " f''[0]= ", 0.65701745266627198, " Terror= ", 0.058349432930894035, " Tprime[0] = ", -0.63814086947759496], Editable->False], TraditionalForm]], "Print"], Cell[BoxData[ FormBox[ InterpretationBox[\("Uerror= "\[InvisibleSpace]\(-0.03571593879669796`\ \)\[InvisibleSpace]" f''[0]= "\[InvisibleSpace]0.6409520362958993`\ \[InvisibleSpace]" Terror= "\[InvisibleSpace]0.008388823852632244`\ \[InvisibleSpace]" Tprime[0] = "\[InvisibleSpace]\(-0.573493809034766`\)\), SequenceForm[ "Uerror= ", -0.03571593879669796, " f''[0]= ", 0.64095203629589925, " Terror= ", 0.0083888238526322444, " Tprime[0] = ", -0.57349380903476599], Editable->False], TraditionalForm]], "Print"], Cell[BoxData[ FormBox[ InterpretationBox[\("Uerror= \ "\[InvisibleSpace]\(-0.0003063076687180981`\)\[InvisibleSpace]" f''[0]= "\ \[InvisibleSpace]0.6413873842297924`\[InvisibleSpace]" Terror= "\ \[InvisibleSpace]0.00006302155579737879`\[InvisibleSpace]" Tprime[0] = "\ \[InvisibleSpace]\(-0.5669086838482106`\)\), SequenceForm[ "Uerror= ", -0.00030630766871809811, " f''[0]= ", 0.64138738422979236, " Terror= ", 6.3021555797378788*^-05, " Tprime[0] = ", -0.56690868384821058], Editable->False], TraditionalForm]], "Print"], Cell[BoxData[ FormBox[ InterpretationBox[\("Uerror= "\[InvisibleSpace]2.7098373300492385`*^-8\ \[InvisibleSpace]" f''[0]= "\[InvisibleSpace]0.641388168222365`\ \[InvisibleSpace]" Terror= "\[InvisibleSpace]\(-7.713707112053856`*^-9\)\ \[InvisibleSpace]" Tprime[0] = \ "\[InvisibleSpace]\(-0.5668529780770206`\)\), SequenceForm[ "Uerror= ", 2.7098373300492385*^-08, " f''[0]= ", 0.641388168222365, " Terror= ", -7.7137071120538559*^-09, " Tprime[0] = ", -0.56685297807702062], Editable->False], TraditionalForm]], "Print"] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell["Here are some plots of the results", "Subsubsection"], Cell[CellGroupData[{ Cell[BoxData[ \(\(Upr1 = \ Plot[y2[x] /. zz\[LeftDoubleBracket]1\[RightDoubleBracket], {x, 0, youter}, PlotRange \[Rule] All, AxesLabel \[Rule] {"\<\[Eta]\>", \*"\"\<\!\(\(u\\\\\\\ x\)\/\(2 \ \[Nu]\)\) \!\(Gr\^\(-\(1\/2\)\)\)\>\""}];\)\)], "Input", AspectRatioFixed->True], 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 0.190476 0.0147152 2.3495 [ [.21429 .00222 -3 -9 ] [.21429 .00222 3 0 ] [.40476 .00222 -3 -9 ] [.40476 .00222 3 0 ] [.59524 .00222 -3 -9 ] [.59524 .00222 3 0 ] [.78571 .00222 -3 -9 ] [.78571 .00222 3 0 ] [.97619 .00222 -3 -9 ] [.97619 .00222 3 0 ] [1.025 .01472 0 -6 ] [1.025 .01472 10 6 ] [.01131 .13219 -24 -4.5 ] [.01131 .13219 0 4.5 ] [.01131 .24967 -18 -4.5 ] [.01131 .24967 0 4.5 ] [.01131 .36714 -24 -4.5 ] [.01131 .36714 0 4.5 ] [.01131 .48462 -18 -4.5 ] [.01131 .48462 0 4.5 ] [.01131 .60209 -24 -4.5 ] [.01131 .60209 0 4.5 ] [.02381 .64303 -24.7188 0 ] [.02381 .64303 24.7188 18.6875 ] [ 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 [(1)] .21429 .00222 0 1 Mshowa .40476 .01472 m .40476 .02097 L s [(2)] .40476 .00222 0 1 Mshowa .59524 .01472 m .59524 .02097 L s [(3)] .59524 .00222 0 1 Mshowa .78571 .01472 m .78571 .02097 L s [(4)] .78571 .00222 0 1 Mshowa .97619 .01472 m .97619 .02097 L s [(5)] .97619 .00222 0 1 Mshowa .125 Mabswid .0619 .01472 m .0619 .01847 L s .1 .01472 m .1 .01847 L s .1381 .01472 m .1381 .01847 L s .17619 .01472 m .17619 .01847 L s .25238 .01472 m .25238 .01847 L s .29048 .01472 m .29048 .01847 L s .32857 .01472 m .32857 .01847 L s .36667 .01472 m .36667 .01847 L s .44286 .01472 m .44286 .01847 L s .48095 .01472 m .48095 .01847 L s .51905 .01472 m .51905 .01847 L s .55714 .01472 m .55714 .01847 L s .63333 .01472 m .63333 .01847 L s .67143 .01472 m .67143 .01847 L s .70952 .01472 m .70952 .01847 L s .74762 .01472 m .74762 .01847 L s .82381 .01472 m .82381 .01847 L s .8619 .01472 m .8619 .01847 L s .9 .01472 m .9 .01847 L s .9381 .01472 m .9381 .01847 L s .25 Mabswid 0 .01472 m 1 .01472 L s gsave 1.025 .01472 -61 -10 Mabsadd m 1 1 Mabs scale currentpoint translate /MISOfy { /newfontname exch def /oldfontname exch def oldfontname findfont dup length dict begin {1 index /FID ne {def} {pop pop} ifelse} forall /Encoding ISOLatin1Encoding def currentdict end newfontname exch definefont pop } def 0 20 translate 1 -1 scale 63.000 13.000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor 0.000 0.000 rmoveto 63.000 13.000 moveto %%IncludeResource: font Math1Mono %%IncludeFont: Math1Mono /Math1Mono findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (h) show 69.000 13.000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor 0.000 0.000 rmoveto 1.000 setlinewidth grestore .02381 .13219 m .03006 .13219 L s [(0.05)] .01131 .13219 1 0 Mshowa .02381 .24967 m .03006 .24967 L s [(0.1)] .01131 .24967 1 0 Mshowa .02381 .36714 m .03006 .36714 L s [(0.15)] .01131 .36714 1 0 Mshowa .02381 .48462 m .03006 .48462 L s [(0.2)] .01131 .48462 1 0 Mshowa .02381 .60209 m .03006 .60209 L s [(0.25)] .01131 .60209 1 0 Mshowa .125 Mabswid .02381 .03821 m .02756 .03821 L s .02381 .06171 m .02756 .06171 L s .02381 .0852 m .02756 .0852 L s .02381 .1087 m .02756 .1087 L s .02381 .15569 m .02756 .15569 L s .02381 .17918 m .02756 .17918 L s .02381 .20268 m .02756 .20268 L s .02381 .22617 m .02756 .22617 L s .02381 .27316 m .02756 .27316 L s .02381 .29666 m .02756 .29666 L s .02381 .32015 m .02756 .32015 L s .02381 .34365 m .02756 .34365 L s .02381 .39064 m .02756 .39064 L s .02381 .41413 m .02756 .41413 L s .02381 .43763 m .02756 .43763 L s .02381 .46112 m .02756 .46112 L s .02381 .50811 m .02756 .50811 L s .02381 .53161 m .02756 .53161 L s .02381 .5551 m .02756 .5551 L s .02381 .5786 m .02756 .5786 L s .25 Mabswid .02381 0 m .02381 .61803 L s gsave .02381 .64303 -85.7188 -4 Mabsadd m 1 1 Mabs scale currentpoint translate /MISOfy { /newfontname exch def /oldfontname exch def oldfontname findfont dup length dict begin {1 index /FID ne {def} {pop pop} ifelse} forall /Encoding ISOLatin1Encoding def currentdict end newfontname exch definefont pop } def 0 26.6875 translate 1 -1 scale 63.000 15.438 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor 0.000 0.000 rmoveto 66.062 9.500 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (u) show 74.688 9.500 moveto (x) show %%IncludeResource: font Math1Mono %%IncludeFont: Math1Mono /Math1Mono findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor 64.750 15.438 moveto (\\200\\200\\200\\200) show 73.250 15.438 moveto (\\200\\200) show 77.500 15.438 moveto (\\200) show 79.688 15.438 moveto (\\200) show 66.062 21.688 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (2) show 72.062 21.688 moveto %%IncludeResource: font Math1Mono %%IncludeFont: Math1Mono /Math1Mono findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor ( ) show 74.688 21.688 moveto %%IncludeResource: font Math1Mono %%IncludeFont: Math1Mono /Math1Mono findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (n) show 83.812 15.438 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (Gr) show 95.812 11.438 moveto %%IncludeResource: font Math1Mono %%IncludeFont: Math1Mono /Math1Mono findfont 7.125 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (-) show 102.188 8.438 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 5.062 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (1) show %%IncludeResource: font Math1Mono %%IncludeFont: Math1Mono /Math1Mono findfont 7.125 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor 101.250 11.438 moveto (\\200\\200) show 104.125 11.438 moveto (\\200) show 104.438 11.438 moveto (\\200) show 102.188 13.375 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 5.062 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (2) show 108.438 15.438 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor 0.000 0.000 rmoveto 1.000 setlinewidth grestore .5 Mabswid .02381 .01472 m .04262 .15227 L .06244 .27396 L .08426 .38213 L .10458 .46066 L .12297 .51492 L .14264 .55708 L .15285 .57309 L .16394 .58632 L .16902 .59102 L .17441 .59512 L .17946 .59818 L .18408 .60032 L .18639 .60117 L .18891 .60193 L .19128 .60249 L .19349 .60288 L .19482 .60306 L .19606 .60318 L .19718 .60326 L .19839 .60331 L .19902 .60332 L .19972 .60332 L .20044 .60331 L .20111 .60329 L .20232 .60322 L .20362 .60311 L .20484 .60297 L .20597 .60281 L .20853 .60234 L .21079 .60181 L .21319 .60112 L .22192 .5976 L .2268 .59498 L .23129 .59218 L .24139 .58462 L .25967 .56703 L .29929 .51619 L .3374 .4579 L .37796 .39293 L .41701 .3323 L .4585 .27317 L .49849 .22296 L .53695 .18139 L .57787 .14425 L .61727 .11482 L .65516 .09172 L .6955 .07192 L .73432 .05679 L .7756 .04419 L Mistroke .81536 .0348 L .85361 .02781 L .89431 .02213 L .93349 .01802 L .97512 .01479 L .97619 .01472 L Mfstroke 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath % End of Graphics MathPictureEnd \ \>"], "Graphics", Evaluatable->False, AspectRatioFixed->True, ImageSize->{281.438, 173.875}, ImageMargins->{{34, 0}, {0, 0}}, ImageRegion->{{0, 1}, {0, 1}}, ImageCache->GraphicsData["Bitmap", "\<\ CF5dJ6E]HGAYHf4PAg9QL6QYHgo6V>o001o02MS_`04001S_f>o000VHkl0 1@00HkmS_f>o0000:F>o00<006>oHkl09F>o00D006>oHkmS_`0001eS_`00CF>o00<006>oHkl0:6>o 00<006>oHkl0:V>o00<006>oHkl096>o1@00:V>o00<006>oHkl06f>o001=Hkl00`00HkmS_`0YHkl0 0`00HkmS_`0WHkl2000WHkl01000HkmS_`00:f>o00<006>oHkl06f>o001=Hkl00`00HkmS_`0VHkl0 1@00HkmS_f>o0000:V>o00<006>oHkl09F>o00<006>o00009f>o10007V>o001o00009V>o00D006>oHkmS_`0002QS_`8002MS_`03001S_f>o01mS_`00CF>o00<006>oHkl0 9f>o0`00:6>o0`00:V>o00<006>oHkl09F>o1@0046>o00<006>oHkl02V>o000SHkl00`00HkmS_`3V Hkl00`00HkmS_`0:Hkl002=S_`03001S_f>o0>=S_`04001S_f>o000o0>=S _`05001S_f>oHkl0000;Hkl001eS_n@000US_`04001S_f>o000;Hkl002=S_`03001S_f>o00ES_`03 001S_f>o00IS_`03001S_f>o00IS_`03001S_f>o00ES_`03001S_f>o00IS_`03001S_f>o00IS_`03 001S_f>o00ES_`03001S_f>o00IS_`03001S_f>o00ES_`03001S_f>o00IS_`03001S_f>o00IS_`03 001S_f>o00ES_`03001S_f>o00IS_`03001S_f>o00IS_`03001S_f>o00ES_`03001S_f>o00IS_`03 001S_f>o00ES_`03001S_f>o00IS_`03001S_f>o00IS_`03001S_f>o00ES_`03001S_f>o00IS_`03 001S_f>o00IS_`03001S_f>o00ES_`04001S_f>oHkl90004Hkl00`00HkmS_`0=Hkl200000f>o001S _`0:Hkl002=S_`03001S_f>o02QS_`03001S_f>o02QS_`03001S_f>o02QS_`03001S_f>o02QS_`03 001S_f>o019S_`T000eS_`03001S_f>o00]S_`<00003Hkl0000000]S_`008f>o00<006>oHkl0]f>o 1`00=F>o000SHkl00`00HkmS_`2bHkl5000lHkl002=S_`800:mS_`@0045S_`008f>o0P00ZV>o1@00 AF>o000SHkl2002WHkl3001:Hkl002=S_`800:AS_`<004eS_`008f>o0P00XF>o0`00D6>o000SHkl2 002NHkl3001CHkl002=S_`8009]S_`<005IS_`008f>o0P00V6>o0`00FF>o000SHkl00`00Hkl0002D Hkl3001LHkl002=S_`03001S_`00099S_`8005mS_`008f>o00<006>o0000T6>o0P00HF>o000SHkl3 002?Hkl00`00HkmS_`1QHkl002=S_`03001S_`0008eS_`8006AS_`008f>o00<006>o0000Rf>o0P00 IV>o000SHkl00`00Hkl00029Hkl2001XHkl002=S_`03001S_`0008MS_`8006YS_`008f>o00@006>o Hkl008ES_`03001S_f>o06YS_`008f>o0P0000=S_`00Hkl0PV>o0P00KF>o000SHkl01000HkmS_`00 PF>o0P00Kf>o000SHkl01000HkmS_`00P6>o00<006>oHkl0Kf>o000:Hkl20004Hkl00`00HkmS_`03 Hkl20003Hkl30005Hkl01000HkmS_`00Of>o00<006>oHkl0L6>o0009Hkl01000HkmS_`0026>o00H0 06>oHkl006>o0003Hkl00`00HkmS_`02Hkl01000HkmS_`00OF>o0P00Lf>o0009Hkl01000HkmS_`00 26>o00@006>oHkl000ES_`03001S_f>o009S_`800003Hkl006>o07]S_`03001S_f>o07=S_`002F>o 00@006>oHkl000QS_`04001S_f>o0005Hkl00`00HkmS_`02Hkl01000HkmS_`00Nf>o00<006>oHkl0 M6>o0009Hkl01000HkmS_`0026>o00D006>oHkl006>o00@000ES_`05001S_f>oHkl0001iHkl00`00 HkmS_`1eHkl000US_`04001S_f>o0008Hkl01P00HkmS_`00Hkl000QS_`05001S_f>oHkl0001gHkl2 001hHkl000YS_`8000YS_`80009S_`D000AS_`05001S_f>oHkl0001fHkl00`00HkmS_`1hHkl002=S _`80009S_`03001S_f>o07=S_`03001S_f>o07US_`008f>o00D006>oHkmS_`0007AS_`03001S_f>o 07YS_`008f>o00D006>oHkmS_`0007=S_`03001S_f>o07]S_`008f>o00D006>oHkmS_`00079S_`03 001S_f>o07aS_`008f>o00<006>oHkl00V>o00<006>oHkl0KV>o00<006>oHkl0OF>o000SHkl00`00 HkmS_`02Hkl00`00HkmS_`1]Hkl00`00HkmS_`1nHkl002=S_`8000=S_`03001S_f>o06aS_`03001S _f>o07mS_`008f>o00<006>oHkl00V>o00<006>oHkl0K6>o00<006>oHkl0Of>o000SHkl00`00HkmS _`02Hkl00`00HkmS_`1[Hkl00`00HkmS_`20Hkl002=S_`03001S_f>o009S_`03001S_f>o06YS_`03 001S_f>o085S_`008f>o00<006>oHkl00V>o00<006>oHkl0JF>o00<006>oHkl0PV>o000SHkl20004 Hkl00`00HkmS_`1WHkl00`00HkmS_`23Hkl002=S_`03001S_f>o00=S_`03001S_f>o06IS_`03001S _f>o08AS_`008f>o00<006>oHkl00f>o00<006>oHkl0IF>o00<006>oHkl0QF>o000SHkl00`00HkmS _`03Hkl00`00HkmS_`1THkl00`00HkmS_`26Hkl002=S_`03001S_f>o00=S_`03001S_f>o06=S_`03 001S_f>o08MS_`008f>o0P0016>o00<006>oHkl0HV>o00<006>oHkl0R6>o000SHkl00`00HkmS_`03 Hkl00`00HkmS_`1RHkl00`00HkmS_`28Hkl002=S_`03001S_f>o00AS_`03001S_f>o061S_`03001S _f>o08US_`0046>o0P0016>o00<006>oHkl00V>o0`001F>o00<006>oHkl016>o00<006>oHkl0Gf>o 00<006>oHkl0RV>o000?Hkl01000HkmS_`002F>o00<006>oHkl016>o00<006>oHkl016>o00<006>o Hkl0GV>o00<006>oHkl0Rf>o000?Hkl01000HkmS_`002F>o00<006>oHkl016>o00<006>oHkl016>o 00<006>oHkl0GF>o00<006>oHkl0S6>o000?Hkl01000HkmS_`002F>o00<006>oHkl016>o0P001F>o 00<006>oHkl0G6>o00<006>oHkl0SF>o000?Hkl01000HkmS_`002F>o00<006>oHkl016>o00<006>o Hkl016>o00<006>oHkl0G6>o00<006>oHkl0SF>o000?Hkl01000HkmS_`0026>o0P001V>o00<006>o Hkl016>o00<006>oHkl0Ff>o00<006>oHkl0SV>o000@Hkl2000:Hkl00`00HkmS_`04Hkl00`00HkmS _`05Hkl00`00HkmS_`1IHkl00`00HkmS_`2?Hkl002=S_`03001S_f>o00ES_`03001S_f>o05QS_`03 001S_f>o091S_`008f>o0P001V>o00<006>oHkl0Ef>o00<006>oHkl0TF>o000SHkl00`00HkmS_`05 Hkl00`00HkmS_`1GHkl00`00HkmS_`2AHkl002=S_`03001S_f>o00ES_`03001S_f>o05IS_`03001S _f>o099S_`008f>o00<006>oHkl01V>o00<006>oHkl0E6>o00<006>oHkl0Tf>o000SHkl00`00HkmS _`06Hkl00`00HkmS_`1DHkl00`00HkmS_`2CHkl002=S_`8000MS_`03001S_f>o05=S_`03001S_f>o 09AS_`008f>o00<006>oHkl01V>o00<006>oHkl0DV>o00<006>oHkl0UF>o000SHkl00`00HkmS_`06 Hkl00`00HkmS_`1AHkl00`00HkmS_`2FHkl002=S_`03001S_f>o00MS_`03001S_f>o051S_`03001S _f>o09IS_`008f>o00<006>oHkl01f>o00<006>oHkl0Cf>o00<006>oHkl0Uf>o000SHkl00`00HkmS _`07Hkl00`00HkmS_`1>Hkl00`00HkmS_`2HHkl002=S_`8000QS_`03001S_f>o04iS_`03001S_f>o 09QS_`008f>o00<006>oHkl01f>o00<006>oHkl0CF>o00<006>oHkl0VF>o000SHkl00`00HkmS_`08 Hkl00`00HkmS_`1;Hkl00`00HkmS_`2JHkl002=S_`03001S_f>o00QS_`03001S_f>o04]S_`03001S _f>o09YS_`008f>o00<006>oHkl026>o00<006>oHkl0BV>o00<006>oHkl0Vf>o000SHkl20009Hkl0 0`00HkmS_`19Hkl00`00HkmS_`2LHkl002=S_`03001S_f>o00QS_`03001S_f>o04US_`03001S_f>o 09aS_`008f>o00<006>oHkl02F>o00<006>oHkl0Af>o00<006>oHkl0WF>o000:Hkl20004Hkl00`00 HkmS_`02Hkl30003Hkl30005Hkl00`00HkmS_`09Hkl00`00HkmS_`16Hkl00`00HkmS_`2NHkl000US _`04001S_f>o0009Hkl01@00HkmS_f>o00000f>o00<006>oHkl00V>o00<006>oHkl02F>o00<006>o Hkl0AV>o00<006>oHkl0WV>o0009Hkl01000HkmS_`002F>o00<006>oHkl01F>o00<006>oHkl00V>o 0P002V>o00<006>oHkl0AF>o00<006>oHkl0Wf>o0009Hkl01000HkmS_`002F>o00<006>oHkl01F>o 00<006>oHkl00V>o00<006>oHkl02F>o00<006>oHkl0AF>o00<006>oHkl0Wf>o0009Hkl01000HkmS _`002F>o00@006>oHkmS_`@000ES_`03001S_f>o00YS_`03001S_f>o04=S_`03001S_f>o0:1S_`00 2F>o00@006>oHkl000QS_`8000=S_`03001S_f>o00IS_`03001S_f>o00YS_`03001S_f>o049S_`03 001S_f>o0:5S_`002V>o0P002V>o00@006>oHkmS_`D000AS_`03001S_f>o00YS_`03001S_f>o049S _`03001S_f>o0:5S_`008f>o00<006>oHkl02V>o00<006>oHkl0@F>o00<006>oHkl0XV>o000SHkl2 000;Hkl00`00HkmS_`10Hkl00`00HkmS_`2SHkl002=S_`03001S_f>o00]S_`03001S_f>o03mS_`03 001S_f>o0:=S_`008f>o00<006>oHkl02f>o00<006>oHkl0?V>o00<006>oHkl0Y6>o000SHkl00`00 HkmS_`0;Hkl00`00HkmS_`0mHkl00`00HkmS_`2UHkl002=S_`03001S_f>o00aS_`03001S_f>o03aS _`03001S_f>o0:ES_`008f>o0P003F>o00<006>oHkl0>f>o00<006>oHkl0YV>o000SHkl00`00HkmS _`0o00eS_`03001S_f>o03US_`03 001S_f>o0:MS_`008f>o00<006>oHkl03F>o00<006>oHkl0>6>o00<006>oHkl0Z6>o000SHkl00`00 HkmS_`0=Hkl00`00HkmS_`0hHkl00`00HkmS_`2XHkl002=S_`8000iS_`03001S_f>o03MS_`03001S _f>o0:US_`008f>o00<006>oHkl03V>o00<006>oHkl0=V>o00<006>oHkl0ZF>o000SHkl00`00HkmS _`0>Hkl00`00HkmS_`0eHkl00`00HkmS_`2ZHkl002=S_`03001S_f>o00iS_`03001S_f>o03AS_`03 001S_f>o0:]S_`008f>o00<006>oHkl03f>o00<006>oHkl0o00<006>oHkl0Zf>o000SHkl00`00 HkmS_`0?Hkl00`00HkmS_`0bHkl00`00HkmS_`2/Hkl002=S_`80011S_`03001S_f>o039S_`03001S _f>o0:aS_`008f>o00<006>oHkl03f>o00<006>oHkl0o00<006>oHkl0[F>o000@Hkl20004Hkl0 1000HkmS_f>o1@0016>o00<006>oHkl046>o00<006>oHkl0;f>o00<006>oHkl0[V>o000?Hkl01000 HkmS_`0026>o00@006>oHkl000AS_`03001S_f>o011S_`03001S_f>o02mS_`03001S_f>o0:iS_`00 3f>o00@006>oHkl000US_`03001S_f>o00AS_`03001S_f>o011S_`03001S_f>o02iS_`03001S_f>o 0:mS_`003f>o00@006>oHkl000YS_`03001S_f>o00=S_`80019S_`03001S_f>o02aS_`03001S_f>o 0;1S_`003f>o00@006>oHkl000MS_`05001S_f>oHkl00004Hkl00`00HkmS_`0AHkl00`00HkmS_`0[ Hkl00`00HkmS_`2aHkl000mS_`04001S_f>o0007Hkl01@00HkmS_f>o000016>o00<006>oHkl04F>o 00<006>oHkl0:f>o00<006>oHkl0/F>o000@Hkl20009Hkl30005Hkl00`00HkmS_`0BHkl00`00HkmS _`0YHkl00`00HkmS_`2bHkl002=S_`03001S_f>o019S_`03001S_f>o02QS_`03001S_f>o0;=S_`00 8f>o0P004f>o00<006>oHkl0:6>o00<006>oHkl0/f>o000SHkl00`00HkmS_`0CHkl00`00HkmS_`0V Hkl00`00HkmS_`2dHkl002=S_`03001S_f>o01=S_`03001S_f>o02ES_`03001S_f>o0;ES_`008f>o 00<006>oHkl04f>o00<006>oHkl09F>o00<006>oHkl0]F>o000SHkl00`00HkmS_`0DHkl00`00HkmS _`0SHkl00`00HkmS_`2fHkl002=S_`03001S_f>o01AS_`03001S_f>o029S_`03001S_f>o0;MS_`00 8f>o0P005V>o00<006>oHkl086>o00<006>oHkl0^6>o000SHkl00`00HkmS_`0FHkl00`00HkmS_`0O Hkl00`00HkmS_`2hHkl002=S_`03001S_f>o01IS_`03001S_f>o01iS_`03001S_f>o0;US_`008f>o 00<006>oHkl05f>o00<006>oHkl076>o00<006>oHkl0^V>o000SHkl00`00HkmS_`0GHkl00`00HkmS _`0KHkl00`00HkmS_`2kHkl002=S_`8001US_`03001S_f>o01YS_`03001S_f>o0;]S_`008f>o00<0 06>oHkl066>o00<006>oHkl06F>o00<006>oHkl0_6>o000SHkl00`00HkmS_`0IHkl00`00HkmS_`0G Hkl00`00HkmS_`2mHkl002=S_`03001S_f>o01US_`03001S_f>o01IS_`03001S_f>o0;iS_`008f>o 00<006>oHkl06V>o00<006>oHkl056>o00<006>oHkl0_f>o000SHkl2000KHkl00`00HkmS_`0CHkl0 0`00HkmS_`30Hkl002=S_`03001S_f>o01]S_`03001S_f>o015S_`03001S_f>o0<5S_`008f>o00<0 06>oHkl076>o0P003f>o0P00a6>o000:Hkl20004Hkl01000HkmS_f>o1@000V>o0`001F>o00<006>o Hkl07V>o00<006>oHkl02f>o00<006>oHkl0a6>o0009Hkl01000HkmS_`0026>o00H006>oHkl006>o 0003Hkl00`00HkmS_`02Hkl00`00HkmS_`0OHkl00`00HkmS_`09Hkl00`00HkmS_`35Hkl000US_`04 001S_f>o0009Hkl00`00HkmS_`05Hkl00`00HkmS_`02Hkl00`00HkmS_`0PHkl30004Hkl40038Hkl0 00US_`04001S_f>o000:Hkl00`00HkmS_`04Hkl00`00HkmS_`02Hkl2000THkl5003;Hkl000US_`04 001S_f>o0007Hkl01P00HkmS_f>o001S_`@000ES_`03001S_f>o0?=S_`002F>o00@006>oHkl000MS _`07001S_f>oHkl006>o000026>o00<006>oHkl0lf>o000:Hkl20009Hkl30002Hkl50004Hkl00`00 HkmS_`3cHkl00?mS_aYS_`00of>o6V>o003oHklJHkl00?mS_aYS_`00of>o6V>o003oHklJHkl00?mS _aYS_`003f>o1@001V>o00<006>oHkl0o6>o000@Hkl01000HkmS_`001V>o00<006>oHkl0o6>o000A Hkl00`00HkmS_`06Hkl2003mHkl0019S_`03001S_f>o00AS_`03001S_`000?eS_`003f>o00D006>o HkmS_`0000AS_`05001S_f>oHkl0003lHkl000mS_`05001S_f>oHkl0003oHkl6Hkl0011S_`<000mS _`<0009S_`<00>mS_`008F>o00D006>oHkmS_`00009S_`03001S_f>o0>iS_`003V>o300000=S_`00 00000P000V>o00D006>oHkmS_`00009S_`03001S_f>o00QS_`800>AS_`008F>o00<006>o00000`00 00=S_`0000002V>o00<006>oHkl0hV>o000QHkl00`00HkmS_`03Hkl200000f>o00000008Hkl00`00 HkmS_`3RHkl0025S_`05001S_f>oHkl00008Hkl300000f>o00000003003RHkl0015S_`800003Hkl0 06>o009S_`800003Hkl0000000ES_`<00?AS_`0046>o00@006>oHkl000ES_`03001S_`0001QS_`03 001S_f>o0>9S_`0046>o00@006>oHkl000IS_`03001S_f>o01MS_`03001S_f>o0>9S_`0046>o00@0 06>oHkl000ES_`03001S_`0001QS_`03001S_f>o0>9S_`003f>o0P0000=S_`00000016>o0P0000=S _`000000o6>o0000\ \>"], ImageRangeCache->{{{0, 280.438}, {172.875, 0}} -> {-0.816041, -0.0222268, \ 0.023172, 0.00187858}}] }, Open ]], Cell["See page 489 of I &D for a comparison.", "Text"], Cell[CellGroupData[{ Cell["\<\ We can plot the temperature to see how linear the profile is. Also \ note where it reaches 0. Do some plots to tell how this scales with Pr. \ \ \>", "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(\(Tpr1 = Plot[y4[x] /. zz\[LeftDoubleBracket]1\[RightDoubleBracket], {x, 0, youter}, PlotRange \[Rule] All, AxesLabel \[Rule] {"\<\[Eta]\>", \ \*"\"\<\!\(\(T\^*\)\)\>\""}];\)\)], "Input", AspectRatioFixed->True], 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 0.190476 0.0147151 0.588604 [ [.21429 .00222 -3 -9 ] [.21429 .00222 3 0 ] [.40476 .00222 -3 -9 ] [.40476 .00222 3 0 ] [.59524 .00222 -3 -9 ] [.59524 .00222 3 0 ] [.78571 .00222 -3 -9 ] [.78571 .00222 3 0 ] [.97619 .00222 -3 -9 ] [.97619 .00222 3 0 ] [1.025 .01472 0 -6 ] [1.025 .01472 10 6 ] [.01131 .13244 -18 -4.5 ] [.01131 .13244 0 4.5 ] [.01131 .25016 -18 -4.5 ] [.01131 .25016 0 4.5 ] [.01131 .36788 -18 -4.5 ] [.01131 .36788 0 4.5 ] [.01131 .4856 -18 -4.5 ] [.01131 .4856 0 4.5 ] [.01131 .60332 -6 -4.5 ] [.01131 .60332 0 4.5 ] [.02381 .64303 -7.625 0 ] [.02381 .64303 7.625 12 ] [ 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 [(1)] .21429 .00222 0 1 Mshowa .40476 .01472 m .40476 .02097 L s [(2)] .40476 .00222 0 1 Mshowa .59524 .01472 m .59524 .02097 L s [(3)] .59524 .00222 0 1 Mshowa .78571 .01472 m .78571 .02097 L s [(4)] .78571 .00222 0 1 Mshowa .97619 .01472 m .97619 .02097 L s [(5)] .97619 .00222 0 1 Mshowa .125 Mabswid .0619 .01472 m .0619 .01847 L s .1 .01472 m .1 .01847 L s .1381 .01472 m .1381 .01847 L s .17619 .01472 m .17619 .01847 L s .25238 .01472 m .25238 .01847 L s .29048 .01472 m .29048 .01847 L s .32857 .01472 m .32857 .01847 L s .36667 .01472 m .36667 .01847 L s .44286 .01472 m .44286 .01847 L s .48095 .01472 m .48095 .01847 L s .51905 .01472 m .51905 .01847 L s .55714 .01472 m .55714 .01847 L s .63333 .01472 m .63333 .01847 L s .67143 .01472 m .67143 .01847 L s .70952 .01472 m .70952 .01847 L s .74762 .01472 m .74762 .01847 L s .82381 .01472 m .82381 .01847 L s .8619 .01472 m .8619 .01847 L s .9 .01472 m .9 .01847 L s .9381 .01472 m .9381 .01847 L s .25 Mabswid 0 .01472 m 1 .01472 L s gsave 1.025 .01472 -61 -10 Mabsadd m 1 1 Mabs scale currentpoint translate /MISOfy { /newfontname exch def /oldfontname exch def oldfontname findfont dup length dict begin {1 index /FID ne {def} {pop pop} ifelse} forall /Encoding ISOLatin1Encoding def currentdict end newfontname exch definefont pop } def 0 20 translate 1 -1 scale 63.000 13.000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor 0.000 0.000 rmoveto 63.000 13.000 moveto %%IncludeResource: font Math1Mono %%IncludeFont: Math1Mono /Math1Mono findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (h) show 69.000 13.000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor 0.000 0.000 rmoveto 1.000 setlinewidth grestore .02381 .13244 m .03006 .13244 L s [(0.2)] .01131 .13244 1 0 Mshowa .02381 .25016 m .03006 .25016 L s [(0.4)] .01131 .25016 1 0 Mshowa .02381 .36788 m .03006 .36788 L s [(0.6)] .01131 .36788 1 0 Mshowa .02381 .4856 m .03006 .4856 L s [(0.8)] .01131 .4856 1 0 Mshowa .02381 .60332 m .03006 .60332 L s [(1)] .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 gsave .02381 .64303 -68.625 -4 Mabsadd m 1 1 Mabs scale currentpoint translate /MISOfy { /newfontname exch def /oldfontname exch def oldfontname findfont dup length dict begin {1 index /FID ne {def} {pop pop} ifelse} forall /Encoding ISOLatin1Encoding def currentdict end newfontname exch definefont pop } def 0 20 translate 1 -1 scale 63.000 13.000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor 0.000 0.000 rmoveto 63.000 13.000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (T) show 69.000 9.000 moveto %%IncludeResource: font Math1Mono %%IncludeFont: Math1Mono /Math1Mono findfont 7.125 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (*) show 74.250 13.000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor 0.000 0.000 rmoveto 1.000 setlinewidth grestore .5 Mabswid .02381 .60332 m .06244 .53569 L .10458 .46259 L .14415 .39593 L .18221 .33525 L .22272 .27599 L .26171 .22536 L .30316 .17916 L .34309 .14215 L .3815 .1131 L .42237 .08849 L .46172 .07004 L .49955 .0563 L .53984 .04514 L .57861 .03706 L .61984 .03067 L .65954 .02614 L .69774 .02291 L .73838 .02038 L .77751 .0186 L .81909 .01722 L .85916 .01626 L .89771 .01559 L .93871 .01506 L .97619 .01472 L s 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath % End of Graphics MathPictureEnd \ \>"], "Graphics", ImageSize->{288, 177.938}, ImageMargins->{{43, 0}, {0, 0}}, ImageRegion->{{0, 1}, {0, 1}}, ImageCache->GraphicsData["Bitmap", "\<\ CF5dJ6E]HGAYHf4PAg9QL6QYHgo8F>o0017Hkl3 000ZHkl5000[Hkl3000/Hkl3000ZHkl3000MHkl004QS_`03001S_f>o02YS_`04001S_f>o000ZHkl0 1@00HkmS_f>o0000;6>o00<006>oHkl0:6>o00D006>oHkmS_`0001aS_`00B6>o00<006>oHkl0:f>o 00<006>oHkl0;V>o00<006>oHkl09f>o1@00;F>o00<006>oHkl06V>o0018Hkl00`00HkmS_`0/Hkl0 0`00HkmS_`0[Hkl2000ZHkl01000HkmS_`00;V>o00<006>oHkl06V>o0018Hkl00`00HkmS_`0YHkl0 1@00HkmS_f>o0000;V>o00<006>oHkl0:6>o00<006>o0000:V>o10007F>o0017Hkl2000[Hkl01@00 HkmS_f>o0000:V>o00D006>oHkmS_`0002]S_`8002YS_`03001S_f>o01iS_`00B6>o00<006>oHkl0 :V>o0`00;6>o0`00;F>o00<006>oHkl0:6>o1@0076>o000KHkl00`00HkmS_`3fHkl00`00HkmS_`09 Hkl001]S_`03001S_f>o0?IS_`03001S_f>o00US_`006f>o00<006>oHkl0lf>o00@006>oHkl000]S _`006f>o00<006>oHkl0lf>o00D006>oHkmS_`0000YS_`005F>om0002F>o00@006>oHkl000YS_`00 6f>o00<006>oHkl01V>o00<006>oHkl01V>o00<006>oHkl01f>o00<006>oHkl01V>o00<006>oHkl0 1V>o00<006>oHkl01V>o00<006>oHkl01f>o00<006>oHkl01V>o00<006>oHkl01V>o00<006>oHkl0 1V>o00<006>oHkl01f>o00<006>oHkl01V>o00<006>oHkl01V>o00<006>oHkl01V>o00<006>oHkl0 1f>o00<006>oHkl01V>o00<006>oHkl01V>o00<006>oHkl01f>o00<006>o00004P001V>o00<006>o Hkl01f>o00<006>oHkl01V>o00<006>oHkl01V>o00<006>oHkl01V>o00<006>oHkl03F>o0P0000=S _`00Hkl02F>o000KHkl00`00HkmS_`2CHklC001o0000000:Hkl001]S_`03001S_f>o 08aS_`L006mS_`006f>o00<006>oHkl0Qf>o1@00MV>o000KHkl00`00HkmS_`22Hkl5001kHkl001]S _`03001S_f>o07eS_`D0081S_`006f>o0P00NF>o1@00QF>o000KHkl00`00HkmS_`1eHkl3002:Hkl0 01]S_`03001S_f>o079S_`<008eS_`006f>o00<006>oHkl0Kf>o0`00T6>o000KHkl00`00HkmS_`1/ Hkl3002CHkl001]S_`03001S_f>o06US_`<009IS_`006f>o00<006>oHkl0IV>o0`00VF>o000KHkl2 001UHkl2002LHkl001]S_`03001S_f>o069S_`8009iS_`006f>o00<006>oHkl0HF>o00<006>oHkl0 WV>o000KHkl00`00HkmS_`1OHkl2002QHkl001]S_`03001S_f>o05eS_`800:=S_`006f>o00<006>o Hkl0Ff>o0P00YF>o000KHkl00`00HkmS_`1JHkl00`00HkmS_`2UHkl001]S_`8005US_`800:QS_`00 6f>o00<006>oHkl0EV>o0P00ZV>o000KHkl00`00HkmS_`1EHkl00`00HkmS_`2ZHkl001]S_`03001S _f>o05=S_`800:eS_`006f>o00<006>oHkl0DV>o00<006>oHkl0[F>o0008Hkl20004Hkl01000HkmS _f>o1@0016>o00<006>oHkl0D6>o0P00/6>o0007Hkl01000HkmS_`0026>o00@006>oHkl000AS_`03 001S_f>o04mS_`03001S_f>o0;1S_`001f>o00@006>oHkl000US_`03001S_f>o00AS_`8004mS_`03 001S_f>o0;5S_`001f>o00@006>oHkl000YS_`03001S_f>o00=S_`03001S_f>o04aS_`800;AS_`00 1f>o00@006>oHkl000MS_`05001S_f>oHkl00004Hkl00`00HkmS_`1;Hkl00`00HkmS_`2dHkl000MS _`04001S_f>o0007Hkl01@00HkmS_f>o000016>o00<006>oHkl0BV>o00<006>oHkl0]F>o0008Hkl2 0009Hkl30005Hkl00`00HkmS_`19Hkl00`00HkmS_`2fHkl001]S_`03001S_f>o04QS_`03001S_f>o 0;MS_`006f>o00<006>oHkl0Af>o00<006>oHkl0^6>o000KHkl20017Hkl00`00HkmS_`2iHkl001]S _`03001S_f>o04ES_`03001S_f>o0;YS_`006f>o00<006>oHkl0A6>o00<006>oHkl0^f>o000KHkl0 0`00HkmS_`13Hkl00`00HkmS_`2lHkl001]S_`03001S_f>o049S_`03001S_f>o0;eS_`006f>o00<0 06>oHkl0@F>o00<006>oHkl0_V>o000KHkl00`00HkmS_`11Hkl00`00HkmS_`2nHkl001]S_`03001S _f>o041S_`03001S_f>o0;mS_`006f>o0P00@6>o00<006>oHkl0`6>o000KHkl00`00HkmS_`0nHkl0 0`00HkmS_`31Hkl001]S_`03001S_f>o03eS_`03001S_f>o0<9S_`006f>o00<006>oHkl0?6>o00<0 06>oHkl0`f>o000KHkl00`00HkmS_`0kHkl00`00HkmS_`34Hkl001]S_`03001S_f>o03YS_`03001S _f>o0o00<006>oHkl0>F>o00<006>oHkl0aV>o000KHkl2000iHkl00`00HkmS_`37Hkl0 01]S_`03001S_f>o03MS_`03001S_f>o0o00<006>oHkl0=f>o00<006>oHkl0b6>o000K Hkl00`00HkmS_`0fHkl00`00HkmS_`39Hkl000QS_`8000AS_`03001S_f>o00=S_`<000AS_`03001S _f>o03ES_`03001S_f>o0o00@006>oHkl000YS_`03001S_f>o00=S_`03001S_f>o03AS _`03001S_f>o0<]S_`001f>o00@006>oHkl000MS_`D000AS_`03001S_f>o03=S_`03001S_f>o0o00@006>oHkl000MS_`04001S_f>o0005Hkl2000cHkl00`00HkmS_`3=Hkl000MS_`04001S _f>o0008Hkl00`00Hkl00005Hkl00`00HkmS_`0bHkl00`00HkmS_`3=Hkl000MS_`04001S_f>o0009 Hkl20005Hkl00`00HkmS_`0aHkl00`00HkmS_`3>Hkl000QS_`8000]S_`03001S_f>o00=S_`03001S _f>o031S_`03001S_f>o0o00<006>oHkl0;f>o00<006>oHkl0d6>o000KHkl00`00HkmS _`0^Hkl00`00HkmS_`3AHkl001]S_`03001S_f>o02eS_`03001S_f>o0=9S_`006f>o0P00;V>o00<0 06>oHkl0dV>o000KHkl00`00HkmS_`0/Hkl00`00HkmS_`3CHkl001]S_`03001S_f>o02]S_`03001S _f>o0=AS_`006f>o00<006>oHkl0:f>o00<006>oHkl0e6>o000KHkl00`00HkmS_`0ZHkl00`00HkmS _`3EHkl001]S_`03001S_f>o02US_`03001S_f>o0=IS_`006f>o00<006>oHkl0:F>o00<006>oHkl0 eV>o000KHkl2000YHkl00`00HkmS_`3GHkl001]S_`03001S_f>o02MS_`03001S_f>o0=QS_`006f>o 00<006>oHkl09f>o00<006>oHkl0f6>o000KHkl00`00HkmS_`0VHkl00`00HkmS_`3IHkl001]S_`03 001S_f>o02ES_`03001S_f>o0=YS_`006f>o00<006>oHkl09F>o00<006>oHkl0fV>o000KHkl00`00 HkmS_`0THkl00`00HkmS_`3KHkl001]S_`8002AS_`03001S_f>o0=aS_`006f>o00<006>oHkl08f>o 00<006>oHkl0g6>o000KHkl00`00HkmS_`0RHkl00`00HkmS_`3MHkl001]S_`03001S_f>o025S_`03 001S_f>o0=iS_`006f>o00<006>oHkl08F>o00<006>oHkl0gV>o0008Hkl20004Hkl00`00HkmS_`02 Hkl30005Hkl00`00HkmS_`0PHkl00`00HkmS_`3OHkl000MS_`04001S_f>o0007Hkl01@00HkmS_f>o 000016>o00<006>oHkl07f>o00<006>oHkl0h6>o0007Hkl01000HkmS_`001f>o00D006>oHkmS_`00 00AS_`03001S_f>o01mS_`03001S_f>o0>1S_`001f>o00@006>oHkl000MS_`@000ES_`8001mS_`03 001S_f>o0>5S_`001f>o00@006>oHkl000MS_`03001S_f>o00IS_`03001S_f>o01iS_`03001S_f>o 0>5S_`001f>o00@006>oHkl000QS_`03001S_f>o00ES_`03001S_f>o01eS_`03001S_f>o0>9S_`00 26>o0P002V>o0P001F>o00<006>oHkl076>o00<006>oHkl0hf>o000KHkl00`00HkmS_`0LHkl00`00 HkmS_`3SHkl001]S_`03001S_f>o01]S_`03001S_f>o0>AS_`006f>o00<006>oHkl06V>o00<006>o Hkl0iF>o000KHkl2000KHkl00`00HkmS_`3UHkl001]S_`03001S_f>o01US_`03001S_f>o0>IS_`00 6f>o00<006>oHkl06F>o00<006>oHkl0iV>o000KHkl00`00HkmS_`0HHkl00`00HkmS_`3WHkl001]S _`03001S_f>o01MS_`03001S_f>o0>QS_`006f>o00<006>oHkl05f>o00<006>oHkl0j6>o000KHkl0 0`00HkmS_`0FHkl00`00HkmS_`3YHkl001]S_`8001MS_`03001S_f>o0>US_`006f>o00<006>oHkl0 5F>o00<006>oHkl0jV>o000KHkl00`00HkmS_`0DHkl00`00HkmS_`3[Hkl001]S_`03001S_f>o01AS _`03001S_f>o0>]S_`006f>o00<006>oHkl04f>o00<006>oHkl0k6>o000KHkl00`00HkmS_`0CHkl0 0`00HkmS_`3/Hkl001]S_`03001S_f>o019S_`03001S_f>o0>eS_`006f>o0P004V>o00<006>oHkl0 kV>o000KHkl00`00HkmS_`0AHkl00`00HkmS_`3^Hkl001]S_`03001S_f>o011S_`03001S_f>o0>mS _`006f>o00<006>oHkl046>o00<006>oHkl0kf>o000KHkl00`00HkmS_`0?Hkl00`00HkmS_`3`Hkl0 00QS_`8000AS_`03001S_f>o009S_`<000ES_`03001S_f>o00mS_`03001S_f>o0?1S_`001f>o00@0 06>oHkl000MS_`05001S_f>oHkl00004Hkl00`00HkmS_`0>Hkl00`00HkmS_`3aHkl000MS_`04001S _f>o0007Hkl01@00HkmS_f>o000016>o0P003V>o00<006>oHkl0lV>o0007Hkl01000HkmS_`0026>o 0`001F>o00<006>oHkl03F>o00<006>oHkl0lV>o0007Hkl01000HkmS_`001f>o00D006>oHkmS_`00 00AS_`03001S_f>o00aS_`03001S_f>o0?=S_`001f>o00@006>oHkl000MS_`05001S_f>oHkl00004 Hkl00`00HkmS_`0o00]S_`03001S _f>o0?AS_`006f>o00<006>oHkl02V>o00<006>oHkl0mF>o000KHkl00`00HkmS_`0:Hkl00`00HkmS _`3eHkl001]S_`8000YS_`03001S_f>o0?IS_`006f>o00<006>oHkl02F>o00<006>oHkl0mV>o000K Hkl00`00HkmS_`08Hkl00`00HkmS_`3gHkl001]S_`03001S_f>o00MS_`03001S_f>o0?QS_`006f>o 00<006>oHkl01f>o00<006>oHkl0n6>o000KHkl00`00HkmS_`06Hkl00`00HkmS_`3iHkl001]S_`03 001S_f>o00IS_`03001S_f>o0?US_`006f>o00<006>oHkl01F>o00<006>oHkl0nV>o000KHkl20006 Hkl00`00HkmS_`3jHkl001]S_`03001S_f>o00AS_`03001S_f>o0?]S_`006f>o00<006>oHkl016>o 00<006>oHkl0nf>o000KHkl00`00HkmS_`03Hkl00`00HkmS_`3lHkl001]S_`03001S_f>o00=S_`03 001S_f>o0?aS_`006f>o00<006>oHkl00V>o00<006>oHkl0oF>o000KHkl01@00HkmS_f>o0000of>o 0F>o000KHkl20002Hkl00`00HkmS_`3nHkl001]S_`04001S_f>o003oHkl2Hkl001]S_`04001S_f>o 003oHkl2Hkl001]S_`03001S_`000?mS_`=S_`004f>o0`001F>o00<006>o0000of>o0f>o000DHkl0 0`00HkmS_`04Hkl2003oHkl4Hkl001AS_`03001S_f>o00AS_`800?mS_`AS_`0056>o00<006>oHkl0 16>o0P00of>o16>o000DHkl00`00HkmS_`04Hkl00`00HkmS_`3oHkl3Hkl001=S_`8000IS_`03001S _f>o0?mS_`=S_`0056>o00<006>oHkl016>o00<006>oHkl0of>o0f>o003oHklQHkl00?mS_b5S_`00 of>o8F>o003oHklQHkl00?mS_b5S_`00of>o8F>o003oHklQHkl00?mS_b5S_`00of>o8F>o000FHkl3 003oHkl8Hkl001MS_`03001S_f>o0?mS_`MS_`005f>o00<006>oHkl0of>o1f>o000GHkl00`00HkmS _`3oHkl7Hkl001ES_`05001S_`00Hkl00002Hkl00`00HkmS_`3oHkl2Hkl001ES_`05001S_`00Hkl0 0003Hkl00`00HkmS_`3oHkl1Hkl001ES_`D00?mS_`MS_`00of>o8F>o003oHklQHkl00001\ \>"], ImageRangeCache->{{{0, 287}, {176.938, 0}} -> {-0.586937, -0.0840689, \ 0.0215746, 0.00698168}}] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[TextData[{ "If we wave the plots with names (the ones that I used have values of ", StyleBox["Pr", FontSlant->"Italic"], " in them), then you can \"Show\" the results together. (I ran this three \ times to make the plots with different names so that I could \"Show\" them.) \ For Pr = 0.1, 1, 10, 100." }], "Text", CellTags->"results"], Cell[CellGroupData[{ Cell[BoxData[ \(\(Show[Upr10, Upr01, \ Upr100, Upr1];\)\)], "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 0.190476 0.0147151 1.51925 [ [.21429 .00222 -3 -9 ] [.21429 .00222 3 0 ] [.40476 .00222 -3 -9 ] [.40476 .00222 3 0 ] [.59524 .00222 -3 -9 ] [.59524 .00222 3 0 ] [.78571 .00222 -3 -9 ] [.78571 .00222 3 0 ] [.97619 .00222 -3 -9 ] [.97619 .00222 3 0 ] [1.025 .01472 0 -6 ] [1.025 .01472 10 6 ] [.01131 .16664 -18 -4.5 ] [.01131 .16664 0 4.5 ] [.01131 .31857 -18 -4.5 ] [.01131 .31857 0 4.5 ] [.01131 .47049 -18 -4.5 ] [.01131 .47049 0 4.5 ] [.02381 .64303 -24.7188 0 ] [.02381 .64303 24.7188 18.6875 ] [ 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 [(1)] .21429 .00222 0 1 Mshowa .40476 .01472 m .40476 .02097 L s [(2)] .40476 .00222 0 1 Mshowa .59524 .01472 m .59524 .02097 L s [(3)] .59524 .00222 0 1 Mshowa .78571 .01472 m .78571 .02097 L s [(4)] .78571 .00222 0 1 Mshowa .97619 .01472 m .97619 .02097 L s [(5)] .97619 .00222 0 1 Mshowa .125 Mabswid .0619 .01472 m .0619 .01847 L s .1 .01472 m .1 .01847 L s .1381 .01472 m .1381 .01847 L s .17619 .01472 m .17619 .01847 L s .25238 .01472 m .25238 .01847 L s .29048 .01472 m .29048 .01847 L s .32857 .01472 m .32857 .01847 L s .36667 .01472 m .36667 .01847 L s .44286 .01472 m .44286 .01847 L s .48095 .01472 m .48095 .01847 L s .51905 .01472 m .51905 .01847 L s .55714 .01472 m .55714 .01847 L s .63333 .01472 m .63333 .01847 L s .67143 .01472 m .67143 .01847 L s .70952 .01472 m .70952 .01847 L s .74762 .01472 m .74762 .01847 L s .82381 .01472 m .82381 .01847 L s .8619 .01472 m .8619 .01847 L s .9 .01472 m .9 .01847 L s .9381 .01472 m .9381 .01847 L s .25 Mabswid 0 .01472 m 1 .01472 L s gsave 1.025 .01472 -61 -10 Mabsadd m 1 1 Mabs scale currentpoint translate /MISOfy { /newfontname exch def /oldfontname exch def oldfontname findfont dup length dict begin {1 index /FID ne {def} {pop pop} ifelse} forall /Encoding ISOLatin1Encoding def currentdict end newfontname exch definefont pop } def 0 20 translate 1 -1 scale 63.000 13.000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor 0.000 0.000 rmoveto 63.000 13.000 moveto %%IncludeResource: font Math1Mono %%IncludeFont: Math1Mono /Math1Mono findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (h) show 69.000 13.000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor 0.000 0.000 rmoveto 1.000 setlinewidth grestore .02381 .16664 m .03006 .16664 L s [(0.1)] .01131 .16664 1 0 Mshowa .02381 .31857 m .03006 .31857 L s [(0.2)] .01131 .31857 1 0 Mshowa .02381 .47049 m .03006 .47049 L s [(0.3)] .01131 .47049 1 0 Mshowa .125 Mabswid .02381 .0451 m .02756 .0451 L s .02381 .07549 m .02756 .07549 L s .02381 .10587 m .02756 .10587 L s .02381 .13626 m .02756 .13626 L s .02381 .19703 m .02756 .19703 L s .02381 .22741 m .02756 .22741 L s .02381 .2578 m .02756 .2578 L s .02381 .28818 m .02756 .28818 L s .02381 .34895 m .02756 .34895 L s .02381 .37933 m .02756 .37933 L s .02381 .40972 m .02756 .40972 L s .02381 .4401 m .02756 .4401 L s .02381 .50087 m .02756 .50087 L s .02381 .53126 m .02756 .53126 L s .02381 .56164 m .02756 .56164 L s .02381 .59203 m .02756 .59203 L s .25 Mabswid .02381 0 m .02381 .61803 L s gsave .02381 .64303 -85.7188 -4 Mabsadd m 1 1 Mabs scale currentpoint translate /MISOfy { /newfontname exch def /oldfontname exch def oldfontname findfont dup length dict begin {1 index /FID ne {def} {pop pop} ifelse} forall /Encoding ISOLatin1Encoding def currentdict end newfontname exch definefont pop } def 0 26.6875 translate 1 -1 scale 63.000 15.438 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor 0.000 0.000 rmoveto 66.062 9.500 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (u) show 74.688 9.500 moveto (x) show %%IncludeResource: font Math1Mono %%IncludeFont: Math1Mono /Math1Mono findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor 64.750 15.438 moveto (\\200\\200\\200\\200) show 73.250 15.438 moveto (\\200\\200) show 77.500 15.438 moveto (\\200) show 79.688 15.438 moveto (\\200) show 66.062 21.688 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (2) show 72.062 21.688 moveto %%IncludeResource: font Math1Mono %%IncludeFont: Math1Mono /Math1Mono findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor ( ) show 74.688 21.688 moveto %%IncludeResource: font Math1Mono %%IncludeFont: Math1Mono /Math1Mono findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (n) show 83.812 15.438 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (Gr) show 95.812 11.438 moveto %%IncludeResource: font Math1Mono %%IncludeFont: Math1Mono /Math1Mono findfont 7.125 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (-) show 102.188 8.438 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 5.062 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (1) show %%IncludeResource: font Math1Mono %%IncludeFont: Math1Mono /Math1Mono findfont 7.125 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor 101.250 11.438 moveto (\\200\\200) show 104.125 11.438 moveto (\\200) show 104.438 11.438 moveto (\\200) show 102.188 13.375 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 5.062 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (2) show 108.438 15.438 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor 0.000 0.000 rmoveto 1.000 setlinewidth grestore .5 Mabswid .02381 .01472 m .04262 .07039 L .06244 .11494 L .08255 .14734 L .09396 .16078 L .10458 .17052 L .11508 .17777 L .12023 .18057 L .125 .18273 L .12924 .18433 L .1339 .18577 L .1386 .1869 L .14357 .18776 L .14606 .18807 L .14737 .18821 L .14877 .18833 L .15011 .18842 L .15133 .18848 L .15255 .18853 L .15371 .18856 L .15484 .18858 L .15603 .18858 L .15705 .18857 L .15816 .18855 L .15937 .18851 L .16066 .18845 L .16298 .1883 L .16567 .18807 L .16859 .18775 L .1739 .18697 L .18386 .18494 L .19326 .18245 L .20363 .17916 L .22464 .17125 L .26485 .15373 L .30353 .13643 L .34467 .11919 L .3843 .10427 L .42241 .09156 L .46297 .07968 L .50201 .0697 L .54351 .06048 L .58349 .05278 L .62196 .04633 L .66288 .04037 L .70228 .03539 L .74017 .03122 L .78051 .02736 L .81933 .02412 L .86061 .02114 L Mistroke .90037 .01866 L .93861 .01658 L .97619 .01479 L Mfstroke .02381 .01472 m .06244 .23533 L .08255 .32682 L .10458 .40985 L .1253 .47257 L .14415 .5176 L .16372 .55319 L .17449 .56829 L .18466 .57985 L .19436 .58854 L .19966 .59239 L .20463 .59543 L .21341 .59954 L .21807 .60109 L .22307 .60229 L .22561 .60272 L .227 .60291 L .22831 .60305 L .22954 .60316 L .23018 .6032 L .23087 .60324 L .23207 .60329 L .23318 .60332 L .23449 .60332 L .23523 .60331 L .2359 .60329 L .2372 .60323 L .23844 .60315 L .23987 .60302 L .2412 .60287 L .24421 .60244 L .24884 .60149 L .25395 .60006 L .2643 .59605 L .27439 .59078 L .28365 .58488 L .3045 .56834 L .34413 .52746 L .38223 .48082 L .4228 .42771 L .46184 .37626 L .50334 .32347 L .54332 .27578 L .58179 .23353 L .62271 .19285 L .66211 .15793 L .7 .12826 L .74034 .10075 L .77916 .0781 L .82044 .05795 L Mistroke .8602 .04216 L .89844 .03016 L .9198 .02475 L .93914 .02063 L .97619 .01472 L Mfstroke .02381 .01472 m .03154 .02756 L .03996 .03897 L .04788 .04743 L .05549 .05372 L .05973 .05651 L .06359 .05864 L .06765 .06049 L .07139 .06185 L .07522 .06294 L .07737 .06342 L .07938 .0638 L .08133 .06408 L .08312 .06429 L .08407 .06438 L .08511 .06445 L .08568 .06449 L .0862 .06452 L .08675 .06454 L .08725 .06456 L .08774 .06457 L .08799 .06458 L .08826 .06458 L .08853 .06459 L .08878 .06459 L .08894 .06459 L .08909 .06459 L .08922 .06459 L .08936 .06459 L .0896 .06459 L .08986 .06459 L .08999 .06459 L .09014 .06459 L .09039 .06458 L .09089 .06458 L .09135 .06456 L .09222 .06453 L .09316 .06449 L .09507 .06436 L .09713 .06417 L .09899 .06394 L .10318 .0633 L .11118 .06156 L .11887 .05935 L .12706 .05653 L .13494 .05346 L .14252 .05025 L .15058 .04662 L .15835 .04296 L .1666 .03895 L Mistroke .17456 .03499 L .1822 .03113 L .19034 .02699 L .19818 .02297 L .20651 .0187 L .21429 .01472 L Mfstroke .02381 .01472 m .04262 .10366 L .06244 .18235 L .08426 .25229 L .10458 .30307 L .12297 .33816 L .14264 .36542 L .15285 .37578 L .16394 .38433 L .16902 .38737 L .17441 .39002 L .17946 .392 L .18408 .39338 L .18639 .39393 L .18891 .39442 L .19128 .39479 L .19349 .39504 L .19482 .39515 L .19606 .39523 L .19718 .39528 L .19839 .39531 L .19902 .39532 L .19972 .39532 L .20044 .39531 L .20111 .3953 L .20232 .39526 L .20362 .39519 L .20484 .3951 L .20597 .39499 L .20853 .39469 L .21079 .39435 L .21319 .3939 L .22192 .39162 L .2268 .38993 L .23129 .38812 L .24139 .38323 L .25967 .37186 L .29929 .33898 L .3374 .30129 L .37796 .25928 L .41701 .22008 L .4585 .18184 L .49849 .14937 L .53695 .12249 L .57787 .09847 L .61727 .07945 L .65516 .06451 L .6955 .0517 L .73433 .04192 L .7756 .03377 L Mistroke .81536 .0277 L .85361 .02318 L .89431 .01951 L .93349 .01685 L .97512 .01476 L .97619 .01472 L Mfstroke 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath % End of Graphics MathPictureEnd \ \>"], "Graphics", ImageSize->{288, 177.938}, ImageMargins->{{0, 0}, {0, 0}}, ImageRegion->{{0, 1}, {0, 1}}, ImageCache->GraphicsData["Bitmap", "\<\ CF5dJ6E]HGAYHf4PAg9QL6QYHgo8F>o001o02QS_`04001S_f>o000XHkl0 1@00HkmS_f>o0000:V>o00<006>oHkl09f>o00D006>oHkmS_`0001iS_`00CF>o00<006>oHkl0:F>o 00<006>oHkl0;6>o00<006>oHkl09F>o1@00;6>o00<006>oHkl076>o001=Hkl00`00HkmS_`0ZHkl0 0`00HkmS_`0YHkl2000XHkl01000HkmS_`00;F>o00<006>oHkl076>o001=Hkl00`00HkmS_`0WHkl0 1@00HkmS_f>o0000;6>o00<006>oHkl09V>o00<006>o0000:F>o10007f>o001o0000:6>o00D006>oHkmS_`0002US_`8002US_`03001S_f>o021S_`00CF>o00<006>oHkl0 :6>o0`00:V>o0`00:f>o00<006>oHkl09f>o1@003f>o00<006>oHkl036>o000RHkl00`00HkmS_`3/ Hkl00`00HkmS_`0o0>US_`04001S_f>o000>Hkl0029S_`03001S_f>o0>US _`05001S_f>oHkl0000=Hkl001aS_nX000US_`04001S_f>o000=Hkl0029S_`03001S_f>o00IS_`03 001S_f>o00IS_`03001S_f>o00ES_`03001S_f>o00IS_`03001S_f>o00=S_`800003Hkl006>o00MS _`03001S_f>o00IS_`03001S_f>o00IS_`03001S_f>o00IS_`03001S_f>o00ES_`03001S_f>o00IS _`03001S_f>o00IS_`03001S_f>o00IS_`03001S_f>o00IS_`03001S_f>o00IS_`03001S_f>o00IS _`03001S_f>o00IS_`03001S_f>o00ES_`03001S_f>o00IS_`03001S_f>o00IS_`03001S_f>o00IS _`03001S_f>o00IS_`03001S_f>o009S_a@0009S_`03001S_f>o00aS_`800003Hkl006>o00aS_`00 8V>o0P009V>o00D006>oHkmS_`0002]S_`03001S_f>o02YS_`03001S_f>o02US_`03001S_f>o00US _a<000=S_`@000MS_`03001S_f>o00YS_`<00003Hkl0000000eS_`008V>o0`008f>o0P00Q6>o4P00 46>o1@00:f>o000RHkl3000QHkl2001mHkl90005Hkl5000DHkl4000`Hkl0029S_`800003Hkl006>o 01eS_`8007IS_`T000YS_`@001IS_`<003AS_`008V>o0`0000=S_`00Hkl06V>o0P00LF>o1`003V>o 1@005f>o0`00=f>o000RHkl300000f>o001S_`0HHkl2001^Hkl5000@Hkl5000JHkl2000jHkl0029S _`<0009S_`03001S_f>o01AS_`8006US_`L0011S_`D001aS_`<003aS_`008V>o10000V>o00<006>o Hkl046>o0`00I6>o1`0056>o0`007f>o0P00?f>o000RHkl01000Hkl000000f>o0P003V>o0P00HV>o 1@0066>o0`007f>o0`00@F>o000RHkl01000Hkl000001F>o0`001f>o1000H6>o10006V>o0`0086>o 0P00A6>o000RHkl01@00Hkl006>o00001f>o1`00Gf>o1@006f>o0`008F>o0P00AV>o000RHkl01@00 Hkl006>o0000JF>o10007V>o0P008F>o0`00B6>o000RHkl00`00Hkl00002001VHkl3000OHkl3000Q Hkl2001;Hkl0029S_`800004Hkl006>o001RHkl3000PHkl2000RHkl2001=Hkl0029S_`06001S_f>o 001S_`00GV>o10008F>o0P008f>o00<006>oHkl0CF>o000RHkl01000HkmS_`000V>o00<006>oHkl0 EV>o1@008f>o0P008f>o0P00D6>o000RHkl01000HkmS_`000V>o00<006>oHkl0Df>o0`009F>o0`00 96>o00<006>oHkl0D6>o000RHkl00`00HkmS_`020002Hkl00`00HkmS_`1?Hkl3000VHkl2000UHkl2 001CHkl0029S_`03001S_f>o0080009S_`03001S_f>o04aS_`<002MS_`8002IS_`03001S_f>o05=S _`008V>o00D006>oHkmS_`00009S_`03001S_f>o04US_`<002QS_`8002IS_`8005IS_`008V>o0P00 0V>o00D006>oHkmS_`0004MS_`<002YS_`03001S_f>o02AS_`8005QS_`008V>o00D006>oHkmS_`00 00=S_`03001S_f>o04=S_`8002]S_`8002IS_`03001S_f>o05QS_`008V>o00@006>oHkmS_`8000=S _`03001S_f>o03mS_`<002]S_`8002IS_`8005]S_`008V>o00@006>oHkmS_`8000=S_`03001S_f>o 03eS_`8002eS_`03001S_f>o02AS_`8005eS_`008V>o00@006>oHkmS_`8000AS_`03001S_f>o03YS _`8002eS_`8002IS_`03001S_f>o05eS_`008V>o00<006>oHkl00V>o00<006>oHkl00V>o00<006>o Hkl0=f>o0`00;V>o00<006>oHkl096>o0P00H6>o000RHkl00`00HkmS_`02Hkl20004Hkl00`00HkmS _`0dHkl2000_Hkl2000VHkl00`00HkmS_`1PHkl0029S_`8000=S_`8000ES_`03001S_f>o035S_`80 031S_`03001S_f>o02ES_`03001S_f>o065S_`008V>o00<006>oHkl00V>o0P001F>o00<006>oHkl0 ;f>o0P00<6>o0P009V>o0P00I6>o000RHkl00`00HkmS_`02Hkl00`00Hkl00005Hkl00`00HkmS_`0[ Hkl3000aHkl00`00HkmS_`0UHkl00`00HkmS_`1THkl0029S_`03001S_f>o009S_`03001S_`0000IS _`03001S_f>o02QS_`80039S_`8002MS_`03001S_f>o06ES_`003f>o0P0016>o00<006>oHkl00V>o 0`001F>o00<006>oHkl00f>o0P001f>o00<006>oHkl09F>o0P00o00<006>oHkl09V>o00<006>o Hkl0IV>o000>Hkl01000HkmS_`002F>o00<006>oHkl016>o00<006>oHkl00f>o0P0026>o00<006>o Hkl08V>o0P00=6>o00<006>oHkl09F>o0P00JF>o000>Hkl01000HkmS_`002F>o00<006>oHkl016>o 00<006>oHkl00f>o00<006>o000026>o00<006>oHkl07V>o0`00=6>o0P009f>o00<006>oHkl0JF>o 000>Hkl01000HkmS_`002F>o00<006>oHkl016>o0P0016>o00<006>o00002F>o00<006>oHkl06f>o 0P00=V>o00<006>oHkl09V>o00<006>oHkl0JV>o000>Hkl01000HkmS_`002F>o00<006>oHkl016>o 00<006>oHkl00f>o00<006>o00002V>o0P006F>o0P00=f>o00<006>oHkl09V>o00<006>oHkl0Jf>o 000>Hkl01000HkmS_`0026>o0P001V>o00<006>oHkl00f>o00<006>o000036>o00<006>oHkl04f>o 0`00>6>o00<006>oHkl09V>o00<006>oHkl0K6>o000?Hkl2000:Hkl00`00HkmS_`04Hkl00`00HkmS _`04Hkl00`00Hkl0000o00AS_`03001S_`0000eS_`@000US_`D003YS_`03001S_f>o02MS_`03001S_f>o06iS_`00 8V>o00<006>oHkl016>o00<006>o00004F>o2@00?V>o00<006>oHkl09f>o00<006>oHkl0Kf>o000R Hkl00`00HkmS_`04Hkl00`00Hkl0001GHkl00`00HkmS_`0WHkl00`00HkmS_`1`Hkl0029S_`8000ES _`04001S_f>o001EHkl00`00HkmS_`0WHkl00`00HkmS_`1aHkl0029S_`03001S_f>o00ES_`03001S _`0005AS_`03001S_f>o02MS_`03001S_f>o079S_`008V>o00<006>oHkl01F>o00<006>o0000Df>o 00<006>oHkl09f>o00<006>oHkl0Lf>o000RHkl00`00HkmS_`05Hkl01000HkmS_`00DF>o00<006>o Hkl09f>o00<006>oHkl0M6>o000RHkl00`00HkmS_`05Hkl01000HkmS_`00D6>o00<006>oHkl09f>o 00<006>oHkl0MF>o000RHkl00`00HkmS_`05Hkl01000HkmS_`00Cf>o00<006>oHkl09f>o00<006>o Hkl0MV>o000RHkl00`00HkmS_`05Hkl01@00HkmS_f>o0000CF>o00<006>oHkl09f>o00<006>oHkl0 Mf>o000RHkl20007Hkl01000HkmS_`00C6>o00<006>oHkl09f>o00<006>oHkl0N6>o000RHkl00`00 HkmS_`06Hkl01000HkmS_`00Bf>o00<006>oHkl09f>o00<006>oHkl0NF>o000RHkl00`00HkmS_`06 Hkl01@00HkmS_f>o0000BF>o00<006>oHkl0:6>o00<006>oHkl0NF>o000RHkl00`00HkmS_`06Hkl0 1@00HkmS_f>o0000B6>o00<006>oHkl0:6>o00<006>oHkl0NV>o000RHkl00`00HkmS_`06Hkl01@00 HkmS_f>o0000Af>o00<006>oHkl0:6>o00<006>oHkl0Nf>o000RHkl00`00HkmS_`06Hkl00`00HkmS _`02Hkl00`00HkmS_`13Hkl00`00HkmS_`0XHkl00`00HkmS_`1lHkl0029S_`03001S_f>o00MS_`05 001S_f>oHkl00014Hkl00`00HkmS_`0XHkl00`00HkmS_`1mHkl0029S_`8000QS_`05001S_f>oHkl0 0013Hkl00`00HkmS_`0XHkl00`00HkmS_`1nHkl0029S_`03001S_f>o00MS_`03001S_f>o009S_`03 001S_f>o041S_`03001S_f>o02MS_`03001S_f>o07mS_`008V>o00<006>oHkl01f>o00<006>oHkl0 0V>o00<006>oHkl0?f>o00<006>oHkl09f>o00<006>oHkl0P6>o000RHkl00`00HkmS_`07Hkl00`00 HkmS_`02Hkl00`00HkmS_`0nHkl00`00HkmS_`0WHkl00`00HkmS_`21Hkl0029S_`03001S_f>o00QS _`03001S_f>o009S_`03001S_f>o03aS_`03001S_f>o02MS_`03001S_f>o089S_`008V>o00<006>o Hkl026>o00<006>oHkl00V>o00<006>oHkl0>f>o00<006>oHkl0:6>o00<006>oHkl0PV>o000RHkl0 0`00HkmS_`08Hkl00`00HkmS_`03Hkl00`00HkmS_`0iHkl00`00HkmS_`0XHkl00`00HkmS_`23Hkl0 029S_`8000US_`03001S_f>o00=S_`03001S_f>o03QS_`03001S_f>o02QS_`03001S_f>o08AS_`00 8V>o00<006>oHkl026>o00<006>oHkl016>o00<006>oHkl0=V>o00<006>oHkl0:6>o00<006>oHkl0 QF>o000RHkl00`00HkmS_`09Hkl00`00HkmS_`03Hkl00`00HkmS_`0eHkl00`00HkmS_`0XHkl00`00 HkmS_`26Hkl0029S_`03001S_f>o00US_`03001S_f>o00=S_`03001S_f>o03AS_`03001S_f>o02QS _`03001S_f>o08MS_`008V>o00<006>oHkl02F>o00<006>oHkl016>o00<006>oHkl0o00<006>o Hkl0:6>o00<006>oHkl0R6>o000?Hkl20004Hkl01000HkmS_f>o1@0016>o00<006>oHkl02F>o00<0 06>oHkl016>o00<006>oHkl0o00<006>oHkl0:6>o00<006>oHkl0RF>o000>Hkl01000HkmS_`00 26>o00@006>oHkl000AS_`03001S_f>o00US_`03001S_f>o00ES_`03001S_f>o02mS_`03001S_f>o 02QS_`03001S_f>o08YS_`003V>o00@006>oHkl000US_`03001S_f>o00AS_`8000]S_`03001S_f>o 00AS_`03001S_f>o02iS_`03001S_f>o02QS_`03001S_f>o08]S_`003V>o00@006>oHkl000YS_`03 001S_f>o00=S_`03001S_f>o00YS_`03001S_f>o00ES_`03001S_f>o02aS_`03001S_f>o02QS_`03 001S_f>o08aS_`003V>o00@006>oHkl000MS_`05001S_f>oHkl00004Hkl00`00HkmS_`0:Hkl00`00 HkmS_`05Hkl00`00HkmS_`0[Hkl00`00HkmS_`0YHkl00`00HkmS_`2o0007 Hkl01@00HkmS_f>o000016>o00<006>oHkl02V>o00<006>oHkl01V>o00<006>oHkl0:F>o00<006>o Hkl0:F>o00<006>oHkl0SF>o000?Hkl20009Hkl30005Hkl00`00HkmS_`0:Hkl00`00HkmS_`06Hkl0 0`00HkmS_`0XHkl00`00HkmS_`0YHkl00`00HkmS_`2>Hkl0029S_`03001S_f>o00]S_`03001S_f>o 00IS_`03001S_f>o02ES_`8002aS_`03001S_f>o08iS_`008V>o00<006>oHkl02f>o00<006>oHkl0 1V>o00<006>oHkl096>o00<006>oHkl0:f>o00<006>oHkl0Sf>o000RHkl2000o00aS_`03001S _f>o00MS_`03001S_f>o021S_`03001S_f>o02]S_`03001S_f>o095S_`008V>o00<006>oHkl036>o 00<006>oHkl01f>o00<006>oHkl07f>o00<006>oHkl0;6>o00<006>oHkl0TF>o000RHkl00`00HkmS _`0o00eS_`03001S_f>o00QS_`03001S_f>o01]S_`03001S_f>o02aS_`03001S_f>o09=S_`00 8V>o00<006>oHkl03F>o00<006>oHkl02F>o00<006>oHkl06F>o00<006>oHkl0;F>o00<006>oHkl0 Tf>o000RHkl00`00HkmS_`0=Hkl00`00HkmS_`0:Hkl00`00HkmS_`0FHkl2000_Hkl00`00HkmS_`2D Hkl0029S_`8000iS_`03001S_f>o00]S_`03001S_f>o019S_`<0031S_`03001S_f>o09ES_`008V>o 00<006>oHkl03V>o00<006>oHkl02f>o00<006>oHkl03f>o0P00o00<006>oHkl0UF>o000RHkl0 0`00HkmS_`0>Hkl00`00HkmS_`0o00iS_`03001S_f>o00mS_`8000QS_`<003AS_`03001S_f>o09MS_`008V>o00<006>oHkl0 3f>o00<006>oHkl046>o2000=V>o00<006>oHkl0V6>o000RHkl00`00HkmS_`0?Hkl00`00HkmS_`1> Hkl00`00HkmS_`2HHkl0029S_`03001S_f>o00mS_`03001S_f>o04eS_`03001S_f>o09US_`008V>o 0P004F>o00<006>oHkl0Bf>o00<006>oHkl0VV>o000RHkl00`00HkmS_`0@Hkl00`00HkmS_`1:Hkl0 0`00HkmS_`2KHkl0029S_`03001S_f>o011S_`03001S_f>o04YS_`03001S_f>o09]S_`008V>o00<0 06>oHkl046>o00<006>oHkl0BF>o00<006>oHkl0W6>o000RHkl00`00HkmS_`0AHkl00`00HkmS_`17 Hkl00`00HkmS_`2MHkl0029S_`03001S_f>o015S_`03001S_f>o04IS_`03001S_f>o09iS_`008V>o 00<006>oHkl04F>o00<006>oHkl0AV>o00<006>oHkl0WV>o000RHkl00`00HkmS_`0BHkl00`00HkmS _`14Hkl00`00HkmS_`2OHkl0029S_`8001=S_`03001S_f>o04=S_`03001S_f>o0:1S_`008V>o00<0 06>oHkl04V>o00<006>oHkl0@V>o00<006>oHkl0XF>o000RHkl00`00HkmS_`0BHkl00`00HkmS_`11 Hkl00`00HkmS_`2RHkl0029S_`03001S_f>o01=S_`03001S_f>o03mS_`03001S_f>o0:=S_`003f>o 0P0016>o00<006>oHkl00V>o0`001F>o00<006>oHkl04f>o00<006>oHkl0?f>o00<006>oHkl0Xf>o 000>Hkl01000HkmS_`001f>o00D006>oHkmS_`0000AS_`03001S_f>o01=S_`03001S_f>o03iS_`03 001S_f>o0:AS_`003V>o00@006>oHkl000]S_`03001S_f>o009S_`03001S_f>o01AS_`03001S_f>o 03aS_`03001S_f>o0:ES_`003V>o00@006>oHkl000US_`8000ES_`8001ES_`03001S_f>o03]S_`03 001S_f>o0:IS_`003V>o00@006>oHkl000]S_`03001S_f>o009S_`03001S_f>o01AS_`03001S_f>o 03YS_`03001S_f>o0:MS_`003V>o00@006>oHkl000MS_`05001S_f>oHkl00004Hkl00`00HkmS_`0E Hkl00`00HkmS_`0hHkl00`00HkmS_`2XHkl000mS_`8000US_`<000ES_`03001S_f>o01ES_`03001S _f>o03QS_`03001S_f>o0:QS_`008V>o00<006>oHkl05V>o00<006>oHkl0=V>o00<006>oHkl0ZF>o 000RHkl00`00HkmS_`0FHkl00`00HkmS_`0eHkl00`00HkmS_`2ZHkl0029S_`03001S_f>o01MS_`03 001S_f>o03=S_`03001S_f>o0:]S_`008V>o0P0066>o00<006>oHkl0o00<006>oHkl0[6>o000R Hkl00`00HkmS_`0GHkl00`00HkmS_`0aHkl00`00HkmS_`2]Hkl0029S_`03001S_f>o01QS_`03001S _f>o031S_`03001S_f>o0:eS_`008V>o00<006>oHkl066>o00<006>oHkl0;f>o00<006>oHkl0[V>o 000RHkl00`00HkmS_`0IHkl00`00HkmS_`0]Hkl00`00HkmS_`2_Hkl0029S_`03001S_f>o01US_`03 001S_f>o02aS_`03001S_f>o0;1S_`008V>o00<006>oHkl06V>o00<006>oHkl0:V>o00<006>oHkl0 /F>o000RHkl2000KHkl00`00HkmS_`0YHkl00`00HkmS_`2bHkl0029S_`03001S_f>o01]S_`03001S _f>o02MS_`03001S_f>o0;=S_`008V>o00<006>oHkl06f>o00<006>oHkl09V>o00<006>oHkl0]6>o 000RHkl00`00HkmS_`0LHkl00`00HkmS_`0THkl00`00HkmS_`2eHkl0029S_`03001S_f>o01aS_`03 001S_f>o02=S_`03001S_f>o0;IS_`008V>o00<006>oHkl07F>o00<006>oHkl08F>o00<006>oHkl0 ]f>o000RHkl00`00HkmS_`0NHkl00`00HkmS_`0OHkl00`00HkmS_`2hHkl0029S_`80021S_`03001S _f>o01eS_`03001S_f>o0;US_`008V>o00<006>oHkl086>o00<006>oHkl06V>o0P00_6>o000RHkl0 0`00HkmS_`0PHkl00`00HkmS_`0IHkl00`00HkmS_`2lHkl0029S_`03001S_f>o025S_`03001S_f>o 01MS_`03001S_f>o0;eS_`008V>o00<006>oHkl08V>o00<006>oHkl05F>o00<006>oHkl0_V>o000R Hkl00`00HkmS_`0SHkl00`00HkmS_`0CHkl00`00HkmS_`2oHkl0029S_`03001S_f>o02AS_`80019S _`03001S_f>o0<1S_`008V>o0P009f>o00<006>oHkl03F>o0P00`f>o000RHkl00`00HkmS_`0WHkl2 000;Hkl20035Hkl0029S_`03001S_f>o02US_`@0009S_`D00o00<006>oHkl0;F>o0`00 bf>o000RHkl00`00HkmS_`3kHkl0029S_`03001S_f>o0?]S_`008V>o00<006>oHkl0nf>o003oHklQ Hkl00?mS_b5S_`00of>o8F>o003oHklQHkl00?mS_b5S_`00of>o8F>o003oHklQHkl000iS_`D000IS _`03001S_f>o0?mS_`ES_`003f>o00@006>oHkl000IS_`03001S_f>o0?mS_`ES_`0046>o00<006>o Hkl01V>o0P00of>o1V>o000AHkl00`00HkmS_`04Hkl00`00Hkl0003oHkl6Hkl000iS_`05001S_f>o Hkl00004Hkl01@00HkmS_f>o0000of>o1F>o000>Hkl01@00HkmS_f>o0000of>o3V>o000?Hkl3000? Hkl30002Hkl3003gHkl0021S_`05001S_f>oHkl00002Hkl00`00HkmS_`3fHkl000eS_``00003Hkl0 00000080009S_`05001S_f>oHkl00002Hkl00`00HkmS_`08Hkl2003/Hkl0021S_`03001S_`0000<0 0003Hkl0000000YS_`03001S_f>o0>YS_`0086>o00<006>oHkl00f>o0P0000=S_`00000026>o00<0 06>oHkl0jV>o000PHkl01@00HkmS_f>o000026>o0`0000=S_`0000001000jF>o000@Hkl200000f>o 001S_`02Hkl200000f>o00000005Hkl3003lHkl000mS_`04001S_f>o0005Hkl00`00Hkl0000HHkl0 0`00HkmS_`3ZHkl000mS_`04001S_f>o0006Hkl00`00HkmS_`0GHkl00`00HkmS_`3ZHkl000mS_`04 001S_f>o0005Hkl00`00Hkl0000HHkl00`00HkmS_`3ZHkl000iS_`800003Hkl0000000AS_`800003 Hkl000000?mS_`ES_`00\ \>"], ImageRangeCache->{{{0, 287}, {176.938, 0}} -> {-0.772438, -0.0336617, \ 0.0225413, 0.00282612}}] }, Open ]], Cell[BoxData[ ButtonBox[\(return\ to\ conclusions\), ButtonData:>"conclusions", ButtonStyle->"Hyperlink"]], "Input"], Cell[CellGroupData[{ Cell[BoxData[ \(\(Show[Tpr10, Tpr01, Tpr100, Tpr1];\)\)], "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 0.190476 0.0147182 0.588601 [ [.21429 .00222 -3 -9 ] [.21429 .00222 3 0 ] [.40476 .00222 -3 -9 ] [.40476 .00222 3 0 ] [.59524 .00222 -3 -9 ] [.59524 .00222 3 0 ] [.78571 .00222 -3 -9 ] [.78571 .00222 3 0 ] [.97619 .00222 -3 -9 ] [.97619 .00222 3 0 ] [1.025 .01472 0 -6 ] [1.025 .01472 10 6 ] [.01131 .13244 -18 -4.5 ] [.01131 .13244 0 4.5 ] [.01131 .25016 -18 -4.5 ] [.01131 .25016 0 4.5 ] [.01131 .36788 -18 -4.5 ] [.01131 .36788 0 4.5 ] [.01131 .4856 -18 -4.5 ] [.01131 .4856 0 4.5 ] [.01131 .60332 -6 -4.5 ] [.01131 .60332 0 4.5 ] [.02381 .64303 -7.625 0 ] [.02381 .64303 7.625 12 ] [ 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 [(1)] .21429 .00222 0 1 Mshowa .40476 .01472 m .40476 .02097 L s [(2)] .40476 .00222 0 1 Mshowa .59524 .01472 m .59524 .02097 L s [(3)] .59524 .00222 0 1 Mshowa .78571 .01472 m .78571 .02097 L s [(4)] .78571 .00222 0 1 Mshowa .97619 .01472 m .97619 .02097 L s [(5)] .97619 .00222 0 1 Mshowa .125 Mabswid .0619 .01472 m .0619 .01847 L s .1 .01472 m .1 .01847 L s .1381 .01472 m .1381 .01847 L s .17619 .01472 m .17619 .01847 L s .25238 .01472 m .25238 .01847 L s .29048 .01472 m .29048 .01847 L s .32857 .01472 m .32857 .01847 L s .36667 .01472 m .36667 .01847 L s .44286 .01472 m .44286 .01847 L s .48095 .01472 m .48095 .01847 L s .51905 .01472 m .51905 .01847 L s .55714 .01472 m .55714 .01847 L s .63333 .01472 m .63333 .01847 L s .67143 .01472 m .67143 .01847 L s .70952 .01472 m .70952 .01847 L s .74762 .01472 m .74762 .01847 L s .82381 .01472 m .82381 .01847 L s .8619 .01472 m .8619 .01847 L s .9 .01472 m .9 .01847 L s .9381 .01472 m .9381 .01847 L s .25 Mabswid 0 .01472 m 1 .01472 L s gsave 1.025 .01472 -61 -10 Mabsadd m 1 1 Mabs scale currentpoint translate /MISOfy { /newfontname exch def /oldfontname exch def oldfontname findfont dup length dict begin {1 index /FID ne {def} {pop pop} ifelse} forall /Encoding ISOLatin1Encoding def currentdict end newfontname exch definefont pop } def 0 20 translate 1 -1 scale 63.000 13.000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor 0.000 0.000 rmoveto 63.000 13.000 moveto %%IncludeResource: font Math1Mono %%IncludeFont: Math1Mono /Math1Mono findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (h) show 69.000 13.000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor 0.000 0.000 rmoveto 1.000 setlinewidth grestore .02381 .13244 m .03006 .13244 L s [(0.2)] .01131 .13244 1 0 Mshowa .02381 .25016 m .03006 .25016 L s [(0.4)] .01131 .25016 1 0 Mshowa .02381 .36788 m .03006 .36788 L s [(0.6)] .01131 .36788 1 0 Mshowa .02381 .4856 m .03006 .4856 L s [(0.8)] .01131 .4856 1 0 Mshowa .02381 .60332 m .03006 .60332 L s [(1)] .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 gsave .02381 .64303 -68.625 -4 Mabsadd m 1 1 Mabs scale currentpoint translate /MISOfy { /newfontname exch def /oldfontname exch def oldfontname findfont dup length dict begin {1 index /FID ne {def} {pop pop} ifelse} forall /Encoding ISOLatin1Encoding def currentdict end newfontname exch definefont pop } def 0 20 translate 1 -1 scale 63.000 13.000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor 0.000 0.000 rmoveto 63.000 13.000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (T) show 69.000 9.000 moveto %%IncludeResource: font Math1Mono %%IncludeFont: Math1Mono /Math1Mono findfont 7.125 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (*) show 74.250 13.000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor 0.000 0.000 rmoveto 1.000 setlinewidth grestore .5 Mabswid .02381 .60332 m .06244 .46446 L .10458 .321 L .14415 .20716 L .18221 .12579 L .20178 .0952 L .22272 .07011 L .24402 .05148 L .2538 .04488 L .26416 .03901 L .27367 .03453 L .28414 .03045 L .30258 .02504 L .31271 .02284 L .32199 .02121 L .33281 .01969 L .34286 .01857 L .35282 .0177 L .36185 .01708 L .37188 .01653 L .38257 .01607 L .3926 .01575 L .40182 .01551 L .41194 .01532 L .42264 .01516 L .43212 .01505 L .44228 .01497 L .45187 .0149 L .46061 .01486 L .47052 .01482 L .4797 .0148 L .48976 .01478 L .50046 .01476 L .50983 .01475 L .51993 .01474 L .52946 .01473 L .53821 .01473 L .54801 .01472 L .55713 .01472 L .56714 .01472 L .57783 .01472 L .58807 .01472 L .59347 .01472 L .59918 .01472 L .60966 .01472 L .61933 .01472 L .62877 .01472 L .6389 .01472 L .64847 .01472 L .65722 .01472 L Mistroke .66709 .01472 L .67624 .01472 L .68628 .01472 L .69186 .01472 L .69698 .01472 L .70729 .01472 L .71271 .01472 L .71843 .01472 L .72806 .01472 L .73352 .01472 L .73861 .01472 L .74909 .01472 L .75424 .01472 L .759 .01472 L .7679 .01472 L .77256 .01472 L .77757 .01472 L .78722 .01472 L .79627 .01472 L .81689 .01472 L .83658 .01472 L .84698 .01472 L .84998 .01472 L .8528 .01472 L .85537 .01472 L .85808 .01472 L .86053 .01472 L .86183 .01472 L .86322 .01472 L .8645 .01472 L .86567 .01472 L .86685 .01472 L .86796 .01472 L .86925 .01472 L .8699 .01472 L .87062 .01472 L .87188 .01472 L .87308 .01472 L .87447 .01472 L .87576 .01472 L .87717 .01472 L .87865 .01472 L .88093 .01472 L .88337 .01472 L .88777 .01472 L .89046 .01472 L .893 .01472 L .89776 .01472 L .90238 .01472 L .90746 .01472 L Mistroke .91775 .01472 L .93619 .01472 L .97619 .01472 L Mfstroke .02381 .60332 m .06244 .57142 L .10458 .53667 L .14415 .50417 L .18221 .47315 L .22272 .44052 L .26171 .40965 L .30316 .37756 L .34309 .34749 L .3815 .31946 L .42237 .29069 L .46172 .26408 L .49955 .23953 L .53984 .21455 L .57861 .19164 L .61984 .16849 L .65954 .14734 L .69774 .12804 L .73838 .10858 L .77751 .09087 L .81909 .07309 L .85916 .05693 L .89771 .04223 L .93871 .02747 L .97619 .01472 L s .02381 .60332 m .03154 .55507 L .03996 .5026 L .04788 .45378 L .05549 .4077 L .06359 .36012 L .07139 .31635 L .07968 .27262 L .08766 .2337 L .09535 .19961 L .10352 .16716 L .11139 .13973 L .11896 .11681 L .12702 .09596 L .13477 .0791 L .14301 .06428 L .15096 .05265 L .1586 .04358 L .16672 .03584 L .17455 .02994 L .18287 .02503 L .19088 .02136 L .19859 .0186 L .20679 .01633 L .21429 .01472 L s .02381 .60332 m .06244 .53569 L .10458 .46259 L .14415 .39593 L .18221 .33525 L .22272 .276 L .26171 .22536 L .30316 .17916 L .34309 .14215 L .3815 .1131 L .42237 .08849 L .46172 .07004 L .49955 .0563 L .53984 .04514 L .57861 .03707 L .61984 .03068 L .65954 .02614 L .69774 .02292 L .73838 .02038 L .77751 .0186 L .81909 .01723 L .85916 .01627 L .89771 .01559 L .93871 .01507 L .97619 .01472 L s 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath % End of Graphics MathPictureEnd \ \>"], "Graphics", ImageSize->{288, 177.938}, ImageMargins->{{0, 0}, {0, 0}}, ImageRegion->{{0, 1}, {0, 1}}, ImageCache->GraphicsData["Bitmap", "\<\ CF5dJ6E]HGAYHf4PAg9QL6QYHgo8F>o0017Hkl3 000ZHkl5000[Hkl3000/Hkl3000ZHkl3000MHkl004QS_`03001S_f>o02YS_`04001S_f>o000ZHkl0 1@00HkmS_f>o0000;6>o00<006>oHkl0:6>o00D006>oHkmS_`0001aS_`00B6>o00<006>oHkl0:f>o 00<006>oHkl0;V>o00<006>oHkl09f>o1@00;F>o00<006>oHkl06V>o0018Hkl00`00HkmS_`0/Hkl0 0`00HkmS_`0[Hkl2000ZHkl01000HkmS_`00;V>o00<006>oHkl06V>o0018Hkl00`00HkmS_`0YHkl0 1@00HkmS_f>o0000;V>o00<006>oHkl0:6>o00<006>o0000:V>o10007F>o0017Hkl2000[Hkl01@00 HkmS_f>o0000:V>o00D006>oHkmS_`0002]S_`8002YS_`03001S_f>o01iS_`00B6>o00<006>oHkl0 :V>o0`00;6>o0`00;F>o00<006>oHkl0:6>o1@0076>o000KHkl00`00HkmS_`3fHkl00`00HkmS_`09 Hkl001]S_`03001S_f>o0?IS_`03001S_f>o00US_`006f>o00<006>oHkl0lf>o00@006>oHkl000]S _`006f>o00<006>oHkl0lf>o00D006>oHkmS_`0000YS_`005F>om0002F>o00@006>oHkl000YS_`00 6f>o00<006>oHkl01V>o00<006>oHkl01V>o00<006>oHkl01f>o00<006>oHkl01V>o00<006>o0000 0`000f>o00<006>oHkl01V>o00<006>oHkl01f>o00<006>oHkl00f>o200016>o00<006>oHkl01V>o 00<006>oHkl01f>o00<006>oHkl01V>o00<006>oHkl01V>o00<006>oHkl01V>o00<006>oHkl01f>o 00<006>oHkl01V>o00<006>oHkl01V>o00<006>oHkl01f>o00<006>o00004P001V>o00<006>oHkl0 1f>o00<006>oHkl01V>o00<006>oHkl01V>o00<006>oHkl00V>o0`0000=S_`00Hkl03V>o0P0000=S _`00Hkl02F>o000KHkl00`00HkmS_`0RHkl2000KHkl5001?HklC000gHkl3000BHkl300000f>o0000 000:Hkl001]S_`03001S_f>o025S_`03001S_f>o01MS_`@004eS_`L004IS_`@002ES_`006f>o00<0 06>oHkl086>o00<006>oHkl05V>o0P00C6>o1@00BV>o0`00:F>o000KHkl00`00HkmS_`0OHkl00`00 HkmS_`0FHkl00`00HkmS_`17Hkl5001o01mS_`03001S_f>o01AS _`8004ES_`D004iS_`<002mS_`006f>o0P007f>o00<006>oHkl04f>o0P00@V>o1@00DF>o0P00o 000KHkl00`00HkmS_`0MHkl00`00HkmS_`0CHkl00`00HkmS_`0oHkl3001CHkl3000dHkl001]S_`03 001S_f>o01aS_`03001S_f>o019S_`8003mS_`<005AS_`8003MS_`006f>o00<006>oHkl076>o00<0 06>oHkl04F>o00<006>oHkl0?6>o0`00EF>o0P00>F>o000KHkl00`00HkmS_`0KHkl00`00HkmS_`0A Hkl00`00HkmS_`0jHkl3001EHkl3000kHkl001]S_`03001S_f>o01YS_`03001S_f>o015S_`03001S _f>o03QS_`<005IS_`8003iS_`006f>o00<006>oHkl06V>o00<006>oHkl046>o00<006>oHkl0=V>o 0`00EV>o0`00@6>o000KHkl2000JHkl00`00HkmS_`0AHkl00`00HkmS_`0dHkl2001GHkl20013Hkl0 01]S_`03001S_f>o01QS_`03001S_f>o015S_`03001S_f>o03=S_`8005IS_`<004ES_`006f>o00<0 06>oHkl066>o00<006>oHkl046>o00<006>oHkl0o00<006>oHkl0E6>o0P00B6>o000KHkl00`00 HkmS_`0GHkl00`00HkmS_`0@Hkl00`00HkmS_`0bHkl2001DHkl3001:Hkl001]S_`03001S_f>o01MS _`03001S_f>o00mS_`03001S_f>o035S_`8005AS_`8004eS_`006f>o00<006>oHkl05V>o00<006>o Hkl03f>o00<006>oHkl0<6>o0P00E6>o0P00Cf>o000KHkl00`00HkmS_`0FHkl00`00HkmS_`0?Hkl0 0`00HkmS_`0_Hkl00`00HkmS_`1AHkl3001AHkl001]S_`8001IS_`03001S_f>o00mS_`03001S_f>o 02iS_`80059S_`8005AS_`006f>o00<006>oHkl05F>o00<006>oHkl03V>o00<006>oHkl0;F>o0P00 DV>o0P00EV>o000KHkl00`00HkmS_`0EHkl00`00HkmS_`0>Hkl00`00HkmS_`0/Hkl00`00HkmS_`1@ Hkl2001HHkl001]S_`03001S_f>o01AS_`03001S_f>o00iS_`03001S_f>o02]S_`80055S_`8005YS _`006f>o00<006>oHkl056>o00<006>oHkl03V>o00<006>oHkl0:V>o00<006>oHkl0Cf>o0P00G6>o 0008Hkl20004Hkl01000HkmS_f>o1@0016>o00<006>oHkl04f>o00<006>oHkl03V>o00<006>oHkl0 :F>o0P00D6>o0P00GV>o0007Hkl01000HkmS_`0026>o00@006>oHkl000AS_`03001S_f>o01=S_`03 001S_f>o00eS_`03001S_f>o02US_`03001S_f>o04iS_`80061S_`001f>o00@006>oHkl000US_`03 001S_f>o00AS_`8001AS_`03001S_f>o00eS_`03001S_f>o02QS_`03001S_f>o04eS_`80069S_`00 1f>o00@006>oHkl000YS_`03001S_f>o00=S_`03001S_f>o019S_`03001S_f>o00eS_`03001S_f>o 02MS_`8004iS_`8006AS_`001f>o00@006>oHkl000MS_`05001S_f>oHkl00004Hkl00`00HkmS_`0B Hkl00`00HkmS_`0=Hkl00`00HkmS_`0VHkl00`00HkmS_`1o0007 Hkl01@00HkmS_f>o000016>o00<006>oHkl04V>o00<006>oHkl036>o00<006>oHkl09V>o00<006>o Hkl0Bf>o0P00J6>o0008Hkl20009Hkl30005Hkl00`00HkmS_`0BHkl00`00HkmS_`0o015S_`03001S_f>o00aS_`03001S_f>o 02ES_`03001S_f>o04US_`8006aS_`006f>o00<006>oHkl04F>o00<006>oHkl036>o00<006>oHkl0 96>o00<006>oHkl0B6>o0P00KV>o000KHkl2000BHkl00`00HkmS_`0;Hkl00`00HkmS_`0THkl00`00 HkmS_`18Hkl00`00HkmS_`1^Hkl001]S_`03001S_f>o011S_`03001S_f>o00aS_`03001S_f>o02=S _`03001S_f>o04MS_`80075S_`006f>o00<006>oHkl046>o00<006>oHkl02f>o00<006>oHkl08f>o 00<006>oHkl0AV>o0P00Lf>o000KHkl00`00HkmS_`0@Hkl00`00HkmS_`0;Hkl00`00HkmS_`0RHkl0 0`00HkmS_`15Hkl2001eHkl001]S_`03001S_f>o011S_`03001S_f>o00]S_`03001S_f>o025S_`03 001S_f>o04ES_`03001S_f>o07ES_`006f>o00<006>oHkl03f>o00<006>oHkl02f>o00<006>oHkl0 8F>o00<006>oHkl0A6>o0P00N6>o000KHkl00`00HkmS_`0?Hkl00`00HkmS_`0;Hkl00`00HkmS_`0Q Hkl00`00HkmS_`12Hkl2001jHkl001]S_`03001S_f>o00mS_`03001S_f>o00YS_`03001S_f>o025S _`03001S_f>o049S_`03001S_f>o07YS_`006f>o0P0046>o00<006>oHkl02V>o00<006>oHkl086>o 00<006>oHkl0@F>o0P00OF>o000KHkl00`00HkmS_`0>Hkl00`00HkmS_`0:Hkl00`00HkmS_`0PHkl0 0`00HkmS_`10Hkl2001oHkl001]S_`03001S_f>o00iS_`03001S_f>o00YS_`03001S_f>o01mS_`03 001S_f>o03mS_`80085S_`006f>o00<006>oHkl03V>o00<006>oHkl02F>o00<006>oHkl07f>o00<0 06>oHkl0?V>o0P00Pf>o000KHkl00`00HkmS_`0>Hkl00`00HkmS_`09Hkl00`00HkmS_`0NHkl00`00 HkmS_`0mHkl20025Hkl001]S_`03001S_f>o00eS_`03001S_f>o00US_`03001S_f>o01iS_`03001S _f>o03aS_`8008MS_`006f>o00<006>oHkl03F>o00<006>oHkl02F>o00<006>oHkl07F>o00<006>o Hkl0>f>o0P00RF>o000KHkl2000>Hkl00`00HkmS_`09Hkl00`00HkmS_`0LHkl00`00HkmS_`0kHkl0 0`00HkmS_`29Hkl001]S_`03001S_f>o00eS_`03001S_f>o00QS_`03001S_f>o01aS_`03001S_f>o 03YS_`8008aS_`006f>o00<006>oHkl036>o00<006>oHkl02F>o00<006>oHkl076>o00<006>oHkl0 >6>o0P00SV>o000KHkl00`00HkmS_`0Hkl000QS_`8000AS_`03001S_f>o00=S_`<000AS_`03001S_f>o00aS_`03001S _f>o00QS_`03001S_f>o01]S_`03001S_f>o03MS_`80095S_`001f>o00@006>oHkl000YS_`03001S _f>o00=S_`03001S_f>o00aS_`03001S_f>o00QS_`03001S_f>o01YS_`03001S_f>o03MS_`03001S _f>o095S_`001f>o00@006>oHkl000MS_`D000AS_`03001S_f>o00]S_`03001S_f>o00US_`03001S _f>o01US_`03001S_f>o03IS_`8009AS_`001f>o00@006>oHkl000MS_`04001S_f>o0005Hkl2000< Hkl00`00HkmS_`08Hkl00`00HkmS_`0IHkl00`00HkmS_`0fHkl00`00HkmS_`2DHkl000MS_`04001S _f>o0008Hkl00`00Hkl00005Hkl00`00HkmS_`0;Hkl00`00HkmS_`08Hkl00`00HkmS_`0IHkl00`00 HkmS_`0dHkl2002GHkl000MS_`04001S_f>o0009Hkl20005Hkl00`00HkmS_`0;Hkl00`00HkmS_`08 Hkl00`00HkmS_`0HHkl00`00HkmS_`0dHkl00`00HkmS_`2GHkl000QS_`8000]S_`03001S_f>o00=S _`03001S_f>o00]S_`03001S_f>o00MS_`03001S_f>o01QS_`03001S_f>o03=S_`8009YS_`006f>o 00<006>oHkl02V>o00<006>oHkl026>o00<006>oHkl05f>o00<006>oHkl0o00<006>oHkl0VV>o 000KHkl00`00HkmS_`0:Hkl00`00HkmS_`07Hkl00`00HkmS_`0GHkl00`00HkmS_`0cHkl00`00HkmS _`2KHkl001]S_`03001S_f>o00YS_`03001S_f>o00MS_`03001S_f>o01IS_`03001S_f>o039S_`80 09iS_`006f>o0P002f>o00<006>oHkl01f>o00<006>oHkl05V>o00<006>oHkl0o00<006>oHkl0 WV>o000KHkl00`00HkmS_`0:Hkl00`00HkmS_`06Hkl00`00HkmS_`0FHkl00`00HkmS_`0aHkl00`00 HkmS_`2OHkl001]S_`03001S_f>o00US_`03001S_f>o00MS_`03001S_f>o01ES_`03001S_f>o035S _`03001S_f>o0:1S_`006f>o00<006>oHkl02F>o00<006>oHkl01f>o00<006>oHkl05F>o00<006>o Hkl0;f>o0P00Xf>o000KHkl00`00HkmS_`09Hkl00`00HkmS_`06Hkl00`00HkmS_`0EHkl00`00HkmS _`0_Hkl00`00HkmS_`2SHkl001]S_`03001S_f>o00US_`03001S_f>o00IS_`03001S_f>o01AS_`03 001S_f>o02iS_`800:IS_`006f>o00<006>oHkl02F>o00<006>oHkl01V>o00<006>oHkl056>o00<0 06>oHkl0;F>o00<006>oHkl0YV>o000KHkl20009Hkl00`00HkmS_`06Hkl00`00HkmS_`0DHkl00`00 HkmS_`0]Hkl00`00HkmS_`2WHkl001]S_`03001S_f>o00QS_`03001S_f>o00IS_`03001S_f>o01=S _`03001S_f>o02aS_`800:YS_`006f>o00<006>oHkl026>o00<006>oHkl01V>o00<006>oHkl04f>o 00<006>oHkl0:f>o00<006>oHkl0ZV>o000KHkl00`00HkmS_`08Hkl00`00HkmS_`05Hkl00`00HkmS _`0CHkl00`00HkmS_`0ZHkl2002]Hkl001]S_`03001S_f>o00QS_`03001S_f>o00ES_`03001S_f>o 019S_`03001S_f>o02YS_`03001S_f>o0:eS_`006f>o00<006>oHkl026>o00<006>oHkl01F>o00<0 06>oHkl04V>o00<006>oHkl0:6>o0P00/6>o000KHkl00`00HkmS_`07Hkl00`00HkmS_`05Hkl00`00 HkmS_`0BHkl00`00HkmS_`0XHkl00`00HkmS_`2`Hkl001]S_`8000QS_`03001S_f>o00ES_`03001S _f>o015S_`03001S_f>o02QS_`03001S_f>o0;5S_`006f>o00<006>oHkl01f>o00<006>oHkl01F>o 00<006>oHkl04F>o00<006>oHkl09V>o0P00]6>o000KHkl00`00HkmS_`07Hkl00`00HkmS_`04Hkl0 0`00HkmS_`0AHkl00`00HkmS_`0VHkl00`00HkmS_`2dHkl001]S_`03001S_f>o00MS_`03001S_f>o 00AS_`03001S_f>o011S_`03001S_f>o02ES_`800;MS_`006f>o00<006>oHkl01V>o00<006>oHkl0 1F>o00<006>oHkl046>o00<006>oHkl096>o00<006>oHkl0]f>o0008Hkl20004Hkl00`00HkmS_`02 Hkl30005Hkl00`00HkmS_`06Hkl00`00HkmS_`05Hkl00`00HkmS_`0?Hkl00`00HkmS_`0THkl00`00 HkmS_`2hHkl000MS_`04001S_f>o0007Hkl01@00HkmS_f>o000016>o00<006>oHkl01V>o00<006>o Hkl016>o00<006>oHkl03f>o00<006>oHkl08f>o0P00^f>o0007Hkl01000HkmS_`001f>o00D006>o HkmS_`0000AS_`03001S_f>o00IS_`03001S_f>o00AS_`03001S_f>o00mS_`03001S_f>o029S_`03 001S_f>o0;]S_`001f>o00@006>oHkl000MS_`@000ES_`8000MS_`03001S_f>o00AS_`03001S_f>o 00iS_`03001S_f>o029S_`03001S_f>o0;aS_`001f>o00@006>oHkl000MS_`03001S_f>o00IS_`03 001S_f>o00IS_`03001S_f>o00=S_`03001S_f>o00mS_`03001S_f>o025S_`03001S_f>o0;eS_`00 1f>o00@006>oHkl000QS_`03001S_f>o00ES_`03001S_f>o00ES_`03001S_f>o00AS_`03001S_f>o 00iS_`03001S_f>o021S_`800<1S_`0026>o0P002V>o0P001F>o00<006>oHkl01F>o00<006>oHkl0 16>o00<006>oHkl03F>o00<006>oHkl086>o00<006>oHkl0`6>o000KHkl00`00HkmS_`05Hkl00`00 HkmS_`04Hkl00`00HkmS_`0=Hkl00`00HkmS_`0OHkl00`00HkmS_`31Hkl001]S_`03001S_f>o00ES _`03001S_f>o00=S_`03001S_f>o00eS_`03001S_f>o01mS_`03001S_f>o0<9S_`006f>o00<006>o Hkl01F>o00<006>oHkl00f>o00<006>oHkl036>o00<006>oHkl07V>o0P00aF>o000KHkl20005Hkl0 0`00HkmS_`04Hkl00`00HkmS_`0o 00AS_`03001S_f>o00=S_`03001S_f>o00aS_`03001S_f>o01eS_`03001S_f>o0o00<0 06>oHkl016>o00<006>oHkl00f>o00<006>oHkl036>o00<006>oHkl076>o00<006>oHkl0af>o000K Hkl00`00HkmS_`04Hkl00`00HkmS_`03Hkl00`00HkmS_`0;Hkl00`00HkmS_`0KHkl2003:Hkl001]S _`03001S_f>o00AS_`03001S_f>o00=S_`03001S_f>o00YS_`03001S_f>o01]S_`03001S_f>o0o00<006>oHkl016>o00<006>oHkl00V>o00<006>oHkl02f>o00<006>oHkl06F>o0P00cF>o 000KHkl00`00HkmS_`04Hkl00`00HkmS_`02Hkl00`00HkmS_`0:Hkl00`00HkmS_`0IHkl00`00HkmS _`3=Hkl001]S_`8000ES_`03001S_f>o009S_`03001S_f>o00YS_`03001S_f>o01QS_`03001S_f>o 0o00<006>oHkl00f>o00<006>oHkl00V>o00<006>oHkl02V>o00<006>oHkl05f>o0P00 dF>o000KHkl00`00HkmS_`03Hkl00`00HkmS_`02Hkl00`00HkmS_`09Hkl00`00HkmS_`0GHkl00`00 HkmS_`3AHkl001]S_`03001S_f>o00=S_`03001S_f>o009S_`03001S_f>o00US_`03001S_f>o01ES _`800=AS_`006f>o00<006>oHkl00f>o00<006>oHkl00V>o00<006>oHkl026>o00<006>oHkl05F>o 00<006>oHkl0e6>o000KHkl00`00HkmS_`03Hkl01@00HkmS_f>o00002f>o00<006>oHkl056>o00<0 06>oHkl0eF>o000KHkl00`00HkmS_`03Hkl01@00HkmS_f>o00002V>o00<006>oHkl056>o00<006>o Hkl0eV>o000KHkl20004Hkl01@00HkmS_f>o00002F>o00<006>oHkl04f>o0P00fF>o000KHkl00`00 HkmS_`03Hkl01000HkmS_`002V>o00<006>oHkl04V>o00<006>oHkl0fF>o000KHkl00`00HkmS_`03 Hkl01000HkmS_`002F>o00<006>oHkl04V>o00<006>oHkl0fV>o000KHkl00`00HkmS_`02Hkl01@00 HkmS_f>o00002F>o00<006>oHkl04F>o00<006>oHkl0ff>o000KHkl00`00HkmS_`02Hkl01@00HkmS _f>o000026>o00<006>oHkl046>o0P00gV>o0008Hkl20004Hkl00`00HkmS_`02Hkl30005Hkl00`00 HkmS_`02Hkl01000HkmS_`002F>o00<006>oHkl03f>o00<006>oHkl0gV>o0007Hkl01000HkmS_`00 1f>o00D006>oHkmS_`0000AS_`03001S_f>o009S_`04001S_f>o0008Hkl00`00HkmS_`0?Hkl00`00 HkmS_`3OHkl000MS_`04001S_f>o0007Hkl01@00HkmS_f>o000016>o0P000f>o00@006>oHkl000MS _`03001S_f>o00mS_`03001S_f>o0>1S_`001f>o00@006>oHkl000QS_`<000ES_`03001S_f>o009S _`04001S_f>o0007Hkl00`00HkmS_`0>Hkl00`00HkmS_`3QHkl000MS_`04001S_f>o0007Hkl01@00 HkmS_f>o000016>o00D006>oHkmS_`00009S_`03001S_f>o00ES_`03001S_f>o00iS_`03001S_f>o 0>9S_`001f>o00@006>oHkl000MS_`05001S_f>oHkl00004Hkl01@00HkmS_f>o00000V>o00<006>o Hkl01F>o00<006>oHkl03F>o00<006>oHkl0hf>o0008Hkl20009Hkl30005Hkl01@00HkmS_f>o0000 0V>o00<006>oHkl016>o00<006>oHkl03F>o00<006>oHkl0i6>o000KHkl01`00HkmS_f>o001S_`00 00IS_`03001S_f>o00aS_`800>MS_`006f>o00L006>oHkmS_`00Hkl00006Hkl00`00HkmS_`0;Hkl0 0`00HkmS_`3WHkl001]S_`80009S_`03001S_`0000ES_`03001S_f>o00]S_`03001S_f>o0>QS_`00 6f>o00@006>oHkl0009S_`03001S_f>o00=S_`03001S_f>o00YS_`03001S_f>o0>US_`006f>o00H0 06>oHkl006>o0005Hkl00`00HkmS_`09Hkl2003/Hkl001]S_`06001S_f>o001S_`0016>o00<006>o Hkl02F>o00<006>oHkl0k6>o000KHkl01P00HkmS_`00Hkl000AS_`03001S_f>o00QS_`03001S_f>o 0>eS_`006f>o00H006>oHkl006>o0003Hkl00`00HkmS_`08Hkl00`00HkmS_`3^Hkl001]S_`03001S _f>o008000AS_`03001S_f>o00IS_`800?5S_`006f>o00<006>oHkl00P000f>o00<006>oHkl01V>o 00<006>oHkl0lF>o000KHkl300000f>o001S_`02Hkl00`00HkmS_`05Hkl00`00HkmS_`3bHkl001]S _`05001S_`00Hkl00002Hkl00`00HkmS_`05Hkl00`00HkmS_`3cHkl001]S_`04001S_`000003Hkl0 0`00HkmS_`03Hkl2003fHkl001]S_`04001S_`000002Hkl00`00HkmS_`03Hkl00`00HkmS_`3fHkl0 01]S_`04001S_`000002Hkl00`00HkmS_`02Hkl00`00HkmS_`3gHkl001]S_`06001S_`00001S_`00 16>o00<006>oHkl0n6>o000KHkl300000f>o001S_`02Hkl2003kHkl001]S_`<00005Hkl006>oHkl0 003mHkl001]S_`@0009S_`03001S_f>o0?aS_`006f>o0P0000AS_`00Hkl00?mS_`006f>o0`0000=S _`00Hkl0of>o000CHkl30005Hkl4003oHkl2Hkl001AS_`03001S_f>o00AS_`<00?mS_`=S_`0056>o 00<006>oHkl016>o0P00of>o16>o000DHkl00`00HkmS_`04Hkl2003oHkl4Hkl001AS_`03001S_f>o 00AS_`03001S_f>o0?mS_`=S_`004f>o0P001V>o00<006>oHkl0of>o0f>o000DHkl00`00HkmS_`04 Hkl00`00HkmS_`3oHkl3Hkl00?mS_b5S_`00of>o8F>o003oHklQHkl00?mS_b5S_`00of>o8F>o003o HklQHkl00?mS_b5S_`00of>o8F>o003oHklQHkl001IS_`<00?mS_`QS_`005f>o00<006>oHkl0of>o 1f>o000GHkl00`00HkmS_`3oHkl7Hkl001MS_`03001S_f>o0?mS_`MS_`005F>o00D006>o001S_`00 009S_`03001S_f>o0?mS_`9S_`005F>o00D006>o001S_`0000=S_`03001S_f>o0?mS_`5S_`005F>o 1@00of>o1f>o003oHklQHkl00?mS_b5S_`00\ \>"], ImageRangeCache->{{{0, 287}, {176.938, 0}} -> {-0.586937, -0.0840746, \ 0.0215746, 0.00698171}}] }, Open ]] }, Open ]] }, Closed]] }, Closed]] }, Closed]], Cell[CellGroupData[{ Cell["Conclusions", "Subtitle"], Cell[TextData[{ "1. For natural convection flows where the ", ButtonBox["Grashof", ButtonData:>"Grashof", ButtonStyle->"Hyperlink"], " number is larger, a ", ButtonBox["boundary layer", ButtonData:>"boundary_layer", ButtonStyle->"Hyperlink"], " can be expected close to a solid surface. \n\n2. For situations where a \ natural convection boundary-layer is occurring, heat transfer will be \ governed by ", ButtonBox["coupled energy and momentum equations", ButtonData:>"govern_eqs", ButtonStyle->"Hyperlink"], ". \n\n3. For transport processes that occur on a semi-infinite domain, \ where there is ", StyleBox["no geometric length scale", FontVariations->{"Underline"->True}], ", it is often possible to define a (dimensionless) ", ButtonBox["similarity", ButtonData:>"similarity_variable", ButtonStyle->"Hyperlink"], " variable that contains the natural length scale. \n\n4. It is often \ possible to", ButtonBox[" reduce PDE's to ODE's", ButtonData:>"pde_to_ode", ButtonStyle->"Hyperlink"], " through simplifications made possible by use of the similarity variable. \ \n\n5. From the scaling identified by the similarity variable, it is often \ possible to ", ButtonBox["predict the macroscopic behavior", ButtonData:>"gr_1/4", ButtonStyle->"Hyperlink"], " without solving the differential equations. In this case Nu ~ ", Cell[BoxData[ \(TraditionalForm\`Gr\^\(1/4\)\)]], ". \n\n6. This prediction agrees with the ", ButtonBox["recommended correlation", ButtonData:>"correlation", ButtonStyle->"Hyperlink"], " for high Gr heat transfer.\n\n7. The coupled, nonlinear ODE's can be \ readily solved with a ", ButtonBox["shooting method", ButtonData:>"shooting", ButtonStyle->"Hyperlink"], ".\n\n8. Note the shape of the ", ButtonBox["temperature profile and the location of the maximum velocity", ButtonData:>"results", ButtonStyle->"Hyperlink"], ". " }], "Text", CellTags->"conclusions"] }, Closed]] }, FrontEndVersion->"4.0 for Macintosh", ScreenRectangle->{{0, 1024}, {0, 748}}, WindowToolbars->{}, CellGrouping->Manual, WindowSize->{582, 692}, WindowMargins->{{Automatic, 63}, {12, Automatic}}, PrintingCopies->1, PrintingPageRange->{1, Automatic}, PrintingOptions->{"PrintingMargins"->{{54, 54}, {72, 72}}, "PrintCellBrackets"->False, "PrintRegistrationMarks"->True, "PrintMultipleHorizontalPages"->False}, PrivateNotebookOptions->{"ColorPalette"->{RGBColor, 128}}, ShowCellLabel->True, CellLabelAutoDelete->True, ShowCellTags->False, SpellingDictionaries->{"CorrectWords"->{ "nondimensionalized.", "Couette", "Tstar", "Ts", "Tfinal", "ufinal", "finit", "deltv", "youter", "alf", "deltvnew", "finitnew", "finitold", "Uerror", "Tprime"}}, RenderingOptions->{"ObjectDithering"->True, "RasterDithering"->False}, CharacterEncoding->"MacintoshAutomaticEncoding", StyleDefinitions -> Notebook[{ Cell[CellGroupData[{ Cell["Style Definitions", "Subtitle"], Cell["\<\ Modify the definitions below to change the default appearance of \ all cells in a given style. Make modifications to any definition using commands in the Format menu.\ \>", "Text"], Cell[CellGroupData[{ Cell["Style Environment Names", "Section"], Cell[StyleData[All, "Working"], ScriptMinSize->9], Cell[StyleData[All, "Printout"], PageWidth->PaperWidth, ShowCellLabel->False, ImageSize->{200, 200}, PrivateFontOptions->{"FontType"->"Outline"}] }, Closed]], Cell[CellGroupData[{ Cell["Notebook Options", "Section"], Cell["\<\ The options defined for the style below will be used at the \ Notebook level.\ \>", "Text"], Cell[StyleData["Notebook"], PageHeaders->{{Cell[ TextData[ { CounterBox[ "Page"]}], "PageNumber"], None, Cell[ TextData[ { ValueBox[ "FileName"]}], "Header"]}, {Cell[ TextData[ { ValueBox[ "FileName"]}], "Header"], None, Cell[ TextData[ { CounterBox[ "Page"]}], "PageNumber"]}}, PageHeaderLines->{True, True}, PrintingOptions->{"FirstPageHeader"->False, "FacingPages"->True}, CellLabelAutoDelete->False, CellFrameLabelMargins->6, StyleMenuListing->None] }, Closed]], Cell[CellGroupData[{ Cell["Styles for Headings", "Section"], Cell[CellGroupData[{ Cell[StyleData["Title"], CellFrame->{{0, 0}, {0, 0.25}}, CellMargins->{{18, 30}, {4, 20}}, CellGroupingRules->{"TitleGrouping", 0}, PageBreakBelow->False, CellFrameMargins->9, LineSpacing->{0.95, 0}, CounterIncrements->"Title", CounterAssignments->{{"Section", 0}, {"Equation", 0}, {"Figure", 0}}, FontSize->36, Background->RGBColor[0.750011, 1, 0.750011]], Cell[StyleData["Title", "Printout"], CellMargins->{{18, 30}, {4, 0}}, CellFrameMargins->4, FontSize->30] }, Open ]], Cell[CellGroupData[{ Cell[StyleData["Subtitle"], CellMargins->{{18, 30}, {0, 10}}, CellGroupingRules->{"TitleGrouping", 10}, PageBreakBelow->False, LineSpacing->{1, 0}, CounterIncrements->"Subtitle", CounterAssignments->{{"Section", 0}, {"Equation", 0}, {"Figure", 0}}, FontSize->24, FontSlant->"Italic", Background->RGBColor[0.8, 0.920012, 0.920012]], Cell[StyleData["Subtitle", "Printout"], CellMargins->{{18, 30}, {0, 10}}, FontSize->18] }, Closed]], Cell[CellGroupData[{ Cell[StyleData["SectionFirst"], CellFrame->{{0, 0}, {0, 3}}, CellMargins->{{18, 30}, {4, 30}}, CellGroupingRules->{"SectionGrouping", 40}, PageBreakBelow->False, CellFrameMargins->3, CounterIncrements->"Section", CounterAssignments->{{"Subsection", 0}, {"Subsubsection", 0}}, FontSize->18, FontWeight->"Bold"], Cell[StyleData["SectionFirst", "Printout"], FontSize->14] }, Closed]], Cell[CellGroupData[{ Cell[StyleData["Section"], CellMargins->{{18, 30}, {4, 30}}, CellGroupingRules->{"SectionGrouping", 40}, PageBreakBelow->False, CounterIncrements->"Section", CounterAssignments->{{"Subsection", 0}, {"Subsubsection", 0}}, FontSize->18, FontWeight->"Bold", Background->RGBColor[0.920012, 0.870024, 0.770016]], Cell[StyleData["Section", "Printout"], FontSize->14] }, Open ]], Cell[CellGroupData[{ Cell[StyleData["Subsection"], CellDingbat->"\[FilledSquare]", CellMargins->{{18, 30}, {4, 20}}, CellGroupingRules->{"SectionGrouping", 50}, PageBreakBelow->False, CounterIncrements->"Subsection", CounterAssignments->{{"Subsubsection", 0}}, FontSize->14, FontWeight->"Bold", Background->RGBColor[0.970001, 0.890013, 0.850004]], Cell[StyleData["Subsection", "Printout"], FontSize->12] }, Closed]], Cell[CellGroupData[{ Cell[StyleData["Subsubsection"], CellDingbat->"\[FilledSmallSquare]", CellMargins->{{18, 30}, {4, 12}}, CellGroupingRules->{"SectionGrouping", 60}, PageBreakBelow->False, CounterIncrements->"Subsubsection", FontSize->12, FontWeight->"Bold", Background->RGBColor[0.870008, 1, 0.960006]], Cell[StyleData["Subsubsection", "Printout"], FontSize->10] }, Closed]] }, Open ]], Cell[CellGroupData[{ Cell["Styles for Body Text", "Section"], Cell[CellGroupData[{ Cell[StyleData["Text"], CellMargins->{{18, 10}, {Inherited, 6}}, TextJustification->1, LineSpacing->{1, 2}, CounterIncrements->"Text", Background->RGBColor[0.900008, 1, 0.940002]], Cell[StyleData["Text", "Printout"], CellMargins->{{18, 30}, {Inherited, 4}}, LineSpacing->{1, 3}, FontSize->10] }, Closed]], Cell[CellGroupData[{ Cell[StyleData["Caption"], CellMargins->{{55, 50}, {5, 5}}, PageBreakAbove->False, FontSize->10], Cell[StyleData["Caption", "Printout"], CellMargins->{{55, 55}, {5, 2}}, FontSize->8] }, Closed]] }, Closed]], Cell[CellGroupData[{ Cell["Styles for Input/Output", "Section"], Cell["\<\ The cells in this section define styles used for input and output \ to the kernel. Be careful when modifying, renaming, or removing these \ styles, because the front end associates special meanings with these style \ names.\ \>", "Text"], Cell[CellGroupData[{ Cell[StyleData["Input"], CellMargins->{{55, 10}, {5, 8}}, Evaluatable->True, CellGroupingRules->"InputGrouping", CellHorizontalScrolling->True, PageBreakWithin->False, GroupPageBreakWithin->False, CellLabelMargins->{{26, Inherited}, {Inherited, Inherited}}, DefaultFormatType->DefaultInputFormatType, AutoItalicWords->{}, FormatType->InputForm, ShowStringCharacters->True, NumberMarks->True, CounterIncrements->"Input", FontSize->12, FontWeight->"Bold"], Cell[StyleData["Input", "Printout"], CellMargins->{{55, 55}, {0, 10}}, ShowCellLabel->False, FontSize->9.5] }, Closed]], Cell[CellGroupData[{ Cell[StyleData["Output"], CellMargins->{{55, 10}, {8, 5}}, CellEditDuplicate->True, CellGroupingRules->"OutputGrouping", CellHorizontalScrolling->True, PageBreakWithin->False, GroupPageBreakWithin->False, GeneratedCell->True, CellAutoOverwrite->True, CellLabelPositioning->Left, CellLabelMargins->{{26, Inherited}, {Inherited, Inherited}}, DefaultFormatType->DefaultOutputFormatType, AutoItalicWords->{}, FormatType->InputForm, CounterIncrements->"Output", Background->RGBColor[0.940017, 0.890013, 0.990005]], Cell[StyleData["Output", "Printout"], CellMargins->{{55, 55}, {10, 10}}, ShowCellLabel->False, FontSize->9.5] }, Open ]], Cell[CellGroupData[{ Cell[StyleData["Message"], CellDingbat->"\[LongDash]", CellMargins->{{55, Inherited}, {Inherited, Inherited}}, CellGroupingRules->"OutputGrouping", PageBreakWithin->False, GroupPageBreakWithin->False, GeneratedCell->True, CellAutoOverwrite->True, ShowCellLabel->False, CellLabelMargins->{{26, Inherited}, {Inherited, Inherited}}, DefaultFormatType->DefaultOutputFormatType, AutoItalicWords->{}, FormatType->InputForm, CounterIncrements->"Message", StyleMenuListing->None, FontSize->10, FontSlant->"Italic"], Cell[StyleData["Message", "Printout"], CellMargins->{{55, 55}, {0, 3}}, FontSize->8] }, Closed]], Cell[CellGroupData[{ Cell[StyleData["Print"], CellMargins->{{55, Inherited}, {Inherited, Inherited}}, CellGroupingRules->"OutputGrouping", CellHorizontalScrolling->True, PageBreakWithin->False, GroupPageBreakWithin->False, GeneratedCell->True, CellAutoOverwrite->True, ShowCellLabel->False, CellLabelMargins->{{26, Inherited}, {Inherited, Inherited}}, DefaultFormatType->DefaultOutputFormatType, TextAlignment->Left, AutoItalicWords->{}, FormatType->InputForm, CounterIncrements->"Print", StyleMenuListing->None], Cell[StyleData["Print", "Printout"], CellMargins->{{54, 72}, {2, 10}}, FontSize->8] }, Closed]], Cell[CellGroupData[{ Cell[StyleData["Graphics"], CellMargins->{{55, Inherited}, {Inherited, Inherited}}, CellGroupingRules->"GraphicsGrouping", CellHorizontalScrolling->True, PageBreakWithin->False, GeneratedCell->True, CellAutoOverwrite->True, ShowCellLabel->False, DefaultFormatType->DefaultOutputFormatType, FormatType->InputForm, CounterIncrements->"Graphics", StyleMenuListing->None, Background->RGBColor[0.780011, 0.920012, 1], ButtonBoxOptions->{Background->RGBColor[0.750011, 0.630014, 0.760021]}], Cell[StyleData["Graphics", "Printout"], CellMargins->{{55, 55}, {0, 15}}, ImageSize->{0.0625, 0.0625}, ImageMargins->{{35, Inherited}, {Inherited, 0}}, FontSize->8] }, Open ]], Cell[CellGroupData[{ Cell[StyleData["CellLabel"], CellMargins->{{9, Inherited}, {Inherited, Inherited}}, StyleMenuListing->None, FontFamily->"Helvetica", FontSize->9, FontSlant->"Oblique"], Cell[StyleData["CellLabel", "Printout"], CellMargins->{{0, Inherited}, {Inherited, Inherited}}, FontSize->8] }, Closed]] }, Open ]], Cell[CellGroupData[{ Cell["Unique Styles", "Section"], Cell[CellGroupData[{ Cell[StyleData["Author"], CellMargins->{{45, Inherited}, {2, 20}}, CellGroupingRules->{"TitleGrouping", 20}, PageBreakBelow->False, CounterAssignments->{{"Section", 0}, {"Equation", 0}, {"Figure", 0}}, FontSize->14, FontWeight->"Bold"], Cell[StyleData["Author", "Printout"], CellMargins->{{36, Inherited}, {2, 30}}, FontSize->12] }, Closed]], Cell[CellGroupData[{ Cell[StyleData["Address"], CellMargins->{{45, Inherited}, {2, 2}}, CellGroupingRules->{"TitleGrouping", 30}, PageBreakBelow->False, LineSpacing->{1, 1}, CounterAssignments->{{"Section", 0}, {"Equation", 0}, {"Figure", 0}}, FontSize->12, FontSlant->"Italic"], Cell[StyleData["Address", "Printout"], CellMargins->{{36, Inherited}, {2, 2}}, FontSize->10] }, Closed]], Cell[CellGroupData[{ Cell[StyleData["Abstract"], CellMargins->{{45, 75}, {Inherited, 30}}, LineSpacing->{1, 0}], Cell[StyleData["Abstract", "Printout"], CellMargins->{{36, 67}, {Inherited, 50}}, FontSize->10] }, Closed]], Cell[CellGroupData[{ Cell[StyleData["Reference"], CellMargins->{{18, 40}, {2, 2}}, TextJustification->1, LineSpacing->{1, 0}], Cell[StyleData["Reference", "Printout"], CellMargins->{{18, 40}, {Inherited, 0}}, FontSize->8] }, Closed]] }, Closed]], Cell[CellGroupData[{ Cell["Styles for Automatic Numbering", "Section"], Cell["\<\ The following styles are useful for numbered equations, figures, \ etc. They automatically give the cell a FrameLabel containing a reference to \ a particular counter, and also increment that counter.\ \>", "Text"], Cell[CellGroupData[{ Cell[StyleData["NumberedEquation"], CellMargins->{{55, 10}, {0, 10}}, CellFrameLabels->{{None, Cell[ TextData[ {"(", CounterBox[ "NumberedEquation"], ")"}]]}, {None, None}}, DefaultFormatType->DefaultInputFormatType, CounterIncrements->"NumberedEquation", FormatTypeAutoConvert->False], Cell[StyleData["NumberedEquation", "Printout"], CellMargins->{{55, 55}, {0, 10}}, FontSize->10] }, Closed]], Cell[CellGroupData[{ Cell[StyleData["NumberedFigure"], CellMargins->{{55, 145}, {2, 10}}, CellHorizontalScrolling->True, CellFrameLabels->{{None, None}, {Cell[ TextData[ {"Figure ", CounterBox[ "NumberedFigure"]}], FontWeight -> "Bold"], None}}, CounterIncrements->"NumberedFigure", FormatTypeAutoConvert->False], Cell[StyleData["NumberedFigure", "Printout"], FontSize->10] }, Closed]], Cell[CellGroupData[{ Cell[StyleData["NumberedTable"], CellMargins->{{55, 145}, {2, 10}}, CellFrameLabels->{{None, None}, {Cell[ TextData[ {"Table ", CounterBox[ "NumberedTable"]}], FontWeight -> "Bold"], None}}, TextAlignment->Center, CounterIncrements->"NumberedTable", FormatTypeAutoConvert->False], Cell[StyleData["NumberedTable", "Printout"], CellMargins->{{18, Inherited}, {Inherited, Inherited}}, FontSize->10] }, Closed]] }, Closed]], Cell[CellGroupData[{ Cell["Formulas and Programming", "Section"], Cell[CellGroupData[{ Cell[StyleData["DisplayFormula"], CellMargins->{{55, 10}, {2, 10}}, CellHorizontalScrolling->True, DefaultFormatType->DefaultInputFormatType, ScriptLevel->0, SingleLetterItalics->True, UnderoverscriptBoxOptions->{LimitsPositioning->True}], Cell[StyleData["DisplayFormula", "Printout"], FontSize->10] }, Closed]], Cell[CellGroupData[{ Cell[StyleData["ChemicalFormula"], CellMargins->{{55, 10}, {2, 10}}, DefaultFormatType->DefaultInputFormatType, AutoSpacing->False, ScriptLevel->1, ScriptBaselineShifts->{0.6, Automatic}, SingleLetterItalics->False, ZeroWidthTimes->True], Cell[StyleData["ChemicalFormula", "Printout"], FontSize->10] }, Closed]], Cell[CellGroupData[{ Cell[StyleData["Program"], CellMargins->{{18, 10}, {Inherited, 6}}, FontFamily->"Courier"], Cell[StyleData["Program", "Printout"], CellMargins->{{18, 30}, {Inherited, 4}}, FontSize->9.5] }, Closed]] }, Closed]] }, Open ]] }], MacintoshSystemPageSetup->"\<\ 01@0004/0B`0000034H9H?o/okD" ] (*********************************************************************** 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->{ "standard"->{ Cell[6987, 194, 456, 8, 61, "Input", CellTags->"standard"]}, "traditional"->{ Cell[9132, 249, 504, 12, 84, "Text", CellTags->"traditional"]}, "Grashof"->{ Cell[13634, 385, 1883, 32, 285, "Text", CellTags->"Grashof"]}, "boundary_layer"->{ Cell[15759, 430, 623, 10, 110, "Text", CellTags->"boundary_layer"]}, "govern_eqs"->{ Cell[19653, 508, 142, 4, 76, "Subsection", CellTags->"govern_eqs"]}, "similarity_variable"->{ Cell[29636, 804, 119, 2, 71, "Output", CellTags->"similarity_variable"]}, "pde_to_ode"->{ Cell[34846, 981, 133, 4, 42, "Subsection", CellTags->"pde_to_ode"]}, "gr_1/4"->{ Cell[71697, 1948, 642, 11, 125, "Text", CellTags->"gr_1/4"]}, "correlation"->{ Cell[72565, 1971, 656, 16, 151, "Text", CellTags->"correlation"]}, "shooting"->{ Cell[81108, 2216, 1369, 43, 262, "Text", CellTags->"shooting"]}, "results"->{ Cell[131571, 3715, 353, 8, 70, "Text", CellTags->"results"]}, "conclusions"->{ Cell[176620, 5176, 2036, 52, 307, "Text", CellTags->"conclusions"]} } *) (*CellTagsIndex CellTagsIndex->{ {"standard", 193956, 5776}, {"traditional", 194046, 5779}, {"Grashof", 194135, 5782}, {"boundary_layer", 194230, 5785}, {"govern_eqs", 194327, 5788}, {"similarity_variable", 194433, 5791}, {"pde_to_ode", 194535, 5794}, {"gr_1/4", 194628, 5797}, {"correlation", 194719, 5800}, {"shooting", 194812, 5803}, {"results", 194902, 5806}, {"conclusions", 194993, 5809} } *) (*NotebookFileOutline Notebook[{ Cell[1717, 49, 106, 3, 150, "Title"], Cell[1826, 54, 895, 18, 323, "Text"], Cell[2724, 74, 142, 3, 64, "Subsubsection"], Cell[CellGroupData[{ Cell[2891, 81, 27, 0, 53, "Subtitle"], Cell[2921, 83, 1121, 19, 166, "Text", Evaluatable->False], Cell[CellGroupData[{ Cell[4067, 106, 96, 4, 58, "Subsection"], Cell[4166, 112, 457, 9, 84, "Text"], Cell[CellGroupData[{ Cell[4648, 125, 39, 0, 48, "Subsubsection"], Cell[4690, 127, 88, 3, 40, "Text"], Cell[CellGroupData[{ Cell[4803, 134, 406, 7, 92, "Input"], Cell[5212, 143, 1550, 38, 69, "Output"] }, Open ]], Cell[6777, 184, 185, 6, 42, "Text"], Cell[CellGroupData[{ Cell[6987, 194, 456, 8, 61, "Input", CellTags->"standard"], Cell[7446, 204, 1550, 38, 69, "Output"] }, Open ]], Cell[9011, 245, 118, 2, 40, "Text"], Cell[9132, 249, 504, 12, 84, "Text", CellTags->"traditional"], Cell[CellGroupData[{ Cell[9661, 265, 1578, 41, 84, "Input", Evaluatable->False], Cell[11242, 308, 1601, 40, 69, "Output", Evaluatable->False] }, Open ]], Cell[12858, 351, 118, 2, 40, "Text"], Cell[12979, 355, 226, 4, 54, "Text"] }, Open ]] }, Closed]] }, Closed]], Cell[CellGroupData[{ Cell[13266, 366, 28, 0, 47, "Subtitle"], Cell[13297, 368, 162, 4, 54, "Text"] }, Closed]], Cell[CellGroupData[{ Cell[13496, 377, 36, 0, 47, "Subtitle"], Cell[CellGroupData[{ Cell[13557, 381, 52, 0, 58, "Subsection"], Cell[CellGroupData[{ Cell[13634, 385, 1883, 32, 285, "Text", CellTags->"Grashof"], Cell[15520, 419, 132, 3, 32, "Input"] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell[15701, 428, 55, 0, 42, "Subsection"], Cell[15759, 430, 623, 10, 110, "Text", CellTags->"boundary_layer"], Cell[16385, 442, 132, 3, 32, "Input"] }, Closed]], Cell[CellGroupData[{ Cell[16554, 450, 44, 0, 42, "Subsection"], Cell[CellGroupData[{ Cell[16623, 454, 2846, 43, 452, "Text"], Cell[19472, 499, 132, 3, 32, "Input"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[19653, 508, 142, 4, 76, "Subsection", CellTags->"govern_eqs"], Cell[19798, 514, 26, 0, 40, "Text"], Cell[19827, 516, 357, 9, 47, "Input", Evaluatable->False], Cell[20187, 527, 24, 0, 40, "Text"], Cell[20214, 529, 655, 14, 50, "Input", Evaluatable->False], Cell[20872, 545, 22, 0, 40, "Text"], Cell[20897, 547, 589, 13, 50, "Input", Evaluatable->False], Cell[21489, 562, 132, 3, 32, "Input"], Cell[21624, 567, 187, 4, 54, "Text"] }, Closed]] }, Closed]], Cell[CellGroupData[{ Cell[21860, 577, 116, 2, 47, "Subtitle", Evaluatable->False], Cell[CellGroupData[{ Cell[22001, 583, 62, 0, 58, "Subsection"], Cell[CellGroupData[{ Cell[22088, 587, 62, 0, 48, "Subsubsection"], Cell[22153, 589, 2495, 46, 418, "Text"] }, Open ]], Cell[CellGroupData[{ Cell[24685, 640, 119, 3, 48, "Subsubsection"], Cell[CellGroupData[{ Cell[24829, 647, 112, 2, 28, "Input"], Cell[24944, 651, 417, 13, 46, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[25398, 669, 201, 4, 44, "Input"], Cell[25602, 675, 772, 21, 46, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[26411, 701, 145, 3, 44, "Input"], Cell[26559, 706, 729, 21, 46, "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[27337, 733, 74, 0, 48, "Subsubsection"], Cell[27414, 735, 266, 5, 68, "Text"], Cell[27683, 742, 324, 9, 62, "Text"], Cell[CellGroupData[{ Cell[28032, 755, 904, 22, 69, "Input"], Cell[28939, 779, 52, 1, 43, "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[29040, 786, 56, 0, 48, "Subsubsection"], Cell[29099, 788, 271, 5, 68, "Text"], Cell[29373, 795, 191, 4, 54, "Text"], Cell[29567, 801, 66, 1, 44, "Input"], Cell[29636, 804, 119, 2, 71, "Output", CellTags->"similarity_variable"], Cell[29758, 808, 115, 2, 40, "Text"], Cell[29876, 812, 49, 0, 40, "Text"], Cell[CellGroupData[{ Cell[29950, 816, 89, 1, 44, "Input"], Cell[30042, 819, 112, 2, 49, "Output"] }, Open ]], Cell[30169, 824, 148, 3, 54, "Text"] }, Open ]], Cell[CellGroupData[{ Cell[30354, 832, 112, 3, 48, "Subsubsection"], Cell[30469, 837, 414, 10, 76, "Text"], Cell[CellGroupData[{ Cell[30908, 851, 68, 1, 28, "Input"], Cell[30979, 854, 446, 11, 49, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[31462, 870, 78, 1, 28, "Input"], Cell[31543, 873, 446, 11, 49, "Output"] }, Open ]], Cell[32004, 887, 75, 2, 40, "Text"], Cell[CellGroupData[{ Cell[32104, 893, 86, 1, 28, "Input"], Cell[32193, 896, 285, 7, 65, "Output"] }, Open ]], Cell[32493, 906, 313, 9, 62, "Text"], Cell[CellGroupData[{ Cell[32831, 919, 54, 1, 28, "Input"], Cell[32888, 922, 786, 20, 66, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[33711, 947, 84, 1, 28, "Input"], Cell[33798, 950, 987, 24, 72, "Output"] }, Open ]] }, Closed]] }, Closed]], Cell[CellGroupData[{ Cell[34846, 981, 133, 4, 42, "Subsection", CellTags->"pde_to_ode"], Cell[34982, 987, 128, 3, 40, "Text"], Cell[35113, 992, 115, 2, 40, "Text"], Cell[CellGroupData[{ Cell[35253, 998, 54, 0, 48, "Subsubsection"], Cell[CellGroupData[{ Cell[35332, 1002, 911, 24, 69, "Input"], Cell[36246, 1028, 4893, 108, 260, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[41176, 1141, 296, 6, 76, "Input"], Cell[41475, 1149, 4003, 85, 341, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[45515, 1239, 75, 1, 28, "Input"], Cell[45593, 1242, 3808, 82, 306, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[49438, 1329, 71, 1, 28, "Input"], Cell[49512, 1332, 1273, 31, 82, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[50822, 1368, 116, 2, 28, "Input"], Cell[50941, 1372, 818, 19, 62, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[51796, 1396, 126, 3, 28, "Input"], Cell[51925, 1401, 877, 20, 82, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[52839, 1426, 162, 4, 44, "Input"], Cell[53004, 1432, 661, 18, 46, "Output"] }, Open ]] }, Closed]], Cell[53692, 1454, 115, 2, 38, "Text"], Cell[CellGroupData[{ Cell[53832, 1460, 52, 0, 48, "Subsubsection"], Cell[CellGroupData[{ Cell[53909, 1464, 1312, 33, 88, "Input"], Cell[55224, 1499, 3222, 77, 170, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[58483, 1581, 296, 6, 76, "Input"], Cell[58782, 1589, 2747, 60, 233, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[61566, 1654, 75, 1, 28, "Input"], Cell[61644, 1657, 2663, 58, 233, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[64344, 1720, 65, 1, 28, "Input"], Cell[64412, 1723, 1406, 32, 106, "Output"] }, Open ]], Cell[65833, 1758, 238, 6, 55, "Text"], Cell[CellGroupData[{ Cell[66096, 1768, 116, 2, 28, "Input"], Cell[66215, 1772, 1232, 30, 106, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[67484, 1807, 119, 3, 28, "Input"], Cell[67606, 1812, 609, 14, 60, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[68252, 1831, 104, 2, 28, "Input"], Cell[68359, 1835, 534, 13, 69, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[68930, 1853, 67, 1, 28, "Input"], Cell[69000, 1856, 405, 10, 44, "Output"] }, Open ]] }, Open ]], Cell[69432, 1870, 115, 2, 40, "Text"] }, Open ]], Cell[CellGroupData[{ Cell[69584, 1877, 56, 0, 58, "Subsection"], Cell[CellGroupData[{ Cell[69665, 1881, 39, 1, 28, "Input"], Cell[69707, 1884, 661, 18, 46, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[70405, 1907, 39, 1, 28, "Input"], Cell[70447, 1910, 405, 10, 44, "Output"] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell[70901, 1926, 41, 0, 42, "Subsection"], Cell[70945, 1928, 533, 9, 173, "Text"] }, Closed]], Cell[CellGroupData[{ Cell[71515, 1942, 61, 0, 42, "Subsection"], Cell[71579, 1944, 115, 2, 40, "Text"], Cell[71697, 1948, 642, 11, 125, "Text", CellTags->"gr_1/4"] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell[72388, 1965, 56, 0, 47, "Subtitle"], Cell[72447, 1967, 115, 2, 40, "Text"], Cell[72565, 1971, 656, 16, 151, "Text", CellTags->"correlation"] }, Closed]], Cell[CellGroupData[{ Cell[73258, 1992, 206, 4, 128, "Subtitle"], Cell[CellGroupData[{ Cell[73489, 2000, 76, 0, 48, "Subsubsection"], Cell[73568, 2002, 6925, 192, 359, "Text"] }, Closed]], Cell[CellGroupData[{ Cell[80530, 2199, 457, 11, 72, "Subsubsection"], Cell[80990, 2212, 115, 2, 40, "Text"], Cell[81108, 2216, 1369, 43, 262, "Text", CellTags->"shooting"] }, Closed]], Cell[CellGroupData[{ Cell[82514, 2264, 52, 0, 42, "Subsection"], Cell[CellGroupData[{ Cell[82591, 2268, 51, 0, 48, "Subsubsection"], Cell[82645, 2270, 979, 23, 108, "Input"], Cell[CellGroupData[{ Cell[83649, 2297, 97, 3, 40, "Text"], Cell[CellGroupData[{ Cell[83771, 2304, 38, 1, 28, "Input"], Cell[83812, 2307, 228, 7, 58, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[84077, 2319, 50, 1, 28, "Input"], Cell[84130, 2322, 55, 1, 43, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[84222, 2328, 50, 1, 28, "Input"], Cell[84275, 2331, 61, 1, 43, "Output"] }, Open ]] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell[84397, 2339, 165, 5, 40, "Subsubsection", Evaluatable->False], Cell[CellGroupData[{ Cell[84587, 2348, 2008, 39, 220, "Input"], Cell[86598, 2389, 523, 9, 22, "Print"] }, Open ]], Cell[87136, 2401, 225, 4, 54, "Text"], Cell[CellGroupData[{ Cell[87386, 2409, 37, 1, 28, "Input"], Cell[87426, 2412, 70, 1, 43, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[87533, 2418, 37, 1, 28, "Input"], Cell[87573, 2421, 75, 1, 43, "Output"] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell[87697, 2428, 57, 0, 40, "Subsubsection"], Cell[CellGroupData[{ Cell[87779, 2432, 8220, 156, 828, "Input"], Cell[96002, 2590, 582, 11, 22, "Print"], Cell[96587, 2603, 592, 11, 22, "Print"], Cell[97182, 2616, 595, 11, 22, "Print"], Cell[97780, 2629, 593, 11, 43, "Print"] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell[98422, 2646, 59, 0, 40, "Subsubsection"], Cell[CellGroupData[{ Cell[98506, 2650, 312, 6, 64, "Input"], Cell[98821, 2658, 18455, 587, 198, 8094, 453, "GraphicsData", "PostScript", \ "Graphics", Evaluatable->False] }, Open ]], Cell[117291, 3248, 54, 0, 40, "Text"], Cell[CellGroupData[{ Cell[117370, 3252, 169, 4, 54, "Text"], Cell[CellGroupData[{ Cell[117564, 3260, 266, 6, 44, "Input"], Cell[117833, 3268, 13689, 441, 202, 5676, 338, "GraphicsData", "PostScript", \ "Graphics"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[131571, 3715, 353, 8, 70, "Text", CellTags->"results"], Cell[CellGroupData[{ Cell[131949, 3727, 72, 1, 28, "Input"], Cell[132024, 3730, 22863, 758, 202, 10367, 599, "GraphicsData", "PostScript", \ "Graphics"] }, Open ]], Cell[154902, 4491, 132, 3, 32, "Input"], Cell[CellGroupData[{ Cell[155059, 4498, 70, 1, 28, "Input"], Cell[155132, 4501, 21369, 664, 202, 8128, 496, "GraphicsData", "PostScript", \ "Graphics"] }, Open ]] }, Open ]] }, Closed]] }, Closed]] }, Closed]], Cell[CellGroupData[{ Cell[176586, 5174, 31, 0, 47, "Subtitle"], Cell[176620, 5176, 2036, 52, 307, "Text", CellTags->"conclusions"] }, Closed]] } ] *) (*********************************************************************** End of Mathematica Notebook file. ***********************************************************************)