之前算DEM反应一直不收敛,后来把计算规模减少到两个固定颗粒,也就是仿照reactive-chemistry和evaperation这两个算例。
目前reactive-chemistry这个算例,把各种输运方程都打开以后,DES_ONEWAY_COUPLED设置为False,就会出现不收敛的情况,而evaperation换成同样的网格和几何,就没有这个情况,于是索性在evaperation上改。
经过尝试发现,反应焓值应该是导致容易发散的一个原因。我设置的反应如下:
如果把异相反应的反应焓设为0,则反应可以顺利进行,不会导致发散(虽然这样并不合理)。
接下来还需要继续寻找解决办法,感觉应该是反应速率编写的不够合理。
时间: 2024-11-04 19:53:38