A new Krylov subspace method based on rational approximation to solve stiff burnup equation. (August 2018)

• Record Type: Journal Article
• Authors: Li, Xuezhong; Cai, Jiejin
• Is Part Of: Annals of nuclear energy. Volume 118(2018)
• Page Start: 99
• Page End: 106
• Publication Date: 2018-08
• Subjects: Burnup equation -- Matrix exponential -- Krylov Subspace Method -- GMRES -- Rational approximation; Nuclear energy -- Periodicals; Nuclear engineering -- Periodicals; 621.4805
• DOI: 10.1016/j.anucene.2018.04.005
• Languages: English
• ISSNs: 0306-4549

Abstract

Highlights: A new Krylov Subspace Method is proposed to solve burnup equation. A rational approximation method is implemented in the new method. Calculation cases show the proposed method has great advantage in computation speed. Abstract: A burnup equation can be solved with matrix exponential method and its solution can be written as n ( t ) = e At n ( 0 ) . In burnup calculation, general Krylov Subspace Method can solute a matrix–vector efficiently in a subspace but fails to keep a high precision. To solve this problem, a new kind of Krylov Subspace Method, Generalized Minimal Residual Method (GMRES) is implemented, based on a rational approximation method. It shows its great advantage in computation speed, which is more than four times faster than the same kind of rational approximation solved in a whole space while its accuracy is also guaranteed. Some optimizations, such as shift-Invariant technique, precondition technique and restart technique, are also implemented on burnup calculation.