用高阶泰勒多项式积分方法求解点堆中子动力学方程

摘要

在求解点堆中子动力学方程组中,对中子密度使用分段全隐式高阶泰勒多项式近似技术,给出一求解点堆中子动力学方程组的数值积分方法,并对该方法进行了修正优化.实例计算并与传统的三阶Hermite插值多项式法的比较表明:该方法能显著消除刚性方程组带来的数值计算不利因素,对给定的反应性输入能够取得较高精确度的数值结果,计算过程简洁,且计算速度快,通过对高阶泰勒多项式的修正,计算精度有了进一步提高,可适宜于反应堆中子动力学控制的设计分析和仿真计算.