对角隐式龙格库塔法在点堆动力学中的应用

来源：核动力工程杂志2016年第1期北京时间：

作者：张知竹, 蔡云, 彭星杰, 李庆, 秦冬

单位：1. 中国核动力研究设计院核反应堆系统设计技术重点实验室, 成都, 610213;；2. 清华大学工程物理系, 北京, 100084；

摘要：点堆动力学方程刚性比较强,采用常规全隐式龙格库塔方法则求解耗时多。对角隐式龙格库塔方法保留全隐式龙格库塔善于求解刚性方程的特点,同时又大大降低计算量。通过嵌入低阶龙格库塔方法,实现自适应时间步技术,提高计算效率。通过计算阶跃、线性、正弦3种反应性变化基准题,计算结果表明该方法和其他方法结果符合很好,而且相对于<i>θ</i>方法能够在相同的计算时间内给出更加精确的解,特别是在快速插入线性反应性的情况下。

关键词：点堆动力学, 对角隐式龙格库塔方法, 自适应时间步