Accuracy considerations for Chebyshev rational approximation method (CRAM) in Burnup calculations
- VTT Technical Research Centre of Finland, P.O. Box 1000, FI-02044 VTT (Finland)
The burnup equations can in principle be solved by computing the exponential of the burnup matrix. However, due to the difficult numerical characteristics of burnup matrices, the problem is extremely stiff and the matrix exponential solution has previously been considered infeasible for an entire burnup system containing over a thousand nuclides. It was recently discovered by the author that the eigenvalues of burnup matrices are generally located near the negative real axis, which prompted introducing the Chebyshev rational approximation method (CRAM) for solving the burnup equations. CRAM can be characterized as the best rational approximation on the negative real axis and it has been shown to be capable of simultaneously solving an entire burnup system both accurately and efficiently. In this paper, the accuracy of CRAM is further studied in the context of burnup equations. The approximation error is analyzed based on the eigenvalue decomposition of the burnup matrix. It is deduced that the relative accuracy of CRAM may be compromised if a nuclide concentration diminishes significantly during the considered time step. Numerical results are presented for two test cases, the first one representing a small burnup system with 36 nuclides and the second one a full a decay system with 1531 nuclides. (authors)
- Research Organization:
- American Nuclear Society, 555 North Kensington Avenue, La Grange Park, IL 60526 (United States)
- OSTI ID:
- 22212772
- Resource Relation:
- Conference: M and C 2013: 2013 International Conference on Mathematics and Computational Methods Applied to Nuclear Science and Engineering, Sun Valley, ID (United States), 5-9 May 2013; Other Information: Country of input: France; 9 refs.; Related Information: In: Proceedings of the 2013 International Conference on Mathematics and Computational Methods Applied to Nuclear Science and Engineering - M and C 2013| 3016 p.
- Country of Publication:
- United States
- Language:
- English
Similar Records
Using Laguerre polynomials to compute the matrix exponential in burnup calculations
Improving the Accuracy of the Chebyshev Rational Approximation Method Using Substeps