Numerical methods for problems of chemical kinetics with diffusion

It is proposed to use the method of lines (MOL) for numerical solving of a problem of chemical kinetics with diffusion and heat conduction. It leads to a stiff system of a very large number of ordinary differential equations (ODEs). The system may be turned to suchform, that its Jacobian matrix has band structure (the width of this band is twice as much as the quantity of the chemical components). Different codes (LSODE, ROS4, RODAS, SEULEX and COCOSE) were compared in order to check which of them is more convenient for solution of such kind of problems. This approach was studied on the problem of radical copolymerization. The accuracy of results was investigated carefully by the nested grids method. The calculations show that ail these codes are stable and reliable and at a modest amount of calculations give good accuracy, that is quite sufficient for purposes of chemical practice.