之前算DEM反应一直不收敛,后来把计算规模减少到两个固定颗粒,也就是仿照reactive-chemistry和evaperation这两个算例. 目前reactive-chemistry这个算例,把各种输运方程都打开以后,DES_ONEWAY_COUPLED设置为False,就会出现不收敛的情况,而evaperation换成同样的网格和几何,就没有这个情况,于是索性在evaperation上改. 经过尝试发现,反应焓值应该是导致容易发散的一个原因.我设置的反应如下: 如果把异相反应的反应焓设为0