气液固三相模拟k-epsilon发散

[青山]

目前在做一个气液固三相的搅拌模拟，我想在搅拌釜中先setFields一个有20%固含率和80%液体的液位，上部是空气，反应器底部有气体鼓入，使用openfoam9的multiphaseEulerFoam，气液湍流使用的是k-e模型，固体湍流使用kineticTheory模型。目前遇到的问题是：当我给上部气相设成纯气体的话，直接运算不起来，第一个时间步都跑不了就结束了。只有给上部气体设一个很小的液体含量才可以正常运行，但是这样的话，计算出来结果是液位上的k.water,k.gas,epsilon.water,epsilon.gas都会出现一块很大的值，计算过程中会出现k和epsilon的bounding,然后导致计算发散，如下图


alpha云图看着到是挺正常，

离散格式也都使用了迎风格式，

使用层流的话计算是没有任何问题的，而且如果不设置上层空气，使用液体填满的话，计算也是没问题的，一旦设置自由液位就有问题。而且为了简化网格中间搅拌桨我使用了一块平板代替，网格是很规整的六面体网格，所以应该也不是网格问题，请问各位大佬有没有遇到过这样的问题，望不吝赐教！

[青山]

会不会是对于这种上层层流，下层湍流的存在较大湍动差异的案例，k-e模型不适用？k-e模型应该使用在完全发展的湍流案例中是吗？

[李东岳]

湍流模型，以及drag模型需要调试一下，这个就是方程奇异的问题。另外fvSolution那面可以不挂预条件试一下