三相搅拌反应器的CFD模拟

waouyouy

各位虫友好，我现在从事三相搅拌反应器的CFD模拟，由于是刚入门的新手，所以很多东西都不知道，因此想请教各位！
1，反应器模型为双层搅拌器，因此利用MRF设置了两个旋转坐标系。
2，因为想模拟反应器的连续工况，选用了E-E-E模型和k-ε湍流模型。
3，模型网格采用非结构网格（旋转坐标系的网格尺寸与外围惯性坐标系网格尺度之比大概为2.5：1到6：1）
采用非稳态方法模拟连续工况，但是计算开始几步就发散了，见下图。在调节欠松弛因子后也没改观，我在想是不是初值的问题和网格质量有问题（在gambit里检查倾斜度没有超出要求）。希望各位道友能给与在下一些建议与支持，在此谢过！

howardD

聽你這樣說應該是沒有太大問題

請先調整time scale的大小，降低conrant number的大小囉！

waouyouy

欧拉-欧拉方法和非稳态好像不能条conrant number。另外，我如果不设置液位，让三相充满整个反应器，收敛是没问题的，但是以patch出液位（上面有部分全是气相），就出现发散的问题，所以有没有可能是我patch初始化的问题了？