(* Content-type: application/mathematica *) (*** Wolfram Notebook File ***) (* http://www.wolfram.com/nb *) (* CreatedBy='Mathematica 7.0' *) (*CacheID: 234*) (* Internal cache information: NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 145, 7] NotebookDataLength[ 17719, 556] NotebookOptionsPosition[ 16570, 520] NotebookOutlinePosition[ 16946, 536] CellTagsIndexPosition[ 16903, 533] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[BoxData[ RowBox[{"vysl", "=", RowBox[{ RowBox[{"T", "[", "x", "]"}], "/.", RowBox[{ RowBox[{"DSolve", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ RowBox[{"T", "''"}], "[", "x", "]"}], "\[Equal]", "0"}], ",", RowBox[{ RowBox[{"T", "[", "0", "]"}], "\[Equal]", "T1"}], ",", RowBox[{ RowBox[{"T", "[", "d", "]"}], "\[Equal]", "T2"}]}], "}"}], ",", RowBox[{"T", "[", "x", "]"}], ",", "x"}], "]"}], "[", RowBox[{"[", "1", "]"}], "]"}]}]}]], "Input", CellChangeTimes->{{3.606116208296875*^9, 3.606116226046875*^9}, { 3.606116256546875*^9, 3.60611630746875*^9}}], Cell[BoxData[ RowBox[{"Needs", "[", "\"\\"", "]"}]], "Input"], Cell[BoxData[ RowBox[{"SetCoordinates", "[", RowBox[{"Cartesian", "[", RowBox[{"x", ",", "y", ",", "z"}], "]"}], "]"}]], "Input", CellChangeTimes->{{3.6061164868125*^9, 3.60611649865625*^9}}], Cell[BoxData[ RowBox[{"Grad", "[", "vysl", "]"}]], "Input", CellChangeTimes->{{3.60611645846875*^9, 3.60611646625*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"\[Lambda]", "*", FractionBox[ RowBox[{"T1", "-", "T2"}], "d"]}], "/.", RowBox[{"{", RowBox[{ RowBox[{"d", "\[Rule]", "0.6"}], ",", RowBox[{"\[Lambda]", "\[Rule]", "0.8"}], ",", RowBox[{"T1", "\[Rule]", "15"}], ",", RowBox[{"T2", "\[Rule]", RowBox[{"-", "15"}]}]}], "}"}]}]], "Input", CellChangeTimes->{{3.606117061859375*^9, 3.60611712915625*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"\[Lambda]", "*", FractionBox[ RowBox[{"T1", "-", "T2"}], "d"]}], "/.", RowBox[{"{", RowBox[{ RowBox[{"d", "\[Rule]", "0.1"}], ",", RowBox[{"\[Lambda]", "\[Rule]", "0.04"}], ",", RowBox[{"T1", "\[Rule]", "15"}], ",", RowBox[{"T2", "\[Rule]", RowBox[{"-", "15"}]}]}], "}"}]}]], "Input", CellChangeTimes->{{3.60611715390625*^9, 3.606117156328125*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{"Rth1", "=", FractionBox["0.45", RowBox[{"0.8", "*", "100"}]]}], ";"}], "\[IndentingNewLine]", RowBox[{"Print", "[", RowBox[{ "\"\\"", ",", "Rth1", ",", "\"\< \!\(\*FractionBox[\(K\), \(\(\\\ \\\ \)\(W\)\)]\)\>\""}], "]"}]}], "Input", CellChangeTimes->{{3.6061176654375*^9, 3.606117666734375*^9}, { 3.6061177085*^9, 3.60611778078125*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{"Rth2", "=", FractionBox["0.12", RowBox[{".04", "*", "100"}]]}], ";"}], "\[IndentingNewLine]", RowBox[{"Print", "[", RowBox[{ "\"\\"", ",", "Rth2", ",", "\"\< \!\(\*FractionBox[\(K\), \(\(\\\ \\\ \)\(W\)\)]\)\>\""}], "]"}]}], "Input", CellChangeTimes->{{3.6061177878125*^9, 3.60611780146875*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{"Q", "=", FractionBox[ RowBox[{"15", "-", RowBox[{"(", RowBox[{"-", "15"}], ")"}]}], RowBox[{"Rth1", "+", "Rth2"}]]}], ";"}], "\[IndentingNewLine]", RowBox[{"Print", "[", RowBox[{"\"\\"", ",", "Q", ",", "\"\< W\>\""}], "]"}]}], "Input", CellChangeTimes->{{3.606117806265625*^9, 3.60611785425*^9}}], Cell[BoxData[ RowBox[{"Solve", "[", RowBox[{ FractionBox[ RowBox[{"15", "-", "Tx"}], "Rth1"], "\[Equal]", "Q"}], "]"}]], "Input", CellChangeTimes->{{3.60611789396875*^9, 3.606117909578125*^9}}], Cell[BoxData[ RowBox[{"SetCoordinates", "[", RowBox[{"Cylindrical", "[", RowBox[{"r", ",", "fi", ",", "z"}], "]"}], "]"}]], "Input", CellChangeTimes->{{3.60611887028125*^9, 3.60611888484375*^9}}], Cell[BoxData[ RowBox[{"0", "==", RowBox[{"Numerator", "[", RowBox[{"Div", "[", RowBox[{"Grad", "[", RowBox[{"T", "[", "r", "]"}], "]"}], "]"}], "]"}]}]], "Input", CellChangeTimes->{{3.6061188930625*^9, 3.60611890746875*^9}, { 3.606118963171875*^9, 3.6061189715*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"res", "=", RowBox[{ RowBox[{"T", "[", "r", "]"}], "/.", RowBox[{ RowBox[{"DSolve", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"0", "==", RowBox[{"Numerator", "[", RowBox[{"Div", "[", RowBox[{"Grad", "[", RowBox[{"T", "[", "r", "]"}], "]"}], "]"}], "]"}]}], ",", RowBox[{ RowBox[{"T", "[", "rmin", "]"}], "\[Equal]", "T1"}], ",", RowBox[{ RowBox[{"T", "[", "rmax", "]"}], "\[Equal]", "T2"}]}], "}"}], ",", RowBox[{"T", "[", "r", "]"}], ",", "r"}], "]"}], "[", RowBox[{"[", "1", "]"}], "]"}]}]}], "\[IndentingNewLine]", "\[IndentingNewLine]"}]], "Input", CellChangeTimes->{{3.606119003828125*^9, 3.60611912140625*^9}, { 3.606119175265625*^9, 3.606119175265625*^9}}], Cell[BoxData[ RowBox[{"Simplify", "[", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", "\[Lambda]"}], "*", RowBox[{ RowBox[{"Grad", "[", "res", "]"}], "[", RowBox[{"[", "1", "]"}], "]"}]}], "/.", RowBox[{"r", "\[Rule]", "rmin"}]}], ")"}], "*", "2", "*", "Pi", "*", "rmin", "*", "delka"}], "]"}]], "Input", CellChangeTimes->{{3.606119179546875*^9, 3.6061191943125*^9}, { 3.6061192298125*^9, 3.606119262125*^9}}], Cell[BoxData[ RowBox[{"Plot", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"Log", "[", "ksi", "]"}], ",", RowBox[{"0.05", RowBox[{"(", RowBox[{ SuperscriptBox[ RowBox[{"(", "ksi", ")"}], "2"], "-", "1"}], ")"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{"ksi", ",", "1", ",", "10"}], "}"}], ",", RowBox[{"AxesLabel", "\[Rule]", RowBox[{"{", "\"\<\!\(\*FractionBox[\(r2\), \(r1\)]\)\>\"", "}"}]}]}], "]"}]], "Input", CellChangeTimes->{{3.606119631921875*^9, 3.606119658671875*^9}, { 3.606119831421875*^9, 3.60611984246875*^9}, {3.6061198856875*^9, 3.606119921796875*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{"c", "=", "4186."}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"\[CapitalDelta]T", "=", "36"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"msteckou", "=", "0.3"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"P", "=", RowBox[{"msteckou", "*", "c", "*", "\[CapitalDelta]T"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"Print", "[", RowBox[{"\"\\"", ",", RowBox[{ SuperscriptBox["10", RowBox[{"-", "3"}]], "P"}], ",", "\"\< kW\>\""}], "]"}]}], "Input", CellChangeTimes->{{3.60612038171875*^9, 3.606120402734375*^9}, { 3.60612044384375*^9, 3.6061204804375*^9}, {3.606120566234375*^9, 3.60612058440625*^9}}], Cell[BoxData[ RowBox[{"(*", "dale", "*)"}]], "Input", CellChangeTimes->{{3.60672044019871*^9, 3.60672044360496*^9}}], Cell[BoxData[ RowBox[{"\[Gamma]", "=", RowBox[{"Fit", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"10", ",", FractionBox["1", "35"]}], "}"}], ",", RowBox[{"{", RowBox[{"50", ",", FractionBox["1", "10"]}], "}"}]}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", RowBox[{"T", "[", "x", "]"}]}], "}"}], ",", RowBox[{"T", "[", "x", "]"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.6067216956205845`*^9, 3.6067217343393345`*^9}, { 3.6067221758080845`*^9, 3.60672217913621*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{"mm", "=", "0.001"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"cm", "=", RowBox[{"10", "mm"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"m", "=", RowBox[{"15", "cm"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"d", "=", RowBox[{"3", "mm"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Vsteckou", "=", RowBox[{"0.3", "*", SuperscriptBox["10", RowBox[{"-", "3"}]]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"vx", "=", FractionBox["Vsteckou", RowBox[{"d", "*", "m"}]]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"\[Rho]", "=", "1000"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"cp", "=", "4186"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"u", "=", "230"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"rce", "=", RowBox[{ RowBox[{"\[Rho]", "*", "cp", "*", RowBox[{ RowBox[{"T", "'"}], "[", "x", "]"}]}], "\[Equal]", RowBox[{"\[Gamma]", "*", SuperscriptBox[ RowBox[{"(", FractionBox["u", "d"], ")"}], "2"]}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"poc", "=", RowBox[{ RowBox[{"T", "[", "0", "]"}], "\[Equal]", "4"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"tepl", "=", RowBox[{"T", "/.", RowBox[{ RowBox[{"NDSolve", "[", RowBox[{ RowBox[{"{", RowBox[{"rce", ",", "poc"}], "}"}], ",", "T", ",", RowBox[{"{", RowBox[{"x", ",", "0", ",", "1"}], "}"}]}], "]"}], "[", RowBox[{"[", "1", "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Tout", "=", "40"}], ";"}], "\[IndentingNewLine]", RowBox[{"Plot", "[", RowBox[{ RowBox[{"{", RowBox[{"Tout", ",", RowBox[{"tepl", "[", "x", "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"x", ",", "0", ",", "1"}], "}"}]}], "]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"xmax", "=", RowBox[{"delka", "/.", RowBox[{"FindRoot", "[", RowBox[{ RowBox[{ RowBox[{"tepl", "[", "delka", "]"}], "\[Equal]", "Tout"}], ",", RowBox[{"{", RowBox[{"delka", ",", "0.1", ",", "1"}], "}"}]}], "]"}]}]}], ";"}]}], "Input", CellChangeTimes->{{3.60672176247996*^9, 3.60672178932371*^9}, { 3.60672182660496*^9, 3.60672185669871*^9}, {3.6067220968393345`*^9, 3.6067221045893345`*^9}, {3.6067221368393345`*^9, 3.60672238841746*^9}, { 3.6067224930893345`*^9, 3.6067225044955845`*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{"poc", "=", RowBox[{ RowBox[{"T", "[", "0", "]"}], "\[Equal]", "12"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"tepl", "=", RowBox[{"T", "/.", RowBox[{ RowBox[{"NDSolve", "[", RowBox[{ RowBox[{"{", RowBox[{"rce", ",", "poc"}], "}"}], ",", "T", ",", RowBox[{"{", RowBox[{"x", ",", "0", ",", "1"}], "}"}]}], "]"}], "[", RowBox[{"[", "1", "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{"tepl", "[", "xmax", "]"}]}], "Input", CellChangeTimes->{{3.6067225699330845`*^9, 3.60672258844871*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{"\[Gamma]MB", "=", RowBox[{"Fit", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"20", ",", FractionBox["1", "72"]}], "}"}], ",", RowBox[{"{", RowBox[{"75", ",", FractionBox["1", "30"]}], "}"}]}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", RowBox[{"T", "[", "x", "]"}]}], "}"}], ",", RowBox[{"T", "[", "x", "]"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"rce", "=", RowBox[{ RowBox[{"\[Rho]", "*", "cp", "*", RowBox[{ RowBox[{"T", "'"}], "[", "x", "]"}]}], "\[Equal]", RowBox[{"\[Gamma]MB", "*", SuperscriptBox[ RowBox[{"(", FractionBox["u", "d"], ")"}], "2"]}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"poc", "=", RowBox[{ RowBox[{"T", "[", "0", "]"}], "\[Equal]", "4"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"tepl", "=", RowBox[{"T", "/.", RowBox[{ RowBox[{"NDSolve", "[", RowBox[{ RowBox[{"{", RowBox[{"rce", ",", "poc"}], "}"}], ",", "T", ",", RowBox[{"{", RowBox[{"x", ",", "0", ",", "1"}], "}"}]}], "]"}], "[", RowBox[{"[", "1", "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{"tepl", "[", "xmax", "]"}]}], "Input", CellChangeTimes->{{3.6067226896205845`*^9, 3.60672274551121*^9}}], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"Remove", "[", "\"\\"", "]"}], ";"}], "\[IndentingNewLine]"}]], "Input", CellChangeTimes->{3.60672355519871*^9}], Cell[BoxData[{ RowBox[{ RowBox[{"Tpoc", "=", "5"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"hodina", "=", "3600"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Tin", "[", "t_", "]"}], ":=", RowBox[{"5", "+", RowBox[{"15", "*", RowBox[{"Tanh", "[", FractionBox["t", RowBox[{"2", "*", "hodina"}]], "]"}]}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Tout", "=", RowBox[{"-", "2"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"a", "=", RowBox[{"5.2", "*", SuperscriptBox["10", RowBox[{"-", "7"}]]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"\[Alpha]out", "=", "20"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"\[Alpha]in", "=", "5"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"d", "=", "0.6"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"\[Lambda]", "=", "0.8"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"rce", "=", RowBox[{ RowBox[{"D", "[", RowBox[{ RowBox[{"T", "[", RowBox[{"x", ",", "t"}], "]"}], ",", "t"}], "]"}], "\[Equal]", RowBox[{"a", "*", RowBox[{"D", "[", RowBox[{ RowBox[{"T", "[", RowBox[{"x", ",", "t"}], "]"}], ",", RowBox[{"{", RowBox[{"x", ",", "2"}], "}"}]}], "]"}]}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"okrout", "=", RowBox[{ RowBox[{ RowBox[{ RowBox[{"-", "\[Lambda]"}], "*", RowBox[{"D", "[", RowBox[{ RowBox[{"T", "[", RowBox[{"x", ",", "t"}], "]"}], ",", "x"}], "]"}]}], "\[Equal]", RowBox[{ RowBox[{"Tanh", "[", "t", "]"}], "*", "\[Alpha]out", "*", RowBox[{"(", RowBox[{ RowBox[{"T", "[", RowBox[{"x", ",", "t"}], "]"}], "-", "Tout"}], ")"}]}]}], "/.", RowBox[{"x", "\[Rule]", "d"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"okrin", "=", RowBox[{ RowBox[{ RowBox[{ RowBox[{"-", "\[Lambda]"}], "*", RowBox[{"D", "[", RowBox[{ RowBox[{"T", "[", RowBox[{"x", ",", "t"}], "]"}], ",", "x"}], "]"}]}], "\[Equal]", RowBox[{ RowBox[{"Tanh", "[", "t", "]"}], "*", "\[Alpha]in", "*", RowBox[{"(", RowBox[{ RowBox[{"Tin", "[", "t", "]"}], "-", RowBox[{"T", "[", RowBox[{"x", ",", "t"}], "]"}]}], ")"}]}]}], "/.", RowBox[{"x", "\[Rule]", "0"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"poc", "=", RowBox[{ RowBox[{"T", "[", RowBox[{"x", ",", "0"}], "]"}], "\[Equal]", "Tpoc"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"dohr", "=", RowBox[{"{", RowBox[{"rce", ",", "okrin", ",", "okrout", ",", "poc"}], "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"den", "=", RowBox[{"24", "*", "hodina"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"tmax", "=", RowBox[{"2", "*", "den"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"res", "=", RowBox[{ RowBox[{"NDSolve", "[", RowBox[{"dohr", ",", RowBox[{"T", "[", RowBox[{"x", ",", "t"}], "]"}], ",", RowBox[{"{", RowBox[{"x", ",", "0", ",", "d"}], "}"}], ",", RowBox[{"{", RowBox[{"t", ",", "0", ",", "tmax"}], "}"}]}], "]"}], "[", RowBox[{"[", "1", "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"Plot3D", "[", RowBox[{ RowBox[{ RowBox[{"T", "[", RowBox[{"x", ",", "t"}], "]"}], "/.", "res"}], ",", RowBox[{"{", RowBox[{"x", ",", "0", ",", "d"}], "}"}], ",", RowBox[{"{", RowBox[{"t", ",", "0", ",", "tmax"}], "}"}]}], "]"}], "\[IndentingNewLine]"}], "Input", CellChangeTimes->{{3.6067235595268345`*^9, 3.6067236139643345`*^9}, { 3.60672398357371*^9, 3.60672399951121*^9}, {3.6067243544018345`*^9, 3.6067243861830845`*^9}, 3.60672442788621*^9, {3.60672450622996*^9, 3.6067245504955845`*^9}, {3.60672460947996*^9, 3.6067247963393345`*^9}, { 3.60672484688621*^9, 3.6067249293705845`*^9}, {3.6071432375073223`*^9, 3.6071432876521907`*^9}}] }, WindowSize->{1003, 628}, WindowMargins->{{0, Automatic}, {Automatic, 0}}, Magnification->2.200000047683716, FrontEndVersion->"8.0 for Microsoft Windows (64-bit) (October 6, 2011)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[545, 20, 683, 19, 178, "Input"], Cell[1231, 41, 80, 1, 63, "Input"], Cell[1314, 44, 200, 4, 63, "Input"], Cell[1517, 50, 121, 2, 63, "Input"], Cell[1641, 54, 423, 12, 93, "Input"], Cell[2067, 68, 424, 12, 93, "Input"], Cell[2494, 82, 408, 11, 165, "Input"], Cell[2905, 95, 364, 10, 165, "Input"], Cell[3272, 107, 374, 10, 134, "Input"], Cell[3649, 119, 205, 5, 93, "Input"], Cell[3857, 126, 204, 4, 63, "Input"], Cell[4064, 132, 288, 7, 63, "Input"], Cell[4355, 141, 859, 23, 254, "Input"], Cell[5217, 166, 473, 13, 102, "Input"], Cell[5693, 181, 646, 18, 139, "Input"], Cell[6342, 201, 699, 18, 222, "Input"], Cell[7044, 221, 118, 2, 63, "Input"], Cell[7165, 225, 572, 17, 93, "Input"], Cell[7740, 244, 2502, 75, 703, "Input"], Cell[10245, 321, 618, 18, 140, "Input"], Cell[10866, 341, 1413, 44, 281, "Input"], Cell[12282, 387, 173, 5, 102, "Input"], Cell[12458, 394, 4108, 124, 1011, "Input"] } ] *) (* End of internal cache information *)