The evaluation of numerical techniques for solution of stiff ordinary differential equations arising from chemical kinetic problems

摘要: Abstract Ordinary differential equations with widely scattered eigenvalues (stiff O.D.E.'s) occur often in the studies of reaction network problems. Five numerical methods, including two methods based on Backward differentiation formulas, a modified Runge-Kutta-Fehlberg method, method PECE Adams and an improved semi-implicit Euler are evaluated by comparing their performance when applied to test systems. The systems represent different combinations linearity nonlinearity, small large dimension, real complex eigenvalues, slightly stiff very relative merits dificiencies discussed.