Repeat problem [9.10] for a plutonium fueled sodium cooled fast reactor with the parameters: Λ = 0.5 × 10^{-6} s , α = -1.8 × 10^{-5} /°C M_ƒc_ƒ = 5.0 × 10^6 J/°C , τ = 4.0 s
β1 := 0.00021 \operatorname{c1}:=\mathbb{\beta 1}\cdot {\frac{\operatorname{P0}}{(\lambda1\cdot\Lambda)}} \mathrm{c}2:=\beta 2.{\frac{\mathrm{P}0}{(\lambda2\cdot\Lambda)}}
β2 := 0.00142
β3 := 0.00128 \operatorname{c3}:=\mathbb{\beta 3}\cdot {\frac{\operatorname{P0}}{(\lambda3\cdot\Lambda)}} \mathrm{c}4:=\beta 4.{\frac{\mathrm{P}0}{(\lambda4\cdot\Lambda)}}
β4 := 0.00257
β5 := 0.00075 \operatorname{c5}:=\mathbb{\beta 5}\cdot {\frac{\operatorname{P0}}{(\lambda5\cdot\Lambda)}} \mathrm{c}6:=\beta 6.{\frac{\mathrm{P}0}{(\lambda6\cdot\Lambda)}}
β6 := 0.00027
Sample calculations shown for $0.10 below using the stiff differential equations solver Radau (Mathcad) to obtain the vector Y.
\mathrm{Rf}:={\frac{\tau}{\mathrm{MfCf}}}=0.8 dollars := 0.1
Rf.WCp = 48 α1 = -1.8 × 10^{-5}
\mathrm y := \begin{pmatrix}P0 \\ c1 \\ c2 \\ c3 \\ c4 \\ c5 \\c6 \\ \mathrm {Tf0} \end{pmatrix}
n := 0 .. 1000