非齐次燃耗方程数值解法

非齐次燃耗方程常用于描述具有显著核素迁移效应的核能系统中核素含量随时间的变化规律.国内外许多燃耗计算程序无法求解方程非齐次项含时的情况.本文在方程非齐次项能够被有限阶关于时间的泰勒展开逼近这一前提下,研究了非齐次项含时情况下方程的两种解法.首先通过Laplace变换推导出了方程基于线性子链方法的解析解形式,然后使用Carathéodory-Fejér方法计算出了方程矩阵级数解的近最佳Chebyshev有理逼近式.将两种方法在燃耗计算程序JBURN中实现,并进行了数值计算,绝大部分计算结果符合很好,部分结果在较长有效数字内仍能保持一致,验证了方法的正确性和精度.同时为求解具有其他非齐次项形式的燃耗方程提供了一种思路.