最近在开展CO/CO2甲烷化反应的模拟,通过UDF植入了多种动力学方程,顺便考察了Species模型中Chemistry Solver选项。在这里分享下自己的体会,正确性不保证。
None-Explicit Source,收敛性差,对于一些动力学方程,开始计算便发散。
Stiff Chemistry Solver,收敛性好。试了几种动力学,都能算的下去,但其收敛时的稳定值与平衡组成偏差较大。
另外有个小技巧,使用Stiff开始计算,待基本达到稳定后,转为Explicit续算。但是,以上两种Solver的稳定值都取决于动力学方程。举个例子,
对于动力学方程A,首先使用Stiff计算,稳定时转化率约为96%,转而Explicit续算,转化率达到99%,且计算过程中不存在结果的波动;若一开始就使用Explicit计算,则发散;若使用Stiff一直算,同样能达到平衡组成。
对于动力学方程B,首先使用Stiff计算,稳定时转为Explicit续算后,仍然未达到平衡组成,且计算结果波动,如下图。两者所得转化率一致。
