Accuracy considerations for Chebyshev rational approximation method (CRAM) in Burnup calculations
Conference
·
OSTI ID:22212772
- 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
- Country of Publication:
- United States
- Language:
- English
Similar Records
An efficient implementation of the Chebyshev Rational Approximation Method (CRAM) for solving the burnup equations
Using Laguerre polynomials to compute the matrix exponential in burnup calculations
Conference
·
Sun Jul 01 00:00:00 EDT 2012
·
OSTI ID:22105637
Using Laguerre polynomials to compute the matrix exponential in burnup calculations
Conference
·
Sun Jul 01 00:00:00 EDT 2012
·
OSTI ID:22105604