FLUENT欧拉-欧拉模型模拟双组份鼓泡流化床中的一些问题
-
尊敬的各位老师，你们好，我是一名研究生，正在做双组份鼓泡流化床的模拟，运用fluent里双欧拉模型。有几点问题想请教各位老师、前辈。
模型描述：二维鼓泡流化床，双组份颗粒密度，直径都不同，底部为速度进口，顶部为压力出口，气相壁面无滑移，颗粒相壁面有滑移。模拟参数设置与颗粒动力学的相关模型均依据已发表的文献，并结合FLUENT中的理论指导，重质颗粒堆积在底部，轻质颗粒在其上，颗粒配比1:1。气固相间曳力模型为Gidaspow，固固相间曳力模型为Syamlal-obrien-symmetric。
出现的问题：计算发散，先是颗粒温度发散，紧接着连续性方程发散。
尝试的解决措施：
1.改变边界表观气速，失败
2.改变曳力模型，失败
3.改变网格数量，失败
4.将颗粒温度求解方程由偏微分方程改为代数方程，成功（想和大家讨论的问题1）请教问题1：关于颗粒温度方程求解采用代数方程还是偏微分方程。代数方程忽略了颗粒温度输运方程中的对流与扩散项，偏微分方程求解的是完整的输运方程，从这个角度来说理应选择偏微分方程，但我无法调收敛，此外，也有为数不多的文献采用代数方程计算，请问各位老师对于双组份鼓泡流化床选择代数方程的结果可信度如何。
请教问题2：颗粒温度方程采用偏微分方程，考虑方程离散格式的选取，二阶迎风、QUICK格式不收敛-发散了，一阶迎风中间有几步没有收敛，后续没有发散（当然，由于一些原因，对于一阶迎风，我没有一直算下去），但对于一阶迎风的精度是否合适，想请教各位老师。
请教问题3：固固相间曳力模型采用Syamlal-obrien-symmetric，此模型中存在一个摩擦因数（i颗粒相与k颗粒相之间的摩擦），在fluent中，官方文档里说明fluent里采用的 摩擦因数值为0，但有文献采取不同的值来进行模拟，不过我并没有找到fluent里如何设置这个值，因此请各位老师指教。
非常欢迎和感激做过相关研究的老师，前辈给予指导，另外我的研究方向是气固两相流，所知十分有限，希望各位做欧拉-欧拉，DEM，MP-PIC的同志一起交流讨论。
第一次发帖，请多多指教。
