(* Content-type: application/vnd.wolfram.mathematica *) (*** Wolfram Notebook File ***) (* http://www.wolfram.com/nb *) (* CreatedBy='Mathematica 10.0' *) (*CacheID: 234*) (* Internal cache information: NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 158, 7] NotebookDataLength[ 13601, 415] NotebookOptionsPosition[ 12889, 390] NotebookOutlinePosition[ 13266, 406] CellTagsIndexPosition[ 13223, 403] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[BoxData[ RowBox[{"(*", RowBox[{"Tepelny", " ", "odpor", " ", "valce"}], "*)"}]], "Input", CellChangeTimes->{{3.670225639903417*^9, 3.6702256519484367`*^9}}], Cell[BoxData[ RowBox[{"Div", "[", RowBox[{ RowBox[{"Grad", "[", RowBox[{ RowBox[{"T", "[", "r", "]"}], ",", " ", RowBox[{"{", RowBox[{"r", ",", "fi", ",", "z"}], "}"}], ",", " ", "\"\\""}], "]"}], ",", " ", RowBox[{"{", RowBox[{"r", ",", "fi", ",", "z"}], "}"}], ",", " ", "\"\\""}], "]"}]], "Input", CellChangeTimes->{{3.6702246181741548`*^9, 3.6702246306660776`*^9}, { 3.670224662686195*^9, 3.6702247022039557`*^9}}], Cell[BoxData[ RowBox[{"resT", "=", RowBox[{ RowBox[{"T", "[", "r", "]"}], "/.", RowBox[{ RowBox[{"DSolve", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"0", "==", RowBox[{"Div", "[", RowBox[{ RowBox[{"Grad", "[", RowBox[{ RowBox[{"T", "[", "r", "]"}], ",", " ", RowBox[{"{", RowBox[{"r", ",", "fi", ",", "z"}], "}"}], ",", " ", "\"\\""}], "]"}], ",", " ", RowBox[{"{", RowBox[{"r", ",", "fi", ",", "z"}], "}"}], ",", " ", "\"\\""}], "]"}]}], ",", RowBox[{ RowBox[{"T", "[", "r1", "]"}], "\[Equal]", "T1"}], ",", RowBox[{ RowBox[{"T", "[", "r2", "]"}], "\[Equal]", "T2"}]}], "}"}], ",", RowBox[{"T", "[", "r", "]"}], ",", "r"}], "]"}], "[", RowBox[{"[", "1", "]"}], "]"}]}]}]], "Input", CellChangeTimes->{{3.6702248306980085`*^9, 3.6702248646452465`*^9}, { 3.6702249008711247`*^9, 3.670224951699763*^9}}], Cell[BoxData[ RowBox[{"Simplify", "[", RowBox[{ RowBox[{"-", "\[Lambda]"}], "*", RowBox[{ RowBox[{"Grad", "[", RowBox[{"resT", ",", " ", RowBox[{"{", RowBox[{"r", ",", "fi", ",", "z"}], "}"}], ",", " ", "\"\\""}], "]"}], "[", RowBox[{"[", "1", "]"}], "]"}], "*", "l", "*", "2", "Pi", "*", "r"}], "]"}]], "Input", CellChangeTimes->{{3.6702250686744223`*^9, 3.670225147576508*^9}}], Cell[BoxData[ RowBox[{"(*", RowBox[{"Prutokovy", " ", "ohrivac"}], "*)"}]], "Input", CellChangeTimes->{{3.6702256778700724`*^9, 3.670225679553935*^9}, { 3.670225732730218*^9, 3.6702257377837243`*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{"msteckou", "=", "0.25"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"\[Rho]", "=", "1000"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"cp", "=", "4186"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Vsteckou", "=", FractionBox["msteckou", "\[Rho]"]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"T1", "=", "15"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"T2", "=", "45"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"m", "=", "0.5"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"cm", "=", "0.01"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"d", "=", RowBox[{"0.5", "cm"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"prikon", "=", RowBox[{"msteckou", "*", "cp", "*", RowBox[{"(", RowBox[{"T2", "-", "T1"}], ")"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"vx", "=", FractionBox["Vsteckou", RowBox[{"d", "*", "m"}]]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"U", "=", "230"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Ecko", "=", FractionBox["U", "d"]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"resab", "=", RowBox[{ RowBox[{"Solve", "[", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ RowBox[{"20", "*", "a"}], "+", "b"}], "==", FractionBox["1", "70."]}], ",", RowBox[{ RowBox[{ RowBox[{"75", "*", "a"}], "+", "b"}], "==", FractionBox["1", "30."]}]}], "}"}], "]"}], "[", RowBox[{"[", "1", "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"(*", RowBox[{"resab", "=", RowBox[{ RowBox[{"Solve", "[", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ RowBox[{"10", "*", "a"}], "+", "b"}], "==", FractionBox["1", "35."]}], ",", RowBox[{ RowBox[{ RowBox[{"50", "*", "a"}], "+", "b"}], "==", FractionBox["1", "10."]}]}], "}"}], "]"}], "[", RowBox[{"[", "1", "]"}], "]"}]}], "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"rce", "=", RowBox[{ RowBox[{ RowBox[{"\[Rho]", "*", "cp", "*", "vx", "*", RowBox[{ RowBox[{"T", "'"}], "[", "x", "]"}]}], "\[Equal]", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"a", "*", RowBox[{"T", "[", "x", "]"}]}], "+", "b"}], ")"}], "*", SuperscriptBox["Ecko", "2"]}]}], "/.", "resab"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"okr", "=", RowBox[{ RowBox[{"T", "[", "0", "]"}], "\[Equal]", "T1"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"res", "=", RowBox[{ RowBox[{"DSolve", "[", RowBox[{ RowBox[{"{", RowBox[{"rce", ",", "okr"}], "}"}], ",", RowBox[{"T", "[", "x", "]"}], ",", "x"}], "]"}], "[", RowBox[{"[", "1", "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"teplota", "=", RowBox[{ RowBox[{"T", "[", "x", "]"}], "/.", "res"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"Plot", "[", RowBox[{ RowBox[{"{", RowBox[{"teplota", ",", "T2"}], "}"}], ",", RowBox[{"{", RowBox[{"x", ",", "0", ",", "0.6"}], "}"}], ",", RowBox[{"AxesLabel", "->", RowBox[{"{", RowBox[{"\"\\"", ",", "\"\\""}], "}"}]}], ",", RowBox[{"PlotStyle", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"Thickness", "[", "0.008", "]"}], ",", RowBox[{"Thickness", "[", "0.008", "]"}]}], "}"}]}]}], "]"}]}], "Input", CellChangeTimes->{ 3.6702259738694324`*^9, {3.6702264682811117`*^9, 3.670226532245496*^9}, { 3.670226593819686*^9, 3.67022669420577*^9}, 3.670226795923707*^9, { 3.670226889318979*^9, 3.6702269089902697`*^9}, {3.670226956811587*^9, 3.670226957033942*^9}, 3.670226988944729*^9, {3.6702270711533813`*^9, 3.6702270861790085`*^9}, {3.670227175907874*^9, 3.670227198434066*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{"Rth1", "=", RowBox[{ FractionBox["1", RowBox[{"5", "*", "100"}]], "+", FractionBox["0.4", RowBox[{"0.5", "*", "100"}]]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Rth2", "=", FractionBox["1", RowBox[{"10", "*", "100"}]]}], ";"}], "\[IndentingNewLine]", RowBox[{"Tepl", "=", RowBox[{ RowBox[{ FractionBox["Rth2", RowBox[{"Rth1", "+", "Rth2"}]], "*", RowBox[{"(", RowBox[{"20", "-", "0"}], ")"}]}], "+", "0"}]}]}], "Input", CellChangeTimes->{{3.6702280093786016`*^9, 3.670228088754363*^9}}], Cell[BoxData[ RowBox[{"Plot", "[", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ FractionBox["1", RowBox[{"alfa", "*", "2", "*", "Pi", "*", "r2", "*", "l"}]], "+", RowBox[{ FractionBox["1", RowBox[{"2", "Pi", "*", "\[Lambda]", "*", "l"}]], RowBox[{"Log", "[", FractionBox["r2", "r1"], "]"}]}]}], ")"}], "/.", RowBox[{"{", RowBox[{ RowBox[{"alfa", "\[Rule]", "5."}], ",", RowBox[{"l", "\[Rule]", "1"}], ",", RowBox[{"\[Lambda]", "\[Rule]", "0.5"}], ",", RowBox[{"r1", "\[Rule]", "0.001"}]}], "}"}]}], ",", RowBox[{"{", RowBox[{"r2", ",", "0", ",", "1"}], "}"}]}], "]"}]], "Input"], Cell[BoxData[{ RowBox[{ RowBox[{"ro", "=", "2400"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Cp", "=", "1011"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Tout", "=", RowBox[{"5", "+", "273"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"ToutS", "=", RowBox[{ RowBox[{"-", "5"}], "+", "273"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Tin", "=", RowBox[{"20", "+", "273"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"sigma", "=", RowBox[{"5.775", "*", SuperscriptBox["10", RowBox[{"-", "8"}]]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"epsilon", "=", "0.93"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"alfa", "=", "15"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"lambda", "=", "0.3"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"tloustka", "=", "0.3"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"tMax", "=", "200000"}], ";"}], "\[IndentingNewLine]"}], "Input", CellChangeTimes->{{3.5612903967481766`*^9, 3.5612905468284235`*^9}, { 3.5612906948522825`*^9, 3.56129073056474*^9}, {3.561290890171415*^9, 3.5612908975190277`*^9}, {3.5612910866295443`*^9, 3.561291189544903*^9}, 3.561291247032189*^9, {3.5613498341450305`*^9, 3.5613498344500475`*^9}, { 3.56136316112529*^9, 3.561363187162779*^9}, {3.5613632253789654`*^9, 3.561363259371909*^9}, {3.5613667382598906`*^9, 3.56136673947696*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{"rce", "=", RowBox[{ RowBox[{"ro", "*", "Cp", "*", RowBox[{"D", "[", RowBox[{ RowBox[{"T", "[", RowBox[{"x", ",", "t"}], "]"}], ",", "t"}], "]"}]}], "\[Equal]", RowBox[{"lambda", "*", RowBox[{"D", "[", RowBox[{ RowBox[{"D", "[", RowBox[{ RowBox[{"T", "[", RowBox[{"x", ",", "t"}], "]"}], ",", "x"}], "]"}], ",", "x"}], "]"}]}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"pocPodminka", "=", RowBox[{ RowBox[{"T", "[", RowBox[{"x", ",", "0"}], "]"}], "\[Equal]", "Tout"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"okrVlevo", "=", RowBox[{ RowBox[{"T", "[", RowBox[{"0", ",", "t"}], "]"}], "\[Equal]", RowBox[{ RowBox[{"Tin", "*", RowBox[{"Tanh", "[", "t", "]"}]}], "+", RowBox[{"Tout", "*", RowBox[{"(", RowBox[{"1", "-", RowBox[{"Tanh", "[", "t", "]"}]}], ")"}]}]}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"okrVpravo", "=", RowBox[{ RowBox[{ RowBox[{ RowBox[{"-", "lambda"}], "*", RowBox[{"D", "[", RowBox[{ RowBox[{"T", "[", RowBox[{"x", ",", "t"}], "]"}], ",", "x"}], "]"}]}], "\[Equal]", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"alfa", "*", RowBox[{"(", RowBox[{ RowBox[{"T", "[", RowBox[{"x", ",", "t"}], "]"}], "-", "Tout"}], ")"}]}], "+", RowBox[{"sigma", "*", "epsilon", "*", RowBox[{"(", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"T", "[", RowBox[{"x", ",", "t"}], "]"}], ")"}], "4"], "-", SuperscriptBox[ RowBox[{"(", "ToutS", ")"}], "4"]}], ")"}]}]}], ")"}], "*", FractionBox["t", RowBox[{"t", "+", "1"}]]}]}], "/.", RowBox[{"x", "\[Rule]", "tloustka"}]}]}], ";"}]}], "Input", CellChangeTimes->{{3.5612905549248376`*^9, 3.561290874758588*^9}, { 3.5612910740003233`*^9, 3.5612911074435773`*^9}, {3.561365654*^9, 3.56136566790625*^9}, {3.56136571203125*^9, 3.561365766703125*^9}, { 3.561365799671875*^9, 3.561365815953125*^9}}], Cell[BoxData[ RowBox[{"res", "=", RowBox[{"NDSolve", "[", RowBox[{ RowBox[{"{", RowBox[{"rce", ",", "pocPodminka", ",", "okrVlevo", ",", "okrVpravo"}], "}"}], ",", RowBox[{"T", "[", RowBox[{"x", ",", "t"}], "]"}], ",", RowBox[{"{", RowBox[{"x", ",", "0", ",", "tloustka"}], "}"}], ",", RowBox[{"{", RowBox[{"t", ",", "0", ",", "tMax"}], "}"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.561290882309001*^9, 3.5612909997893972`*^9}}], Cell[BoxData[ RowBox[{"Plot3D", "[", RowBox[{ RowBox[{ RowBox[{ RowBox[{"T", "[", RowBox[{"x", ",", "t"}], "]"}], "-", "273"}], "/.", "res"}], ",", RowBox[{"{", RowBox[{"x", ",", "0", ",", "tloustka"}], "}"}], ",", RowBox[{"{", RowBox[{"t", ",", "0", ",", "tMax"}], "}"}]}], "]"}]], "Input", CellChangeTimes->{{3.5612910091378117`*^9, 3.561291047137076*^9}, { 3.5613631979713974`*^9, 3.561363199031458*^9}}] }, WindowSize->{1280, 685}, WindowMargins->{{-8, Automatic}, {Automatic, -8}}, Magnification:>2.4 Inherited, FrontEndVersion->"10.0 for Microsoft Windows (64-bit) (September 9, 2014)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[558, 20, 166, 3, 70, "Input"], Cell[727, 25, 498, 13, 114, "Input"], Cell[1228, 40, 1062, 28, 286, "Input"], Cell[2293, 70, 445, 12, 144, "Input"], Cell[2741, 84, 207, 4, 70, "Input"], Cell[2951, 90, 3982, 118, 1206, "Input"], Cell[6936, 210, 601, 19, 315, "Input"], Cell[7540, 231, 691, 20, 208, "Input"], Cell[8234, 253, 1450, 35, 576, "Input"], Cell[9687, 290, 2252, 69, 453, "Input"], Cell[11942, 361, 489, 13, 144, "Input"], Cell[12434, 376, 451, 12, 144, "Input"] } ] *) (* End of internal cache information *)