Solution of the Equations of Nuclide Kinetics

• Published in: Atomic energy (New York, N.Y.) N.Y.), 2021-09, Vol.130 (5), p.262-266
• Main Authors: Seleznev, E. F., Drobyshev, Yu. Yu, Karpov, S. A., Bukhtiyarov, I. R.
• Format: Article
• Language: English
• Subjects: Approximation; Approximation method; Chebyshev approximation; Computer applications; Energy industry; Exact solutions; Hadrons; Heavy Ions; Iterative methods; Kinetics; Mathematical analysis; Nuclear Chemistry; Nuclear Energy; Nuclear Physics; Nuclear power plants; Physics; Physics and Astronomy; Uncertainty analysis
• ISSN: 1063-4258; 1573-8205
• DOI: 10.1007/s10512-021-00806-8

Many methods and computational programs are now available for solving burnup equations. Three fundamentally different approaches are combined in a single control program OpenBPS: solution of the nuclide kinetics equation by expansion of a matrix exponential, an iterative method with the possibility of analyzing uncertainties, and a direct analytical solution. The matrix exponential solution is obtained using the modern Chebyshev rational approximation method CRAM. An iterative method taking into account the errors of the initial data and data on decay and cross sections, made possible by an exponential representation of the increment of the nuclide concentration at the BPSE computational step, affords the possibility of analyzing the uncertainty of the nuclear concentration. An analytical computational method based on the use of modified Bateman functions, which is the basis of the ASBE computational code, was developed for the accelerated solution of the burnup equations. Using this program, a solution is obtained more than 100 times faster than the solution speed of programs based on the iterative method with high solution accuracy.