GEAR算法在随机轨道模型计算中的应用

本文对随机轨道模型中颗粒相常微分方程组的刚性问题进行了分析，结果表明：当采用常规算法如四阶Runge－kutta法求解方程组时，方程组的刚性是导致某些情况下计算发散或计算时间过长的原因。为此，本文将适用于求解刚性方程组的Gear算法应用于随机轨道模型的计算中，取得了良好的效果．