Higher-order Chebyshev rational approximation method and application to burnup equations

Maria Pusa

Research output: Contribution to journalArticleScientificpeer-review

15 Citations (Scopus)

Abstract

The burnup equations can, in principle, be solved by computing the exponential of the burnup matrix. However, the problem is extremely stiff, and the matrix exponential solution was long considered infeasible for entire burnup systems containing short-lived nuclides. After discovering that the eigenvalues of burnup matrices are confined to the vicinity of the negative real axis, the Chebyshev rational approximation method (CRAM) was introduced for solving the burnup equations and it was shown to be capable of providing accurate and efficient solutions without the need to exclude the short-lived nuclides. The main difficulty in using CRAM is determining the coefficients of the rational approximant for a given approximation order, with the previously published coefficients enabling only approximations up to order 16 for computing the matrix exponential. In this paper, a Remez-type method is presented for the computation of higher-order CRAM approximations. The optimal form of CRAM for the solution of burnup equations is discussed, and the method of incomplete partial fractions is proposed for this purpose. The CRAM coefficients based on this factorization are provided for approximation orders 4, 8, 12, . . ., 48. The accuracy of the method is demonstrated by applying it to large burnup and decay systems. It is shown that higher-order CRAM can be used to solve the burnup equations accurately for time steps of the order of 1 million years.

Original languageEnglish
Pages (from-to)297-318
Number of pages22
JournalNuclear Science and Engineering
Volume182
Issue number3
DOIs
Publication statusPublished - 1 Jan 2016
MoE publication typeA1 Journal article-refereed

Fingerprint

Isotopes
Factorization

Keywords

  • Burnup equations
  • CRAM
  • Matrix exponential

Cite this

@article{89826c3e5fe948acbb54f52f5dc62b13,
title = "Higher-order Chebyshev rational approximation method and application to burnup equations",
abstract = "The burnup equations can, in principle, be solved by computing the exponential of the burnup matrix. However, the problem is extremely stiff, and the matrix exponential solution was long considered infeasible for entire burnup systems containing short-lived nuclides. After discovering that the eigenvalues of burnup matrices are confined to the vicinity of the negative real axis, the Chebyshev rational approximation method (CRAM) was introduced for solving the burnup equations and it was shown to be capable of providing accurate and efficient solutions without the need to exclude the short-lived nuclides. The main difficulty in using CRAM is determining the coefficients of the rational approximant for a given approximation order, with the previously published coefficients enabling only approximations up to order 16 for computing the matrix exponential. In this paper, a Remez-type method is presented for the computation of higher-order CRAM approximations. The optimal form of CRAM for the solution of burnup equations is discussed, and the method of incomplete partial fractions is proposed for this purpose. The CRAM coefficients based on this factorization are provided for approximation orders 4, 8, 12, . . ., 48. The accuracy of the method is demonstrated by applying it to large burnup and decay systems. It is shown that higher-order CRAM can be used to solve the burnup equations accurately for time steps of the order of 1 million years.",
keywords = "Burnup equations, CRAM, Matrix exponential",
author = "Maria Pusa",
year = "2016",
month = "1",
day = "1",
doi = "10.13182/NSE15-26",
language = "English",
volume = "182",
pages = "297--318",
journal = "Nuclear Science and Engineering",
issn = "0029-5639",
publisher = "American Nuclear Society ANS",
number = "3",

}

Higher-order Chebyshev rational approximation method and application to burnup equations. / Pusa, Maria.

In: Nuclear Science and Engineering, Vol. 182, No. 3, 01.01.2016, p. 297-318.

Research output: Contribution to journalArticleScientificpeer-review

TY - JOUR

T1 - Higher-order Chebyshev rational approximation method and application to burnup equations

AU - Pusa, Maria

PY - 2016/1/1

Y1 - 2016/1/1

N2 - The burnup equations can, in principle, be solved by computing the exponential of the burnup matrix. However, the problem is extremely stiff, and the matrix exponential solution was long considered infeasible for entire burnup systems containing short-lived nuclides. After discovering that the eigenvalues of burnup matrices are confined to the vicinity of the negative real axis, the Chebyshev rational approximation method (CRAM) was introduced for solving the burnup equations and it was shown to be capable of providing accurate and efficient solutions without the need to exclude the short-lived nuclides. The main difficulty in using CRAM is determining the coefficients of the rational approximant for a given approximation order, with the previously published coefficients enabling only approximations up to order 16 for computing the matrix exponential. In this paper, a Remez-type method is presented for the computation of higher-order CRAM approximations. The optimal form of CRAM for the solution of burnup equations is discussed, and the method of incomplete partial fractions is proposed for this purpose. The CRAM coefficients based on this factorization are provided for approximation orders 4, 8, 12, . . ., 48. The accuracy of the method is demonstrated by applying it to large burnup and decay systems. It is shown that higher-order CRAM can be used to solve the burnup equations accurately for time steps of the order of 1 million years.

AB - The burnup equations can, in principle, be solved by computing the exponential of the burnup matrix. However, the problem is extremely stiff, and the matrix exponential solution was long considered infeasible for entire burnup systems containing short-lived nuclides. After discovering that the eigenvalues of burnup matrices are confined to the vicinity of the negative real axis, the Chebyshev rational approximation method (CRAM) was introduced for solving the burnup equations and it was shown to be capable of providing accurate and efficient solutions without the need to exclude the short-lived nuclides. The main difficulty in using CRAM is determining the coefficients of the rational approximant for a given approximation order, with the previously published coefficients enabling only approximations up to order 16 for computing the matrix exponential. In this paper, a Remez-type method is presented for the computation of higher-order CRAM approximations. The optimal form of CRAM for the solution of burnup equations is discussed, and the method of incomplete partial fractions is proposed for this purpose. The CRAM coefficients based on this factorization are provided for approximation orders 4, 8, 12, . . ., 48. The accuracy of the method is demonstrated by applying it to large burnup and decay systems. It is shown that higher-order CRAM can be used to solve the burnup equations accurately for time steps of the order of 1 million years.

KW - Burnup equations

KW - CRAM

KW - Matrix exponential

UR - http://www.scopus.com/inward/record.url?scp=84975044074&partnerID=8YFLogxK

U2 - 10.13182/NSE15-26

DO - 10.13182/NSE15-26

M3 - Article

VL - 182

SP - 297

EP - 318

JO - Nuclear Science and Engineering

JF - Nuclear Science and Engineering

SN - 0029-5639

IS - 3

ER -