基于interFoam的新求解器,加入蒸发模型后不收敛

大家好奥,
我目前的课题是基于OpemFOAM进行液滴蒸发的研究,VOF模型中的表面张力模型常常会引起非物理的spurious currents出现在液滴的界面附近,导致计算不稳定等一系列问题,我前期在interFoam中植入了Sharp Surface-tension Force (SSF)和Filtering Surface-tension Force (FSF),发现在没有相变的时候,可以在很短的时间内,比如0.001s内将Spurious Currents降到1.0E-8。

但是,一旦加入蒸发模型后,求解就变得不稳定了,模拟不收敛,而且模拟过程中,液滴逐渐被spurious currents拖拽偏离平衡位置。具体的大家可以参考我的两个图片,是只有spurious currents的速度场,欢迎大家的讨论。
01