1. Computing the Matrix Exponential in Burnup Calculations
- Author
-
Jaakko Leppänen and Maria Pusa
- Subjects
State-transition matrix ,Approximation theory ,010308 nuclear & particles physics ,Computation ,MathematicsofComputing_NUMERICALANALYSIS ,0211 other engineering and technologies ,Context (language use) ,02 engineering and technology ,Nuclear reactor ,01 natural sciences ,law.invention ,Nuclear Energy and Engineering ,law ,ComputingMethodologies_SYMBOLICANDALGEBRAICMANIPULATION ,0103 physical sciences ,Applied mathematics ,021108 energy ,Matrix exponential ,Algorithm ,Eigenvalues and eigenvectors ,Mathematics ,Burnup - 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.
- Published
- 2010
- Full Text
- View/download PDF