求刚性中子动力学方程数值解的新方法 (2003年)

中子动力学方程是刚性方程,数值计算很费机时,为减少计算机时提出了许多算法,常用的有隐式差分法、多步法和自动变步长法等,著名的方法有吉尔方法和埃米特方法。但这些方法的计算程序复杂且总计算时间减少不明显。消去刚性法是一个全新的算法,稳定步长可达到 20 s,计算时间大大下降。