Computing the Matrix Exponential in Burnup Calculations

Maria Pusa, Jaakko Leppänen

    Research output: Contribution to journalArticleScientificpeer-review

    124 Citations (Scopus)

    Abstract

    The topic of this paper is the computation of the matrix exponential in the context of burnup equations. The established matrix exponential methods are introduced briefly. The eigenvalues of the burnup matrix are important in choosing the matrix exponential method, and their characterization is considered. Based on the characteristics of the burnup matrix, the Chebyshev rational approximation method (CRAM) and its interpretation as a numeric contour integral are discussed in detail. The introduced matrix exponential methods are applied to two test cases representing an infinite pressurized water reactor pin-cell lattice, and the numerical results are presented. The results suggest that CRAM is capable of providing a robust and accurate solution to the burnup equations with a very short computation time.
    Original languageEnglish
    Pages (from-to)140-150
    Number of pages11
    JournalNuclear Science and Engineering
    Volume164
    Issue number2
    DOIs
    Publication statusPublished - 2010
    MoE publication typeA1 Journal article-refereed

      Fingerprint

    Cite this