- Introduce
- Editorial Board
- Journal
- Special Source
- Excellent Articles
- Authors
- Reviewers
- Publication Statement
- 个人中心
- 退出登录
浏览全部资源
扫码关注微信
- Introduce
- Editorial Board
- Journal
- Special Source
- Excellent Articles
- Authors
- Reviewers
- Publication Statement
Comparison and analysis of several typical algorithms in solving depletion equations
- NUCLEAR TECHNIQUES Vol. 41, Issue 3, Article number: 030501(2018)
- 作者机构：
水力机械过渡过程教育部重点实验室（武汉大学） 武汉 430072
- 作者简介：
- 基金信息：National Natural Science Foundation of China
DOI：10.11889/j.0253-3219.2018.hjs.41.030501
CLC： TL329.2Received：21 August 2017，
Revised：28 September 2017，
Published：10 March 2018
- 稿件说明：
移动端阅览
Jie TAN, Peng ZHANG. Comparison and analysis of several typical algorithms in solving depletion equations[J]. Nuclear techniques, 2018, 41(3): 030501
Jie TAN, Peng ZHANG. Comparison and analysis of several typical algorithms in solving depletion equations[J]. Nuclear techniques, 2018, 41(3): 030501 DOI： 10.11889/j.0253-3219.2018.hjs.41.030501.
摘要
为求解点堆燃耗方程组，通过自主编程燃耗计算模块，分别采用Chebyshev有理近似方法（Chebyshev Rational Approximation Method，CRAM）、Padé近似方法、Krylov Subspace方法以及线性子链法（Transmutation Trajectory Analysis，TTA）对不同规模的点堆燃耗系统进行燃耗计算，并与高精度下16阶系数的CRAM的结果进行对比分析。结果表明，CRAM能够高效地获得精确的计算结果，而且其系数阶数越高，精度越高，步长稳定性越高，因此目前燃耗计算多采用这种方法；Padé近似方法在双精度下对于中小规模燃耗系统能够快速获得精确结果，受步长影响小，而对于大规模燃耗系统需要采用高精度计算，且受步长影响也小；Krylov Subspace方法在双精度下可准确求解中小规模燃耗矩阵，但不适用于大规模燃耗求解；TTA适用于求解各规模燃耗方程，但由于线性链的简化，相对参考值精度较低，求解速度相对较慢。本文所获结果对于不同燃耗计算系统的开发及应用具有较好的借鉴意义。
Abstract
Background
2
It is very important to trace fuel concentration changing with high precision and efficiency in the current core simulation calculations
where an essential part is to solve depletion equations precisely.
Purpose
2
This study aims to solve point-core depletion equations in different sizes with various methods and find out the characteristics of these methods by a self-programming burnup module.
Methods
2
Four methods including Chebyshev rational approximation method (CRAM)
Padé approximation method
Krylov Subspace method and transmutation trajectory analysis (TTA) were used to realize burnup calculations and comparative analysis with the results by CRAM with 16 orders in high precision for different-scale burnup systems.
Results
2
The results suggested that CRAM could provide very accurate solutions in a very short computing time
and the higher the coefficient order
the higher the precision and step stability. Padé approximation method could get accurate results affected marginally by step length with double precision in small and medium systems quickly
but high precision computation should be adopted for large systems while little affected by step length. Krylov Subspace method had accurate solutions in small and medium systems with double floats. However
it can't be applied to large cases. TTA was applicable to burnup equations at any scale with a weak accuracy and a slower speed because of linear chain simplification.
Conclusion
2
The four methods can be adapted to the depletion calculations of the cores in different cases
and the responding results can be treated as reference for the development and application of different burnup calculation systems.
关键词
Keywords
references
J Cetenar . General solution of bateman equations for nuclear transmutations . Annuals of Nuclear Energy , 2006 . 33 640 - 645 . DOI: 10.1016/j.anucene.2006.02.004 http://doi.org/10.1016/j.anucene.2006.02.004 .
M Pusa , J Leppänen . Computing the matrix exponential in burnup calculations . Nuclear Science and Engineering , 2010 . 164 ( Pt23 ): 140 - 150 . DOI: 10.13182/NSE09-14 http://doi.org/10.13182/NSE09-14 .
A E Isotalo , P A Aarnio . Comparison of depletion algorithms for large systems of nuclides . Annals of Nuclear Energy , 2011 . 38 ( 2 ): 261 - 268 . DOI: 10.1016/j.anucene.2010.10.019 http://doi.org/10.1016/j.anucene.2010.10.019 .
范 文玎 , 孙 光耀 , 张 彬航 , . 基于切比雪夫有理逼近方法的蒙特卡罗燃耗计算研究与验证 . 核技术 , 2016 . 39 ( 4 ): 040501 DOI: 10.11889/j.0253-3219.2016.hjs.39.040501 http://doi.org/10.11889/j.0253-3219.2016.hjs.39.040501 http://www.j.sinap.ac.cn/hejishu/CN/abstract/abstract570.shtml .
Wending FAN , Guangyao SUN , Binhang ZHANG , . Research and verification of Monte Carlo burnup calculations based on Chebyshev rational approximation method . Nuclear Techniques , 2016 . 39 ( 4 ): 040501 DOI: 10.11889/j.0253-3219.2016.hjs.39.040501 http://doi.org/10.11889/j.0253-3219.2016.hjs.39.040501 http://www.j.sinap.ac.cn/hejishu/CN/abstract/abstract570.shtml .
吴 明宇 , 王 事喜 , 杨 勇 , . 回溯算法在燃耗计算中的应用 . 原子能科学技术 , 2013 . 47 ( 7 ): 1127 - 1132 . DOI: 10.7538/yzk.2013.47.07.1127 http://doi.org/10.7538/yzk.2013.47.07.1127 .
Mingyu WU , Shixi WANG , Yong YANG , . Application of backtracking algorithm to depletion calculations . Atomic Energy Science and Technology , 2013 . 47 ( 7 ): 1127 - 1132 . DOI: 10.7538/yzk.2013.47.07.1127 http://doi.org/10.7538/yzk.2013.47.07.1127 .
W J Cody , G Meinardus , R S Varga . Chebyshev rational approximations to exp (-x) in[0, +∞) and applications to heat-conduction problems . Journal of Approximation Theory , 1969 . 2 ( 1 ): 50 - 65 . DOI: 10.1016/0021-9045(69)90030-6 http://doi.org/10.1016/0021-9045(69)90030-6 .
Schmelzer T. Carathéodory-Fejér approximation, MATLAB central[CP/OL]. 2008. http://www.mathworks.com/matlabcentral http://www.mathworks.com/matlabcentral .
A J Carpenter , A Ruttan , R S Varga . Extended numerical computations on the"1/9"conjecture in rational approximation theory . Lecture Notes in Mathematics , 1984 . 1105 ( 4 ): 383 - 411 . DOI: 10.1007/BFb0072427 http://doi.org/10.1007/BFb0072427 .
Pusa M. Correction to partial fraction decomposition coefficients for Chebyshev rational approximation on the negative real axis[J/OL]. arXiv: 1206. 2880[math. NA], 2012. http://www.oalib.com/paper/4122717
张 竞宇 , 马 亚栋 , 陈 义学 , . CRAM在放射性核素存量计算中的应用 . 核技术 , 2017 . 40 ( 8 ): 080502 DOI: 10.11889/j.0253-3219.2017.hjs.40.080502 http://doi.org/10.11889/j.0253-3219.2017.hjs.40.080502 http://www.j.sinap.ac.cn/hejishu/CN/abstract/abstract799.shtml .
Yadong ZHANG , Yadong MA , Yixue CHEN , . Research and verification of Monte Carlo burnup calculations based on Chebyshev rational approximation method . Nuclear Techniques , 2017 . 40 ( 8 ): 080502 DOI: 10.11889/j.0253-3219.2017.hjs.40.080502 http://doi.org/10.11889/j.0253-3219.2017.hjs.40.080502 http://www.j.sinap.ac.cn/hejishu/CN/abstract/abstract799.shtml .
C Moler , C F Loan . Nineteen dubious ways to compute the exponential of a matrix, twenty-five years later . SIAM Review , 2003 . 45 ( 1 ): 3 - 49 . DOI: 10.1137/S00361445024180 http://doi.org/10.1137/S00361445024180 .
N J Higham . The scaling and squaring method for the matrix exponential revisited . SIAM Review , 2009 . 51 ( 4 ): 747 - 764 . DOI: 10.1137/090768539 http://doi.org/10.1137/090768539 .
R B Sidje . Expokit:a software package for computing matrix exponentials . ACM Transactions on Mathematical Software , 1998 . 24 ( 1 ): 130 - 156 . DOI: 10.1145/285861.285868 http://doi.org/10.1145/285861.285868 .
A Yamamoto , M Tatsumi , N Sugimura . Numerical solution of stiff burnup equation with short half lived nuclides by the Krylov Subspace method . Journal of Nuclear Science and Technology , 2007 . 44 ( 2 ): 147 - 154 . DOI: 10.1080/18811248.2007.9711268 http://doi.org/10.1080/18811248.2007.9711268 .
R B James , J R Donald . Sparse matrix computations , : New York Academic Press , 1976 . 3 - 22 . DOI: 10.1016/B978-0-12-141050-6.50006-4 http://doi.org/10.1016/B978-0-12-141050-6.50006-4 .
M Pusa , J Leppänen . Solving linear systems with sparse gaussian elimination in the Chebyshev rational approximation method . Nuclear Science and Engineering , 2013 . 175 ( 3 ): 250 - 258 . DOI: 10.13182/nse12-52 http://doi.org/10.13182/nse12-52 .
0
Views
22
下载量
1
CSCD
- L-PDFDownload
- Meta-XMLDownload
- BibTeXDownload
- EndnoteDownload
- RefMan Download
- NoteExpressDownload
- NoteFirstDownload
Publicity Resources
Related Articles
Related Author
Related Institution
AI问答- Address：2019 Jialuo Road, Jiading, Shanghai, China Postal code：201800
- Tel：021-39194977 Email：LHB@sinap.ac.cn
- Technical support is provided by Beijing Founder electronics co., LTD 沪ICP备05005479号-1
京公网安备11010802024621 - It is recommended to read the content of this site in Chrome&IE9+. Please switch to extreme mode in browser 360.
- Cookies We use cookies to help provide and enhance our service and tailor content. By continuing, you agree to the use of cookies.
- Total Visits：880473
