液固多相模拟遇到的奇怪问题

waouyouy

各位好！
本人利用双层搅拌器模拟气液固三相搅拌，暂时先不通气体，只patch出了固体床层高度和液位高度，然后进行非稳态计算。
所选模型为E-E，SKE湍流模型，并且没加额外源项，算法采用压力-速度耦合的SIMPLE算法，对流项的离散格式均为一阶迎风，时间步长0.001和0.0005都尝试过，能收敛（能算到2.2s左右，此时的固相分布及水的流场图1，2），之后就会出现solid相w方向上速度达不到收敛残差限（见数据图3，4），并且在计算一小段时间后就发散了，软件提示如下图。
让我很奇怪的是，第一个时间步长内，固相的体积分数残差相vf一直维持在10e-1这个级别（见图5），到第二个时间步长内反而直接从10e-1降到10e-5，我不知道这是不是一个反常的现象。
想咨询下，这种情况是什么引起的，如何能解决这种情况？
拜谢各位！！！

[ 本帖最后由 waouyouy 于 2013-12-7 13:22 编辑 ]

waouyouy

各位大神速速现身啊，小弟纠结这个问题两周了，一直没弄清楚，求指教求指教！！

waouyouy

没有人做过类似的模拟吗？急急！！