Using Laguerre polynomials to compute the matrix exponential in burnup calculations
Abstract
An essential part of burnup analysis is to solve the burnup equations. The burnup equations can be regarded as a first-order linear system and solved by means of matrix exponential methods. Because of its large spectrum, it is difficult to compute the exponential of the burnup matrix. Conventional methods of computing the matrix exponential, such as the truncated Taylor expansion and the Pade approximation, are not applicable to burnup calculations. Recently the Chebyshev Rational Approximation Method (CRAM) has been applied to solve burnup matrix exponential and shown to be robust and accurate. However, the main defect of CRAM is that its coefficients are not easy to obtain. In this paper, an orthogonal polynomial expansion method, called Laguerre Polynomial Approximation Method (LPAM), is proposed to compute the matrix exponential in burnup calculations. The polynomial sequence of LPAM can be easily computed in any order and thus LPAM is quite convenient to be utilized into burnup codes. Two typical test cases with the decay and cross-section data taken from the standard ORIGEN 2.1 libraries are calculated for validation, against the reference results provided by CRAM of 14 order. Numerical results show that, LPAM is sufficiently accurate for burnup calculations. The influences ofmore »
- Authors:
-
- Dept. of Engineering Physics, Tsinghua Univ., Beijing, 100084 (China)
- Publication Date:
- Research Org.:
- American Nuclear Society, Inc., 555 N. Kensington Avenue, La Grange Park, Illinois 60526 (United States)
- OSTI Identifier:
- 22105604
- Resource Type:
- Conference
- Resource Relation:
- Conference: PHYSOR 2012: Conference on Advances in Reactor Physics - Linking Research, Industry, and Education, Knoxville, TN (United States), 15-20 Apr 2012; Other Information: Country of input: France; 14 refs.
- Country of Publication:
- United States
- Language:
- English
- Subject:
- 97 MATHEMATICAL METHODS AND COMPUTING; 22 GENERAL STUDIES OF NUCLEAR REACTORS; BURNUP; CROSS SECTIONS; EQUATIONS; EXPANSION; LAGUERRE POLYNOMIALS; MATRICES; NUCLEAR DATA COLLECTIONS; PADE APPROXIMATION
Citation Formats
She, D., Zhu, A., and Wang, K. Using Laguerre polynomials to compute the matrix exponential in burnup calculations. United States: N. p., 2012.
Web.
She, D., Zhu, A., & Wang, K. Using Laguerre polynomials to compute the matrix exponential in burnup calculations. United States.
She, D., Zhu, A., and Wang, K. 2012.
"Using Laguerre polynomials to compute the matrix exponential in burnup calculations". United States.
@article{osti_22105604,
title = {Using Laguerre polynomials to compute the matrix exponential in burnup calculations},
author = {She, D. and Zhu, A. and Wang, K.},
abstractNote = {An essential part of burnup analysis is to solve the burnup equations. The burnup equations can be regarded as a first-order linear system and solved by means of matrix exponential methods. Because of its large spectrum, it is difficult to compute the exponential of the burnup matrix. Conventional methods of computing the matrix exponential, such as the truncated Taylor expansion and the Pade approximation, are not applicable to burnup calculations. Recently the Chebyshev Rational Approximation Method (CRAM) has been applied to solve burnup matrix exponential and shown to be robust and accurate. However, the main defect of CRAM is that its coefficients are not easy to obtain. In this paper, an orthogonal polynomial expansion method, called Laguerre Polynomial Approximation Method (LPAM), is proposed to compute the matrix exponential in burnup calculations. The polynomial sequence of LPAM can be easily computed in any order and thus LPAM is quite convenient to be utilized into burnup codes. Two typical test cases with the decay and cross-section data taken from the standard ORIGEN 2.1 libraries are calculated for validation, against the reference results provided by CRAM of 14 order. Numerical results show that, LPAM is sufficiently accurate for burnup calculations. The influences of the parameters on the convergence of LPAM are also discussed. (authors)},
doi = {},
url = {https://www.osti.gov/biblio/22105604},
journal = {},
number = ,
volume = ,
place = {United States},
year = {Sun Jul 01 00:00:00 EDT 2012},
month = {Sun Jul 01 00:00:00 EDT 2012}
}