An efficient implementation of the Chebyshev Rational Approximation Method (CRAM) for solving the burnup equations
- VTT Technical Research Centre of Finland, P.O. Box 1000, FI-02044 VTT (Finland)
The Chebyshev Rational Approximation Method (CRAM) has been recently introduced by the authors for solving the burnup equations with excellent results. This method has been shown to be capable of simultaneously solving an entire burnup system with thousands of nuclides both accurately and efficiently. The method was prompted by an analysis of the spectral properties of burnup matrices and it can be characterized as the best rational approximation on the negative real axis. The coefficients of the rational approximation are fixed and have been reported for various approximation orders. In addition to these coefficients, implementing the method only requires a linear solver. This paper describes an efficient method for solving the linear systems associated with the CRAM approximation. The introduced direct method is based on sparse Gaussian elimination where the sparsity pattern of the resulting upper triangular matrix is determined before the numerical elimination phase. The stability of the proposed Gaussian elimination method is discussed based on considering the numerical properties of burnup matrices. Suitable algorithms are presented for computing the symbolic factorization and numerical elimination in order to facilitate the implementation of CRAM and its adoption into routine use. The accuracy and efficiency of the described technique are demonstrated by computing the CRAM approximations for a large test case with over 1600 nuclides. (authors)
- Research Organization:
- American Nuclear Society, Inc., 555 N. Kensington Avenue, La Grange Park, Illinois 60526 (United States)
- OSTI ID:
- 22105637
- 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; 8 refs.
- Country of Publication:
- United States
- Language:
- English
Similar Records
Using Laguerre polynomials to compute the matrix exponential in burnup calculations
Solving the Bateman equations in CASMO5 using implicit ode numerical methods for stiff systems