International Journal of Systems Science and Applied Mathematics
Volume 5, Issue 3, September 2020, Pages: 27-31
Received: Dec. 18, 2019;
Accepted: Dec. 31, 2019;
Published: Sep. 14, 2020
Views 112 Downloads 34
Michael Erle Lord, Intera Inc., Boulder, CO, USA
Many engineering applications involving radioactive materials requires the time history of the radioactivity of nuclides within a decay chain. The system of differential equations with initial conditions or initial value problem describing radioactive decay of a parent and daughter nuclides was posited by Ernest Rutherford who was awarded a Noble Prize in 1910 for this work. Harry Bateman (1910) provided an analytic solution to the radioactive decay chain problem with the constraint that initial inventory of all daughter elements is zero. Required data for the decay chain calculation consists of the parent and all daughters’ radioactive half-life. The half-life for essentially all radionuclides has been established and is available from multiple sources. Solutions other than Bateman’s (non-zero initial conditions) can be computed analytically but become unwieldy for longer decay chains. For this reason, many applications use a numerical solution. However, a numerical solution can require constraints on the time step size. The proposed method of false rates provides a unique algorithm for the decay chain activities. The method treats the decay chain with arbitrary initial conditions and the calculation is analytic or exact. The method is unexpectedly simplistic. An example decay chain calculation compares the solutions by the method of false rates with a numerical method. The comparison is a verification of the method of false rates calculation. The method of false rates is easily coded as a stand-alone application or as a sub-module of a more general code such as a contaminant transport model.
Michael Erle Lord,
Method of False Rates for Radioactive Decay Chain Calculation, International Journal of Systems Science and Applied Mathematics.
Vol. 5, No. 3,
2020, pp. 27-31.
Copyright © 2020 Authors retain the copyright of this article.
This article is an open access article distributed under the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0/
) which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Goldsim Technology Group, GoldSim Simulation Software.
Argonne National Laboratory, Managing Radioactive Material Inventories, LMS-PROC_45, Rev. 11, 08/04/2014.
Oak Ridge Nation Laboratory, Radiological Toolbox, v3.0.0/5/1/2004. https://crpk.oml.gov/software.
Radiological Toolbox User’s Guide, NUREG/CR-7166 2013.
Bateman, H., Solution of a system of differential equations occurring in the theory of radioactive transformations, Proc. Cambridge Philos. Soc., 15, 1910.
Moral, L., Pacheco, A. F., Algebraic approach to the radioactive decay equations, American Journal of Physics, 71, 2003.
Bakin, R. I., Kiselev, A. A., Shvedov, A. M. et al., Computational Errors in the Calculation of Long Radioactive Decay Chains. At Energy 123, 406–411, 2018, doi: 10.1007/s10512-018-0360-2.
Massey, F., Prentis, J., Power law approximations for radioactive decay chains, Applied Mathematics and Computation, 245, October 2014.
Krane, K., Introductory Nuclear Physics, John Wiley and Sons Inc., 1987.
International Commission of Radiological Protection, 2008. Nuclear Decay Data for Dosimetric Calculations. ICRP Publication 107. Ann. ICRP 38 (3).
Burden, R., Faires, D., Burden, A., Numerical Analysis, 10 ed. Cengage Learning, 2016.