核电计算中刚性常微分方程的数值解

在核电计算中,核岛中的反应堆动态学方程和常规岛中的湿蒸汽汽轮机甩负荷时各级压力变化方程都是典型的刚性常微分方程。本文介绍用通用性强的Gear方法求解这些刚性方程的有效性。在解点堆动态学方程时,Gear方法的计算效率与分段多项式逼近法相当;在对矩阵求逆的算法作了改进后,该方法的计算效率就高于分段多项式逼近法。