(* 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[ 9707, 300] NotebookOptionsPosition[ 9121, 278] NotebookOutlinePosition[ 9497, 294] CellTagsIndexPosition[ 9454, 291] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[BoxData[ RowBox[{"(*", RowBox[{"valcovy", " ", "odpor"}], "*)"}]], "Input", CellChangeTimes->{{3.670219545442359*^9, 3.670219553369952*^9}}], Cell[BoxData[ RowBox[{"r", "*", RowBox[{"Div", "[", RowBox[{ RowBox[{"Grad", "[", RowBox[{ RowBox[{"T", "[", "r", "]"}], ",", " ", RowBox[{"{", RowBox[{"r", ",", "fi", ",", "z"}], "}"}], ",", " ", "\"\\""}], "]"}], ",", " ", RowBox[{"{", RowBox[{"r", ",", "fi", ",", "z"}], "}"}], ",", " ", "\"\\""}], "]"}]}]], "Input", CellChangeTimes->{{3.670218323401511*^9, 3.670218396725069*^9}, { 3.670219103317941*^9, 3.6702191036863976`*^9}}], Cell[BoxData[{ RowBox[{"resT", "=", RowBox[{ RowBox[{"T", "[", "r", "]"}], "/.", RowBox[{ RowBox[{"DSolve", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"0", "\[Equal]", RowBox[{"r", "*", 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", "]"}], "]"}]}]}], "\n", RowBox[{"Simplify", "[", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", "\[Lambda]"}], "*", RowBox[{ RowBox[{"Grad", "[", RowBox[{"resT", ",", " ", RowBox[{"{", RowBox[{"r", ",", "fi", ",", "z"}], "}"}], ",", " ", "\"\\""}], "]"}], "[", RowBox[{"[", "1", "]"}], "]"}]}], "/.", RowBox[{"r", "\[Rule]", "rmin"}]}], ")"}], "*", "2", "*", "Pi", "*", "rmin", "*", "l"}], "]"}]}], "Input", CellChangeTimes->{{3.6702184634956775`*^9, 3.670218517761662*^9}, { 3.670218567406943*^9, 3.6702185938448076`*^9}, {3.6702187085525646`*^9, 3.670218711206176*^9}, 3.670218859240677*^9, {3.670219134527495*^9, 3.670219135249816*^9}}], Cell[BoxData[ RowBox[{"Plot", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"Log", "[", "ksi", "]"}], ",", RowBox[{"0.05", "*", RowBox[{"(", RowBox[{ SuperscriptBox["ksi", "2"], "+", "1"}], ")"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{"ksi", ",", "1", ",", "10"}], "}"}]}], "]"}]], "Input", CellChangeTimes->{{3.670219344691554*^9, 3.6702194098658915`*^9}}], Cell[BoxData[ RowBox[{"(*", RowBox[{"prutokovy", " ", "ohrivac"}], "*)"}]], "Input", CellChangeTimes->{{3.6702195401063766`*^9, 3.670219574562686*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{"msteckou", "=", "0.2"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"\[Rho]", "=", "1000"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"cp", "=", "4186"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Vsteckou", "=", FractionBox["msteckou", "\[Rho]"]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"T1", "=", "4"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"T2", "=", "45"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"prikon", "=", RowBox[{"msteckou", "*", "cp", "*", RowBox[{"(", RowBox[{"T2", "-", "T1"}], ")"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"vx", "=", FractionBox["Vsteckou", RowBox[{"m", "*", "d"}]]}], ";"}], "\[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[{"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[{"U", "=", "230"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"m", "=", "0.5"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"cm", "=", "0.01"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"d", "=", RowBox[{"0.5", "cm"}]}], ";"}], "\n", RowBox[{ RowBox[{"Ecko", "=", FractionBox["U", "d"]}], ";"}], "\[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[{"poc", "=", RowBox[{ RowBox[{"T", "[", "0", "]"}], "\[Equal]", "T1"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"res", "=", RowBox[{ RowBox[{"DSolve", "[", RowBox[{ RowBox[{"{", RowBox[{"rce", ",", "poc"}], "}"}], ",", 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.5"}], "}"}], ",", RowBox[{"AxesLabel", "->", RowBox[{"{", RowBox[{"\"\\"", ",", "\"\\""}], "}"}]}], ",", RowBox[{"PlotStyle", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"Thickness", "[", "0.008", "]"}], ",", RowBox[{"Thickness", "[", "0.008", "]"}]}], "}"}]}]}], "]"}]}], "Input", CellChangeTimes->{ 3.6702163187762327`*^9, 3.6702200563450613`*^9, {3.6702201936732864`*^9, 3.670220197237046*^9}, 3.6702203973419476`*^9, {3.670220687105709*^9, 3.6702207376788387`*^9}, {3.670220788565041*^9, 3.6702207958265963`*^9}, { 3.670221173888477*^9, 3.670221187839735*^9}, {3.6702212417086706`*^9, 3.670221288862954*^9}, 3.6702264620634813`*^9}], Cell[BoxData[ RowBox[{"(*", "alfa", "*)"}]], "Input", CellChangeTimes->{{3.6702228487000985`*^9, 3.6702228515544257`*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{"Rth1", "=", RowBox[{ FractionBox["1", RowBox[{"5", "*", "50"}]], "+", FractionBox["0.4", RowBox[{"0.5", "*", "50"}]]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Rth2", "=", FractionBox["1", RowBox[{"10", "*", "50"}]]}], ";"}], "\[IndentingNewLine]", RowBox[{"Tepl", "=", RowBox[{ RowBox[{ FractionBox["Rth2", RowBox[{"Rth1", "+", "Rth2"}]], "*", RowBox[{"(", RowBox[{"20", "-", "0"}], ")"}]}], "+", "0"}]}]}], "Input", CellChangeTimes->{{3.6702221687021832`*^9, 3.670222250267418*^9}, { 3.670222285201171*^9, 3.6702222893703003`*^9}, {3.6702223517654285`*^9, 3.670222367082036*^9}, {3.670222418804846*^9, 3.6702224190677724`*^9}, 3.670819892650345*^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", CellChangeTimes->{{3.670222773706421*^9, 3.6702227822210307`*^9}}] }, WindowSize->{1280, 685}, WindowMargins->{{-8, Automatic}, {Automatic, 0}}, Magnification:>2.2 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, 150, 3, 66, "Input"], Cell[711, 25, 527, 14, 109, "Input"], Cell[1241, 41, 1729, 46, 403, "Input"], Cell[2973, 89, 410, 12, 72, "Input"], Cell[3386, 103, 156, 3, 66, "Input"], Cell[3545, 108, 3904, 117, 1090, "Input"], Cell[7452, 227, 124, 2, 66, "Input"], Cell[7579, 231, 776, 22, 271, "Input"], Cell[8358, 255, 759, 21, 167, "Input"] } ] *) (* End of internal cache information *)