化学非平衡流的求解问题
最近在编写求解化学非平衡流的程序，遇见一个问题，我求解氢氧燃烧的时候，处理化学源项采用流动与化学反应解耦的算法，然后在单独求解化学源项的时候冻结流动，化学时间步采用流动时间步的1/n，但是利用不同反应机理的时候，发现采用7组分8反应机理程序能够运行，但是换成9组分19步反应的时候，前面所用的1/n就不能用了，在第3个时间步，出现了密度为负的非物理解，这是不是由于不用反应机理中某些组分的化学反应特征时间不同造成的？我试着将n变为原来的5倍，算一段时间才出现这种状况，是不是这是由于上述特征时间短造成的严重的刚性问题导致的？ 好高深
页:
[1]