Xenon Concentration

The Xenon buildup in fission fuel is an important parameter to consider during burnup of fuel. Xenon is a radioactive gas and causes many problems within the fuel cycle. Besides these issues, calculation of Xenon buildup is also an illustration of how to use the Bateman equations, which has many other uses in nuclear engineering.

To start, we know that Xenon can be produced directly in the fission process, at the rate of \(\gamma_{Xe}\), and can also come from the decay of Iodine that was produced by fission (at the rate of \(\gamma_{I}\)). Thus, we must define two equations

\[\frac{dN_{I}}{dt}=\underbrace{\gamma_{I}F\left(t\right)}_{\text{production}} -\underbrace{\overline{\lambda}_{I}N_{I}\left(t\right)}_{\text{loss}}\]
\[ \begin{align}\begin{aligned}\frac{dN_{Xe}}{dt}=\underbrace{\gamma_{Xe}F\left(t\right)+\lambda_{I}N_{I}\left(t\right)}_{\text{production}} -\underbrace{\overline{\lambda}_{Xe}N_{Xe}\left(t\right)}_{\text{loss}}\\with :math:`\sigma_{.}` the cross section for neutron absorption,\end{aligned}\end{align} \]

\(\phi\) the neutron flux, \(F\left( t \right)\) the fission rate, and

\[\overline{\lambda}_{I}=\lambda_{I}+\sigma_{I}\phi\left(t\right)\]
\[ \begin{align}\begin{aligned}\overline{\lambda}_{Xe}=\lambda_{Xe}+\sigma_{Xe}\phi\left(t\right)\\with :math:`\lambda_{x}` the decay constant.\end{aligned}\end{align} \]

Conceptually, both of these equations are simple: they both read that the change in the concentration of the atom is the production of that atom (whether directly from fission or from decay of Iodine in the case of xenon) less the loss of that atom. The loss of that atom is the sum of the natural decay and the neutron absorption probability (which will transmute the atom).

To solve these equations, the integrating factor is used, and these steps are shown below

\[\frac{d}{dt}\left[N_{I}\left(t\right)e^{\overline{\lambda}_{I}t}\right]=\gamma_{I}Fe^{\overline{\lambda}_{I}t}\]
\[N_{I}\left(t\right)=N_{I}^{0}e^{\overline{\lambda}_{I}t}+\frac{\gamma_{I}F}{\lambda_{I}}\left(1-e^{\overline{\lambda}_{I}t}\right)\]
\[\frac{d}{dt}\left[N_{Xe}e^{\overline{\lambda}_{Xe}t}\right]=\left[\gamma_{Xe}F+\lambda_{I}N_{I}\left(t\right)\right]e^{\overline{\lambda}_{Xe}t}\]
\[\begin{split}\begin{multline*} N_{Xe}\left(t\right)=N_{Xe}^{0}e^{\overline{\lambda}_{Xe}t}+\frac{\gamma_{Xe}F}{\lambda_{Xe}}\left(1-e^{-\overline{\lambda}_{Xe}t}\right)+\frac{\lambda_{I}N_{I}^{0}}{\overline{\lambda}_{Xe}-\overline{\lambda}_{I}}\left(e^{-\overline{\lambda}_{I}t}-e^{-\overline{\lambda}_{Xe}t}\right)\\ +\frac{\lambda_{I}\gamma_{I}F}{\overline{\lambda}_{I}}\left[\frac{1}{\overline{\lambda}_{Xe}}\left(1-e^{-\overline{\lambda}_{Xe}t}\right)-\frac{1}{\overline{\lambda}_{Xe}-\overline{\lambda}_{I}}\left(e^{-\overline{\lambda}_{I}t}-e^{-\overline{\lambda}_{Xe}t}\right)\right] \end{multline*}\end{split}\]

It should be noted that this can be more easily solved in matrix form, and modern computing makes it almost trivial to solve. The program ORIGEN is often used for this purpose.