(* Content-type: application/mathematica *) (*** Wolfram Notebook File ***) (* http://www.wolfram.com/nb *) (* CreatedBy='Mathematica 6.0' *) (*CacheID: 234*) (* Internal cache information: NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 145, 7] NotebookDataLength[ 62141, 1609] NotebookOptionsPosition[ 58978, 1496] NotebookOutlinePosition[ 59391, 1512] CellTagsIndexPosition[ 59348, 1509] WindowFrame->Normal ContainsDynamic->False*) (* Beginning of Notebook Content *) Notebook[{ Cell[CellGroupData[{ Cell["Reactor Kinetics", "Title", CellChangeTimes->{{3.393089088911792*^9, 3.393089092292838*^9}}], Cell["\<\ Nuclear Reactor Physics J.A. McNeil April 2008\ \>", "Text", CellChangeTimes->{{3.4173669863642883`*^9, 3.417367015513379*^9}}], Cell[CellGroupData[{ Cell["Prompt neutron lifetime", "Section", CellChangeTimes->{{3.393089099923243*^9, 3.393089106580451*^9}, { 3.393415895798856*^9, 3.393415898836241*^9}, {3.393635453992888*^9, 3.393635457056673*^9}}], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"vavg", "=", RowBox[{"2.200", " ", RowBox[{"10", "^", "5"}]}]}], " ", RowBox[{"(*", " ", RowBox[{"cm", "/", "s"}], " ", "*)"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"\[Sigma]ScattU", "=", ".397"}], ";", RowBox[{"\[Sigma]AbsU238", "=", ".363"}], ";", RowBox[{"\[Nu]value", "=", "2.07"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"\[Sigma]AbsU235", "=", "687"}], ";", RowBox[{"\[Sigma]FissU", "=", "587"}], ";", RowBox[{"PercentU235", "=", ".0072"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"\[Sigma]AbsNet", "=", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"1", "-", "PercentU235"}], ")"}], " ", "\[Sigma]AbsU238"}], "+", RowBox[{"PercentU235", " ", "\[Sigma]AbsU235"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"\[Sigma]Fiss", "=", RowBox[{"PercentU235", " ", "\[Sigma]FissU"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"PromptLifetime", "=", RowBox[{ RowBox[{"1", "/", RowBox[{"(", RowBox[{"vavg", " ", "\[Sigma]AbsNet"}], ")"}]}], " ", RowBox[{"10", "^", "6"}]}]}], " ", ";", RowBox[{"(*", " ", RowBox[{"micro", " ", "sec"}], "*)"}], "\[IndentingNewLine]", RowBox[{"\[CapitalLambda]value", "=", RowBox[{ RowBox[{"1", "/", RowBox[{"(", RowBox[{"\[Nu]value", " ", "vavg", " ", "\[Sigma]Fiss"}], " ", ")"}]}], " ", RowBox[{"10", "^", "6"}]}]}], " ", ";", RowBox[{"(*", " ", RowBox[{"micro", " ", "sec"}], "*)"}], "\[IndentingNewLine]", RowBox[{"Print", "[", RowBox[{ "\"\< Natural Uranium: \[Tau](prompt) = \>\"", ",", "PromptLifetime", ",", "\"\< \[Mu]s\>\""}], "]"}]}], "\[IndentingNewLine]", RowBox[{"Print", "[", RowBox[{ "\"\< Natural Uranium: \[CapitalLambda](prompt) = \>\"", ",", "\[CapitalLambda]value", ",", "\"\< \[Mu]s\>\""}], "]"}]}], "Input", CellChangeTimes->{{3.393089129882485*^9, 3.393089167788682*^9}, { 3.3930892060690536`*^9, 3.3930892583736677`*^9}, {3.39308943017686*^9, 3.3930894328737383`*^9}, {3.3930896658750677`*^9, 3.3930896918475533`*^9}, {3.393091066368746*^9, 3.393091169887536*^9}, { 3.393091235968754*^9, 3.393091357778775*^9}, {3.39309168206713*^9, 3.393091687309807*^9}, {3.3930918464295883`*^9, 3.393091863588975*^9}, { 3.3930920848410597`*^9, 3.393092114209614*^9}, {3.393092150700025*^9, 3.393092160225504*^9}, {3.3934157656975927`*^9, 3.3934157674990683`*^9}, { 3.393415800587804*^9, 3.3934158685771103`*^9}, {3.393415911029256*^9, 3.393415912418807*^9}, {3.3935382065697117`*^9, 3.393538238041498*^9}, { 3.393635466991721*^9, 3.3936355130687723`*^9}, {3.3936355468776817`*^9, 3.3936355557257338`*^9}, 3.39363559722965*^9, {3.417184374928616*^9, 3.417184377747621*^9}, {3.4171879940985622`*^9, 3.417188067575696*^9}}], Cell[CellGroupData[{ Cell[BoxData[ InterpretationBox[ RowBox[{"\<\" Natural Uranium: \[Tau](prompt) = \"\>", "\[InvisibleSpace]", "0.856536178930161`", "\[InvisibleSpace]", "\<\" \[Mu]s\"\>"}], SequenceForm[ " Natural Uranium: \[Tau](prompt) = ", 0.856536178930161, " \[Mu]s"], Editable->False]], "Print", CellChangeTimes->{3.393635516450933*^9, 3.393635598555071*^9, 3.393672462448592*^9, 3.417184336820006*^9, 3.4171844134133463`*^9, 3.4171880725144978`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\" Natural Uranium: \[CapitalLambda](prompt) = \"\>", "\[InvisibleSpace]", "0.5195607990462693`", "\[InvisibleSpace]", "\<\" \[Mu]s\"\>"}], SequenceForm[ " Natural Uranium: \[CapitalLambda](prompt) = ", 0.5195607990462693, " \[Mu]s"], Editable->False]], "Print", CellChangeTimes->{3.393635516450933*^9, 3.393635598555071*^9, 3.393672462448592*^9, 3.417184336820006*^9, 3.4171844134133463`*^9, 3.4171880725184813`*^9}] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"\[CapitalSigma]aGSTR", "=", ".12"}], ";", RowBox[{"\[Tau]prompt", "=", RowBox[{ RowBox[{"1", "/", RowBox[{"(", RowBox[{"vavg", " ", "\[CapitalSigma]aGSTR"}], ")"}]}], " ", RowBox[{"10", "^", "6"}]}]}], " ", ";"}], RowBox[{"(*", " ", RowBox[{"micro", " ", "sec"}], "*)"}]}], "\[IndentingNewLine]", RowBox[{"Print", "[", RowBox[{ "\"\< GSTR: \[Tau](prompt) = \>\"", ",", "\[Tau]prompt", ",", "\"\< \[Mu]s\>\""}], "]"}]}], "Input", CellChangeTimes->{{3.393538250451482*^9, 3.393538305766047*^9}, { 3.393628252374881*^9, 3.393628254708067*^9}, {3.393628442630764*^9, 3.3936284443239737`*^9}, {3.393628480878504*^9, 3.393628482644068*^9}, { 3.393635562729759*^9, 3.3936356028359413`*^9}, {3.393672442855397*^9, 3.393672454436557*^9}}], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\" GSTR: \[Tau](prompt) = \"\>", "\[InvisibleSpace]", "37.878787878787875`", "\[InvisibleSpace]", "\<\" \[Mu]s\"\>"}], SequenceForm[" GSTR: \[Tau](prompt) = ", 37.878787878787875`, " \[Mu]s"], Editable->False]], "Print", CellChangeTimes->{{3.3936355855010653`*^9, 3.39363560344289*^9}, 3.3936724700789757`*^9, 3.4171843369558287`*^9, 3.4171844135089827`*^9}] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["Delayed neutrons - one delayed neutron group", "Section", CellChangeTimes->{{3.3930891134288683`*^9, 3.393089115588193*^9}, { 3.393415901098483*^9, 3.393415902148323*^9}, {3.3935413309091043`*^9, 3.393541337583991*^9}, 3.3935438745287447`*^9}], Cell[CellGroupData[{ Cell["Mean lifetime of delayed neutron precursors", "Subsection", CellChangeTimes->{{3.393635418312327*^9, 3.393635440664928*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{ RowBox[{"U235", " ", "-", " ", RowBox[{"List", " ", RowBox[{"format", ":", " ", RowBox[{"(", RowBox[{"\[Lambda]i", ",", "\[Beta]i"}], ")"}]}]}]}], ",", " ", RowBox[{ RowBox[{ RowBox[{"Data", " ", "from", " ", "Lamarsh"}], " ", "&"}], " ", "Barrata"}], ",", " ", RowBox[{"Table", " ", "3.5"}], ",", " ", RowBox[{"p", ".", " ", "88"}]}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"\[Lambda]DelayedList", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{".0124", ",", ".000215"}], "}"}], ",", RowBox[{"{", RowBox[{".0305", ",", ".001412"}], "}"}], ",", RowBox[{"{", RowBox[{".111", ",", ".001274"}], "}"}], ",", RowBox[{"{", RowBox[{".301", ",", ".002568"}], "}"}], ",", RowBox[{"{", RowBox[{"1.14", ",", ".000748"}], "}"}], ",", RowBox[{"{", RowBox[{"3.01", ",", ".000273"}], "}"}]}], "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"\[Tau]DelayedList", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"1", "/", RowBox[{"\[Lambda]DelayedList", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}]}], ",", RowBox[{"\[Lambda]DelayedList", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "6"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"\[Beta]Total", "=", RowBox[{"Sum", "[", RowBox[{ RowBox[{"\[Tau]DelayedList", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], " ", ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "6"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"\[Tau]Avg", "=", RowBox[{ RowBox[{"Sum", "[", RowBox[{ RowBox[{ RowBox[{"\[Tau]DelayedList", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], RowBox[{"\[Tau]DelayedList", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}]}], " ", ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "6"}], "}"}]}], "]"}], "/", "\[Beta]Total"}]}], ";"}], " ", "\[IndentingNewLine]", RowBox[{ RowBox[{"\[Lambda]Avg", "=", RowBox[{"1", "/", "\[Tau]Avg"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"Print", "[", RowBox[{ "\"\<| \[Beta] (total) = \>\"", ",", "\[Beta]Total", ",", "\"\< | \[Lambda] (avg) = \>\"", ",", "\[Lambda]Avg", ",", "\"\< s^-1 | \[Tau] (avg) = \>\"", ",", "\[Tau]Avg", ",", "\"\< s |\>\""}], "]"}]}]}]], "Input", CellChangeTimes->{{3.39353833097665*^9, 3.393538518251545*^9}, { 3.393543894609749*^9, 3.3935439139062223`*^9}, {3.393610553578022*^9, 3.393610575824182*^9}, {3.393610678226407*^9, 3.393610689750658*^9}, { 3.393612278863268*^9, 3.3936122859429693`*^9}, {3.393628969103298*^9, 3.393629087781427*^9}, {3.393629233780505*^9, 3.393629239441525*^9}, { 3.393630315024229*^9, 3.393630413609267*^9}, {3.393630479795793*^9, 3.393630482993236*^9}, {3.3936307506304893`*^9, 3.3936309235969477`*^9}, { 3.3936310062621717`*^9, 3.393631383111895*^9}, {3.417184669630856*^9, 3.417184674869697*^9}, {3.417184739884676*^9, 3.4171848291210613`*^9}, { 3.4171849028849783`*^9, 3.417184955046471*^9}, {3.417184999113188*^9, 3.4171849996787167`*^9}, {3.417185041549592*^9, 3.417185205885105*^9}}], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"| \[Beta] (total) = \"\>", "\[InvisibleSpace]", "0.00649`", "\[InvisibleSpace]", "\<\" | \[Lambda] (avg) = \"\>", "\[InvisibleSpace]", "0.07690514854436212`", "\[InvisibleSpace]", "\<\" s^-1 | \[Tau] (avg) = \"\>", "\[InvisibleSpace]", "13.003030602342026`", "\[InvisibleSpace]", "\<\" s |\"\>"}], SequenceForm[ "| \[Beta] (total) = ", 0.00649, " | \[Lambda] (avg) = ", 0.07690514854436212, " s^-1 | \[Tau] (avg) = ", 13.003030602342026`, " s |"], Editable->False]], "Print", CellChangeTimes->{ 3.393630896135035*^9, 3.3936313235049*^9, {3.3936313675232553`*^9, 3.393631384006197*^9}, 3.3936725569214897`*^9, 3.417184337562768*^9, 3.417184831149588*^9, 3.417184956454279*^9, 3.417185000672927*^9, { 3.417185086632812*^9, 3.417185141847331*^9}, {3.417185200116458*^9, 3.4171852086005907`*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{ "Data", " ", "from", " ", "Measurement", " ", "of", " ", "Reactor", " ", "Period", " ", "handout"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"\[Beta]value", "=", ".00642"}], ";"}], "\[IndentingNewLine]", RowBox[{"\[Tau]DelayedList", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"55.72", ",", ".00021"}], "}"}], ",", RowBox[{"{", RowBox[{"14.1", ",", ".00141"}], "}"}], ",", RowBox[{"{", RowBox[{"6.22", ",", ".00126"}], "}"}], ",", RowBox[{"{", RowBox[{"2.3", ",", ".00253"}], "}"}], ",", RowBox[{"{", RowBox[{".61", ",", ".00074"}], "}"}], ",", RowBox[{"{", RowBox[{".05", ",", ".00027"}], "}"}]}], "}"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"\[Tau]Avg", "=", RowBox[{ RowBox[{ RowBox[{"Sum", "[", RowBox[{ RowBox[{ RowBox[{"\[Tau]DelayedList", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], RowBox[{"\[Tau]DelayedList", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}]}], " ", ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "6"}], "}"}]}], "]"}], "/", "\[Beta]value"}], "/", RowBox[{"Log", "[", "2", "]"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"\[Beta]alt", "=", RowBox[{"Sum", "[", RowBox[{ RowBox[{"\[Tau]DelayedList", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], " ", ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "6"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"\[Lambda]Avg", "=", RowBox[{"1", "/", "\[Tau]Avg"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"Print", "[", RowBox[{ "\"\<| \[Beta] (GSR) = \>\"", ",", "\[Beta]alt", ",", "\"\< | \[Lambda] (GSR) = \>\"", ",", "\[Lambda]Avg", ",", "\"\< s^-1 | \[Tau] (GSR) = \>\"", ",", "\[Tau]Avg", ",", "\"\< s |\>\""}], "]"}]}]}]], "Input", CellChangeTimes->{{3.39353833097665*^9, 3.393538518251545*^9}, { 3.393543894609749*^9, 3.3935439139062223`*^9}, {3.393610553578022*^9, 3.393610575824182*^9}, {3.393610678226407*^9, 3.393610689750658*^9}, { 3.393612173199638*^9, 3.3936122760309134`*^9}, {3.3936285884803743`*^9, 3.3936288651831636`*^9}, 3.3936289010122023`*^9, {3.3936289324901237`*^9, 3.3936289332872047`*^9}, {3.393629197171206*^9, 3.3936292214246273`*^9}, { 3.417185222920745*^9, 3.4171852523544483`*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"55.72`", ",", "0.00021`"}], "}"}], ",", RowBox[{"{", RowBox[{"14.1`", ",", "0.00141`"}], "}"}], ",", RowBox[{"{", RowBox[{"6.22`", ",", "0.00126`"}], "}"}], ",", RowBox[{"{", RowBox[{"2.3`", ",", "0.00253`"}], "}"}], ",", RowBox[{"{", RowBox[{"0.61`", ",", "0.00074`"}], "}"}], ",", RowBox[{"{", RowBox[{"0.05`", ",", "0.00027`"}], "}"}]}], "}"}]], "Output", CellChangeTimes->{3.393612251895975*^9, 3.3936122912502937`*^9, 3.393628811704967*^9, 3.393628867661894*^9, 3.393628903154868*^9, 3.393628935007777*^9, 3.41718433768976*^9, 3.417185258136026*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"| \[Beta] (GSR) = \"\>", "\[InvisibleSpace]", "0.00642`", "\[InvisibleSpace]", "\<\" | \[Lambda] (GSR) = \"\>", "\[InvisibleSpace]", "0.09736725573853201`", "\[InvisibleSpace]", "\<\" s^-1 | \[Tau] (GSR) = \"\>", "\[InvisibleSpace]", "10.270393187267997`", "\[InvisibleSpace]", "\<\" s |\"\>"}], SequenceForm[ "| \[Beta] (GSR) = ", 0.00642, " | \[Lambda] (GSR) = ", 0.09736725573853201, " s^-1 | \[Tau] (GSR) = ", 10.270393187267997`, " s |"], Editable->False]], "Print", CellChangeTimes->{3.4171852581438923`*^9}] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["Laplace transform method", "Subsection", CellChangeTimes->{{3.3935413715044317`*^9, 3.3935413773747873`*^9}}], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{"Clear", "[", RowBox[{"\[Lambda]", ",", "\[Beta]", ",", "\[CapitalLambda]"}], "]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"\[Nu]", "=", "2.07"}], ";", RowBox[{"vavg", "=", RowBox[{"2.2", " ", RowBox[{"10", "^", "5"}]}]}], ";", RowBox[{"\[CapitalSigma]fGSTR", "=", "0.083"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"params", "=", RowBox[{"{", RowBox[{ RowBox[{"\[Lambda]", "\[Rule]", "\[Lambda]Avg"}], ",", RowBox[{"\[Beta]", "\[Rule]", "\[Beta]alt"}], ",", RowBox[{"\[CapitalLambda]", "\[Rule]", RowBox[{"1", "/", RowBox[{"(", RowBox[{"vavg", " ", "\[Nu]", " ", "\[CapitalSigma]fGSTR"}], ")"}]}]}]}], "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"sEqn", "[", "s_", "]"}], "=", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"s", "^", "2"}], "+", RowBox[{ RowBox[{"(", RowBox[{"\[Lambda]", "-", RowBox[{ RowBox[{"(", RowBox[{"\[Rho]", "-", "\[Beta]"}], ")"}], "/", "\[CapitalLambda]"}]}], ")"}], "s"}], "-", RowBox[{"\[Rho]", " ", RowBox[{"\[Lambda]", "/", "\[CapitalLambda]"}]}]}], ")"}], "/.", RowBox[{"\[Rho]", "\[Rule]", RowBox[{"\[Delta]", " ", "\[Beta]"}]}]}], ")"}], "/.", "params"}]}]}], "Input", CellChangeTimes->{{3.393091439857071*^9, 3.39309158990304*^9}, { 3.393091626002818*^9, 3.393091674419024*^9}, {3.393091705304225*^9, 3.393091810047282*^9}, {3.393091867708343*^9, 3.393091941740755*^9}, { 3.393091994769335*^9, 3.393092081594578*^9}, {3.393092204859825*^9, 3.393092358425973*^9}, {3.3930927524572687`*^9, 3.393092755413435*^9}, { 3.3930928763237762`*^9, 3.393092921771482*^9}, {3.393092954262705*^9, 3.393093010427066*^9}, {3.3930942584125757`*^9, 3.393094271104239*^9}, { 3.3931723258951893`*^9, 3.393172329712378*^9}, 3.39341577043896*^9, { 3.393415907756563*^9, 3.393416026084511*^9}, {3.393416058735971*^9, 3.393416091062639*^9}, {3.393543389256345*^9, 3.393543505830901*^9}, { 3.3935436010737877`*^9, 3.3935436016689177`*^9}, {3.3936726758770943`*^9, 3.393672676922543*^9}, {3.393672802844356*^9, 3.39367284866545*^9}, { 3.417185281051614*^9, 3.417185304152109*^9}}], Cell[BoxData[ RowBox[{ SuperscriptBox["s", "2"], "+", RowBox[{"s", " ", RowBox[{"(", RowBox[{"0.09736725573853201`", "\[InvisibleSpace]", "-", RowBox[{"37798.200000000004`", " ", RowBox[{"(", RowBox[{ RowBox[{"-", "0.00642`"}], "+", RowBox[{"0.00642`", " ", "\[Delta]"}]}], ")"}]}]}], ")"}]}], "-", RowBox[{"23.627570977596683`", " ", "\[Delta]"}]}]], "Output", CellChangeTimes->{ 3.3930916060173273`*^9, {3.393092130215004*^9, 3.3930921324742413`*^9}, 3.3930921671419573`*^9, {3.3930923101348257`*^9, 3.3930923302104483`*^9}, 3.393092361614689*^9, 3.393092925722363*^9, {3.3930929647733183`*^9, 3.393092968534107*^9}, 3.3930930124947443`*^9, 3.393094273548036*^9, 3.393171555890235*^9, {3.3931723118587227`*^9, 3.393172330603929*^9}, 3.393416030410925*^9, 3.393416093985003*^9, {3.393543499898193*^9, 3.3935435074311867`*^9}, 3.393672678922257*^9, {3.393672835396709*^9, 3.393672850583949*^9}, 3.417184337752564*^9, 3.4171853365404863`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"s1", "[", "\[Delta]_", "]"}], "=", RowBox[{"s", "/.", RowBox[{"Part", "[", RowBox[{ RowBox[{"Solve", "[", RowBox[{ RowBox[{ RowBox[{"sEqn", "[", "s", "]"}], "\[Equal]", "0"}], ",", "s"}], "]"}], ",", "1"}], "]"}]}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"s2", "[", "\[Delta]_", "]"}], "=", RowBox[{"s", "/.", RowBox[{"Part", "[", RowBox[{ RowBox[{"Solve", "[", RowBox[{ RowBox[{ RowBox[{"sEqn", "[", "s", "]"}], "\[Equal]", "0"}], ",", "s"}], "]"}], ",", "2"}], "]"}]}]}]}], "Input", CellChangeTimes->{{3.393092371858533*^9, 3.393092452745001*^9}, { 3.393416121927808*^9, 3.393416123562792*^9}, {3.393543520404553*^9, 3.3935435251765823`*^9}}], Cell[BoxData[ RowBox[{"0.5`", " ", RowBox[{"(", RowBox[{ RowBox[{"-", "242.76181125573856`"}], "+", RowBox[{"242.66444400000003`", " ", "\[Delta]"}], "-", RowBox[{"242.66444400000003`", " ", SqrtBox[ RowBox[{"1.000802645728128`", "\[InvisibleSpace]", "-", RowBox[{"1.9991975152673083`", " ", "\[Delta]"}], "+", RowBox[{"1.`", " ", SuperscriptBox["\[Delta]", "2"]}]}]]}]}], ")"}]}]], "Output", CellChangeTimes->{ 3.3930923890168858`*^9, 3.3930924544791327`*^9, 3.393093015656516*^9, 3.3930942783402777`*^9, 3.393171556106493*^9, 3.393172314952239*^9, { 3.3934161082163887`*^9, 3.3934161246470823`*^9}, {3.3935435109152107`*^9, 3.393543526381467*^9}, 3.3936728563600597`*^9, 3.417184337854116*^9, 3.417185341008685*^9}], Cell[BoxData[ RowBox[{"0.5`", " ", RowBox[{"(", RowBox[{ RowBox[{"-", "242.76181125573856`"}], "+", RowBox[{"242.66444400000003`", " ", "\[Delta]"}], "+", RowBox[{"242.66444400000003`", " ", SqrtBox[ RowBox[{"1.000802645728128`", "\[InvisibleSpace]", "-", RowBox[{"1.9991975152673083`", " ", "\[Delta]"}], "+", RowBox[{"1.`", " ", SuperscriptBox["\[Delta]", "2"]}]}]]}]}], ")"}]}]], "Output", CellChangeTimes->{ 3.3930923890168858`*^9, 3.3930924544791327`*^9, 3.393093015656516*^9, 3.3930942783402777`*^9, 3.393171556106493*^9, 3.393172314952239*^9, { 3.3934161082163887`*^9, 3.3934161246470823`*^9}, {3.3935435109152107`*^9, 3.393543526381467*^9}, 3.3936728563600597`*^9, 3.417184337854116*^9, 3.41718534101639*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"\[Delta]value", "=", RowBox[{"-", ".2"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"t1value", "=", RowBox[{"1", "/", RowBox[{"s1", "[", "\[Delta]value", "]"}]}]}], "\[IndentingNewLine]", RowBox[{"\[Tau]2value", "=", RowBox[{"1", "/", RowBox[{"s2", "[", "\[Delta]value", "]"}]}]}], "\[IndentingNewLine]"}]}]], "Input", CellChangeTimes->{{3.393092463442194*^9, 3.393092647919981*^9}, { 3.393094288960079*^9, 3.3930943125141172`*^9}, {3.393543541565152*^9, 3.393543571027337*^9}, {3.417187429543401*^9, 3.41718746287964*^9}, { 3.417187500869237*^9, 3.417187534076241*^9}}], Cell[BoxData[ RowBox[{"-", "0.003433140722677435`"}]], "Output", CellChangeTimes->{ 3.393543571649393*^9, 3.393672861493699*^9, 3.417184337971614*^9, 3.417185344506762*^9, {3.4171874478651*^9, 3.41718746408388*^9}, 3.417187501735918*^9, 3.4171875365368423`*^9}], Cell[BoxData[ RowBox[{"-", "61.639530567083185`"}]], "Output", CellChangeTimes->{ 3.393543571649393*^9, 3.393672861493699*^9, 3.417184337971614*^9, 3.417185344506762*^9, {3.4171874478651*^9, 3.41718746408388*^9}, 3.417187501735918*^9, 3.417187536545512*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"\[Delta]value", "=", ".2"}], ";"}], "\[IndentingNewLine]", RowBox[{"t1value", "=", RowBox[{"1", "/", RowBox[{"s1", "[", "\[Delta]value", "]"}]}]}], "\[IndentingNewLine]", RowBox[{"t2value", "=", RowBox[{"1", "/", RowBox[{"s2", "[", "\[Delta]value", "]"}]}]}], "\[IndentingNewLine]"}]}]], "Input", CellChangeTimes->{{3.393092463442194*^9, 3.393092647919981*^9}, { 3.393094288960079*^9, 3.3930943125141172`*^9}, {3.393543541565152*^9, 3.393543571027337*^9}, {3.417187429543401*^9, 3.417187446452249*^9}, 3.417187496805028*^9, {3.417187542148982*^9, 3.417187545867958*^9}}], Cell[BoxData[ RowBox[{"-", "0.005147919044050122`"}]], "Output", CellChangeTimes->{ 3.393543571649393*^9, 3.393672861493699*^9, 3.417184337971614*^9, 3.417185344506762*^9, {3.4171874478651*^9, 3.417187467801383*^9}, 3.417187505102186*^9, 3.417187547337769*^9}], Cell[BoxData["41.10732525245341`"], "Output", CellChangeTimes->{ 3.393543571649393*^9, 3.393672861493699*^9, 3.417184337971614*^9, 3.417185344506762*^9, {3.4171874478651*^9, 3.417187467801383*^9}, 3.417187505102186*^9, 3.417187547346258*^9}] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["Direct solution", "Subsection", CellChangeTimes->{{3.3935413828047523`*^9, 3.393541385661656*^9}}], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{"Clear", "[", RowBox[{ "\[Lambda]", ",", "\[Beta]", ",", "\[CapitalLambda]", ",", "\[Rho]"}], "]"}], "\[IndentingNewLine]", RowBox[{"eq1", "=", RowBox[{ RowBox[{ RowBox[{"n", "'"}], "[", "t", "]"}], "\[Equal]", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"\[Rho]", "-", "\[Beta]"}], ")"}], RowBox[{ RowBox[{"n", "[", "t", "]"}], "/", "\[CapitalLambda]"}]}], " ", "+", " ", RowBox[{"\[Lambda]", " ", RowBox[{"c", "[", "t", "]"}]}]}]}]}], "\[IndentingNewLine]", RowBox[{"eq2", "=", RowBox[{ RowBox[{ RowBox[{"c", "'"}], "[", "t", "]"}], "\[Equal]", RowBox[{ RowBox[{"\[Beta]", " ", RowBox[{ RowBox[{"n", "[", "t", "]"}], "/", "\[CapitalLambda]"}]}], "-", RowBox[{"\[Lambda]", " ", RowBox[{"c", "[", "t", "]"}]}]}]}]}]}], "Input", CellChangeTimes->{{3.393541397927334*^9, 3.393541503091761*^9}, { 3.393541786833552*^9, 3.393541815745105*^9}}], Cell[BoxData[ RowBox[{ RowBox[{ SuperscriptBox["n", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Equal]", RowBox[{ RowBox[{"\[Lambda]", " ", RowBox[{"c", "[", "t", "]"}]}], "+", FractionBox[ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"-", "\[Beta]"}], "+", "\[Rho]"}], ")"}], " ", RowBox[{"n", "[", "t", "]"}]}], "\[CapitalLambda]"]}]}]], "Output", CellChangeTimes->{3.393541508409251*^9, 3.393541851486508*^9, 3.3936316792933826`*^9, 3.393672912395813*^9, 3.417184338365992*^9, 3.417185365559346*^9}], Cell[BoxData[ RowBox[{ RowBox[{ SuperscriptBox["c", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Equal]", RowBox[{ RowBox[{ RowBox[{"-", "\[Lambda]"}], " ", RowBox[{"c", "[", "t", "]"}]}], "+", FractionBox[ RowBox[{"\[Beta]", " ", RowBox[{"n", "[", "t", "]"}]}], "\[CapitalLambda]"]}]}]], "Output", CellChangeTimes->{3.393541508409251*^9, 3.393541851486508*^9, 3.3936316792933826`*^9, 3.393672912395813*^9, 3.417184338365992*^9, 3.417185365568185*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{"Clear", "[", RowBox[{"\[Lambda]", ",", "\[Beta]", ",", "\[CapitalLambda]"}], "]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"\[Nu]", "=", "2.07"}], ";", RowBox[{"vavg", "=", RowBox[{"2.2", " ", RowBox[{"10", "^", "5"}]}]}], ";", RowBox[{"\[CapitalSigma]fGSTR", "=", "0.083"}], ";"}], "\[IndentingNewLine]", RowBox[{"params", "=", RowBox[{"{", RowBox[{ RowBox[{"\[Lambda]", "\[Rule]", "\[Lambda]Avg"}], ",", RowBox[{"\[Beta]", "\[Rule]", "\[Beta]alt"}], ",", RowBox[{"\[CapitalLambda]", "\[Rule]", RowBox[{"1", "/", RowBox[{"(", RowBox[{"vavg", " ", "\[Nu]", " ", "\[CapitalSigma]fGSTR"}], ")"}]}]}]}], "}"}]}]}], "Input", CellChangeTimes->{{3.393541523045928*^9, 3.393541611685073*^9}, { 3.393541643723878*^9, 3.393541750542674*^9}, {3.3935418249720707`*^9, 3.393541834559133*^9}, 3.3935418663600283`*^9, {3.3935425999180527`*^9, 3.393542600463711*^9}, {3.3935428098846903`*^9, 3.393542901738736*^9}, { 3.393543090765023*^9, 3.3935430929211206`*^9}, {3.3936316943982077`*^9, 3.393631703532329*^9}, {3.393672920062644*^9, 3.3936729207405148`*^9}, { 3.417185388443492*^9, 3.417185389921506*^9}, 3.417188130445341*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"\[Lambda]", "\[Rule]", "0.09736725573853201`"}], ",", RowBox[{"\[Beta]", "\[Rule]", "0.00642`"}], ",", RowBox[{"\[CapitalLambda]", "\[Rule]", "0.000026456286278182555`"}]}], "}"}]], "Output", CellChangeTimes->{3.417188132016727*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"soln", "=", RowBox[{"Flatten", "[", RowBox[{"DSolve", "[", RowBox[{ RowBox[{"{", RowBox[{"eq1", ",", "eq2", ",", RowBox[{ RowBox[{"n", "[", "0", "]"}], "\[Equal]", "1"}], ",", RowBox[{ RowBox[{"c", "[", "0", "]"}], "\[Equal]", RowBox[{"\[Beta]", " ", "/", RowBox[{"(", RowBox[{"\[Lambda]", " ", "\[CapitalLambda]"}], ")"}]}]}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"n", "[", "t", "]"}], ",", RowBox[{"c", "[", "t", "]"}]}], "}"}], ",", "t"}], "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"nsoln", "[", "t_", "]"}], "=", RowBox[{"Simplify", "[", RowBox[{ RowBox[{"n", "[", "t", "]"}], "/.", "soln"}], "]"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"csoln", "[", "t_", "]"}], "=", RowBox[{ RowBox[{"c", "[", "t", "]"}], "/.", "soln"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"nplot", "[", RowBox[{"\[Delta]_", ",", "t_"}], "]"}], "=", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"nsoln", "[", "t", "]"}], "/.", RowBox[{"\[Rho]", "\[Rule]", RowBox[{"\[Delta]", " ", "\[Beta]"}]}]}], ")"}], "/.", "params"}]}], ";"}]}], "Input", CellChangeTimes->{{3.393541879112661*^9, 3.393541930697509*^9}, { 3.393542063691494*^9, 3.393542104880645*^9}, {3.3935421421236467`*^9, 3.393542260465369*^9}, {3.3935423224069233`*^9, 3.393542373946813*^9}, 3.393631672495837*^9, {3.3936322074496737`*^9, 3.393632224479371*^9}, 3.393633173215163*^9}], Cell[BoxData[ RowBox[{ FractionBox["1", RowBox[{"2", " ", SqrtBox[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"\[Beta]", "+", RowBox[{"\[Lambda]", " ", "\[CapitalLambda]"}], "-", "\[Rho]"}], ")"}], "2"], "+", RowBox[{ "4", " ", "\[Lambda]", " ", "\[CapitalLambda]", " ", "\[Rho]"}]}]]}]], RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"-", FractionBox[ RowBox[{"t", " ", RowBox[{"(", RowBox[{"\[Beta]", "+", RowBox[{"\[Lambda]", " ", "\[CapitalLambda]"}], "-", "\[Rho]", "+", SqrtBox[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"\[Beta]", "+", RowBox[{"\[Lambda]", " ", "\[CapitalLambda]"}], "-", "\[Rho]"}], ")"}], "2"], "+", RowBox[{ "4", " ", "\[Lambda]", " ", "\[CapitalLambda]", " ", "\[Rho]"}]}]]}], ")"}]}], RowBox[{"2", " ", "\[CapitalLambda]"}]]}]], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "\[Beta]"}], "+", RowBox[{ SuperscriptBox["\[ExponentialE]", FractionBox[ RowBox[{"t", " ", SqrtBox[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"\[Beta]", "+", RowBox[{"\[Lambda]", " ", "\[CapitalLambda]"}], "-", "\[Rho]"}], ")"}], "2"], "+", RowBox[{ "4", " ", "\[Lambda]", " ", "\[CapitalLambda]", " ", "\[Rho]"}]}]]}], "\[CapitalLambda]"]], " ", "\[Beta]"}], "-", RowBox[{"\[Lambda]", " ", "\[CapitalLambda]"}], "+", RowBox[{ SuperscriptBox["\[ExponentialE]", FractionBox[ RowBox[{"t", " ", SqrtBox[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"\[Beta]", "+", RowBox[{"\[Lambda]", " ", "\[CapitalLambda]"}], "-", "\[Rho]"}], ")"}], "2"], "+", RowBox[{ "4", " ", "\[Lambda]", " ", "\[CapitalLambda]", " ", "\[Rho]"}]}]]}], "\[CapitalLambda]"]], " ", "\[Lambda]", " ", "\[CapitalLambda]"}], "-", "\[Rho]", "+", RowBox[{ SuperscriptBox["\[ExponentialE]", FractionBox[ RowBox[{"t", " ", SqrtBox[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"\[Beta]", "+", RowBox[{"\[Lambda]", " ", "\[CapitalLambda]"}], "-", "\[Rho]"}], ")"}], "2"], "+", RowBox[{ "4", " ", "\[Lambda]", " ", "\[CapitalLambda]", " ", "\[Rho]"}]}]]}], "\[CapitalLambda]"]], " ", "\[Rho]"}], "+", SqrtBox[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"\[Beta]", "+", RowBox[{"\[Lambda]", " ", "\[CapitalLambda]"}], "-", "\[Rho]"}], ")"}], "2"], "+", RowBox[{ "4", " ", "\[Lambda]", " ", "\[CapitalLambda]", " ", "\[Rho]"}]}]], "+", RowBox[{ SuperscriptBox["\[ExponentialE]", FractionBox[ RowBox[{"t", " ", SqrtBox[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"\[Beta]", "+", RowBox[{"\[Lambda]", " ", "\[CapitalLambda]"}], "-", "\[Rho]"}], ")"}], "2"], "+", RowBox[{ "4", " ", "\[Lambda]", " ", "\[CapitalLambda]", " ", "\[Rho]"}]}]]}], "\[CapitalLambda]"]], " ", SqrtBox[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"\[Beta]", "+", RowBox[{"\[Lambda]", " ", "\[CapitalLambda]"}], "-", "\[Rho]"}], ")"}], "2"], "+", RowBox[{ "4", " ", "\[Lambda]", " ", "\[CapitalLambda]", " ", "\[Rho]"}]}]]}]}], ")"}]}]}]], "Output", CellChangeTimes->{ 3.393542123999763*^9, {3.393542191470846*^9, 3.393542201690546*^9}, 3.393542276528651*^9, {3.393542355632977*^9, 3.3935423754406443`*^9}, 3.393542814438076*^9, 3.393542853060609*^9, {3.393542887629794*^9, 3.39354290838708*^9}, 3.393543104806591*^9, 3.39363179805155*^9, 3.393632232644978*^9, 3.3936331753563137`*^9, 3.393672993270234*^9, 3.417184339676201*^9, 3.417185397300538*^9}] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["Step negative reactivity insertion", "Subsection", CellChangeTimes->{{3.3936327230078278`*^9, 3.393632730182086*^9}}], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"\[Delta]value", "=", RowBox[{"-", ".2"}]}], ";", RowBox[{"\[CapitalDelta]t", "=", "10"}], ";", RowBox[{"tvalue", "=", "40"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Period", "=", RowBox[{"\[CapitalDelta]t", "/", RowBox[{"(", RowBox[{ RowBox[{"Log", "[", RowBox[{"nplot", "[", RowBox[{"\[Delta]value", ",", RowBox[{"tvalue", "+", "\[CapitalDelta]t"}]}], "]"}], "]"}], "-", RowBox[{"Log", "[", RowBox[{"nplot", "[", RowBox[{"\[Delta]value", ",", "tvalue"}], "]"}], "]"}]}], ")"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", RowBox[{"\"\< Insertion of -$\>\"", ",", RowBox[{"-", "\[Delta]value"}], ",", "\"\< of reactivity gives a period of \>\"", ",", "Period", ",", "\"\< sec\>\""}], "]"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"horizaxis", " ", "=", " ", RowBox[{"StyleForm", "[", RowBox[{"\"\\"", ",", RowBox[{"FontFamily", "\[Rule]", "\"\\""}], ",", RowBox[{"FontColor", "\[Rule]", "Blue"}], ",", RowBox[{"FontWeight", "\[Rule]", "Bold"}], ",", RowBox[{"FontSize", "\[Rule]", "12"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"vertaxis", " ", "=", " ", RowBox[{"StyleForm", "[", RowBox[{"\"\\"", ",", RowBox[{"FontFamily", "\[Rule]", "\"\\""}], ",", RowBox[{"FontColor", "\[Rule]", "Blue"}], ",", RowBox[{"FontWeight", "\[Rule]", "Bold"}], ",", RowBox[{"FontSize", "\[Rule]", "12"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"plotname", " ", "=", " ", RowBox[{"StyleForm", "[", RowBox[{"\"\\"", ",", RowBox[{"FontFamily", "\[Rule]", "\"\\""}], ",", RowBox[{"FontColor", "\[Rule]", "Black"}], ",", RowBox[{"FontWeight", "\[Rule]", "Bold"}], ",", RowBox[{"FontSize", "\[Rule]", "14"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"Plot", "[", RowBox[{ RowBox[{"Log", "[", RowBox[{"nplot", "[", RowBox[{"\[Delta]value", ",", "t"}], "]"}], "]"}], ",", RowBox[{"{", RowBox[{"t", ",", "0", ",", "100"}], "}"}], ",", RowBox[{"PlotStyle", "\[Rule]", RowBox[{"{", RowBox[{"{", RowBox[{"Red", ",", RowBox[{"Thickness", "[", ".005", "]"}]}], "}"}], "}"}]}], ",", RowBox[{"Frame", "\[Rule]", "True"}], ",", RowBox[{"GridLines", "\[Rule]", "Automatic"}], ",", RowBox[{"PlotLabel", "\[Rule]", "plotname"}], ",", "\[IndentingNewLine]", RowBox[{"FrameLabel", "\[Rule]", RowBox[{"{", RowBox[{"horizaxis", ",", "vertaxis"}], "}"}]}], ",", RowBox[{"ImageSize", "\[Rule]", "600"}], ",", RowBox[{"Background", "\[Rule]", "LightYellow"}], ",", RowBox[{"PlotRange", "\[Rule]", "All"}]}], "]"}]}], "Input", CellChangeTimes->{{3.3935422852800093`*^9, 3.393542315281653*^9}, { 3.3935423827055283`*^9, 3.393542576600267*^9}, {3.393542929905593*^9, 3.393542987447997*^9}, {3.3935431330409718`*^9, 3.393543322264104*^9}, 3.39363269309023*^9, {3.393632779306052*^9, 3.393632804826092*^9}, 3.39363298466875*^9, 3.417187560982374*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\" Insertion of -$\"\>", "\[InvisibleSpace]", "0.2`", "\[InvisibleSpace]", "\<\" of reactivity gives a period of \"\>", "\[InvisibleSpace]", RowBox[{"-", "61.6395305678390921412`10.718264126891784"}], "\[InvisibleSpace]", "\<\" sec\"\>"}], SequenceForm[ " Insertion of -$", 0.2, " of reactivity gives a period of ", -61.639530568`10.718264126891784, " sec"], Editable->False]], "Print", CellChangeTimes->{ 3.393543325048053*^9, {3.393632666794676*^9, 3.3936326952401543`*^9}, { 3.393632785535231*^9, 3.393632806159863*^9}, 3.393632985299943*^9, 3.393673116383761*^9, 3.417184339904735*^9, 3.417185409197756*^9, 3.417187561708706*^9}], Cell[BoxData[ GraphicsBox[{{}, {}, {RGBColor[1, 0, 0], Thickness[0.005], LineBox[CompressedData[" 1:eJwd03s4VHkYB/CRy9SWddmyaSuDWillpKyEc8iEGCFqkkpuM8h9qqkMSfea LsZKW25dN7UhuunhNauawUpJNyqXRBgMY4xmhtnf8cd5zvN5vr/3/b3nfZ7j Hx9qGkEikaahJ4x4OwsW5l2bOGcGxjI312b3xVg+K9bEcvojEIbOy3ARuWLF rCXX73EEIEympW3SHcIendUN7rLjgzCDH20f/goDlg4zkfoChPwO+o3qXqym QP1YPScFhNetbQw/fcOEPJUzLTMbhOXc2WnjnVgjU+bRdEYEwrdzW4Ns+rDX RZyeWUckIPwWWSX06sHeFMgti2u4IJSVX10V1Y+950nvjeHZIDLyjfq5sAdr Y4qfZDjUgmjL8bEaPTHWXrw69RfHOyBivm2hWvdhHUXs/PDWYRBxzCtz3fux rwUjndPnpYIoB47sS+vFenhDTJ+LF0H0UW64bEiKDTG/J7ZeqIPa4MjlWW/6 MIkNvzY75wTU7i430JD2YJJi2yCMZQW1XI3R+Jl92HBR0sEL5CSozc196unV j0kLJFWraJNQ++Wtx0TZGCbnDaw7UDMX6na5hUfyBzA1s9uHXLkB6iPNLttv GsHUEiNZqNoG6vc1ZoWGDeAkR/OXOylnoP5kCu9M2CBOehawpyKUDfV336W1 Jw/hGs0PniV0u0G99FT4Cb4U15Rxdn0e6IT/Do0s/1DTi0+3m7jyaIICDZcF 1fvNB3D97Voz7HQboPF1yLeKdYP43I23qmnnXOH1eDOJXDiKU3RtdBXWEfBm 8nLEtEIZvsjdaDjCkg7vqLT0vzpk+BL72j3c5vvwwZGts8B0DLdK0ypeargW WthHXR3TZTiVcSM3U8mET1cdNlu5jOEr6cW2PS/i4UtRcpO2QIav2j2xWCfo FLTHe0yrwhW4Hcsm8qAwDDr25vjGFShx+60xfsZkAXRyKk6RKSrcgd2bfcjA Br6ef2kSh+y428fKSRoCXdnm3+tClLgzTz7/eV4hfCuN0R1tV+C4f1Orw9ky 6BYeMSG3K3EX71xfpwo29NS1hesVKHDXbSdHbtQBfO+y0upG9esYUZVr3KOg t6/iYymyW3TMA5u/r0A/ee1PXmge2r6E9ErpBRAv5P3eEKLA1yex9u4p6oYB c9N9bih3P3w4J3P9Mxh0umV+CtV7HDrX2bJRA4Zo8ZrPKUrcc3+BxkMXEkgW aTGlEi6+4XgTXf7rCZCUeCrFJam4V4b837CgBBi2pd5uk6Ti3n9i4avSO2BY oK9dkcDF6YUnWy7xZ8OIN2PLCVTvc1m4f6n+TRh56VLjh7zxjgXNJEkA0m3m n6eh874BTvMLKHkg/XJL/x8q8isWSZFSCaM7+Q5e6LzfFnMV2csURvs5Bk0o 92v+eNqysRxksZOa/ij3306hTr44AzJlytWnKPf/FBLdm8yHsbQrwQuQN201 01uI/l+5xuKZbMLt7Dpty9sg591klaL5A3YIjDVGdWBcX9TWh+YJ6MlLLaXe gvG8Lu0ZKA+MnfApCrWEH0tZlAUlXDywP/D+E89W+HGXF7gMnd+cjC0/cK0J FGvu2Jui/ptHzBz9ysZBURm92gDlW+KOPUh6rAalc6aFjLAshbae6gXKuvui atSPkRwTZb1dCCo6O+MYYWX/tbqJ46BqXSG2RftmKJUy46wAZOzzAcKqWetZ yyjIvi+r0fczVFa92owK5KRSuoTI41ZgZWLkh3tZRD+V5HF5tC+oPjmTcon6 CVlD7vt5oPpMn0MmPKmzUBzbjbxDh070mzRKWKtVhhwnzyLqJ/8w/LDSG/l8 i+lUzmEY8tKQm/MdEwgrOo+v6wLVl2DLRCrqp55p+7alBDl23mO0D4b6t6OL ErnIqTMnifvUjhaCfCPk/MHTU+buVio8kb+W3ST6qXlVGy4Q+eh1MXGfOk/v kkUnqNq0s22JedRV9+0CU5AtOFP7UzdqHhF7INtHk6fcFvAmYw6y57apfaiH bpoZE/VB3llT96nlCSX3kGOcWiTc/wHZ7VqA "]]}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->True, AxesOrigin->{0, 0}, Background->RGBColor[1, 1, 0.85], Frame->True, FrameLabel->{ FormBox[ TagBox[ StyleBox[ "\"t (seconds)\"", FontFamily -> "Tahoma", FontColor -> RGBColor[0, 0, 1], FontWeight -> Bold, FontSize -> 12], StyleForm[#, FontFamily -> "Tahoma", FontColor -> RGBColor[0, 0, 1], FontWeight -> Bold, FontSize -> 12]& ], TraditionalForm], FormBox[ TagBox[ StyleBox[ "\"Log[P/P0]\"", FontFamily -> "Tahoma", FontColor -> RGBColor[0, 0, 1], FontWeight -> Bold, FontSize -> 12], StyleForm[#, FontFamily -> "Tahoma", FontColor -> RGBColor[0, 0, 1], FontWeight -> Bold, FontSize -> 12]& ], TraditionalForm]}, GridLines->Automatic, ImageSize->600, PlotLabel->FormBox[ TagBox[ StyleBox[ "\"Log[P/P0] versus Time \"", FontFamily -> "Tahoma", FontColor -> GrayLevel[0], FontWeight -> Bold, FontSize -> 14], StyleForm[#, FontFamily -> "Tahoma", FontColor -> GrayLevel[0], FontWeight -> Bold, FontSize -> 14]& ], TraditionalForm], PlotRange->{All, All}, PlotRangeClipping->True, PlotRangePadding->{Automatic, Automatic}]], "Output", CellChangeTimes->{ 3.393542413262844*^9, {3.393542554287809*^9, 3.393542577596322*^9}, 3.393542816589963*^9, 3.393542855208921*^9, 3.393542889860279*^9, 3.393542937161759*^9, 3.393542988881764*^9, 3.39354310873623*^9, 3.3935432457099867`*^9, 3.393543325083208*^9, {3.393632667302384*^9, 3.3936326952823877`*^9}, {3.3936327855924683`*^9, 3.3936328062178392`*^9}, 3.393632985356455*^9, 3.393673117114114*^9, 3.4171843402382*^9, 3.4171854093048162`*^9, 3.417187561777461*^9}] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["Step positive reactivity insertion (sub prompt critical)", "Subsection", CellChangeTimes->{{3.3936327230078278`*^9, 3.3936327487743063`*^9}, { 3.3936331839866257`*^9, 3.393633193897419*^9}}], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"\[Delta]value", "=", RowBox[{"+", ".25"}]}], ";", RowBox[{"\[CapitalDelta]t", "=", "10"}], ";", RowBox[{"tvalue", "=", "40"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Period", "=", RowBox[{"\[CapitalDelta]t", "/", RowBox[{"(", RowBox[{ RowBox[{"Log", "[", RowBox[{"nplot", "[", RowBox[{"\[Delta]value", ",", RowBox[{"tvalue", "+", "\[CapitalDelta]t"}]}], "]"}], "]"}], "-", RowBox[{"Log", "[", RowBox[{"nplot", "[", RowBox[{"\[Delta]value", ",", "tvalue"}], "]"}], "]"}]}], ")"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", RowBox[{ "\"\< Insertion of +$\>\"", ",", "\[Delta]value", ",", "\"\< of reactivity gives a period of \>\"", ",", "Period", ",", "\"\< sec\>\""}], "]"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"horizaxis", " ", "=", " ", RowBox[{"StyleForm", "[", RowBox[{"\"\\"", ",", RowBox[{"FontFamily", "\[Rule]", "\"\\""}], ",", RowBox[{"FontColor", "\[Rule]", "Blue"}], ",", RowBox[{"FontWeight", "\[Rule]", "Bold"}], ",", RowBox[{"FontSize", "\[Rule]", "12"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"vertaxis", " ", "=", " ", RowBox[{"StyleForm", "[", RowBox[{"\"\\"", ",", RowBox[{"FontFamily", "\[Rule]", "\"\\""}], ",", RowBox[{"FontColor", "\[Rule]", "Blue"}], ",", RowBox[{"FontWeight", "\[Rule]", "Bold"}], ",", RowBox[{"FontSize", "\[Rule]", "12"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"plotname", " ", "=", " ", RowBox[{"StyleForm", "[", RowBox[{"\"\\"", ",", RowBox[{"FontFamily", "\[Rule]", "\"\\""}], ",", RowBox[{"FontColor", "\[Rule]", "Black"}], ",", RowBox[{"FontWeight", "\[Rule]", "Bold"}], ",", RowBox[{"FontSize", "\[Rule]", "14"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"Plot", "[", RowBox[{ RowBox[{"Log", "[", RowBox[{"nplot", "[", RowBox[{"\[Delta]value", ",", "t"}], "]"}], "]"}], ",", RowBox[{"{", RowBox[{"t", ",", "0", ",", "100"}], "}"}], ",", RowBox[{"PlotStyle", "\[Rule]", RowBox[{"{", RowBox[{"{", RowBox[{"Red", ",", RowBox[{"Thickness", "[", ".005", "]"}]}], "}"}], "}"}]}], ",", RowBox[{"Frame", "\[Rule]", "True"}], ",", RowBox[{"GridLines", "\[Rule]", "Automatic"}], ",", RowBox[{"PlotLabel", "\[Rule]", "plotname"}], ",", "\[IndentingNewLine]", RowBox[{"FrameLabel", "\[Rule]", RowBox[{"{", RowBox[{"horizaxis", ",", "vertaxis"}], "}"}]}], ",", RowBox[{"ImageSize", "\[Rule]", "600"}], ",", RowBox[{"Background", "\[Rule]", "LightYellow"}], ",", RowBox[{"PlotRange", "\[Rule]", "All"}]}], "]"}]}], "Input", CellChangeTimes->{{3.3935422852800093`*^9, 3.393542315281653*^9}, { 3.3935423827055283`*^9, 3.393542576600267*^9}, {3.393542929905593*^9, 3.393542987447997*^9}, {3.3935431330409718`*^9, 3.393543322264104*^9}, 3.39363269309023*^9, {3.39363275165596*^9, 3.3936327523434772`*^9}, 3.393632816425988*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\" Insertion of +$\"\>", "\[InvisibleSpace]", "0.25`", "\[InvisibleSpace]", "\<\" of reactivity gives a period of \"\>", "\[InvisibleSpace]", "30.8331538693915819654`11.124425262328362", "\[InvisibleSpace]", "\<\" sec\"\>"}], SequenceForm[ " Insertion of +$", 0.25, " of reactivity gives a period of ", 30.8331538694`11.124425262328362, " sec"], Editable->False]], "Print", CellChangeTimes->{ 3.393543325048053*^9, {3.393632666794676*^9, 3.3936326952401543`*^9}, 3.3936327533072767`*^9, 3.393632817744987*^9, 3.3936732362934227`*^9, 3.4171843403004923`*^9, 3.4171854186789017`*^9, 3.417187584442231*^9}], Cell[BoxData[ GraphicsBox[{{}, {}, {RGBColor[1, 0, 0], Thickness[0.005], LineBox[CompressedData[" 1:eJwd0HtUzGkYB/AZ3TXRImJRyp5MG0o2Yup5XEpKV0WyVopmQowuFE10kVIj idJZJSdCUlYo2iaJbmdpzYxKWb+KaiPbzVRzqdl3+uN3fudzvu/7Pt/39Toa sOQAjUabRr5A1d++enFu/kSaCRjbOhfXLq2A65wQI6Z2GQj3xxt+10qDEs6y m8WR1SBMMMjM/yKFsgt6v362yQDhZa77hqRgqOJoso9Z1oLw5p/c51H9UJOn LFcaRIPwsXYCFiDU8RX2DpcyQfjKO6v0QC80sSVOwtR6EHb3C642yuFtYWQv I2EQhKNrhboxJ0CUN8YsqeGBSDOhO8bqOLTwR4pHMRNEZosYgVdkQLH7n8av awARx9VvxT4ZdJT8EjObdQ9EkdkhN+ZIobMw/Pr+9iEQJXXHzn4lh095w13a C2JAdJd3W2oxCr38AbZbVhaIvtyX1IzKYID977H29EYQh8xI35U6AYNWGQ2Z V5NAzPO7+RdKYbDE2g84FiDmF5TbD0thqDD0VLpWKIiL7amlfgoYyRsUrHaY BPHAUYsBEwmM8b9tOlljCO+ONdUnPBsHJbvHTavSGZqPX1SW9MhBOThXEqC0 guYkFr01To40lumbvcap0Jzdq043GUfaS++IZwHh0Fxpp+0VIEW6+PFLbs9m aFHr0//eJkU1SeS+f751QUv6BhPbNzLUtpm4VjZhDK33hjbXPFGg/h51HRu9 19DW6ZnckihFQ/fbzx3SNgJFN7U1eyNHYz0rPdnKA9CZaBT6cR4Pl26ZO3SA 6QqfylKMuqt4uGxtQwRP/BC6a07vja3jocVp9RLzWeuhV+4Z0b6Gh5a+t3Iu ydnwZVn9+eayaFzlWmLdW3sU+i3rVwVl8XD14YmfNP3Ow3/7bz0obeGhDccq 6FRdIAzanu1orBDg2l2HPOdrVcPQ+hROaoYA14X3ZZ75wQqGt59uv29ahazD bhZ2I/4w4lt+dd4nAdrzxxa+yr0B38M6oDNQgOglbF93oRQkyW1BcpJv2Jbj YfcsHEb5hi3+5LyNu5OHbzVWwViBo7XpiAA3+QZX2m4JhvGislareAFuPnjo sdWdayBtMOTk1FShwwlubOVIOsg6XCJ3mwrQMZRzPKKwB+SfX7wNIOu3xMVd veT4EhRKpkG9SIBOZ9K62tzpMKkpyo+yrsKtUXn0JxtooNTxeHCSzHc+J3Qd m5eENP2P3iG7KHSJH3sR6MdFWsaTiI1RFG67AvtXx3YifUmibp0Oha43ktuy M+YgvUjQOSKl0O33uihz/QKctiL7Qy2TQvd7Zg5GodU4reKw65a7FHp42y3M M85FNfuOriMk9/ibQ5NFV6JalaDMbiuFnjtNFVouS1Dd5TK9ooFY/D6F2fQI 1d/ocilbCr32GFtO1qaixu6TPg/Jfq8P/gf7wjJQg0rZa068fZfJzMUG0ai5 T+OOczJxR3ijBvMuan47wV5A+nv/Vj2f/l0Ttbgpzmnkft69uTF/WN5GrcnM /NKnFPqETLgVBjBR+7wu6wxZ7/PV5+HTre2ow1giHSbeEQbLT+YLUSfbcvls cp8dwyYsz9JxnL6wI+F9H4U7jyQ+Di1X4vQc2UIn0n+nJNrB0dIFdc11mWzS zzfsUPDKPXWoWyTYbR5MLP+a3zhxDhnrgqfl6Kgsl8y/7E0cNUwn5/sqGI6c n42Jz3cGqXKFRZ+G7zPioucrp3xkBZT2Ew/ynvepPFj+6KAHMtafGP+s2j8h eZ3TsgAZrMQBC/K+vpOai/tDeoizPl5sUHkud716KfGd1xJVv8k1s1pXbSNu LBJMOdJ3Fv80MuxmHPSYsqzr3KbPxFe6w1TzlLrW79oeEBeIW5NV/vHs0mM8 4rKXrKmcZVZ9fS7x+3yNKfMOy2VbkWG/yD9L1U/JFzink9x+hbtC1U+ZOzPb rIsYwH8qFzy08Ykm3reYqeqrbFJL6HciDtW7oOqjpLxF8QbEcYrhqfkDBSbz Vfsz+neq3kepHOM+KCa+2V5xl/of6ygIsw== "]]}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->True, AxesOrigin->{0, 0}, Background->RGBColor[1, 1, 0.85], Frame->True, FrameLabel->{ FormBox[ TagBox[ StyleBox[ "\"t (seconds)\"", FontFamily -> "Tahoma", FontColor -> RGBColor[0, 0, 1], FontWeight -> Bold, FontSize -> 12], StyleForm[#, FontFamily -> "Tahoma", FontColor -> RGBColor[0, 0, 1], FontWeight -> Bold, FontSize -> 12]& ], TraditionalForm], FormBox[ TagBox[ StyleBox[ "\"Log[P/P0]\"", FontFamily -> "Tahoma", FontColor -> RGBColor[0, 0, 1], FontWeight -> Bold, FontSize -> 12], StyleForm[#, FontFamily -> "Tahoma", FontColor -> RGBColor[0, 0, 1], FontWeight -> Bold, FontSize -> 12]& ], TraditionalForm]}, GridLines->Automatic, ImageSize->600, PlotLabel->FormBox[ TagBox[ StyleBox[ "\"Log[P/P0] versus Time \"", FontFamily -> "Tahoma", FontColor -> GrayLevel[0], FontWeight -> Bold, FontSize -> 14], StyleForm[#, FontFamily -> "Tahoma", FontColor -> GrayLevel[0], FontWeight -> Bold, FontSize -> 14]& ], TraditionalForm], PlotRange->{All, All}, PlotRangeClipping->True, PlotRangePadding->{Automatic, Automatic}]], "Output", CellChangeTimes->{ 3.393542413262844*^9, {3.393542554287809*^9, 3.393542577596322*^9}, 3.393542816589963*^9, 3.393542855208921*^9, 3.393542889860279*^9, 3.393542937161759*^9, 3.393542988881764*^9, 3.39354310873623*^9, 3.3935432457099867`*^9, 3.393543325083208*^9, {3.393632667302384*^9, 3.3936326952823877`*^9}, 3.393632753349968*^9, 3.3936328177674637`*^9, 3.3936732363526707`*^9, 3.4171843403580503`*^9, 3.417185418755378*^9, 3.417187584515891*^9}] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["Step positive reactivity insertion (super prompt critical)", \ "Subsection", CellChangeTimes->{{3.3936327230078278`*^9, 3.3936327487743063`*^9}, { 3.393633201328497*^9, 3.39363321053627*^9}}], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"\[Delta]value", "=", RowBox[{"+", "2.25"}]}], ";", RowBox[{"\[CapitalDelta]t", "=", ".10"}], ";", RowBox[{"tvalue", "=", ".40"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Period", "=", RowBox[{"\[CapitalDelta]t", "/", RowBox[{"(", RowBox[{ RowBox[{"Log", "[", RowBox[{"nplot", "[", RowBox[{"\[Delta]value", ",", RowBox[{"tvalue", "+", "\[CapitalDelta]t"}]}], "]"}], "]"}], "-", RowBox[{"Log", "[", RowBox[{"nplot", "[", RowBox[{"\[Delta]value", ",", "tvalue"}], "]"}], "]"}]}], ")"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", RowBox[{ "\"\< Insertion of +$\>\"", ",", "\[Delta]value", ",", "\"\< of reactivity gives a period of \>\"", ",", "Period", ",", "\"\< sec\>\""}], "]"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"horizaxis", " ", "=", " ", RowBox[{"StyleForm", "[", RowBox[{"\"\\"", ",", RowBox[{"FontFamily", "\[Rule]", "\"\\""}], ",", RowBox[{"FontColor", "\[Rule]", "Blue"}], ",", RowBox[{"FontWeight", "\[Rule]", "Bold"}], ",", RowBox[{"FontSize", "\[Rule]", "12"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"vertaxis", " ", "=", " ", RowBox[{"StyleForm", "[", RowBox[{"\"\\"", ",", RowBox[{"FontFamily", "\[Rule]", "\"\\""}], ",", RowBox[{"FontColor", "\[Rule]", "Blue"}], ",", RowBox[{"FontWeight", "\[Rule]", "Bold"}], ",", RowBox[{"FontSize", "\[Rule]", "12"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"plotname", " ", "=", " ", RowBox[{"StyleForm", "[", RowBox[{"\"\\"", ",", RowBox[{"FontFamily", "\[Rule]", "\"\\""}], ",", RowBox[{"FontColor", "\[Rule]", "Black"}], ",", RowBox[{"FontWeight", "\[Rule]", "Bold"}], ",", RowBox[{"FontSize", "\[Rule]", "14"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"Plot", "[", RowBox[{ RowBox[{"Log", "[", RowBox[{"nplot", "[", RowBox[{"\[Delta]value", ",", "t"}], "]"}], "]"}], ",", RowBox[{"{", RowBox[{"t", ",", "0", ",", "1"}], "}"}], ",", RowBox[{"PlotStyle", "\[Rule]", RowBox[{"{", RowBox[{"{", RowBox[{"Red", ",", RowBox[{"Thickness", "[", ".005", "]"}]}], "}"}], "}"}]}], ",", RowBox[{"Frame", "\[Rule]", "True"}], ",", RowBox[{"GridLines", "\[Rule]", "Automatic"}], ",", RowBox[{"PlotLabel", "\[Rule]", "plotname"}], ",", "\[IndentingNewLine]", RowBox[{"FrameLabel", "\[Rule]", RowBox[{"{", RowBox[{"horizaxis", ",", "vertaxis"}], "}"}]}], ",", RowBox[{"ImageSize", "\[Rule]", "600"}], ",", RowBox[{"Background", "\[Rule]", "LightYellow"}], ",", RowBox[{"PlotRange", "\[Rule]", "All"}]}], "]"}]}], "Input", CellChangeTimes->{{3.3935422852800093`*^9, 3.393542315281653*^9}, { 3.3935423827055283`*^9, 3.393542576600267*^9}, {3.393542929905593*^9, 3.393542987447997*^9}, {3.3935431330409718`*^9, 3.393543322264104*^9}, 3.39363269309023*^9, {3.39363275165596*^9, 3.3936327523434772`*^9}, 3.393632816425988*^9, {3.393633214770175*^9, 3.393633218296982*^9}, { 3.393633249908114*^9, 3.393633260089766*^9}, {3.3936332956363993`*^9, 3.393633321241653*^9}, {3.393633545931552*^9, 3.393633548007856*^9}, { 3.393634976719101*^9, 3.393635004078205*^9}, {3.393673300355991*^9, 3.393673307105048*^9}, {3.39367340264501*^9, 3.3936734113301888`*^9}}], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\" Insertion of +$\"\>", "\[InvisibleSpace]", "2.25`", "\[InvisibleSpace]", "\<\" of reactivity gives a period of \"\>", "\[InvisibleSpace]", "0.003295887613385802`", "\[InvisibleSpace]", "\<\" sec\"\>"}], SequenceForm[ " Insertion of +$", 2.25, " of reactivity gives a period of ", 0.003295887613385802, " sec"], Editable->False]], "Print", CellChangeTimes->{3.3936350052255783`*^9, 3.3936733170828323`*^9, 3.39367341248453*^9, 3.4171843404712143`*^9, 3.417185429847183*^9}], Cell[BoxData[ GraphicsBox[{{}, {}, {RGBColor[1, 0, 0], Thickness[0.005], LineBox[CompressedData[" 1:eJwVznk4FHgcBnBXKK2rpOFx5a5pYxIRfX9Fj9KUdUauXDPJkZm1RU/KiFLo 2FyhJWmFEIbqoYhJT3KFolwZ0wyGGTZT65z97R/v8z6fv97X9WyQXqiEhIQU jsf/vd9jqvnDZKjTfp6/RZcWeQIsCGPRGvYywKqJv8IY9AQHQm1+TsFu6Ffa cygkYhbcCfdrBjqPwoS0rfEUxR5CCH++/BlyAuaTcssSFXhAyW99pmHPgu+B wsgamQgI01l4aqdvBiKL9ursjZUQZeRZlMiWh3+bic4XIoMghpBuf23XCIiv vSulyg/AHz6GzVsLMpBE1gOSiokinMt/ZVum5IQkqRRuo8kCxOkI93QI6pG0 9JDlvthOuGTkYqJacRPJEnsmSrNSIIlw46qzWiHa6Bm+SdavBDII141r5RlI Y3uhyRuzXiien6t64O2FTMU59+LdJqFK29GXIVZAu+W0kp/G2wIztXVS4GKI 9s6ZpFkyHOHFQg5Yd51H4E2L/nXCHJokJUdVmOroYNk9/ZPzj6DFtrvDnl2P Dk3mxJ/qPw6sGW/N+2da0BH1mE4ixRDeRlYo2RmJEFkpfQscOwfthprK5YtS yNmRdXxzdgh0VLrVnxJbIheDvKacPAp0b2/rskjqQK73W/vWt1RDj6B0kX9b B7lbLvsNa5Kh1+zj6NhPPvJg9nxz1S6GvpuK7Q3TxeiEax3H1M0J+jkcYYmW IvKaWu83EZACnwwKKXVhaehk8FL/gS1XYODyaU45+oF81pZfi9c6YLB1C5NF ikF+l3KfBN3phS96H1JZlWMoQPuIopGnMQxdMM3rP7OMTvmQRJS4DTDcQFmT O7ADBRJuBynk+8OoejVVr8gRBfJENwK9fGDMd3xbgNZ5FPS4tFp+iglfjw0o eMwQUXBoppoe9wt8XWLJzyltQiGkFAv+xccwnq9fldckg0Lm291keR+BfXT9 5waLcBT6ejOxl18C7B+vCoOOlCNKWqaDso8jTGRxlTX4IkR1CvUuKHQAjkOG tKNbHTq91W77WPw0cLjMguTDl9HpweAlmtksfEt+OOg0oofCMlrkYpwbgUsi ZCRZ2aAzgcY9JnIk4PZJLnZr+6JwneHsnRmywGNEs9PbRShccI8w8SIRJg1c Y/xlclBEVXEun+kPkyyporrRTSgygucdqikNUzRrd0b8OhRlY953S8MFplXZ Vf6lEihq8T15g2AVpltU06ty6ejsu4gEK9/fgB+lIp+y4yWKdiqKnc59DjOK iUS6kwqKLq03Y9N+wkxjXF/PzmFEU+5sq8m8DrP+HAP46ymi0TYOlT98D7Nr WiqWVX6I1k295Uf9BwRlw4wHyQmIvqdNlxcbCUKyYp7WZAWiZ5LNc3/PAuHS 272KKi6IvvJJ17BBDXtYU2rdDWx2wnPRW+y5VdEiC1swTjaLBeGyRusI2xLR V+WK6I8GsaOcK2o1sK0Nmm7ngXBFjXrMnY2db+pN0QXhanBWWk4Moq8FkdJV dmLHxiWk4f21qLu6R0ex031jErjYF77XXr2F/Wybb5gX9l3m5+U5EK5tqCTa 2GG/sTLiMrGr37wflkV0scm+5kZr7LbHTT2AbZHv9mMaeyi1lhWLjVZ45rkg FK9zyXsyje396pe/V7C9RsLju7BT0ck7TdiRzQF0OezswtmOs9iJD90oCLtY kiGH/4ufhO07XoPd2FJykYH9mrzrIB/7nb7NMxL2p136VgbYH690zrOx+arq O/yxxzmBxAxssUhBJxtb4LBAPfQfqlCh4g== "]]}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->True, AxesOrigin->{0, 0}, Background->RGBColor[1, 1, 0.85], Frame->True, FrameLabel->{ FormBox[ TagBox[ StyleBox[ "\"t (seconds)\"", FontFamily -> "Tahoma", FontColor -> RGBColor[0, 0, 1], FontWeight -> Bold, FontSize -> 12], StyleForm[#, FontFamily -> "Tahoma", FontColor -> RGBColor[0, 0, 1], FontWeight -> Bold, FontSize -> 12]& ], TraditionalForm], FormBox[ TagBox[ StyleBox[ "\"Ln[P/P0]\"", FontFamily -> "Tahoma", FontColor -> RGBColor[0, 0, 1], FontWeight -> Bold, FontSize -> 12], StyleForm[#, FontFamily -> "Tahoma", FontColor -> RGBColor[0, 0, 1], FontWeight -> Bold, FontSize -> 12]& ], TraditionalForm]}, GridLines->Automatic, ImageSize->600, PlotLabel->FormBox[ TagBox[ StyleBox[ "\"Ln[P/P0] versus Time \"", FontFamily -> "Tahoma", FontColor -> GrayLevel[0], FontWeight -> Bold, FontSize -> 14], StyleForm[#, FontFamily -> "Tahoma", FontColor -> GrayLevel[0], FontWeight -> Bold, FontSize -> 14]& ], TraditionalForm], PlotRange->{All, All}, PlotRangeClipping->True, PlotRangePadding->{Automatic, Automatic}]], "Output", CellChangeTimes->{3.393635005235808*^9, 3.393673317126027*^9, 3.393673412486953*^9, 3.417184340513402*^9, 3.417185429905602*^9}] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["Subcritical multiplication", "Subsection", CellChangeTimes->{{3.3936356601946507`*^9, 3.393635664973283*^9}}], Cell[BoxData[""], "Input", CellChangeTimes->{{3.393635814623437*^9, 3.393635818729594*^9}}], Cell["\<\ Suppose k < 1; therefore each neutron in the parent generation produces k \ offspring neutrons. If a source is present that delivers S neutrons/sec, \ then in each neutron lifetime Ns = S \[Tau] neutrons are delivered by the \ source. If the source is introduced at time t=0, then the number of neutrons \ after one lifetime is Ns, after two lifetimes, Ns + k Ns, after three \ lifetimes, Ns + k (Ns + k Ns). In other words after one additional lifetime, \ the number of neutrons is increased by the number delivered by the source \ plus k x the number already present. Thus, after many generations: \ N=Ns(1+k+k^2+k^3+...)=Ns/(1-k). Thus, introducing a neutron source into a \ subcritical (k<1) reactor causes the number of neutrons to be multiplied by \ the factor: 1/(1-k). This can be re-expressed in terms of the reactivity: \ \[Rho] = (k-1)/k or k = 1/(1-\[Rho]), giving the subcritical multiplication \ factor as (1-\[Rho])/\[Rho].\ \>", "Text", CellChangeTimes->{{3.393635678641704*^9, 3.3936362824059153`*^9}, { 3.39363673800054*^9, 3.3936369502719088`*^9}, {3.393636996690373*^9, 3.393637012954454*^9}, {3.393637049241707*^9, 3.3936371448697157`*^9}, { 3.417185453765007*^9, 3.4171854575661497`*^9}}, FontSize->14] }, Open ]] }, Open ]] }, Open ]] }, WindowSize->{1024, 905}, WindowMargins->{{182, Automatic}, {Automatic, 27}}, FrontEndVersion->"6.0 for Mac OS X PowerPC (32-bit) (March 13, 2008)", StyleDefinitions->FrontEnd`FileName[{"Creative"}, "PastelColor.nb", CharacterEncoding -> "UTF-8"] ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[CellGroupData[{ Cell[590, 23, 99, 1, 66, "Title"], Cell[692, 26, 138, 5, 56, "Text"], Cell[CellGroupData[{ Cell[855, 35, 206, 3, 72, "Section"], Cell[CellGroupData[{ Cell[1086, 42, 2859, 66, 158, "Input"], Cell[CellGroupData[{ Cell[3970, 112, 459, 9, 20, "Print"], Cell[4432, 123, 489, 11, 20, "Print"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[4970, 140, 841, 20, 53, "Input"], Cell[5814, 162, 418, 7, 20, "Print"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[6281, 175, 257, 3, 72, "Section"], Cell[CellGroupData[{ Cell[6563, 182, 131, 1, 34, "Subsection"], Cell[CellGroupData[{ Cell[6719, 187, 3678, 93, 128, "Input"], Cell[10400, 282, 897, 18, 20, "Print"] }, Open ]], Cell[CellGroupData[{ Cell[11334, 305, 2638, 67, 128, "Input"], Cell[13975, 374, 676, 17, 37, "Output"], Cell[14654, 393, 602, 13, 20, "Print"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[15305, 412, 116, 1, 34, "Subsection"], Cell[CellGroupData[{ Cell[15446, 417, 2324, 54, 83, "Input"], Cell[17773, 473, 1028, 20, 40, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[18838, 498, 783, 23, 53, "Input"], Cell[19624, 523, 794, 17, 55, "Output"], Cell[20421, 542, 793, 17, 55, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[21251, 564, 699, 16, 98, "Input"], Cell[21953, 582, 273, 5, 37, "Output"], Cell[22229, 589, 269, 5, 37, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[22535, 599, 699, 15, 98, "Input"], Cell[23237, 616, 272, 5, 37, "Output"], Cell[23512, 623, 252, 4, 37, "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[23813, 633, 105, 1, 34, "Subsection"], Cell[CellGroupData[{ Cell[23943, 638, 972, 29, 68, "Input"], Cell[24918, 669, 571, 16, 54, "Output"], Cell[25492, 687, 510, 14, 54, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[26039, 706, 1242, 27, 68, "Input"], Cell[27284, 735, 296, 7, 37, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[27617, 747, 1653, 48, 83, "Input"], Cell[29273, 797, 4277, 119, 137, "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[33599, 922, 124, 1, 34, "Subsection"], Cell[CellGroupData[{ Cell[33748, 927, 3260, 77, 173, "Input"], Cell[37011, 1006, 724, 16, 20, "Print"], Cell[37738, 1024, 3914, 77, 419, "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[41701, 1107, 199, 2, 34, "Subsection"], Cell[CellGroupData[{ Cell[41925, 1113, 3225, 77, 173, "Input"], Cell[45153, 1192, 683, 13, 20, "Print"], Cell[45839, 1207, 3892, 77, 423, "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[49780, 1290, 200, 3, 34, "Subsection"], Cell[CellGroupData[{ Cell[50005, 1297, 3565, 81, 173, "Input"], Cell[53573, 1380, 548, 11, 20, "Print"], Cell[54124, 1393, 3301, 67, 420, "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[57474, 1466, 116, 1, 34, "Subsection"], Cell[57593, 1469, 92, 1, 37, "Input"], Cell[57688, 1472, 1250, 19, 131, "Text"] }, Open ]] }, Open ]] }, Open ]] } ] *) (* End of internal cache information *)