On a point kinetic model for nuclear reactors considering the variation in fuel composition

• Published in: Progress in nuclear energy (New series) Vol. 118; p. 103134
• Main Authors: Paganin, T.M., Bodmann, B.E.J., Vilhena, M.T.
• Format: Journal Article
• Language: English
• Published: Oxford Elsevier Ltd 01.01.2020; Elsevier BV
• Subjects: Adomian decomposition; Algorithms; Chemical composition; Computer simulation; Decomposition; Decomposition reactions; Fuels; Infinite series; Kinetics; Linear equations; Matrix methods; Neutron poisons; Nonlinear analysis; Nonlinear equations; Nuclear capture; Nuclear fuels; Nuclear reactors; Organic chemistry; Point kinetics; Poisons; Polynomials; Reaction kinetics; Reactivity; Reactivity decomposition; Stiffness; Time; Transuranium; Reactivity decomposition; Neutron poisons; Point kinetics; Adomian decomposition; Transuranium
• ISSN: 0149-1970; 1878-4224
• DOI: 10.1016/j.pnucene.2019.103134

The present discussion is an extension of a new model recently developed for neutron point kinetics, where the reactivity is decomposed in terms of short and long timescales. The short timescale corresponds to the operational control of the reactor, while the long timescale corresponds to the effects of the chemical composition variation in the nuclear fuel as a consequence of fission and neutron capture processes. This work considers not only the effects of the main neutron poisons but contains also a first analysis of the effects of some transuranium. The proposed model consists of a system of non-linear equations, which is solved by the use of the Adomian decomposition method. This method expands the solution into an infinite series, whose contributions are obtained by the solution of a recursive equation system where the non-linear terms are treated as source terms. These are calculated by so-called Adomian polynomials and are constructed from the solutions obtained in previous recursion steps. An approach is proposed, that diagonalizes the linear matrix in order to circumvent the existing stiffness problem, which is also responsible for instabilities in the computational algorithm. Two case studies are presented, analysing the non-linear behaviour of nuclear reactor kinetics with constant reactivity using input parameters for which the algorithm was stable. Further, a standard point kinetics case with periodic reactivity is performed in order to compare the present model to a standard point kinetics simulation from the literature. •We present an extension of a model recently developed for neutron point kinetics.•Reactivity is decomposed in terms of short and long time scales.•The model takes into account effects of variations in the nuclear fuel composition.•The system of nonlinear equations is solved by the Adomian decomposition method.•An analysis of the solution consistency and convergence is proposed.