A new Krylov subspace method based on rational approximation to solve stiff burnup equation

作者：Li, Xuezhong; Cai, Jiejin^*

来源：Annals of Nuclear Energy, 2018, 118: 99-106.

DOI：10.1016/j.anucene.2018.04.005

• Burnup equation
• Matrix exponential
• Krylov Subspace Method
• GMRES
• Rational approximation

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.

• 出版日期2018-8
• 单位华南理工大学