[求助]化学反应动力学问题！各位请进！！！急急急

shhj

本人FORTRAN 编程计算乙醇燃烧的化学反应动力学过程，有372个基元反应，57种物质，用GEAR算法，但是不收敛，请问是什么问题，如果是JCAOBI矩阵不可逆造成的吗？有什么办法让他可逆？
谢谢

shyang5

建议用龙格-库塔试一试

sunmoon

此乃化学反应刚性所致，R-K方法除非步长极小，否则无用，而步长极小往往导致计算过缓。其实，Gear算法是解决此类问题的最好方法，不知是你自己写的吗？一般来说，自己难以写出较好的Gear算法程序，建议你去网上搜搜，也许有收获。