K.A. Bala, M.R. Omar, John Y.H. Soo, W.M.H. Wan Mokhtar
{"title":"CNUCTRAN: A program for computing final nuclide concentrations using a direct simulation approach","authors":"K.A. Bala, M.R. Omar, John Y.H. Soo, W.M.H. Wan Mokhtar","doi":"10.1016/j.cpc.2024.109258","DOIUrl":null,"url":null,"abstract":"<div><p>It is essential to precisely determine the evolving concentrations of radioactive nuclides within transmutation problems. It is also a crucial aspect of nuclear physics with widespread applications in nuclear waste management and energy production. This paper introduces <span>CNUCTRAN</span>, a novel computer program that employs a probabilistic approach to estimate nuclide concentrations in transmutation problems. <span>CNUCTRAN</span> directly simulates nuclei transformations arising from various nuclear reactions, diverging from the traditional deterministic methods that solve the Bateman equation using matrix exponential approximation. This approach effectively addresses numerical challenges associated with solving the Bateman equations, therefore, circumventing the need for matrix exponential approximations that risk producing nonphysical concentrations. Our sample calculations using <span>CNUCTRAN</span> shows that the concentration predictions of <span>CNUCTRAN</span> have a relative error of less than 0.001% compared to the state-of-the-art method, CRAM, in different test cases. This makes <span>CNUCTRAN</span> a valuable alternative tool for transmutation analysis.</p></div><div><h3>Program summary</h3><p><em>Program Title:</em> <span>CNUCTRAN</span></p><p><em>CPC Library link to program files:</em> <span>https://doi.org/10.17632/b484w2vx52.1</span><svg><path></path></svg></p><p><em>Developer's repository link:</em> <span>https://github.com/rabieomar92/cnuctran/releases</span><svg><path></path></svg></p><p><em>Licensing provisions:</em> MIT</p><p><em>Programming language:</em> C++</p><p><em>Nature of problem:</em> <span>CNUCTRAN</span> simulates the transmutation of various nuclides such as decays, fissions, and neutron induced reactions using a direct simulation approach. It has the capability of predicting the final concentration of a large system of nuclides altogether after a specified time step, <span><math><msub><mrow><mi>t</mi></mrow><mrow><mi>f</mi></mrow></msub></math></span>.</p><p><em>Solution method:</em> <span>CNUCTRAN</span> works based on the novel probabilistic method such that it does not compute the final nuclide concentrations by solving Bateman equations. Instead, it statistically tracks nuclide transformations into one another in a transmutation problem. The technique encapsulates various possible nuclide transformations into a sparse transfer matrix, <span><math><mi>T</mi></math></span>, whose elements are made up of various nuclear reaction probabilities. Next, <span><math><mi>T</mi></math></span> serves as a matrix operator acting on the initial nuclide concentrations, <span><math><mi>y</mi><mo>(</mo><mn>0</mn><mo>)</mo></math></span>, producing the final nuclide concentrations, <strong>y</strong>.</p></div>","PeriodicalId":285,"journal":{"name":"Computer Physics Communications","volume":null,"pages":null},"PeriodicalIF":7.2000,"publicationDate":"2024-05-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Computer Physics Communications","FirstCategoryId":"101","ListUrlMain":"https://www.sciencedirect.com/science/article/pii/S0010465524001814","RegionNum":2,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q1","JCRName":"COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS","Score":null,"Total":0}
引用次数: 0
Abstract
It is essential to precisely determine the evolving concentrations of radioactive nuclides within transmutation problems. It is also a crucial aspect of nuclear physics with widespread applications in nuclear waste management and energy production. This paper introduces CNUCTRAN, a novel computer program that employs a probabilistic approach to estimate nuclide concentrations in transmutation problems. CNUCTRAN directly simulates nuclei transformations arising from various nuclear reactions, diverging from the traditional deterministic methods that solve the Bateman equation using matrix exponential approximation. This approach effectively addresses numerical challenges associated with solving the Bateman equations, therefore, circumventing the need for matrix exponential approximations that risk producing nonphysical concentrations. Our sample calculations using CNUCTRAN shows that the concentration predictions of CNUCTRAN have a relative error of less than 0.001% compared to the state-of-the-art method, CRAM, in different test cases. This makes CNUCTRAN a valuable alternative tool for transmutation analysis.
Program summary
Program Title:CNUCTRAN
CPC Library link to program files:https://doi.org/10.17632/b484w2vx52.1
Nature of problem:CNUCTRAN simulates the transmutation of various nuclides such as decays, fissions, and neutron induced reactions using a direct simulation approach. It has the capability of predicting the final concentration of a large system of nuclides altogether after a specified time step, .
Solution method:CNUCTRAN works based on the novel probabilistic method such that it does not compute the final nuclide concentrations by solving Bateman equations. Instead, it statistically tracks nuclide transformations into one another in a transmutation problem. The technique encapsulates various possible nuclide transformations into a sparse transfer matrix, , whose elements are made up of various nuclear reaction probabilities. Next, serves as a matrix operator acting on the initial nuclide concentrations, , producing the final nuclide concentrations, y.
期刊介绍:
The focus of CPC is on contemporary computational methods and techniques and their implementation, the effectiveness of which will normally be evidenced by the author(s) within the context of a substantive problem in physics. Within this setting CPC publishes two types of paper.
Computer Programs in Physics (CPiP)
These papers describe significant computer programs to be archived in the CPC Program Library which is held in the Mendeley Data repository. The submitted software must be covered by an approved open source licence. Papers and associated computer programs that address a problem of contemporary interest in physics that cannot be solved by current software are particularly encouraged.
Computational Physics Papers (CP)
These are research papers in, but are not limited to, the following themes across computational physics and related disciplines.
mathematical and numerical methods and algorithms;
computational models including those associated with the design, control and analysis of experiments; and
algebraic computation.
Each will normally include software implementation and performance details. The software implementation should, ideally, be available via GitHub, Zenodo or an institutional repository.In addition, research papers on the impact of advanced computer architecture and special purpose computers on computing in the physical sciences and software topics related to, and of importance in, the physical sciences may be considered.