Численные методы решения задач химической кинетики с диффузией

Предлагается использовать метод прямых (МП) для численного решения задач химической кинетики с диффузией и теплопроводностью. При таком подходе получается система обыкновенных дифференциальных уравнений (ОДУ), жесткая и достаточно большой размерности. Матрица Якоби этой системы имеет ленточную структуру (половина ширины которой равна количеству химических компонент). Для выбора более подходящего алгоритма решения такого типа задач было проведено сравнение наиболее известных программ (LSODE, ROS4, RODAS, SEULEX и COCOSE), ориентированных на решение жестких систем ОДУ. Конкретно была взята задача полной сополимеризации. Точность результатов была тщательно исследована с помощью метода вложенных сеток. Тестирование показало, что все эти программы устойчивы и надежны, и при разумных вычислительных затратах дают решение с хорошей точностью, вполне достаточной для практики.