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

Application of Diagonally Implicit Runge Kutta in Point Kinetics

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

机构地区：[1]中国核动力研究设计院核反应堆系统设计技术重点实验室,成都610213 [2]清华大学工程物理系,北京100084

出　　处：《核动力工程》2016年第1期13-17,共5页Nuclear Power Engineering

摘　　要：点堆动力学方程刚性比较强,采用常规全隐式龙格库塔方法则求解耗时多。对角隐式龙格库塔方法保留全隐式龙格库塔善于求解刚性方程的特点,同时又大大降低计算量。通过嵌入低阶龙格库塔方法,实现自适应时间步技术,提高计算效率。通过计算阶跃、线性、正弦3种反应性变化基准题,计算结果表明该方法和其他方法结果符合很好,而且相对于θ方法能够在相同的计算时间内给出更加精确的解,特别是在快速插入线性反应性的情况下。The stiffness of point kinetics is very strong, while it usually takes a lot of time to solve by using full implicit Runge Kutta （FIRK）. Diagonally Implicit Runge Kutta （DIRK） is a useful tool like FIRK to solve the stiff differential equations, while it greatly reduces the computation compared with FIRK. By embedded lower order Runge Kutta, the time step-size adaptation technique is implemented, which improves the computation efficiency of DIRK. Through three typical cases with step, ramp and sinusoidal reactivity insertions, it shows that the results obtained by DIRK are in good agreement with other available results and DIRK can give far more accurate results than θ method at the same computation cost, especially in the case of fast ramp reactivity insertion.

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

分 类 号：TL329[核科学技术—核技术及应用]